応用数学 第6回 (4) 定数係数同次線形微分方程式の解
定数係数同次線形微分方程式の解
f(X) を前ページの通りとするとき、定数係数
n 階同次線形微分方程式
f(D)y=0 ⋯⋯ (∗)
の解の公式を作ります。
Def.7 n 次方程式
f(X)=0
を (∗) の特性方程式と言い、その根を特性根と言う。
Th.8
- y=eαx が (∗) の解であること ⇔ α が特性根であること
- α が特性方程式の k 重根ならば、
eαx, xeαx, ⋯, xk−1eαx は全て (∗) の解である
- f(X) の全ての根 α について (2) の形の関数を書き出すと、それらはちょうど n 個あり、(∗) の解空間の基底を成す。
証明 (1)
Th.6 (1) より
f(D)eαx=0
⇔ f(α)eαx=0
⇔ f(α)=0.
- f(X)=g(X)(X−α)k と書けますので、Th.6 (2) より
f(D){eαxxj}=eαxf(D+α)xj=eαxg(D+α)Dkxj.
j<k ならば Dkxj=0 ゆえ右辺は 0 になります。
- (2) の形の関数が n 個あることは根の重複度を考えればわかります。
それらが一次独立であることは少し長くなるので省略します。( 付録参照 )
Rem.9 特性根が複素数
α=λ+iμ のときは、共役複素数
ˉα=λ−iμ も特性根で、
それらの重複度を
k とすると
- eλxcos(μx), xeλxcos(μx), ⋯, xk−1eλxcos(μx)
- eλxsin(μx), xeλxsin(μx), ⋯, xk−1eλxsin(μx)
たちが
2k 次元分の基底になります。
やってみよう
(1) y‴−5y″+6y′=0
解
f(D)=D3−5D2+6D の場合で、特性方程式は
X3−5X2+6X=X(X−2)(X−3)=0
特性根は
0,
2,
3 ( 全て単根 ) で、一般解は
y=Ae0x+Be2x+Ce3x=A+Be2x+Ce3x
(2) (D−1)4(D+3)2y=0
解
特性根は 4 重根
1 と 2 重根
−3 なので、一般解は
y=(A+Bx+Cx2+Dx3)ex+(E+Fx)e−3x
(3) (D2+4D+13)2y=0
解
特性方程式は
(X2+4X+13)2=(X+2+3i)2(X+2−3i)2=0
特性根は
−2±3i ( 共に 2 重根 ) なので、一般解は
y=(A+Bx)e(−2+3i)x+(C+Dx)e(−2−3i)x=(A+Bx)e−2xei(3x)+(C+Dx)e−2xei(−3x)=(A+Bx)e−2x{cos(3x)+isin(3x)}+(C+Dx)e−2x{cos(−3x)+isin(−3x)}=e−2x{((A+C)+(B+D)x)cos(3x)+i((A−C)+(B−D)x)sin(3x)}=e−2x{(A′+B′x)cos(3x)+(C′+D′x)sin(3x)}