第 11 讲:降阶法与幂级数解

高阶非线性方程通常没有统一方法,但方程若缺少某个变量,就提示存在可利用的不变性;线性方程在普通函数难以求解时,还可以在一点附近直接构造级数解。

缺少 yy:令 p=y′p=y'

若方程不显含 yy,形如

F(x,y′,y′′,… )=0,F(x,y',y'',\dots)=0,

令 p(x)=y′(x)p(x)=y'(x),则

y′′=p′,y′′′=p′′,…y''=p',\quad y'''=p'',\dots

方程阶数降低一阶。先解 pp,再积分 y=∫p,dx+Cy=\int p,dx+C。

缺少 xx:把 y′y' 看成 yy 的函数

若方程不显含 xx,令

p(y)=dydx.p(y)=\frac{dy}{dx}.

链式法则给

y′′=dpdx=dpdydydx=pdpdy.y''=\frac{dp}{dx} =\frac{dp}{dy}\frac{dy}{dx} =p\frac{dp}{dy}.

三阶导还要继续按 d/dx=p,d/dyd/dx=p,d/dy 计算。解出 p(y)p(y) 后,再由 dy/dx=p(y)dy/dx=p(y) 分离变量。

除以 pp 前要检查 p=0p=0 对应的常值解。

已知一个解的线性降阶

二阶齐次线性方程已知解 y1y_1 时,令 y=vy1y=v y_1 得第二个解

y2=y1∫e−∫pdxy12dx.y_2=y_1\int\frac{e^{-\int p dx}}{y_1^2}dx.

更高阶方程若已知一个非零解,也可因式分解相应微分算子或令 y=y1uy=y_1u 降低阶数。

常点附近的幂级数解

考虑

y′′+p(x)y′+q(x)y=0.y''+p(x)y'+q(x)y=0.

若 p,qp,q 在 x0x_0 附近解析,x0x_0 是常点。设

y=∑n=0∞an(x−x0)n,y=\sum_{n=0}^{\infty}a_n(x-x_0)^n,

逐项求导并代入,统一成相同幂次,比较系数得到递推关系。a0,a1a_0,a_1 通常是两个自由常数,对应二阶解空间的两个维度。

计算步骤:

  1. 写出 y,y′,y′′y,y',y'';
  2. 所有求和指标平移到 (x−x0)n(x-x_0)^n;
  3. 处理单独出现的低次项;
  4. 比较每个 nn 的系数;
  5. 用递推式生成若干项并识别熟悉函数;
  6. 给出收敛区间或至少说明受最近奇点限制。

正则奇点与 Frobenius 法

若 x0x_0 不是常点,但

(x−x0)p(x),(x−x0)2q(x)(x-x_0)p(x),\qquad (x-x_0)^2q(x)

在 x0x_0 附近解析,则 x0x_0 是正则奇点。尝试

y=(x−x0)r∑n=0∞an(x−x0)n,a0≠0.y=(x-x_0)^r \sum_{n=0}^{\infty}a_n(x-x_0)^n,\qquad a_0\ne0.

最低次幂给出指标方程,先解 rr;其余系数给递推。两个指标根差不是整数时,通常得到两个独立 Frobenius 解;重根或相差整数时,第二个解可能含 ln⁡∣x−x0∣\ln|x-x_0|。

课件的幂级数部分强调:定理保证的是在相应解析条件下存在这种形式的解,不是所有奇点都能用普通幂级数处理。

常见错误

  • 缺 xx 的降阶中把 y′′y'' 错写成 dp/dydp/dy,漏掉因子 pp;
  • 除以 pp 丢掉常值解;
  • 级数换指标后起始下标不一致;
  • 忘记单独处理 n=0,1n=0,1 项;
  • Frobenius 法没有先解指标方程;
  • 写出形式级数却不说明适用点是常点还是正则奇点。

评论