python贝塞尔函数-Python中使用大指数的Bessel函数.pdf

python贝塞尔函数_Python中使⽤⼤指数的Bessel函数

mpmath是⼀个⾮常棒的库，是进⾏⾼精度计算的径。值得注意的是，这些函数可以从它们更基本的组成部分计算出来。因此，您不必遵

守scipy的限制，您可以使⽤不同的⾼精度库。下⾯是最简单的例⼦：importnumpyasnp

fromscipy.specialimport*

X=np.random.random(3)

v=2.000000000

printBesselFunctionJ

printjn(v,X)

printModifiedBesselFunction,Iv

print((1j**(-v))*jv(v,1j*X)).real

printiv(v,X)

printModifiedBesselFunctionofthesecondkind,Kv

print(iv(-v,X)-iv(v,X))*(np.pi/(2*sin(v*pi)))

printkv(v,X)

printModifiedsphericalBesselFunction,in

printnp.sqrt(np.pi/(2*X))*iv(v+0.5,X)

print[sph_in(floor(v),x)[0][-1]forxinX]

printModifiedsphericalBesselFunction,kn

printnp.sqrt(np.pi/(2*X))*kv(v+0.5,X)

print[sph_kn(floor(v),x)[0][-1]forxinX]

printModifiedsphericalBesselFunction,in

printnp.sqrt(np.pi/(2*X))*iv(v+0.5,X)

print[sph_in(floor(v),x)[0][-1]forxinX]

这就提供了：BesselFunctionJ

[0.018870980.001842020

ModifiedBesselFunction,Iv

[0.019358080.001846560

[0.019358080.001846560

ModifiedBesselFunctionofthesecondkind,Kv

[12.614948641352

[12.614948651352

ModifiedsphericalBesselFunction,in

[0.010305600

[0.010305631072943869,0.00098466280846548084,0.050033450286650107]

ModifiedsphericalBesselFunction,kn

[76.867386312622.982284116.99803515]

[76.867205587011171,2622.9730878542782,6.998023749439338]

ModifiedsphericalBesselFunction,in

[0.010305600

[0.010305631072943869,0.00098466280846548084,0.050033450286650107]

除⾮底层数据具有⾼精度，否则对于要查找的⼤值，这将失败。