物理数学I 線形代数のソースを表示

<small> [[物理数学I]] > 線形代数</small>	
----
==線形代数==


===行列の定義と特別な行列===

====行列の定義====


数値を何らかの仕方で組み合わせたものを行列と呼ぶ。
ただし、縦の長さと横の長さを、各行と列でそろえなくてはならない。

例えば、
:<math>
\begin{pmatrix}
a _{11}&a _{12}\\
 a _{21}&a _{22}\\
\end{pmatrix}
</math>
は行列である。
高校までの範囲では、行列は3*3までしか扱わなかった。
しかし、実際には行列はm*n行列が存在し、(m.nは正の整数。)
全てにおいて和、積などの演算を行なうことが出来る。


====行列の和,積====


行列の和は各要素ごとに和を取ることによって定義される。
このことは、行列の和が可換であり、結合則を満たすことを保証する。
実数倍は、各要素に実数を書けることによって
定義する。この演算は行列に単位行列の定数倍ををかける演算と
等しいことに注意。n*n行列の単位行列はすぐ後に定義する。
これらの操作が可能なことを、行列の線形性と呼ぶ。


行列の積は、
:<math>
\begin{pmatrix}
a _{11}& \cdots & a _{1n} \\
&\cdots &\\
a _{n1}& \cdots & a _{nn} 
\end{pmatrix}
 \times 
\begin{pmatrix}
b _{11}& \cdots & b _{1n} \\
&\cdots &\\
b _{n1}& \cdots & b _{nn} 
\end{pmatrix}
=
\begin{pmatrix}
\sum  _i a _{1i} b _{i1} & \cdots & \sum  _i a _{1i} b _{in}\\
&\cdots &\\
\sum  _i a _{ni} b _{i1} & \cdots & \sum  _i a _{ni} b _{in}\\
\end{pmatrix}
</math> 
で与えられる。これらは2*2,3*3などの行列の演算の
拡張となっている。
この演算は短く
:<math>
(a b ) _{ij} = \Sigma  _k (a _{ik} b _{kj} )
</math>
と書かれることがある。
重要な事は、行列の積は
n*m行列とm*l行列のように(n,m,lは正の整数。)前の行列の列の数と後の行列の行の数が
一致していないと計算できないことが挙げられる。
このことからもわかるように、行列の積は一般に非可換である。
このことを式で書くと、
:<math>
AB \ne BA
</math>
となる。
 

====単位行列====

n*nのように行と列の数が等しい行列のことを正方行列と呼ぶ。
正方行列について、行と列の番号が等しい成分のことを対角成分と呼ぶ。
対角成分でない成分を非対角成分と呼ぶ。

n*n正方行列において、対角成分が全て1であり、非対角成分が
全て0である行列を、n*n単位行列と呼ぶ。
この行列は直接計算を行なうことで、
:<math>
EA = AE = A
</math>
を満たすことが分かる。(Aは任意の行列。)
このことは、短く
:<math>
\Sigma  _k \delta _{ik} a _{kj} = a _{ij}
</math>
と書かれることがある。
ここで、
:<math>
\delta _{ij}
</math>
は、ijが等しい値を持つとき1を返し、ijが等しくない値を持つとき0を返す行列であり
ちょうど単位行列に対応するものになっている。
この
:<math>
\delta
</math>
のことを歴史的理由によりクロネッカーのデルタと呼ぶ。


====複素共役行列====

ここまで、行列の要素がどのような数で与えられているかを指定していなかった。
実用的には多くの場合、行列の中身は実数か複素数のどちらかが与えられる。
とくに行列の複素共役を行列のうちで全ての要素の複素行列を取ったもので
与え、これを行列<math>A</math>に対して<math>A^*</math>と書く。(これらの記号は本によって違いがあるので
注意が必要である。場合によっては、<math>A^\dagger</math>、<math>\bar A</math>と書かれることが
ある。)
:<math>
(A^*)  _{ij} = a _{ij}^*
</math>
また、Aが実数だけで与えられる行列であるなら、その複素共役は
Aと一致する。


====転置行列====

転置行列とは、ある行列Aを
:<math>
A = a _{ij}
</math>
で与えるとき、
:<math>
A^t=a _{ji}
</math>
で与えられる。つまり、行列の行と列を入れ換えることで定義される
行列である。この記号も本によって異なるので、それぞれの本をよく読み比べることが
必要となる。

====随伴行列====

随伴行列はある行列Aに対して、
:<math>
A^\dagger = (A^t)^*=(A^*)^t
</math>
で与えられる。つまり、複素共役を取ってしかも転置した行列を
随伴行列と呼ぶのである。
ここでは、<math>\dagger</math>を用いたが、
この記号は前のものと同様、本によって統一されていないので注意が
必要となる。


====エルミート行列====

エルミート行列とは、その行列の随伴行列がその行列自身と一致する行列である。
エルミート行列は行列が実数で与えられているとき、転置が等しい行列に
帰着するが、この行列を対称行列と呼ぶ。
例えば、単位行列はエルミート行列である。

====ユニタリ行列====

ユニタリ行列(Uで表わされることが多い。)とは、その行列の随伴行列が元の行列の
逆行列となっている行列のことである。
:<math>
U U^{\dagger} = E
</math>
ユニタリ行列はもとの行列が実数で与えられているとき、
その行列の転置が、元の行列の逆行列に等しくなっている行列に
帰着するが、この行列を直交行列と呼ぶ。

(
TODO
ノルムを導入した場合の結論?
)

===逆行列の一般形===

====逆行列と一次方程式====


一次方程式は、


:<math>
a _{11}x _1 + \cdots a _{1n}x _n = b _1 
</math>

:<math>
a _{n1}x _1 + \cdots a _{nn}x _n = b _n 
</math>


で表わされる方程式である。(<math>a _i</math>, <math>b _i</math>は、定数。）
これらの一般解を求める。


上の連立方程式は、行列記法では

:<math>
Ax = b
</math>
と書ける。仮に、Aが逆行列を持つなら、
この式の一般解は、
:<math>
x = A^{-1} b
</math>
となる。
よって、1次方程式を解くのは、行列の逆行列を
求めることに等しい。

ここでは、逆行列を具体的に計算する方法を考察し、
その結果を用いて1次方程式を解く方法を考える。
そのために、いくらかの用語を導入する必要がある。

====行列式====

行列式は、
行列
:<math>
A = 
\begin{pmatrix}
a _{11} & \cdots & a _{1n} \\
\cdots & \cdots &\cdots\\
a _{n1} & \cdots & a _{nn} 
\end{pmatrix}
</math>
に対して、

:<math>
\sum _{\sigma} \textrm{sgn} \sigma a _{1\sigma _1} \cdot \cdot a _{n\sigma _n}
</math>
で与えられる数である。
これだけでは何を言っているか分からないかもしれないが、
順に説明を追っていって欲しい。

ここで、
:<math>
\sigma
</math>
は、1からnまでの整数の置換のうちのどれかを表わしている。

ここで整数の置換とは、ある整数の集合を取ったとき、それらが
互いに重複しないように
ある値をある値に対応させたものである。
例えば、整数の組1,2,3をとったとき
:<math>
1 \rightarrow 3, 2 \rightarrow 1, 3 \rightarrow 2
</math>
は1つの置換である。
1,2,3がそれぞれ1,2,3のうちの別の数に移っていることに
注目して欲しい。

置換の数はn個の整数の組を用いたときn!個ある。
例えば、
3個の整数の組では、
:<math>
1\rightarrow 1 , 2 \rightarrow 2 ,3\rightarrow 3
</math>,
:<math>
1\rightarrow 1 , 2 \rightarrow 3 ,3\rightarrow 2
</math>,
:<math>
1\rightarrow 2 , 2 \rightarrow 1 ,3\rightarrow 3
</math>,
:<math>
1\rightarrow 2 , 2 \rightarrow 3 ,3\rightarrow 1
</math>,
:<math>
1\rightarrow 3 , 2 \rightarrow 1 ,3\rightarrow 2
</math>,
:<math>
1\rightarrow 3 , 2 \rightarrow 2 ,3\rightarrow 1
</math>
の(6=3!)組となる。

上で挙げた6つの置換の中で
1番上の、全ての整数が変化しない
置換を、単位置換と呼ぶ。単位置換から偶数回だけの
変更を行なって得られる置換を隅置換、奇数回だけの
変更を行なって得られる置換を奇置換と呼ぶ。
ここで言う変更とは、置換の後にある整数が移る数を、別の整数が移る数と
互いに入れ換えることを言う。これが1度だけ起これば奇置換であり、
2回だけ起これば隅置換である。厳密には例えば、1と2という結果を何度も移し変えて
、単位置換自身が<math>2^n</math>(nは正の整数。)回置換に対応するということもできる。
しかし、この場合でも単位置換が隅置換である、という主張は変化していない
ことがわかる。(ただし0回は偶数なので0置換は隅置換であるとした。)
実はこの結果はより一般的なものであり、
ある置換について様々な入れ換えを行なっても、その置換が
隅置換か奇置換かということには影響しないことが知られている。
*TODO
導出?

つまり、結局のところある置換が隅置換か奇置換かどうかを決めたいときには、
最も簡単な仕方で1つ結果を得れば、全てのやり方で同じ結果が得られる
ことがわかる。

例として上に挙げた6つの置換について、それらが隅置換か奇置換かを判別する。

上でいう2,3,6番目の置換は、それぞれ
単位置換から1度だけ対応する値を交換して得られるので
奇置換である。(それぞれ2と3、1と2、3と1)を交換して
得られる。 
単位置換は0階の交換で得られるので隅置換であり
残った2つはそれぞれ2組の数値の入れ換えをすることで
得られるので、隅置換である。

行列式の定義の式
:<math>
\sum _{\sigma} \textrm{sgn} \sigma a _{1\sigma _1} \cdot \cdot a _{n\sigma _n}
</math>
で、
:<math>
\textrm{sgn} \sigma
</math>
は、
:<math>
\sigma
</math>
が隅置換であるとき、+1、奇置換であるとき-1となる。
つまり、単位置換をして得られる、
:<math>
a _{11} \cdot \cdot a _{nn}
</math>
の様な項に対しては、1をかけ算し、
ある奇置換
:<math>
a _{12}a _{21}a _{33}a _{44} \cdots \alpha _{nn}
</math>
については-1をかけ、
そのような値を全て足し合わせるということがこの計算の主旨となっている。

実際には、置換の数はn!の速さで増えるので、これらの項数は急激に増える。
そのため、n=2,3,4ぐらいのときを除いて、計算機を使わないで
値を得ることは困難になる。



*計算例

例えば、2次行列について計算する。このときは置換の数は
2つであるので、計算は簡単である。

2*2行列Aを
:<math>
A = 
\begin{pmatrix}
a&b\\
c&d
\end{pmatrix}
</math>
で与える。行列式の定義にしたがって計算すると、 
:<math>
\begin{matrix}
\textrm{det} A \\
=& a _{11}a _{22} - a _{12}a _{21}\\
=& a\times d - c \times b\\
=& ad -bc
\end{matrix}
</math>
となる。この値は高校では<math>\Delta</math>と呼ばれていた量であるが
この量は行列の行列式であったわけである。

3次の行列式では、
:<math>
\begin{matrix}
A = 
\begin{pmatrix}
a&b&c\\
d&e&f\\
g&h&i\\
\end{pmatrix}
\\
=a _{11}a _{22}a _{33}+
a _{12}a _{23}a _{31}+
a _{13}a _{21}a _{32}-
a _{11}a _{23}a _{32}-
a _{13}a _{22}a _{31}-
a _{12}a _{21}a _{33} \\
=aei + bfg + cdh - afh - bdi -ceg
\end{matrix}
</math>
となる。
これは、斜めに数を掛け合わせていったものに等しいことに注意。
例えば、第1項aeiは、1行1列のaから、3行3列のiまでを右下に向かって
順に書けていったものに等しい。また、次のbfgは、1行2列のbから始めて、
右下に向かってかけ算していったものに等しい。2行3列のfのあとは
端を突き抜けて、3行1列のgにいたることに注意。
4から6番目の項は、右下に向かってではなく左下に向かって
書けていった値となり同時にかけ算した値に(-1)をかける必要がある。
このような計算法を
サラスの公式と呼ぶことがある。

一方、4 * 4 以降の行列ではこのような簡単な計算法は
得られないので定義に従ってじょじょに計算していくことが求められる。


行列式は重要な性質として、
:<math>
\det AB = \det A \det B
</math>
という性質を持っている。
*TODO
(:導出?)

====小行列式====

ここでは、いきなり定義を与える。
:<math>
A = 
\begin{pmatrix}
a _{11} & \cdots & a _{1n} \\
\cdots & \cdots &\cdots \\
a _{n1} & \cdots & a _{nn} 
\end{pmatrix}
</math>
に対して、
m行とr列を除いて得られる(n-1)*(n-1)行列の行列式を、
行列Aのm行r列に関する小行列式と呼ぶ。

つまり、ある列と行を1本ずつ除いて、出来た行列の行列式を取ると
それが小行列式となっているわけである。

====行列式の展開====

行列式の計算を簡単に行なうため、
式の展開を導入する。この方法を用いると、
多くの行と列を持つ行列の行列式が比較的簡単に計算できる。
ただし、実際にはこの方法は基本行列による変形と組み合わせて使われることが
多いので、これだけを見てもそれほど便利とは感じないかも知れない。

*TODO
基本行列の導入
:永年行列式の計算のためにも...

行列式の展開とは、
:<math>
A = 
\begin{pmatrix}
a_{11} & \cdots & a_{1n} \\
\cdots & \cdots & \cdots \\
a_{n1} & \cdots & a_{nn} 
\end{pmatrix}
</math>
に対して、
:<math>
\det A = \sum  _{k = 1} ^n (-1) ^{m+k} a _{mk} b _{mk}
</math>,または
:<math>
\det A = \sum  _{k = 1} ^n (-1) ^{r+k} a _{kr} b _{kr}
</math>
が成り立つ。
ここで、
:<math>
b _{mk}
</math>
は、
m行k列に関する行列Aの小行列式である。

この式はn*nの行列の行列式を
n個の(n-1)*(n-1)行列の行列式の和によって表わしている。
このことから、この操作を行列式の展開と呼ぶのである。


*導出

例えば、n*nの行列式において
mk項を含む値は、
他に、m行または、k列に含まれる項を含んでいてはならない。
(これは、行列式に含まれる値がそれぞれn個の整数の置換であり、
その中でm行またはk列を表わす数は、1度しか含まれていないことによる。)
またそれ以外の項は、全体の置換が隅置換であったら前の符合が1になるように、
奇置換であったら前の符合が-1になるように計算されるが、
これはまさしく、m行とk列を除いた(n-1)*(n-1)行列の行列式、
すなわちm行k列に関する行列Aの小行列式 に
他ならない。
同様の考察を行列中の他の項についても繰りかえすと、
行列式の展開の式を得る。

例えば、3行3列の行列式の計算を行なうとき、
行列式の展開を使うと、
:<math>
\det 
\begin{pmatrix}
a&b&c\\
d&e&f\\
k&j&i
\end{pmatrix}
= a 
\det
\begin{pmatrix}
e&f\\
j&i
\end{pmatrix}
- b 
\det
\begin{pmatrix}
d&f\\
k&i
\end{pmatrix}
+c
\det
\begin{pmatrix}
d&e\\
k&j
\end{pmatrix}
</math>
となり、2次の行列の行列式の計算に帰着する。
それぞれの小行列式の前につく符合がi行j列であったら
<math>(-1)^{i+j}</math>に比例することに注意。
実際にはこのままだとあまり計算量が減っている感じがせず、
実際そうなのだが、特にある行または列にただ1つだけ0でない数が
入っており、残りは全て0という行や列が存在したとき、この行列式の
展開によってn*n行列はただ1つの(n-1)*(n-1)行列に帰着し、非常に
計算が楽になる。

====余因子行列====

ここでも無味乾燥な定義が続くが、がんばってもらいたい。
実はこの量は線形代数が終わるとほとんど出て来なくなるので
定義などしなくてもよさそうなものなのだが、とはいえこの量を
置かないと、次の逆行列の記述が非常に大変になるので、
ここでは導入するのである。

ある行列n*n行列Aに対して、 
:<math>
c _{ij} = (-1)^{i+j} b _{ji}
</math>
で定義される行列Cを行列Aの余因子行列と呼ぶ。
ここで<math>b _{ji} </math>は、
行列Aの行j列iに関する小行列式である。

例えば、2*2行列でこの値を計算してみる。
実際には3*3行列以降では、これを得るために
n*n個の行列式の計算をせねばならなくなるので、
計算はかなり大変になる。そのため、ここでは2*2行列を
扱うのである。

2*2行列では余因子行列も2*2なので、4つの小行列式を計算すればよい。
実際には2*2の行列式の小行列式は1*1の行列の行列式、というか
1*1の行列はただの数であるので
余因子行列は簡単に求められる。
:<math>
A =
\begin{pmatrix}
a&b\\
c&d
\end{pmatrix}
</math>
として実際に計算すると、
:<math>
c _{11} = b _{11}=a _{22} = d
</math>
:<math>
c _{12} = -b _{21}= -a _{12} = -b
</math>
:<math>
c _{21} = -b _{12}= -a _{21} = -c
</math>
:<math>
c _{22} = b _{22}= a _{11} = a
</math>
が得られる。よってAの余因子行列Cは、
:<math>
C = 
\begin{pmatrix}
d & -b\\
-c & a
\end{pmatrix}
</math>
となる。この形は実は2*2行列の逆行列の形と同じ形をしている。
つまり、この行列はAの逆行列に比例しているのである。
実はこれは一般的な結果で、ある行列の逆行列は
その行列の余因子行列を元の行列の行列式で割ったものになる。
ただし、Aの行列式が0に等しいときは例外であり、このとき行列A
は逆行列を持たない。
次のセクションではそのことの一般的な証明を
与える。 

====逆行列の一般型====

まずは、先ほどの結果をまとめる。

ある行列Aの逆行列は、 
:<math>
A^{-1} = \frac 1 {\det A} C
</math>
で書かれる。
ここでCは、Aの余因子行列である。
ただし、det A が0に等しいとき、行列Aは逆行列を持たない。

*導出

ここでは具体的に行列ACを計算し、
その結果が
:<math>
(\det A)  E
</math>
に等しいことを示す。(Eは単位行列。)
一般的に、第l行について考える。(l = 1 , ... , nとする。)
このとき、ACのll要素を考えると、
:<math>
\sum  _{m=1} ^ n a _{lm} c _{ml}
</math>
:<math>
=\sum  _{m=1} ^ n a _{lm} (-1)^{m+ l} b _{lm}
</math>,
(<math>b _{lm}</math>は、行列Aの行l、列mに関する小行列式。)
:<math>
=\det A
</math>
(式の展開の逆)
となり、ここでは確かに
:<math>
 (\det A)  E
</math>
と一致した結果になる。

次に、l行で、i列(i = 1, ... , n : l 以外) について
ACを考える。ここでは0になってくれるとよい。
:<math>
\sum  _{m=1} ^ n a _{lm} c _{mi}
</math>
:<math>
\sum  _{m=1} ^ n a _{lm} (-1)^{m+ i} b _{im}
</math>
これは、行列Aで、i行目をl行目で置き換えた行列の行列式に等しい。
行列式で行列のうちのある行か、ある列が他の行か他の列と一致する場合、
その2つの行または列からの寄与は必ず打ち消しあう。
よってi列からの寄与は0に等しい。
よって求める行列
ACは、
:<math>
(\det A ) E
</math>
となり、
:<math>
\frac 1 {\det A} C
</math>
は、(CはAの余因子行列)
Aの逆行列に等しいことが分る。

これによって、ある正方行列が与えられたとき、その行列の
逆行列を求める一般的な方法が得られたわけである。

しかし、実際にはこの計算は多くの計算量を必要とするので
実用的な計算には用いられない。
実用的な計算にはガウスの消去法が
用いられることが多い。
ガウスの消去法は計算機科学か数学の線形代数で扱われる。

====クラメルの公式====
{{main|クラメルの公式}}
ある一次方程式
:<math>
A \vec x = \vec b 
</math>
が与えられたとき、
Aの行列式が0でないとき、その解は
:<math>
b  _i = \frac 1 {\det A}
\begin{vmatrix}
a _{11} & \cdots & a _{1,i-1}&b _1 & a _{1,i+1}&\cdots &a _{1n}\\
&&&\cdots&&&\\
a _{n1} & \cdots & a _{n,i-1}&b _n & a _{n,i+1}&\cdots &a _{nn}
\end{vmatrix}
</math>
で与えられる。

*導出

:<math>
\vec x = A^{-1} \vec b
</math>
となる。
このとき、<math>A^{-1}</math>について
一般的な表式を用いると、
:<math>
\begin{matrix}
\textrm{rhs} =& \frac 1 {\det A} C \vec b\\
=& \frac 1 {\det A}\sum  _ {j} C  _{ij}  b _j\\
=& \frac 1 {\det A}\sum  _ {j} (-1)^{i+j} B  _{ji}  b _j\\
=& \frac 1 {\det A} |A _{\textrm{replaced} } |  _i (?)
\end{matrix}
</math>
ここで、
<math>B _{ij}</math>はAのi行j列に関する小行列式であり、<math>A _{\textrm{replaced} }</math>は、上で示した<math>i</math>番目の要素については行列<math>A</math>の<math>i</math>列目を列ベクトル<math> \vec b</math>で置き換えたものである。

これを用いて一般的に1次方程式の解が得られることが分る。

つまり、どんな1次方程式が現われても瞬時に答が得られるということである。

実際には計算が大変な事も多いので、この公式に頼らない方法も計算機科学の分野を中心に多く知られている。

一般的な方法としては基本行列を使った方法が良く知られているがいつその方法が書かれるかは今の時点では分からない。

==== determinant = 0 の幾何学的意味 ====

determinant が0であるとき、その行列には
逆行列が存在しない。
例えば、2次の行列A では、det A = ad - bc= 0  のとき、その行列には
逆行列が存在しない。

これは実際には、a/b = c/d に対応し、２つの直線が
平行である場合に対応する。


===行列の対角化===

====対角化の具体例====

ここではベクトルを幾何的なベクトルとして取ってみる。
つまり、
:<math>
\vec e _1 = \vec e _x=
\begin{pmatrix}
1\\
0
\end{pmatrix}
</math>
:<math>
\vec e _2 = \vec e _y=
\begin{pmatrix}
0\\
1
\end{pmatrix}
</math>
としてみる。このとき一般的な2*2行列は、あるベクトルを
別のベクトルに移す働きをするということが出来る。
例えば、
行列
:<math>
\begin{pmatrix}
a&b\\
c&d
\end{pmatrix}
</math>
は、
:<math>
e _1 = (a,c)
</math>
:<math>
e _2= (b,d)
</math>
の変更をする。

一般に多くの行列は、行列はもっともよい基底をとった場合、
対角化されてみえる。対角化とはどういうことかというと、
:<math>
A = \alpha  _i \delta  _{ij}
</math>
の形で書かれる行列のことである。
ただし、ここでiの和を取ってはいけない。
ここで<math>\delta _{ij}</math>はクロネッカのデルタであるので、
この式は、対角化された行列は対角要素だけを持つ行列であることが
わかる。つまり、ある行列はある基底をとることで
対角化されて見えるのである。

例えば
:<math>
A=
\begin{pmatrix}
1&1\\
1&1
\end{pmatrix}
</math>
のときを考えてみる。
この行列は
ベクトル(1 1)を2倍の長さにし、
(1 -1)を長さ0のベクトルにする。
全てのベクトルはこれら2つのベクトルを用いて張ることが出来るので
全てのベクトルに対するAの動作が判ったことになる。
例えば、
(2 0) = (1 1) + (1 -1)
を用いれば、
A(2 0) = A((1 1) + (1 -1)) = 2(1 1) + 0 (1 -1) = (2 2)
となる。
また、
<math>A^n</math> (2 0) = <math>2^n</math> (1 1) = (<math>2^n</math> <math>2^n</math>)
となる。
つまり、<math>e _1</math>,<math>e _2</math>という行列を取ったときには
その行列を互いに混ぜ合わせたり長さを変えたり更には向きを変えたりと
複雑なマッピングをしていたようにみえた行列が、非常に扱い易くなっている
のである。このように行列が最も扱い易くなる
ベクトルを探すことが、このセクションの主題である。

====対角化の方法====

====固有値と固有ベクトル====

行列Aについてあるベクトルxに対して
:<math>
Ax = \lambda x
</math>
が成り立つときこれを固有ベクトルと呼ぶ。
また<math>\lambda</math>を固有値と呼ぶ。
これは先ほどでいうと、行列の動作をただの定数倍にしてしまうベクトルのことである。
このようなベクトルを見つけられると非常に都合が良い。

この方程式は実際に解くことが出来る。
値は、
:<math>
(A - \lambda E) x=0
</math>
となる。
ここで、
:<math>
\lambda x = \lambda E x
</math>
となることを用いた。
(Eは単位行列。)
これは一見奇異に思えるかも知れないが、元々1次方程式が要素ごとの
連立方程式だったことをふまえて、要素ごとの記述に戻って
計算しても、同じ値が得られる。

さて、この方程式は常にx=0という解を持っている。
しかし、今必要なのはx = 0という解ではなく、0でないxである。
ここで、仮に
<math>A - \lambda E</math>
が、
:<math>
\det (A - \lambda E) \ne 0
</math>
を満たすとすると、この行列はただ1つの解を持ち、
その解はx=0である。
これはこのときには
<math>A - \lambda E</math>が逆行列を持つことからすぐにわかる。

よって、仮に行列Aがある0でない固有ベクトルを持つとすれば、
すくなくとも<math>\det (A- \lambda E) = 0</math> とならなければならない。
いいかえれば、そのことの必要条件は <math>\det (A- \lambda E) = 0</math> となる。
この方程式を固有方程式と呼ぶ。この方程式は
<math>\lambda</math>に関するn次の代数方程式であり、
nが大きいときには一般にとじた形で解けるとは限らない。
ただし、この式は必ず複素数の範囲でn個の解を持つことが
知られている。

==== 固有方程式の計算例 ====

固有方程式には様々な場合がある。
例として、全ての解が単根であるときや、いくつかの重根があるときが
ある。
一般に全ての解が単根であるときには、それぞれの固有値に対応する
固有ベクトルがただ1つずつ存在し、元の行列は完全に対角化
されることが知られている。このとき対角成分に
現われる数は、その固有ベクトルに対応する固有値である。
また、一般に
:<math>
A^\dagger A = A A^\dagger
</math>
を満たす行列は対角化されることが知られている。
重要な例としてユニタリ行列、エルミート行列、直交行列、対称行列などは
対角化可能である。実際物理などででてくるのは多くがこの場合であり、
特にエルミート行列の対角化は、量子力学の定式化において
重要な位置を占める。


いくつかの重根があるときには、行列は必ずしも対角化されるとは限らない。
例えば、
:<math>
\begin{pmatrix}
1&0\\
0&1
\end{pmatrix}
</math>
は対角化可能である、というより既に対角化されている。
しかし、
:<math>
\begin{pmatrix}
1&1\\
0&1
\end{pmatrix}
</math>
のような行列を対角化することは出来ない。
このような場合には一般にジョルダン標準型と呼ばれる
最も簡単な形に帰着させることが慣習的になっている。
ジョルダン標準型についてはおそらく数学の線形代数で
扱われるが、後にも少し扱う。
上の対角化できない行列は実際には既にジョルダン標準型の形になっている。

ここでは典型的な場合をいくつか解析してみる。
例えば、
:<math>
A =
\begin{pmatrix}
1&1\\
1&1
\end{pmatrix}
</math>
では、
固有方程式は
:<math>
(\lambda-1)^2 -1 = 0
</math>
:<math>
\lambda = 0,2
</math>
となる。
ここでこの固有方程式は全ての根が単根であるので、
これらの固有値にはそれぞれ1つずつの固有ベクトルが対応し、
それらによって行列が対角化されることが期待される。
実際、固有ベクトルが2本あるなら、その方向に関しては
その行列が良い振舞いをするという方向が2本分かっていることになり、
そして、ここでは2次元のベクトルを考えているので、
これらの2本だけで作り得る全てのベクトルを作ることが出来るので、
元の行列が全てのベクトルについて良い振舞いをすることはある意味で
当然であると思われる。

さて、得られた固有値を 
もともとの式
:<math>
(A -\lambda E )x = 0
</math>
または、
:<math>
A x = \lambda x
</math>
に代入することで固有ベクトルを求められる。
この操作は少しトリッキーに思えるかも知れないが、
こういうものだと思ってもらいたい。

ここでは、
<math>\lambda = 0</math>に対して、上の方程式は
:<math>
x _1 + x _2 = 0
</math>
となる。2つの方程式が同じ方程式になることに注意。
これは、この<math>\lambda</math>に対して
<math>A -\lambda E</math>の行列式が0になることからの帰結と見ることが
出来る。

この条件を満たすベクトルは全てこの固有値に関する固有ベクトルとなる。
ベクトルの長さは自由に選んで良い。
例えば、おそらく最も簡明なのは
:<math>
\begin{pmatrix}
1\\
-1
\end{pmatrix}
</math>
であり、これは<math>\lambda=0</math>に関する
固有ベクトルとなる。
実際には、ベクトルの長さを1にするために規格化を行ない、
:<math>
\frac 1 {\sqrt{2}}
\begin{pmatrix}
1\\
-1
\end{pmatrix}
</math>
としたものもよく用いられる。

<math>\lambda = 2</math> に対しては、例えば
:<math>
\begin{pmatrix}
1\\
1
\end{pmatrix}
</math>
が得られる。
これでこの行列は全てのベクトルに対して良い振舞いを
することが出来ることが分かった。


==== 対角化が出来ない場合 ====

先ほども述べた通り、
対角化は全ての行列について出来るとは限らない。
例えば、
:<math>
\begin{pmatrix}
1&0\\
0&1
\end{pmatrix}
</math>
は対角化可能であるが、
:<math>
\begin{pmatrix}
1&1\\
0&1
\end{pmatrix}
</math>
は対角化不可能である。
問題は、
固有値<math>\lambda</math>が重解を持ったとき
固有ベクトルの張るベクトルの次元が
充分な自由度を持たないことがありうるということにある。
この問題では、固有ベクトルをきめる
方程式は、<math>x_2 = 0 </math>に帰着するが、これを満たすベクトルは
(a 0) (aは任意の定数)であり、1次元でしか無い。
そのため、固有ベクトルが充分な数だけとれないということになっている。
このような場合は対角化は原理的に不可能なのである。
このときには、代わりにジョルダン標準型という形式が使われることが
多い。

===2次形式===

ここでは行列の対角化の応用として、2次形式の計算を扱う。 
実際の計算では2次式の計算を扱うことが良くある。
例えば、2次元の放物線型ポテンシャルを考えるとき、ある方向を取ることで
放物線が非常に簡単に見えることがある。そのようなことのために
行列の対角化の計算が使われるのである。

一般にn*n個の文字<math>x _1, \cdots x _n</math>を用いて書かれる2次式は、
:<math>
 {}^{t}{x} A x
</math>
の形で書ける。ただし、A,xともに実数しか取り得ないものとする。
ここで、便宜上Aは対称な行列とする。
例えば、
<math>x^2 + xy + y^2</math> に対しては、
:<math>
A = 
\begin{pmatrix}
1&1/2\\
1/2&1
\end{pmatrix}
</math>
と選べばよい。このように、うまく数値を選ぶことで、2次形式の
行列は必ず対称に選べる。

対称行列の対角化を用いると、
ある直交行列Oを用いて、 
:<math>
{}^t{x} A x = ^t {x} O^t B O x
</math>
と書ける。ここでBは対角行列である。
Ox = y と定義することで、
:<math>
{}^t{x} A x = \sum _{i=1}^n  \lambda _i y _i^2
</math>
が得られる。
右辺の形を2次形式の標準型と呼ぶ。


例えば、
２つの2次形式
<math>x^2 - y^2</math> と、<math>2xy</math>は、同一の標準型をもつ。
直交行列Oは、
:<math>
\begin{pmatrix}
\cos \pi/4&-\sin \pi/4\\
\sin \pi/4&\cos \pi/4
\end{pmatrix}
</math>
また、
<math>x^2 - y^2 = 1</math> のグラフは、
2xy = 1 のグラフを <math>\pi/4</math>回転させてみたものに
等しいこともこれに対応する。
結局<math>2xy</math>を見るときにはグラフを傾けてみるのが分かりやすく解析する
のに良い方法であり、そのよい方法は、実は行列の対角化と対応していた
ということである。

{{DEFAULTSORT:せんけいだいすう}}
[[カテゴリ:線形代数学]]