常微分方程

考点定位

  • 大纲位置:高等数学第四章,解方程"程序化"程度最高——识别类型 → 套固定流程 → 算准
  • 常见题型:一阶方程识别与求解、二阶常系数齐次/非齐次通解、可降阶方程、含变限积分的方程(求导脱限)、由几何/物理条件列方程
  • 考查频率:高。解答题几乎每年一道,且常与变限积分、几何应用综合
  • 定位:全部考法收敛为"类型识别表 + 特解设法表"两张表;综合题套路是"求导脱限 → 带初始条件的方程 → 回代验证"

卷种差异(卷种确认前按数一超集学,确认后可划走):可降阶方程数一数二大纲均含;伯努利方程、全微分方程、欧拉方程仅数一考。

核心讲解

基本概念:微分方程的阶、解、通解、特解

微分方程:含有未知函数 y(x)y(x) 及其导数(或微分)的等式。阶:方程中出现的最高阶导数的阶数(如 y′′+p(x)y′+q(x)y=f(x)y''+p(x)y'+q(x)y = f(x) 是二阶)。

解:代入方程后使之成为恒等式的函数。通解:含有独立任意常数、且任意常数的个数等于方程的阶数的解;特解:不含任意常数的解,通常由初始条件定出。初始条件与初值问题:nn 阶方程一般需要 nn 个条件 y(x0), y′(x0), …, y(n−1)(x0)y(x_0),\ y'(x_0),\ \ldots,\ y^{(n-1)}(x_0) 把通解中的常数逐一锁定,方程加初始条件合称初值问题。

线性方程:未知函数 yy 及其各阶导数都是一次、且系数只依赖 xx 的方程,一般形式为

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)

右端 f(x)≡0f(x)\equiv 0 时称齐次线性方程,否则称非齐次线性方程。

一词两义警报(命题人的措辞陷阱):本章"齐次"有两个互不相关的含义——① 线性方程语境下的齐次(右端为 0,如 y′+P(x)y=0y'+P(x)y=0);② 一阶方程 y′=φ(yx)y'=\varphi\left(\dfrac{y}{x}\right) 形态的齐次(方程各项关于 x,yx,y 的总次数相同)。看到"齐次"先判断语境属于哪一种,二者解法完全不同。

类型识别总表(一阶)

拿到一阶方程先化成标准形态(把 y′y' 解出来或写成 P dx+Q dy=0P\,dx + Q\,dy = 0),再对照:

形态特征类型解法入口
y′=f(x)g(y)y' = f(x)g(y)可分离变量∫dyg(y)=∫f(x) dx+C\int\dfrac{dy}{g(y)} = \int f(x)\,dx + C
y′=φ(yx)y' = \varphi\left(\dfrac{y}{x}\right)齐次令 u=yxu = \dfrac{y}{x}(即 y=uxy = ux,y′=u+xu′y' = u + xu')化为可分离
y′+P(x)y=Q(x)y' + P(x)y = Q(x)一阶线性通解公式(见下)
y′+P(x)y=Q(x)yny' + P(x)y = Q(x)y^n(n≠0,1n\ne 0,1)伯努利(仅数一)除以 yny^n,令 z=y1−nz = y^{1-n} 化为一阶线性

识别技巧:看各项关于 x,yx, y 的"总次数"是否相同(齐次的代数特征:x2+xyx2−y2\dfrac{x^2+xy}{x^2-y^2} 型每项次数一致)。

一阶线性方程与常数变易法推导

对 y′+P(x)y=Q(x)y' + P(x)y = Q(x),先解齐次方程 y′+Py=0y' + Py = 0(可分离):

y=Ce−∫P(x)dxy = Ce^{-\int P(x)dx}

常数变易:把常数 CC 升级为待定函数 C(x)C(x),设 y=C(x)e−∫P dxy = C(x)e^{-\int P\,dx},代回原方程:

C′(x)e−∫Pdx−P⋅C(x)e−∫Pdx+P⋅C(x)e−∫Pdx=Q(x)C'(x)e^{-\int Pdx} - P\cdot C(x)e^{-\int Pdx} + P\cdot C(x)e^{-\int Pdx} = Q(x)

中间两项相消,剩下 C′(x)=Q(x)e∫P dxC'(x) = Q(x)e^{\int P\,dx},积分得 C(x)=∫Qe∫Pdxdx+C1C(x) = \int Qe^{\int Pdx}dx + C_1。于是通解公式:

y=e−∫Pdx(∫Q(x)e∫Pdxdx+C)y = e^{-\int Pdx}\left(\int Q(x)e^{\int Pdx}dx + C\right)

使用公式前必须化成标准形态(y′y' 的系数为 1),否则 PP 取错全盘皆输。理解推导比背公式可靠——考场上现推 30 秒,比记错公式强。

伯努利方程(仅数一)

y′+P(x)y=Q(x)yny' + P(x)y = Q(x)y^n:两边除以 yny^n(y≠0y\ne 0),得 y−ny′+Py1−n=Qy^{-n}y' + Py^{1-n} = Q。注意 y1−ny^{1-n} 的导数恰是 (1−n)y−ny′(1-n)y^{-n}y',令 z=y1−nz = y^{1-n} 即化为一阶线性 z′+(1−n)Pz=(1−n)Qz' + (1-n)Pz = (1-n)Q。别忘检查 y≡0y\equiv 0 是否为解(除以 yny^n 时丢掉了它)。

可降阶的二阶方程(数一数二均考)

类型特征降阶手段
y(n)=f(x)y^{(n)} = f(x)右端只含 xx直接积分 nn 次(每次补一个 CC)
y′′=f(x,y′)y'' = f(x, y')不含 yy令 p(x)=y′p(x) = y',y′′=p′y'' = p',降为一阶方程解 pp 再积分
y′′=f(y,y′)y'' = f(y, y')不含 xx令 p=p(y)p = p(y),y′′=dpdx=dpdy⋅dydx=pdpdyy'' = \dfrac{dp}{dx} = \dfrac{dp}{dy}\cdot\dfrac{dy}{dx} = p\dfrac{dp}{dy},降为 pp 关于 yy 的一阶方程

第三类的 y′′=pdpdyy'' = p\dfrac{dp}{dy} 是本节最易错点——pp 被视为 yy 的函数,不是 xx 的。区分口诀:缺谁不代谁(缺 yy 则设 p=p(x)p = p(x);缺 xx 则设 p=p(y)p = p(y))。

二阶常系数齐次:特征方程法

y′′+py′+qy=0y'' + py' + qy = 0(p,qp,q 为常数),猜 y=erxy = e^{rx} 代入得特征方程 r2+pr+q=0r^2 + pr + q = 0。三种情形:

判别式特征根通解
Δ>0\Delta > 0两不等实根 r1≠r2r_1 \ne r_2y=C1er1x+C2er2xy = C_1e^{r_1x} + C_2e^{r_2x}
Δ=0\Delta = 0二重根 rry=(C1+C2x)erxy = (C_1 + C_2x)e^{rx}(重根补 xx)
Δ<0\Delta < 0共轭复根 α±βi\alpha\pm\beta iy=eαx(C1cos⁡βx+C2sin⁡βx)y = e^{\alpha x}(C_1\cos\beta x + C_2\sin\beta x)

结构定理背景:两个线性无关解的线性组合即通解。Δ<0\Delta<0 情形的实解由欧拉公式 e(α+βi)xe^{(\alpha+\beta i)x} 的实虚部拆出——记住结论即可,但要记得写 eαxe^{\alpha x} 这个因子。

二阶常系数非齐次:特解设法规律表

y′′+py′+qy=f(x)y'' + py' + qy = f(x),通解 = 齐次通解 + 非齐次特解。特解不用积分,用待定系数猜:

右端 f(x)f(x)特解设法kk 的取法
Pm(x)eλxP_m(x)e^{\lambda x}(多项式×指数)y∗=xkQm(x)eλxy^* = x^k Q_m(x)e^{\lambda x},QmQ_m 为系数待定的同次多项式λ\lambda 不是特征根 k=0k=0;单根 k=1k=1;二重根 k=2k=2
eαx[P(x)cos⁡βx+Q(x)sin⁡βx]e^{\alpha x}\left[P(x)\cos\beta x + Q(x)\sin\beta x\right]y∗=xkeαx[Rmcos⁡βx+Smsin⁡βx]y^* = x^k e^{\alpha x}\left[R_m\cos\beta x + S_m\sin\beta x\right],m=max⁡{deg⁡P,deg⁡Q}m = \max\{\deg P, \deg Q\}α±βi\alpha\pm\beta i 不是特征根 k=0k=0;是特征根 k=1k=1

三条铁律:

  1. 共振必补 xkx^k:λ\lambda 落在特征根上时,不补 xkx^k 则待定方程无解(猜出的项已被齐次解吸收)
  2. cos 与 sin 成对设:右端哪怕只有 cos⁡βx\cos\beta x,特解也必须同时含 cos⁡\cos 与 sin⁡\sin(求导会互相产生)
  3. 多项式要设满次:右端是 xx 就设 ax+bax+b,不能只设 axax

叠加原理:右端是几项之和,特解就分别设再相加。

综合题套路:微分方程 × 变限积分 × 几何应用

含变限积分的方程(最高频综合形态)

方程形如 f(x)=g(x)+∫0xf(t) dtf(x) = g(x) + \int_0^x f(t)\,dt(或 ∫0xf(t)dt\int_0^x f(t)dt 的更复杂组合):

  1. 求导脱限:两边对 xx 求导,∫0xf(t)dt\int_0^x f(t)dt 变成 f(x)f(x),得到微分方程
  2. 定初值:在原方程中令 x=0x = 0(积分下限),得 f(0)f(0)——初值藏在原方程里,不求白丢
  3. 解初值问题,重要的话最后代回原方程检验(求导可能放大解集)

变体:∫0xxf(t) dt\int_0^x xf(t)\,dt、∫0xf(x−t) dt\int_0^x f(x-t)\,dt 这类被积函数含 xx 的,先拆项提出 / 换元(u=x−tu = x-t)再求导——衔接一元函数积分学变限求导节。

几何应用列方程

翻译表:切线斜率 →y′\to y';曲线在 (x,y)(x,y) 处切线纵截距/法线条件 →\to 关于 x,y,y′x,y,y' 的方程;"过点 (x0,y0)(x_0,y_0)" →\to 初始条件。列完照旧解方程。物理应用(功、压力、质心,仅数一)同理把物理定律翻译成 y′y' 与 x,yx,y 的关系。

例题

例题 1(模拟题)一阶线性:公式或凑微分

解方程 y′+ycos⁡x=e−sin⁡xy' + y\cos x = e^{-\sin x}。

解析:标准形态已就位,P=cos⁡xP = \cos x,Q=e−sin⁡xQ = e^{-\sin x}。观察 e−sin⁡xe^{-\sin x} 恰是 esin⁡xe^{\sin x} 的倒数,用凑微分更快:两边乘 esin⁡xe^{\sin x}:

esin⁡xy′+esin⁡xcos⁡x y=1  ⟹  (yesin⁡x)′=1e^{\sin x}y' + e^{\sin x}\cos x\, y = 1 \implies (y e^{\sin x})' = 1

积分得 yesin⁡x=x+Cy e^{\sin x} = x + C,即 y=(x+C)e−sin⁡xy = (x+C)e^{-\sin x}。

两种殊途:这就是通解公式里"积分因子 e∫Pdxe^{\int Pdx}"的手工版——公式忘了现场能重推,这是底线保障。

例题 2(模拟题)齐次方程

解 (x2+y2) dx−xy dy=0(x^2+y^2)\,dx - xy\,dy = 0。

解析:化为 y′=x2+y2xy=1u+uy' = \dfrac{x^2+y^2}{xy} = \dfrac{1}{u} + u 型(分子分母同除 x2x^2 即见 yx\frac{y}{x} 结构),令 y=uxy = ux:

u+xdudx=1u+u  ⟹  xdudx=1u  ⟹  u du=dxxu + x\frac{du}{dx} = \frac{1}{u} + u \implies x\frac{du}{dx} = \frac{1}{u} \implies u\,du = \frac{dx}{x}

积分 u22=ln⁡∣x∣+C1\dfrac{u^2}{2} = \ln|x| + C_1,回代 u=yxu = \dfrac{y}{x}:

y2=2x2ln⁡∣x∣+Cx2y^2 = 2x^2\ln|x| + Cx^2

坑:回代后结果可整理成隐式通解,不必强求显式;x=0x = 0 处的奇性由 ln⁡∣x∣\ln|x| 自然体现。

例题 3(模拟题)二阶常系数非齐次:多项式×指数共振

求 y′′−3y′+2y=xexy'' - 3y' + 2y = xe^x 的通解。

解析:

齐次部分:r2−3r+2=0  ⟹  r1=1, r2=2r^2 - 3r + 2 = 0 \implies r_1 = 1,\ r_2 = 2,齐次通解 Y=C1ex+C2e2xY = C_1e^x + C_2e^{2x}。

特解:右端 λ=1\lambda = 1 恰是单特征根,k=1k=1;Pm(x)=xP_m(x) = x(一次),设 y∗=(ax2+bx)exy^* = (ax^2+bx)e^x。代入时用多项式整体代换(exe^x 因子提出)最省力:记 P=ax2+bxP = ax^2+bx,则 y′′−3y′+2y=ex(P′′−3P′+2P)y'' - 3y' + 2y = e^x(P'' - 3P' + 2P),其中 P′′=2aP'' = 2a,P′=2ax+bP' = 2ax+b:

2a−3(2ax+b)+2(ax2+bx)=2ax2+(2b−6a)x+(2a−3b)=!x2a - 3(2ax+b) + 2(ax^2+bx) = 2ax^2 + (2b-6a)x + (2a-3b) \stackrel{!}{=} x

比较系数:2a=0⇒a=02a = 0 \Rightarrow a = 0;2b=1⇒b=122b = 1 \Rightarrow b = \dfrac{1}{2}。故 y∗=x2exy^* = \dfrac{x}{2}e^x,通解:

y=C1ex+C2e2x+x2exy = C_1e^x + C_2e^{2x} + \frac{x}{2}e^x

技巧复述:exe^x 型右端代入时保持 exe^x 不动、只算多项式组合 P′′−3P′+2PP''-3P'+2P,比三次求导快得多且不易错。

例题 4(模拟题)可降阶:缺 x 型

解 yy′′=y′2yy'' = y'^2(y>0y > 0)。

解析:方程不含 xx,令 p=p(y)p = p(y),y′′=pdpdyy'' = p\dfrac{dp}{dy}:

y pdpdy=p2  ⟹  dpp=dyy  ⟹  p=C1y  ⟹  y′=C1yy\,p\frac{dp}{dy} = p^2 \implies \frac{dp}{p} = \frac{dy}{y} \implies p = C_1y \implies y' = C_1y

再积分 dyy=C1dx  ⟹  y=C2eC1x\dfrac{dy}{y} = C_1dx \implies y = C_2e^{C_1x}。

坑:第一步除以 pp 时若 p=0p = 0(即 y≡y\equiv 常数)也是解,被 C1=0C_1 = 0 情形吸收,故通解完整。

例题 5(模拟题)综合:变限积分方程

设 f(x)f(x) 连续,且满足 f(x)=ex+∫0xf(t) dtf(x) = e^x + \int_0^x f(t)\,dt,求 f(x)f(x)。

解析:

求导脱限:f′(x)=ex+f(x)f'(x) = e^x + f(x),即 f′−f=exf' - f = e^x(一阶线性)。

定初值:原方程令 x=0x = 0:f(0)=e0+0=1f(0) = e^0 + 0 = 1。

解初值问题:积分因子 e−xe^{-x}:(fe−x)′=1  ⟹  fe−x=x+C(fe^{-x})' = 1 \implies fe^{-x} = x + C,由 f(0)=1f(0) = 1 得 C=1C = 1:

f(x)=(x+1)exf(x) = (x+1)e^x

检验:∫0x(t+1)et dt=[tet]0x=xex\int_0^x (t+1)e^t\,dt = \left[te^t\right]_0^x = xe^x,代回右边 ex+xex=(x+1)ex=f(x)e^x + xe^x = (x+1)e^x = f(x),成立。

套路总结:脱限 → 取初值 → 解 → 回代检验,四步缺一不可;初值必须从原方程取(求导后的方程已丢失这一信息)。

常见坑

  1. 非齐次特解不查 λ\lambda 是否特征根:共振不补 xkx^k,待定系数方程直接矛盾(例题 3 若不补 xx 将解不出 a,ba,b)
  2. 右端只有 cos⁡\cos 就只设 cos⁡\cos:特解必须 cos⁡\cos、sin⁡\sin 成对。对照:y′′+y=cos⁡xy''+y=\cos x 的特解是 −x2sin⁡x-\dfrac{x}{2}\sin x(共振且由 sin⁡\sin 承担),设 acos⁡xa\cos x 毫无胜算
  3. 一阶线性不化标准形态就取 P,QP,Q,符号错全盘错
  4. 可降阶"缺 x 型"把 pp 当 xx 的函数:y′′=f(y,y′)y''=f(y,y') 必须 y′′=pdpdyy'' = p\dfrac{dp}{dy}
  5. 变限积分方程求导后忘记取初值,或初值不从原方程取
  6. Δ<0\Delta < 0 忘写 eαxe^{\alpha x} 因子:y′′+4y=0y''+4y=0 的解是 C1cos⁡2x+C2sin⁡2xC_1\cos 2x + C_2\sin 2x(α=0\alpha=0 隐去),但 y′′+2y′+5y=0y''+2y'+5y=0 的解必须带 e−xe^{-x}
  7. 通解常数个数不足:nn 阶方程 nn 个独立常数;可降阶每积一次分补一个
  8. 伯努利/欧拉方程当公共考点:二者仅数一考(卷种确认前按超集学)
  9. 特解多项式没设满次:右端 xx 要设 ax+bax+b,常数项不能省

自测

  1. 解 y′=ex−yy' = e^{x-y}
  2. 解 (1+x)y′+y=ex(1+x)y' + y = e^x
  3. 求 y′′−2y′+y=exy'' - 2y' + y = e^x 的通解(注意共振)
  4. 求 y′′+4y=3cos⁡2xy'' + 4y = 3\cos 2x 的通解
  5. 曲线过点 (0,1)(0,1),其上任一点 (x,y)(x,y) 处切线斜率为 2x+y2x + y,求曲线方程

自测答案与简析

  1. 可分离:ey dy=ex dx  ⟹  ey=ex+Ce^y\,dy = e^x\,dx \implies e^y = e^x + C,即 y=ln⁡(ex+C)y = \ln(e^x + C)
  2. 化标准形 y′+y1+x=ex1+xy' + \dfrac{y}{1+x} = \dfrac{e^x}{1+x},积分因子 1+x1+x:[(1+x)y]′=ex[(1+x)y]' = e^x,故 (1+x)y=ex+C(1+x)y = e^x + C
  3. Δ=0\Delta = 0,r=1r=1 二重根,齐次通解 (C1+C2x)ex(C_1+C_2x)e^x;右端 λ=1\lambda = 1 是二重根,k=2k=2:设 y∗=ax2exy^* = ax^2e^x,代入得 2aex=ex2ae^x = e^x,a=12a = \dfrac{1}{2}。通解 y=(C1+C2x)ex+x22exy = (C_1+C_2x)e^x + \dfrac{x^2}{2}e^x
  4. r=±2ir = \pm 2i,齐次通解 C1cos⁡2x+C2sin⁡2xC_1\cos 2x + C_2\sin 2x;±2i\pm 2i 是特征根,k=1k=1:设 y∗=x(acos⁡2x+bsin⁡2x)y^* = x(a\cos 2x + b\sin 2x),代入得 y′′+4y=−4asin⁡2x+4bcos⁡2x=3cos⁡2xy''+4y = -4a\sin 2x + 4b\cos 2x = 3\cos 2x,故 a=0a = 0、b=34b = \dfrac{3}{4}。通解 y=C1cos⁡2x+C2sin⁡2x+3x4sin⁡2xy = C_1\cos 2x + C_2\sin 2x + \dfrac{3x}{4}\sin 2x
  5. 列方程 y′=2x+yy' = 2x + y 且 y(0)=1y(0) = 1,即 y′−y=2xy' - y = 2x。积分因子 e−xe^{-x}:(ye−x)′=2xe−x(ye^{-x})' = 2xe^{-x},积分得 ye−x=−2xe−x−2e−x+Cye^{-x} = -2xe^{-x} - 2e^{-x} + C,即 y=Cex−2x−2y = Ce^x - 2x - 2;代入 y(0)=1y(0)=1 得 C=3C = 3:y=3ex−2x−2y = 3e^x - 2x - 2

本章关联