整理自张宇30讲

前置:解的结构

对于二阶常系数线性微分方程y′′+py′+qy=f(x)y''+py'+qy=f(x)y′′+py+qy=f(x)

  • f(x)≡0f(x) \equiv 0f(x)0,则该方程为齐次方程
  • f(x)f(x)f(x)不恒等于0,则该方程为非齐次方程

齐次方程通解:由两个线性无关y1(x)y_1(x)y1(x)y2(x)y_2(x)y2(x)组合成齐次方程的通解y(x)y(x)y(x)

y(x)=C1y1(x)+C2y2(x)y(x)=C_1y_1(x)+C_2y_2(x)y(x)=C1y1(x)+C2y2(x)


特解记为y∗(x)y^*(x)y(x)

齐次方程通解=齐次方程通解+齐次方程特解

齐次方程通解为y(x)+y∗(x)y(x)+y^*(x)y(x)+y(x)


当有y′′+py′+qy=P1(x)y''+py'+qy=P_1(x)y′′+py+qy=P1(x)y′′+py′+qy=P2(x)y''+py'+qy=P_2(x)y′′+py+qy=P2(x)时,且y1∗(x)y_1^*(x)y1(x)y2∗(x)y_2^*(x)y2(x)分别是两式的解,则y′′+py′+qy=P1(x)+P2(x)y''+py'+qy=P_1(x)+P_2(x)y′′+py+qy=P1(x)+P2(x)的解为y1∗(x)+y2∗(x)y_1^*(x)+y_2^*(x)y1(x)+y2(x)

二阶常系数非齐次线性微分方程的特解

二阶常系数非齐次线性微分方程求通解的整体思路都如下:

  • 先设特解
  • 求出特解
  • 和齐次方程通解组合成为非齐次方程的通解

其中特解按自由项f(x)f(x)f(x)分为两种情况:

  • 自由项f(x)=Pn(x)eαxf(x)=P_n(x)e^{\alpha x}f(x)=Pn(x)eαx
  • 自由项f(x)=eαx[Pm(x)cos⁡βx+Pn(x)sin⁡βx]f(x)=e^{\alpha x}[P_m(x) \cos \beta x + P_n(x) \sin \beta x]f(x)=eαx[Pm(x)cosβx+Pn(x)sinβx]

确定xkx^kxk:一看(自由项中的α±βi\alpha \pm \beta iα±βi)二算(齐次方程中的特征根)三比较(两者是否相等)

自由项f(x)=Pn(x)eαxf(x)=P_n(x)e^{\alpha x}f(x)=Pn(x)eαx

当自由项为f(x)=Pn(x)eαxf(x)=P_n(x)e^{\alpha x}f(x)=Pn(x)eαx,将特解设为eαxQn(x)xke^{\alpha x}Q_n(x)x^keαxQn(x)xk

其中:

  • eαxe^{\alpha x}eαx照抄
  • Qn(x)Q_n(x)Qn(x)为x的nnn次多项式
  • kkk的取值依据a+bia+bia+bi是否是特征根
    • α≠λ1,2\alpha \neq \lambda_{1,2}α=λ1,2k=0k=0k=0(不是特征根)
    • α=λ1\alpha = \lambda_{1}α=λ1α=λ2\alpha = \lambda_{2}α=λ2k=1k=1k=1(是单特征根)
    • α=λ1=λ2\alpha = \lambda_{1}=\lambda_{2}α=λ1=λ2k=2k=2k=2(是二重特征根)

Qn(x)Q_n(x)Qn(x)为x的nnn次多项式,意思是2x2x2x这种写为Ax+BAx+BAx+B,而x3−1x^3-1x31这种写为ax3+bx2+cax^3+bx^2+cax3+bx2+c


【例】求y′′−2y′+5y=exy''-2y'+5y=e^xy′′2y+5y=ex的通解

自由项符合f(x)=Pn(x)eαxf(x)=P_n(x)e^{\alpha x}f(x)=Pn(x)eαx形式,故特解要设为eαxQn(x)xke^{\alpha x}Q_n(x)x^keαxQn(x)xk形式。

由题意知α=1\alpha = 1α=1,而特征根λ1,2=1±2i\lambda_{1,2}= 1 \pm 2iλ1,2=1±2i,故kkk写0。

因为Qn(x)=1Q_n(x)=1Qn(x)=1是0次多项式,所以写成aaa就行,故特解y∗=aexy^*=ae^xy=aex

y∗y^*y带入y′′−2y′+5y=exy''-2y'+5y=e^xy′′2y+5y=ex,即aex−2aex+5aex=exae^x-2ae^x+5ae^x=e^xaex2aex+5aex=ex,可解出a=14a=\frac14a=41

故特解为y∗=14exy^*=\frac14e^xy=41ex

但是题目的式子是非齐次的,其解的结构构成是:齐次方程通解+齐次方程特解

所以先得求出y′′−2y′+5y=0y''-2y'+5y=0y′′2y+5y=0这个齐次方程的通解y(x)y(x)y(x),再和特解y∗y^*y组合即可,y(x)+y∗y(x)+y^*y(x)+y是最终答案

关于y′′−2y′+5y=0y''-2y'+5y=0y′′2y+5y=0,带入eλxe^{\lambda x}eλx得到特征方程:

λ2−2λ+5=0\lambda^2-2\lambda+5=0λ22λ+5=0

由于Δ=b2−4a=−16<0\Delta=b^2-4a=-16<0Δ=b24a=16<0,故有一对共轭复根,即

λ1,2=2±Δi2=1±2i\lambda_{1,2}=\frac{2 \pm \sqrt{\Delta}i}{2}=1\pm 2iλ1,2=22±Δ i=1±2i
对于共轭复根α±βi\alpha \pm \beta iα±βi可知通解为y=eαx(C1cos⁡βx+C2sin⁡βx)y=e^{\alpha x}(C_1\cos \beta x + C_2 \sin \beta x)y=eαx(C1cosβx+C2sinβx),具体到本例则通解为:

y=ex(C1cos⁡2x+C2sin⁡2x)y=e^x(C_1\cos 2 x + C_2 \sin 2 x)y=ex(C1cos2x+C2sin2x)

则该非齐次方程的通解由y(x)+y∗y(x)+y^*y(x)+y可得如下:

ex(C1cos⁡2x+C2sin⁡2x)+14exe^x(C_1\cos 2 x + C_2 \sin 2 x)+\frac14e^xex(C1cos2x+C2sin2x)+41ex

自由项f(x)=eαx[Pm(x)cos⁡βx+Pn(x)sin⁡βx]f(x)=e^{\alpha x}[P_m(x) \cos \beta x + P_n(x) \sin \beta x]f(x)=eαx[Pm(x)cosβx+Pn(x)sinβx]

当自由项f(x)=eαx[Pm(x)cos⁡βx+Pn(x)sin⁡βx]f(x)=e^{\alpha x}[P_m(x) \cos \beta x + P_n(x) \sin \beta x]f(x)=eαx[Pm(x)cosβx+Pn(x)sinβx]时,特解应当设为

y∗=eαx[Ql(1)(x)cos⁡βx+Ql(2)(x)sin⁡βx]xky^*=e^{\alpha x}[Q_l^{(1)}(x) \cos \beta x + Q_l^{(2)}(x)\sin \beta x]x^ky=eαx[Ql(1)(x)cosβx+Ql(2)(x)sinβx]xk

其中:

  • eαxe^{\alpha x}eαx照抄
  • l=max{m,n}l=max\{m, n\}l=max{m,n}Ql(1)(x)Q_l^{(1)}(x)Ql(1)(x)Ql(2)(x)Q_l^{(2)}(x)Ql(2)(x)分别为xxx的两个不同的lll次多项式
  • kkk的取值依据a+bia+bia+bi是否是特征根
    • 不是则k=0k=0k=0
    • 是则k=1k=1k=1

【注】上面提到的lll次多项式,如2x22x^22x2222,则l=2l=2l=2,分别写为a1x2+b1x+c1a_1x^2+b_1x+c_1a1x2+b1x+c1a2x2+b2x+c2a_2x^2+b_2x+c_2a2x2+b2x+c2


例题:

y′′−2y′+5y=−excos⁡2xy'' - 2y' + 5y = -e^x\cos 2xy′′2y+5y=excos2x

首先右侧不是多项式形式,有点像复根形式,所以得进行变换,得到下式

y′′−2y′+5y=−ex[(−1)cos⁡2x+0⋅sin⁡2x]y'' - 2y' + 5y = -e^x[(-1)\cos 2x+0·\sin2x]y′′2y+5y=ex[(1)cos2x+0sin2x]

此处α=1\alpha = 1α=1l=1l=1l=1,所以QQQ是一次多项式,将特解设为

y∗=ex[Acos⁡2x+Bsin⁡2x]x1y*=e^x[A\cos 2x+B\sin2x]x^1y=ex[Acos2x+Bsin2x]x1

xkx^kxk处是由于齐次方程$y′′−2y′+5y=0y'' - 2y' + 5y =0y′′2y+5y=0特征根为1±2i1 \pm 2 i1±2i与上面式子里的α±βi\alpha \pm \beta iα±βi带入值相等,所以k=1k=1k=1

y∗y^*y带入y′′−2y′+5y=−excos⁡2xy'' - 2y' + 5y = -e^x\cos 2xy′′2y+5y=excos2x可求出AAABBB

其他

具有
解出
微分方程
特征方程
特征根

注意关系,做题时也可能逆着推回去。尤其是在高阶微分方程的情况中,利用特解反推特征根,一步一步反推回微分方程的可能形式。

Logo

有“AI”的1024 = 2048,欢迎大家加入2048 AI社区

更多推荐