
微积分笔记(27)——常系数线性方程(1)
高阶线性方程解的结构(续)
推广:n 阶线性方程
y(n)+a1(x)y(n–1)+⋯+an–1(x)y′+an(x)y=f(x) (x∈I)
通解:
y=y∗+c1y1+⋯+cnyn
其中 y1,⋯,yn 是齐次方程的 n 个线性无关解,y∗ 是非齐次方程的一个解。
常系数线性方程
常系数线性方程的特征方程
y(n)+a1y(n–1)+⋯+an–1y′+any=f(x) (1)
方法:求形如 y=eλx 的解(齐次方程 f≡0)。
代入方程:
(λn+a1λn–1+⋯+an–1λ+an)eλ=0∴λn+a1λn–1+⋯+an–1λ+an=0 (2)
方程 (2) 称为方程 (1) 的特征方程。
求解方法之三:特征方程法——适用常系数情况
以 n=2 为例:
y′′+ay′+by=0 (3)
特征方程:
λ2+aλ+b=0 (4)
分情况讨论:
- 若特征方程 (4) 有 2 个不同的实根 λ1,λ2,则 y1=eλ1x,y2=eλ2x 都是 (3) 的解,W[y1,y2]=(λ2–λ1)e(λ1+λ2)x≠0,从而通解 y=c1eλ1x+c2eλ2x。
若特征方程 (4) 只有 1 个 2 重根 λ=–a2,则 y1=eλx,y2=xeλx 是两个线性无关解。
验证:y′2=(1+λx)eλx,y′′2=(2λ+λ2x)eλx。
代入方程得:
y′′2+ay′2+by2=[(2λ+λ2x)+a(1+λx)+bx]eλx=[x(λ2+aλ+b)+(a+2λ)]eλx=0
W[y1,y2]=eλx≠0。∴ 方程 (3) 的通解 y=(c1+c2x)e−a2x。
若特征方程 (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 的情况:
- 若 λ 为特征方程 (2) 的 k 重实根,则:
y1=eλx,⋯,yk=xk–1eλx 为方程 (1) 的 k 个线性无关解。
若 λ=α±iβ 是特征方程 (2) 的 k 重共轭复根,则 y1=eαxcos(βx),⋯,yk=xk–1eαxcos(βx),~y1=eαxsin(βx),⋯,~yk=xk–1eαxsin(βx) 为方程 (1) 的 2k 个线性无关解。
非齐次方程(常系数)
y(n)+a1y(n–1)+⋯+an–1y′+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)
No Comments