LaTeX在线编辑器:Equation Editor
[toc]
一、定义
二阶常系数线性齐次微分方程:
y′′(x)+py′(x)+qy(x)=0
二阶常系数线性非齐次微分方程:
y′′(x)+py′(x)+qy(x)=g(x)
二阶常系数线性非齐次微分方程的通解:
y(x)=齐次通解+非齐次特解=y0(x)+y∗(x)
二、齐次通解
特征方程为
r2+pr+q=0
根据特征方程的根r1,r2的情况,设通解为
y0(x)=⎩⎨⎧C1er1x+C2er2x,(C1+C2x)er1x,eαx[C1cos(βx)+C2sin(βx)],r1=r2r1=r2r=α±iβ
三、非齐次特解
1. α不是特征根
| g(x)的形式 |
特解y∗(x)的形式 |
| a |
A |
| ax+b |
Ax+B |
| ax2+bx+c |
Ax2+Bx+C |
| aeαx |
Aeαx |
| (ax+b)eαx |
(Ax+B)eαx |
| (ax2+bx+c)eαx |
(Ax2+Bx+C)eαx |
2. α是特征单重根
(1)α=0时
| g(x)的形式 |
特解y∗(x)的形式 |
| aeαx |
xAeαx |
| (ax+b)eαx |
x(Ax+B)eαx |
| (ax2+bx+c)eαx |
x(Ax2+Bx+C)eαx |
(2)α=0时
| g(x)的形式 |
特解y∗(x)的形式 |
| a |
xA |
| ax+b |
x(Ax+B) |
| ax2+bx+c |
x(Ax2+Bx+C) |
3. α是特征二重根
(1)α=0时
| g(x)的形式 |
特解y∗(x)的形式 |
| aeαx |
x2Aeαx |
| (ax+b)eαx |
x2(Ax+B)eαx |
| (ax2+bx+c)eαx |
x2(Ax2+Bx+C)eαx |
(2)α=0时
| g(x)的形式 |
特解y∗(x)的形式 |
| a |
x2A |
| ax+b |
x2(Ax+B) |
| ax2+bx+c |
x2(Ax2+Bx+C) |
4. α+iβ不是特征根
| g(x)的形式 |
特解y∗(x)的形式 |
| eαx[acos(βx)+bsin(βx)] |
eαx[Acos(βx)+Bsin(βx)] |
| eαx[(ax+b)cos(βx)+(cx+d)sin(βx)] |
eαx[(Ax+B)cos(βx)+(Cx+D)sin(βx)] |
| eαx[(ax2+bx+c)cos(βx)+(dx2+ex+f)sin(βx)] |
eαx[(Ax2+Bx+C)cos(βx)+(Dx2+Ex+F)sin(βx)] |
5. α+iβ是特征根
| g(x)的形式 |
特解y∗(x)的形式 |
| eαx[acos(βx)+bsin(βx)] |
xeαx[Acos(βx)+Bsin(βx)] |
| eαx[(ax+b)cos(βx)+(cx+d)sin(βx)] |
xeαx[(Ax+B)cos(βx)+(Cx+D)sin(βx)] |
| eαx[(ax2+bx+c)cos(βx)+(dx2+ex+f)sin(βx)] |
xeαx[(Ax2+Bx+C)cos(βx)+(Dx2+Ex+F)sin(βx)] |
6. 一个小技巧稍微降低运算量
求解二阶常系数线性非齐次微分方程时,需要对特解进行二次求导,若特解是两个含x项的乘积,则求导时需要不断使用求导乘法法则,且需代入方程求解系数,这个过程的计算量比较大。有一个小技巧可以稍微降低运算量:
若特解是两个含x项的乘积,可将两项分别记为y1和y2此时特解为:y∗=y1y2先对特解求一阶导可得:(y∗)′=y1′y2+y1y2′再对特解求二阶导可得:(y∗)′′=y1′′y2+2y1′y2′+y1y2′′计算y1′,y1′′,y2′,y2′′,代入以上两式,整理并化简最后代入原方程:y′′(x)+py′(x)+qy(x)=0
建议将(y∗)′和(y∗)′′的式子记下来,以加快运算速度。
【例】y′′−3y′+2y=10e−xsinx
特征方程:r2−3r+2=0由此设齐次通解:y0=C1ex+C2e2x设非齐次特解:y∗=e−x(Asinx+Bcosx)记:y1=e−x,y2=Asinx+Bcosx,则:y∗=y1y2有:y1′=−e−x,y2′=Acosx−Bsinx有:y1′′=e−x,y2′′=−Asinx−Bcosx所以有:(y∗)′=y1′y2+y1y2′=e−x[(A−B)cosx−(A+B)sinx]所以有:(y∗)′′=y1′′y2+2y1′y2′+y1y2′′=e−x[2Bsinx−2Acosx]最后代入原方程化简:(A+B)sinx+(B−A)cosx=2sinx求出A=1,B=1所以:y∗=e−x(sinx+cosx)