1350 字
7 分钟
关于特殊微分方程的一些笔记

方程引入#

物理问题中,我们经常需要求解拉普拉斯方程、泊松方程、波动方程,这一些方程都有一个共同点,都有拉普拉斯算子,也就是 ∇2\nabla^{2},一般的直角坐标系的拉普拉斯算子可以直接使用分离变量法进行求解,那在柱坐标、球坐标中也是如此:

先来看球坐标的拉普拉斯算子:

∇2=1r2∂∂r(r2∂∂r)+1r2sin⁡θ∂∂θ(sin⁡θ∂∂θ)+1r2sin⁡2θ∂2∂φ2\nabla^{2} = \frac{1}{r^{2}}\frac{ \partial }{ \partial r } \left( r^{2}\frac{ \partial }{ \partial r } \right) + \frac{1}{r^{2}\sin \theta}\frac{ \partial }{ \partial \theta } \left( \sin \theta \frac{ \partial }{ \partial \theta } \right)+\frac{1}{r^{2}\sin^{2}\theta}\frac{ \partial ^{2} }{ \partial \varphi^{2} }

柱坐标的拉普拉斯算子:

∇2=1r∂∂r(r∂∂r)+1r∂2∂θ2+∂2∂z2\nabla ^{2} = \frac{1}{r}\frac{ \partial }{ \partial r } \left( r\frac{ \partial }{ \partial r } \right) + \frac{1}{r}\frac{ \partial^{2} }{ \partial \theta^{2} } + \frac{ \partial^{2} }{ \partial z^{2} }

拿拉普拉斯方程举例子,球坐标的拉普拉斯方程可以表示为,设解为 u(r,θ,φ)=R(r)Θ(θ)Φ(φ)u(r,\theta,\varphi) = R(r)\Theta(\theta)\Phi(\varphi)

1r2R∂∂r(r2∂R∂r)+1r2sin⁡θ Θ∂∂θ(sin⁡θ∂Θ∂θ)+1r2sin⁡2θ Φ∂2Φ∂φ2=0\frac{1}{r^{2}R}\frac{ \partial }{ \partial r } \left( r^{2}\frac{ \partial R }{ \partial r } \right) + \frac{1}{r^{2}\sin \theta \,\Theta}\frac{ \partial }{ \partial \theta } \left( \sin \theta \frac{ \partial \Theta }{ \partial \theta } \right)+\frac{1}{r^{2}\sin^{2}\theta \, \Phi}\frac{ \partial ^{2} \Phi }{ \partial \varphi^{2}} = 0

设:

∂2Φ∂φ2=−m2Φ∂∂r(r2∂R∂r)=l(l+1)R\begin{align} \frac{ \partial ^{2}\Phi }{ \partial \varphi^{2} } &= -m^{2}\Phi \\ \frac{ \partial }{ \partial r } \left( r^{2}\frac{ \partial R }{ \partial r } \right) &= l(l+1) R \end{align}

得到极角方程:

1sin⁡θ∂∂θ(sin⁡θ∂Θ∂θ)=(1sin⁡2θm2−l(l+1))Θ\frac{1}{\sin \theta}\frac{ \partial }{ \partial \theta } \left( \sin \theta \frac{ \partial \Theta }{ \partial \theta } \right) = \left(\frac{1}{\sin^{2}\theta}m^{2} - l(l+1)\right)\Theta

这就是球坐标下的分离变量法得到的三个分量方程,我们称极角方程为连带勒让德方程,当 m=0m = 0 方程简化为勒让德方程:

1sin⁡θ∂∂θ(sin⁡θ∂Θ∂θ)+l(l+1)Θ=0\frac{1}{\sin \theta}\frac{ \partial }{ \partial \theta } \left( \sin \theta \frac{ \partial \Theta }{ \partial \theta } \right) +l(l+1)\Theta = 0

径向方向可以化简为:

r2∂2R∂r2+2r∂R∂r−l(l+1)R=0r^{2}\frac{ \partial ^{2}R }{ \partial r^{2} } + 2r\frac{ \partial R }{ \partial r } - l(l+1)R = 0

这是欧拉方程,可以直接写出解:

Rl(r)=Al rl+Bl r−(l+1)R_{l}(r) = A_{l}\,r^{l} + B_{l}\,r^{-(l+1)}

方位角部分是一个简单的谐振子方程,可以直接写出解:

Φm(φ)=Amcos⁡mφ+Bmsin⁡mφ\Phi_{m}(\varphi) = A_{m}\cos m\varphi + B_{m}\sin m\varphi

另外方位角部分有一个天然的周期限制条件:

Φ(φ+2π)=Φ(φ)\Phi(\varphi + 2\pi) = \Phi(\varphi)

这要求 mm 必须是一个整数,记住这个限制条件,后续会用到的。

要解这个拉普拉斯方程,相当于解极角部分的连带勒让德方程。

勒让德方程#

要解连带勒让德方程,我们先来解简单一些的勒让德方程,也就是:

1sin⁡θ∂∂θ(sin⁡θ∂Θ∂θ)+l(l+1)Θ=0 \frac{1}{\sin \theta}\frac{ \partial }{ \partial \theta } \left( \sin \theta \frac{ \partial \Theta }{ \partial \theta } \right) +l(l+1)\Theta = 0

换元: x=cos⁡θ,y=Θ(θ)x= \cos \theta,\quad y = \Theta(\theta),得到:

11−x2y′′−2xy′+l(l+1)y=0\frac{1}{1-x^{2}}y^{\prime\prime} - 2 x y^{\prime} + l(l+1)y = 0

设级数解:y=∑CkxC+ky = \sum C_{k}x^{C+k},代入得到系数递推方程:

C=0Ck+2=(−1)(l−k)(l+k+1)(k+2)(k+1)Ck\begin{align} C & = 0 \\ C_{k+2} &=(-1) \frac{(l-k)(l+k+1)}{(k+2)(k+1)}C_{k} \end{align}

发现偶数列和奇数列无关,那么解就可以表示成奇数列和偶数列的线性组合,对于任何一个 ll ,都存在数列满足递归。但是这里有一个问题,一般来看,如果这是一个无穷项级数,那么可以证明的是在 x=1x = 1 处级数趋向正无穷,这显然和物理事实不相符,那要让级数是一个有界值,必须让系数在某一点截断,也就是在此之后的系数递归均为0,也就是级数变成了一个多项式。要想实现这样的效果,必须要求 ll 为一个整数。比如 l=2l = 2 ,那么偶数列会在 C2C_{2} 处截断,即 C4=C6=⋯=0C_{4} = C_{6} = \dots= 0

那么我们就确定了要想获得在物理意义上合理的数学解,就需要 ll 为一个整数,当 ll 为偶数,那为了保证级数有界,奇数列必须全为 00 。

比如我们可以写出前面几个解:

P0=1P1=xP2=3x2−1P3=5x3−3x\begin{alignat}{2} P_{0} &=1 \qquad \qquad \qquad & &P_{1} =x\\ P_{2} &=3x^{2}-1 & &P_{3} = 5x^{3} - 3x \end{alignat}

如果告知了最高次幂的系数,那么整一个数列都是已知的了。那为了简化表达式,我们有一个规定:

Cl=(2l)!2l(l!)2C_{l} = \frac{(2l)!}{2^{l}(l!)^{2}}

后续计算出来的多项式满足 Pl(1)=1P_{l}(1) = 1,满足有界性。我们称这样子获得的多项式称为勒让德多项式:

P0=1P1=xP2=12(3x2−1)P3=12(5x3−3x)\begin{alignat}{2} P_{0} &= 1 \qquad \qquad\qquad & P_{1} &= x \\ P_{2} &= \frac{1}{2}(3x^{2}-1) & P_{3} &= \frac{1}{2}(5x^{3}-3x) \end{alignat}

当我们定义内积为:f⋅g=∫−11f(x)g(x)dxf \cdot g = \int_{-1}^{1}f(x)g(x)\mathrm{d}x 时,勒让德多项式是成交的:

∫−11PnPm dx=22n+1δmn\int _{-1}^{1}P_{n}P_{m}\, dx = \sqrt{ \frac{2}{2n+1} }\delta_{mn}

另外我们也可以证明勒让德多项式也是一个完备的函数集,可以用勒让德函数来表示任何平方可积函数,类似傅里叶级数。

回到勒让德方程,我们写出解:

Θ(θ)=AlPl(cos⁡θ)\Theta(\theta) = A_{l}P_{l}(\cos \theta)

连带勒让德函数#

我们原本的方程为

1sin⁡θ∂∂θ(sin⁡θ∂Θ∂θ)=(1sin⁡2θm2−l(l+1))Θ\frac{1}{\sin \theta}\frac{ \partial }{ \partial \theta } \left( \sin \theta \frac{ \partial \Theta }{ \partial \theta } \right) = \left(\frac{1}{\sin^{2}\theta}m^{2} - l(l+1)\right)\Theta

这是连带勒让德方程,其中 mm 是一个整数。利用相同的换元得到化简之后的方程:

11−x2y′′−2xy′+(l(l+1)−m21−x2)y=0\frac{1}{1-x^{2}}y^{\prime\prime} - 2 x y^{\prime} + \left( l(l+1) - \frac{m^{2}}{1-x^{2}} \right)y = 0

为了解这个方程,我们需要先对勒让德方程求 mm 阶导数,然后进行换元,最后得到解:

Plm(x)=(−1)∣m∣(1−x2)∣m∣/2P l(m)(x)P_{l}^{m}(x) = (-1)^{|m|}(1-x^{2})^{|m|/2}P_{\text{ }l}^{(m)}(x)

带回到原本的方程中:

Θ(θ)=AlPlm(cos⁡θ)\Theta(\theta) = A_{l}P_{l}^{m}(\cos \theta)

球谐函数#

将方位角和极角两部分结合起来就得到了球谐函数 Ylm(θ,φ)Y_{l}^{m}(\theta ,\varphi)

Ylm(θ,φ)=εAlmPlm(cos⁡θ)eimφY^{m}_{l}(\theta, \varphi) = \varepsilon A_{lm}P_{l}^{m}(\cos \theta)e^{im\varphi}

其中 AlmA_{lm} 是归一化系数,ε\varepsilon 是相位规定系数

拉普拉斯方程的解#

最后和径向部分结合起来得到拉普拉斯方程的通解:

y(r,θ,φ)=∑l,m(Al rl+Bl r−(l+1))Ylm(θ,φ)y(r,\theta,\varphi) = \sum_{l,m}(A_{l}\,r^{l}+B_{l}\,r^{-(l+1)})Y_{lm}(\theta,\varphi)
关于特殊微分方程的一些笔记
https://fuwari.bianwenzhao.workers.dev/posts/贝塞尔函数和勒让德函数/
作者
Siesta
发布于
2026-10-12
许可协议
CC BY-NC-SA 4.0