跳转至

Sturm-Liouville 本征值理论

Sturm–Liouville 特征值理论

本页数学主线使用「特征值、特征函数、微分算子」;量子力学应用段使用「本征值、本征态、算符」.两套叫法对应不同语境,不改变特征方程本身.

物理问题引入:为什么所有特殊函数都具有正交完备性?

在学习数理方法时,学生通常会被迫背诵一大批看似孤立的函数正交公式:

  • 傅里叶三角函数的正交性:∫0Lsin⁡(nπxL)sin⁡(mπxL)dx∝δnmint_0^L sinleft(npi x/Lright)sinleft(mpi x/Lright)dx propto delta_nm;
  • 勒让德多项式的正交性:∫−11Pn(x)Pm(x)dx∝δnmint_-1^1 P_n(x) P_m(x)dx propto delta_nm;
  • 贝塞尔函数的带权正交性:∫0aJν(knr)Jν(kmr)rdr∝δnmint_0^a J_nu(k_n r) J_nu(k_m r) , r , dr propto delta_nm;
  • 厄米多项式与拉盖尔多项式的加权正交性……

** 为什么大自然在不同坐标系下分离变量得到的每一个方程,其解都奇迹般地能够作为「基底」展开任意物理场?** 夏尔·施图姆(Charles Sturm) 与约瑟夫·刘维尔(Joseph Liouville) 在 19 世纪给出了终极答案:它们全都是同一种具备「自共轭性」的微分算子的特征函数!


1. Sturm–Liouville 算子与标准方程

定义在有限或无穷区间 [a,b][a, b] 上的标准 Sturm–Liouville (S-L) 特征值方程 为:

Ly(x)=λw(x)y(x)L y(x) = lambda w(x) y(x)

其中微分算子 LL 定义为:

L=−ddx[p(x)ddx]+q(x)L = -ddxleft[ p(x) ddx right] + q(x)
  • p(x)>0,w(x)>0p(x) > 0, w(x) > 0:实值连续函数;w(x)w(x) 称为 权函数(Weight Function);
  • q(x)q(x):实值连续位势项;
  • λlambda:待确定的 特征值;y(x)y(x):对应的非零 特征函数.

2. 拉格朗日恒等式与自共轭性(Self-Adjointness)

考虑两个任意光滑函数 u(x)u(x) 与 v(x)v(x),直接计算双线性差式:

uLv−vLu=−ddx[p(x)(udvdx−vdudx)]u Lv - v Lu = -ddxleft[ p(x)left( u dvdx - v dudx right) right]

在区间 [a,b][a, b] 上积分,得到著名的 格林公式/拉格朗日恒等式:

∫ab(uLv−vLu)dx=−p(x)(uv′−vu′)|abint_a^b left( u Lv - v Lu right) dx = -left. p(x) left( u v' - v u' right) right|_a^b

自共轭边界条件

若施加的边界条件使得上式右端的边界项为零,即:

p(x)W(u,v)(x)|ab=0left. p(x) W(u, v)(x) right|_a^b = 0

则算子 LL 称为在给定位形空间中是 自共轭(自伴/Hermitian) 的!满足此条件的常见物理边界包括:

  1. 齐次边界条件(Dirichlet、Neumann 或 Robin):α1y(a)+β1y′(a)=0alpha_1 y(a) + beta_1 y'(a) = 0 且 α2y(b)+β2y′(b)=0alpha_2 y(b) + beta_2 y'(b) = 0;
  2. 周期性边界条件:y(a)=y(b)y(a) = y(b) 且 y′(a)=y′(b)y'(a) = y'(b),且 p(a)=p(b)p(a) = p(b);
  3. 自然奇点边界:在端点处 p(a)=0p(a) = 0 或 p(b)=0p(b) = 0(如勒让德方程在 x=±1x = +/- 1 处 p(x)=1−x2=0p(x) = 1-x^2 = 0),只要解在该点保持有界,边界项自动消失!

3. Sturm–Liouville 定理的三大皇冠推论

当 LL 自共轭时,数学分析定理给出了物理学家最梦寐以求的三大结论:

3.1 特征值必为纯实数

设特征函数为复函数 y(x)y(x),对应的特征值为 λlambda:Ly=λwyLy = lambda w y.取共轭 Ly∗=λ∗wy∗Ly^* = lambda^* w y^*. 代入自共轭拉格朗日恒等式:

(λ−λ∗)∫ab|y(x)|2w(x)dx=0(lambda - lambda^*) int_a^b |y(x)|^2 w(x) dx = 0

因为 w(x)>0w(x) > 0 且 y(x)≢0y(x) notequiv 0,积分严格为正实数,故强制要求:

λ=λ∗⟹λ∈Rlambda = lambda^* implies lambda in R

3.2 不同特征值对应的特征函数严格带权正交

设 Lyn=λnwynLy_n = lambda_n w y_n,Lym=λmwymLy_m = lambda_m w y_m(λn≠λmlambda_n != lambda_m):

(λn−λm)∫abyn(x)ym(x)w(x)dx=0(lambda_n - lambda_m) int_a^b y_n(x) y_m(x) w(x) dx = 0

因为 λn≠λmlambda_n != lambda_m,必须满足:

∫abyn(x)ym(x)w(x)dx=Nnδnmint_a^b y_n(x) y_m(x) w(x) dx = N_n , delta_nm

其中归一化模长平方 Nn=∫ab[yn(x)]2w(x)dxN_n = int_a^b [y_n(x)]^2 w(x)dx.

3.3 特征函数系的完备性(Completeness)

所有的特征函数 {yn(x)}n=1∞y_n(x) _n=1^inf 构成了函数空间 Lw2[a,b]L_w^2[a, b] 的一组 完备正交基! 区间内任意满足狄利克雷条件的物理态 f(x)f(x),都可以唯一展开为广义傅里叶级数:

f(x)=∑n=1∞cnyn(x),cn=1Nn∫abf(x)yn(x)w(x)dxf(x) = sum_n=1^inf c_n y_n(x), quad c_n = 1N_n int_a^b f(x) y_n(x) w(x) dx

4. 全系特殊函数在 S-L 框架下的统一大对照

经典方程区间 [a,b][a, b]算子参数 p(x)p(x)位势项 q(x)q(x)权函数 w(x)w(x)特征函数系 yn(x)y_n(x)
傅里叶级数[0,L][0, L]110011三角函数 sin⁡(nπxL)sinleft(npi x/Lright)
勒让德多项式[−1,1][-1, 1]1−x21 - x^20011勒让德多项式 Pl(x)P_l(x)
连带勒让德[−1,1][-1, 1]1−x21 - x^2m21−x2m^21-x^211连带勒让德函数 Plm(x)P_l^m(x)
贝塞尔方程[0,R][0, R]xxν2xnu^2xxx贝塞尔函数 Jν(knx)J_nu(k_n x)
切比雪夫多项式[−1,1][-1, 1]1−x21-x^20011−x211-x^2切比雪夫多项式 Tn(x)T_n(x)
厄米多项式(量子谐振子)(−∞,∞)(-inf, inf)e−x2e^-x^200e−x2e^-x^2厄米多项式 Hn(x)H_n(x)
拉盖尔多项式(氢原子)[0,∞)[0, inf)xe−xx e^-x00e−xe^-x连带拉盖尔多项式 Lnk(x)L_n^k(x)

5. 跨入现代量子力学的数学接口

在有限维线性代数中,实对称矩阵 A=ATA = A^T 的特征值为实数,且不同特征值的特征向量正交; 在无限维希尔伯特空间中,自共轭 Sturm–Liouville 微分算子 LL 正是实对称矩阵向 连续函数空间的无限维自然推广! 量子力学中:

  • 哈密顿算符 H^=−ℏ22m∇2+V(r)H = -hbar^22mnabla^2 + V(r) 是自共轭(Hermitian)的;
  • 定态薛定谔方程 H^ψn=EnψnHpsi_n = E_npsi_n 正是 S-L 本征值问题;
  • 能量本征值 EnE_n 为可观测的实数,波函数正交归一化 ∫ψn∗ψmdV=δnmint psi_n^* psi_m dV = delta_nm,任意量子态叠加 |Ψ⟩=∑cn|ψn⟩|Psirangle = sum c_n |psi_nrangle——这一切全部直接继承自 Sturm–Liouville 理论!

6. 学习衔接

  • 下一节:利用特征函数展开法求解任意非齐次点源方程,进入 格林函数方法;
  • 特殊函数各论:在 特殊函数族导学 中细化各方程的具体解法与几何性质.


评论