対角行列に似た行列
線型代数 において 、 正方行列は、 対角行列 に 相似で ある場合に、 対角化可能 または 非欠陥で ある と呼ばれます 。つまり、 となる 可逆行列 と対角行列が存在する場合です 。これは と同等です 。 (このような 、 は一意ではありません。) この特性は任意の線型写像に対して存在します。 有限次元 ベクトル空間に対して 、 線型写像は 、 の 固有ベクトル から なる の 順序付き基底 が存在する場合、 対角化可能で ある と呼ばれます 。これらの定義は同等です。 が上記のように 行列 表現 される場合 、 の列ベクトルは の固有ベクトルからなる基底を形成し 、 の対角要素 は の 対応する 固有値です 。 この固有ベクトル基底に関して、 は によって表されます 。
あ
{\displaystyle A}
ポ
{\displaystyle P}
だ
{\displaystyle D}
ポ
−
1
あ
ポ
=
だ
{\displaystyle P^{-1}AP=D}
あ
=
ポ
だ
ポ
−
1
{\displaystyle A=PDP^{-1}}
ポ
{\displaystyle P}
だ
{\displaystyle D}
五
{\displaystyle V}
T
:
五
→
五
{\displaystyle T:V\to V}
五
{\displaystyle V}
T
{\displaystyle T}
T
{\displaystyle T}
あ
=
ポ
だ
ポ
−
1
{\displaystyle A=PDP^{-1}}
ポ
{\displaystyle P}
T
{\displaystyle T}
だ
{\displaystyle D}
T
{\displaystyle T}
T
{\displaystyle T}
だ
{\displaystyle D}
対角化は 、上記の とを 見つけるプロセスであり 、これにより、その後の多くの計算が容易になります。 対角行列を べき乗するには、対角要素をそのべき乗するだけです。対角行列の 行列式は 、すべての対角要素の積です。このような計算は、 に簡単に一般化できます 。
ポ
{\displaystyle P}
だ
{\displaystyle D}
だ
{\displaystyle D}
あ
=
ポ
だ
ポ
−
1
{\displaystyle A=PDP^{-1}}
対角化可能な行列によって表される幾何学的変換は、 非同次膨張 (または 異方性スケーリング )です。つまり、異なる方向に異なる量で空間を スケーリング できます。各固有ベクトルの方向は、対応する固有値によって指定された係数によってスケーリングされます。
対角化できない正方行列は、 欠陥行列 と呼ばれます。実数 要素を持つ 行列が 実数上で欠陥行列になることがあります。これは、実数要素を持つ任意の 逆行列 および対角行列では不可能であるが、 複素数 要素 では可能であるため、 複素数上で対角化可能であることを意味します。たとえば、一般的な 回転行列 がこれに該当します。
あ
{\displaystyle A}
あ
=
ポ
だ
ポ
−
1
{\displaystyle A=PDP^{-1}}
ポ
{\displaystyle P}
だ
{\displaystyle D}
あ
{\displaystyle A}
対角化可能行列の結果の多くは、 代数的に閉じた体 (複素数など)上でのみ成り立ちます。この場合、対角化可能行列はすべての行列の空間で 稠密であり、これは、欠陥のある行列を小さな 摂動 で対角化可能行列に変形できることを意味します 。また、 ジョルダン・シュヴァレー分解は、任意の 行列が対角化可能行列とべき 零行列の一意の和であることを示します。代数的に閉じた体上では、対角化可能行列は 半単純行列 と同等です 。
意味
体 に要素を持つ 正方 行列は、逆行列(つまり、 一般線型群 GL n ( F ) の元) が存在し、 が 対角行列で
ある場合 、 対角化可能 または 非欠陥である と呼ばれます。
ん
×
ん
{\displaystyle n\times n}
あ
{\displaystyle A}
ふ
{\displaystyle F}
ん
×
ん
{\displaystyle n\times n}
ポ
{\displaystyle P}
ポ
−
1
あ
ポ
{\displaystyle P^{-1}AP}
特徴づけ
対角化可能な写像と行列に関する基本的な事実は次のように表現されます。
体上の 行列 が 対角化可能であるのは、その固有空間の 次元 の合計 が に等しい 場合であり 、これはの固有ベクトルからなる の 基底 が存在する場合のみです。そのような基底が見つかった場合は、 これらの 基底ベクトルを 列として持つ行列 を形成することができ 、 は の固有値を対角要素とする対角行列になります 。この行列は の モード行列 として知られています 。
ん
×
ん
{\displaystyle n\times n}
あ
{\displaystyle A}
ふ
{\displaystyle F}
ん
{\displaystyle n}
ふ
ん
{\displaystyle F^{n}}
あ
{\displaystyle A}
ポ
{\displaystyle P}
ポ
−
1
あ
ポ
{\displaystyle P^{-1}AP}
あ
{\displaystyle A}
ポ
{\displaystyle P}
あ
{\displaystyle A}
線型写像は、その固有空間の 次元 の合計が に等しい 場合のみ対角化可能であり 、 これはの固有ベクトルからなる の基底が存在する場合のみ当てはまります 。このような基底に関して、 は対角行列によって表されます。この行列の対角要素は の固有値です 。
T
:
五
→
五
{\displaystyle T:V\to V}
薄暗い
(
五
)
{\displaystyle \dim(V)}
五
{\displaystyle V}
T
{\displaystyle T}
T
{\displaystyle T}
T
{\displaystyle T}
次の十分条件(ただし必要条件ではない)は、多くの場合役立ちます。
行列 が 体上で対角化可能である のは、行列が において異なる固有値を持つ場合 、 つまり、その 特性多項式が において異なる根を 持つ場合です 。 ただし、その逆は偽である可能性があります。 が固有値 1、2、2 (すべて異なるわけではない) を持ち、対角形式 ( と 同様 ) と 基底行列の変更 によって対角化可能であることを考えます。が 1 より大きい次元の固有空間を持つ 場合、その逆は失敗します。この例では、 固有値 2 に関連付けられた の固有空間は次元 2 です。
ん
×
ん
{\displaystyle n\times n}
あ
{\displaystyle A}
ふ
{\displaystyle F}
ん
{\displaystyle n}
ふ
{\displaystyle F}
ん
{\displaystyle n}
ふ
{\displaystyle F}
[
−
1
3
−
1
−
3
5
−
1
−
3
3
1
]
、
{\displaystyle {\begin{bmatrix}-1&3&-1\\-3&5&-1\\-3&3&1\end{bmatrix}},}
あ
{\displaystyle A}
[
1
0
0
0
2
0
0
0
2
]
{\displaystyle {\begin{bmatrix}1&0&0\\0&2&0\\0&0&2\end{bmatrix}}}
ポ
{\displaystyle P}
[
1
1
−
1
1
1
0
1
0
3
]
。
{\displaystyle {\begin{bmatrix}1&1&-1\\1&1&0\\1&0&3\end{bmatrix}}.}
あ
{\displaystyle A}
あ
{\displaystyle A}
を持つ 線型写像は、異なる固有値 を持つ場合 、つまりその特性多項式が において異なる根を持つ場合に対角化可能です 。
T
:
五
→
五
{\displaystyle T:V\to V}
ん
=
薄暗い
(
五
)
{\displaystyle n=\dim(V)}
ん
{\displaystyle n}
ん
{\displaystyle n}
ふ
{\displaystyle F}
を 上の行列と します 。 が対角化可能な場合 、 の任意のべき乗も対角化可能です。逆に、 が可逆で、 代数的に閉じており、 の特性の整数倍でない に対して対角化可能な 場合 、 は対角化可能です。証明: が 対角化可能な場合、 は ある多項式 によって消滅します が、 は重根を持たず ( であるため ) 、 の最小多項式で割り切れます 。
あ
{\displaystyle A}
ふ
{\displaystyle F}
あ
{\displaystyle A}
あ
{\displaystyle A}
ふ
{\displaystyle F}
あ
ん
{\displaystyle A^{n}}
ん
{\displaystyle n}
ふ
{\displaystyle F}
あ
{\displaystyle A}
あ
ん
{\displaystyle A^{n}}
あ
{\displaystyle A}
(
x
ん
−
λ
1
)
⋯
(
x
ん
−
λ
け
)
{\displaystyle \left(x^{n}-\lambda _{1}\right)\cdots \left(x^{n}-\lambda _{k}\right)}
λ
j
≠
0
{\displaystyle \lambda _{j}\neq 0}
A
{\displaystyle A}
複素数 上では、ほぼすべての行列が対角化可能です。より正確には、 上で対角化でき ない 複素行列 の集合は 、 の サブセット として考えられ 、 ルベーグ測度は ゼロです。また、対角化可能な行列は、 ザリスキー位相 に関して稠密なサブセットを形成するとも言えます 。つまり、対角化できない行列は、 特性多項式の 判別式の 消失集合 の内部にあり、これは 超曲面です。このことから、 ノルム によって与えられる通常の( 強い )位相における密度も導かれます 。 上では同じことは当てはまりません 。
C
{\displaystyle \mathbb {C} }
n
×
n
{\displaystyle n\times n}
C
{\displaystyle \mathbb {C} }
C
n
×
n
{\displaystyle \mathbb {C} ^{n\times n}}
R
{\displaystyle \mathbb {R} }
ジョルダン ・シュヴァレー分解は、 演算子を半単純 (つまり、対角化可能な) 部分とべき 零 部分の合計として表現します。したがって、行列は、べき零部分がゼロの場合にのみ対角化可能です。言い換えると、 ジョルダン形式 の各ブロックにべき零部分がない場合は、行列は対角化可能です。つまり、各「ブロック」は 1 行 1 列の行列です。
対角化
次の 2 つの任意の基底 とを考えます。 基底 E に関して記述される
行列によって表される線形変換が存在するとします。また、次の固有方程式が存在するとします。
E
=
{
e
i
|
∀
i
∈
[
n
]
}
{\displaystyle E=\{{{\boldsymbol {e}}_{i}|\forall i\in [n]}\}}
F
=
{
α
i
|
∀
i
∈
[
n
]
}
{\displaystyle F=\{{{\boldsymbol {\alpha }}_{i}|\forall i\in [n]}\}}
A
E
{\displaystyle A_{E}}
A
E
α
E
,
i
=
λ
i
α
E
,
i
{\displaystyle A_{E}{\boldsymbol {\alpha }}_{E,i}=\lambda _{i}{\boldsymbol {\alpha }}_{E,i}}
アルファ固有ベクトルは、E 基底に関しても記述されます。集合 F は行列 A の固有ベクトルの集合であり、任意のベクトル空間にまたがるため、 に類似した対角行列である行列が存在すると言えます 。言い換えると、 行列 が基底 F で記述されている場合、 は対角化可能な行列です。基底の変更計算は 、次のように基底を E から F に変更する遷移行列 を使用して実行します。
D
F
{\displaystyle D_{F}}
A
E
{\displaystyle A_{E}}
A
E
{\displaystyle A_{E}}
S
{\displaystyle S}
D
F
=
S
E
F
A
E
S
E
−
1
F
{\displaystyle D_{F}=S_{E}^{F}\ A_{E}\ S_{E}^{-1F}}
、
ここで、 はE 基底から F 基底への遷移行列です。逆行列は、 基底を F から E に変更する
新しい遷移行列と等しくすることができ、次の関係が得られます。
S
E
F
{\displaystyle S_{E}^{F}}
P
{\displaystyle P}
S
E
−
1
F
=
P
F
E
{\displaystyle S_{E}^{-1F}=P_{F}^{E}}
遷移行列と遷移行列は 両方とも 可逆です。したがって、次の方法で行列を操作できます。 行列は と表され 、これは依然として E 基底にあります。同様に、対角行列は F 基底にあります。
S
{\displaystyle S}
P
{\displaystyle P}
D
=
S
A
E
S
−
1
D
=
P
−
1
A
E
P
{\displaystyle {\begin{aligned}D=S\ A_{E}\ S^{-1}\\D=P^{-1}\ A_{E}\ P\end{aligned}}}
A
E
{\displaystyle A_{E}}
A
{\displaystyle A}
対称行列の対角化は、軸を回転させて固有ベクトルに合わせることとして解釈できます。
行列が 対角化できる場合、つまり、
A
{\displaystyle A}
P
−
1
A
P
=
[
λ
1
0
⋯
0
0
λ
2
⋯
0
⋮
⋮
⋱
⋮
0
0
⋯
λ
n
]
=
D
,
{\displaystyle P^{-1}AP={\begin{bmatrix}\lambda _{1}&0&\cdots &0\\0&\lambda _{2}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &\lambda _{n}\end{bmatrix}}=D,}
それから:
A
P
=
P
[
λ
1
0
⋯
0
0
λ
2
⋯
0
⋮
⋮
⋱
⋮
0
0
⋯
λ
n
]
.
{\displaystyle AP=P{\begin{bmatrix}\lambda _{1}&0&\cdots &0\\0&\lambda _{2}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &\lambda _{n}\end{bmatrix}}.}
遷移行列 S には、基底 F に書き込まれた列として E 基底ベクトルがあります。逆に、逆遷移行列 P には、 基底 E に書き込まれた F 基底ベクトルがあるため、次のようにして P をブロック行列形式で表すことができます。
α
i
{\displaystyle {\boldsymbol {\alpha }}_{i}}
P
=
[
α
E
,
1
α
E
,
2
⋯
α
E
,
n
]
,
{\displaystyle P={\begin{bmatrix}{\boldsymbol {\alpha }}_{E,1}&{\boldsymbol {\alpha }}_{E,2}&\cdots &{\boldsymbol {\alpha }}_{E,n}\end{bmatrix}},}
結果として、次のように書くことができます。
A
[
α
E
,
1
α
E
,
2
⋯
α
E
,
n
]
=
[
α
E
,
1
α
E
,
2
⋯
α
E
,
n
]
D
.
{\displaystyle {\begin{aligned}A{\begin{bmatrix}{\boldsymbol {\alpha }}_{E,1}&{\boldsymbol {\alpha }}_{E,2}&\cdots &{\boldsymbol {\alpha }}_{E,n}\end{bmatrix}}={\begin{bmatrix}{\boldsymbol {\alpha }}_{E,1}&{\boldsymbol {\alpha }}_{E,2}&\cdots &{\boldsymbol {\alpha }}_{E,n}\end{bmatrix}}D.\end{aligned}}}
ブロック行列形式では、A 行列は 1x1 次元の行列であるのに対し、P は 1xn 次元の行列であると考えることができます。D 行列は、すべての対角要素を含む nxn 次元の行列として完全な形式で記述できます。
A
[
α
E
,
1
α
E
,
2
⋯
α
E
,
n
]
=
[
α
E
,
1
α
E
,
2
⋯
α
E
,
n
]
[
λ
1
0
⋯
0
0
λ
2
⋯
0
⋮
⋮
⋱
⋮
0
0
⋯
λ
n
]
.
{\displaystyle A{\begin{bmatrix}{\boldsymbol {\alpha }}_{E,1}&{\boldsymbol {\alpha }}_{E,2}&\cdots &{\boldsymbol {\alpha }}_{E,n}\end{bmatrix}}={\begin{bmatrix}{\boldsymbol {\alpha }}_{E,1}&{\boldsymbol {\alpha }}_{E,2}&\cdots &{\boldsymbol {\alpha }}_{E,n}\end{bmatrix}}{\begin{bmatrix}\lambda _{1}&0&\cdots &0\\0&\lambda _{2}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &\lambda _{n}\end{bmatrix}}.}
上記の行列乗算を実行すると、次の結果が得られます。 ブロック行列の各要素を両側で個別に取ると、次のようになります。
A
[
α
1
α
2
⋯
α
n
]
=
[
λ
1
α
1
λ
2
α
2
⋯
λ
n
α
n
]
{\displaystyle {\begin{aligned}A{\begin{bmatrix}{\boldsymbol {\alpha }}_{1}&{\boldsymbol {\alpha }}_{2}&\cdots &{\boldsymbol {\alpha }}_{n}\end{bmatrix}}={\begin{bmatrix}\lambda _{1}{\boldsymbol {\alpha }}_{1}&\lambda _{2}{\boldsymbol {\alpha }}_{2}&\cdots &\lambda _{n}{\boldsymbol {\alpha }}_{n}\end{bmatrix}}\end{aligned}}}
A
α
i
=
λ
i
α
i
(
i
=
1
,
2
,
…
,
n
)
.
{\displaystyle A{\boldsymbol {\alpha }}_{i}=\lambda _{i}{\boldsymbol {\alpha }}_{i}\qquad (i=1,2,\dots ,n).}
したがって、 の列ベクトルはの 右固有ベクトル で あり 、 対応する対角要素は対応する 固有値 です。 の可逆性は、 固有ベクトルが 線形独立で あり、 の基底を形成することも示唆しています 。 これは、対角化可能性と対角化の標準的なアプローチの必要かつ十分な条件です。 の 行ベクトル はの 左固有ベクトル です 。
P
{\displaystyle P}
A
{\displaystyle A}
P
{\displaystyle P}
F
n
{\displaystyle F^{n}}
P
−
1
{\displaystyle P^{-1}}
A
{\displaystyle A}
複素行列が エルミート行列 (またはより一般的には 正規行列 ) である場合、 の固有ベクトルは の 正規直交基底 を形成するように選択でき 、 は ユニタリ行列 になるように選択できます 。 さらに、 が 実 対称行列 である場合、その固有ベクトルは の正規直交基底になるように選択でき 、 は 直交行列 になるように選択できます 。
A
∈
C
n
×
n
{\displaystyle A\in \mathbb {C} ^{n\times n}}
A
{\displaystyle A}
C
n
{\displaystyle \mathbb {C} ^{n}}
P
{\displaystyle P}
A
∈
R
n
×
n
{\displaystyle A\in \mathbb {R} ^{n\times n}}
R
n
{\displaystyle \mathbb {R} ^{n}}
P
{\displaystyle P}
実際の作業のほとんどでは、行列はコンピューター ソフトウェアを使用して数値的に対角化されます。これを実現するための アルゴリズムは多数 存在します。
同時対角化
行列の集合が同時対角化可能 であるとは、集合内の 任意の行列に対して が対角行列となる ような 可逆行列が 1 つだけ存在する場合を言う 。同時対角化可能な行列を特徴付ける定理は次の通りである: 対角化可能な 行列の集合が可換と なるのは、集合が同時対角化可能である場合のみである。 [1] : p. 64
P
{\displaystyle P}
P
−
1
A
P
{\displaystyle P^{-1}AP}
A
{\displaystyle A}
対角化可能な行列( 上 ) 全体の集合は 同時に 対角化可能ではない。例えば、行列
n
×
n
{\displaystyle n\times n}
C
{\displaystyle \mathbb {C} }
n
>
1
{\displaystyle n>1}
[
1
0
0
0
]
and
[
1
1
0
0
]
{\displaystyle {\begin{bmatrix}1&0\\0&0\end{bmatrix}}\quad {\text{and}}\quad {\begin{bmatrix}1&1\\0&0\end{bmatrix}}}
は対角化可能ですが、可換ではないため同時に対角化することはできません。
集合が可換な 正規行列から構成される場合、かつその集合が ユニタリ行列 によって同時に対角化可能である場合に 限ります。つまり、 集合内の
任意の に対して が対角となる ようなユニタリ行列が存在します。
U
{\displaystyle U}
U
∗
A
U
{\displaystyle U^{*}AU}
A
{\displaystyle A}
リー理論 の言語では 、同時に対角化可能な行列の集合は トーラル・リー代数 を生成します。
例
対角化可能な行列
反転は 実数(および特性が 2 でない任意の体)上で対角化可能であり、対角線上には ±1 があります。
有限位数の自己準 同型は 、対角線上に単位根を持つ (または体の特性が自己準同型の位数を割り切れない任意の代数的に閉じた体) 上 で対角化可能である。これは、単位根が異なるため、 最小多項式が 分離可能であることから導かれる。
C
{\displaystyle \mathbb {C} }
投影は 対角化可能であり、対角線上に 0 と 1 が配置されます。
実 対称行列は 直交行列 で対角化できます 。つまり、実対称行列 が与えられた 場合 、 何らかの直交行列 に対して は対角行列になります 。 より一般的には、行列 が ユニタリ行列で対角化できるのは、それらが 正規 で ある場合のみです。実対称行列の場合、 となるため 、 明らかにが成り立ちます。正規行列の例には、実対称 (または 歪対称 ) 行列 (共分散行列など) や エルミート行列 (または歪エルミート行列) があります。 無限次元ベクトル空間への一般化については、 スペクトル定理を参照してください。
A
{\displaystyle A}
Q
T
A
Q
{\displaystyle Q^{\mathrm {T} }AQ}
Q
{\displaystyle Q}
A
=
A
T
{\displaystyle A=A^{\mathrm {T} }}
A
A
T
=
A
T
A
{\displaystyle AA^{\mathrm {T} }=A^{\mathrm {T} }A}
対角化できない行列
一般に、 回転行列は 実数体上では対角化できませんが、すべての 回転行列は 複素数体上では対角化できます。行列が対角化できない場合でも、常に「最善を尽くして」、主対角線上に固有値があり、上対角線上に 1 または 0 がある同じ特性を持つ行列を見つけることができます。これは、 ジョルダン正規形 として知られています。
いくつかの行列は、どの体でも対角化できません。最も顕著なのは、非零冪 零行列です。これは、より一般的には、固有値の 代数的重複度と幾何的 重複度が一致しない場合に発生します 。たとえば、
C
=
[
0
1
0
0
]
.
{\displaystyle C={\begin{bmatrix}0&1\\0&0\end{bmatrix}}.}
この行列は対角化できません。つまり、が 対角行列となる ような行列は存在しません。実際、 には 1 つの固有値 (つまり 0) があり、この固有値の代数的重複度は 2、幾何学的重複度は 1 です。
U
{\displaystyle U}
U
−
1
C
U
{\displaystyle U^{-1}CU}
C
{\displaystyle C}
実数行列の中には、実数上で対角化できないものがある。例えば、次の行列を考える。
B
=
[
0
1
−
1
0
]
.
{\displaystyle B=\left[{\begin{array}{rr}0&1\\\!-1&0\end{array}}\right].}
この行列には 実固有値がないので、 対角行列と なるような 実行 列は存在しません。しかし、 複素数を許せば対角化することができます。実際、
B
{\displaystyle B}
Q
{\displaystyle Q}
Q
−
1
B
Q
{\displaystyle Q^{-1}BQ}
B
{\displaystyle B}
Q
=
[
1
i
i
1
]
,
{\displaystyle Q={\begin{bmatrix}1&i\\i&1\end{bmatrix}},}
は対角行列です。反 時計回りに角度だけ回転する回転行列である ことは簡単にわかります。
Q
−
1
B
Q
{\displaystyle Q^{-1}BQ}
B
{\displaystyle B}
θ
=
−
π
2
{\textstyle \theta =-{\frac {\pi }{2}}}
上記の例は、対角化可能な行列の合計が対角化可能である必要がないことを示していることに注意してください。
行列を対角化する方法
行列を対角化するということは、固有ベクトルが基底を形成する場合、
その 固有値と固有ベクトルを 求めるのと同じプロセスです。たとえば、次の行列を考えます。
A
=
[
0
1
−
2
0
1
0
1
−
1
3
]
.
{\displaystyle A=\left[{\begin{array}{rrr}0&1&\!\!\!-2\\0&1&0\\1&\!\!\!-1&3\end{array}}\right].}
特性多項式 の根は 固有値 です 。 この線形方程式を解くと、 固有ベクトル と が得られ 、 は と なります 。 つまり、 の 場合 となります 。 これらのベクトルは の基底を形成するため 、 基底変換 行列 の列ベクトルとして組み立てると 、次の式が得られます。
この方程式を変換で表すと、次のようになります。 は 標準基底 を固有基底 に変換するため 、 次の式が得られます。
したがって、 は 標準基底 を固有ベクトルとして持ち、これが の定義特性となります 。
p
(
λ
)
=
det
(
λ
I
−
A
)
{\displaystyle p(\lambda )=\det(\lambda I-A)}
λ
1
=
1
,
λ
2
=
1
,
λ
3
=
2
{\displaystyle \lambda _{1}=1,\lambda _{2}=1,\lambda _{3}=2}
(
I
−
A
)
v
=
0
{\displaystyle \left(I-A\right)\mathbf {v} =\mathbf {0} }
v
1
=
(
1
,
1
,
0
)
{\displaystyle \mathbf {v} _{1}=(1,1,0)}
v
2
=
(
0
,
2
,
1
)
{\displaystyle \mathbf {v} _{2}=(0,2,1)}
(
2
I
−
A
)
v
=
0
{\displaystyle \left(2I-A\right)\mathbf {v} =\mathbf {0} }
v
3
=
(
1
,
0
,
−
1
)
{\displaystyle \mathbf {v} _{3}=(1,0,-1)}
A
v
i
=
λ
i
v
i
{\displaystyle A\mathbf {v} _{i}=\lambda _{i}\mathbf {v} _{i}}
i
=
1
,
2
,
3
{\displaystyle i=1,2,3}
V
=
R
3
{\displaystyle V=\mathbb {R} ^{3}}
P
{\displaystyle P}
P
−
1
A
P
=
[
1
0
1
1
2
0
0
1
−
1
]
−
1
[
0
1
−
2
0
1
0
1
−
1
3
]
[
1
0
1
1
2
0
0
1
−
1
]
=
[
1
0
0
0
1
0
0
0
2
]
=
D
.
{\displaystyle P^{-1}AP=\left[{\begin{array}{rrr}1&0&1\\1&2&0\\0&1&\!\!\!\!-1\end{array}}\right]^{-1}\left[{\begin{array}{rrr}0&1&\!\!\!-2\\0&1&0\\1&\!\!\!-1&3\end{array}}\right]\left[{\begin{array}{rrr}1&\,0&1\\1&2&0\\0&1&\!\!\!\!-1\end{array}}\right]={\begin{bmatrix}1&0&0\\0&1&0\\0&0&2\end{bmatrix}}=D.}
P
{\displaystyle P}
P
e
i
=
v
i
{\displaystyle P\mathbf {e} _{i}=\mathbf {v} _{i}}
P
−
1
A
P
e
i
=
P
−
1
A
v
i
=
P
−
1
(
λ
i
v
i
)
=
λ
i
e
i
,
{\displaystyle P^{-1}AP\mathbf {e} _{i}=P^{-1}A\mathbf {v} _{i}=P^{-1}(\lambda _{i}\mathbf {v} _{i})=\lambda _{i}\mathbf {e} _{i},}
P
−
1
A
P
{\displaystyle P^{-1}AP}
D
{\displaystyle D}
における固有ベクトルの優先順序は存在しないことに注意する 。における 固有ベクトル の順序を変更すると、 の対角化された形式における 固有値 の順序が変わるだけである 。 [ 2]
P
{\displaystyle P}
P
{\displaystyle P}
A
{\displaystyle A}
行列関数への応用
対角化は行列の累乗を効率的に計算するために使用できます 。
A
=
P
D
P
−
1
{\displaystyle A=PDP^{-1}}
A
k
=
(
P
D
P
−
1
)
k
=
(
P
D
P
−
1
)
(
P
D
P
−
1
)
⋯
(
P
D
P
−
1
)
=
P
D
(
P
−
1
P
)
D
(
P
−
1
P
)
⋯
(
P
−
1
P
)
D
P
−
1
=
P
D
k
P
−
1
,
{\displaystyle {\begin{aligned}A^{k}&=\left(PDP^{-1}\right)^{k}=\left(PDP^{-1}\right)\left(PDP^{-1}\right)\cdots \left(PDP^{-1}\right)\\&=PD\left(P^{-1}P\right)D\left(P^{-1}P\right)\cdots \left(P^{-1}P\right)DP^{-1}=PD^{k}P^{-1},\end{aligned}}}
後者は対角行列の累乗のみを扱うので計算が簡単です。例えば、 上の例の
固有値を持つ行列の場合、次のように計算します。
A
{\displaystyle A}
λ
=
1
,
1
,
2
{\displaystyle \lambda =1,1,2}
A
k
=
P
D
k
P
−
1
=
[
1
0
1
1
2
0
0
1
−
1
]
[
1
k
0
0
0
1
k
0
0
0
2
k
]
[
1
0
1
1
2
0
0
1
−
1
]
−
1
=
[
2
−
2
k
−
1
+
2
k
2
−
2
k
+
1
0
1
0
−
1
+
2
k
1
−
2
k
−
1
+
2
k
+
1
]
.
{\displaystyle {\begin{aligned}A^{k}=PD^{k}P^{-1}&=\left[{\begin{array}{rrr}1&\,0&1\\1&2&0\\0&1&\!\!\!\!-1\end{array}}\right]{\begin{bmatrix}1^{k}&0&0\\0&1^{k}&0\\0&0&2^{k}\end{bmatrix}}\left[{\begin{array}{rrr}1&\,0&1\\1&2&0\\0&1&\!\!\!\!-1\end{array}}\right]^{-1}\\[1em]&={\begin{bmatrix}2-2^{k}&-1+2^{k}&2-2^{k+1}\\0&1&0\\-1+2^{k}&1-2^{k}&-1+2^{k+1}\end{bmatrix}}.\end{aligned}}}
このアプローチは、行列指数関数 や、べき級数として定義できる他の 行列関数 に一般化できます 。たとえば、 を定義すると 、 次のようになります。
exp
(
A
)
=
I
+
A
+
1
2
!
A
2
+
1
3
!
A
3
+
⋯
{\textstyle \exp(A)=I+A+{\frac {1}{2!}}A^{2}+{\frac {1}{3!}}A^{3}+\cdots }
exp
(
A
)
=
P
exp
(
D
)
P
−
1
=
[
1
0
1
1
2
0
0
1
−
1
]
[
e
1
0
0
0
e
1
0
0
0
e
2
]
[
1
0
1
1
2
0
0
1
−
1
]
−
1
=
[
2
e
−
e
2
−
e
+
e
2
2
e
−
2
e
2
0
e
0
−
e
+
e
2
e
−
e
2
−
e
+
2
e
2
]
.
{\displaystyle {\begin{aligned}\exp(A)=P\exp(D)P^{-1}&=\left[{\begin{array}{rrr}1&\,0&1\\1&2&0\\0&1&\!\!\!\!-1\end{array}}\right]{\begin{bmatrix}e^{1}&0&0\\0&e^{1}&0\\0&0&e^{2}\end{bmatrix}}\left[{\begin{array}{rrr}1&\,0&1\\1&2&0\\0&1&\!\!\!\!-1\end{array}}\right]^{-1}\\[1em]&={\begin{bmatrix}2e-e^{2}&-e+e^{2}&2e-2e^{2}\\0&e&0\\-e+e^{2}&e-e^{2}&-e+2e^{2}\end{bmatrix}}.\end{aligned}}}
これは、フィボナッチ数列 などの 線形再帰シーケンス の項の閉じた形式の表現を見つけるときに特に役立ちます 。
特定の用途
たとえば、次の行列を考えてみましょう。
M
=
[
a
b
−
a
0
b
]
.
{\displaystyle M={\begin{bmatrix}a&b-a\\0&b\end{bmatrix}}.}
のさまざまな累乗を計算すると、 驚くべきパターンが明らかになります。
M
{\displaystyle M}
M
2
=
[
a
2
b
2
−
a
2
0
b
2
]
,
M
3
=
[
a
3
b
3
−
a
3
0
b
3
]
,
M
4
=
[
a
4
b
4
−
a
4
0
b
4
]
,
…
{\displaystyle M^{2}={\begin{bmatrix}a^{2}&b^{2}-a^{2}\\0&b^{2}\end{bmatrix}},\quad M^{3}={\begin{bmatrix}a^{3}&b^{3}-a^{3}\\0&b^{3}\end{bmatrix}},\quad M^{4}={\begin{bmatrix}a^{4}&b^{4}-a^{4}\\0&b^{4}\end{bmatrix}},\quad \ldots }
上記の現象は を対角化することで説明できます 。 これを実現するには、の固有ベクトルからなる の基底が必要です 。 そのような固有ベクトル基底の1つは次のように与えられます。
M
{\displaystyle M}
R
2
{\displaystyle \mathbb {R} ^{2}}
M
{\displaystyle M}
u
=
[
1
0
]
=
e
1
,
v
=
[
1
1
]
=
e
1
+
e
2
,
{\displaystyle \mathbf {u} ={\begin{bmatrix}1\\0\end{bmatrix}}=\mathbf {e} _{1},\quad \mathbf {v} ={\begin{bmatrix}1\\1\end{bmatrix}}=\mathbf {e} _{1}+\mathbf {e} _{2},}
ここで e i は R n の標準基底を表す 。基底の逆変換は次のように表される。
e
1
=
u
,
e
2
=
v
−
u
.
{\displaystyle \mathbf {e} _{1}=\mathbf {u} ,\qquad \mathbf {e} _{2}=\mathbf {v} -\mathbf {u} .}
単純な計算でわかるのは
M
u
=
a
u
,
M
v
=
b
v
.
{\displaystyle M\mathbf {u} =a\mathbf {u} ,\qquad M\mathbf {v} =b\mathbf {v} .}
したがって、 a と bはそれぞれ u と v に対応する固有値です 。行列の乗算の線形性により、次の式が得られます。
M
n
u
=
a
n
u
,
M
n
v
=
b
n
v
.
{\displaystyle M^{n}\mathbf {u} =a^{n}\mathbf {u} ,\qquad M^{n}\mathbf {v} =b^{n}\mathbf {v} .}
標準基準に戻ると、
M
n
e
1
=
M
n
u
=
a
n
e
1
,
M
n
e
2
=
M
n
(
v
−
u
)
=
b
n
v
−
a
n
u
=
(
b
n
−
a
n
)
e
1
+
b
n
e
2
.
{\displaystyle {\begin{aligned}M^{n}\mathbf {e} _{1}&=M^{n}\mathbf {u} =a^{n}\mathbf {e} _{1},\\M^{n}\mathbf {e} _{2}&=M^{n}\left(\mathbf {v} -\mathbf {u} \right)=b^{n}\mathbf {v} -a^{n}\mathbf {u} =\left(b^{n}-a^{n}\right)\mathbf {e} _{1}+b^{n}\mathbf {e} _{2}.\end{aligned}}}
前述の関係を行列形式で表すと、
M
n
=
[
a
n
b
n
−
a
n
0
b
n
]
,
{\displaystyle M^{n}={\begin{bmatrix}a^{n}&b^{n}-a^{n}\\0&b^{n}\end{bmatrix}},}
それによって上記の現象が説明されます。
量子力学の応用
量子力学 および 量子化学 計算において、 行列の対角化は最も頻繁に適用される数値処理の 1 つです。基本的な理由は、無限次元 ヒルベルト空間上のほとんどの物理的状況ではあるが、時間に依存しない シュレーディンガー方程式 が固有値方程式であるためです 。
非常に一般的な近似は、ヒルベルト空間を有限次元に切り捨てる(または投影する)ことです。その後、シュレーディンガー方程式を実対称行列または複素エルミート行列の固有値問題として定式化できます。正式には、この近似は 変分原理 に基づいており、下から制限されるハミルトニアンに有効です。
一次摂動論は 縮退状態の行列固有値問題にもつながります。
参照
注記
参考文献
^ Horn, Roger A.; Johnson, Charles R. (2013). マトリックス分析、第2版 。ケンブリッジ大学出版局 。ISBN 9780521839402 。
^ Anton, H.; Rorres, C. (2000 年 2 月 22 日)。 初等線形代数 (応用版) (第 8 版)。John Wiley & Sons。ISBN 978-0-471-17052-5 。