正交曲线坐标系中的分离变量
这一节介绍了柱坐标系及球坐标系下亥姆霍兹方程如何进行分离变量，通过变量的分离，我们顺理成章地引入了勒让德函数、贝塞尔函数、球谐函数等数理方程常用函数。读者可以先行浏览一遍了解方程推导过程，待特殊函数掌握之后再来仔细回忆回味。
柱坐标系下的亥姆霍兹方程
柱坐标系下 { ∇ 2 u + λ u = 0 u = R ( ρ ) Φ ( φ ) Z ( z ) \left\{\begin{array}{l}\nabla^2u+\lambda u=0\\u=R(\rho)\Phi(\varphi)Z(z)\end{array}\right. {∇2u+λu=0u=R(ρ)Φ(φ)Z(z)
在数学物理方程相关介绍中我们讨论过，亥姆霍兹方程中的常数 λ \lambda λ与时间有关。
带入柱坐标系中
∇
2
u
\nabla^2 u
∇2u的表达式则有
1
ρ
R
d
d
ρ
(
ρ
d
R
d
ρ
)
+
1
ρ
2
(
1
Φ
d
2
Φ
d
φ
2
)
+
λ
=
−
1
Z
d
2
Z
d
z
2
=
μ
\frac{1}{\rho R} \frac{\mathrm{d}}{\mathrm{d} \rho}\left(\rho \frac{\mathrm{d} R}{\mathrm{~d} \rho}\right)+\frac{1}{\rho^{2}}\left(\frac{1}{ \Phi} \frac{\mathrm{d}^{2} \Phi}{\mathrm{d} \varphi^{2}}\right)+\lambda=-\frac{1}{Z} \frac{\mathrm{d}^{2} Z}{\mathrm{~d} z^{2}}=\mu
ρR1dρd(ρ dρdR)+ρ21(Φ1dφ2d2Φ)+λ=−Z1 dz2d2Z=μ
由于自然边界条件限制，
1
Φ
d
2
Φ
d
φ
2
\frac{1}{ \Phi} \frac{\mathrm{d}^{2} \Phi}{\mathrm{d} \varphi^{2}}
Φ1dφ2d2Φ 只能取
−
n
2
-n^2
−n2 (
n
n
n 为正整数)，令
λ
−
μ
=
k
2
\lambda-\mu=k^2
λ−μ=k2，则
ρ
2
R
′
′
+
ρ
R
′
+
(
k
2
ρ
2
−
n
2
)
R
=
0
\rho^{2} R^{\prime \prime}+\rho R^{\prime}+\left(k^{2} \rho^{2}-n^{2}\right) R=0
ρ2R′′+ρR′+(k2ρ2−n2)R=0
综上所述，通过在柱坐标中分离变量，解偏微分方程
∇
2
u
+
λ
u
=
0
\nabla^2u+\lambda u=0
∇2u+λu=0 的问题，就化为了解三个常微分方程
{
Z
′
′
+
μ
Z
=
0
Φ
′
′
+
n
2
Φ
=
0
ρ
2
R
′
′
+
ρ
R
′
+
(
k
2
ρ
2
−
n
2
)
R
=
0
径
向
方
程
\left\{\begin{array}{l} Z^{\prime \prime}+\mu Z=0 \\ \Phi^{\prime \prime}+n^{2} \Phi=0 \\ \rho^{2} R^{\prime \prime}+\rho R^{\prime}+\left(k^{2} \rho^{2}-n^{2}\right) R=0\quad 径向方程 \end{array}\right.
⎩⎨⎧Z′′+μZ=0Φ′′+n2Φ=0ρ2R′′+ρR′+(k2ρ2−n2)R=0径向方程
的问题，其中
μ
、
n
2
、
k
2
\mu 、 n^{2} 、 k^{2}
μ、n2、k2 都是在分离变量过程中所引人的常数，与弦的振动问题一样，它们不能任意取值，而要根据边界条件取某些特定的值，它们称为上述方程的本征值。
方程
Z
′
′
+
μ
Z
=
0
Z^{\prime \prime}+\mu Z=0
Z′′+μZ=0 和
Φ
′
′
+
n
2
Φ
=
0
\Phi^{\prime \prime}+n^{2} \Phi=0
Φ′′+n2Φ=0 是常系数常微分方程，其解易于得到。方程
ρ
2
R
′
′
+
ρ
R
′
+
(
k
2
ρ
2
−
n
2
)
R
=
0
\rho^{2} R^{\prime \prime}+\rho R^{\prime}+\left(k^{2} \rho^{2}-n^{2}\right) R=0
ρ2R′′+ρR′+(k2ρ2−n2)R=0 是变系数常微分方程，若作变换
x
=
k
ρ
,
y
(
x
)
=
R
(
ρ
)
x=k \rho, y(x)=R(\rho)
x=kρ,y(x)=R(ρ)，则可化为
x
2
y
′
′
+
x
y
′
+
(
k
2
−
n
2
)
y
=
0
x^{2} y^{\prime \prime}+x y^{\prime}+\left(k^{2}-n^{2}\right) y=0
x2y′′+xy′+(k2−n2)y=0
称之为
n
n
n 阶贝塞尔 (Bessel) 方程，Bessel 方程的解为
y
n
(
x
)
=
c
n
J
n
(
x
)
+
d
n
N
n
(
x
)
y_n(x)=c_nJ_n(x)+d_nN_n(x)
yn(x)=cnJn(x)+dnNn(x)。它的求解和解的性质将在特殊函数中讨论。
本征值 k 2 = λ − μ k^2=\lambda-\mu k2=λ−μ， λ \lambda λ是与时间有关的系数， λ = 0 \lambda=0 λ=0表示方程与时间无关，代表稳定场； μ \mu μ与 z z z轴有关， μ = 0 \mu=0 μ=0表示方程是二维问题或具有 z z z方向平移不变性。若在二维稳定场内，或 λ − μ = k 2 \lambda-\mu=k^2 λ−μ=k2恰好为0，则径向方程变为欧拉方程 x 2 y ′ ′ + x y ′ − n 2 y = 0 x^{2} y^{\prime \prime}+x y^{\prime}-n^{2} y=0 x2y′′+xy′−n2y=0；若 k 2 > 0 k^2>0 k2>0，则方程为Bessel 方程；若 k 2 < 0 k^2
