1 数学背景
Bessel方程
贝塞尔方程的形式为:
其中, v v v为参数,y为因变量,z为自变量。Bessel方程的解就是Bessel函数,Bessel函数有两种:
第一类Bessel函数,J = besselj(nu,z,scale)
y = { J ν ( z ) = ( z 2 ) ν ∑ ( k = 0 ) ∞ ( − z 2 4 ) k k ! Γ ( ν + k + 1 ) J − ν ( z ) = ( z 2 ) − ν ∑ ( k = 0 ) ∞ ( − z 2 4 ) k k ! Γ ( − ν + k + 1 ) y=\begin{cases} J_{\nu}(z)=\left(\frac{z}{2}\right)^{\nu} \sum_{(k=0)}^{\infty} \frac{\left(\frac{-z^{2}}{4}\right)^{k}}{k ! \Gamma(\nu+k+1)} \\ J_{-\nu}(z)=\left(\frac{z}{2}\right)^{-\nu} \sum_{(k=0)}^{\infty} \frac{\left(\frac{-z^{2}}{4}\right)^{k}}{k ! \Gamma(-\nu+k+1)} \end{cases} y=⎩⎪⎨⎪⎧Jν(z)=(2z)ν∑(k=0)∞k!Γ(ν+k+1)(4−z2)kJ−ν(z)=(2z)−ν∑(k=0)∞k!Γ(−ν+k+1)(4−z2)k
第一类Bessel修正函数,I = besseli(nu,z,scale)
y = { I ν ( z ) = ( z 2 ) ν ∑ ( k = 0 ) ∞ ( z 2 4 ) k k ! Γ ( ν + k + 1 ) I − ν ( z ) = ( z 2 ) − ν ∑ ( k = 0 ) ∞ ( z 2 4 ) k k ! Γ ( ν + k + 1 ) y=\begin{cases} I_{\nu}(z)=\left(\frac{z}{2}\right)^{\nu} \sum_{(k=0)}^{\infty} \frac{\left(\frac{z^{2}}{4}\right)^{k}}{k ! \Gamma(\nu+k+1)} \\ I_{-\nu}(z)=\left(\frac{z}{2}\right)^{-\nu} \sum_{(k=0)}^{\infty} \frac{\left(\frac{z^{2}}{4}\right)^{k}}{k ! \Gamma(\nu+k+1)} \end{cases} y=⎩⎪⎨⎪⎧Iν(z)=(2z)ν∑(k=0)∞k!Γ(ν+k+1)(4z2)k