高阶非线性方程通常没有统一方法,但方程若缺少某个变量,就提示存在可利用的不变性;线性方程在普通函数难以求解时,还可以在一点附近直接构造级数解。
缺少 y:令 p=y′
若方程不显含 y,形如
F(x,y′,y′′,…)=0,
令 p(x)=y′(x),则
y′′=p′,y′′′=p′′,…
方程阶数降低一阶。先解 p,再积分 y=∫p,dx+C。
缺少 x:把 y′ 看成 y 的函数
若方程不显含 x,令
p(y)=dxdy.
链式法则给
y′′=dxdp=dydpdxdy=pdydp.
三阶导还要继续按 d/dx=p,d/dy 计算。解出 p(y) 后,再由 dy/dx=p(y) 分离变量。
除以 p 前要检查 p=0 对应的常值解。
已知一个解的线性降阶
二阶齐次线性方程已知解 y1 时,令 y=vy1 得第二个解
y2=y1∫y12e−∫pdxdx.
更高阶方程若已知一个非零解,也可因式分解相应微分算子或令 y=y1u 降低阶数。
常点附近的幂级数解
考虑
y′′+p(x)y′+q(x)y=0.
若 p,q 在 x0 附近解析,x0 是常点。设
y=n=0∑∞an(x−x0)n,
逐项求导并代入,统一成相同幂次,比较系数得到递推关系。a0,a1 通常是两个自由常数,对应二阶解空间的两个维度。
计算步骤:
- 写出 y,y′,y′′;
- 所有求和指标平移到 (x−x0)n;
- 处理单独出现的低次项;
- 比较每个 n 的系数;
- 用递推式生成若干项并识别熟悉函数;
- 给出收敛区间或至少说明受最近奇点限制。
正则奇点与 Frobenius 法
若 x0 不是常点,但
(x−x0)p(x),(x−x0)2q(x)
在 x0 附近解析,则 x0 是正则奇点。尝试
y=(x−x0)rn=0∑∞an(x−x0)n,a0=0.
最低次幂给出指标方程,先解 r;其余系数给递推。两个指标根差不是整数时,通常得到两个独立 Frobenius 解;重根或相差整数时,第二个解可能含 ln∣x−x0∣。
课件的幂级数部分强调:定理保证的是在相应解析条件下存在这种形式的解,不是所有奇点都能用普通幂级数处理。
常见错误
- 缺 x 的降阶中把 y′′ 错写成 dp/dy,漏掉因子 p;
- 除以 p 丢掉常值解;
- 级数换指标后起始下标不一致;
- 忘记单独处理 n=0,1 项;
- Frobenius 法没有先解指标方程;
- 写出形式级数却不说明适用点是常点还是正则奇点。