logo

正方行列のシューア分解 📂行列代数

正方行列のシューア分解

定義

あるユニタリー行列 QQ上三角行列 TT について、A=QTQA = Q T Q^{\ast} ならばAAシュア分解schur Factorizationを持つという。

定理

すべての正方行列 ACm×mA \in \mathbb{C}^{ m \times m} はシュア分解を持つ。

説明

固有値対角化 の欠点は、A=SΛS1A = S \Lambda S^{-1} で分解された場合でもなおS1S^{-1} を見つける労が必要であることだ。冪乗の計算時間が劇的に減るのは事実だが、行列代数のすべての問題が逆行列を求めることに帰結すると考えれば、冪乗が非常に多く必要な問題以外では、治療よりも害が大きくなる可能性がある。

一方、シュア分解はTTが上三角行列として与えられるため、冪乗には少し役に立たないかもしれないが、QQ^{\ast} を見つけることが非常に簡単であるため、TT を適切に扱えば固有値対角化よりもずっと速く問題を解決する可能性がある。上三角行列自体がそんなに扱いづらい行列ではないことを考えると、その汎用性はかなりのものと言えるだろう。そして何よりも、まずシュア分解は正方行列であるだけで十分なほど条件が緩い。

証明

有限次元 mm に対して数学的帰納法を使って証明する。AmA_{m} をあるm×m m \times m 行列とする。

m=1m = 1 ならばA1A_{1}スカラーであり、自明である。

m2m \ge 2 に対してAm1A_{m-1}がシュア分解を持つと仮定する。

x\mathbf{x}を固有値λ\lambdaに対応する単位固有ベクトルとし、U:=[xU1]U := \begin{bmatrix} \mathbf{x} & U_{1} \end{bmatrix}ユニタリー行列と定義すればxAmx=λ\mathbf{x}^{\ast} A_{m} \mathbf{x} = \lambdaであり、U1x=0U_{1}^{\ast} \mathbf{x} = 0

UAmU=[xU1]Am[xU1]=[xAmxxAmU1U1AmxU1AmU1]=[λxAmU1U1λxAm1]=[λB0Am1] \begin{align*} U^{\ast} A_{m} U =& \begin{bmatrix} \mathbf{x}^{\ast} \\ U_{1}^{\ast} \end{bmatrix} A_{m} \begin{bmatrix} \mathbf{x} & U_{1} \end{bmatrix} \\ =& \begin{bmatrix} \mathbf{x}^{\ast} A_{m} \mathbf{x} & \mathbf{x}^{\ast} A_{m} U_{1} \\ U_{1}^{\ast} A_{m} \mathbf{x} & U_{1}^{\ast} A_{m} U_{1} \end{bmatrix} \\ =& \begin{bmatrix} \lambda & \mathbf{x}^{\ast} A_{m} U_{1} \\ U_{1}^{\ast} \lambda \mathbf{x} & A_{m-1} \end{bmatrix} \\ =& \begin{bmatrix} \lambda & B \\ \mathbb{0} & A_{m-1} \end{bmatrix} \end{align*}

BCm1B \in \mathbb{C}^{m-1}は任意のベクトルであり、仮定によりAm1A_{m-1}はシュア分解を持つ。

したがって、Am1A_{m-1}はあるユニタリー行列 VV上三角行列 TT に対してAm1=VTVA_{m-1} = V T V^{\ast}として表される。

それならば、ユニタリー行列 Q:=[100V]Q := \begin{bmatrix} 1 & 0 \\ 0 & V \end{bmatrix}上三角行列 [λBV0T]\begin{bmatrix} \lambda & B V \\ 0 & T \end{bmatrix}に対して

Q[λBV0T]Q=[100V][λBV0T][100V]=[100V][λB0TV]=[λB0VTV]=[λB0Am1] \begin{align*} Q \begin{bmatrix} \lambda & B V \\ 0 & T \end{bmatrix} Q^{\ast} &= \begin{bmatrix} 1 & 0 \\ 0 & V \end{bmatrix} \begin{bmatrix} \lambda & B V \\ 0 & T \end{bmatrix} \begin{bmatrix} 1 & 0 \\ 0 & V^{\ast} \end{bmatrix} \\ =& \begin{bmatrix} 1 & 0 \\ 0 & V \end{bmatrix} \begin{bmatrix} \lambda & B \\ 0 & T V^{\ast} \end{bmatrix} \\ =& \begin{bmatrix} \lambda & B \\ 0 & V T V^{\ast} \end{bmatrix} \\ =& \begin{bmatrix} \lambda & B \\ 0 & A_{m-1} \end{bmatrix} \end{align*}

すなわち

Am=UQ[λBV0T]QU=(UQ)[λBV0T](UQ) \begin{align*} A_{m} =& U Q \begin{bmatrix} \lambda & B V \\ 0 & T \end{bmatrix} Q^{\ast} U^{\ast} \\ =& ( U Q ) \begin{bmatrix} \lambda & B V \\ 0 & T \end{bmatrix} ( U Q )^{\ast} \end{align*}

Am1A_{m-1}がシュア分解を持てばAmA_{m}もシュア分解を持つので、数学的帰納法により、すべての正方行列はシュア分解を持つ。