幂级数解法与Frobenius方法

幂级数解法与 Frobenius 方法

物理问题引入：为什么球坐标和柱坐标总会引出变系数方程？

在直角笛卡尔坐标系中，拉普拉斯算子是简单的常系数二阶偏导之和 $\nabla^{2} = \frac{\partial^{2}}{\partial x^{2}} + \frac{\partial^{2}}{\partial y^{2}} + \frac{\partial^{2}}{\partial z^{2}}$ nabla^2 = d^2/d x^2 + d^2/d y^2 + d^2/d z^2．然而，真实物理系统往往具备球对称性（如行星引力场、氢原子库仑场）或柱对称性（如圆柱波导、同轴电缆）．当我们在曲线正交坐标系下分离变量时，几何面积与体积元的坐标依赖性使得方程 不可避免地带有变量系数：

球坐标角向方程（勒让德方程）： $(1 - x^{2}) y^{″} - 2 x y^{'} + l (l + 1) y = 0$ (1 - x^2) y'' - 2x y' + l(l+1)y = 0；
柱坐标径向方程（贝塞尔方程）： $x^{2} y^{″} + x y^{'} + (x^{2} - ν^{2}) y = 0$ x^2 y'' + x y' + (x^2 - nu^2)y = 0．

面对这些无法写出初等函数解析解的方程，幂级数展开是探索其精确解的基石工具．

1. 常点与奇点的严格判别

将二阶齐次线性常微分方程化为标准形式：

y^{″} + P (x) y^{'} + Q (x) y = 0

y'' + P(x) y' + Q(x) y = 0

1.1 常点（Ordinary Point）

若系数函数 $P (x)$ P(x) 和 $Q (x)$ Q(x) 在点 $x_{0}$ x_0 处解析（即能在 $x_{0}$ x_0 展开为收敛的泰勒级数），则称 $x_{0}$ x_0 为方程的常点．

1.2 奇点与正则奇点（Regular Singular Point）

若 $x_{0}$ x_0 处 $P (x)$ P(x) 或 $Q (x)$ Q(x) 不解析，则称 $x_{0}$ x_0 为奇点．若奇点的发散程度被严格约束在一定阶数以内，即：

(x - x_{0}) P (x) 和 (x - x_{0})^{2} Q (x)

(x - x_0) P(x) quad 和 quad (x - x_0)^2 Q(x)

在 $x_{0}$ x_0 处均解析，则称 $x_{0}$ x_0 为 正则奇点（Regular Singular Point）；否则称为非正则奇点．

物理经典方程	标准形式系数	关键点判别
勒让德方程	$P (x) = \frac{- 2 x}{1 - x^{2}}, Q (x) = \frac{l (l + 1)}{1 - x^{2}}$ P(x) = -2x1-x^2, quad Q(x) = l(l+1)1-x^2	$x = 0$ x=0 为常点； $x = \pm 1$ x = +/- 1（对应物理球面的南北极点 $θ = 0, π$ theta = 0, pi）为正则奇点
贝塞尔方程	$P (x) = \frac{1}{x}, Q (x) = \frac{x^{2} - ν^{2}}{x^{2}}$ P(x) = 1x, quad Q(x) = x^2 - nu^2x^2	$x = 0$ x=0（柱坐标对称轴原点）为正则奇点

2. 常点处的普通幂级数法

在常点 $x_{0}$ x_0（通常取 $x_{0} = 0$ x_0 = 0）附近，方程必定存在两个线性无关的解析解，可直接设为幂级数：

y (x) = \sum_{n = 0}^{\infty} a_{n} x^{n}

y(x) = sum_n=0^inf a_n x^n

逐项求导代入方程：

y^{'} (x) = \sum_{n = 1}^{\infty} n a_{n} x^{n - 1}, y^{″} (x) = \sum_{n = 2}^{\infty} n (n - 1) a_{n} x^{n - 2}

y'(x) = sum_n=1^inf n a_n x^n-1, quad y''(x) = sum_n=2^inf n(n-1) a_n x^n-2

通过统一 $x^{k}$ x^k 的幂次指标，令每一阶幂次的合并系数为零，即可得到关于展开系数 $a_{n}$ a_n 的 递推关系（Recurrence Relation）．

3. 弗罗贝尼乌斯方法（Frobenius Method）

在正则奇点 $x_{0} = 0$ x_0 = 0 附近，普通幂级数可能无法求出全部通解．弗罗贝尼乌斯法引入广义幂级数形式（假定 $a_{0} \neq 0$ a_0 != 0）：

y (x) = x^{r} \sum_{n = 0}^{\infty} a_{n} x^{n} = \sum_{n = 0}^{\infty} a_{n} x^{n + r}

y(x) = x^r sum_n=0^inf a_n x^n = sum_n=0^inf a_n x^n+r

其中指标 $r$ r 是待定的实数或复数．

3.1 指标方程（Indicial Equation）

将广义级数代入微分方程，考察最低幂次项 $x^{r}$ x^r 的系数，令其为零得到关于 $r$ r 的二次代数方程——指标方程：

I (r) = 0 ⟹ r_{1}, r_{2}

I(r) = 0 implies r_1, r_2

根据两根之差 $r_{1} - r_{2}$ r_1 - r_2 的性质分为三种情形：

两根不相差整数：直接得到两个线性无关的 Frobenius 级数解；
两根相等 ( $r_{1} = r_{2}$ r_1 = r_2)：第二个解必定包含对数项： $y_{2} (x) = y_{1} (x) \ln x + x^{r_{1}} \sum b_{n} x^{n}$ y_2(x) = y_1(x)ln x + x^r_1sum b_n x^n；
两根相差非零整数：较小根的解可能带有对数因子，需通过对指标参数求导法确定．

对数项的物理图景：奇异场源的本征体现

贝塞尔方程在 $x = 0$ x=0 处的指标方程两根相差整数或相等，由此解出包含 $\ln x$ ln x 奇异项的第二类贝塞尔函数（诺伊曼函数 $Y_{ν} (x)$ Y_nu(x)）．在无源自由介质中由于能量有限性要求通常舍去对数解；而在对称轴存在线电荷或线电流源时，对数解正是静电线势 $ϕ \propto \ln r$ phi propto ln r 的精确数学写照！

4. 核心物理思想：级数截断与量子化

在物理学中，求解勒让德方程或一维谐振子薛定谔方程时，级数递推关系往往给出无穷级数：

形式分析表明：当 $n \to \infty$ n -> inf 时，无穷级数的渐近行为类似于 $\frac{1}{1 - x^{2}}$ 1/1-x^2 或 $e^{x^{2}}$ e^x^2，使得解在物理边界处（如球面的南北极 $x = \pm 1$ x = +/- 1 或无穷远处 $x \to \infty$ x -> inf）发生 严重的发散；
物理现实（概率有限、波函数归一化、场强有限）强制要求解在全空间保持有界；
唯一的自救之道：递推关系分子中必须存在某一项恰好等于零，使得从 $n = l$ n = l 阶起的所有更高阶系数全部为零（级数在第 $l$ l 项发生截断，退化为有限项多项式）！

这一「为了边界有界而强制截断级数」的要求，直接决定了方程中分离常数只能取分立整数 $l = 0, 1, 2, \dots$ l = 0, 1, 2, dots！这就是为什么经典电磁场会出现分立多极矩，而量子微观世界会出现能级量子化的数学根本原理．

5. 学习衔接

下一节：研究二维和三维偏微分方程的空间与时间分离，进入三类经典物理偏微分方程；
成果应用：在勒让德多项式与贝塞尔函数中检验幂级数截断的具体成果．

本页面最近更新：2026/10/1 10:49:01，更新历史
发现错误？想一起完善？在 GitHub 上编辑此页！
本页面贡献者：Leafuke, Physics Learning Wiki Team
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供，附加条款亦可能应用

幂级数解法与Frobenius方法