贝塞尔函数 - 码农知识堂

贝塞尔函数
- 贝塞尔方程的导出
  - 二维热传导模型
- 贝塞尔方程的广义幂级数解
  - $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.$
相关阅读:
【LeetCode】No.88. Merge Sorted Array -- Java Version
MySQL主从复制
 云原生-FRP内网穿透（详解）使用云服务器将内网集群服务暴露至公网（二）
QT之QScriptEngine的用法介绍
 Linux学习-53-定时任务管理(at、anacron命令)
SSM“点点通”餐饮点餐小程序-计算机毕业设计源码11264
2023最新版JavaSE教程——第7天：面向对象编程(进阶)
还不知道图片去水印怎么去？看看这篇文章就知道了
 PPT架构师架构技能图
 Maven 使用IDEA环境操作
原文地址：https://blog.csdn.net/Chandler_river/article/details/126104201