LaTeX在线编辑器:Equation Editor

[toc]

一、定义

二阶常系数线性齐次微分方程:

y(x)+py(x)+qy(x)=0y''(x)+py'(x)+qy(x)=0

二阶常系数线性非齐次微分方程:

y(x)+py(x)+qy(x)=g(x)y''(x)+py'(x)+qy(x)=g(x)

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

y(x)=齐次通解+非齐次特解=y0(x)+y(x)y(x)=齐次通解+非齐次特解=y_0(x)+y^*(x)

二、齐次通解

特征方程为

r2+pr+q=0r^2+pr+q=0

根据特征方程的根r1,r2r_1,r_2的情况,设通解为

y0(x)={C1er1x+C2er2x,r1r2(C1+C2x)er1x,r1=r2eαx[C1cos(βx)+C2sin(βx)],r=α±iβy_0(x)= \begin{cases} C_1e^{r_1x}+C_2e^{r_2x}, & r_1 \neq r_2 \\ (C_1+C_2x)e^{r_1x}, & r_1 = r_2 \\ e^{\alpha x}[C_1\cos(\beta x) + C_2\sin(\beta x)], & r = \alpha \pm \mathrm{i} \beta \end{cases}

三、非齐次特解

1. α\alpha不是特征根

g(x)g(x)的形式 特解y(x)y^*(x)的形式
aa AA
ax+bax+b Ax+BAx+B
ax2+bx+cax^2+bx+c Ax2+Bx+CAx^2+Bx+C
aeαxae^{\alpha x} AeαxAe^{\alpha x}
(ax+b)eαx(ax+b)e^{\alpha x} (Ax+B)eαx(Ax+B)e^{\alpha x}
(ax2+bx+c)eαx(ax^2+bx+c)e^{\alpha x} (Ax2+Bx+C)eαx(Ax^2+Bx+C)e^{\alpha x}

2. α\alpha是特征单重根

(1)α0\alpha \neq 0

g(x)g(x)的形式 特解y(x)y^*(x)的形式
aeαxae^{\alpha x} xAeαxxAe^{\alpha x}
(ax+b)eαx(ax+b)e^{\alpha x} x(Ax+B)eαxx(Ax+B)e^{\alpha x}
(ax2+bx+c)eαx(ax^2+bx+c)e^{\alpha x} x(Ax2+Bx+C)eαxx(Ax^2+Bx+C)e^{\alpha x}

(2)α=0\alpha=0

g(x)g(x)的形式 特解y(x)y^*(x)的形式
aa xAxA
ax+bax+b x(Ax+B)x(Ax+B)
ax2+bx+cax^2+bx+c x(Ax2+Bx+C)x(Ax^2+Bx+C)

3. α\alpha是特征二重根

(1)α0\alpha \neq 0

g(x)g(x)的形式 特解y(x)y^*(x)的形式
aeαxae^{\alpha x} x2Aeαxx^2Ae^{\alpha x}
(ax+b)eαx(ax+b)e^{\alpha x} x2(Ax+B)eαxx^2(Ax+B)e^{\alpha x}
(ax2+bx+c)eαx(ax^2+bx+c)e^{\alpha x} x2(Ax2+Bx+C)eαxx^2(Ax^2+Bx+C)e^{\alpha x}

(2)α=0\alpha=0

g(x)g(x)的形式 特解y(x)y^*(x)的形式
aa x2Ax^2A
ax+bax+b x2(Ax+B)x^2(Ax+B)
ax2+bx+cax^2+bx+c x2(Ax2+Bx+C)x^2(Ax^2+Bx+C)

4. α+iβ\alpha+\mathrm{i}\beta不是特征根

g(x)g(x)的形式 特解y(x)y^*(x)的形式
eαx[acos(βx)+bsin(βx)]e^{\alpha x}[a\cos(\beta x) + b\sin(\beta x)] eαx[Acos(βx)+Bsin(βx)]e^{\alpha x}[A\cos(\beta x) + B\sin(\beta x)]
eαx[(ax+b)cos(βx)+(cx+d)sin(βx)]e^{\alpha x}[(ax+b)\cos(\beta x) + (cx+d)\sin(\beta x)] eαx[(Ax+B)cos(βx)+(Cx+D)sin(βx)]e^{\alpha x}[(Ax+B)\cos(\beta x) + (Cx+D)\sin(\beta x)]
eαx[(ax2+bx+c)cos(βx)+(dx2+ex+f)sin(βx)]e^{\alpha x}[(ax^2+bx+c)\cos(\beta x) + (dx^2+ex+f)\sin(\beta x)] eαx[(Ax2+Bx+C)cos(βx)+(Dx2+Ex+F)sin(βx)]e^{\alpha x}[(Ax^2+Bx+C)\cos(\beta x) + (Dx^2+Ex+F)\sin(\beta x)]

5. α+iβ\alpha+\mathrm{i}\beta是特征根

g(x)g(x)的形式 特解y(x)y^*(x)的形式
eαx[acos(βx)+bsin(βx)]e^{\alpha x}[a\cos(\beta x) + b\sin(\beta x)] xeαx[Acos(βx)+Bsin(βx)]xe^{\alpha x}[A\cos(\beta x) + B\sin(\beta x)]
eαx[(ax+b)cos(βx)+(cx+d)sin(βx)]e^{\alpha x}[(ax+b)\cos(\beta x) + (cx+d)\sin(\beta x)] xeαx[(Ax+B)cos(βx)+(Cx+D)sin(βx)]xe^{\alpha x}[(Ax+B)\cos(\beta x) + (Cx+D)\sin(\beta x)]
eαx[(ax2+bx+c)cos(βx)+(dx2+ex+f)sin(βx)]e^{\alpha x}[(ax^2+bx+c)\cos(\beta x) + (dx^2+ex+f)\sin(\beta x)] xeαx[(Ax2+Bx+C)cos(βx)+(Dx2+Ex+F)sin(βx)]xe^{\alpha x}[(Ax^2+Bx+C)\cos(\beta x) + (Dx^2+Ex+F)\sin(\beta x)]

6. 一个小技巧稍微降低运算量

求解二阶常系数线性非齐次微分方程时,需要对特解进行二次求导,若特解是两个含xx项的乘积,则求导时需要不断使用求导乘法法则,且需代入方程求解系数,这个过程的计算量比较大。有一个小技巧可以稍微降低运算量:

若特解是两个含x项的乘积,可将两项分别记为y1y2此时特解为:y=y1y2先对特解求一阶导可得:(y)=y1y2+y1y2再对特解求二阶导可得:(y)=y1y2+2y1y2+y1y2计算y1y1y2y2,代入以上两式,整理并化简最后代入原方程:y(x)+py(x)+qy(x)=0\begin{aligned} &若特解是两个含x项的乘积,可将两项分别记为y_1和y_2 \\ &此时特解为:y^*=y_1y_2 \\ &先对特解求一阶导可得:(y^*)'=y_1'y_2+y_1y_2' \\ &再对特解求二阶导可得:(y^*)''=y_1''y_2+2y_1'y_2'+y_1y_2'' \\ &计算y_1',y_1'',y_2',y_2'',代入以上两式,整理并化简 \\ &最后代入原方程:y''(x)+py'(x)+qy(x)=0 \\ \end{aligned}

建议将(y)(y^*)'(y)(y^*)''的式子记下来,以加快运算速度。

【例】y3y+2y=10exsinxy''-3y'+2y = 10e^{-x} \sin x

特征方程:r23r+2=0由此设齐次通解:y0=C1ex+C2e2x设非齐次特解:y=ex(Asinx+Bcosx)记:y1=exy2=Asinx+Bcosx,则:y=y1y2有:y1=exy2=AcosxBsinx有:y1=exy2=AsinxBcosx所以有:(y)=y1y2+y1y2=ex[(AB)cosx(A+B)sinx]所以有:(y)=y1y2+2y1y2+y1y2=ex[2Bsinx2Acosx]最后代入原方程化简:(A+B)sinx+(BA)cosx=2sinx求出A=1,B=1所以:y=ex(sinx+cosx)\begin{aligned} &特征方程:r^2-3r+2=0 \\ &由此设齐次通解:y_0=C_1e^x+C_2e^{2x} \\ &设非齐次特解:y^*=e^{-x}(A \sin x+B \cos x) \\ &记:y_1=e^{-x},y_2=A \sin x+B \cos x,则:y^*=y_1y_2 \\ &有:y_1'=-e^{-x},y_2'=A \cos x-B \sin x \\ &有:y_1''=e^{-x},y_2''=-A \sin x-B \cos x \\ &所以有:(y^*)'=y_1'y_2+y_1y_2'=e^{-x}[(A-B)\cos x-(A+B) \sin x] \\ &所以有:(y^*)''=y_1''y_2+2y_1'y_2'+y_1y_2''=e^{-x}[2B\sin x-2A \cos x] \\ &最后代入原方程化简:(A+B)\sin x+(B-A)\cos x=2 \sin x \\ &求出A=1,B=1 \\ &所以:y^*=e^{-x}(\sin x+\cos x) \end{aligned}