二阶常系数非齐次微分方程求解

二阶常系数非齐次微分方程求解

Last updated on


385 字 约 2 分钟

设有二阶常系数非齐次微分方程

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

考虑非齐次方程的通解由两部分组成:

y=yh+yy = y_{h} + y^*

其中 yhy_{h} 为齐次方程 y+py+y=0y'' + py' + y = 0 的通解,yy^* 为非齐次特解。


首先求齐次方程的通解:

设有特征方程

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

对特征根的存在情况进行分类:

判别式 Δ=p24q\Delta = p^{2} - 4q特征根情况齐次通解 yhy_{h}
Δ>0\Delta > 0两个不同实根 r1,r2r_{1}, r_{2}C1er1x+C2er2xC_{1} e^{r_{1}x} + C_{2} e^{r_{2}x}
Δ=0\Delta = 0两个相同实根 rr(C1+C2x)erx(C_{1} + C_{2} x) e^{rx}
Δ<0\Delta < 0一对共轭复根 α±βi\alpha \pm \beta ieαx(C1cosβx+C2sinβx)e^{\alpha x}(C_{1} cos \beta x + C_{2} sin \beta x)

随后求非齐次方程的特解,我们根据 f(x)f(x) 的形式进行预期设解。

f(x)f(x)yy^* 预设形式
Pn(x)P_{n}(x)xkQn(x)x^{k} Q_{n} (x)
Pn(x)eλxP_{n}(x) e^{\lambda x}xkQn(x)eλxx^{k} Q_{n}(x) e^{\lambda x}
eλx[Pl(x)cosωx+Pm(x)sinωx]e^{\lambda x} [P_{l}(x) cos \omega x + P_{m}(x) sin \omega x]xkeλx[Rn(x)cosωx+Sn(x)sinωx]x^{k} e^{\lambda x} [R_{n}(x) cos \omega x + S_{n}(x) sin \omega x]

其中对于 kk ,即为特征根与 f(x)f(x) 中指数部分 (λ+ωi)(\lambda + \omega i) 的重合次数。

  • 如果 λ+ωi\lambda + \omega i 不是特征方程的根, k=0k = 0
  • 如果 λ+ωi\lambda + \omega i 是特征方程的单根, k=1k = 1
  • 如果 λ+ωi\lambda + \omega i 是特征方程的重根, k=2k = 2

其中对于所谓指数部分,即

eλ+ωi=eλeωie^{\lambda + \omega i} = e^{\lambda} \cdot e^{\omega i}

有欧拉公式:

eλ(cosω+isinω)e^{\lambda} (cos \omega + i sin \omega)

欧拉公式的本质,是将微分这个高成本的操作,转化为乘法这个低成本的操作。

对于方程

y+py+qy=eλxcos(ωx)y'' + py' + qy = e^{\lambda x} cos (\omega x)

将其视为复指数函数的实部

y+py+qy=Re[e(λ+ωi)x]y'' + py' + qy = \text{Re} [e^{(\lambda + \omega i) x}]
🐙