速成 · 微分方程

微分方程不是“见到导数就积分”,而是根据方程结构选择能保留下来的量、合适的变量和解空间。考试最重要的能力是:先分类,再动手。

先按方程分类

模块要解决的问题主要工具
一阶显式方程y′=f(x,y)y'=f(x,y) 能否化成标准类型?变量分离、线性方程、Bernoulli、恰当方程、积分因子
一阶隐式方程方程没有直接解出 y′y' 怎么办?令 p=y′p=y'、Clairaut、参数法、可解出 xx 或 yy
高阶线性方程如何描述全部解,而不只找一个解?Wronski 行列式、特征方程、待定系数、常数变易
Laplace 变换初值、分段输入和冲激怎样一次处理?导数变换、部分分式、卷积、移位
线性方程组多个状态量联立变化怎么解?基解矩阵、矩阵指数、特征向量、Jordan 链
定性理论不求显式解,能否保证解存在、唯一、可延拓?方向场、Picard 迭代、Lipschitz、比较定理、连续依赖

一阶方程识别顺序

拿到 y′=f(x,y)y'=f(x,y),按下面顺序试:

  1. 能否分离变量:y′=X(x)Y(y)y'=X(x)Y(y);
  2. 是否齐次型:y′=F(y/x)y'=F(y/x),令 v=y/xv=y/x;
  3. 是否一阶线性:y′+p(x)y=q(x)y'+p(x)y=q(x);
  4. 是否 Bernoulli:y′+py=qyny'+py=qy^n,令 z=y1−nz=y^{1-n};
  5. 写成 Mdx+Ndy=0Mdx+Ndy=0,检查 My=NxM_y=N_x;
  6. 不恰当时看能否找到只依赖 xx 或 yy 的积分因子;
  7. 若没有解出 y′y',令 p=y′p=y' 检查 Clairaut、Lagrange 或参数结构。

分离变量时,除以 Y(y)Y(y) 之前先检查 Y(y)=0Y(y)=0 是否给出常值解;否则会把解丢掉。

一阶线性方程的积分因子

μ(x)=e∫p(x)dx\mu(x)=e^{\int p(x)dx}

把左边变成 (μy)′(\mu y)':

y=e−∫pdx(C+∫q(x)e∫pdxdx).y=e^{-\int p dx} \left(C+\int q(x)e^{\int p dx}dx\right).

恰当方程 Mdx+Ndy=0Mdx+Ndy=0 对应一个势函数 FF:

Fx=M,Fy=N,F(x,y)=C.F_x=M,\qquad F_y=N,\qquad F(x,y)=C.

高阶线性方程的骨架

标准形式

y(n)+a1(x)y(n−1)+⋯+an(x)y=f(x).y^{(n)}+a_1(x)y^{(n-1)}+\cdots+a_n(x)y=f(x).

齐次方程的解构成 nn 维线性空间。找齐 nn 个线性无关解 y1,…,yny_1,\dots,y_n 后,

yh=C1y1+⋯+Cnyn.y_h=C_1y_1+\cdots+C_ny_n.

非齐次通解是

y=yh+yp.y=y_h+y_p.

常系数齐次方程用特征多项式。根 r=α±iβr=\alpha\pm i\beta、重数 mm 对应

eαxxkcos⁡βx,eαxxksin⁡βx,k=0,…,m−1.e^{\alpha x}x^k\cos\beta x,\quad e^{\alpha x}x^k\sin\beta x,\qquad k=0,\dots,m-1.

待定系数法只适用于多项式、指数、正弦余弦及其有限组合。试探式若与齐次解重合,要乘足够高次的 xx 消除共振。一般右端则用常数变易法。

Laplace 变换

L{f(t)}=F(s)=∫0∞e−stf(t)dt.\mathcal L\{f(t)\}=F(s)=\int_0^\infty e^{-st}f(t)dt.

导数公式自动带入初值:

L{y(n)}=snY−sn−1y(0)−⋯−y(n−1)(0).\mathcal L\{y^{(n)}\} =s^nY-s^{n-1}y(0)-\cdots-y^{(n-1)}(0).

因此流程是:对方程整体变换,代初值,代数求 Y(s)Y(s),部分分式,再逆变换。分段输入优先写成 Heaviside 阶跃函数;卷积对应变换域乘法。

线性方程组

x′=A(t)x+f(t).\mathbf x'=A(t)\mathbf x+\mathbf f(t).

齐次组的基解矩阵 Φ\Phi 满足

Φ′=AΦ,det⁡Φ≠0.\Phi'=A\Phi,\qquad \det\Phi\ne0.

非齐次组通解为

x(t)=Φ(t)[c+∫t0tΦ−1(τ)f(τ)dτ].\mathbf x(t)=\Phi(t) \left[ \mathbf c+\int_{t_0}^{t}\Phi^{-1}(\tau)\mathbf f(\tau)d\tau \right].

常系数时 Φ=eAt\Phi=e^{At}。若 AA 可对角化,eAt=PeΛtP−1e^{At}=Pe^{\Lambda t}P^{-1};不可对角化则用 Jordan 块,块中的幂零部分会产生 t,t2/2!,…t,t^2/2!,\dots 因子。

存在、唯一和定性理论

对初值问题

y′=f(x,y),y(x0)=y0,y'=f(x,y),\qquad y(x_0)=y_0,

ff 连续通常保证局部存在;对 yy 满足局部 Lipschitz 条件进一步保证局部唯一。fyf_y 连续是常用的 Lipschitz 充分条件。

区别必须说清:

  • 存在不等于唯一;
  • 局部存在不等于解能延拓到全实轴;
  • 方程右端在整个平面连续,也不自动排除有限时间爆炸;
  • 初值和参数小变化导致解在共同区间内小变化,需要连续依赖定理的条件。

考场检查清单

  • 除法、换元前检查是否丢掉常值解或特殊解;
  • 通解应含正确数量的任意常数;
  • 初值题最后代回初值,不要只写通解;
  • 线性非齐次解始终分成齐次通解和一个特解;
  • 矩阵解要检查基解矩阵是否可逆;
  • 存在唯一性题要写区域和 Lipschitz 条件,不只写定理名称;
  • 得到解后直接代回原方程,是发现符号错误最快的方法。

评论