微分方程不是“见到导数就积分”,而是根据方程结构选择能保留下来的量、合适的变量和解空间。考试最重要的能力是:先分类,再动手。
先按方程分类
| 模块 | 要解决的问题 | 主要工具 |
|---|
| 一阶显式方程 | y′=f(x,y) 能否化成标准类型? | 变量分离、线性方程、Bernoulli、恰当方程、积分因子 |
| 一阶隐式方程 | 方程没有直接解出 y′ 怎么办? | 令 p=y′、Clairaut、参数法、可解出 x 或 y |
| 高阶线性方程 | 如何描述全部解,而不只找一个解? | Wronski 行列式、特征方程、待定系数、常数变易 |
| Laplace 变换 | 初值、分段输入和冲激怎样一次处理? | 导数变换、部分分式、卷积、移位 |
| 线性方程组 | 多个状态量联立变化怎么解? | 基解矩阵、矩阵指数、特征向量、Jordan 链 |
| 定性理论 | 不求显式解,能否保证解存在、唯一、可延拓? | 方向场、Picard 迭代、Lipschitz、比较定理、连续依赖 |
一阶方程识别顺序
拿到 y′=f(x,y),按下面顺序试:
- 能否分离变量:y′=X(x)Y(y);
- 是否齐次型:y′=F(y/x),令 v=y/x;
- 是否一阶线性:y′+p(x)y=q(x);
- 是否 Bernoulli:y′+py=qyn,令 z=y1−n;
- 写成 Mdx+Ndy=0,检查 My=Nx;
- 不恰当时看能否找到只依赖 x 或 y 的积分因子;
- 若没有解出 y′,令 p=y′ 检查 Clairaut、Lagrange 或参数结构。
分离变量时,除以 Y(y) 之前先检查 Y(y)=0 是否给出常值解;否则会把解丢掉。
一阶线性方程的积分因子
μ(x)=e∫p(x)dx
把左边变成 (μy)′:
y=e−∫pdx(C+∫q(x)e∫pdxdx).
恰当方程 Mdx+Ndy=0 对应一个势函数 F:
Fx=M,Fy=N,F(x,y)=C.
高阶线性方程的骨架
标准形式
y(n)+a1(x)y(n−1)+⋯+an(x)y=f(x).
齐次方程的解构成 n 维线性空间。找齐 n 个线性无关解 y1,…,yn 后,
yh=C1y1+⋯+Cnyn.
非齐次通解是
y=yh+yp.
常系数齐次方程用特征多项式。根 r=α±iβ、重数 m 对应
eαxxkcosβx,eαxxksinβx,k=0,…,m−1.
待定系数法只适用于多项式、指数、正弦余弦及其有限组合。试探式若与齐次解重合,要乘足够高次的 x 消除共振。一般右端则用常数变易法。
Laplace 变换
L{f(t)}=F(s)=∫0∞e−stf(t)dt.
导数公式自动带入初值:
L{y(n)}=snY−sn−1y(0)−⋯−y(n−1)(0).
因此流程是:对方程整体变换,代初值,代数求 Y(s),部分分式,再逆变换。分段输入优先写成 Heaviside 阶跃函数;卷积对应变换域乘法。
线性方程组
x′=A(t)x+f(t).
齐次组的基解矩阵 Φ 满足
Φ′=AΦ,detΦ=0.
非齐次组通解为
x(t)=Φ(t)[c+∫t0tΦ−1(τ)f(τ)dτ].
常系数时 Φ=eAt。若 A 可对角化,eAt=PeΛtP−1;不可对角化则用 Jordan 块,块中的幂零部分会产生 t,t2/2!,… 因子。
存在、唯一和定性理论
对初值问题
y′=f(x,y),y(x0)=y0,
f 连续通常保证局部存在;对 y 满足局部 Lipschitz 条件进一步保证局部唯一。fy 连续是常用的 Lipschitz 充分条件。
区别必须说清:
- 存在不等于唯一;
- 局部存在不等于解能延拓到全实轴;
- 方程右端在整个平面连续,也不自动排除有限时间爆炸;
- 初值和参数小变化导致解在共同区间内小变化,需要连续依赖定理的条件。
考场检查清单
- 除法、换元前检查是否丢掉常值解或特殊解;
- 通解应含正确数量的任意常数;
- 初值题最后代回初值,不要只写通解;
- 线性非齐次解始终分成齐次通解和一个特解;
- 矩阵解要检查基解矩阵是否可逆;
- 存在唯一性题要写区域和 Lipschitz 条件,不只写定理名称;
- 得到解后直接代回原方程,是发现符号错误最快的方法。