当前位置：网站首页>贝塞尔函数

贝塞尔函数

2022-08-09 14:52:00 【Chandler_river】

贝塞尔方程的导出
- 二维热传导模型
贝塞尔方程的广义幂级数解
- $y(x)=\sum_{k=0}^{\infty}a_kx^{s+k}(a_0 \neq 0)$
- 代入解出 $\left\{\begin{matrix} (s^2-n^2)a_0=0\\ [(s+1)^2-n^2]a_1=0\\ [(s+k)^2-n^2]a_k+a_{k-2}=0(k=2,3,...) \end{matrix}\right.$
- 讨论1 n不为整数(包括0)和半奇数,则s1-s2=2n也不为整数
  - s1=n
    - $\left\{\begin{matrix} a_1=0\\ a_k=-\frac{a_{k-2}}{k(2n+k)}(k=2,3,...) \\a_{2m}=(-1)^m\frac{a_0}{2^{2m}m!(n+1)(n+2)...(n+m)} \end{matrix}\right.$
    - 因为a_0是任意常数,取
      - $a_{2m}=(-1)^m\frac{1}{2^{n+2m}m!\Gamma(n+m+1)}$
    - 得到方程的一个特解Jn(x)(n阶第一类贝塞尔函数)
      - $J_n(x)=\sum_{m=0}^\infty a_{2m}x^{n+2m}=\sum_{m=0}^{\infty}(-1)^m\frac{x^{n+2m}}{2^{n+2m}m!\Gamma(n+m+1)}$
      - 由达朗贝尔判别法(lim u(m+1)/um )可知级数在实轴上绝对收敛
  - s2=-n
    - $\left\{\begin{matrix} a_1=0\\ a_k=-\frac{a_{k-2}}{k(k-2n)}\\ a_{2m}=(-1)^m\frac{a_0}{2^{2m}m!(1-n)(2-n)...(m-n)} \end{matrix}\right.$
    - 取 $a_0=\frac{1}{a^{-n}\Gamma(1-n)}$
    - 得到方程另外一个特解,记作(-n阶第一类贝塞尔函数)
      - $J_{-n}(x)=\sum_{m=0}^\infty (-1)^m\frac{x^{2m-n}}{2^{2m-n}m!\Gamma(1+m-n)}$
  - 因为n \neq -n 所以Jn(x) 与 J-n(x) 线性无关,方程的通解可以表示为
    - 如果令,可得一与Jn(x)线性无关的特解(第二类贝塞尔函数/诺伊曼函数)
      - $Y_n(x)=\frac{J_n(x)cosn\pi-J_{-n}(x)}{sinn\pi}$
- 讨论2 n为整数(包括0) 则s1-s2=2n也为整数同样可得方程的两个特解
  - $J_n(x)=\sum_{m=0}^\infty a_{2m}x^{n+2m}=\sum_{m=0}^{\infty}(-1)^m\frac{x^{n+2m}}{2^{n+2m}m!\Gamma(n+m+1)}$
  - $J_{-n}(x)=\sum_{m=0}^\infty (-1)^m\frac{x^{2m-n}}{2^{2m-n}m!\Gamma(1+m-n)}$
  - 利用Γ函数的递推公式,可推知此时J(-n)(x) 与 J(n)(x) 线性相关
    - 定义
      - 其中n为整数而a不为整数,可推出
      - 与J(n)(x)线性无关的Y(n)(x) $Y_n(x)=\frac{2}{\pi}J_n(x)(ln(\frac{x}{2})+C)-\frac{1}{\pi}\sum_{m=0}^{n-1}\frac{(n-m-1)!}{m!}(\frac{x}{2})^{2m-n}-\frac{1}{\pi}\sum_{m=0}^\infty(-1)^m\frac{1}{m!(n+m)!}(\frac{x}{2})^{n+2m}(\sum_{k=0}^{n+m-1}\frac{1}{1+k}+\sum_{k=0}^{m-1}\frac{1}{1+k})(n=1,2,3...)$
- 讨论3 n为半奇数时贝塞尔函数可以用初等函数表示
  - $J_{\frac{1}{2}}(x)=\sqrt{\frac{2}{\pi x}}sinx$
  - $\left\{\begin{matrix} J_{\frac{2m+1}{2}}=(-1)^m\sqrt{\frac{2}{\pi}}x^{m+\frac{1}{2}}(\frac{1}{x}\frac{d}{dx})^m(\frac{sinx}{x})\\ J_{-\frac{2m+1}{2}}=\sqrt{\frac{2}{\pi}}x^{m+\frac{1}{2}}(\frac{1}{x}\frac{d}{dx})^m(\frac{cosx}{x}) \end{matrix}\right.$
贝塞尔函数的递推公式
- $J_n(x)=\sum_{m=0}^{\infty}(-1)^m\frac{x^{n+2m}}{2^{n+2m}m!\Gamma(n+m+1)}$
- $Y_n(x\frac{J_n(x)cosn\pi-J_{-n}(x)}{sinn\pi}$
- $\left\{\begin{matrix} J_{n-1}(x)+J_{n+1}(x)=\frac{2n}{x}J_n(x)\\ J_{n-1}(x)-J_{n+1}(x)=2J'_n(x) \end{matrix}\right.$
- $\left\{\begin{matrix} Y_{n-1}(x)+Y_{n+1}(x)=\frac{2n}{x}Y_n(x)\\ Y_{n-1}(x)-Y_{n+1}(x)=2Y'_n(x) \end{matrix}\right.$

版权声明
本文为[Chandler_river]所创，转载请带上原文链接，感谢
https://blog.csdn.net/Chandler_river/article/details/126104201

边栏推荐

猜你喜欢

随机推荐