Processing math: 100%

数値解析 第13回 (7) 付録:線形問題が特性方程式で解ける仕組み

定数係数線形微分方程式  y+ay+by=0 ()

 これを固有値・固有ベクトルの考え方で解きましょう。
  1. () の解は
    • 足してもスカラー倍しても解になる
    • 2階微分方程式なので、一般解は2個の積分定数を含む
    ことから、2次元のベクトル空間 V を成すことがわかります。
  2. 「微分をする」という写像を
    D:VV,  Dy=y
    と書くことにすれば、DV の一次変換になります。
  3. D の固有値・固有ベクトル λ, y とは
    Dy=λy
    を満たすスカラーと非零関数 y のことですから、
    0=y+ay+b=D2y+aDy+by=(λ2+aλ+b)y
    が成り立ちます。すなわち λ は特性方程式
    λ2+aλ+b=0
    の根、すなわち特性根です。
  4. λ に対する固有ベクトルは Dy=y=λy の解ですから y=Aeλx です。
  5. 特性方程式が異なる特性根 λ, μ を持つときは、以上の議論から
    y=Aeλx+Beμx
    () の一般解になります。
行列は表に出てきませんが、定数係数線形微分方程式の解法にはこのように固有値・固有ベクトルの考え方が潜んでいます。 ( 特性根が重根のときはもう少し議論が必要ですが。)

線形漸化式  xn+2+axn+1+bxn=0 ()

 これも定数係数線形微分方程式と全く同じストーリーで解けます。
  1. () の解は
    • 足してもスカラー倍しても解になる
    • x0, x1 を決めると解がただ一通りに決まるので、一般式は2個のパラメータを含む
    ことから、2次元のベクトル空間 V を成すことがわかります。
  2. 「番号 n1 進める」という写像を
    D:VV,  Dxn=xn+1
    と書くことにすれば、DV の一次変換になります。
  3. D の固有値・固有ベクトル λ, {xn} とは
    Dxn=λxn
    を満たすスカラーと非零数列 {xn} のことですから、
    0=xn+2+axn+1+bxn=D2xn+aDxn+bxn=(λ2+aλ+b)xn
    が成り立ちます。すなわち λ は特性方程式
    λ2+aλ+b=0
    の根、すなわち特性根です。
  4. λ に対する固有ベクトルは Dxn=xn+1=λxn の解ですから等比数列 xn=Aλn です。
  5. 特性方程式が異なる特性根 λ, μ を持つときは、以上の議論から
    xn=Aλn+Bμn
    () の一般解になります。