1500字范文,内容丰富有趣,写作好帮手!
1500字范文 > 常系数齐次线性微分方程

常系数齐次线性微分方程

时间:2024-05-30 06:30:45

相关推荐

常系数齐次线性微分方程

先讨论二阶常系数齐次线性微分方程的解法,再把二阶方程的解法推广到n阶方程。

在二阶齐次线性微分方程

y′′+P(x)y′+Q(x)y=0(1)y''+P(x)y'+Q(x)y=0 \tag{1} y′′+P(x)y′+Q(x)y=0(1)

中,如果y′,yy',yy′,y的系数P(x)P(x)P(x),Q(x)Q(x)Q(x)均为常数,即(1)式成为

y′′+py′+qy=0(2)y''+py'+qy=0 \tag{2} y′′+py′+qy=0(2)

其中p,q是常数,那么称(2)为二阶常系数齐次线性微分方程。如果p,q不全为常数,称(1)为二阶变系数齐次线性微分方程。

由之前讨论可知,要找微分方程(2)的通解,可以先求出它的两个解y1,y2y_1,y_2y1​,y2​,如果它们之比不为常数,即y1y_1y1​与y2y_2y2​线性无关,那么y=C1y1+C2y2y=C_1y_1+C_2y_2y=C1​y1​+C2​y2​就是方程(2)的通解。

当r为常数时,指数函数y=erxy=e^{rx}y=erx和它的各阶导数都只相差一个常熟因子。由于指数函数有这个特点,因此用y=erxy=e^{rx}y=erx来尝试,看能否选取适当的常数r,使y=erxy=e^{rx}y=erx满足方程(2)

将y=erxy=e^{rx}y=erx求导,得到

y′=rerx,y′′=r2erxy'=re^{rx},\quad y''=r^2e^{rx} y′=rerx,y′′=r2erx

把y,y′y,y'y,y′和y′′y''y′′代入方程(2),得

(r2+pr+q)erx=0(r^2+pr+q)e^{rx}=0 (r2+pr+q)erx=0

由于erx≠0e^{rx}\neq 0erx​=0,所以

r2+pr+1=0(3)r^2+pr+1=0 \tag{3} r2+pr+1=0(3)

由此可见,只要r满足代数方程(3),函数r=erxr=e^{rx}r=erx就是微分方程(2)的解,我们把代数方程(3)叫做微分方程(2)的特征方程。

特征方程(3)是一个二次代数方程,其中r2,rr^2,rr2,r的系数及常数项恰好依次是微分方程(2)中y′′,y′y'',y'y′′,y′及yyy的系数。

特征方程(3)的两个根r1,r2r_1,r_2r1​,r2​可以用公式

r1,2=−p±p2−4q2r_{1,2}=\frac{-p\pm \sqrt{p^2-4q}}{2} r1,2​=2−p±p2−4q​​

求出。它们有三种不同的情形:

(i)当P2−4q>0P^2-4q>0P2−4q>0时,r1,r2r_1,r_2r1​,r2​是两个不相等的实根

r1=−p+p2−4q2,r2=−p−p2−4q2r_1=\frac{-p+\sqrt{p^2-4q}}{2},\quad r_2=\frac{-p-\sqrt{p^2-4q}}{2} r1​=2−p+p2−4q​​,r2​=2−p−p2−4q​​

(ii)当p2−4q=0p^2-4q=0p2−4q=0时,r1,r2r_1,r_2r1​,r2​是两个相等的实根

r1=r2=−p2r_1=r_2=-\frac{p}{2} r1​=r2​=−2p​

(iii)当p2−4q<0p^2-4q<0p2−4q<0时,r1,r2r_1,r_2r1​,r2​是一对共轭复根

r1=α+βi,r2=α−βir_1=\alpha+\beta i,\quad r_2=\alpha-\beta i r1​=α+βi,r2​=α−βi

其中

α=−p2,β=4q−p22\alpha=-\frac{p}{2},\quad \beta=\frac{\sqrt{4q-p^2}}{2} α=−2p​,β=24q−p2​​

相应地,微分方程(2)的通解也有三种不同的情形。分别讨论如下:

(i)特征方程有两个不相等的实根:r1≠r2r_1\neq r_2r1​​=r2​

由上面讨论知道,y1=er1x,y2=er2xy_1=e^{r_1x},y_2=e^{r_2x}y1​=er1​x,y2​=er2​x是微分方程(2)的两个解,并且y2y1=er2xer1x=e(r2−r1)x\frac{y_2}{y_1}=\frac{e^{r_2x}}{e^{r_1x}}=e^{(r_2-r_1)x}y1​y2​​=er1​xer2​x​=e(r2​−r1​)x不是常数,因此微分方程(2)的通解为

y=C1er1x+C2er2xy=C_1e^{r_1x}+C_2e^{r_2x} y=C1​er1​x+C2​er2​x

(ii)特征方程有两个相等的实根:r1=r2r_1=r_2r1​=r2​

这时,只得到微分方程(2)的一个解

y1=er1xy_1=e^{r_1x} y1​=er1​x

为了得出微分方程(2)的通解,还需求出另一个解y2y_2y2​,并且要求y2y1\frac{y_2}{y_1}y1​y2​​不是常数。设y2y1=u(x)\frac{y_2}{y_1}=u(x)y1​y2​​=u(x),即y2=er1xu(x)y_2=e^{r_1x}u(x)y2​=er1​xu(x)。下面来求u(x)u(x)u(x)。将y2y_2y2​求导,得

y2′=er1x(u′+r1u)y2′′=er1x(u′′+2r1u′+r12u)y_2'=e^{r_1x}(u'+r_1u) \\ y_2''=e^{r_1x}(u''+2r_1u'+r_1^2u) y2′​=er1​x(u′+r1​u)y2′′​=er1​x(u′′+2r1​u′+r12​u)

将y2,y2′y_2,y_2'y2​,y2′​和y2′′y_2''y2′′​代入微分方程(2),得

er1x[(u′′+2r1u′+r12u)+p(u′+r1u)+qu]=0e^{r_1x}[(u''+2r_1u'+r_1^2u)+p(u'+r_1u)+qu]=0 er1​x[(u′′+2r1​u′+r12​u)+p(u′+r1​u)+qu]=0

约去er1xe^{r_1x}er1​x,并合并同类项,得

u′′+(2r1+p)u′+(r12+pr1+q)u=0u''+(2r_1+p)u'+(r_1^2+pr_1+q)u=0 u′′+(2r1​+p)u′+(r12​+pr1​+q)u=0

由于r1r_1r1​是特征方程(3)的二重根。因此r12+pr1+q=0r_1^2+pr_1+q=0r12​+pr1​+q=0,且2r1+p=02r_1+p=02r1​+p=0,于是得

u′′=0u''=0 u′′=0

因为这里只要得到一个不为常数的解,所以不妨选取u=x,由此得到微分方程(2)的另一个解

y2=xer1xy_2=xe^{r_1x} y2​=xer1​x

从而微分方程(2)的通解为

y=C1er1x+C2xer1xy=C_1e^{r_1x}+C_2xe^{r_1x} y=C1​er1​x+C2​xer1​x

y=(C1+C2x)er1xy=(C_1+C_2x)e^{r_1x} y=(C1​+C2​x)er1​x

(iii)特征方程有一对共轭复根:r1=α+βi,r2=α−βi(β≠0)r_1=\alpha+\beta i,r_2=\alpha-\beta i(\beta \neq 0)r1​=α+βi,r2​=α−βi(β​=0)

这时,y1=e(α+βi)x,y2=e(α−βi)xy_1=e^{(\alpha+\beta i)x},y_2=e^{(\alpha-\beta i)x}y1​=e(α+βi)x,y2​=e(α−βi)x是微分方程(2)的两个解,但它们是复值函数形式。为了得出实值函数形式的解,先利用欧拉公式eiθ=cosθ+isinθe^{i\theta}=cos\theta+isin\thetaeiθ=cosθ+isinθ把y1,y2y_1,y_2y1​,y2​改写为

y1=e(α+βi)x=eax⋅eβxi=eax(cosβx+isinβx)y2=e(α−βi)x=eax⋅e−βxi=eax(cosβx−isinβx)y_1=e^{(\alpha+\beta i)x}=e^{ax}·e^{\beta xi}=e^{ax}(cos\beta x+isin\beta x) \\ y_2=e^{(\alpha-\beta i)x}=e^{ax}·e^{-\beta xi}=e^{ax}(cos\beta x-isin\beta x) y1​=e(α+βi)x=eax⋅eβxi=eax(cosβx+isinβx)y2​=e(α−βi)x=eax⋅e−βxi=eax(cosβx−isinβx)

由于复值函数y1y_1y1​与y2y_2y2​之间成共轭关系,因此,取它们的和除以2就得到它们的实部,取它们的差除以2i就得到它们的虚部。由于方程(2)的解复合叠加原理,所以实值函数

y^1=12(y1+y2)=eaxcosβx,y^2=12i(y1−y2)=eaxsinβx\hat y_1=\frac{1}{2}(y_1+y_2)=e^{ax}cos\beta x ,\\ \hat y_2=\frac{1}{2i}(y_1-y_2)=e^{ax}sin\beta x y^​1​=21​(y1​+y2​)=eaxcosβx,y^​2​=2i1​(y1​−y2​)=eaxsinβx

还是微分方程(2)的解,且y^1y^2=eaxcosβxeaxsinβx=cotβx\frac{\hat y_1}{\hat y_2}=\frac{e^{ax}cos\beta x}{e^{ax}sin\beta x}=cot \,\beta xy^​2​y^​1​​=eaxsinβxeaxcosβx​=cotβx不是常数,所以微分方程(2)的通解为

y=eax(C1cosβx+C2sinβx)y=e^{ax}(C_1cos\beta x+C_2sin\beta x) y=eax(C1​cosβx+C2​sinβx)

综上所述,求二阶常系数齐次线性微分方程

y′′+py′+qy=0y''+py'+qy =0 y′′+py′+qy=0

的通解的步骤如下:

第一步:写出微分方程(2)的特征方程

r2+pr+q=0(3)r^2+pr+q=0 \tag{3} r2+pr+q=0(3)

第二步:求出特征方程(3)的两个根r1,r2r_1,r_2r1​,r2​。

第三步:根据特征方程(3)的两个根的不同情形,按照下列表格写出微分方程(2)的通解:

本内容不代表本网观点和政治立场,如有侵犯你的权益请联系我们处理。
网友评论
网友评论仅供其表达个人看法,并不表明网站立场。