重根と複素数根(multiple roots and complex roots)

前節では特性方程式が異なる実数根をもつ場合を考えました．ここでは重根と複素数根をもつ場合を考えます．もし複素数$a + bi$が微分方程式 ${\textbf X}^{\prime} = A{\textbf X}$の特性方程式の根ならば, ${\textbf X}^{\prime} = A{\textbf X}$は $c{\textbf Q}e^{(a+bi)t}$の形をした解をもっているはずです．ただし, ${\textbf Q}$は複素固有ベクトル,$c$は任意の複素定数です．

${\textbf X}^{\prime} = A{\textbf X}$の$A$が実行列のとき,複素指数や複素ベクトルの存在は扱いにくいものです．とくに初期値 ${\textbf X}(0) = {\textbf X}_{0}$が実数のときはなおさらです．そこで解を実数と実関数だけで表わす必要があります．

$A$を実行列, $\lambda = a+bi$が複素固有値で ${\textbf C} = {\textbf C}_{1} + i{\textbf C}_{2} \ ({\textbf C}_{1},{\textbf C}_{2}$   実数ベクトル$)$が$\lambda$に対する固有ベクトルだとすると固有方程式(eigenequation)

$$(A - \lambda I){\textbf C} = {\textbf 0}$$

よりその共役も

$$\overline{(A - \lambda I){\textbf C}} = (\bar{A} - \bar{\lambda}\bar{I})\bar{\textbf C} = {\textbf 0}$$

を満たします．よって $\bar{\lambda} = a - bi$も固有値となり,対応する固有ベクトル $\bar{\textbf C}$をもちます．したがって, ${\textbf X}_{1} = {\textbf C}e^{\lambda t}$と ${\textbf X}_{2} = \bar{\textbf C}e^{\bar{\lambda} t}$はともに ${\textbf X}^{\prime} = A{\textbf X}$の解になります．もちろん ${\textbf X}_{1}$と ${\textbf X}_{2}$の一次結合も解です．とくにその中でも

$$\Re{\textbf C}e^{\lambda t} = \Re{\textbf X}_{1} = \frac{{\textbf... ...tbf X}_{2}}{2} = {\textbf C}_{1}e^{at}\cos{bt} - {\textbf C}_{2}e^{at}\sin{bt}$$

と

$$\Im{\textbf C}e^{\lambda t} = \Im{\textbf X}_{1} = \frac{{\textbf... ...bf X}_{2}}{2i} = {\textbf C}_{1}e^{at}\sin{bt} - {\textbf C}_{2}e^{at}\cos{bt}$$

は一次独立な実数値解です．ここにでてきた $\Re{\textbf C}e^{\lambda t}$と $\Im{\textbf C}e^{\lambda t}$のことを,複素数 ${\textbf C}e^{\lambda t}$の実部(real part)と虚部(imaginary part)といいます．

例題 3..3   ${\textbf X}^{\prime} = \left(\begin{array}{rrr} 2&-1&-1\\ 2&1&-1\\ 0&-1&1 \end{array}\right){\textbf X}$を解け．

解

$$\det(A - \lambda I) = \det\left(\begin{array}{rrr} 2-\lambda&-1&-... ...&-1&1-\lambda \end{array}\right) = (2-\lambda)(\lambda^{2} - 2\lambda + 2) = 0$$

より固有値 $\lambda = 2, 1 \pm i$を得る．

まず,固有値$2$に対応する固有ベクトルを求めると

$$A - 2I = \left(\begin{array}{rrr} 0&-1&-1\\ 2&-1&-1\\ 0&-1&-1 \end{array... ...rrow} \left(\begin{array}{rrr} 1&-1/2&-1/2\\ 0&1&1\\ 0&0&0 \end{array}\right)$$

$$\stackrel{\frac{R_{2}}{2}+ R_{1}\to R_{1}} {\longrightarrow} \left(\begin{array}{rrr} 1&0&0\\ 0&1&1\\ 0&0&0 \end{array}\right)$$

より固有ベクトルは ${\textbf C} = \left(\begin{array}{r} 0\\ -1\\ 1 \end{array}\right)$である．

次に固有値$1+i$に対応する固有ベクトルを求める．

$$A - (1+i) I = \left(\begin{array}{rrr} 1-i&-1&-1\\ 2&-i&-1\\ 0&-1&-i \end{array}\right)$$

$$\stackrel{\begin{array}{c} {}_{\frac{1}{1-i} \times R_{1}\to R_{1}}\\ {}_{-(1+i)R_{1} + R_{2}\to R_{2}} \end{array}}{\longrightarrow} \left(\begin{array}{rrr} 1&-1/1-i&-1/1-i\\ 0&1&i\\ 0&-1&-i \end{array}\right)$$

$$\stackrel{\begin{array}{c} {}_{\frac{R_{2}}{1-i} + R_{1}\to R_{1}}\\ {}_{R_{2} + R_{3}\to R_{3}} \end{array}}{\longrightarrow} \left(\begin{array}{rrr} 1&0&-1\\ 0&1&i\\ 0&0&0 \end{array}\right)$$

より固有ベクトルは ${\textbf C} = \left(\begin{array}{r} 1\\ -i\\ 1 \end{array}\right)$である．これより ${\textbf C}e^{\lambda t}$の実部と虚部を求めると,

$${\textbf C}e^{\lambda t} = \left(\begin{array}{c} 1\\ -i\\ 1 \end{array}\right)e^{(1+i)t} = \left(\begin{array}{c} 1\\ -i\\ 1 \end{array}\right)e^{t}(\cos{t} + i\sin{t})$$

$$= \left(\begin{array}{c} e^{t}\cos{t} + ie^{t}\sin{t}\\ e^{t}\sin{... ...}\\ - e^{t}\cos{t}\\ e^{t}\sin{t} \end{array}\right)}_{\mbox{imaginary part}}$$

となる．よって一般解は次のように表わすことができる．

$${\textbf X}(t) = c_{1}\left(\begin{array}{c} 0\\ e^{2t}\\ -e^{2... ... e^{t}\cos{t}\\ e^{t}\sin{t} \end{array}\right). \ensuremath{\ \blacksquare}$$

さて,行列$A$の特性方程式が重根をもつ場合はどうでしょう．この場合も,$A$が対角化可能(つまり,重複度と同数の一次独立な固有ベクトルをもつ)ならば,単根のときと同じように一般解を求めることができます．では$A$が対角化不可能の場合はどうでしょう．こんな場合に便利な方法として行列の指数関数を用いる方法があります．

${\textbf X}^{\prime} = A{\textbf X}$の解として, $y = ce^{at}$が $y^{\prime} = ay$の解であったように, ${\textbf X} = {\textbf C}e^{At}$を考えます．ただし$A$は$n$次の正方行列, ${\textbf C}$は定数ベクトルです．これが可能なためには$e^{At}$に意味づけを行なう必要があります．そこで$e^{at}$のMaclaurin級数

$$e^{at} = 1 + at + \frac{(at)^{2}}{2!} + \frac{(at)^{3}}{3!} + \cdots$$

より

$$e^{At} = I + tA + \frac{t^{2}}{2!}A^{2} + \frac{t^{3}}{3!}A^{3} + \cdots$$

で定義します．この展開はすべての$t$で収束することが知られています．これを行列$A$の指数行列(exponential matrix)といいます．また

$$\frac{d}{dt}e^{At} = A + tA^{2} + \frac{t^{2}}{2!}A^{3} + \cdots = Ae^{At}$$

より, ${\textbf X} = {\textbf C}e^{At}$は ${\textbf X}^{\prime} = A{\textbf X}$の解といえます．

ここでひとつ問題があります．それは,どうやって ${\textbf C}e^{At}$を計算するかということです．定義を使うと無限級数の計算をすることになりますが,じつは無限級数の計算をせずに,解を求める方法があるのです．まず ${\textbf C}e^{(A - \lambda I)t}$の展開を考えて下さい．つまり,

$${\textbf C}e^{(A - \lambda I)t} = {\textbf C} + t(A - \lambda I){\textbf C} + \frac{t^{2}}{2!}(A - \lambda I)^{2}{\textbf C}$$

$$+ \frac{t^{3}}{3!}(A - \lambda I)^{3}{\textbf C} + \cdots + \frac{t^{k}}{k!}(A - \lambda I)^{k}{\textbf C} + \cdots$$

ここである$k$で $(A - \lambda I)^{k}{\textbf C} = 0$が成り立てば,上の展開で$k$の後はすべて0になります．よってこのときは,有限級数の計算をすればよいことになります．

例題 3..4   $\left\{\begin{array}{rl} x_{1}^{\prime} =& 4x_{1} - x_{2} + 2x_{3}\\ x_{2}^{\prime} =& 4x_{2} + 2x_{3}\\ x_{3}^{\prime} =& 4x_{3} \end{array}\right .$を解け．

解

$$\det(A - \lambda I) = \det\left(\begin{array}{rrr} 4-\lambda&-1&2\\ 0&4-\lambda&2\\ 0&0&4-\lambda \end{array}\right) = (4 - \lambda)^{3}$$

より,固有値 $\lambda = 4$である．次に固有値$4$に対する固有ベクトルCを求める．

$$A - 4I = \left(\begin{array}{rrr} 0&-1&2\\ 0&0&2\\ 0&0&0 \end{a... ...ghtarrow} \left(\begin{array}{rrr} 0&1&0\\ 0&0&1\\ 0&0&0 \end{array}\right)$$

より ${\textbf C} = \left(\begin{array}{c} 1\\ 0\\ 0 \end{array}\right)$となり,ひとつの解

$$\left(\begin{array}{c} 1\\ 0\\ 0 \end{array}\right)e^{4t}$$

を得る．

$A$は三次の正方行列なので一般解を表わすには,$3$つの一次独立な解が必要である．そこで

$$(A - 4I)^{2}{\textbf C} = {\textbf 0},\ (A - 4I){\textbf C} \neq {\textbf 0}$$

を満たす ${\textbf C}$を見つける．

$$(A - 4I)^{2} = \left(\begin{array}{rrr} 0&0&-2\\ 0&0&0\\ 0&0&0 \end{array}\right)$$

より $c_{1},c_{2}$は任意の値を取ることができるので, ${\textbf C} = \left(\begin{array}{c} \alpha\\ \beta\\ 0 \end{array}\right)$とおくと, $(A - 4I)^{2}{\textbf C} = {\textbf 0}$を満たす．次にもうひとつの条件 $(A - 4I){\textbf C} \neq {\textbf 0}$を満たすように $\alpha ,\beta$を選ぶと,

$$(A - 4I){\textbf C} = \left(\begin{array}{rrr} 0&-1&2\\ 0&0&2\\ ... ...end{array}\right) = \left(\begin{array}{c} -\beta\\ 0\\ 0 \end{array}\right)$$

より$\beta = 1$．よって ${\textbf C} = \left(\begin{array}{c} 0\\ 1\\ 0 \end{array}\right)$となり,2つめの解

$$e^{At}{\textbf C} = e^{4t}e^{(A - 4I)t}{\textbf C} = e^{4t}[{\textbf C} + t(A - 4I){\textbf C}]$$

$$= e^{4t}\{\left(\begin{array}{c} 0\\ 1\\ 0 \end{array}\right) + t... ...0&0&0 \end{array}\right)\left(\begin{array}{c} 0\\ 1\\ 0 \end{array}\right)\}$$

$$= e^{4t}\{\left(\begin{array}{c} 1\\ 1\\ 0 \end{array}\right) + \... ...array}\right)\} = e^{4t}\left(\begin{array}{c} 1-t\\ 1\\ 0 \end{array}\right)$$

を得る．3つめの解を求めるために,

$$(A - 4I)^{3}{\textbf C} = {\textbf 0},\ (A - 4I)^{2}{\textbf C} \neq {\textbf 0}$$

を満たす ${\textbf C}$を見つける． $(A - 4I)^{3}$は零行列なので ${\textbf C} = \left(\begin{array}{c} \alpha\\ \beta\\ \gamma \end{array}\right)$は $(A - 4I)^{3}{\textbf C} = {\textbf 0}$を満たす．ここでもうひとつの条件を用い,前に選んだ ${\textbf C}$と一次独立になるように, $\alpha,\beta,\gamma$を選ぶと, ${\textbf C} = \left(\begin{array}{c} 0\\ 0\\ 1 \end{array}\right)$ となる．これより3つめの解を計算すると

$$e^{At}{\textbf C} = e^{4t}e^{(A - 4I)t}{\textbf C}$$

$$= e^{4t}[{\textbf C} + t(A - 4I){\textbf C} + \frac{t^{2}}{2!}(A - 4I)^{2}{\textbf C}]$$

$$= e^{4t}[\left(\begin{array}{c} 0\\ 0\\ 1 \end{array}\right) + t\... ... 0&0&0 \end{array}\right)\left(\begin{array}{c} 0\\ 0\\ 1 \end{array}\right)]$$

$$= e^{4t}\left(\begin{array}{c} 2t - t^{2}\\ 2t\\ 1 \end{array}\right)$$

である．ここで $\left(\begin{array}{c} 1\\ 0\\ 0 \end{array}\right), \left(\begin{array}{c} ... ...array}\right), \left(\begin{array}{c} 2t - t^{2}\\ 2t\\ 1 \end{array}\right)$は一次独立なので,一般解は

$${\textbf X} = e^{4t}[c_{1}\left(\begin{array}{c} 1\\ 0\\ 0 \end... ...right) + c_{3}\left(\begin{array}{c} 2t - t^{2}\\ 2t\\ 1 \end{array}\right)]$$

で与えられる． $\blacksquare$