亥姆霍兹方程是如何推导出球面波表达式的?
首先写下球面坐标系下的亥姆霍兹方程:
\frac{1}{r^2}\frac{\partial}{\partial r}\left(r^2+\frac{\partial u}{\partial r}\right)+\frac{1}{r^2\sin \theta}\frac{\partial}{\partial \theta}\left(\sin \theta \cdot \frac{\partial u}{\partial \theta}\right)+\frac{1}{r^2\sin^2\theta}\frac{\partial^2u}{\partial\varPhi^2}+k^2u=0
由于是球面坐标系,利用球谐函数分离变量作试探解,即代入下式到方程中:
u=R(r)Y ^m _l (\theta,\text{\o})
得到径向的方程为:
r21drd(r2drdR(r))+[k2−r2l(l+1)]R(r)=0
作一次标度变换,x=kr, y(x)=R(r) ,得到球贝塞尔方程:
x21dxd(xdxdy)+[1−x2l(l+1)]y=0
再作变换 y(x)=√xv(x) ,代回球贝塞尔方程,得到:
x1dxd(xdxdv)+[1−x2(l+1/2)2]v=0
这就是柱坐标和平面极坐标系下常见的贝塞尔方程,不过在柱坐标下常见的是整数阶的贝塞尔函数,而这里是 l+21 阶的贝塞尔方程,
显然可以定义:
jl(x)=√2xπJl+1/2(x)
nl(1,2)(x)=√2xπNl+1/2(x)
注:此函数在 x=0 处是发散的。
hl(1,2)(x)=jl(x)∓i⋅nl(x)
注:贝塞尔函数 J ,诺依曼函数 N 都是贝塞尔函数方程的解,可以通过级数展开来获得级数解,对于 J ,直接在原点处展开就可以,对于 N 要通过 J 进行构造。这两者是贝塞尔方程的两个线性无关解。
由此,亥姆霍兹方程的一般解就是:
u=\mathrm{\sum^\infty_{l=0}\sum^l_{m=-l}}A_{lm}j_l(kr)Y_l^m(\theta,\text{\o})+\mathrm{\sum^\infty_{l=0}\sum^l_{m=-l}}B_{lm}n_l(kr)Y_l^m(\theta,\text{\o})
A, B 由方程的边界条件和初始条件给定。这种展开的完备性由斯图姆刘维尔定理保证。特别地,对于 l=0 的情况,可以验证 j0(x)=xsinx, n0(x)=−xcosx,又因为 Y00=1 ,此时球汉克尔函数对应的解就是\displaystyle u = \frac{e^{\plusmn\mathrm{i}kr}}{r}这个最常见的形式。
参考:
- 通过亥姆霍兹方程是如何推导出球面波表达式的? - 知乎
- 用到的部分 Katex语法
$\displaystyle$
x=0∑100x
$\mathrm{d}$
dxdy
$\left( 1 + \frac{a}{b} \right)$
(1+ba)
$\plusmn$
\plusmn
$\cdot$
⋅
$\varPhi$
\varPhi
$\text{\o}$
\text{\o}
$\begin{aligned} a & = b +c \\ &= d \end{aligned}$
a=b+c=d