问题标签 [bessel-functions]

我想估计涉及贝塞尔函数和积分的函数参数。但是，当我尝试运行它时，我收到一条消息“f(x, ...) 中的错误：找不到函数“BesselI””。我不知道要解决它，并会感谢任何相关的建议。

Newton-Raphson 最大化迭代次数：0 返回代码：100 初始值超出范围。

我得到“有 50 个或更多警告（使用警告（）查看前 50 个）”，当我使用警告（）时，以下是警告。

希望得到一些关于我的代码的指导，我对 python 还是很陌生，并且在我的代码中修改后的贝塞尔函数 k0 时遇到了困难。实际代码分为两部分

1）第一节使用sympy求解t，返回t=0时x0的值。

2) 第二部分使用这个值 x0 来求解不同 x 值的 t。

在第二步，我遇到了困难，scipy.special.k0 函数似乎不适用于所述已解决的 x0 值，我不知道为什么

当我运行代码时，它会导致以下错误

非常感谢您的任何帮助，并在此先感谢您。

在 2.7 中学习 Pythonic。有没有办法避免显式循环？答案=[5, 4, 4, 3, 3, 2]

假设圆柱体被加热到非常高的温度并被淬火成淬火剂，其中周围温度为非零温度，即T0，圆柱体的初始温度为Ti。

该等式是针对给定的导热系数和传热系数，在给定的时间 t 找出表面和中心之间任意一点的温度。

J0 是第一类贝塞尔函数，J1=-d(J0)/dz

其中特征值 beta(m) 是从超越方程的根中获得的

如何使用上述等式绘制冷却曲线？

我是 matlab 新手，这完全不符合我的要求。在其他任何问题中，我都没有看到类似的情况。谢谢你的帮助。

我进行了以下计算：

可以理解的是，如果不是更多的话，这需要一些时间 > 10 秒（没有计时）。

随后的调用速度明显更快，这让我认为它缓存了结果。我想知道还有哪些其他功能/计算 sympy cahces？

这是否意味着如果我想计算mp.besseljzero(n, m)n < N 和 m < M 最好计算mp.besseljzero(N-1, M-1)然后访问其他结果？

提前致谢。

我应该补充一点，我在 sympy FAQ 上找到了这个信息。似乎支持某种缓存。如果有更多文档，您可以指出我会有所帮助！

https://github.com/sympy/sympy/wiki/Faq

方程是I ₀ (a*x)=b，其中a,b是常数，I ₀ (y)是零阶贝塞尔函数。

I ₀ (x)=1/(2*pi)*integral(x*cos(t)) dt [从 0 ->2*pi]

我想在给出 a 和 b 时得到 x 的值。我只想要结果，所以不必全部通过Matlab来实现求解过程，一个近似的答案也可以。

我正在尝试使用零阶的MATLAB Besselbesselj(0,Z)函数 ---> 。我注意到的是：它可以处理 n 次方的实部（例如），但如果它在复杂参数中e100大于 ( ) 则不能处理虚部。e2

例如：(2+2*i)它将处理的复杂参数。但它不能处理(20000+20000*i)----->给出的答案为“无穷大”。

这就是我得到的：

有人可以解释一下吗？虚部的功率有限制吗？

我注意到scipy.specialn 阶贝塞尔函数和参数 xjv(n,x)在 x 中被矢量化：

In [14]: import scipy.special as sp In [16]: sp.jv(1, range(3)) # n=1, [x=0,1,2] Out[16]: array([ 0., 0.44005059, 0.57672481])

但是没有相应的球面贝塞尔函数的矢量化形式sp.sph_jn：

此外，球面贝塞尔函数一次计算 N 的所有阶数。因此，如果我想要n=5Bessel 函数作为参数x=10，它会返回 n=1,2,3,4,5。它实际上一次返回 jn 及其导数：

为什么 API 中存在这种不对称性，有没有人知道一个库会返回向量化的球形贝塞尔函数，或者至少更快（即在 cython 中）？

我正在尝试在 Python 中绘制 FM 信号的 Bessel 等效项的频域。

我遇到的第一个问题是提高 FFT 的分辨率以查看较小的带宽 - 我相信我已经通过将贝塞尔波函数的总和与一个大小为阵列的两倍的大型零填充数组连接起来解决了这个问题保存波形数据。

但是，这会对 Y 轴的缩放产生影响——当我增加零填充的大小时，Y 轴的大小会下降。我可以将 Y 轴乘以什么因素来反转这个？到目前为止的实验使我陷入了死胡同...

非常感谢！

我正在使用 scipy 库来查找 Bessel 函数的零点。它似乎适用于低阶贝塞尔函数

但是，在切断后，它挂起

上面的行永远不会返回。这是 scipy 的预期行为吗？