数値解析 第13回 (7) 付録:線形問題が特性方程式で解ける仕組み
定数係数線形微分方程式 y″+ay′+by=0 ⋯⋯ (∗)
これを固有値・固有ベクトルの考え方で解きましょう。
- (∗) の解は
- 足してもスカラー倍しても解になる
- 2階微分方程式なので、一般解は2個の積分定数を含む
ことから、2次元のベクトル空間 V を成すことがわかります。
- 「微分をする」という写像を
D:V→V, Dy=y′
と書くことにすれば、D は V の一次変換になります。
- D の固有値・固有ベクトル λ, y とは
Dy=λy
を満たすスカラーと非零関数 y のことですから、
0=y″+ay′+b=D2y+aDy+by=(λ2+aλ+b)y
が成り立ちます。すなわち λ は特性方程式
λ2+aλ+b=0
の根、すなわち特性根です。
- λ に対する固有ベクトルは Dy=y′=λy の解ですから y=Aeλx です。
- 特性方程式が異なる特性根 λ, μ を持つときは、以上の議論から
y=Aeλx+Beμx
が (∗) の一般解になります。
行列は表に出てきませんが、定数係数線形微分方程式の解法にはこのように固有値・固有ベクトルの考え方が潜んでいます。
( 特性根が重根のときはもう少し議論が必要ですが。)
線形漸化式 xn+2+axn+1+bxn=0 ⋯⋯ (♯)
これも定数係数線形微分方程式と全く同じストーリーで解けます。
- (♯) の解は
- 足してもスカラー倍しても解になる
- x0, x1 を決めると解がただ一通りに決まるので、一般式は2個のパラメータを含む
ことから、2次元のベクトル空間 V を成すことがわかります。
- 「番号 n を 1 進める」という写像を
D:V→V, Dxn=xn+1
と書くことにすれば、D は V の一次変換になります。
- D の固有値・固有ベクトル λ, {xn} とは
Dxn=λxn
を満たすスカラーと非零数列 {xn} のことですから、
0=xn+2+axn+1+bxn=D2xn+aDxn+bxn=(λ2+aλ+b)xn
が成り立ちます。すなわち λ は特性方程式
λ2+aλ+b=0
の根、すなわち特性根です。
- λ に対する固有ベクトルは Dxn=xn+1=λxn の解ですから等比数列 xn=Aλn です。
- 特性方程式が異なる特性根 λ, μ を持つときは、以上の議論から
xn=Aλn+Bμn
が (♯) の一般解になります。