高阶线性微分方程

区别于初等积分法用积分找解。实际上n>=2就算高阶了

就拿二阶举例子吧

y+py+qy=f(x)

pq是函数,那么就是变系数。有f(x)是非齐次

高等数学的语境里,高阶只包含了常系数

常系数高阶线性微分方程

解的结构

两个解,比值不是常数,就叫成了线性无关的解,实际上有线性代数的影子。

解的结构实际上也和线性代数方程组的解的结构很相似。也是齐通+非齐特

更一般的

y(n)+an1y(n1)++a1y˙+a0y=f(x)

和对应的齐次线性方程

y(n)+an1y(n1)++a1y˙+a0y=0

借助微分方程组的结论,上面这个齐次方程等价于

dxdt=Ax A=[010000100001a0a1a2an1]

这时候就可以有这个矩阵的特征行列式和特征方程了,虽然是通过这个矩阵定义出来的,但是这个特征方程也被叫做了微分方程的特征方程,解叫作特征根。

特征根对应的模态就是齐次方程组的解,一堆模态搞起来作为基本解组。(说的太含糊了)

有复数根的基本解组使用EUler公式处理。有重根….

很严谨,但是搞复杂了。

还是拿二阶举例子,用有代表性的二阶研究,这里pq是常数。

y+py+qy=0

很明显满足这个式子的是指数函数。

特征式λ2+pλ+q=0

解出来就是e上面的指数。

  • 两个不同根,就是两个模态的线性组合
y=C1eλ1x+C2eλ2x
  • 二重根,模态就线性相关了,xeλx也是个解,而且无关
y=(C1+C2x)eλx
  • 共轭复根λ1,2=α±βi
y=eαx(C1cosβx+C2sinβx)
补充

控制系统,二阶系统,共轭复根,复数引入震动(三角函数),实部包络线,震动幅值衰减还是增加的问题。

非齐次的问题,找特解。

特解是设出来的。。很痛苦

y+py+qy=Pn(x)eαx

当自由项是Pn(x)eαx,设特解为eαxQn(x)xk,其中Qn(x)为待定系数的同阶多项式,k取值0、1、2,

y+py+qy=Pn(x)eαx 特征根 通解 特解
y4y+4y=3xe2x 2,2 (C1+C2x)e2x e2x(ax+b)x2
y4y+3y=xe3x 1,3 C1ex+C2e3x e3x(ax+b)x1
  • 更痛苦的,当自由项出现了三角函数了
y+py+qy=[Pm(x)cosβx+Pn(x)sinβx]eαx

设特解为eαx[Ql(1)(x)cosβx+Ql(1)(x)sinβx]xk,其中Ql(x)为待定系数的同阶多项式,取mn的最大值,k取值0、1,

y+py+qy=[Pm(x)cosβx+Pn(x)sinβx]eαx 特征根 通解 特解
y4y+4y=e2xcosβx 2,2 (C1+C2x)e2x e2x(ax+b)x2
y4y+3y=xe3x 1,3 C1ex+C2e3x e3x(ax+b)x1