二阶常系数非齐次线性微分方程的解
(防止忘了记在这里)二阶常系数线性微分方程的通解由对应的齐次方程的通解+该非齐次方程的特解构成,所以关键在于求此类方程的特解。
整理自张宇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,2则k=0k=0k=0(不是特征根)
- α=λ1\alpha = \lambda_{1}α=λ1或α=λ2\alpha = \lambda_{2}α=λ2则k=1k=1k=1(是单特征根)
- α=λ1=λ2\alpha = \lambda_{1}=\lambda_{2}α=λ1=λ2则k=2k=2k=2(是二重特征根)
Qn(x)Q_n(x)Qn(x)为x的nnn次多项式,意思是2x2x2x这种写为Ax+BAx+BAx+B,而x3−1x^3-1x3−1这种写为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^xaex−2aex+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λ2−2λ+5=0
由于Δ=b2−4a=−16<0\Delta=b^2-4a=-16<0Δ=b2−4a=−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(C1cos2x+C2sin2x)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(C1cos2x+C2sin2x)+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^22x2和222,则l=2l=2l=2,分别写为a1x2+b1x+c1a_1x^2+b_1x+c_1a1x2+b1x+c1和a2x2+b2x+c2a_2x^2+b_2x+c_2a2x2+b2x+c2
例题:
y′′−2y′+5y=−excos2xy'' - 2y' + 5y = -e^x\cos 2xy′′−2y′+5y=−excos2x
首先右侧不是多项式形式,有点像复根形式,所以得进行变换,得到下式
y′′−2y′+5y=−ex[(−1)cos2x+0⋅sin2x]y'' - 2y' + 5y = -e^x[(-1)\cos 2x+0·\sin2x]y′′−2y′+5y=−ex[(−1)cos2x+0⋅sin2x]
此处α=1\alpha = 1α=1,l=1l=1l=1,所以QQQ是一次多项式,将特解设为
y∗=ex[Acos2x+Bsin2x]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=−excos2xy'' - 2y' + 5y = -e^x\cos 2xy′′−2y′+5y=−excos2x可求出AAA和BBB
其他
注意关系,做题时也可能逆着推回去。尤其是在高阶微分方程的情况中,利用特解反推特征根,一步一步反推回微分方程的可能形式。
更多推荐
所有评论(0)