贝塞尔函数 - 码农知识堂

贝塞尔函数
- 贝塞尔方程的导出
  - 二维热传导模型
- 贝塞尔方程的广义幂级数解
  - $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_0=\frac{1}{2^n\Gamma(n+1)}$
    $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)}$
    得到方程另外一个特解,记作 $J_{-n}(x)$ (-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) 线性无关,方程的通解可以表示为 $y=AJ_n(x)+BJ_{-n}(x)$
    如果令 $A=cotn\pi,B=-cscn\pi$ ,可得一与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) 线性相关
    定义 $Y_n(x)=\lim_{a->n}\frac{J_a(x)cosa\pi-J_{-a}(x)}{sina\pi}$
    其中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.$
相关阅读:
基于Python爬虫的股票成交量数据抓取分析系统
 数据库事务，JDBC操作和数据类型
 ASEMI整流桥KBPC3510，KBPC3510封装，KBPC3510应用
 2023.10.17 关于 wait 和 notify 的使用
 《nginx》二、nginx反向代理
 .NET Aspire 预览版 6 发布
 java毕业设计传统文化网站mybatis+源码+调试部署+系统+数据库+lw
javaee SpringMVC中json的使用
 2022/11/21[指针] 多维数组与指针的联系
 英国开发者年龄歧视为29岁，女程序员幸福指数略高于男性 | 全球开发者幸福指数报告
原文地址：https://blog.csdn.net/Chandler_river/article/details/126104201