ここでは固有値と固有ベクトルを解説します。行列の基本事項は当たり前のように現れるので、自信のない場合は、本シリーズの中学でもわかる⑬「行列」から読み始めるとよいでしょう。最後にケイリー・ハミルトンの定理について解説します。 まずは、特に必要な「逆行列」と、「行列を用いた連立1次方程式の解法」を復習します(中学でも分かる⑭「逆行列」)。
<逆行列の定義>(復習) 行列 A = ( a b c d ) A=\begin{pmatrix} a & b \\ c & d \end{pmatrix} A = ( a c b d ) について
∣ A ∣ = a d − b c ≠ 0 |A|=ad-bc\ne0 ∣ A ∣ = a d − b c = 0 のとき、A A A には逆行列 A − 1 A^{-1} A − 1 が存在して
A − 1 = 1 ∣ A ∣ ( d − b − c a ) = 1 a d − b c ( d − b − c a ) \begin{align*}
A^{-1}&=\dfrac{1}{|A|}\begin{pmatrix*}[r] d & -b \\ -c & a \end{pmatrix*}\\[10pt]
&=\dfrac{1}{ad-bc}\begin{pmatrix*}[r] d & -b \\ -c & a \end{pmatrix*}
\end{align*} A − 1 = ∣ A ∣ 1 ( d − c − b a ) = a d − b c 1 ( d − c − b a )
∣ A ∣ = a d − b c = 0 |A|=ad-bc=0 ∣ A ∣ = a d − b c = 0 のとき、A A A には逆行列は存在しない。
***
∣ A ∣ |A| ∣ A ∣ を行列式といいます。下図のように「たすきがけ」を実行します。
行列式は、det ( A ) \det(A) det ( A ) と表すこともあります。つまり
det ( A ) = a d − b c \begin{align*}
\det(A)=ad-bc
\end{align*} det ( A ) = a d − b c
”det \det det ” は、行列式を意味する英単語 determinant(デターミナント)に由来します。 また、行列式は成分を明示して ∣ a b c d ∣ \begin{vmatrix} a & b \\ c & d \end{vmatrix} a c b d と表すこともあります。つまり
∣ a b c d ∣ = a d − b c \begin{align*}
\begin{vmatrix} a & b \\ c & d \end{vmatrix}=ad-bc
\end{align*} a c b d = a d − b c
具体的な成分計算を行う場面では、成分が明示されている方が分かりやすいこともあります。
次に、行列を用いた連立1次方程式の解法を復習します。A X = 0 AX=0 A X = 0 の形をした連立方程式です。
行列を用いた連立1次方程式の解法(復習) 復習として、次の例題を与えます。後ほど固有値、固有ベクトルを求めるときにこの解法を用います。
例題1 次の式を満たす ( x y ) \begin{pmatrix} x \\ y \end{pmatrix} ( x y ) を求めなさい。 (1) ( 2 − 1 4 3 ) ( x y ) = ( 0 0 ) \begin{pmatrix*}[r] 2 & -1 \\ 4 & 3 \end{pmatrix*}\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix}\raisebox{3ex}{ } ( 2 4 − 1 3 ) ( x y ) = ( 0 0 ) (2) ( 2 1 4 2 ) ( x y ) = ( 0 0 ) \begin{pmatrix*}[c] 2 & 1 \\ 4 & 2 \end{pmatrix*}\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix}\raisebox{3ex}{ } ( 2 4 1 2 ) ( x y ) = ( 0 0 )
指針 (1) 2次正方行列 ( 2 − 1 4 3 ) \begin{pmatrix*}[r] 2 & -1 \\ 4 & 3 \end{pmatrix*} ( 2 4 − 1 3 ) に、逆行列があるかどうかを確認します。 そのためには、この行列の行列式を計算して、(case1) 行列式が 0 0 0 でなければ、逆行列が存在し、その逆行列を左● か● ら● 掛けます。(case2) 行列式が 0 0 0 ならば、逆行列は存在しないので、そのまま行列の連立1次方程式を計算します。
(2) も同様です。具体的に解いていきましょう。
解答 (1) A = ( 2 − 1 4 3 ) A=\begin{pmatrix*}[r] 2 & -1 \\ 4 & 3 \end{pmatrix*} A = ( 2 4 − 1 3 ) とおくと、与式は、
A ( x y ) = ( 0 0 ) A\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix} A ( x y ) = ( 0 0 )
A A A の行列式は ∣ A ∣ = 2 ⋅ 3 − ( − 1 ) ⋅ 4 = 6 − ( − 4 ) = 6 + 4 = 10 ≠ 0 |A|=2\cdot3-(-1)\cdot4=6-(-4)=6+4=10\ne0 ∣ A ∣ = 2 ⋅ 3 − ( − 1 ) ⋅ 4 = 6 − ( − 4 ) = 6 + 4 = 10 = 0 より、A A A には逆行列 A − 1 A^{-1} A − 1 が存在して
A − 1 = 1 ∣ A ∣ ( 3 1 − 4 2 ) = 1 10 ( 3 1 − 4 2 ) \begin{align*}
A^{-1}&=\dfrac{1}{|A|}\begin{pmatrix*}[r] 3 & 1 \\ -4 & 2 \end{pmatrix*}\\
&=\dfrac{1}{10}\begin{pmatrix*}[r] 3 & 1 \\ -4 & 2 \end{pmatrix*}
\end{align*} A − 1 = ∣ A ∣ 1 ( 3 − 4 1 2 ) = 10 1 ( 3 − 4 1 2 )
与式 A ( x y ) = ( 0 0 ) A\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix} A ( x y ) = ( 0 0 ) の両辺に、左● か● ら● A − 1 A^{-1} A − 1 を掛けて
A − 1 A ( x y ) = A − 1 ( 0 0 ) \begin{align*}
A^{-1}A\begin{pmatrix} x \\ y \end{pmatrix}=A^{-1}\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} A − 1 A ( x y ) = A − 1 ( 0 0 )
行列の計算では交換法則が成り立たないので、左から掛けるか右から掛けるかをあらかじめ決めて、掛ける方向を統一します(中学でもわかる⑬「行列」、⑭「逆行列」)。
A − 1 A = E = ( 1 0 0 1 ) A^{-1}A=E=\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} A − 1 A = E = ( 1 0 0 1 ) より
E ( x y ) = A − 1 ( 0 0 ) \begin{align*}
E\begin{pmatrix} x \\ y \end{pmatrix}=A^{-1}\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} E ( x y ) = A − 1 ( 0 0 )
単位行列 E E E は掛けても変わらないので
( x y ) = A − 1 ( 0 0 ) \begin{align*}
\begin{pmatrix} x \\ y \end{pmatrix}=A^{-1}\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( x y ) = A − 1 ( 0 0 )
よって、解は
( x y ) = A − 1 ( 0 0 ) = 1 10 ( 3 1 − 4 2 ) ( 0 0 ) = 1 10 ( 3 ⋅ 0 + 1 ⋅ 0 − 4 ⋅ 0 + 2 ⋅ 0 ) = 1 10 ( 0 0 ) = ( 0 0 ) ⋯ (答) \begin{align*}
\begin{pmatrix} x \\ y \end{pmatrix}&=A^{-1}\begin{pmatrix} 0 \\ 0 \end{pmatrix}\\
&=\dfrac{1}{10}\begin{pmatrix*}[r] 3 & 1 \\ -4 & 2 \end{pmatrix*}\begin{pmatrix} 0 \\ 0 \end{pmatrix}\\
&=\dfrac{1}{10}\begin{pmatrix*}[r] 3\cdot0+1\cdot0 \\ -4\cdot0+2\cdot0 \end{pmatrix*}\\
&=\dfrac{1}{10}\begin{pmatrix*}[r] 0 \\ 0 \end{pmatrix*}\\
&=\begin{pmatrix*}[r] 0 \\ 0 \end{pmatrix*} \cdots\small{(答)}
\end{align*} ( x y ) = A − 1 ( 0 0 ) = 10 1 ( 3 − 4 1 2 ) ( 0 0 ) = 10 1 ( 3 ⋅ 0 + 1 ⋅ 0 − 4 ⋅ 0 + 2 ⋅ 0 ) = 10 1 ( 0 0 ) = ( 0 0 ) ⋯ (答)
■ \blacksquare ■
零ゼロ ベクトル ( 0 0 ) \begin{pmatrix} 0 \\ 0 \end{pmatrix} ( 0 0 ) に逆行列 A − 1 A^{-1} A − 1 を掛けているので、明らかに解は零ベクトル ( 0 0 ) \begin{pmatrix} 0 \\ 0 \end{pmatrix} ( 0 0 ) になります。
(グラフで確認) この解について、グラフで考えてみましょう。与式 ( 2 − 1 4 3 ) ( x y ) = ( 0 0 ) \begin{pmatrix*}[r] 2 & -1 \\ 4 & 3 \end{pmatrix*}\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix} ( 2 4 − 1 3 ) ( x y ) = ( 0 0 ) について、行列の積を計算すると
( 2 x − y 4 x + 3 y ) = ( 0 0 ) \begin{pmatrix*}[c]
2x-y \\ 4x+3y
\end{pmatrix*}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix} ( 2 x − y 4 x + 3 y ) = ( 0 0 )
対応する成分を比較して
{ 2 x − y = 0 4 x + 3 y = 0 \begin{cases}
2x-y=0\\
4x+3y=0
\end{cases} { 2 x − y = 0 4 x + 3 y = 0
この2つの直線 2 x − y = 0 2x-y=0 2 x − y = 0 と 4 x + 3 y = 0 4x+3y=0 4 x + 3 y = 0 の交点の座標が、この連立方程式の解になります。 すると、第1式より
y = 2 x y=2x y = 2 x
第2式より
3 y = − 4 x y = − 4 3 x \begin{align*}
3y&=-4x\\
y&=-\dfrac{4}{3}x
\end{align*} 3 y y = − 4 x = − 3 4 x
それぞれのグラフを描くと
この2つの直線の交点は、上図より原点 ( 0 , 0 ) (0, 0) ( 0 , 0 ) なので
( x y ) = ( 0 0 ) ⋯ (答) \begin{align*}
\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix} \cdots\small{(答)}
\end{align*} ( x y ) = ( 0 0 ) ⋯ (答)
2つの直線は原点のみで交わることから、解は ( x y ) = ( 0 0 ) \begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix} ( x y ) = ( 0 0 ) のみで、他に解は存在しません。
(2) A = ( 2 1 4 2 ) A=\begin{pmatrix*}[r] 2 & 1 \\ 4 & 2 \end{pmatrix*} A = ( 2 4 1 2 ) とおくと、A A A の行列式は ∣ A ∣ = 2 ⋅ 2 − 1 ⋅ 4 = 4 − 4 = 0 |A|=2\cdot2-1\cdot4=4-4=0 ∣ A ∣ = 2 ⋅ 2 − 1 ⋅ 4 = 4 − 4 = 0 より、A A A には逆行列が存在しません。 このときは、与式 ( 2 1 4 2 ) ( x y ) = ( 0 0 ) \begin{pmatrix*}[r] 2 & 1 \\ 4 & 2 \end{pmatrix*}\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix} ( 2 4 1 2 ) ( x y ) = ( 0 0 ) について、行列の積を計算すると
( 2 x + y 4 x + 2 y ) = ( 0 0 ) \begin{pmatrix*}[c]
2x+y \\ 4x+2y
\end{pmatrix*}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix} ( 2 x + y 4 x + 2 y ) = ( 0 0 )
対応する成分を比較して
{ 2 x + y = 0 4 x + 2 y = 0 \begin{cases}
2x+y=0\\
4x+2y=0
\end{cases} { 2 x + y = 0 4 x + 2 y = 0
第2式の両辺を 2 2 2 で割ると
{ 2 x + y = 0 2 x + y = 0 \begin{cases}
2x+y=0\\
2x+y=0
\end{cases} { 2 x + y = 0 2 x + y = 0
すると、同じ式 2 x + y = 0 2x+y=0 2 x + y = 0 の連立方程式となるので、この 2 x + y = 0 2x+y=0 2 x + y = 0 を満たすすべての ( x , y ) (x, y) ( x , y ) が、この連立方程式の解になります。 2 x + y = 0 2x+y=0 2 x + y = 0 より
y = − 2 x \begin{align*}
y=-2x
\end{align*} y = − 2 x
C C C を任意の実数として x = C x=C x = C とすると
y = − 2 C \begin{align*}
y=-2C
\end{align*} y = − 2 C
よって、解は C C C を任意の実数として
( x y ) = ( C − 2 C ) = C ( 1 − 2 ) ⋯ (答) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}&=\begin{pmatrix*}[c] C \\ -2C \end{pmatrix*}\\
&=C\begin{pmatrix*}[r] 1 \\ -2 \end{pmatrix*} \cdots(答)
\end{align*} ( x y ) = ( C − 2 C ) = C ( 1 − 2 ) ⋯ (答)
■ \blacksquare ■
(グラフで確認) この解について、グラフで考えてみましょう。 2 x + y = 0 2x+y=0 2 x + y = 0 と 4 x + 2 y = 0 4x+2y=0 4 x + 2 y = 0 は同じ式で、どちらも直線 y = − 2 x y=-2x y = − 2 x を表すので、この直線上のすべての点が解となります。
この直線 y = − 2 x y=-2x y = − 2 x 上のすべての点の座標は、x = C x=C x = C (C C C は任意の実数)とすると、y = − 2 C y=-2C y = − 2 C となるので
( x y ) = ( C − 2 C ) = C ( 1 − 2 ) ⋯ (答) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}&=\begin{pmatrix*}[c] C \\ -2C \end{pmatrix*}\\
&=C\begin{pmatrix*}[r] 1 \\ -2 \end{pmatrix*} \cdots(答)
\end{align*} ( x y ) = ( C − 2 C ) = C ( 1 − 2 ) ⋯ (答)
C C C は任意の実数をとるので、それに対応して ( x y ) \begin{pmatrix} x \\ y \end{pmatrix} ( x y ) は無数の解をもちます。
一般に、A X = O AX=O A X = O の形をした連立一次方程式について、次のことがいえます。
<AX=O の連立一次方程式の解> ( a b c d ) ( x y ) = ( 0 0 ) \begin{pmatrix*}[r] a & b \\ c & d \end{pmatrix*}
\begin{pmatrix} x \\ y \end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix} ( a c b d ) ( x y ) = ( 0 0 )
の解 ( x y ) \begin{pmatrix} x \\ y \end{pmatrix} ( x y ) について、次のことが成り立つ。
(case1) ( a b c d ) \begin{pmatrix*}[r] a & b \\ c & d \end{pmatrix*} ( a c b d ) が逆行列をもつときは
( x y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( x y ) = ( 0 0 )
のみ解をもつ(他に解は存在しない)。
<グラフの視点> a x + b y = 0 ax+by=0 a x + b y = 0 と c x + d y = 0 cx+dy=0 c x + d y = 0 の2つの直線は、原点 ( 0 , 0 ) (0, 0) ( 0 , 0 ) のみで交わるので、解は ( x y ) = ( 0 0 ) \begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix} ( x y ) = ( 0 0 ) のみ。
(case2) ( a b c d ) \begin{pmatrix*}[r] a & b \\ c & d \end{pmatrix*} ( a c b d ) が逆行列をもたないときは、a x + b y = 0 ax+by=0 a x + b y = 0 をみたす無数の解をもつ。
<グラフの視点> a x + b y = 0 ax+by=0 a x + b y = 0 と c x + d y = 0 cx+dy=0 c x + d y = 0 は原点を通る同じ直線となるので、その直線上のすべての点が解となる。
***
以上のように、A X = O AX=O A X = O の解 X X X は、行列 A A A に逆行列が存在するかどうかよって2通りの場合に分かれます。このことを踏まえて、再び例題を与えます。
例題2 次の式を満たす ( x y ) \begin{pmatrix} x \\ y \end{pmatrix} ( x y ) を求めなさい。 (1) ( 3 2 − 1 4 ) ( x y ) = ( 0 0 ) \begin{pmatrix*}[r] 3 & 2 \\ -1 & 4 \end{pmatrix*}\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix}\raisebox{3ex}{ } ( 3 − 1 2 4 ) ( x y ) = ( 0 0 ) (2) ( 3 2 6 4 ) ( x y ) = ( 0 0 ) \begin{pmatrix*}[c] 3 & 2 \\ 6 & 4 \end{pmatrix*}\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix}\raisebox{3ex}{ } ( 3 6 2 4 ) ( x y ) = ( 0 0 )
解答 (1) A = ( 3 2 − 1 4 ) A=\begin{pmatrix*}[r] 3 & 2 \\ -1 & 4 \end{pmatrix*} A = ( 3 − 1 2 4 ) とおくと、与式は
A ( x y ) = ( 0 0 ) \begin{align*}
A\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} A ( x y ) = ( 0 0 )
A A A の行列式は ∣ A ∣ = 3 ⋅ 4 − 2 ⋅ ( − 1 ) = 12 + 2 = 14 ≠ 0 |A|=3\cdot4-2\cdot(-1)=12+2=14\ne0 ∣ A ∣ = 3 ⋅ 4 − 2 ⋅ ( − 1 ) = 12 + 2 = 14 = 0 より、A A A には逆行列 A − 1 A^{-1} A − 1 が存在するので
A − 1 A ( x y ) = A − 1 ( 0 0 ) E ( x y ) = A − 1 ( 0 0 ) \begin{align*}
A^{-1}A\begin{pmatrix} x \\ y \end{pmatrix}&=A^{-1}\begin{pmatrix} 0 \\ 0 \end{pmatrix}\\[10pt]
E\begin{pmatrix} x \\ y \end{pmatrix}&=A^{-1}\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} A − 1 A ( x y ) E ( x y ) = A − 1 ( 0 0 ) = A − 1 ( 0 0 )
単位行列 E E E は掛けても変わらないので
( x y ) = A − 1 ( 0 0 ) = ( 0 0 ) ⋯ (答) \begin{align*}
\begin{pmatrix} x \\ y \end{pmatrix}&=A^{-1}\begin{pmatrix} 0 \\ 0 \end{pmatrix}\\
&=\begin{pmatrix*}[r] 0 \\ 0 \end{pmatrix*} \cdots\small{(答)}\\
\end{align*} ( x y ) = A − 1 ( 0 0 ) = ( 0 0 ) ⋯ (答)
すでに述べた通り、零ゼロ ベクトル ( 0 0 ) \begin{pmatrix} 0 \\ 0 \end{pmatrix} ( 0 0 ) に逆行列を掛けているので、計算しなくても明らかに解は零ベクトル ( 0 0 ) \begin{pmatrix} 0 \\ 0 \end{pmatrix} ( 0 0 ) です。
(2) A = ( 3 2 6 4 ) A=\begin{pmatrix*}[r] 3 & 2 \\ 6 & 4 \end{pmatrix*} A = ( 3 6 2 4 ) とおくと、A A A の行列式は ∣ A ∣ = 3 ⋅ 4 − 2 ⋅ 6 = 12 − 12 = 0 |A|=3\cdot4-2\cdot6=12-12=0 ∣ A ∣ = 3 ⋅ 4 − 2 ⋅ 6 = 12 − 12 = 0 より、A A A には逆行列は存在しません。 よって、与式 ( 3 2 6 4 ) ( x y ) = ( 0 0 ) \begin{pmatrix*}[r] 3 & 2 \\ 6 & 4 \end{pmatrix*}\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix} ( 3 6 2 4 ) ( x y ) = ( 0 0 ) について、行列の積を計算すると
( 3 x + 2 y 6 x + 4 y ) = ( 0 0 ) \begin{pmatrix*}[c]
3x+2y \\ 6x+4y
\end{pmatrix*}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix} ( 3 x + 2 y 6 x + 4 y ) = ( 0 0 )
対応する成分を比較して
{ 3 x + 2 y = 0 6 x + 4 y = 0 \begin{cases}
3x+2y=0\\
6x+4y=0
\end{cases} { 3 x + 2 y = 0 6 x + 4 y = 0
第2式を 2 2 2 で割ると
{ 3 x + 2 y = 0 3 x + 2 y = 0 \begin{cases}
3x+2y=0\\
3x+2y=0
\end{cases} { 3 x + 2 y = 0 3 x + 2 y = 0
すると、同じ式の連立方程式となるので、この 3 x + 2 y = 0 3x+2y=0 3 x + 2 y = 0 を満たすすべての ( x , y ) (x, y) ( x , y ) が、この連立方程式の解となります。 3 x + 2 y = 0 3x+2y=0 3 x + 2 y = 0 より
2 y = − 3 x y = − 3 2 x \begin{align*}
2y&=-3x\\
y&=-\dfrac{3}{2}x
\end{align*} 2 y y = − 3 x = − 2 3 x
C C C を任意の実数として x = 2 C x=2C x = 2 C とすると
y = − 3 2 ⋅ 2 C = − 3 C \begin{align*}
\\[-14pt]
y=-\dfrac{3}{\cancel{2}}\cdot\cancel{2}C=-3C
\end{align*} y = − 2 3 ⋅ 2 C = − 3 C
よって、解は C C C を任意の実数として
( x y ) = ( 2 C − 3 C ) = C ( 2 − 3 ) ⋯ (答) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}&=\begin{pmatrix*}[c] 2C \\ -3C \end{pmatrix*}\\
&=C\begin{pmatrix*}[r] 2 \\ -3 \end{pmatrix*} \cdots(答)
\end{align*} ( x y ) = ( 2 C − 3 C ) = C ( 2 − 3 ) ⋯ (答)
■ \blacksquare ■
ここから、本記事のテーマである固有値と固有ベクトルに入ります。その前に、1次変換の復習です(中学でも分かる⑮「一次変換」)。
1次変換の定義(復習) a , b , c , d a, b, c, d a , b , c , d を定数として
x ′ = a x + b y y ′ = c x + d y \begin{align*}
x'&=ax+by\\
y'&=cx+dy
\end{align*} x ′ y ′ = a x + b y = c x + d y
によって点 P( x , y ) (x, y) ( x , y ) と点 Q( x ′ , y ′ ) (x', y') ( x ′ , y ′ ) の関係を定めると、座標平面上の1つの変換が得られます。 一般に、この形の式によって表される変換
f : ( x , y ) ⟶ ( x ′ , y ′ ) f:(x, y) \longrightarrow (x', y') f : ( x , y ) ⟶ ( x ′ , y ′ )
を1次変換 といいます。1次変換 f f f は、行列を用いると
( x ′ y ′ ) = ( a b c d ) ( x y ) ⋯ ( ∗ 1 ) \begin{align*}
\begin{pmatrix}
x' \\ y'
\end{pmatrix}
=\begin{pmatrix}
a & b \\ c & d
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix} \cdots (*1)
\end{align*} ( x ′ y ′ ) = ( a c b d ) ( x y ) ⋯ ( ∗ 1 )
と表すことができます。つまり、変換 f f f は行列
( a b c d ) \begin{align*}
\begin{pmatrix}
a & b \\ c & d
\end{pmatrix}
\end{align*} ( a c b d )
によって定まるので、この変換 f f f を、行列 ( a b c d ) \begin{pmatrix} a & b \\ c & d \end{pmatrix} ( a c b d ) で表される1次変換 といいます。また、この行列のことを変換行列 とよびます。
f : P ( x , y ) → ( a b c d ) 1次変換 Q ( x ′ , y ′ ) \begin{align*}
\\[-12pt]
f:\mathrm{P}(x, y) \xrightarrow[\scriptsize \begin{pmatrix} a & b \\ c & d \end{pmatrix}]{1次変換 } \mathrm{Q}(x', y')
\end{align*} f : P ( x , y ) 1次変換 ( a c b d ) Q ( x ′ , y ′ )
***
一次変換をベクトルを移す変換とみなす この一次変換
( x ′ y ′ ) = ( a b c d ) ( x y ) ⋯ ( ∗ 1 ) \begin{align*}
\begin{pmatrix}
x' \\ y'
\end{pmatrix}
=\begin{pmatrix}
a & b \\ c & d
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix} \cdots(*1)
\end{align*} ( x ′ y ′ ) = ( a c b d ) ( x y ) ⋯ ( ∗ 1 )
について、上の定義では ( x y ) \begin{pmatrix}x \\ y\end{pmatrix} ( x y ) は平面上の1点の座標を表していましたが、これをベクトル O P → = ( x y ) \overrightarrow{\mathrm{OP}}=\begin{pmatrix}x \\ y\end{pmatrix} OP = ( x y ) として考えることができます。なお、O P → = ( x y ) \overrightarrow{\mathrm{OP}}=\begin{pmatrix}x \\ y\end{pmatrix} OP = ( x y ) は、高校で習うベクトル表記 O P → = ( x , y ) \overrightarrow{\mathrm{OP}}=(x, y) OP = ( x , y ) を縦書きにし、列ベクトルで表したものです(中学でも分かる⑬「行列」)。ベクトルの本質である向き・大きさ・成分は変わらず、行列の計算と相性をよくするために表記を変えているだけです。 ベクトルの始点を原点に置いた場合、そのベクトルの終点の座標がベクトルの成分となります。これにより、ある点 P \mathrm{P} P の座標 ( x y ) \begin{pmatrix}x \\ y\end{pmatrix} ( x y ) は、原点 O \mathrm{O} O から P \mathrm{P} P へのベクトル O P → = ( x y ) \overrightarrow{\mathrm{OP}}= \begin{pmatrix}x \\ y\end{pmatrix} OP = ( x y ) と同一視することができます(中学でも分かる⑫「ベクトル」)。
つまり、上記の ( ∗ 1 ) (*1) ( ∗ 1 ) 式
( x ′ y ′ ) = ( a b c d ) ( x y ) ⋯ ( ∗ 3 ) \begin{align*}
\begin{pmatrix}
x' \\ y'
\end{pmatrix}
=\begin{pmatrix}
a & b \\ c & d
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix} \cdots(*3)
\end{align*} ( x ′ y ′ ) = ( a c b d ) ( x y ) ⋯ ( ∗ 3 )
は、「ベクトル O P → = ( x y ) \overrightarrow{\mathrm{OP}}=\begin{pmatrix}x \\ y\end{pmatrix} OP = ( x y ) が、変換行列 ( a b c d ) \begin{pmatrix} a & b \\ c & d\end{pmatrix} ( a c b d ) によって、ベクトル O Q → = ( x ′ y ′ ) \overrightarrow{\mathrm{OQ}}=\begin{pmatrix}x' \\ y' \end{pmatrix} OQ = ( x ′ y ′ ) に変換される」とみることができます。一次変換は、「ベクトルをベクトルに写す」という新しいものの見方です。
f : O P → = ( x y ) → ( a b c d ) 1次変換 O Q → = ( x ′ y ′ ) \begin{align*}
\\[-10pt]
f:\overrightarrow{\mathrm{OP}}=\begin{pmatrix}x \\ y\end{pmatrix} \xrightarrow[\scriptsize \begin{pmatrix} a & b \\ c & d \end{pmatrix}]{1次変換 } \overrightarrow{\mathrm{OQ}}=\begin{pmatrix} x' \\ y' \end{pmatrix}
\end{align*} f : OP = ( x y ) 1次変換 ( a c b d ) OQ = ( x ′ y ′ )
ベクトルは、向きと大きさが同じであれば、その位置に関係なく同じベクトルとみなします。したがって、どのような位置にあるベクトルも、始点が原点 O \mathrm{O} O になるように平行移動して考えることができます。このとき、原点 O \mathrm{O} O を始点とするベクトルを位置ベクトル といいます。 このようにみれば、一次変換は「座標を他の座標に変換する操作 」であると同時に、「位置ベクトルを他の位置ベクトルに変換する操作 」、より一般には「ベクトルを他のベクトルに変換する操作 」と考えることができます。
これが、「固有ベクトル」を理解するための最初のステップです。 以上を踏まえ、固有値と固有ベクトルの解説です。
固有値と固有ベクトル まずは、行列の固有値、固有ベクトルがどのようなものか、具体例を通してみていきます。
固有値、固有ベクトルの具体例 例えば、1次変換の変換行列として、次のような2次正方行列 A A A があるとします。
A = ( 1 2 − 1 4 ) \begin{align*}
A=\begin{pmatrix*}[r]
1 & 2 \\ -1 & 4
\end{pmatrix*}
\end{align*} A = ( 1 − 1 2 4 )
この変換行列 A A A を、ある位置ベクトル O P → \overrightarrow{\mathrm{OP}} OP に掛けて、別の位置ベクトル O Q → \overrightarrow{\mathrm{OQ}} OQ に変換することを考えます。すると
(例1) O P → = ( 2 3 ) \overrightarrow{\mathrm{OP}}=\begin{pmatrix} 2 \\ 3 \end{pmatrix} OP = ( 2 3 ) のとき
O Q → = A O P → = ( 1 2 − 1 4 ) ( 2 3 ) = ( 1 ⋅ 2 + 2 ⋅ 3 ( − 1 ) ⋅ 2 + 4 ⋅ 3 ) = ( 2 + 6 − 2 + 12 ) = ( 8 10 ) \begin{align*}
\overrightarrow{\mathrm{OQ}}&=A\,\overrightarrow{\mathrm{OP}}\\
&=\begin{pmatrix*}[r]
1 & 2 \\ -1 & 4
\end{pmatrix*}
\begin{pmatrix}
2 \\ 3
\end{pmatrix}\\
&=\begin{pmatrix}
1\cdot2+2\cdot3 \\ (-1)\cdot2+4\cdot3
\end{pmatrix}\\
&=\begin{pmatrix*}[c]
2+6 \\ -2+12
\end{pmatrix*}\\
&=\begin{pmatrix}
8 \\ 10
\end{pmatrix}
\end{align*} OQ = A OP = ( 1 − 1 2 4 ) ( 2 3 ) = ( 1 ⋅ 2 + 2 ⋅ 3 ( − 1 ) ⋅ 2 + 4 ⋅ 3 ) = ( 2 + 6 − 2 + 12 ) = ( 8 10 )
より
( 1 2 − 1 4 ) ( 2 3 ) = ( 8 10 ) \begin{align*}
\begin{pmatrix*}[r]
1 & 2 \\ -1 & 4
\end{pmatrix*}
\begin{pmatrix}
2 \\ 3
\end{pmatrix}
=\begin{pmatrix}
8 \\ 10
\end{pmatrix}
\end{align*} ( 1 − 1 2 4 ) ( 2 3 ) = ( 8 10 )
上図のように、ベクトル O P → = ( 2 3 ) \overrightarrow{\mathrm{OP}}=\begin{pmatrix} 2 \\ 3 \end{pmatrix} OP = ( 2 3 ) は、変換行列 ( 1 2 − 1 4 ) \begin{pmatrix*}[r] 1 & 2 \\ -1 & 4 \end{pmatrix*} ( 1 − 1 2 4 ) によって、別のベクトル O Q → = ( 8 10 ) \overrightarrow{\mathrm{OQ}}=\begin{pmatrix} 8 \\ 10 \end{pmatrix} OQ = ( 8 10 ) に写されます。
(例2)
O P → = ( 2 − 1 ) \overrightarrow{\mathrm{OP}}=\begin{pmatrix*}[r] 2 \\ -1 \end{pmatrix*} OP = ( 2 − 1 ) のとき
O Q → = A O P → = ( 1 2 − 1 4 ) ( 2 − 1 ) = ( 1 ⋅ 2 + 2 ⋅ ( − 1 ) ( − 1 ) ⋅ 2 + 4 ⋅ ( − 1 ) ) = ( 2 − 2 − 2 − 4 ) = ( 0 − 6 ) \begin{align*}
\overrightarrow{\mathrm{OQ}}&=A\,\overrightarrow{\mathrm{OP}}\\
&=\begin{pmatrix*}[r]
1 & 2 \\ -1 & 4
\end{pmatrix*}
\begin{pmatrix*}[r]
2 \\ -1
\end{pmatrix*}\\
&=\begin{pmatrix}
1\cdot2+2\cdot(-1) \\ (-1)\cdot2+4\cdot(-1)
\end{pmatrix}\\
&=\begin{pmatrix}
2-2 \\ -2-4
\end{pmatrix}\\
&=\begin{pmatrix*}[r]
0 \\ -6
\end{pmatrix*}
\end{align*} OQ = A OP = ( 1 − 1 2 4 ) ( 2 − 1 ) = ( 1 ⋅ 2 + 2 ⋅ ( − 1 ) ( − 1 ) ⋅ 2 + 4 ⋅ ( − 1 ) ) = ( 2 − 2 − 2 − 4 ) = ( 0 − 6 )
より
( 1 2 − 1 4 ) ( 2 − 1 ) = ( 0 − 6 ) \begin{align*}
\begin{pmatrix*}[r]
1 & 2 \\ -1 & 4
\end{pmatrix*}
\begin{pmatrix}
2 \\ -1
\end{pmatrix}
=\begin{pmatrix}
0 \\ -6
\end{pmatrix}
\end{align*} ( 1 − 1 2 4 ) ( 2 − 1 ) = ( 0 − 6 )
上図のように、ベクトル O P → = ( 2 − 1 ) \overrightarrow{\mathrm{OP}}=\begin{pmatrix*}[r] 2 \\ -1 \end{pmatrix*} OP = ( 2 − 1 ) は、変換行列 ( 1 2 − 1 4 ) \begin{pmatrix*}[r] 1 & 2 \\ -1 & 4 \end{pmatrix*} ( 1 − 1 2 4 ) によって、別のベクトル O Q → = ( 0 − 6 ) \overrightarrow{\mathrm{OQ}}=\begin{pmatrix*}[r] 0 \\ -6 \end{pmatrix*} OQ = ( 0 − 6 ) に写されます。
上記の例では、どちらも変換後のベクトルの「向き」と「大きさ(長さ)」は変わっています。しかし、中には
(例3)
O P → = ( 2 1 ) \overrightarrow{\mathrm{OP}}=\begin{pmatrix*}[r] 2 \\ 1 \end{pmatrix*} OP = ( 2 1 ) のとき
A p ⃗ = ( 1 2 − 1 4 ) ( 2 1 ) = ( 1 ⋅ 2 + 2 ⋅ 1 ( − 1 ) ⋅ 2 + 4 ⋅ 1 ) = ( 2 + 2 − 2 + 4 ) = ( 4 2 ) = 2 ( 2 1 ) \begin{align*}
A\vec{p}&=\begin{pmatrix*}[r]
1 & 2 \\ -1 & 4
\end{pmatrix*}
\begin{pmatrix}
2 \\ 1
\end{pmatrix}\\
&=\begin{pmatrix}
1\cdot2+2\cdot1 \\ (-1)\cdot2+4\cdot1
\end{pmatrix}\\
&=\begin{pmatrix}
2+2 \\ -2+4
\end{pmatrix}\\
&=\begin{pmatrix}
4 \\ 2
\end{pmatrix}\\
&=2\begin{pmatrix}
2 \\ 1
\end{pmatrix}
\end{align*} A p = ( 1 − 1 2 4 ) ( 2 1 ) = ( 1 ⋅ 2 + 2 ⋅ 1 ( − 1 ) ⋅ 2 + 4 ⋅ 1 ) = ( 2 + 2 − 2 + 4 ) = ( 4 2 ) = 2 ( 2 1 )
より
( 1 2 − 1 4 ) ( 2 1 ) = 2 ( 2 1 ) \begin{align*}
\begin{pmatrix*}[r]
1 & 2 \\ -1 & 4
\end{pmatrix*}
\begin{pmatrix}
2 \\ 1
\end{pmatrix}
=2\begin{pmatrix}
2 \\ 1
\end{pmatrix}
\end{align*} ( 1 − 1 2 4 ) ( 2 1 ) = 2 ( 2 1 )
この場合、O Q → \overrightarrow{\mathrm{OQ}} OQ は O P → \overrightarrow{\mathrm{OP}} OP と向きが同じで、大きさだけが 2 2 2 倍になります。
(例4)
A p ⃗ = ( 1 2 − 1 4 ) ( 1 1 ) = ( 1 ⋅ 1 + 2 ⋅ 1 ( − 1 ) ⋅ 1 + 4 ⋅ 1 ) = ( 1 + 2 − 1 + 4 ) = ( 3 3 ) = 3 ( 1 1 ) \begin{align*}
A\vec{p}&=\begin{pmatrix*}[r]
1 & 2 \\ -1 & 4
\end{pmatrix*}
\begin{pmatrix}
1 \\ 1
\end{pmatrix}\\
&=\begin{pmatrix}
1\cdot1+2\cdot1 \\ (-1)\cdot1+4\cdot1
\end{pmatrix}\\
&=\begin{pmatrix}
1+2 \\ -1+4
\end{pmatrix}\\
&=\begin{pmatrix}
3 \\ 3
\end{pmatrix}\\
&=3\begin{pmatrix}
1 \\ 1
\end{pmatrix}
\end{align*} A p = ( 1 − 1 2 4 ) ( 1 1 ) = ( 1 ⋅ 1 + 2 ⋅ 1 ( − 1 ) ⋅ 1 + 4 ⋅ 1 ) = ( 1 + 2 − 1 + 4 ) = ( 3 3 ) = 3 ( 1 1 )
より
( 1 2 − 1 4 ) ( 1 1 ) = 3 ( 1 1 ) \begin{align*}
\begin{pmatrix*}[r]
1 & 2 \\ -1 & 4
\end{pmatrix*}
\begin{pmatrix}
1 \\ 1
\end{pmatrix}
=3\begin{pmatrix}
1 \\ 1
\end{pmatrix}
\end{align*} ( 1 − 1 2 4 ) ( 1 1 ) = 3 ( 1 1 )
この場合、O Q → \overrightarrow{\mathrm{OQ}} OQ は O P → \overrightarrow{\mathrm{OP}} OP と向きが同じで、大きさだけが 3 3 3 倍になります。
これらのように、「大きさ」が実数倍されるだけで、変換しても「向き」が変わらないベクトルがあります。 このように、ある変換行列に対して、ベクトル を掛けても「向き」が変わらない特別なベクトルを固有ベクトル といいます。また、そのときのベクトルの大きさの変化を表す値(何倍するのかの値)を固有値 といいます。 (例3) の
( 1 2 − 1 4 ) ( 2 1 ) = 2 ( 2 1 ) \begin{align*}
\begin{pmatrix*}[r]
1 & 2 \\ -1 & 4
\end{pmatrix*}
\begin{pmatrix}
2 \\ 1
\end{pmatrix}
=2\begin{pmatrix}
2 \\ 1
\end{pmatrix}
\end{align*} ( 1 − 1 2 4 ) ( 2 1 ) = 2 ( 2 1 )
については、固有ベクトルは ( 2 1 ) \begin{pmatrix} 2 \\ 1 \end{pmatrix} ( 2 1 ) 、固有値は 2 2 2 となります。 (例4) の
( 1 2 − 1 4 ) ( 1 1 ) = 3 ( 1 1 ) \begin{align*}
\begin{pmatrix*}[r]
1 & 2 \\ -1 & 4
\end{pmatrix*}
\begin{pmatrix}
1 \\ 1
\end{pmatrix}
=3\begin{pmatrix}
1 \\ 1
\end{pmatrix}
\end{align*} ( 1 − 1 2 4 ) ( 1 1 ) = 3 ( 1 1 )
については、固有ベクトルは ( 1 1 ) \begin{pmatrix} 1 \\ 1 \end{pmatrix} ( 1 1 ) 、固有値は 3 3 3 となります。
固有ベクトルは、「変換行列によって向きが変わらない」という、極めて特別な性質をもちます。これは偶然ではなく、行列そのものの性質を浮き彫りにする、ある種の「骨格」のような存在です。
(注意1)固有ベクトルは零ベクトルを含まない なお、 固有ベクトルは零ベクトル 0 ⃗ \vec{0} 0 を含みません。例えば C C C を任意の実数として
(例5)
( 1 2 − 1 4 ) ( 0 0 ) = C ( 0 0 ) \begin{align*}
\begin{pmatrix*}[r]
1 & 2 \\ -1 & 4
\end{pmatrix*}
\begin{pmatrix}
0 \\ 0
\end{pmatrix}
=C\begin{pmatrix}
0 \\ 0
\end{pmatrix}
\end{align*} ( 1 − 1 2 4 ) ( 0 0 ) = C ( 0 0 )
という式はつねに成り立ちますが、このとき、零ベクトル ( 0 0 ) \begin{pmatrix} 0 \\ 0 \end{pmatrix} ( 0 0 ) は固有ベクトルではありません。
固有ベクトルとは、「変換行列を掛けたときに、方向は変わらず、大きさだけが伸び縮みするベクトル」のことです。そのため、固有ベクトルは「大きさがあり、特定の方向をもつ」ことが不可欠なので、大きさも方向ももたない零ベクトルは固有ベクトルになりません。
あらためて、固有ベクトルと固有値の教科書的な定義を与えましょう。
<固有値と固有値ベクトルの定義> ある正方行列 A A A に対して、零ではない ベクトル v ⃗ \vec{v} v と、あるスカラー λ \lambda λ が存在し、以下の関係が成り立つとき、v ⃗ \vec{v} v を固有ベクトル 、 λ \lambda λ を固有値 という。
A v ⃗ = λ v ⃗ A\vec{v}=\lambda\vec{v} A v = λ v
***
なお、ここでは行列の成分、および固有値は、実数の範囲で考えることとします。量子力学などの特定の分野では、行列の成分、および固有値は、複素数まで拡張されることがあります。
行列 A A A を2次正方行列 ( a b c d ) \begin{pmatrix}a & b \\ c & d\end{pmatrix} ( a c b d ) で考えると、上の定義は次のようになります。
( a b c d ) ( x y ) = λ ( x y ) \begin{align*}
\begin{pmatrix}
a & b \\ c & d
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\lambda
\begin{pmatrix}
x \\ y
\end{pmatrix}
\end{align*} ( a c b d ) ( x y ) = λ ( x y )
上記の関係式が成り立つとき、( x y ) ( ≠ 0 ⃗ ) \begin{pmatrix} x \\ y \end{pmatrix}(\ne{\vec{0}}) ( x y ) ( = 0 ) が固有ベクトル、λ \lambda λ が固有値となります。
固有値と固有ベクトルを求める計算 ある行列に対して、その固有値と固有ベクトルを求めることを考えましょう。まずは次の例題から。
例題3 行列 A = ( 2 1 1 2 ) A=\begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix} A = ( 2 1 1 2 ) の固有値と固有ベクトルを求めなさい。
解答 2次正方行列 A = ( 2 1 1 2 ) A=\begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix} A = ( 2 1 1 2 ) に対して、固有値を λ \lambda λ 、固有ベクトルを ( x y ) ( ≠ 0 ⃗ ) \begin{pmatrix} x \\ y \end{pmatrix}(\ne\vec{0}) ( x y ) ( = 0 ) とすると、以下の関係式が成り立ちます。
( 2 1 1 2 ) ( x y ) = λ ( x y ) \begin{align*}
\begin{pmatrix}
2 & 1 \\ 1 & 2
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\lambda
\begin{pmatrix}
x \\ y
\end{pmatrix}
\end{align*} ( 2 1 1 2 ) ( x y ) = λ ( x y )
ここで、単位行列 E = ( 1 0 0 1 ) E=\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} E = ( 1 0 0 1 ) を掛けても変わらないので
( 2 1 1 2 ) ( x y ) = λ ( 1 0 0 1 ) ( x y ) \begin{align*}
\begin{pmatrix}
2 & 1 \\ 1 & 2
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\lambda
\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
\end{align*} ( 2 1 1 2 ) ( x y ) = λ ( 1 0 0 1 ) ( x y )
より
( 2 1 1 2 ) ( x y ) = ( λ 0 0 λ ) ( x y ) \begin{align*}
\begin{pmatrix}
2 & 1 \\ 1 & 2
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} \lambda & 0 \\ 0 & \lambda \end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
\end{align*} ( 2 1 1 2 ) ( x y ) = ( λ 0 0 λ ) ( x y )
左辺から右辺を引いて
( 2 1 1 2 ) ( x y ) − ( λ 0 0 λ ) ( x y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
2 & 1 \\ 1 & 2
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
-\begin{pmatrix} \lambda & 0 \\ 0 & \lambda \end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( 2 1 1 2 ) ( x y ) − ( λ 0 0 λ ) ( x y ) = ( 0 0 )
この式の右辺は、零ベクトル ( 0 0 ) \begin{pmatrix} 0 \\ 0 \end{pmatrix} ( 0 0 ) になることに注意しましょう。左辺を計算すると2行1列の行列(列ベクトル)になるので、右辺も同じ型の零ベクトルになります。等式の両辺は同じ型でなければならないので、実数の 0 0 0 にはならないことに注意しましょう。 左辺を ( x y ) \begin{pmatrix} x \\ y \end{pmatrix} ( x y ) でくくって
{ ( 2 1 1 2 ) − ( λ 0 0 λ ) } ( x y ) = ( 0 0 ) \begin{align*}
\left\{\begin{pmatrix}
2 & 1 \\ 1 & 2
\end{pmatrix}
-\begin{pmatrix} \lambda & 0 \\ 0 & \lambda \end{pmatrix}\right\}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} { ( 2 1 1 2 ) − ( λ 0 0 λ ) } ( x y ) = ( 0 0 )
行列のひき算を計算して
( 2 − λ 1 − 0 1 − 0 2 − λ ) ( x y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
2-\lambda & 1-0 \\ 1-0 & 2-\lambda
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( 2 − λ 1 − 0 1 − 0 2 − λ ) ( x y ) = ( 0 0 )
より
( 2 − λ 1 1 2 − λ ) ( x y ) = ( 0 0 ) ⋯ ( ∗ 2 ) \begin{align*}
\begin{pmatrix}
2-\lambda & 1 \\ 1 & 2-\lambda
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} \cdots(*2) ( 2 − λ 1 1 2 − λ ) ( x y ) = ( 0 0 ) ⋯ ( ∗ 2 )
後は、この方程式を解きます。 これには、この記事ですでにやった<AX= O の連立1次方程式の解> が役立ちます。
<AX= O の連立1次方程式の解>(再掲)
( a b c d ) ( x y ) = ( 0 0 ) \begin{pmatrix*}[r] a & b \\ c & d \end{pmatrix*}
\begin{pmatrix} x \\ y \end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix} ( a c b d ) ( x y ) = ( 0 0 )
の解 ( x y ) \begin{pmatrix} x \\ y \end{pmatrix} ( x y ) について、次のことが成り立つ。
(case1) ( a b c d ) \begin{pmatrix*}[r] a & b \\ c & d \end{pmatrix*} ( a c b d ) が逆行列をもつときは
( x y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( x y ) = ( 0 0 )
のみ解をもつ(他に解は存在しない)。
(case2) ( a b c d ) \begin{pmatrix*}[r] a & b \\ c & d \end{pmatrix*} ( a c b d ) が逆行列をもたないときは、a x + b y = 0 ax+by=0 a x + b y = 0 をみたす無数の解をもつ。
***
すると、行列
( 2 − λ 1 1 2 − λ ) \begin{align*}
\begin{pmatrix}
2-\lambda & 1 \\ 1 & 2-\lambda
\end{pmatrix}
\end{align*} ( 2 − λ 1 1 2 − λ )
が逆行列をもつときは、(case1) より、( ∗ 2 ) (*2) ( ∗ 2 ) 式
( 2 − λ 1 1 2 − λ ) ( x y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
2-\lambda & 1 \\ 1 & 2-\lambda
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( 2 − λ 1 1 2 − λ ) ( x y ) = ( 0 0 )
の解は零ベクトル
( x y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( x y ) = ( 0 0 )
となりますが、零ベクトルは固有ベクトルにはならないので不適となります。 よって、(case2) のとき、つまり、行列
( 2 − λ 1 1 2 − λ ) \begin{align*}
\begin{pmatrix}
2-\lambda & 1 \\ 1 & 2-\lambda
\end{pmatrix}
\end{align*} ( 2 − λ 1 1 2 − λ )
が逆行列をもたないときを考えます。 この行列が逆行列をもたないための条件は、最初にやった<逆行列の定義> より、行列式
∣ 2 − λ 1 1 2 − λ ∣ \begin{align*}
\begin{vmatrix}
2-\lambda & 1 \\ 1 & 2-\lambda
\end{vmatrix}
\end{align*} 2 − λ 1 1 2 − λ
が 0 0 0 になることなので
∣ 2 − λ 1 1 2 − λ ∣ = ( 2 − λ ) 2 − 1 = 2 2 − 2 ⋅ 2 λ + λ 2 − 1 = λ 2 − 4 λ + 4 − 1 = λ 2 − 4 λ + 3 = ( λ − 1 ) ( λ − 3 ) = 0 \begin{align*}
\begin{vmatrix}
2-\lambda & 1 \\ 1 & 2-\lambda
\end{vmatrix}&=(2 - \lambda)^2-1\\
&=2^2-2\cdot2\lambda+\lambda^2-1\\
&=\lambda^2-4\lambda+4-1\\
&=\lambda^2-4\lambda+3\\
&=(\lambda-1)(\lambda-3)\\
&=0
\end{align*} 2 − λ 1 1 2 − λ = ( 2 − λ ) 2 − 1 = 2 2 − 2 ⋅ 2 λ + λ 2 − 1 = λ 2 − 4 λ + 4 − 1 = λ 2 − 4 λ + 3 = ( λ − 1 ) ( λ − 3 ) = 0
よって
λ − 1 = 0 または λ − 3 = 0 \begin{align*}
\lambda-1=0 または \lambda-3=0
\end{align*} λ − 1 = 0 または λ − 3 = 0
より
λ = 1 または λ = 3 \begin{align*}
\lambda=1 または \lambda=3
\end{align*} λ = 1 または λ = 3
(ⅰ) λ = 1 \lambda=1 λ = 1 のとき
( ∗ 2 ) (*2) ( ∗ 2 ) に代入して
( 2 − 1 1 1 2 − 1 ) ( x y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
2-1 & 1 \\ 1 & 2-1
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( 2 − 1 1 1 2 − 1 ) ( x y ) = ( 0 0 )
より
( 1 1 1 1 ) ( x y ) = ( 0 0 ) ( x + y x + y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
1 & 1 \\ 1 & 1
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
&=\begin{pmatrix} 0 \\ 0 \end{pmatrix}\\
\begin{pmatrix}
x+y \\ x+y
\end{pmatrix}
&=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( 1 1 1 1 ) ( x y ) ( x + y x + y ) = ( 0 0 ) = ( 0 0 )
対応する成分は等しいので
{ x + y = 0 x + y = 0 \begin{cases}
x+y=0\\
x+y=0
\end{cases} { x + y = 0 x + y = 0
すると、同じ式の連立方程式となるので、この連立方程式の解は x + y = 0 x+y=0 x + y = 0 を満たすすべての実数が解となります。すると
y = − x \begin{align*}
y=-x
\end{align*} y = − x
より、C C C を任意の実数として x = C x=C x = C とすると
y = − C y=-C y = − C
したがって、この連立方程式の解は、C C C を任意の実数として
( x y ) = ( 3 C − 2 C ) = C ( 3 − 2 ) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}&=\begin{pmatrix*}[c] 3C \\ -2C \end{pmatrix*}\\
&=C\begin{pmatrix*}[r] 3 \\ -2 \end{pmatrix*}
\end{align*} ( x y ) = ( 3 C − 2 C ) = C ( 3 − 2 )
ここで、特にことわりの無い場合は、固有値はもっとも簡単な値を選択します。もっとも簡単な値は C = 1 C=1 C = 1 のときなので、λ = 1 \lambda=1 λ = 1 のときの固有値は、C = 1 C=1 C = 1 のときを選んで
( x y ) = 1 ⋅ ( 3 − 2 ) = ( 3 − 2 ) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}
&=1\cdot\begin{pmatrix*}[r] 3 \\ -2 \end{pmatrix*}\\
&=\begin{pmatrix*}[r] 3 \\ -2 \end{pmatrix*}
\end{align*} ( x y ) = 1 ⋅ ( 3 − 2 ) = ( 3 − 2 )
(ⅱ) λ = 3 \lambda=3 λ = 3 のとき
( ∗ 2 ) (*2) ( ∗ 2 ) に代入して
( 2 − 3 1 1 2 − 3 ) ( x y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
2-3 & 1 \\ 1 & 2-3
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( 2 − 3 1 1 2 − 3 ) ( x y ) = ( 0 0 )
より
( − 1 1 1 − 1 ) ( x y ) = ( 0 0 ) ( − x + y x − y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
-1 & 1 \\ 1 & -1
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
&=\begin{pmatrix} 0 \\ 0 \end{pmatrix}\\
\begin{pmatrix}
-x+y \\ x-y
\end{pmatrix}
&=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( − 1 1 1 − 1 ) ( x y ) ( − x + y x − y ) = ( 0 0 ) = ( 0 0 )
対応する成分は等しいので
{ − x + y = 0 x − y = 0 \begin{cases}
-x+y=0\\
x-y=0
\end{cases} { − x + y = 0 x − y = 0
第1式を − 1 -1 − 1 倍すると
{ x − y = 0 x − y = 0 \begin{cases}
x-y=0\\
x-y=0
\end{cases} { x − y = 0 x − y = 0
すると、同じ式の連立方程式となるので、この連立方程式の解は x − y = 0 x-y=0 x − y = 0 を満たすすべての実数となります。すると
y = x \begin{align*}
y=x
\end{align*} y = x
より、C C C を任意の実数として x = C x=C x = C とすると
y = C y=C y = C
したがって、この連立方程式の解は、C C C を任意の実数として
( x y ) = ( C C ) = C ( 1 1 ) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}&=\begin{pmatrix*}[c] C \\ C \end{pmatrix*}\\
&=C\begin{pmatrix*}[r] 1 \\ 1 \end{pmatrix*}
\end{align*} ( x y ) = ( C C ) = C ( 1 1 )
よって、λ = 3 \lambda=3 λ = 3 のときの固有ベクトルは、最も簡単な C = 1 C=1 C = 1 のときを選んで
( x y ) = 1 ⋅ ( 1 1 ) = ( 1 1 ) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}
&=1\cdot\begin{pmatrix*}[r] 1 \\ 1 \end{pmatrix*}\\
&=\begin{pmatrix*}[r] 1 \\ 1 \end{pmatrix*}
\end{align*} ( x y ) = 1 ⋅ ( 1 1 ) = ( 1 1 )
以上より、
固有値が λ = 1 \lambda=1 λ = 1 のときの固有ベクトルは
( x y ) = ( 1 − 1 ) ⋯ (答) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}&=\begin{pmatrix*}[r] 1 \\ -1 \end{pmatrix*} \cdots(答)
\end{align*} ( x y ) = ( 1 − 1 ) ⋯ (答)
固有値が λ = 3 \lambda=3 λ = 3 のときの固有ベクトルは
( x y ) = ( 1 1 ) ⋯ (答) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}&=\begin{pmatrix*}[r] 1 \\ 1 \end{pmatrix*} \cdots(答)
\end{align*} ( x y ) = ( 1 1 ) ⋯ (答)
■ \blacksquare ■
後のために、上記でやった固有値、固有ベクトルを求め方を、文字を用いて一般化しましょう。
固有値、固有ベクトルの求め方の一般化 2次正方行列 A = ( a b c d ) A=\begin{pmatrix}a & b \\ c & d\end{pmatrix} A = ( a c b d ) に対して、固有ベクトルを ( x y ) ( ≠ 0 ⃗ ) \begin{pmatrix} x \\ y \end{pmatrix}(\ne\vec{0}) ( x y ) ( = 0 ) 、固有値を λ \lambda λ とすると、以下の関係が成り立ちます。
( a b c d ) ( x y ) = λ ( x y ) \begin{align*}
\begin{pmatrix}
a & b \\ c & d
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\lambda
\begin{pmatrix}
x \\ y
\end{pmatrix}
\end{align*} ( a c b d ) ( x y ) = λ ( x y )
これに、単位行列 ( 1 0 0 1 ) \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} ( 1 0 0 1 ) を掛けて
( a b c d ) ( x y ) = λ ( 1 0 0 1 ) ( x y ) \begin{align*}
\begin{pmatrix}
a & b \\ c & d
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\lambda\begin{pmatrix}
1 & 0 \\ 0 & 1
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
\end{align*} ( a c b d ) ( x y ) = λ ( 1 0 0 1 ) ( x y )
より
( a b c d ) ( x y ) = ( λ 0 0 λ ) ( x y ) \begin{align*}
\begin{pmatrix}
a & b \\ c & d
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} \lambda & 0 \\ 0 & \lambda \end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
\end{align*} ( a c b d ) ( x y ) = ( λ 0 0 λ ) ( x y )
左辺から右辺を引いて
( a b c d ) ( x y ) − ( λ 0 0 λ ) ( x y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
a & b \\ c & d
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
-\begin{pmatrix}
\lambda & 0 \\ 0 & \lambda
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( a c b d ) ( x y ) − ( λ 0 0 λ ) ( x y ) = ( 0 0 )
( x y ) \begin{pmatrix} x \\ y \end{pmatrix} ( x y ) でくくって
{ ( a b c d ) − ( λ 0 0 λ ) } ( x y ) = ( 0 0 ) \begin{align*}
\left\{\begin{pmatrix}
a & b \\ c & d
\end{pmatrix}
-\begin{pmatrix} \lambda & 0 \\ 0 & \lambda \end{pmatrix}\right\}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} { ( a c b d ) − ( λ 0 0 λ ) } ( x y ) = ( 0 0 )
行列のひき算を計算して
( a − λ b − 0 c − 0 d − λ ) ( x y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
a-\lambda & b-0 \\ c-0 & d-\lambda
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( a − λ c − 0 b − 0 d − λ ) ( x y ) = ( 0 0 )
より
( a − λ b c d − λ ) ( x y ) = ( 0 0 ) ⋯ ( ∗ 3 ) \begin{align*}
\begin{pmatrix}
a-\lambda & b \\ c & d-\lambda
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix} \cdots(*3)
\end{align*} ( a − λ c b d − λ ) ( x y ) = ( 0 0 ) ⋯ ( ∗ 3 )
このとき、固有ベクトルに掛けられている行列 ( a − λ b c d − λ ) \begin{pmatrix} a-\lambda & b \\ c & d-\lambda \end{pmatrix} ( a − λ c b d − λ ) の行列式
∣ a − λ b c d − λ ∣ = ( a − λ ) ( d − λ ) − b c \begin{align*}
\begin{vmatrix}
a-\lambda & b \\ c & d-\lambda
\end{vmatrix}=(a-\lambda)(d-\lambda)-bc
\end{align*} a − λ c b d − λ = ( a − λ ) ( d − λ ) − b c
について、次のことが成り立ちます。
(case1) ∣ a − λ b c d − λ ∣ ≠ 0 \begin{vmatrix} a-\lambda & b \\ c & d-\lambda \end{vmatrix}\ne0 a − λ c b d − λ = 0 のとき
このときは、( x y ) = ( 0 0 ) \begin{pmatrix} x \\ y \end{pmatrix}=\begin{pmatrix} 0 \\ 0 \end{pmatrix} ( x y ) = ( 0 0 ) となるので、A A A には固有値、および固有ベクトルは存在しません。
(case2) ∣ a − λ b c d − λ ∣ = 0 \begin{vmatrix} a-\lambda & b \\ c & d-\lambda \end{vmatrix}=0 a − λ c b d − λ = 0 のとき
このときは、A A A には固有値、および固有ベクトルは存在し、この
∣ a − λ b c d − λ ∣ = ( a − λ ) ( d − λ ) − b c = 0 \begin{align*}
\begin{vmatrix} a-\lambda & b \\ c & d-\lambda \end{vmatrix}=(a-\lambda)(d-\lambda)-bc=0
\end{align*} a − λ c b d − λ = ( a − λ ) ( d − λ ) − b c = 0
を解くことによって、固有値 λ \lambda λ が求まります。この
( a − λ ) ( d − λ ) − b c = 0 \begin{align*}
(a-\lambda)(d-\lambda)-bc=0
\end{align*} ( a − λ ) ( d − λ ) − b c = 0
を特性方程式 、または固有方程式 といいます。 この特性方程式を解いて固有値が求まれば、それを元々の ( ∗ 3 ) (*3) ( ∗ 3 ) 式
( a − λ b c d − λ ) ( x y ) = ( 0 0 ) ⋯ ( ∗ 3 ) \begin{align*}
\begin{pmatrix}
a-\lambda & b \\ c & d-\lambda
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix} \cdots(*3)
\end{align*} ( a − λ c b d − λ ) ( x y ) = ( 0 0 ) ⋯ ( ∗ 3 )
に代入し、この(行列で表された)連立一次方程式を解くことによって、固有ベクトル ( x y ) \begin{pmatrix} x \\ y \end{pmatrix} ( x y ) が求まります。
以上のことをまとめると、次の必要十分条件が成り立ちます。
<固有値をもつための必要十分条件> ある行列 A = ( a b c d ) A=\begin{pmatrix} a & b \\ c & d \end{pmatrix} A = ( a c b d ) に固有値 λ \lambda λ 、および固有ベクトルが存在するための必要十分条件は
( a − λ b c d − λ ) \begin{align*}
\begin{pmatrix} a-\lambda & b \\ c & d-\lambda \end{pmatrix}
\end{align*} ( a − λ c b d − λ )
の行列式が 0 0 0 になること、つまり
∣ a − λ b c d − λ ∣ = ( a − λ ) ( d − λ ) − b c = 0 \begin{align*}
\begin{vmatrix} a-\lambda & b \\ c & d-\lambda \end{vmatrix}=(a-\lambda)(d-\lambda)-bc=0
\end{align*} a − λ c b d − λ = ( a − λ ) ( d − λ ) − b c = 0
である。
***
証明しましょう。 まず、(case1) の命題
∣ a − λ b c d − λ ∣ ≠ 0 \begin{vmatrix} a-\lambda & b \\ c & d-\lambda \end{vmatrix}\ne0 a − λ c b d − λ = 0 → ならば \xrightarrow[ならば ]{} ならば 固有値 λ \lambda λ が存在しない
について、元の命題が真であれば、その対偶も真となるので、(case1) の対偶
固有値 λ \lambda λ が存在する → ならば \xrightarrow[ならば ]{} ならば ∣ a − λ b c d − λ ∣ = 0 \begin{vmatrix} a-\lambda & b \\ c & d-\lambda \end{vmatrix}=0 a − λ c b d − λ = 0
も真となります(補足) 。
次に、∣ a − λ b c d − λ ∣ = 0 \begin{vmatrix} a-\lambda & b \\ c & d-\lambda \end{vmatrix}=0 a − λ c b d − λ = 0 を仮定すると、(case2) より、固定値 λ \lambda λ 、および固有ベクトルが存在するので
固有値 λ \lambda λ が存在する ← ならば \xleftarrow[ならば ]{} ならば ∣ a − λ b c d − λ ∣ = 0 \begin{vmatrix} a-\lambda & b \\ c & d-\lambda \end{vmatrix}=0 a − λ c b d − λ = 0
よって、”→ \rightarrow → ” と ”← \leftarrow ← ” の両方の向きの ”ならば” が真であるので
固有値 λ \lambda λ が存在する ⇄ \rightleftarrows ⇄ ∣ a − λ b c d − λ ∣ = 0 \begin{vmatrix} a-\lambda & b \\ c & d-\lambda \end{vmatrix}=0 a − λ c b d − λ = 0
となり、これで固有値 λ \lambda λ 、および固有ベクトルが存在するための必要十分条件は
∣ a − λ b c d − λ ∣ = 0 \begin{align*}
\begin{vmatrix} a-\lambda & b \\ c & d-\lambda \end{vmatrix}=0
\end{align*} a − λ c b d − λ = 0
であることが証明されます。
では、この<固有値をもつための必要十分条件> を用いて、固有値と固有ベクトルを求めてみましょう。
例題4 行列 A = ( 1 2 2 4 ) A=\begin{pmatrix} 1 & 2 \\ 2 & 4 \end{pmatrix} A = ( 1 2 2 4 ) の固有値と固有ベクトルを求めなさい。
指針 まずは、この行列の特性方程式を解いて固有値を求めます。そして、その固有値を(行列で表された)連立一次方程式に代入し、それを解くことで固有ベクトルが得られます。
解答 2次正方行列 A = ( 1 2 2 4 ) A=\begin{pmatrix} 1 & 2 \\ 2 & 4 \end{pmatrix} A = ( 1 2 2 4 ) に対して、固有ベクトルを ( x y ) ( ≠ 0 ⃗ ) \begin{pmatrix} x \\ y \end{pmatrix}(\ne\vec{0}) ( x y ) ( = 0 ) 、固有値を λ \lambda λ とすると、以下の関係が成り立ちます。
( 1 2 2 4 ) ( x y ) = λ ( x y ) \begin{align*}
\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\lambda
\begin{pmatrix}
x \\ y
\end{pmatrix}
\end{align*} ( 1 2 2 4 ) ( x y ) = λ ( x y )
これより
( 1 2 2 4 ) ( x y ) = λ ( x y ) ↔ ( 1 2 2 4 ) ( x y ) = λ ( 1 0 0 1 ) ( x y ) ↔ ( 1 2 2 4 ) ( x y ) = ( λ 0 0 λ ) ( x y ) ↔ ( 1 2 2 4 ) ( x y ) − ( λ 0 0 λ ) ( x y ) = ( 0 0 ) ↔ { ( 1 2 2 4 ) − ( λ 0 0 λ ) } ( x y ) = ( 0 0 ) ↔ ( 1 − λ 2 2 4 − λ ) ( x y ) = ( 0 0 ) ⋯ ( ∗ 4 ) \begin{align*}
&\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\lambda
\begin{pmatrix}
x \\ y
\end{pmatrix}\\
\leftrightarrow\,&\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\lambda
\begin{pmatrix}
1 & 0 \\ 0 & 1
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}\\
\leftrightarrow\,&\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix}
\lambda & 0 \\ 0 & \lambda
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}\\
\leftrightarrow\,&\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
-\begin{pmatrix}
\lambda & 0 \\ 0 & \lambda
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix}
0 \\ 0
\end{pmatrix}\\
\leftrightarrow\,&\left\{\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
-\begin{pmatrix} \lambda & 0 \\ 0 & \lambda \end{pmatrix}\right\}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}\\
\leftrightarrow\,&\begin{pmatrix}
1-\lambda & 2 \\ 2 & 4-\lambda
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix} \cdots(*4)
\end{align*} ↔ ↔ ↔ ↔ ↔ ( 1 2 2 4 ) ( x y ) = λ ( x y ) ( 1 2 2 4 ) ( x y ) = λ ( 1 0 0 1 ) ( x y ) ( 1 2 2 4 ) ( x y ) = ( λ 0 0 λ ) ( x y ) ( 1 2 2 4 ) ( x y ) − ( λ 0 0 λ ) ( x y ) = ( 0 0 ) { ( 1 2 2 4 ) − ( λ 0 0 λ ) } ( x y ) = ( 0 0 ) ( 1 − λ 2 2 4 − λ ) ( x y ) = ( 0 0 ) ⋯ ( ∗ 4 )
ここでは解説用に計算しましたが、最初から特性方程式を考えればよいので、慣れたら ( ∗ 4 ) (*4) ( ∗ 4 ) から計算を始めてもかまいません。 すると、特性方程式
∣ 1 − λ 2 2 4 − λ ∣ = 0 \begin{align*}
\begin{vmatrix} 1-\lambda & 2 \\ 2 & 4-\lambda \end{vmatrix}=0
\end{align*} 1 − λ 2 2 4 − λ = 0
より
( 1 − λ ) ( 4 − λ ) − 2 2 = 0 ( 1 − λ ) ( 4 − λ ) − 4 = 0 \begin{align*}
&(1-\lambda)(4-\lambda)-2^2=0\\
&(1-\lambda)(4-\lambda)-4=0
\end{align*} ( 1 − λ ) ( 4 − λ ) − 2 2 = 0 ( 1 − λ ) ( 4 − λ ) − 4 = 0
展開して整理すると
1 ⋅ 4 − λ − 4 λ + λ 2 − 4 = 0 4 − 5 λ + λ 2 − 4 = 0 λ 2 − 5 λ = 0 λ ( λ − 5 ) = 0 \begin{align*}
1\cdot4-\lambda-4\lambda+\lambda^2-4&=0\\
4-5\lambda+\lambda^2-4&=0\\
\lambda^2-5\lambda&=0\\
\lambda(\lambda-5)&=0
\end{align*} 1 ⋅ 4 − λ − 4 λ + λ 2 − 4 4 − 5 λ + λ 2 − 4 λ 2 − 5 λ λ ( λ − 5 ) = 0 = 0 = 0 = 0
よって
λ = 0 または λ − 5 = 0 \begin{align*}
\lambda=0 または \lambda-5=0
\end{align*} λ = 0 または λ − 5 = 0
より
λ = 0 または λ = 5 \begin{align*}
\lambda=0 または \lambda=5
\end{align*} λ = 0 または λ = 5
(ⅰ) λ = 0 \lambda=0 λ = 0 のとき
( ∗ 4 ) (*4) ( ∗ 4 ) に代入して
( 1 − 0 2 3 4 − 0 ) ( x y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
1-0 & 2 \\ 3 & 4-0
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( 1 − 0 3 2 4 − 0 ) ( x y ) = ( 0 0 )
より
( 1 2 2 4 ) ( x y ) = ( 0 0 ) ( x + 2 y 2 x + 4 y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
&=\begin{pmatrix} 0 \\ 0 \end{pmatrix}\\
\begin{pmatrix}
x+2y \\ 2x+4y
\end{pmatrix}
&=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( 1 2 2 4 ) ( x y ) ( x + 2 y 2 x + 4 y ) = ( 0 0 ) = ( 0 0 )
対応する成分は等しいので
{ x + 2 y = 0 2 x + 4 y = 0 \begin{cases}
x+2y=0\\
2x+4y=0
\end{cases} { x + 2 y = 0 2 x + 4 y = 0
第2式の両辺を 2 2 2 で割ると
{ x + 2 y = 0 x + 2 y = 0 \begin{cases}
x+2y=0\\
x+2y=0
\end{cases} { x + 2 y = 0 x + 2 y = 0
この(同じ式の)連立方程式の解は x + 2 y = 0 x+2y=0 x + 2 y = 0 を満たすすべての実数となるので、C C C を任意の実数として y = − C y=-C y = − C とすると
x + 2 ( − C ) = 0 x − 2 C = 0 \begin{align*}
x+2(-C)&=0\\
x-2C&=0
\end{align*} x + 2 ( − C ) x − 2 C = 0 = 0
より
x = 2 C x=2C x = 2 C
したがって、この連立方程式の解は、C C C を任意の実数として
( x y ) = ( 2 C − C ) = C ( 2 − 1 ) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}&=\begin{pmatrix*}[c] 2C \\ -C \end{pmatrix*}\\
&=C\begin{pmatrix*}[r] 2 \\ -1 \end{pmatrix*}
\end{align*} ( x y ) = ( 2 C − C ) = C ( 2 − 1 )
よって、λ = 0 \lambda=0 λ = 0 のときの固有ベクトルは、最も簡単な C = 1 C=1 C = 1 のときを選んで
( x y ) = ( 2 − 1 ) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}=\begin{pmatrix*}[r] 2 \\ -1 \end{pmatrix*}
\end{align*} ( x y ) = ( 2 − 1 )
(ⅱ) λ = 5 \lambda=5 λ = 5 のとき
( ∗ 4 ) (*4) ( ∗ 4 ) に代入して
( 1 − 5 2 3 4 − 5 ) ( x y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
1-5 & 2 \\ 3 & 4-5
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( 1 − 5 3 2 4 − 5 ) ( x y ) = ( 0 0 )
より
( − 4 2 2 − 1 ) ( x y ) = ( 0 0 ) ( − 4 x + 2 y 2 x − y ) = ( 0 0 ) \begin{align*}
\begin{pmatrix}
-4 & 2 \\ 2 & -1
\end{pmatrix}
\begin{pmatrix}
x \\ y
\end{pmatrix}
&=\begin{pmatrix} 0 \\ 0 \end{pmatrix}\\
\begin{pmatrix}
-4x+2y \\ 2x-y
\end{pmatrix}
&=\begin{pmatrix} 0 \\ 0 \end{pmatrix}
\end{align*} ( − 4 2 2 − 1 ) ( x y ) ( − 4 x + 2 y 2 x − y ) = ( 0 0 ) = ( 0 0 )
対応する成分は等しいので
{ − 4 x + 2 y = 0 2 x − y = 0 \begin{cases}
-4x+2y =0\\
2x-y=0
\end{cases} { − 4 x + 2 y = 0 2 x − y = 0
第1式の両辺を − 2 -2 − 2 で割ると
{ 2 x − y = 0 2 x − y = 0 \begin{cases}
2x-y=0\\
2x-y=0
\end{cases} { 2 x − y = 0 2 x − y = 0
この連立方程式の解は 2 x − y = 0 2x-y=0 2 x − y = 0 を満たすすべての実数となるので、C C C を任意の実数として x = C x=C x = C とすると
2 C − y = 0 2C-y=0 2 C − y = 0
より
y = 2 C y=2C y = 2 C
したがって、この連立方程式の解は、C C C を任意の実数として
( x y ) = ( C 2 C ) = C ( 1 2 ) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}&=\begin{pmatrix*}[c] C \\ 2C \end{pmatrix*}\\
&=C\begin{pmatrix*}[r] 1 \\ 2 \end{pmatrix*}
\end{align*} ( x y ) = ( C 2 C ) = C ( 1 2 )
よって、λ = 5 \lambda=5 λ = 5 のときの固有ベクトルは、最も簡単な C = 1 C=1 C = 1 のときを選んで
( x y ) = ( 1 2 ) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}=\begin{pmatrix*}[r] 1 \\ 2 \end{pmatrix*}
\end{align*} ( x y ) = ( 1 2 )
以上より、
固有値が λ = 0 \lambda=0 λ = 0 のときの固有ベクトルは
( x y ) = ( 2 − 1 ) ⋯ (答) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}&=\begin{pmatrix*}[r] 2 \\ -1 \end{pmatrix*} \cdots(答)
\end{align*} ( x y ) = ( 2 − 1 ) ⋯ (答)
固有値が λ = 5 \lambda=5 λ = 5 のときの固有ベクトルは
( x y ) = ( 1 2 ) ⋯ (答) \begin{align*}
\begin{pmatrix*}[r] x \\ y \end{pmatrix*}&=\begin{pmatrix*}[r] 1 \\ 2 \end{pmatrix*} \cdots(答)
\end{align*} ( x y ) = ( 1 2 ) ⋯ (答)
■ \blacksquare ■
結果をまとめると、行列 A = ( 1 2 2 4 ) A=\begin{pmatrix} 1 & 2 \\ 2 & 4 \end{pmatrix} A = ( 1 2 2 4 ) について、以下の2つの式が成り立ちます(要確認)。
( 1 2 2 4 ) ( 2 − 1 ) = 0 ( 2 − 1 ) ( 1 2 2 4 ) ( 1 2 ) = 5 ( 1 2 ) \begin{align*}
&\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
\begin{pmatrix}
2 \\ -1
\end{pmatrix}
=0
\begin{pmatrix}
2 \\ -1
\end{pmatrix}\\[10pt]
&\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
\begin{pmatrix}
1 \\ 2
\end{pmatrix}
=5
\begin{pmatrix}
1 \\ 2
\end{pmatrix}
\end{align*} ( 1 2 2 4 ) ( 2 − 1 ) = 0 ( 2 − 1 ) ( 1 2 2 4 ) ( 1 2 ) = 5 ( 1 2 )
前者より、固有値は 0 0 0 、固有ベクトルは ( 2 − 1 ) \begin{pmatrix*}[r] 2 \\ -1 \end{pmatrix*} ( 2 − 1 ) となります。 後者より、固有値は 5 5 5 、固有ベクトルは ( 1 2 ) \begin{pmatrix*}[r] 1 \\ 2 \end{pmatrix*} ( 1 2 ) となります。
(注意2)固有値は 0 になる場合がある 固有ベクトルは零ベクトルを含めませんが、固有値は 0 0 0 になる場合があります。そのような場合にも、固有値 0 0 0 に対応する固有ベクトルは必ず存在します。 固有値は、その固有ベクトルの大きさ(長さ)を何倍にするかを表しますが、固有値が 0 0 0 のときは、その大きさを「0 0 0 倍にする」場合だといえます。すなわち、その固有ベクトルの方向は、変換行列によってベクトルの大きさが 0 0 0 になる特別な方向だと考えることができます。言い換えると、「その方向のどのようなベクトルも、変換によって必ず零ベクトルにしてしまう」特別な方向だということです。 例えば、例題4 の
( 1 2 2 4 ) ( 2 − 1 ) = 0 ( 2 − 1 ) \begin{align*}
&\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
\begin{pmatrix}
2 \\ -1
\end{pmatrix}
=0
\begin{pmatrix}
2 \\ -1
\end{pmatrix}
\end{align*} ( 1 2 2 4 ) ( 2 − 1 ) = 0 ( 2 − 1 )
のケースを考えると、変換行列 ( 1 2 2 4 ) \begin{pmatrix} 1 & 2 \\ 2 & 4 \end{pmatrix} ( 1 2 2 4 ) によって、( 2 − 1 ) \begin{pmatrix} 2 \\ -1 \end{pmatrix} ( 2 − 1 ) と同じ向きの任意のベクトル(大きさは何でもよい)は、零ベクトルになることが分かります。このことは、逆向きの任意のベクトルについてもいえます。
例 ( 2 − 1 ) \begin{pmatrix} 2 \\ -1 \end{pmatrix} ( 2 − 1 ) と同じ向きのベクトル ( 3 − 1.5 ) \begin{pmatrix} 3 \\ -1.5 \end{pmatrix} ( 3 − 1.5 ) 、または逆向きのベクトル ( − 4 2 ) \begin{pmatrix} -4 \\ 2 \end{pmatrix} ( − 4 2 ) について、次の式が成り立ちます(要確認)。
( 1 2 2 4 ) ( 3 − 1.5 ) = 0 ( 3 − 1.5 ) ( 1 2 2 4 ) ( − 4 2 ) = 0 ( − 4 2 ) \begin{align*}
&\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
\begin{pmatrix}
3 \\ -1.5
\end{pmatrix}
=0
\begin{pmatrix}
3 \\ -1.5
\end{pmatrix}\\[10pt]
&\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
\begin{pmatrix}
-4 \\ 2
\end{pmatrix}
=0
\begin{pmatrix}
-4 \\ 2
\end{pmatrix}
\end{align*} ( 1 2 2 4 ) ( 3 − 1.5 ) = 0 ( 3 − 1.5 ) ( 1 2 2 4 ) ( − 4 2 ) = 0 ( − 4 2 )
一般に、k k k を 0 0 0 以外の任意の実数として、( 2 − 1 ) \begin{pmatrix} 2 \\ -1 \end{pmatrix} ( 2 − 1 ) を k k k 倍したベクトル k ( 2 − 1 ) = ( 2 k − k ) k\begin{pmatrix} 2 \\ -1 \end{pmatrix}=\begin{pmatrix} 2k \\ -k \end{pmatrix} k ( 2 − 1 ) = ( 2 k − k ) について、次の式はつねに成り立ちます(要確認)。
( 1 2 2 4 ) ( 2 k − k ) = 0 ( 2 k − k ) \begin{align*}
&\begin{pmatrix}
1 & 2 \\ 2 & 4
\end{pmatrix}
\begin{pmatrix}
2k \\ -k
\end{pmatrix}
=0
\begin{pmatrix}
2k \\ -k
\end{pmatrix}
\end{align*} ( 1 2 2 4 ) ( 2 k − k ) = 0 ( 2 k − k )
固有値が 0 0 0 のときは、そのときの固有ベクトルの大きさを「0 0 0 倍にする」、つまり、その方向の大きさを「完全につぶしてしまう」特別な場合だと考えられます。
ケイリー・ハミルトンの定理 まずは、先にやった固有値をもつための必要十分条件を確認しましょう。
<固有値をもつための必要十分条件> ある行列 A = ( a b c d ) A=\begin{pmatrix} a & b \\ c & d \end{pmatrix} A = ( a c b d ) に固有値 λ \lambda λ 、および固有ベクトルが存在するための必要十分条件は
( a − λ b c d − λ ) \begin{align*}
\begin{pmatrix} a-\lambda & b \\ c & d-\lambda \end{pmatrix}
\end{align*} ( a − λ c b d − λ )
の行列式が 0 0 0 になること、つまり
∣ a − λ b c d − λ ∣ = ( a − λ ) ( d − λ ) − b c = 0 \begin{align*}
\begin{vmatrix} a-\lambda & b \\ c & d-\lambda \end{vmatrix}=(a-\lambda)(d-\lambda)-bc=0
\end{align*} a − λ c b d − λ = ( a − λ ) ( d − λ ) − b c = 0
である。この方程式を特定方程式 、または固有方程式 という。
***
この特性方程式を展開して、λ \lambda λ について整理すると
( a − λ ) ( d − λ ) − b c = 0 a d − a λ − d λ + λ 2 − b c = 0 λ 2 − a λ − d λ + a d − b c = 0 λ 2 − ( a + b ) λ + a d − b c = 0 \begin{align*}
(a-\lambda)(d-\lambda)-bc=0\\
ad-a\lambda-d\lambda+\lambda^2-bc=0\\
\lambda^2-a\lambda-d\lambda+ad-bc=0\\
\lambda^2-(a+b)\lambda+ad-bc=0
\end{align*} ( a − λ ) ( d − λ ) − b c = 0 a d − aλ − d λ + λ 2 − b c = 0 λ 2 − aλ − d λ + a d − b c = 0 λ 2 − ( a + b ) λ + a d − b c = 0
この2次方程式を解くことによって、固有値 λ \lambda λ が求まります。 この λ \lambda λ を、もとの行列 A A A に置き換えた次の式
A 2 − ( a + b ) A + ( a d − b c ) E = O \begin{align*}
A^2-(a+b)A+(ad-bc)E=O
\end{align*} A 2 − ( a + b ) A + ( a d − b c ) E = O
が、かつて高校でも習ったケーリーハミルトンの定理 とよばれる関係式です。この定理の証明自体は、中学でも分かる⑬「行列」で行いました。 なお、a d − b c ad-bc a d − b c に、2次の単位行列 E E E を掛けていることに注意しましょう。E E E を掛ける理由は、行列として扱うためです。つまり、この式の行列を成分で表すと
( a b c d ) 2 − ( a + b ) ( a b c d ) + ( a d − b c ) ( 1 0 0 1 ) = ( 0 0 0 0 ) \begin{align*}
\begin{pmatrix} a & b \\ c & d \end{pmatrix}^2-(a+b)\begin{pmatrix} a & b \\ c & d \end{pmatrix}+(ad-bc)\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}=\begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix}
\end{align*} ( a c b d ) 2 − ( a + b ) ( a c b d ) + ( a d − b c ) ( 1 0 0 1 ) = ( 0 0 0 0 )
となります。仮に E E E を掛けないと
( a b c d ) 2 − ( a + b ) ( a b c d ) + ( a d − b c ) = ( 0 0 0 0 ) \begin{align*}
\begin{pmatrix} a & b \\ c & d \end{pmatrix}^2-(a+b)\begin{pmatrix} a & b \\ c & d \end{pmatrix}+(ad-bc)=\begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix}
\end{align*} ( a c b d ) 2 − ( a + b ) ( a c b d ) + ( a d − b c ) = ( 0 0 0 0 )
と、行列に実数 a d − b c ad-bc a d − b c をたすことになり、意味のない式になります。 また、O O O は零行列 ( 0 0 0 0 ) \begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix} ( 0 0 0 0 ) であって、実数値の 0 0 0 ではないことにも注意しましょう。左辺全体は2行2列の行列になるので、右辺も同じ型の、2行2列の零行列でなければなりません。
このケイリー・ハミルトンの定理は、次のように n n n 次正方行列に拡張されます。この一般化された主張が、線形代数における重要な「ケイリー・ハミルトンの定理」 です。
<ケイリー・ハミルトンの定理(一般)> n × n n\times n n × n 行列 A A A の特性多項式 P A ( λ ) P_A(λ) P A ( λ ) が、次のように与えられたとします。これは固有多項式ともよばれます。
P A ( λ ) = ∣ A − λ E ∣ = c n λ n + c n − 1 λ n − 1 + ⋯ + c 1 λ + c 0 \begin{align*}
P_A(λ)=|A−λE|=c_nλ^n+c_{n−1}λ^{n−1}+\cdots+c_1λ+c_0
\end{align*} P A ( λ ) = ∣ A − λ E ∣ = c n λ n + c n − 1 λ n − 1 + ⋯ + c 1 λ + c 0
c 1 , c 2 , c 3 , ⋯ , c n c_1, c_2, c_3, \cdots, c_n c 1 , c 2 , c 3 , ⋯ , c n は、行列式の計算によって得られる、実数の値をもつ係数です。E E E は n × n n\times n n × n の単位行列です。 このとき、ケイリー・ハミルトンの定理は、 「特性多項式 P A ( λ ) P_A(λ) P A ( λ ) の固有値 λ λ λ に行列 A A A を代入すると、零行列 O O O になる」 ことを主張します。つまり、以下の式が成り立ちます。
P A ( A ) = c n A n + c n − 1 A n − 1 + ⋯ + c 1 A + c 0 E = O ⋯ ( ∗ 5 ) \begin{align*}
P_A(A)=c_nA^n+c_{n−1}A^{n−1}+⋯+c_1A+c_0E=O \cdots{(*5)}
\end{align*} P A ( A ) = c n A n + c n − 1 A n − 1 + ⋯ + c 1 A + c 0 E = O ⋯ ( ∗ 5 )
なお、定数項 c 0 c_0 c 0 には、行列にするために n × n n\times n n × n の単位行列 E E E を掛けます。また、O O O は n × n n\times n n × n の零行列です。
この ( ∗ 5 ) (*5) ( ∗ 5 ) について、n = 2 n=2 n = 2 として、c 2 = a + d , c 3 = a d − b c c_2=a+d, c_3=ad-bc c 2 = a + d , c 3 = a d − b c としたものが、2次正方行列におけるケイリー・ハミルトンの定理です。
A 2 − ( a + b ) A + ( a d − b c ) E = O \begin{align*}
A^2-(a+b)A+(ad-bc)E=O
\end{align*} A 2 − ( a + b ) A + ( a d − b c ) E = O
これは3次正方行列、4次正方行列、さらには一般の n n n 次正方行列まで拡張されます。拡張するにあたって、各 c n c_n c n の式は複雑なものになりますが、その過程でさまざまな定理や性質が明らかになります。 その内容はこの記事の範囲を超えるので、より深く学びたい場合は、大学の線形代数の教科書などをご参照ください。
固有値の歴史 フランスの数学者、オーギュスタン=ルイ・コーシー(Augustin Louis Cauchy)は、固有値の研究において重要な役割を果たした数学者の一人です。
Augustin Louis Cauchy(1789年 - 1857年) 19世紀初頭、道具としての「行列 」はすでに存在していました。しかし、行列の性質やその応用など、行列そのものを研究の対象とする、抽象的な理論としての「行列理論 」は、まだ確立されてはいませんでした。 その中で彼は、「二次形式の標準化 」という問題に取り組む中で、行列の「固有値 」や「固有ベクトル 」の重要性に着目し、それを深く研究しました。
コーシーがとりわけ関心を寄せていたのは、A 2 + B x y + C y 2 = D A^2+Bxy+Cy^2=D A 2 + B x y + C y 2 = D のような二次形式とよばれる式を、座標変換によって A ′ x 2 + C ′ y 2 = D A'x^2+C'y^2=D A ′ x 2 + C ′ y 2 = D という単純な形、いわゆる標準形に変換することでした。この操作を、二次形式の標準化といいます。 コーシーは、この標準化による変換後の係数 A ′ , C ′ A', C' A ′ , C ′ が、ある特定の方程式、今日の用語でいえば「特性方程式 」の解として得られることを発見しました。 この解こそが、現代でいう「固有値 」に他なりません。 さらにコーシーは、この固有値が、物体がどのくらい速く振動するか、どのくらい回転しやすいかといった、力学における運動の特徴を表していることにも気づいていました。 固有値や固有ベクトルは、行列の性質をとらえる上で非常に重要な概念です。コーシーは二次形式の標準化を研究する中で、行列の背後にある構造を明らかにし、その性質を引き出すための手法を切り開いていたのです。
コーシーの研究を土台として、ケイリーやシルベスター、ハミルトンといった後続の数学者たちによって、行列理論の本格的な研究が加速されていきました。そしてそれは今日の「線形代数 」という学問へと受け継がれていきます。コーシーがいなければ、現在の線形代数の発展は遅れていたといっても過言ではないでしょう。
(補足)必要十分条件と対偶 必要十分条件と対偶について、簡単に解説します。詳しくは高校1年で習います。
命題の真と偽、反例 2つの条件 p , q p, q p , q に対して「p p p ならば q q q 」という命題を考えます。これは「ならば」を意味する記号 → \rightarrow → や ⇒ \Rightarrow ⇒ を用いて
p → q \begin{align*}
p \rightarrow q
\end{align*} p → q
または
p ⇒ q \begin{align*}
p \Rightarrow q
\end{align*} p ⇒ q
と表します。このとき、p p p をこの命題の仮定 、q q q を結論 といいます。 例えば「x = 1 x=1 x = 1 ならば x 2 = 1 x^2=1 x 2 = 1 」という命題は
x = 1 → x 2 = 1 \begin{align*}
x=1 \rightarrow x^2=1
\end{align*} x = 1 → x 2 = 1
と表せます。このとき x = 1 x=1 x = 1 は仮定で、x 2 = 1 x^2=1 x 2 = 1 は結論です。
また、命題が正しいとき、その命題は真 であるといい、正しくないとき、その命題は偽 であるといいいます。 例えば、 x = 1 x=1 x = 1 ならば x 2 = 1 x^2=1 x 2 = 1 という命題は真です。x = 1 x=1 x = 1 を x 2 = 1 x^2=1 x 2 = 1 に代入すれば成り立ちます。 一方、その命題の仮定と結論を入れ換えた x 2 = 1 x^2=1 x 2 = 1 ならば x = 1 x=1 x = 1 という命題は、正しくないので偽です。
ある命題が偽であることを示すには、その命題が成り立たない例を一つあげます、そのような例を反例 といいます。反例とは、仮定を満たすが結論を満たさない例のことです。 すると、先ほどの「x 2 = 1 x^2=1 x 2 = 1 ならば x = 1 x=1 x = 1 」の反例は x = − 1 x=-1 x = − 1 となります。x = − 1 x=-1 x = − 1 は、仮定 x 2 = 1 x^2=1 x 2 = 1 を満たしますが、結論 x = 1 x=1 x = 1 を満たしません。
必要条件と十分条件 一般に、2つの条件 p , q p, q p , q に対して、命題 p → q p \rightarrow q p → q が真であるとき、 p p p は q q q であるための十分条件 q q q は p p p であるための必要条件 という。
p 十分条件 → q 必要条件 \begin{align*}
\underset{\tiny 十分条件}{p}\hspace{-6pt}\rightarrow\hspace{-6pt}\underset{\tiny 必要条件}{q}
\end{align*} 十分条件 p → 必要条件 q
***
「必要」や「十分」という日常語の意味から理解しようとすると、多くの場合はかえって混乱するので、上記のように機械的に覚えればよいでしょう。
<覚え方> 「~は」の主語だけに着目して、 ・「~は」が矢印の根元にあれば十分条件 ・「~は」が矢印の先端にあれば必要条件 と覚えます。下図のように、矢印の根元は「十」と書けるので十分条件、矢印の先端は「必」とかけるので必要条件、と覚えればよいでしょう。
例えば、先ほどの「x = 1 x=1 x = 1 ならば x 2 = 1 x^2=1 x 2 = 1 」という命題を考えます。
x = 1 → x 2 = 1 \begin{align*}
x=1 \rightarrow x^2=1
\end{align*} x = 1 → x 2 = 1
この命題は真です。このとき、 x = 1 x=1 x = 1 は x 2 = 1 x^2=1 x 2 = 1 であるための十分条件 であり x 2 = 1 x^2=1 x 2 = 1 は x = 1 x=1 x = 1 であるための必要条件 となります。
x = 1 十分条件 → x 2 = 1 必要条件 \begin{align*}
\underset{\tiny 十分条件}{x=1} \rightarrow \underset{\tiny 必要条件}{x^2=1}
\end{align*} 十分条件 x = 1 → 必要条件 x 2 = 1
十分条件や必要条件は、2つの条件の間の論理的な関係を表す言葉なので、なにかある1つの条件だけを取り出して「~は十分条件である」というような言い方は、論理的には適切ではありません。 ただし、文脈から「何のための条件か」がはっきりしているときには、「p は十分条件である」「q は必要条件である」のように、「~であるための」を省略していうこともよくあります。
例 x = 1 x=1 x = 1 は十分条件 、x 2 = 1 x^2=1 x 2 = 1 は必要条件 である。
次に、この命題の仮定と結論を入れ換えた「x 2 = 1 x^2=1 x 2 = 1 ならば x = 1 x=1 x = 1 」という命題
x 2 = 1 → x = 1 \begin{align*}
x^2=1 \rightarrow x=1
\end{align*} x 2 = 1 → x = 1
は偽でした(反例は x = − 1 x=-1 x = − 1 )。なお、十分条件や必要条件といえるのは、その命題「p → q p \rightarrow q p → q 」が真である場合だけです。命題が偽である場合には、必要条件や十分条件という枠組みは適用できません。
必要十分条件 「p p p ならば q q q 」と「q q q ならば p p p 」の両方の命題を、両向きの「ならば」を意味する記号 ⇄ \rightleftarrows ⇄ や ↔ \leftrightarrow ↔ や ⇔ \Leftrightarrow ⇔ を用いて
p ⇄ q \begin{align*}
p \rightleftarrows q
\end{align*} p ⇄ q
または
p ↔ q \begin{align*}
p \leftrightarrow q
\end{align*} p ↔ q
または
p ⇔ q \begin{align*}
p \Leftrightarrow q
\end{align*} p ⇔ q
と表す。 「p p p ならば q q q 」と「q q q ならば p p p 」がともに真であるとき、すなわち p ↔ q p \leftrightarrow q p ↔ q が成り立つとき、 p p p は q q q であるための必要十分条件 という。この場合、 q q q は p p p であるための必要十分条件 でもある。
p 必要十分 条件 ↔ q 必要十分 条件 \begin{align*}
\underset{必要十分 \atop 条件}{p}\hspace{-7pt}\leftrightarrow\hspace{-7pt}\underset{必要十分 \atop 条件}{q}
\end{align*} 条件 必要十分 p ↔ 条件 必要十分 q
またこのとき、 p p p と q q q は同値 である という。
***
例えば、「x = 0 x=0 x = 0 ならば x 2 = 0 x^2=0 x 2 = 0 」という命題を考えます。
x = 0 → x 2 = 0 \begin{align*}
x=0 \rightarrow x^2=0
\end{align*} x = 0 → x 2 = 0
この命題は真です。x = 0 x=0 x = 0 を x 2 = 0 x^2=0 x 2 = 0 に代入すると成り立ちます。 次に、この命題の仮定を結論を入れ換えた「x 2 = 0 x^2=0 x 2 = 0 ならば x = 0 x=0 x = 0 」という命題を考えます。
x 2 = 0 → x = 0 \begin{align*}
x^2=0 \rightarrow x=0
\end{align*} x 2 = 0 → x = 0
すると、この命題も真です。x 2 = 0 x^2=0 x 2 = 0 を満たす x x x は、x = 0 x=0 x = 0 のみだからです。よって、 x = 0 x=0 x = 0 は x 2 = 0 x^2=0 x 2 = 0 であるための必要十分条件 であり、 x 2 = 0 x^2=0 x 2 = 0 は x = 0 x=0 x = 0 であるための必要十分条件 です。
x = 0 必要十分 条件 ↔ x 2 = 0 必要十分 条件 \begin{align*}
\underset{必要十分 \atop 条件}{x=0} \leftrightarrow \underset{必要十分 \atop 条件}{x^2=0}
\end{align*} 条件 必要十分 x = 0 ↔ 条件 必要十分 x 2 = 0
つまり、 x 2 = 0 x^2=0 x 2 = 0 と x = 0 x=0 x = 0 は同値 となります。
逆・裏・対偶 p → q p \rightarrow q p → q という命題について、仮定 p p p と結論 q q q を入れ換えた命題
q → p \begin{align*}
q \rightarrow p
\end{align*} q → p
を逆 という。 また、仮定 p p p と結論 q q q の両方を否定にした「p p p でないならば q q q でない」という命題
p ‾ → q ‾ \begin{align*}
\overline{p} \rightarrow \overline{q}
\end{align*} p → q
を裏 という。条件 p p p に対して、「p p p でない」という条件を p p p の否定 といい、条件の記号の上に横棒をつけた p ‾ \overline{p} p という記号で表す。
さらに、この裏を逆にした
p ‾ → q ‾ \begin{align*}
\overline{p} \rightarrow \overline{q}
\end{align*} p → q
を対偶 という。
***
以上のことを、まとめておきます。
<逆・裏・対偶> 命題 p → q p \rightarrow q p → q について、 逆は「q → p q \rightarrow p q → p 」 裏は「p ‾ → q ‾ \overline{p} \rightarrow \overline{q} p → q 」 対偶は「p ‾ → q ‾ \overline{p} \rightarrow \overline{q} p → q 」
***
逆と裏は互いに対偶の関係にあり、逆の対偶は裏、裏の対偶は逆です。また、対偶の対偶はもとの命題になります。 以上を図式化すると次のようになります。
例えば、先に述べた「x = 1 x=1 x = 1 ならば x 2 = 1 x^2=1 x 2 = 1 」という命題について 逆は「x 2 = 1 x^2=1 x 2 = 1 ならば x = 1 x=1 x = 1 」 裏は「x ≠ 1 x\ne1 x = 1 ならば x 2 ≠ 1 x^2\ne1 x 2 = 1 」 対偶は「x 2 = 1 x^2=1 x 2 = 1 ならば x = 1 x=1 x = 1 」 となります。
ここで、次のような重要な定理があります。
「命題とその対偶の真偽は一致する。 」
「x = 1 x=1 x = 1 ならば x 2 = 1 x^2=1 x 2 = 1 」という命題について、この命題は真です。この命題の逆、裏、対偶の真偽は次のようになります。
命題「x = 1 x=1 x = 1 ならば x 2 = 1 x^2=1 x 2 = 1 」は真 逆「x 2 = 1 x^2=1 x 2 = 1 ならば x = 1 x=1 x = 1 」は偽 裏「x ≠ 1 x\ne1 x = 1 ならば x 2 ≠ 1 x^2\ne1 x 2 = 1 」は偽 対偶「x 2 ≠ 1 x^2\ne1 x 2 = 1 ならば x ≠ 1 x\ne1 x = 1 」は真
すると、元の命題と対偶はどちらも真であり、同じく対偶の関係にある逆と裏はどちらも偽であり、「もとの命題とその対偶の真偽は一致する」ことが分かります。
対偶による証明 ある命題を証明するとき、その命題を直接証明するよりも、その対偶を証明するほうが簡単な場合があります。そのようなときは、対偶を証明します。
例題 ある整数 n n n について、次の命題を証明しなさい。 「n 2 n^2 n 2 が偶数であるとき、n n n は偶数である。」証明 この命題の対偶をとると 「n n n は偶数でないとき、n 2 n^2 n 2 は偶数ではない。」 偶数ではない数は奇数なので、さらに次のように表せる。 「n n n は奇数のとき、n 2 n^2 n 2 は奇数である。」 まず、n n n は奇数なので、k k k を整数 (⋯ , − 2 , − 1 , 0 , 1 , 2 , ⋯ \cdots, -2, -1, 0, 1, 2, \cdots ⋯ , − 2 , − 1 , 0 , 1 , 2 , ⋯ ) として
n = 2 k + 1 ⋯ ( ∗ ) n=2k+1 \cdots(*) n = 2 k + 1 ⋯ ( ∗ )
とおける。ここで展開公式 ( a + b ) 2 = a 2 + 2 a b + b 2 (a+b)^2=a^2+2ab+b^2 ( a + b ) 2 = a 2 + 2 ab + b 2 を用いると
n 2 = ( 2 k + 1 ) 2 = ( 2 k ) 2 + 2 ⋅ 2 k ⋅ 1 + 1 2 = 4 k 2 + 4 k + 1 = 2 ( 2 k 2 + 2 k ) + 1 = 2 l + 1 (ただし l = 2 k 2 + 2 k とおいた) \begin{align*}
n^2&=(2k+1)^2\\
&=(2k)^2+2\cdot2k\cdot1+1^2\\
&=4k^2+4k+1\\
&=2(2k^2+2k)+1\\
&=2l+1 (ただし l=2k^2+2k とおいた)
\end{align*} n 2 = ( 2 k + 1 ) 2 = ( 2 k ) 2 + 2 ⋅ 2 k ⋅ 1 + 1 2 = 4 k 2 + 4 k + 1 = 2 ( 2 k 2 + 2 k ) + 1 = 2 l + 1 (ただし l = 2 k 2 + 2 k とおいた)
すると、k k k は整数なので、l ( = 2 k 2 + 2 k ) l(=2k^2+2k) l ( = 2 k 2 + 2 k ) も整数となり 2 l + 1 2l+1 2 l + 1 は奇数。よって、n 2 n^2 n 2 は奇数となる。
■ \blacksquare ■
( ∗ ) (*) ( ∗ ) 例えば、奇数 − 5 , − 3 , − 1 , 1 , 3 , 5 , ⋯ -5, -3, -1, 1, 3, 5, \cdots − 5 , − 3 , − 1 , 1 , 3 , 5 , ⋯ は − 5 = 2 × ( − 3 ) + 1 -5=2\times(-3)+1 − 5 = 2 × ( − 3 ) + 1 − 3 = 2 × ( − 2 ) + 1 -3=2\times(-2)+1 − 3 = 2 × ( − 2 ) + 1 − 1 = 2 × ( − 1 ) + 1 -1=2\times(-1)+1 − 1 = 2 × ( − 1 ) + 1 1 = 2 × 0 + 1 1=2\times0+1 1 = 2 × 0 + 1 3 = 2 × 1 + 1 3=2\times1+1 3 = 2 × 1 + 1 5 = 2 × 2 + 1 5=2\times2+1 5 = 2 × 2 + 1 ⋮ \hspace{30pt}\vdots ⋮ より、2 k + 1 2k+1 2 k + 1 (k k k は整数)の形で書けます。
なお、確認として、2 2 2 で割り切れる整数を偶数 といいます。また、2 2 2 で割り切れない整数を奇数 といいます。よって、偶数ではない数は奇数です。 一般に、偶数と奇数は負の場合も含まれます。例えば − 4 -4 − 4 は偶数、− 3 -3 − 3 は奇数です。また 0 0 0 は 2 2 2 で割り切れるので偶数です。
バックナンバー 理数系小説『シュレーディンガーの Φ』