Loading [MathJax]/jax/output/HTML-CSS/jax.js

微积分笔记(27)——常系数线性方程(1)

微积分笔记(27)——常系数线性方程(1)

高阶线性方程解的结构(续)

推广:n 阶线性方程

y(n)+a1(x)y(n1)++an1(x)y+an(x)y=f(x)    (xI)

通解:
y=y+c1y1++cnyn


其中 y1,,yn 是齐次方程的 n 个线性无关解,y 是非齐次方程的一个解。

常系数线性方程

常系数线性方程的特征方程

y(n)+a1y(n1)++an1y+any=f(x)    (1)

方法:求形如 y=eλx 的解(齐次方程 f0)。

代入方程:
(λn+a1λn1++an1λ+an)eλ=0λn+a1λn1++an1λ+an=0    (2)


方程 (2) 称为方程 (1)特征方程

求解方法之三:特征方程法——适用常系数情况

n=2 为例:
y+ay+by=0    (3)


特征方程:
λ2+aλ+b=0    (4)

分情况讨论:

  1. 若特征方程 (4)2 个不同的实根 λ1,λ2,则 y1=eλ1x,y2=eλ2x 都是 (3) 的解,W[y1,y2]=(λ2λ1)e(λ1+λ2)x0,从而通解 y=c1eλ1x+c2eλ2x

  2. 若特征方程 (4) 只有 12 重根 λ=a2,则 y1=eλx,y2=xeλx 是两个线性无关解。

    验证y2=(1+λx)eλx,y2=(2λ+λ2x)eλx

    代入方程得:
    y2+ay2+by2=[(2λ+λ2x)+a(1+λx)+bx]eλx=[x(λ2+aλ+b)+(a+2λ)]eλx=0


    W[y1,y2]=eλx0

    方程 (3) 的通解 y=(c1+c2x)ea2x

  3. 若特征方程 (4) 只有一对共轭复根,λ=α±iβ(β0),则 y±=eα±iβ(3)2 个解。

    回忆:Euler 公式(定义)eiθ=cosθ+isinθ。(导数相等)
    y±=eαx(cos(βx)±isin(βx))


    为了得到实值的解,利用线性性质,y1=12(y++y)=eαxcos(βx),y2=12i(y+y)=eαxsin(βx),得到 (4) 的两个解,W[y1,y2]=βeαx,方程 (1) 的通解 y=eαx(c1cos(βx)+c2sin(βx))

回到一般 n 的情况:

  1. λ 为特征方程 (2)k 重实根,则:

    y1=eλx,,yk=xk1eλx 为方程 (1)k 个线性无关解。

  2. λ=α±iβ 是特征方程 (2)k 重共轭复根,则 y1=eαxcos(βx),,yk=xk1eαxcos(βx),~y1=eαxsin(βx),,~yk=xk1eαxsin(βx) 为方程 (1)2k 个线性无关解。

非齐次方程(常系数)

y(n)+a1y(n1)++an1y+any=f(x)    (1)

假设 f(x)=Pn(x)eμxμ 可以是复数。

μ=α+iβeμx=eαx(cos(βx)+isin(βx))

方程化为:
L(y)=Pn(x)eαxcos(βx)+iPn(x)eαxsin(βx)=f1(x)+if2(x)


y=u(x)+iv(x),则:
L(y)=L(u)+iL(v)=f1(x)+if2(x)L(u)=f1(x),L(v)=f2(x)

 

点赞 0

No Comments

Add your comment

规矩是阻碍,谁定的规矩谁去守吧!