G
N
I
D
A
O
L

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

LaTeX在线编辑器:Equation Editor

一、定义

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

\[y''(x)+py'(x)+qy(x)=0 \]

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

\[y''(x)+py'(x)+qy(x)=g(x) \]

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

\[y(x)=齐次通解+非齐次特解=y_0(x)+y^*(x) \]

二、齐次通解

特征方程为

\[r^2+pr+q=0 \]

根据特征方程的根\(r_1,r_2\)的情况,设通解为

\[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)\)的形式 特解\(y^*(x)\)的形式
\(a\) \(A\)
\(ax+b\) \(Ax+B\)
\(ax^2+bx+c\) \(Ax^2+Bx+C\)
\(ae^{\alpha x}\) \(Ae^{\alpha x}\)
\((ax+b)e^{\alpha x}\) \((Ax+B)e^{\alpha x}\)
\((ax^2+bx+c)e^{\alpha x}\) \((Ax^2+Bx+C)e^{\alpha x}\)

2. \(\alpha\)是特征单重根

(1)\(\alpha \neq 0\)

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

(2)\(\alpha=0\)

\(g(x)\)的形式 特解\(y^*(x)\)的形式
\(a\) \(xA\)
\(ax+b\) \(x(Ax+B)\)
\(ax^2+bx+c\) \(x(Ax^2+Bx+C)\)

3. \(\alpha\)是特征二重根

(1)\(\alpha \neq 0\)

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

(2)\(\alpha=0\)

\(g(x)\)的形式 特解\(y^*(x)\)的形式
\(a\) \(x^2A\)
\(ax+b\) \(x^2(Ax+B)\)
\(ax^2+bx+c\) \(x^2(Ax^2+Bx+C)\)

4. \(\alpha+\mathrm{i}\beta\)不是特征根

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

5. \(\alpha+\mathrm{i}\beta\)是特征根

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

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

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

\[\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''-3y'+2y = 10e^{-x} \sin x\)

\[\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} \]

posted @ 2023-06-25 11:34  漫舞八月(Mount256)  阅读(660)  评论(0编辑  收藏  举报