Struve 函数：Struve 函数 H0(z)、H1(z)、H0(z)-Y0(z)、H1(z)-Y1(z) 和修改的 Str...

在计算贝塞尔函数积分时通常会遇到 Struve 函数。 在本次提交中，Struve 函数 H0(z) 和 H1(z) 的前两个阶数是针对向量/矩阵复数参数计算的。 另外，提供了用于相应函数H0（z）-Y0（z）和H1（z）-Y1（z）的例程，这些例程对大z有用，并且对修改后的函数L0（z）和L1（z）也有用。 这些是高精度例程。 他们针对 mfun 进行了测试，并且 Matlab 的超几何函数和一致性为 14 位有效数字。 关于计算时间，对于大型矩阵参数，它们要快几个数量级。 与Mathematica结果进行比较时，请注意，在某些参数区域中，不同的Mathematica版本会给出不同的结果。 所使用的方法是 abs(z)<=16 的 Chebyshev 展开式和 abs(z)>16 的有理近似，后者在右半复平面上进行映射。 使用了 Matlab 例程 bessely 和 besselh。