769 字
4 分钟
常系数(非)齐次线性微分方程

齐次线性微分方程求通解#

y(n)+an1y(n1)++a1y+a0y=0\large y^{(n)}+a_{n-1}y^{(n-1)}+\cdots +a_{1}y^{\prime}+a_{0}y=0为常系数齐次线性微分方程
λn+an1λn1++a1λ+a0=0\large \lambda ^{n}+a_{n-1}\lambda ^{n-1}+\cdots +a_{1}\lambda +a_{0}=0为它的特征方程

特征值四种情况#

  • λ\large \lambda是单实特征根,则y=eλx\large y=e^{\lambda x} (基础解系)
  • λ\large \lambdak\large k重实特征根,则基础解系中的k\large k个解y1=eλx,y2=xeλx,,yk=xk1eλx\large y_1=e^{\lambda x},y_2=xe^{\lambda x},\cdots,y_k=x^{k-1}e^{\lambda x}
  • λ=α±iβ\large \lambda=\alpha \pm i\beta是单复特征根,则方程基础解系中对应的两个解y1=eαxcosβx,y2=eαxsinβx\large y_1=e^{\alpha x}\cos {\beta x},y_2=e^{\alpha x}\sin{\beta x}
  • λ=α±iβ\large \lambda=\alpha \pm i\betak\large k重复特征根,则对应2k\large 2k个解
y1=eαxcosβx,y2=eαxsinβxxeαxcosβx,xeαxsinβx,xk1eαxcosβx,xk1eαxsinβx\begin{aligned} \large y_1=e^{\alpha x}\cos{\beta x}&\large,y_2=e^{\alpha x}\sin{\beta x} \\ \large xe^{\alpha x}\cos{\beta x}&\large,xe^{\alpha x}\sin{\beta x} \\ \large \cdots&\large,\cdots \\ \large x^{k-1}e^{\alpha x}\cos{\beta x}&\large,x^{k-1}e^{\alpha x}\sin{\beta x} \end{aligned}
TIP

由欧拉公式eiθ=cosθ±isinθ\large e^{i\theta}=\cos{\theta} \pm i\sin{\theta} 得出

NOTE

例题:
y+3y+2y=0\large y^{\prime \prime} +3y^{\prime}+2y=0的通解

解:特征方程λ2+3λ+2=0\large \lambda ^{2}+3\lambda +2=0
特征根为λ1=1,λ2=2\large \lambda_{1}=-1,\lambda_2=-2
故通解为y=C1ex+C2e2x\large y=C_{1}e^{-x}+C_{2}e^{-2x}

NOTE

例题:
y+2y+y=0\large y^{\prime \prime} +2y^{\prime}+y=0的通解

解:特征方程λ2+2λ+1=0\large \lambda ^{2}+2\lambda +1=0
特征根为λ1=λ2=1\large \lambda_{1}=\lambda_2=-1
故通解为y=C1ex+C2xex\large y=C_{1}e^{-x}+C_{2}xe^{-x}

NOTE

例题:
y+y=0\large y^{\prime \prime} +y=0的通解

解:特征方程λ2+1=0\large \lambda ^{2}+1=0
特征根为λ1=i,λ2=i\large \lambda_{1}=i,\lambda_2=-i
(α=0,β=1\large \alpha = 0,\beta = 1 代入公式)
故通解为y=C1cosx+C2sinx\large y=C_{1}\cos{x}+C_{2}\sin{x}

非齐次#

y(n)+an1y(n1)++a1y+a0y=f(x)\large y^{(n)}+a_{n-1}y^{(n-1)}+\cdots +a_{1}y^{\prime}+a_{0}y=f(x)为常系数非齐次线性微分方程
区别就是右边的0变成了关于x的函数

求通解#

方法:非齐次方程特解加对应的齐次方程通解(跟线代里学的线性方程组一样)

求特解#

  • f(x)\large f(x)多项式型

f(x)=Pn(x)\large f(x)=P_n(x)

k\large k为0作为该方程特征根的重数,0不是特征根则重数为0,k=0\large k=0

可设特解y=xkQn(x)\large y^{\ast}=x^{k}Q_n(x)

Qn(x)=anxn++a1x+a0\large Q_n(x)=a_nx^{n}+\cdots+a_1x+a_0

NOTE

例题:
yy=x2\large y^{\prime \prime} -y^{\prime}=x^{2}特解

解:设y=ax3+bx2+cx+d,y=3ax2+2bx+c,y=6ax+2b\large y^{\ast}=ax^{3}+bx^{2}+cx+d, y^{\ast \prime}=3ax^{2}+2bx+c,y^{\ast \prime \prime}=6ax+2b
(设ax3\large ax^{3}是因为等号右侧最高次为x2\large x^{2},左侧最低阶导数为1阶,ax3\large ax^{3}求导一次得3ax2\large 3ax^{2})
代入方程,6ax+2b3ax22bxc=x2\large 6ax+2b-3ax^{2}-2bx-c=x^{2}
解出a,b,c\large a,b,c 得特解

  • f(x)\large f(x)非多项式型

一般可设f(x)=Pn(x)eαxcosβx+Qm(x)eαxsinβx\large f(x)=P_{n}(x)e^{\alpha x}\cos{\beta x}+Q_{m}(x)e^{\alpha x}\sin{\beta x}

Pn(x)\large P_{n}(x)Qm(x)\large Q_{m}(x)为多项式

特解为y=xk[P~t(x)eαxcosβx+Q~t(x)eαxsinβx]\large y^{\ast}=x^{k}\left [ \tilde{P} _{t}(x)e^{\alpha x}\cos{\beta x}+\tilde{Q} _{t}(x)e^{\alpha x}\sin{\beta x}\right ]

k\large kα±iβ\large \alpha \pm i\beta作为特征根时的重数且t=max{n,m}\large t=\max \left \{ n,m \right \}

k\large k分为两种情况

  • α±iβ\large \alpha \pm i\beta为方程的n\large n重特征根,则k=n\large k=n,特解中要多乘xk\large x^{k}
  • α±iβ\large \alpha \pm i\beta不是特征根则重数为0,则k=0\large k=0,也就是乘x0=1\large x^{0}=1

所以解特解的思路就变成了先将f(x)\large f(x)化成Pn(x)eαxcosβx+Qm(x)eαxsinβx\large P_{n}(x)e^{\alpha x}\cos{\beta x}+Q_{m}(x)e^{\alpha x}\sin{\beta x}的形式,进而就求得了α,β\large \alpha,\beta的值,得k\large k的值为α±iβ\large \alpha \pm i\beta作为方程特征根时的重数

NOTE

例题:
y+y=x2+3e2x\large y^{\prime \prime} +y\prime=x-2+3e^{2x}特解

其实是分成了两个方程x2\large x-23e2x\large 3e^{2x}分别求特解再相加
解:特征方程λ2+λ=0\large \lambda ^{2}+\lambda=0
特征根为λ1=0,λ2=1\large \lambda_1=0,\lambda_2=-1
对于x2\large x-2的部分
α±iβ=0±i0=0\large \alpha \pm i\beta=0\pm i0=0是1重特征根,所以k1=1\large k_1 = 1
所以设y1=xk1(ax+b)=x(ax+b)\large y^{\ast}_1=x^{k_1}\left ( ax+b \right )=x\left ( ax+b \right )
对于方程3e2x\large 3e^{2x}的部分
α±iβ=2±i0=2\large \alpha \pm i\beta=2\pm i0=2,不是3e2x\large 3e^{2x}的特征根,所以重数为0,k2=0\large k_2=0
所以设y2=xk2[ce2xcos0x+ae2xsin0x]=ce2x\large y^{\ast}_2=x^{k_2}\left [ ce^{2x}\cos{0x}+ae^{2x}\sin{0x} \right ]=ce^{2x}
因此y=y1+y2=x(ax+b)+ce2x\large y^{\ast}=y_1^{\ast}+y_2^{\ast}=x(ax+b)+ce^{2x}
代入求解a,b,c\large a,b,c即可

常系数(非)齐次线性微分方程
https://a1kari8.github.io/posts/const_coefficient_diff_equa/
作者
A1kari8
发布于
2025-01-09
许可协议
CC BY-NC-SA 4.0