微分方程的幂级数解法

用幂级数表示微分方程的解的方法。当方程无法用初等函数显式求解时（如变系数方程），可假设解在 $x_0$ 附近展开为幂级数。

方法步骤

1. 设幂级数解：假设解在 $x_0$ 邻域内可展开为

$$y = \sum_{n=0}^{\infty} a_n (x-x_0)^n$$

2. 逐项微分：代入原方程，得到 $y', y'', \dots$ 的级数表达式。
3. 比较系数：合并同次幂，由恒等式系数相等得到 $a_n$ 的递推关系。
4. 求通项：解递推关系确定各项系数，得级数形式的解。

示例

求 $y'' + xy = 0$ 在 $x=0$ 附近的级数解。

设 $y = \sum_{n=0}^{\infty} a_n x^n$，则 $y'' = \sum_{n=2}^{\infty} n(n-1)a_n x^{n-2}$。代入得：

$$\sum_{n=2}^{\infty} n(n-1)a_n x^{n-2} + \sum_{n=0}^{\infty} a_n x^{n+1} = 0$$

比较系数：

• $x^0$: $2a_2 = 0 \Rightarrow a_2 = 0$
• $x^1$: $6a_3 = 0 \Rightarrow a_3 = 0$
• $x^k\ (k\ge 2)$: $(k+2)(k+1)a_{k+2} + a_{k-1} = 0 \Rightarrow a_{k+2} = -\frac{a_{k-1}}{(k+2)(k+1)}$

适用与局限

• 适用于解析系数（$P(x), Q(x)$ 在 $x_0$ 处解析）的线性齐次微分方程。
• 解在 $[x_0-R, x_0+R]$ 内收敛，$R$ 为系数函数的收敛半径。
• 涉及正则奇点时需用推广的 Frobenius 方法（含非整数幂项）。