问题标签 [bessel-functions]

我想在python中计算修改后的贝塞尔函数的积分或导数。我想计算无限积分（无限制）。最近我找到了一种方法来做到这一点。您可以在下面看到一个简单函数 (x**2) 的示例：

结果是： x^3/3 。但是当我放置一个修改过的贝塞尔函数而不是 x**2 时：

在这种情况下，我收到此错误：

AttributeError：“模块”对象没有属性“iv”

需要说明的是，第一类贝塞尔函数的积分是零类贝塞尔函数。我希望得到：iv(0,x)

我怎么能在python中做到这一点？

我制作了一个使用以下在线库的 .cpp 代码：

Complex_Bessel_functions。

我想用我的代码制作一个 .mex 文件。当我打字时：

Matlab 给出以下错误：

这些错误发生的我的功能之一是：

关键是这个库将一些 Fortran 代码与 C++ 连接起来。例如，该sph_besselY函数将调用适当的 Fortran 代码zbesy_wrap。

我想以某种方式编译 Fortran 函数吗？我以与网页相同的方式安装库。另外，我正在将库导入到我的代码中，如下所示：

meijerG 函数是否可能包含负值（即 is {-1,0,0}）？我尝试了 Mathematica 和 Matlab 来计算这个 meijerG 函数，但是它们产生了一个错误，这meijerG不是为给定的参数定义的。

这是我的代码：

这里我还附上了文本中的函数图像。

我正在尝试使用以下等式计算 Ricean Fading PDF。Ricean 褪色 PDF . 其中“y”是归一化包络，“gamma”是 SNR

如果 K 值很大，那么

exp 产生大浮点数（eq 1.01e-5088）。在 python 中，它将显示 '0.0' 作为值

Bessel 函数的值显示大 int 值（eq 7.78e+5092）。在 python 中它会显示 '**inf**' 值

如何在 python 中存储大整数和浮点值并计算 pdf？

谢谢。

作为一个 SymPy 新手，我正在考虑以下 SymPy 表达式（它应该等于 1 if 0<x<1）：

lambda_m的m第-个零在哪里besselj(0,x)。

现在，mpmath有一个函数besseljzero(0,m,0)可以精确计算这个。

不幸的是，如果我在上面的公式中替换（手动，我的意思是......）lambda_m，besseljzero(0,m,0)SymPy 会给我一个错误，因为m它不是整数......

我想通过创建一个函数来解决这个问题，该函数应该：

• j0(m)如果m不计算为整数则返回
• besseljzero(0,m,0)如果确实如此，请返回

但我不知道如何进行。

这是一个好主意，有人可以帮助我吗？

我正在尝试绘制一个特殊（例如贝塞尔）函数的集成，我的最小代码如下。

但是，这给了我几条错误消息，说明该函数很容易在Mathematicaquadpack.error: Supplied function does not return a valid float中绘制。Python的贝塞尔函数有限制吗？

我已使用此文档进行绘图。

我试图在表面上积分以下方程：强度=（2*J1（z）/z）^2，z=A*sqrt（（x-mu1）^2+（y-mu2）^2）， A(L) x 和 y 的常数和 J1 一阶贝塞尔函数。为此，我使用 dblquad 函数，如下所示：

这里唯一重要的参数是极坐标中的 r 和 phi（其他参数取决于此处不重要的其他参数），其中 x=r cos(phi) 和 y=r sin(phi) 但是当我尝试整合函数时，我得到了这个信息：

C:\pyzo2013b\Lib\pyzo-packages\scipy\integrate\quadpack.py:289：UserWarning：已达到最大细分数（50）。如果增加限制没有改善，建议分析被积函数以确定困难。如果可以确定局部难度的位置（奇点、不连续性），则可能会从拆分区间并调用子范围上的积分器中获益。也许应该使用专用的积分器。
警告。警告（味精）

然后是一个完全不准确的结果，然后是：

C:\pyzo2013b\Lib\pyzo-packages\scipy\integrate\quadpack.py:289：UserWarning：积分可能是发散的，或者是慢慢收敛的。警告。警告（味精）

我确实理解此消息的含义，但我有两个问题：

• 除了将我的积分区间划分为更小的部分之外，我还有什么方法可以避免这种细分错误（我想通过将它们与无限域上的范数进行比较来检查我的其他结果，我将无法做到所以如果我不能正确整合无限域）？也许有一个特殊用途的集成商？但我不知道它是什么或如何使用它们。
• 为什么当 z 收敛到零时（就像窦基数一样），我知道 J1(z)/(z) 收敛到 1 时会收到关于发散积分或奇点的警告？

有人有答案吗？

这是完整的有用代码行（所有其他参数都以其他方式定义）：

（为了更好地理解该功能，我根据 gboffi 的建议对其进行了修改。）

我需要集成这种振荡函数：

我有一个正在振荡的贝塞尔函数，而 F 不是很振荡。我正在寻找在 C++ 中执行此操作的最精确/准确的方法。希望这应该是一个已经存在的实现，我可以通过例如库等使用它......

GSL 库可能是一种选择，但即使在这种情况下，您能否向我推荐许多可用的例程中哪一个对我最有用？

编辑：我知道这个 可能重复的问题的存在

但我在那里没有看到明确的答案。此外，除非有某种包装器，否则 Fortran 库对我不利。

我正在尝试对由多个 Bessel 函数（第一类和第二类）组成的实值被积函数进行数值评估。被积函数是振荡衰减的，需要在 0 到 +∞ 之间求值。到目前为止，我使用 scipy.integrate 子包（quad 和 fixed_quad）的尝试都没有成功。当实际上它应该是平滑的时，评估值会跳跃。对于某些参数值集，我还会收到警告：“IntegrationWarning：积分可能发散或缓慢收敛。” （已知是收敛的）或“IntegrationWarning：已达到最大细分数（50）。”

方程式来自：http ://dx.doi.org/10.1029/WR003i001p00241

它也可以在这里找到：http ://www.aqtesolv.com/papadopu.htm

感谢您在 python 中繁琐函数的数值集成方面提供的任何帮助...

代码示例

我是 julia 的新手，所以我欢迎一些建议来改进以下功能，

基本上，我不太确定在这两个参数（x 和 nu）上获取矩阵的推荐方法是什么。该文档没有提供太多信息，但我知道底层的 fortran 例程在内部循环 nu，所以为了性能，我宁愿不再这样做。

编辑： 有人问我目标；它是针对 $x$ 和 $\nu$ 的多个值计算 Riccati-Bessel 函数 $j_1(x,\nu),h_1(x,\nu)$。

我已经从原始版本中删除了风格问题，以专注于这个核心问题。