LaTeX在线编辑器:Equation Editor
一、定义
二阶常系数线性齐次微分方程:
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,r1≠r2(C1+C2x)er1x,r1=r2eαx[C1cos(βx)+C2sin(βx)],r=α±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βα+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βα+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∗)′=y′1y2+y1y′2再对特解求二阶导可得:(y∗)′′=y′′1y2+2y′1y′2+y1y′′2计算y′1,y′′1,y′2,y′′2,代入以上两式,整理并化简最后代入原方程: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有:y′1=−e−x,y′2=Acosx−Bsinx有:y′′1=e−x,y′′2=−Asinx−Bcosx所以有:(y∗)′=y′1y2+y1y′2=e−x[(A−B)cosx−(A+B)sinx]所以有:(y∗)′′=y′′1y2+2y′1y′2+y1y′′2=e−x[2Bsinx−2Acosx]最后代入原方程化简:(A+B)sinx+(B−A)cosx=2sinx求出A=1,B=1所以:y∗=e−x(sinx+cosx)
本文作者:漫舞八月(Mount256)
本文链接:https://www.cnblogs.com/Mount256/p/17502532.html
版权声明:本作品采用CC 4.0 BY-SA许可协议进行许可。
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 无需6万激活码!GitHub神秘组织3小时极速复刻Manus,手把手教你使用OpenManus搭建本
· Manus爆火,是硬核还是营销?
· 终于写完轮子一部分:tcp代理 了,记录一下
· 别再用vector<bool>了!Google高级工程师:这可能是STL最大的设计失误
· 单元测试从入门到精通