第8章行列と行列式
第7章までの微分積分学では,つねに実数(あるいは1個か2個の変数)が主役だった.本章から始まる第II部「線形代数学」では,主役が行列(matrix)——数を長方形に並べたもの——に変わる.行列は,見た目こそ「数の表」にすぎないが,連立1次方程式,平面や空間の回転・拡大縮小,統計データの分析,量子力学の状態,コンピュータグラフィックスの座標変換など,理系のほぼすべての分野で姿を変えて現れる,驚くほど汎用性の高い道具である.
本章はその出発点として,望月泰英『数学ノート 線形代数学』の最初の8ページ(付録:転置・べき乗のパターン・置換・行列式の定義・余因子・行列式の性質・行列式の幾何学的意味)を丁寧にたどる.とりわけ中心的な問いは「行列式(determinant)とは何か,そして何の役に立つのか」である.行列式は,一見すると「$ad-bc$ のような謎の計算規則」に見えるかもしれない.しかし本章を読み終える頃には,行列式が(1)連立1次方程式が一意に解けるかどうかの判定条件であり(次章で詳しく使う),(2)ベクトルが作る平行四辺形の面積や平行六面体の体積そのものであり,(3)置換という「並べ替え」の符号を数え上げる操作から自然に生まれる量である,という3つの顔を持つことが分かるはずである.
- 行列の基本演算(和・スカラー倍・積)の確認と,転置行列 ${}^{t}A$ の定義
- 列ベクトルに分割して行列を書く技法と,対角行列・上三角行列・回転行列という3種類の行列のべき乗パターン
- 置換(permutation)の定義,逆置換,偶置換・奇置換と符号 $\sgn$
- 行列式の定義(置換の和による一般の $n$ 次の定義)と,$3\times3$ 行列式の展開(サラス(Sarrus)の公式との対応)
- 余因子 $A_{ij}$ の定義と,余因子展開(行または列に関する行列式の展開)
- 行列式の4つの基本性質(スカラー倍・行の入れ替え・行の加法・多重線形性)とそれを使った効率的な計算
- 行列式の幾何学的意味:2次元では平行四辺形の面積,3次元では平行六面体の体積
もとにしたノート:望月泰英『数学ノート 線形代数学』 pp. 1–8.
8.1 行列の基本操作とべき乗のパターン
8.1.1 行列とは何か —— 行列の演算の確認
数を長方形に並べたものを行列(matrix)という.たとえば
$$ A=\begin{bmatrix}1&2&-1\\0&3&2\end{bmatrix} $$は,横に2行,縦に3列並んだ行列である.一般に,横に $m$ 行,縦に $n$ 列並んだ行列を$(m,n)$ 型の行列といい,$i$ 行目・$j$ 列目にある数を$(i,j)$ 成分といって $a_{ij}$ と書く($i$ を行の添字,$j$ を列の添字という.「添字」(そえじ,subscript)とは,文字の右下・右上に小さく書き添える番号のことである).上の例なら $a_{11}=1$,$a_{23}=2$ である.行列は本書では $A,B,C,\ldots$ のように大文字で表す.
定義8.1 行列(matrix)
$m\times n$ 個の数 $a_{11},a_{12},\ldots,a_{mn}$ を,横 $m$ 行・縦 $n$ 列の長方形に並べたもの
$$ A=\begin{bmatrix}a_{11}&a_{12}&\cdots&a_{1n}\\a_{21}&a_{22}&\cdots&a_{2n}\\\vdots&\vdots&&\vdots\\a_{m1}&a_{m2}&\cdots&a_{mn}\end{bmatrix} $$を,$(m,n)$ 型の行列という.とくに $n=1$ のとき(1列しかないとき)$A$ を列ベクトル(column vector),$m=1$ のとき(1行しかないとき)行ベクトル(row vector)という.また $m=n$ のとき$n$ 次正方行列(square matrix)という.
行列どうしの計算は,高校数学で学んだ数ベクトルの計算を素直に拡張したものである.型が同じ2つの行列 $A=(a_{ij})$,$B=(b_{ij})$ に対して,和とスカラー倍を次のように定める.
定義8.2 行列の和とスカラー倍
$A=(a_{ij})$,$B=(b_{ij})$ を同じ $(m,n)$ 型の行列,$\lambda$ を実数(スカラー,scalar)とするとき,
$$ A+B=(a_{ij}+b_{ij}),\qquad \lambda A=(\lambda a_{ij}) $$と定める.すなわち,和は同じ位置の成分どうしを足し,スカラー倍はすべての成分を $\lambda$ 倍する.(型の異なる行列どうしの和は定義しない.)
問題は積(product)である.行列の積は「同じ位置の成分どうしをかける」のではなく,もっと込み入った規則に従う.まず具体的な数値で確認しよう.$2\times3$ 行列 $A$ と $3\times2$ 行列 $B$ を
$$ A=\begin{bmatrix}1&2&-1\\0&3&2\end{bmatrix},\qquad B=\begin{bmatrix}2&0\\1&-1\\3&1\end{bmatrix} $$とする.積 $AB$ の $(1,1)$ 成分は,$A$ の第1行 $(1,2,-1)$ と $B$ の第1列 $\begin{bmatrix}2\\1\\3\end{bmatrix}$ を「同じ順番の成分どうしかけて全部足す」ことで作る:$1\times2+2\times1+(-1)\times3=2+2-3=1$.同様に $(1,2)$ 成分は $A$ の第1行と $B$ の第2列 $\begin{bmatrix}0\\-1\\1\end{bmatrix}$ から $1\times0+2\times(-1)+(-1)\times1=0-2-1=-3$.$(2,1)$ 成分は $A$ の第2行 $(0,3,2)$ と $B$ の第1列から $0\times2+3\times1+2\times3=9$,$(2,2)$ 成分は $A$ の第2行と $B$ の第2列から $0\times0+3\times(-1)+2\times1=-1$.まとめると
$$ AB=\begin{bmatrix}1&-3\\9&-1\end{bmatrix} $$という $2\times2$ 行列になる.「$A$ の行」と「$B$ の列」を組にして計算する,という規則を一般化したのが次の定義である.
定義8.3 行列の積
$A=(a_{ij})$ を $(m,n)$ 型,$B=(b_{jk})$ を $(n,p)$ 型(=$A$ の列の数と $B$ の行の数が一致)とするとき,積 $AB$ を $(m,p)$ 型の行列
$$ AB=(c_{ik}),\qquad c_{ik}=\sum_{j=1}^{n}a_{ij}b_{jk}=a_{i1}b_{1k}+a_{i2}b_{2k}+\cdots+a_{in}b_{nk} $$で定める.すなわち $(i,k)$ 成分は,$A$ の第 $i$ 行と $B$ の第 $k$ 列を「対応する成分どうしかけて足し合わせた」ものである.$A$ の列の数と $B$ の行の数が一致しないときは,積 $AB$ は定義されない.
注意:行列の積は交換法則が成り立たない
実数の掛け算では $ab=ba$ が常に成り立つが,行列の積では一般に $AB\neq BA$ である.実際,上の例の $A$($2\times3$型)と $B$($3\times2$型)に対し,$AB$ は$2\times2$型だが,$BA$ は$3\times3$型になり,そもそも型が違う.型が同じ正方行列どうしでも,$AB$ と $BA$ は一般に一致しない.行列の積を扱うときは,かける順序を勝手に入れ替えてはいけない.
8.1.2 転置行列
行列の行と列を入れ換える操作を考えよう.
定義8.4 転置行列(transpose)
$(m,n)$ 型の行列 $A=(a_{ij})$ の行と列を入れ換えた $(n,m)$ 型の行列を,$A$ の転置行列といい,${}^{t}A$ と書く.すなわち,${}^{t}A$ の $(j,i)$ 成分は $A$ の $(i,j)$ 成分 $a_{ij}$ である.たとえば
\begin{equation} A=\begin{bmatrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\\a_{31}&a_{32}&a_{33}\end{bmatrix} \qquad {}^{t}A=\begin{bmatrix}a_{11}&a_{21}&a_{31}\\a_{12}&a_{22}&a_{32}\\a_{13}&a_{23}&a_{33}\end{bmatrix} \label{eq:8-transpose-def} \end{equation}である($A$ の第1行 $(a_{11},a_{12},a_{13})$ が ${}^{t}A$ の第1列に,$A$ の第2行が ${}^{t}A$ の第2列に,……という具合に,行がそのまま列になる).
正方行列の場合,転置は「左上から右下への対角線(対角成分,$a_{11},a_{22},a_{33}$ のこと)を軸に鏡映しにする」操作だと考えると分かりやすい.対角線上にある成分はそもそも動かないので,転置しても値は変わらない.図8.1にこの様子を示す.
転置には,あとの章でも繰り返し使う便利な性質がある.とくに,2つの行列の積を転置すると,順序が逆になってそれぞれの転置の積になる.これはノートには書かれていないが,行列を扱ううえで非常によく使う事実なので,ここで確認しておく.
定理8.1 転置の性質
(i) ${}^{t}({}^{t}A)=A$(2回転置すると元に戻る).
(ii) $A$ が $(m,n)$ 型,$B$ が $(n,p)$ 型のとき,${}^{t}(AB)={}^{t}B\,{}^{t}A$(積の転置は,転置の積を逆順にしたものに等しい).
証明
(i) 転置は行と列を入れ換える操作だから,2回行えば元の位置に戻る.成分で確認すると,${}^{t}A$ の $(j,i)$ 成分は $A$ の $(i,j)$ 成分 $a_{ij}$ であり,${}^{t}({}^{t}A)$ の $(i,j)$ 成分は ${}^{t}A$ の $(j,i)$ 成分,すなわちまた $a_{ij}$ である.よって ${}^{t}({}^{t}A)=A$.
(ii) $AB$ は $(m,p)$ 型だから,${}^{t}(AB)$ は $(p,m)$ 型.一方 ${}^{t}B$ は $(p,n)$ 型,${}^{t}A$ は $(n,m)$ 型だから,${}^{t}B\,{}^{t}A$ も $(p,m)$ 型で型は合っている.成分を比べる.$AB$ の $(i,k)$ 成分は定義8.3より $\sum_j a_{ij}b_{jk}$ だから,${}^{t}(AB)$ の $(k,i)$ 成分はこれと同じ $\sum_j a_{ij}b_{jk}$ である.他方,${}^{t}B$ の $(k,j)$ 成分は $b_{jk}$,${}^{t}A$ の $(j,i)$ 成分は $a_{ij}$ だから,定義8.3をそのまま適用すると,${}^{t}B\,{}^{t}A$ の $(k,i)$ 成分は $\sum_j b_{jk}a_{ij}=\sum_j a_{ij}b_{jk}$ となり,両辺の $(k,i)$ 成分が一致する.よって ${}^{t}(AB)={}^{t}B\,{}^{t}A$.
(証明終わり)
例題8.1 積の転置 ${}^{t}(AB)={}^{t}B\,{}^{t}A$ の確認
$A=\begin{bmatrix}1&2&-1\\0&3&2\end{bmatrix}$,$B=\begin{bmatrix}2&0\\1&-1\\3&1\end{bmatrix}$ に対して,$AB$,${}^{t}(AB)$,${}^{t}B\,{}^{t}A$ をそれぞれ計算し,定理8.1(ii)が成り立つことを確認せよ.
解答 8.1.1節ですでに計算したように $AB=\begin{bmatrix}1&-3\\9&-1\end{bmatrix}$ である.よって
$$ {}^{t}(AB)=\begin{bmatrix}1&9\\-3&-1\end{bmatrix} $$次に ${}^{t}B\,{}^{t}A$ を計算する.${}^{t}B=\begin{bmatrix}2&1&3\\0&-1&1\end{bmatrix}$($2\times3$型),${}^{t}A=\begin{bmatrix}1&0\\2&3\\-1&2\end{bmatrix}$($3\times2$型)だから,積 ${}^{t}B\,{}^{t}A$ は $2\times2$ 型になる.$(1,1)$ 成分は ${}^{t}B$ の第1行 $(2,1,3)$ と ${}^{t}A$ の第1列 $(1,2,-1)$ から $2\times1+1\times2+3\times(-1)=2+2-3=1$.$(1,2)$ 成分は ${}^{t}B$ の第1行と ${}^{t}A$ の第2列 $(0,3,2)$ から $2\times0+1\times3+3\times2=0+3+6=9$.$(2,1)$ 成分は ${}^{t}B$ の第2行 $(0,-1,1)$ と ${}^{t}A$ の第1列から $0\times1+(-1)\times2+1\times(-1)=-2-1=-3$.$(2,2)$ 成分は ${}^{t}B$ の第2行と ${}^{t}A$ の第2列から $0\times0+(-1)\times3+1\times2=-3+2=-1$.まとめると
$$ {}^{t}B\,{}^{t}A=\begin{bmatrix}1&9\\-3&-1\end{bmatrix} $$となり,たしかに ${}^{t}(AB)={}^{t}B\,{}^{t}A$ が成り立っている.(sympyで検算済み.)
8.1.3 列ベクトルに分解する —— 計算の技法
行列を扱うとき,「列ごとにベクトルとして取り出す」という見方をすると計算の見通しがよくなることがある.次の2つの式の違いを考えよう.$A=\begin{bmatrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\\a_{31}&a_{32}&a_{33}\end{bmatrix}$ の各列を,列ベクトル $\vec{x}_1=\begin{bmatrix}a_{11}\\a_{21}\\a_{31}\end{bmatrix}$,$\vec{x}_2=\begin{bmatrix}a_{12}\\a_{22}\\a_{32}\end{bmatrix}$,$\vec{x}_3=\begin{bmatrix}a_{13}\\a_{23}\\a_{33}\end{bmatrix}$ とおく.このとき,実数 $\lambda_1,\lambda_2,\lambda_3$ を使って作る次の2つの式は,見た目は似ているが,意味も型もまったく違う.
(記法について:ここでの矢印付き $\vec{x}_1,\vec{x}_2,\vec{x}_3$ は「数を縦に並べた列ベクトル」を表す記号として,本章ではこの形を採用する.第II部(線形代数学)では,同じ「列ベクトル」という対象でも,章によって記号の飾り方が異なる場合がある——たとえば第9章以降の一部の章では太字 $\bm{x}$ を用いる.いずれも指している対象($\R^n$ の要素)は同じであり,章をまたぐときの記号の違いは,そのつど本文中で断ることにする.)
注意:行列の成分にベクトルが入り混じっても良い
$[\vec{x}_1\ \ \vec{x}_2\ \ \vec{x}_3]$ のように,列ベクトルを横に3つ並べたものは,そのまま $3\times3$ 行列 $A$ そのものである.このように「行列の成分」の代わりに「列ベクトル」を並べたと考えても,何も矛盾は起きない.
数学ノート:$[\lambda_1\vec{x}_1+\lambda_2\vec{x}_2+\lambda_3\vec{x}_3]$ と $[\lambda_1\vec{x}_1\ \ \lambda_2\vec{x}_2\ \ \lambda_3\vec{x}_3]$ の違い
1つ目の $[\lambda_1\vec{x}_1+\lambda_2\vec{x}_2+\lambda_3\vec{x}_3]$ は,3本の列ベクトルを$\lambda_1,\lambda_2,\lambda_3$倍してから1本のベクトルに足し合わせたものであり,型は $(3,1)$(列ベクトル1本)である.これは,$A$ の列ベクトルを横に並べた行列 $[\vec{x}_1\ \ \vec{x}_2\ \ \vec{x}_3]$($=A$ そのもの)に,縦ベクトル $\begin{bmatrix}\lambda_1\\\lambda_2\\\lambda_3\end{bmatrix}$ を右からかけたものに等しい:
\begin{equation} [\lambda_1\vec{x}_1+\lambda_2\vec{x}_2+\lambda_3\vec{x}_3] = [\vec{x}_1\ \ \vec{x}_2\ \ \vec{x}_3]\begin{bmatrix}\lambda_1\\\lambda_2\\\lambda_3\end{bmatrix}=\begin{bmatrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\\a_{31}&a_{32}&a_{33}\end{bmatrix}\begin{bmatrix}\lambda_1\\\lambda_2\\\lambda_3\end{bmatrix} \label{eq:8-colvec-sum} \end{equation}(右辺が定義8.3の行列の積になっていることを確かめよ.$A$ の第1行 $(a_{11},a_{12},a_{13})$ と縦ベクトル $(\lambda_1,\lambda_2,\lambda_3)$ の積 $a_{11}\lambda_1+a_{12}\lambda_2+a_{13}\lambda_3$ が,まさに左辺のベクトルの第1成分 $\lambda_1 a_{11}+\lambda_2a_{12}+\lambda_3a_{13}$ に一致している.)
2つ目の $[\lambda_1\vec{x}_1\ \ \lambda_2\vec{x}_2\ \ \lambda_3\vec{x}_3]$ は,3本の列ベクトルを$\lambda_1,\lambda_2,\lambda_3$倍して横に並べたままにしたものであり,型は $(3,3)$(元の $A$ と同じ)である.これは,$A$ の右から対角行列 $\begin{bmatrix}\lambda_1&0&0\\0&\lambda_2&0\\0&0&\lambda_3\end{bmatrix}$ をかけたものに等しい:
\begin{equation} [\lambda_1\vec{x}_1\ \ \lambda_2\vec{x}_2\ \ \lambda_3\vec{x}_3] = [\vec{x}_1\ \ \vec{x}_2\ \ \vec{x}_3]\begin{bmatrix}\lambda_1&0&0\\0&\lambda_2&0\\0&0&\lambda_3\end{bmatrix}=\begin{bmatrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\\a_{31}&a_{32}&a_{33}\end{bmatrix}\begin{bmatrix}\lambda_1&0&0\\0&\lambda_2&0\\0&0&\lambda_3\end{bmatrix} \label{eq:8-colvec-diag} \end{equation}(右辺の積の第1列は,$A$ とこの対角行列の第1列 $(\lambda_1,0,0)$ との積だから,$a_{11}\lambda_1+a_{12}\times0+a_{13}\times0=\lambda_1a_{11}$,……となり,ちょうど $\lambda_1\vec{x}_1$ の各成分に一致する.第2列・第3列も同様である.)
このように,行列を列ベクトルに分割して考えると,「まとめて1本にする(左から縦ベクトルをかける)」のか「別々のまま並べる(右から対角行列をかける)」のかが,かける相手の行列の形の違いとしてすっきり整理できる.この見方は,次の8.1.4項で対角行列のべき乗を考えるときにも役立つ.
8.1.4 行列のべき乗のパターン
正方行列 $A$ に対して,$A^2=AA$,$A^3=AAA$,……と,自然数 $n$ 個の $A$ の積を $A^n$ と書く($A^0=E$,単位行列とする.単位行列の記号 $E$ は,第9章以降でも一貫して使う).一般の行列では $A^n$ の規則的な形を求めるのは難しいが,特別な形の $2\times2$ 行列については,$A^n$ が驚くほどきれいなパターンに従う.ここでは3つの代表的なパターンを見る.
公式8.1 対角行列のべき乗
$A=\begin{bmatrix}\alpha&0\\0&\beta\end{bmatrix}$(対角行列,diagonal matrix:対角成分以外がすべて $0$ の正方行列)のとき,任意の自然数 $n$ に対して
\begin{equation} A^{n}=\begin{bmatrix}\alpha^{n}&0\\0&\beta^{n}\end{bmatrix} \label{eq:8-diag-power} \end{equation}導出:$n=2,3$ で確かめてから,数学的帰納法で証明する
まず $n=2,3$ で具体的に計算してみる.
$$ A^{2}=\begin{bmatrix}\alpha&0\\0&\beta\end{bmatrix}\begin{bmatrix}\alpha&0\\0&\beta\end{bmatrix}=\begin{bmatrix}\alpha^2+0&\alpha\cdot0+\beta\cdot0\\\alpha\cdot0+\beta\cdot0&0+\beta^2\end{bmatrix}=\begin{bmatrix}\alpha^{2}&0\\0&\beta^{2}\end{bmatrix} $$ $$ A^{3}=A^{2}\cdot A=\begin{bmatrix}\alpha^2&0\\0&\beta^2\end{bmatrix}\begin{bmatrix}\alpha&0\\0&\beta\end{bmatrix}=\begin{bmatrix}\alpha^3+0&\alpha\cdot0+\beta^2\cdot0\\\alpha^2\cdot0+\beta\cdot0&0+\beta^3\end{bmatrix}=\begin{bmatrix}\alpha^{3}&0\\0&\beta^{3}\end{bmatrix} $$(積の非対角成分がつねに $0\times(\text{何か})$ の形になって消えることに注目する.)たしかに $n=2,3$ でパターン通りになっている.これを一般の $n$ で証明するには,数学的帰納法(mathematical induction,「$n=1$ で成り立ち,$n=k$ で成り立てば $n=k+1$ でも成り立つ」ことを示せば,すべての自然数 $n$ で成り立つとする証明法)を使う.$n=1$ のとき $A^1=A=\begin{bmatrix}\alpha^1&0\\0&\beta^1\end{bmatrix}$ で成り立つ.$n=k$ で $A^{k}=\begin{bmatrix}\alpha^{k}&0\\0&\beta^{k}\end{bmatrix}$ が成り立つと仮定すると,
$$ A^{k+1}=A^{k}A=\begin{bmatrix}\alpha^{k}&0\\0&\beta^{k}\end{bmatrix}\begin{bmatrix}\alpha&0\\0&\beta\end{bmatrix}=\begin{bmatrix}\alpha^{k}\alpha+0\cdot0&\alpha^{k}\cdot0+0\cdot\beta\\0\cdot\alpha+\beta^{k}\cdot0&0\cdot0+\beta^{k}\beta\end{bmatrix}=\begin{bmatrix}\alpha^{k+1}&0\\0&\beta^{k+1}\end{bmatrix} $$となり,$n=k+1$ でも成り立つ.よってすべての自然数 $n$ について式\eqref{eq:8-diag-power}が成り立つ.
(証明終わり)
公式8.2 上三角行列のべき乗
$A=\begin{bmatrix}1&\alpha\\0&1\end{bmatrix}$(上三角行列,upper triangular matrix:対角線より左下がすべて $0$ の正方行列)のとき,任意の自然数 $n$ に対して
\begin{equation} A^{n}=\begin{bmatrix}1&n\alpha\\0&1\end{bmatrix} \label{eq:8-uppertri-power} \end{equation}導出:$n=2,3$ で確かめてから,数学的帰納法で証明する
$$ A^{2}=\begin{bmatrix}1&\alpha\\0&1\end{bmatrix}\begin{bmatrix}1&\alpha\\0&1\end{bmatrix}=\begin{bmatrix}1+\alpha\cdot0&\alpha+\alpha\\1\cdot0+1\cdot0&\alpha\cdot0+1\end{bmatrix}=\begin{bmatrix}1&2\alpha\\0&1\end{bmatrix} $$ $$ A^{3}=A^{2}\cdot A=\begin{bmatrix}1&2\alpha\\0&1\end{bmatrix}\begin{bmatrix}1&\alpha\\0&1\end{bmatrix}=\begin{bmatrix}1+2\alpha\cdot0&\alpha+2\alpha\\1\cdot0+0\cdot1&\alpha\cdot0+1\end{bmatrix}=\begin{bmatrix}1&3\alpha\\0&1\end{bmatrix} $$よって $n=2,3$ ではパターン $A^n=\begin{bmatrix}1&n\alpha\\0&1\end{bmatrix}$ が成り立っている.一般の $n$ を数学的帰納法で示す.$n=1$ のとき $A^1=A=\begin{bmatrix}1&\alpha\\0&1\end{bmatrix}=\begin{bmatrix}1&1\times\alpha\\0&1\end{bmatrix}$ で成り立つ.$n=k$ で $A^{k}=\begin{bmatrix}1&k\alpha\\0&1\end{bmatrix}$ が成り立つと仮定すると,
$$ A^{k+1}=A^{k}A=\begin{bmatrix}1&k\alpha\\0&1\end{bmatrix}\begin{bmatrix}1&\alpha\\0&1\end{bmatrix}=\begin{bmatrix}1\cdot1+k\alpha\cdot0&1\cdot\alpha+k\alpha\cdot1\\0\cdot1+1\cdot0&0\cdot\alpha+1\cdot1\end{bmatrix}=\begin{bmatrix}1&(k+1)\alpha\\0&1\end{bmatrix} $$となり,$n=k+1$ でも成り立つ.よってすべての自然数 $n$ について式\eqref{eq:8-uppertri-power}が成り立つ.
(証明終わり)
公式8.3 回転行列のべき乗
$A=\begin{bmatrix}\cos\theta&-\sin\theta\\\sin\theta&\cos\theta\end{bmatrix}$(回転行列,rotation matrix:平面上のベクトルを原点周りに角 $\theta$ だけ回転させる行列.高校数学の複素数平面や図形と方程式で登場した回転の式と同じものである)のとき,任意の自然数 $n$ に対して
\begin{equation} A^{n}=\begin{bmatrix}\cos n\theta&-\sin n\theta\\\sin n\theta&\cos n\theta\end{bmatrix} \label{eq:8-rotation-power} \end{equation}すなわち,角 $\theta$ の回転を $n$ 回繰り返すことは,角 $n\theta$ の回転1回と同じである.
導出:加法定理を使って $n=2$ を確かめ,数学的帰納法で証明する
$$ A^{2}=\begin{bmatrix}\cos\theta&-\sin\theta\\\sin\theta&\cos\theta\end{bmatrix}\begin{bmatrix}\cos\theta&-\sin\theta\\\sin\theta&\cos\theta\end{bmatrix}=\begin{bmatrix}\cos^{2}\theta-\sin^{2}\theta&-2\sin\theta\cos\theta\\2\sin\theta\cos\theta&\cos^{2}\theta-\sin^{2}\theta\end{bmatrix} $$ここで,高校数学で学んだ加法定理(高校数学 第5章:三角関数の加法定理と諸公式)から
$$ \cos(\theta+\theta)=\cos^{2}\theta-\sin^{2}\theta,\qquad \sin(\theta+\theta)=2\sin\theta\cos\theta $$が成り立つので,これを代入すると
$$ A^{2}=\begin{bmatrix}\cos2\theta&-\sin2\theta\\\sin2\theta&\cos2\theta\end{bmatrix} $$となり,たしかに $n=2$ のパターンが確認できた.一般の $n$ を数学的帰納法で示す.$n=1$ のとき $A^1=A=\begin{bmatrix}\cos\theta&-\sin\theta\\\sin\theta&\cos\theta\end{bmatrix}=\begin{bmatrix}\cos1\theta&-\sin1\theta\\\sin1\theta&\cos1\theta\end{bmatrix}$ で成り立つ.$n=k$ で $A^{k}=\begin{bmatrix}\cos k\theta&-\sin k\theta\\\sin k\theta&\cos k\theta\end{bmatrix}$ が成り立つと仮定すると,
$$ A^{k+1}=A^{k}A=\begin{bmatrix}\cos k\theta&-\sin k\theta\\\sin k\theta&\cos k\theta\end{bmatrix}\begin{bmatrix}\cos\theta&-\sin\theta\\\sin\theta&\cos\theta\end{bmatrix}=\begin{bmatrix}\cos k\theta\cos\theta-\sin k\theta\sin\theta&-(\cos k\theta\sin\theta+\sin k\theta\cos\theta)\\\sin k\theta\cos\theta+\cos k\theta\sin\theta&-\sin k\theta\sin\theta+\cos k\theta\cos\theta\end{bmatrix} $$ここでも加法定理 $\cos k\theta\cos\theta-\sin k\theta\sin\theta=\cos(k+1)\theta$,$\sin k\theta\cos\theta+\cos k\theta\sin\theta=\sin(k+1)\theta$ を使うと,
$$ A^{k+1}=\begin{bmatrix}\cos(k+1)\theta&-\sin(k+1)\theta\\\sin(k+1)\theta&\cos(k+1)\theta\end{bmatrix} $$となり,$n=k+1$ でも成り立つ.よってすべての自然数 $n$ について式\eqref{eq:8-rotation-power}が成り立つ.
(証明終わり)
例題8.2 3種類のべき乗パターンの数値計算
次の行列の指定されたべき乗を求めよ.(a) $A=\begin{bmatrix}2&0\\0&-1\end{bmatrix}$ の $A^{5}$.(b) $B=\begin{bmatrix}1&3\\0&1\end{bmatrix}$ の $B^{4}$.(c) $C=\begin{bmatrix}\cos\frac{\pi}{6}&-\sin\frac{\pi}{6}\\\sin\frac{\pi}{6}&\cos\frac{\pi}{6}\end{bmatrix}$ の $C^{4}$.
解答 (a) 公式8.1($\alpha=2,\beta=-1,n=5$)より $A^{5}=\begin{bmatrix}2^{5}&0\\0&(-1)^{5}\end{bmatrix}=\begin{bmatrix}32&0\\0&-1\end{bmatrix}$.
(b) 公式8.2($\alpha=3,n=4$)より $B^{4}=\begin{bmatrix}1&4\times3\\0&1\end{bmatrix}=\begin{bmatrix}1&12\\0&1\end{bmatrix}$.
(c) 公式8.3($\theta=\frac{\pi}{6},n=4$)より $C^{4}$ は角 $4\times\frac{\pi}{6}=\frac{2\pi}{3}$($120°$)の回転行列である.$\cos\frac{2\pi}{3}=-\frac{1}{2}$,$\sin\frac{2\pi}{3}=\frac{\sqrt{3}}{2}$(高校数学の三角比の値)だから,
$$ C^{4}=\begin{bmatrix}-\dfrac{1}{2}&-\dfrac{\sqrt{3}}{2}\\[4pt]\dfrac{\sqrt{3}}{2}&-\dfrac{1}{2}\end{bmatrix} $$である.(3問ともsympyで検算済み.)
イメージ:なぜこの3つの行列だけ簡単にべき乗できるのか
対角行列・上三角行列・回転行列のべき乗がこれほどきれいな形になるのは偶然ではない.じつは,多くの正方行列 $A$ は,適切な行列 $P$ を使って $P^{-1}AP$ が対角行列になるように「作り替える」ことができる(これを対角化,diagonalizationという).対角化できてしまえば,$A^{n}=P(P^{-1}AP)^{n}P^{-1}$ という関係を使って,対角行列のべき乗(公式8.1)に帰着させられるので,どんな複雑な行列のべき乗も簡単に求まる.回転行列は複素数の目を通して見ると対角化に近い扱いができ,上三角行列 $\begin{bmatrix}1&\alpha\\0&1\end{bmatrix}$ は対角化できない特殊な例(対角成分が重なっているため)として登場する.この「対角化」の話は第12章(固有値・固有ベクトルと行列の対角化)で本格的に扱う.そこでは,本節で見た3つのパターンが,じつはすべて同じ1つの原理(固有値と固有ベクトル)から導かれることが分かるはずである.
8.2 置換とその符号
次の8.3節で行列式を定義するには,「並べ替え」という操作を数学的に扱う準備が要る.たとえば $1,2,3$ という3つの数を $2,3,1$ の順に並べ替える,といった操作を置換(permutation)という.本節では,置換の定義,逆向きの置換(逆置換),そして置換に「符号」(プラスかマイナスか)を割り当てる方法(偶置換・奇置換)を順に見ていく.
8.2.1 置換の定義
定義8.5 置換(permutation)
$1$ から $n$ までの自然数を,同じ $1$ から $n$ までの自然数のいずれかに1対1に対応させる変換を,$n$ 次の置換という.置換は次のように上下2段の記法で表す.
\begin{equation} \begin{pmatrix}1&2&3&4&\cdots&n\\i_{1}&i_{2}&i_{3}&i_{4}&\cdots&i_{n}\end{pmatrix} \label{eq:8-perm-def} \end{equation}上段(1行目)には $1$ から $n$ までを小さい順に並べ,下段(2行目)には,上段の数 $k$ に対応する行き先 $i_{k}$ を書く.すなわち,この記法は「$1$ は $i_1$ に,$2$ は $i_2$ に,……,$n$ は $i_n$ に対応する」という意味である($i_1,\ldots,i_n$ も $1,\ldots,n$ を並べ替えたものであり,同じ数が2回現れることはない).
例8.1 $n=2$ の置換の例
$n=2$ で $\begin{pmatrix}1\\2\end{pmatrix}\begin{pmatrix}2\\1\end{pmatrix}$ と略記される置換 $\begin{pmatrix}1&2\\2&1\end{pmatrix}$ は,「$1$ は $2$ に,$2$ は $1$ に対応させる」という置換である.
もう少し要素数の多い $n=4$ の置換で,上下2段の対応を矢印で表した図を見てみよう.図8.2に,置換 $\begin{pmatrix}1&2&3&4\\3&1&2&4\end{pmatrix}$(「$1$は$3$に,$2$は$1$に,$3$は$2$に,$4$は$4$に対応」)の対応を矢印で表した図を示す.
$n$ の値によって,置換が何通りあるかを数えてみよう.
$n=2$ のとき,$\begin{pmatrix}1&2\\1&2\end{pmatrix}$,$\begin{pmatrix}1&2\\2&1\end{pmatrix}$ の $2$ 通り($2!=2\times1=2$ 通り)の置換がある.
$n=3$ のとき,$\begin{pmatrix}1&2&3\\1&2&3\end{pmatrix}$,$\begin{pmatrix}1&2&3\\1&3&2\end{pmatrix}$,$\begin{pmatrix}1&2&3\\2&1&3\end{pmatrix}$,$\begin{pmatrix}1&2&3\\2&3&1\end{pmatrix}$,$\begin{pmatrix}1&2&3\\3&1&2\end{pmatrix}$,$\begin{pmatrix}1&2&3\\3&2&1\end{pmatrix}$ の $6$ 通り($3!=3\times2\times1=6$ 通り)の置換がある.
一般に,$1$ の行き先の選び方が $n$ 通り,残りの $n-1$ 個から $2$ の行き先を選ぶ方法が $n-1$ 通り,……と数えていくと,$n$ 次の置換はちょうど $n!$($n$ の階乗,$n$ factorial)通りあることが分かる.これは高校数学の順列の考え方(異なる $n$ 個のものを1列に並べる並べ方の総数)とまったく同じである.
8.2.2 逆置換
定義8.6 逆置換(inverse permutation)
ある置換 $\sigma$ に対して,$\sigma$ とちょうど逆の対応を行う置換を,$\sigma$ の逆置換という.すなわち,$\sigma$ が「$k$ を $i_k$ に対応させる」ものなら,$\sigma$ の逆置換は「$i_k$ を $k$ に対応させる」ものである.
例として,$\sigma=\begin{pmatrix}1&2&3&4\\3&1&2&4\end{pmatrix}$($1$ は $3$ に,$2$ は $1$ に,$3$ は $2$ に,$4$ は $4$ に対応)の逆置換を考える.「A は B に対応している」というのを「B は A に対応している」に変えると,$3$ は $1$ に,$1$ は $2$ に,$2$ は $3$ に,$4$ は $4$ に対応することになる.よって $\sigma$ の逆置換は $\begin{pmatrix}1&2&3&4\\2&3&1&4\end{pmatrix}$ となる.
逆置換を求める手順は,次の2段階に整理できる.
(i) 上段と下段を Exchange(交換)する:$\begin{pmatrix}1&2&3&4\\3&1&2&4\end{pmatrix}\ \longrightarrow\ \begin{pmatrix}3&1&2&4\\1&2&3&4\end{pmatrix}$
(ii) 上段の自然数を小さい順に並べ替える(下段も一緒に動かす):$\begin{pmatrix}3&1&2&4\\1&2&3&4\end{pmatrix}\ \longrightarrow\ \begin{pmatrix}1&2&3&4\\2&3&1&4\end{pmatrix}$
図8.3に,この2段階の手順を示す.
8.2.3 偶置換と奇置換
置換のうち2つの数の対応だけを入れ換える操作を互換(transposition)という.たとえば $\begin{pmatrix}1&2&3&4\\1&2&3&4\end{pmatrix}$(恒等置換,identity permutation:すべての数がそれ自身に対応する,何も変えない置換)から出発して,下段の2つの数字を入れ換える操作を何回か繰り返すと,任意の置換を作ることができる.
定義8.7 偶置換と奇置換
ある置換が,恒等置換から偶数回の互換で作れるとき,その置換を偶置換(even permutation)という.奇数回の互換で作れるとき,奇置換(odd permutation)という.(同じ置換でも互換の仕方は何通りもあるが,その回数の偶奇はどの作り方でも変わらないことが知られている.)
例として,$\begin{pmatrix}1&2&3&4\\3&1&2&4\end{pmatrix}$ が偶置換か奇置換かを調べよう.恒等置換 $\begin{pmatrix}1&2&3&4\\1&2&3&4\end{pmatrix}$ から出発し,下段の数字を互換で少しずつ目標の並びに近づけていく.
まず下段の $3$ と $1$ を互換すると $\begin{pmatrix}1&2&3&4\\1&2&3&4\end{pmatrix}\to\begin{pmatrix}1&2&3&4\\3&2&1&4\end{pmatrix}$.次に下段の $2$ と $1$ を互換すると $\begin{pmatrix}1&2&3&4\\3&2&1&4\end{pmatrix}\to\begin{pmatrix}1&2&3&4\\3&1&2&4\end{pmatrix}$.これで目標の置換に到達し,互換の回数は $2$ 回(偶数回)である.よって $\begin{pmatrix}1&2&3&4\\3&1&2&4\end{pmatrix}$ は偶置換である.
例題8.3 5次の置換の逆置換と偶奇の判定
置換 $\sigma=\begin{pmatrix}1&2&3&4&5\\4&1&3&5&2\end{pmatrix}$ について,(a) 逆置換 $\sigma^{-1}$ を求めよ.(b) $\sigma$ が偶置換か奇置換かを判定せよ.
解答 (a) 定義8.6の手順に従う.(i) 上下を交換:$\begin{pmatrix}4&1&3&5&2\\1&2&3&4&5\end{pmatrix}$.(ii) 上段を小さい順($1,2,3,4,5$)に並べ替え,下段も一緒に動かす.上段で $1$ は元の第2列にあったので下段も第2列の $2$ が付いてくる.同様に $2$ は元の第5列にあったので下段は $5$,$3$ は第3列のままなので下段は $3$,$4$ は元の第1列にあったので下段は $1$,$5$ は元の第4列にあったので下段は $4$.よって
$$ \sigma^{-1}=\begin{pmatrix}1&2&3&4&5\\2&5&3&1&4\end{pmatrix} $$(検算:$\sigma$ で $1\to4$ なら $\sigma^{-1}$ で $4\to1$ になっているはずである.実際 $\sigma^{-1}$ の表で上段 $4$ の下は $1$ であり,一致している.同様に $\sigma$ の $2\to1,3\to3,4\to5,5\to2$ に対し,$\sigma^{-1}$ では $1\to2,3\to3,5\to4,2\to5$ となっていて,すべて矛盾なく対応している.)
(b) $\sigma$ の下段の並び $(4,1,3,5,2)$ を恒等置換 $(1,2,3,4,5)$ から互換で作る代わりに,並びの中で「大きい数が小さい数より先(左)に来ている組」(転倒数,inversion number という.転倒数の偶奇と互換回数の偶奇は一致することが知られている)を数える方法を使う.$(4,1,3,5,2)$ の中で,位置の左右を問わず値の大小が逆転している組をすべて挙げると,$(4,1)$,$(4,3)$,$(4,2)$,$(3,2)$,$(5,2)$ の $5$ 組(たとえば $4$ は $1,3,2$ より前にあるのにこれらより大きいので3組,$5$ は $2$ より前にあるのに大きいので1組,さらに $3$ も $2$ より前にあるのに大きいので1組,合計 $5$ 組)である.転倒数は $5$(奇数)だから,$\sigma$ は奇置換である.(sympyで転倒数と符号を検算し,どちらも一致することを確認済み.)
8.3 行列式の定義
8.3.1 符号 $\sgn$ とΣ記号
置換に,プラス($+1$)かマイナス($-1$)かの「符号」を割り当てる関数を定義する.
定義8.8 置換の符号 $\sgn$(signature)
置換 $\sigma=\begin{pmatrix}1&2&\cdots&n\\i_{1}&i_{2}&\cdots&i_{n}\end{pmatrix}$ に対して,
$$ \sgn\sigma=\begin{cases}+1&(\sigma\text{ が偶置換のとき})\\-1&(\sigma\text{ が奇置換のとき})\end{cases} $$と定める.$\sgn$ は signature(符号)の略で,「符号」と呼ぶ.8.2.3項の例では $\sgn\begin{pmatrix}1&2&3&4\\3&1&2&4\end{pmatrix}=+1$ である.
また,本節以降で頻繁に使う $\Sigma$(シグマ,Sigma.ギリシャ文字で,和の記号)についても確認しておく.$\displaystyle\sum_{k=1}^{n}c_k$ は「$c_1$ から $c_n$ まで,$k$ を $1$ から $n$ まで動かしながら全部足す」という意味である.次の定義8.9では,添字を省略して単に $\sum_{\sigma}$(「すべての置換 $\sigma$ にわたって足す」の意味)と書く.
8.3.2 行列式の定義
定義8.9 行列式(determinant)
$n$ 次正方行列 $A=\begin{bmatrix}a_{11}&a_{12}&\cdots&a_{1n}\\a_{21}&a_{22}&\cdots&a_{2n}\\\vdots&\vdots&&\vdots\\a_{n1}&a_{n2}&\cdots&a_{nn}\end{bmatrix}$ の行列式(determinant)$|A|$($\det A$ とも書く)を,
\begin{equation} |A|=\sum_{\sigma}\sgn\begin{pmatrix}1&2&3&\cdots&n\\i_{1}&i_{2}&i_{3}&\cdots&i_{n}\end{pmatrix}\,a_{1i_{1}}a_{2i_{2}}a_{3i_{3}}\cdots a_{ni_{n}} \label{eq:8-det-def} \end{equation}で定める.ここで和 $\sum_\sigma$ は,$n$ 次のすべての置換 $\sigma=\begin{pmatrix}1&2&\cdots&n\\i_1&i_2&\cdots&i_n\end{pmatrix}$(全部で $n!$ 通り)にわたってとる.すなわち行列式とは,「各行から1個ずつ,かつ各列から1個ずつ成分を選んで作った積」を,選び方(=置換)ごとに符号 $\sgn\sigma$ を付けて全部足し合わせたものである.
この定義は,17世紀のドイツの数学者ライプニッツ(Gottfried Wilhelm Leibniz,1646–1716)に由来することからライプニッツの公式(Leibniz formula)とも呼ばれる(実際にはライプニッツ自身がこの形で発表したわけではないが,行列式の考え方を最初に用いた一人とされる).「各行・各列から1個ずつ」という条件は,将棋やチェスのルーク(飛車)を $n\times n$ の盤に,互いの利き筋(同じ行・同じ列)が重ならないように $n$ 個置く置き方に対応する.図8.5に,$n=3$ の場合の $3!=6$ 通りすべての置き方(=置換)と,それぞれの符号・対応する項を示す.
この定義を使って,一般の $3$ 次正方行列 $A=\begin{bmatrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\\a_{31}&a_{32}&a_{33}\end{bmatrix}$ の行列式 $|A|$ を実際に書き下してみよう.
導出:$3$次の行列式(サラスの公式)を定義から書き下す
$n=3$ の置換は8.2.1項で数えたとおり $3!=6$ 通りある.それぞれの置換について $\sgn$ を求めると(8.2.3項の方法で1つずつ確かめられる),
$$ |A|=\sum\sgn\begin{pmatrix}1&2&3\\i_1&i_2&i_3\end{pmatrix}a_{1i_1}a_{2i_2}a_{3i_3} $$ \begin{equation} \begin{aligned} &=\underbrace{\sgn\begin{pmatrix}1&2&3\\1&2&3\end{pmatrix}}_{+1}a_{11}a_{22}a_{33}+\underbrace{\sgn\begin{pmatrix}1&2&3\\1&3&2\end{pmatrix}}_{-1}a_{11}a_{23}a_{32}\\ &\quad+\underbrace{\sgn\begin{pmatrix}1&2&3\\2&1&3\end{pmatrix}}_{-1}a_{12}a_{21}a_{33}+\underbrace{\sgn\begin{pmatrix}1&2&3\\2&3&1\end{pmatrix}}_{+1}a_{12}a_{23}a_{31}\\ &\quad+\underbrace{\sgn\begin{pmatrix}1&2&3\\3&1&2\end{pmatrix}}_{+1}a_{13}a_{21}a_{32}+\underbrace{\sgn\begin{pmatrix}1&2&3\\3&2&1\end{pmatrix}}_{-1}a_{13}a_{22}a_{31} \end{aligned} \label{eq:8-det3-leibniz} \end{equation}(各項の符号は,8.2.3項で $\begin{pmatrix}1&2&3\\3&1&2\end{pmatrix}$ について行ったのと同じように,恒等置換から何回互換すれば到達できるかを数えれば求まる.たとえば $\begin{pmatrix}1&2&3\\2&3&1\end{pmatrix}$ は,$1\leftrightarrow2$ の互換で $\begin{pmatrix}1&2&3\\2&1&3\end{pmatrix}$ にし,続けて $1\leftrightarrow3$(下段の1と3)の互換で $\begin{pmatrix}1&2&3\\2&3&1\end{pmatrix}$ に到達するので,2回の互換=偶置換で $\sgn=+1$ となる.)これらの符号をまとめて代入すると,
\begin{equation} |A| = a_{11}a_{22}a_{33} - a_{11}a_{23}a_{32} - a_{12}a_{21}a_{33} + a_{12}a_{23}a_{31} + a_{13}a_{21}a_{32} - a_{13}a_{22}a_{31} \label{eq:8-det3-sarrus} \end{equation}を得る.これが,高校数学や物理でもよく使われる「サラスの公式」(Sarrus's rule:$3$次行列式を,右に列を2つ複製して斜めの積の和と差で計算する覚え方)の中身である.本章の定義8.9は,このサラスの公式を「なぜそう計算してよいのか」の根っこから説明する,より一般的な枠組みになっている.
例題8.4 $3$次行列式の計算
$A=\begin{bmatrix}2&-1&3\\0&4&1\\5&2&-2\end{bmatrix}$ の行列式 $|A|$ を,式\eqref{eq:8-det3-sarrus}を使って求めよ.
解答 $a_{11}=2,a_{12}=-1,a_{13}=3,a_{21}=0,a_{22}=4,a_{23}=1,a_{31}=5,a_{32}=2,a_{33}=-2$ を式\eqref{eq:8-det3-sarrus}に代入する.
$$ \begin{aligned} |A| &= a_{11}a_{22}a_{33} - a_{11}a_{23}a_{32} - a_{12}a_{21}a_{33} + a_{12}a_{23}a_{31} + a_{13}a_{21}a_{32} - a_{13}a_{22}a_{31}\\ &= 2\times4\times(-2) - 2\times1\times2 - (-1)\times0\times(-2) + (-1)\times1\times5 + 3\times0\times2 - 3\times4\times5\\ &= -16 - 4 - 0 - 5 + 0 - 60\\ &= -85 \end{aligned} $$よって $|A|=-85$ である.(sympyで直接 $\det$ を計算し,同じ $-85$ になることを確認済み.)
8.4 余因子と行列式の展開
定義8.9のライプニッツの公式は,行列式が何であるかを明確に定めてくれるが,$n$ が大きくなると項の数が $n!$ 個に爆発的に増えるため,実際の計算には向かない.そこで,$n$ 次の行列式を,より小さな $(n-1)$ 次の行列式の組み合わせに帰着させる方法——余因子展開(cofactor expansion)——を用意する.
定義8.10 余因子(cofactor)
$n$ 次正方行列 $A=(a_{ij})$ に対して,$A$ から第 $i$ 行と第 $j$ 列を取り除いてできる $(n-1)$ 次の行列の行列式に $(-1)^{i+j}$ を掛けたものを,$A$ の $(i,j)$ 余因子(cofactor)といい,$A_{ij}$ と書く:
\begin{equation} A_{ij}=(-1)^{i+j}\begin{vmatrix}a_{11}&\cdots&\widehat{a_{1j}}&\cdots&a_{1n}\\ \vdots&&\vdots&&\vdots\\ \widehat{a_{i1}}&\cdots&\widehat{a_{ij}}&\cdots&\widehat{a_{in}}\\ \vdots&&\vdots&&\vdots\\ a_{n1}&\cdots&\widehat{a_{nj}}&\cdots&a_{nn}\end{vmatrix} \label{eq:8-cofactor-def} \end{equation}ここで $\widehat{\phantom{a}}$ を付けた行・列は「取り除く」ことを表す記号であり,実際の行列式には第 $i$ 行と第 $j$ 列を除いた残りの $(n-1)\times(n-1)$ 個の成分だけが並ぶ.図8.6に,この「第 $i$ 行・第 $j$ 列を取り除く」操作を図示する.
この定義が正しいことを確かめるために,$A_1=\begin{bmatrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\\a_{31}&a_{32}&a_{33}\end{bmatrix}$ の行列式を求めてみよう.
証明:$3$次の行列式を余因子で書き直す
式\eqref{eq:8-det3-sarrus}(サラスの公式)より,
$$ |A_1| = a_{11}a_{22}a_{33} + a_{12}a_{23}a_{31} + a_{13}a_{21}a_{32} - a_{13}a_{22}a_{31} - a_{12}a_{21}a_{33} - a_{11}a_{23}a_{32} $$(項の順序を並べ替えただけで,式\eqref{eq:8-det3-sarrus}と同じ6項である.)これを $a_{11},a_{12},a_{13}$ でくくり出すと,
$$ |A_1| = a_{11}(a_{22}a_{33} - a_{23}a_{32}) - a_{12}(a_{23}a_{31} - a_{21}a_{33}) + a_{13}(a_{21}a_{32} - a_{22}a_{31}) $$($a_{12}$ の項の符号がマイナスであることに注意.括弧の中を $a_{23}a_{31}-a_{21}a_{33}$ の順に書いたのは,直後の $2$ 次行列式の並びに合わせるためである.)ここで,右辺の3つの括弧の中身は,$A_1$ からそれぞれ第1行と第1列,第2列,第3列を取り除いた $2$ 次の行列式に一致している:
$$ \begin{vmatrix}a_{22}&a_{23}\\a_{32}&a_{33}\end{vmatrix}=a_{22}a_{33}-a_{23}a_{32},\qquad \begin{vmatrix}a_{21}&a_{23}\\a_{31}&a_{33}\end{vmatrix}=a_{21}a_{33}-a_{23}a_{31},\qquad \begin{vmatrix}a_{21}&a_{22}\\a_{31}&a_{32}\end{vmatrix}=a_{21}a_{32}-a_{22}a_{31} $$定義8.10に従えば,これらはそれぞれ $A_{11}=(-1)^{1+1}\begin{vmatrix}a_{22}&a_{23}\\a_{32}&a_{33}\end{vmatrix}$,$A_{12}=(-1)^{1+2}\begin{vmatrix}a_{21}&a_{23}\\a_{31}&a_{33}\end{vmatrix}$,$A_{13}=(-1)^{1+3}\begin{vmatrix}a_{21}&a_{22}\\a_{31}&a_{32}\end{vmatrix}$ である.$(-1)^{1+1}=1$,$(-1)^{1+2}=-1$,$(-1)^{1+3}=1$ なので,$A_{11}=a_{22}a_{33}-a_{23}a_{32}$,$A_{12}=-(a_{21}a_{33}-a_{23}a_{31})=a_{23}a_{31}-a_{21}a_{33}$,$A_{13}=a_{21}a_{32}-a_{22}a_{31}$ であり,これはまさに上の3つの括弧の中身と一致する.よって
\begin{equation} |A_1|=a_{11}A_{11}+a_{12}A_{12}+a_{13}A_{13} \label{eq:8-cofactor-row1} \end{equation}が成り立つ.すなわち,$|A_1|$ は第1行の成分 $a_{11},a_{12},a_{13}$ と,それぞれに対応する余因子 $A_{11},A_{12},A_{13}$ の積の和として書ける.
(証明終わり)
式\eqref{eq:8-cofactor-row1}は「第1行に関する余因子展開」だが,同じことは任意の行・任意の列について成り立つ.これを一般の $n$ 次の場合にまとめたのが,次の定理である.
定理8.2 余因子展開(cofactor expansion)
$n$ 次正方行列 $A=(a_{ij})$ の行列式は,任意の行 $i$ に関して
\begin{equation} |A| = a_{i1}A_{i1}+a_{i2}A_{i2}+a_{i3}A_{i3}+\cdots+a_{in}A_{in} \label{eq:8-cofactor-expansion-row} \end{equation}と展開できる(第 $i$ 行に関する展開).また,任意の列 $j$ に関しても同様に
\begin{equation} |A| = a_{1j}A_{1j}+a_{2j}A_{2j}+a_{3j}A_{3j}+\cdots+a_{nj}A_{nj} \label{eq:8-cofactor-expansion-col} \end{equation}と展開できる(第 $j$ 列に関する展開).
なぜ一般の $n$ でも成り立つのか(証明の考え方)
上の証明で見たように,$n=3$ の場合は,サラスの公式の6項を「第1行のどの成分を含むか」で3つのグループ($a_{11}$ を含む2項,$a_{12}$ を含む2項,$a_{13}$ を含む2項)に仕分けし,各グループが対応する余因子に一致することを確かめた.一般の $n$ でも考え方は同じで,定義8.9のライプニッツの公式($n!$ 個の項の和)を「置換 $\sigma$ が第 $i$ 行の成分をどの列に対応させるか($\sigma(i)=j$ を固定する)」で $n$ 個のグループに仕分けると,各グループの和がちょうど $a_{ij}A_{ij}$ に一致することが示せる.この仕分けの計算は,添字の付け替えが煩雑になるため本書では省略するが,$n=3$ の場合で確認した「くくり出し」の操作を,より多くの項について辛抱強く繰り返すだけである,という点をつかんでおけば十分である.
例題8.5 余因子展開による$4$次行列式の計算
$M=\begin{bmatrix}1&2&0&-1\\0&1&3&2\\0&0&1&4\\2&-1&0&1\end{bmatrix}$ の行列式 $|M|$ を,余因子展開を使って求めよ.
解答 第3行 $(0,0,1,4)$ に $0$ が2個あるので,第3行に関して展開する(式\eqref{eq:8-cofactor-expansion-row}で $i=3$).$0$ の係数の項は計算しなくてよいので,
$$ |M| = 0\cdot A_{31}+0\cdot A_{32}+1\cdot A_{33}+4\cdot A_{34} $$だけを計算すればよい.$A_{33}$ は,$M$ から第3行・第3列を取り除いた$3$次行列式に $(-1)^{3+3}=1$ を掛けたもの:
$$ A_{33}=(-1)^{3+3}\begin{vmatrix}1&2&-1\\0&1&2\\2&-1&1\end{vmatrix}=\begin{vmatrix}1&2&-1\\0&1&2\\2&-1&1\end{vmatrix} $$サラスの公式(式\eqref{eq:8-det3-sarrus})で計算すると,$1\times1\times1+2\times2\times2+(-1)\times0\times(-1)-(-1)\times1\times2-2\times0\times1-1\times2\times(-1)=1+8+0+2-0+2=13$.よって $A_{33}=13$.
$A_{34}$ は,$M$ から第3行・第4列を取り除いた$3$次行列式に $(-1)^{3+4}=-1$ を掛けたもの:
$$ A_{34}=(-1)^{3+4}\begin{vmatrix}1&2&0\\0&1&3\\2&-1&0\end{vmatrix}=-\begin{vmatrix}1&2&0\\0&1&3\\2&-1&0\end{vmatrix} $$サラスの公式で $1\times1\times0+2\times3\times2+0\times0\times(-1)-0\times1\times2-2\times0\times0-1\times3\times(-1)=0+12+0-0-0+3=15$.よって $A_{34}=-15$.
以上より,
$$ |M| = 1\times13+4\times(-15)=13-60=-47 $$である.(sympyで直接 $\det$ を計算し,同じ $-47$ になることを確認済み.)
8.5 行列式の性質
行列式には,計算を大幅に楽にしてくれる基本的な性質が4つある.ここでは $2$ 次の行列式で確かめるが,ノートの欄外にある注記のとおり,これらの性質は列についても全く同じように成り立ち,また$2$次に限らず$n$次の行列式でも成り立つ(証明の考え方は8.4節の余因子展開と同様,$n$次の場合に一般化できる).
定理8.3 行列式の基本性質
(i) 1つの行を $c$ 倍すると,行列式も $c$ 倍になる.
\begin{equation} \begin{vmatrix}ca_{11}&ca_{12}\\da_{21}&da_{22}\end{vmatrix}=cd\begin{vmatrix}a_{11}&a_{12}\\a_{21}&a_{22}\end{vmatrix} \label{eq:8-prop-scalar} \end{equation}(ii) 2つの行を入れ替えると,行列式の符号が変わる.
\begin{equation} \begin{vmatrix}a_{11}&a_{12}\\a_{21}&a_{22}\end{vmatrix}=-\begin{vmatrix}a_{21}&a_{22}\\a_{11}&a_{12}\end{vmatrix} \label{eq:8-prop-swap} \end{equation}(iii) ある行に,他の行の定数倍を加えても,行列式は変わらない.
\begin{equation} \begin{vmatrix}a_{11}&a_{12}\\a_{21}&a_{22}\end{vmatrix}=\begin{vmatrix}a_{11}&a_{12}\\a_{21}+ca_{11}&a_{22}+ca_{12}\end{vmatrix} \label{eq:8-prop-addrow} \end{equation}($c$ は任意の実数.ノートに書かれている $a_{21}\pm a_{11}$ の形は,$c=\pm1$ の特別な場合にあたる.)
(iv) ある行を2つのベクトルの和に分けると,行列式もその2つの行列式の和に分けられる(多重線形性).
\begin{equation} \begin{vmatrix}a_{11}&a_{12}\\a_{21}+b_{21}&a_{22}+b_{22}\end{vmatrix}=\begin{vmatrix}a_{11}&a_{12}\\a_{21}&a_{22}\end{vmatrix}+\begin{vmatrix}a_{11}&a_{12}\\b_{21}&b_{22}\end{vmatrix} \label{eq:8-prop-additive} \end{equation}証明:(i) 行のスカラー倍
$$ \begin{aligned} \begin{vmatrix}ca_{11}&ca_{12}\\da_{21}&da_{22}\end{vmatrix} &= ca_{11}\cdot da_{22} - ca_{12}\cdot da_{21}\\ &= cd\cdot a_{11}a_{22} - cd\cdot a_{12}a_{21}\qquad(\text{$c,d$ をくくり出した})\\ &= cd(a_{11}a_{22}-a_{12}a_{21})\\ &= cd\begin{vmatrix}a_{11}&a_{12}\\a_{21}&a_{22}\end{vmatrix} \end{aligned} $$(証明終わり)
証明:(ii) 行の入れ替え
$$ \begin{aligned} \begin{vmatrix}a_{11}&a_{12}\\a_{21}&a_{22}\end{vmatrix} &= a_{11}a_{22} - a_{12}a_{21}\\ &= -(a_{12}a_{21} - a_{11}a_{22})\qquad(\text{$-1$ をくくり出した})\\ &= -\begin{vmatrix}a_{21}&a_{22}\\a_{11}&a_{12}\end{vmatrix} \end{aligned} $$(最後の行列式は,1行目に $a_{21},a_{22}$,2行目に $a_{11},a_{12}$ を並べたものであり,その値を定義通り計算すると $a_{21}a_{12}-a_{22}a_{11}=-(a_{11}a_{22}-a_{12}a_{21})$ の符号を反転したもの,すなわち $a_{12}a_{21}-a_{11}a_{22}$ に等しいことを使っている.)
(証明終わり)
証明:(iii) 他の行の定数倍を加えても不変
任意の実数 $c$ で確かめる(ノートにある $c=\pm1$ の場合は,この一般の証明の特別な場合である).
$$ \begin{aligned} \begin{vmatrix}a_{11}&a_{12}\\a_{21}+ca_{11}&a_{22}+ca_{12}\end{vmatrix} &= a_{11}(a_{22}+ca_{12}) - a_{12}(a_{21}+ca_{11})\\ &= a_{11}a_{22} + \underline{ca_{11}a_{12}} - \underline{ca_{12}a_{11}} - a_{12}a_{21}\qquad(\text{下線部は同じ量なので,$c$ の値によらず互いに打ち消し合う})\\ &= a_{11}a_{22} - a_{12}a_{21}\\ &= \begin{vmatrix}a_{11}&a_{12}\\a_{21}&a_{22}\end{vmatrix} \end{aligned} $$(証明終わり)
証明:(iv) 行の分割と和(多重線形性)
右辺を先に展開する.
$$ \begin{aligned} (\text{右辺}) &= (a_{11}a_{22} - a_{12}a_{21}) + (a_{11}b_{22} - a_{12}b_{21})\\ &= a_{11}(a_{22}+b_{22}) - a_{12}(a_{21}+b_{21})\qquad(a_{11},a_{12}\text{ をくくり出した})\\ &= \begin{vmatrix}a_{11}&a_{12}\\a_{21}+b_{21}&a_{22}+b_{22}\end{vmatrix}\qquad(\text{定義通り書き直した})\\ &= (\text{左辺}) \end{aligned} $$(証明終わり)
これらの性質を組み合わせると,大きな行列式でも,高校数学の連立方程式の消去法と同じ要領で「上三角形」の形(対角線より左下がすべて $0$)に変形してから,対角成分の積を計算するだけで求められる.実際,上三角行列の行列式は対角成分の積に等しい——これは定義8.9のライプニッツの公式で,下三角部分がすべて $0$ である以上,恒等置換 $\sigma=\mathrm{id}$ 以外のどの置換を選んでも積 $a_{1i_1}\cdots a_{ni_n}$ のどこかに $0$ が現れてしまう(下三角の成分を選ばざるを得なくなる)ことから確かめられる,覚えておくと便利な事実である.
例題8.6 行の性質を使った$4$次行列式の効率的な計算
$P=\begin{bmatrix}2&4&2&0\\1&3&1&2\\3&5&4&1\\1&1&2&3\end{bmatrix}$ の行列式 $|P|$ を,定理8.3の性質を使って求めよ.
解答 まず第1行の成分がすべて偶数なので,性質(i)を使って $2$ をくくり出す.$R_1$で第1行を表すと,$R_1=(2,4,2,0)=2\times(1,2,1,0)$ だから,
$$ |P| = 2\begin{vmatrix}1&2&1&0\\1&3&1&2\\3&5&4&1\\1&1&2&3\end{vmatrix} $$次に,性質(iii)を繰り返し使って,第1列の第1行より下をすべて $0$ にする:$R_2\leftarrow R_2-R_1$,$R_3\leftarrow R_3-3R_1$,$R_4\leftarrow R_4-R_1$(いずれも「他の行の定数倍を加える」操作なので行列式は変わらない).
$$ 2\begin{vmatrix}1&2&1&0\\1&3&1&2\\3&5&4&1\\1&1&2&3\end{vmatrix} = 2\begin{vmatrix}1&2&1&0\\0&1&0&2\\0&-1&1&1\\0&-1&1&3\end{vmatrix} $$続けて,第2列の第2行より下を $0$ にする:$R_3\leftarrow R_3+R_2$,$R_4\leftarrow R_4+R_2$(第2行の $(2,2)$ 成分は $1$ なので,そのまま足せばよい).
$$ 2\begin{vmatrix}1&2&1&0\\0&1&0&2\\0&-1&1&1\\0&-1&1&3\end{vmatrix} = 2\begin{vmatrix}1&2&1&0\\0&1&0&2\\0&0&1&3\\0&0&1&5\end{vmatrix} $$最後に,第3列の第3行より下を $0$ にする:$R_4\leftarrow R_4-R_3$.
$$ 2\begin{vmatrix}1&2&1&0\\0&1&0&2\\0&0&1&3\\0&0&1&5\end{vmatrix} = 2\begin{vmatrix}1&2&1&0\\0&1&0&2\\0&0&1&3\\0&0&0&2\end{vmatrix} $$ここまでで使った操作はすべて性質(i)(の逆=くくり出し)と性質(iii)(他の行の定数倍を加える)だけであり,行の入れ替え(性質(ii))は使っていないので符号は変わらない.最後の行列式は上三角行列だから,対角成分の積 $1\times1\times1\times2=2$ に等しい.よって
$$ |P| = 2\times2 = 4 $$である.(sympyで直接 $\det$ を計算し,同じ $4$ になることを確認済み.)
8.6 行列式の幾何学的意味
ここまで行列式を,置換の符号を使った代数的な計算規則として扱ってきた.しかし行列式には,もう1つの顔——図形の大きさを表すという顔——がある.行列式が $0$ であることと $0$ でないことは,いったい何を意味するのだろうか.具体的な数値とグラフで確かめてみよう.
8.6.1 2次元:面積
ベクトル $\vec{a}=\begin{bmatrix}2\\1\end{bmatrix}$,$\vec{b}=\begin{bmatrix}2\\4\end{bmatrix}$ を考える.$\vec{a}$ と $\vec{b}$ を列として並べると,行列 $A=\begin{bmatrix}2&2\\1&4\end{bmatrix}$ になる.
行列式を計算すると,
\begin{equation} |A| = 2\times4-2\times1 = 8-2 = 6 \qquad\cdots① \label{eq:8-area-detA} \end{equation}他方,$\vec{a},\vec{b}$ が作る平行四辺形(図8.7左,原点 $O$,$\vec{a}$ の先端 $(2,1)$,$\vec{a}+\vec{b}$ の先端 $(4,5)$,$\vec{b}$ の先端 $(2,4)$ を頂点とする)の面積 $S_1$ を,これらの点を囲む長方形から余分な部分を引く方法で求めてみよう.4点をちょうど対角の2頂点に持つ長方形は,$x$ 座標の範囲 $[0,4]$,$y$ 座標の範囲 $[0,5]$ なので,面積は $4\times5=20$ である.この長方形から平行四辺形の部分を除いた「余り」は,$\vec a=(a_1,a_2)=(2,1)$ が作る小さな直角三角形2つ(合わせて長方形 $a_1\times a_2$ 分)と,$\vec b=(b_1,b_2)=(2,4)$ が作る直角三角形2つ(合わせて長方形 $b_1\times b_2$ 分),そして幅 $b_1$・高さ $a_2$ の長方形2つ分に分かれる(一般に $a_1,b_1>0$ かつ $0\lt a_2\lt b_2$ のとき成り立つ分解であり,代数的には $(a_1+b_1)(a_2+b_2)=(a_1a_2+b_1b_2)+2a_2b_1+(a_1b_2-a_2b_1)$ という恒等式に対応する):
\begin{equation} S_1 = 20-(4+10)=6 \qquad\cdots② \label{eq:8-area-S1} \end{equation}ここで $10=a_1a_2+b_1b_2=2\times1+2\times4=2+8$(2組の直角三角形の分),$4=2a_2b_1=2\times1\times2$(2つの長方形の分)である.①=②より,$|A|$ は平行四辺形の面積 $S_1$ を示していたことが分かる.
もう一組,$\vec{c}=\begin{bmatrix}2\\1\end{bmatrix}$,$\vec{d}=\begin{bmatrix}4\\2\end{bmatrix}$ を考える.$\vec c$ と $\vec d$ を合わせると $B=\begin{bmatrix}2&4\\1&2\end{bmatrix}$ となる.行列式は
$$ |B| = 4-4 = 0 $$である.図8.7右に示すように,$\vec d=2\vec c$(同じ向きで長さだけ2倍)なので,$\vec c,\vec d$ は同一直線上にあり,平行四辺形ができない(つぶれて線分になる).よって面積は $0$ である.
このことから,2次元グラフでは,行列式が $0$ でなければ面積の値を示し,$0$ であれば2つのベクトルが重なり,面積が$0$である線が表れる,とまとめられる.これを一般の $2\times2$ 行列で証明しておこう.
定理8.4 行列式の幾何学的意味(2次元)
ベクトル $\vec a=(a_1,a_2)$,$\vec b=(b_1,b_2)$ を列として並べた行列を $A=\begin{bmatrix}a_1&b_1\\a_2&b_2\end{bmatrix}$ とする.このとき,$\vec a,\vec b$ が作る平行四辺形の面積は $|A|=|a_1b_2-a_2b_1|$(絶対値)に等しい.とくに $|A|=0$ であることと,$\vec a,\vec b$ が同一直線上にある(一方が他方の実数倍である)こととは同値である.
証明
$\vec a,\vec b$ を,原点 $O$ から測った角を $\alpha,\beta$($x$ 軸から測った角),長さを $r=\|\vec a\|,\ s=\|\vec b\|$ として $\vec a=(r\cos\alpha,r\sin\alpha)$,$\vec b=(s\cos\beta,s\sin\beta)$ と書く.平行四辺形の面積は「底辺 $\times$ 高さ」で求められ,底辺を $\vec a$(長さ $r$)に取ると,高さは $\vec b$ から $\vec a$ の方向への垂線の長さ,すなわち $\vec b$ と $\vec a$ のなす角を $\theta=\beta-\alpha$ として $s|\sin\theta|$ である.よって面積は
$$ S = rs|\sin\theta| = rs|\sin(\beta-\alpha)| = rs|\sin\beta\cos\alpha-\cos\beta\sin\alpha| $$(高校数学の加法定理 $\sin(\beta-\alpha)=\sin\beta\cos\alpha-\cos\beta\sin\alpha$ を使った).ここで $r\cos\alpha=a_1$,$r\sin\alpha=a_2$,$s\cos\beta=b_1$,$s\sin\beta=b_2$ を代入すると,
$$ rs(\sin\beta\cos\alpha-\cos\beta\sin\alpha) = (s\sin\beta)(r\cos\alpha)-(s\cos\beta)(r\sin\alpha) = b_2a_1-b_1a_2 = a_1b_2-a_2b_1 $$となり,$S=|a_1b_2-a_2b_1|=|A|$(絶対値をつけたもの)が示された.また,$\vec a,\vec b$ が同一直線上にある($\vec b=k\vec a$ となる実数 $k$ が存在する,または $\vec a=\vec 0$)ことと,$\theta=0$ または $\pi$(あるいは $\vec a,\vec b$ の一方が $\vec 0$)となることは同値であり,これは $\sin\theta=0$,すなわち $S=0$,すなわち $|A|=0$ と同値である.
(証明終わり)
例題8.7 平行四辺形の面積
$\vec a=(3,1)$,$\vec b=(1,4)$ が作る平行四辺形の面積を求めよ.
解答 定理8.4より,面積は $|a_1b_2-a_2b_1|=|3\times4-1\times1|=|12-1|=11$ である.(sympyで検算済み.)
8.6.2 3次元:体積
次に3次元で考える.ベクトル $\vec{f}=\begin{bmatrix}1\\2\\3\end{bmatrix}$,$\vec{g}=\begin{bmatrix}2\\-1\\1\end{bmatrix}$,$\vec{h}=\begin{bmatrix}1\\1\\2\end{bmatrix}$ を考える.この $\vec f,\vec g,\vec h$ は,じつはすべて同じ平面 $S$ 上にある——すなわち線形従属(linearly dependent,どれか1本が他の2本の実数倍の和で書けること.厳密な定義は第10章で学ぶ)である.実際,$\vec h=\dfrac{3}{5}\vec f+\dfrac{1}{5}\vec g$ という関係が成り立つ.$\vec f,\vec g,\vec h$ を列として合わせると,$C=\begin{bmatrix}1&2&1\\2&-1&1\\3&1&2\end{bmatrix}$ となる.行列式は
$$ |C| = -2+6+2+3-8-1 = 0 $$である(サラスの公式で計算.各自確かめよ).図8.8に,$\vec f,\vec g,\vec h$ が同一平面 $S$($\vec f,\vec g$ が張る平面)上にあり,原点 $O$ からの3本のベクトルが平行六面体を作らずにつぶれてしまう様子を示す.
3次元グラフでは,行列式が $0$ でなければ平行六面体の体積の値を示し,$0$ であれば,3本のベクトルが同一平面上に収まってしまい,体積が$0$である面が表れる.2次元のときと同様,一般の場合を証明しておこう.高校数学のベクトル(空間ベクトル)で登場する内積(dot product,$\vec u\cdot\vec v=u_1v_1+u_2v_2+u_3v_3$)を使う.
定理8.5 行列式の幾何学的意味(3次元)
ベクトル $\vec f=(f_1,f_2,f_3)$,$\vec g=(g_1,g_2,g_3)$,$\vec h=(h_1,h_2,h_3)$ を列として並べた行列を $C$ とする.このとき,$\vec f,\vec g,\vec h$ が作る平行六面体の体積は $|C|$ の絶対値に等しい.とくに $|C|=0$ であることと,$\vec f,\vec g,\vec h$ が同一平面上にある(線形従属である)こととは同値である.
証明(概略)
平行六面体の体積は,底面($\vec f,\vec g$ が作る平行四辺形)の面積 $S$ と,高さ($\vec h$ のうち底面に垂直な方向の成分の大きさ)の積である.空間ベクトルには,2本のベクトル $\vec f,\vec g$ の両方に垂直で,長さが底面積 $S$ に等しいベクトル $\vec f\times\vec g$(外積,cross product)が存在することが知られている(成分では $\vec f\times\vec g=(f_2g_3-f_3g_2,\ f_3g_1-f_1g_3,\ f_1g_2-f_2g_1)$).このとき,体積は $\vec h$ を $\vec f\times\vec g$ の方向に正射影した長さと $S$ の積,すなわち内積を使って $V=|\vec h\cdot(\vec f\times\vec g)|$ と書ける($\vec f\times\vec g$ の長さがすでに $S$ を含んでいるため,正射影の長さ×$S$ は内積そのものになる).そこで $\vec h\cdot(\vec f\times\vec g)$ を成分で計算すると,
$$ \vec h\cdot(\vec f\times\vec g) = h_1(f_2g_3-f_3g_2)+h_2(f_3g_1-f_1g_3)+h_3(f_1g_2-f_2g_1) $$となる.他方,$|C|$ を第3列($\vec h$ の列)に関して余因子展開(定理8.2)すると,
$$ |C| = h_1\begin{vmatrix}f_2&g_2\\f_3&g_3\end{vmatrix}-h_2\begin{vmatrix}f_1&g_1\\f_3&g_3\end{vmatrix}+h_3\begin{vmatrix}f_1&g_1\\f_2&g_2\end{vmatrix} = h_1(f_2g_3-f_3g_2)-h_2(f_1g_3-f_3g_1)+h_3(f_1g_2-f_2g_1) $$であり,$-h_2(f_1g_3-f_3g_1)=h_2(f_3g_1-f_1g_3)$ に書き直せば,これは上の $\vec h\cdot(\vec f\times\vec g)$ の式とまったく同じである.よって $V=|\vec h\cdot(\vec f\times\vec g)|=|C|$.また,$\vec f,\vec g,\vec h$ が同一平面上にある(線形従属である)ことと,高さが $0$,すなわち $V=0$,すなわち $|C|=0$ となることは同値である.
(証明終わり)
例題8.8 平行六面体の体積と線形従属性の判定
$\vec u=(1,2,0)$,$\vec v=(0,1,3)$,$\vec w=(2,0,1)$ は線形従属か.また,これらが作る平行六面体の体積を求めよ.
解答 $\vec u,\vec v,\vec w$ を列として並べた行列式を計算する.
$$ \begin{vmatrix}1&0&2\\2&1&0\\0&3&1\end{vmatrix} = 1\times(1\times1-0\times3)-0\times(2\times1-0\times0)+2\times(2\times3-1\times0) = 1-0+12=13 $$$|{\,}|=13\neq0$ なので,定理8.5より $\vec u,\vec v,\vec w$ は線形従属ではない(線形独立である).また,これらが作る平行六面体の体積は $|13|=13$ である.(sympyで検算済み.)
8.7 まとめと演習
8.7.1 まとめ
- 行列の和・スカラー倍・積(定義8.2,定義8.3)を確認した.積は交換法則が成り立たない($AB\neq BA$).
- 転置行列 ${}^{t}A$(定義8.4)は行と列を入れ換えたもので,${}^{t}({}^{t}A)=A$,${}^{t}(AB)={}^{t}B\,{}^{t}A$(定理8.1,積の転置は順序が逆になる).
- 対角行列・上三角行列・回転行列の3種類は,べき乗が数学的帰納法によってきれいな形(公式8.1〜8.3)で求まる.この理由は第12章の対角化で明らかになる.
- 置換(定義8.5)は $1,\ldots,n$ の並べ替えで,全部で $n!$ 通りある.逆置換(定義8.6)は上下の段を交換して上段を昇順に並べ直せば求まる.偶置換・奇置換(定義8.7)は,恒等置換からの互換の回数の偶奇で決まり,符号 $\sgn$(定義8.8)は偶置換で $+1$,奇置換で $-1$ である.
- 行列式(定義8.9)は,各行・各列から1つずつ成分を選ぶすべての選び方(置換)について,符号付きの積を全部足したもの(ライプニッツの公式).$3$次の場合はサラスの公式に一致する(式\eqref{eq:8-det3-sarrus}).
- 余因子 $A_{ij}$(定義8.10)を使うと,行列式は好きな行・好きな列に関して余因子展開(定理8.2)で計算できる.
- 行列式の基本性質(定理8.3)——(i)行のスカラー倍,(ii)行の入れ替えで符号反転,(iii)他の行の定数倍を加えても不変,(iv)多重線形性——を使うと,行を上三角形に変形して対角成分の積を求めるだけで,大きな行列式も効率的に計算できる.
- 行列式の幾何学的意味(定理8.4,定理8.5):2次元では2本のベクトルが作る平行四辺形の面積,3次元では3本のベクトルが作る平行六面体の体積に等しい.行列式が $0$ になることは,ベクトルたちが「潰れて」次元が落ちる(線形従属になる)ことを意味する.
例題8.9 行列の性質のまとめ確認——転置・行列式・余因子展開・幾何学的意味
$A=\begin{bmatrix}1&2&-1\\0&3&2\\1&0&4\end{bmatrix}$ について,(a) 転置行列 ${}^{t}A$ を求めよ.(b) サラスの公式(式\eqref{eq:8-det3-sarrus})を使って行列式 $|A|$ を求めよ.(c) 第1行に関する余因子展開(定理8.2,式\eqref{eq:8-cofactor-expansion-row})で $|A|$ を計算し,(b)の値と一致することを確認せよ.(d) $A$ の3つの列ベクトルが作る平行六面体の体積を求め,これら3本のベクトルが線形従属か線形独立かを判定せよ.
解答 (a) 定義8.4(行と列を入れ換える)に従うと,
$$ {}^{t}A=\begin{bmatrix}1&0&1\\2&3&0\\-1&2&4\end{bmatrix} $$である($A$ の第1行 $(1,2,-1)$ が ${}^{t}A$ の第1列に,第2行 $(0,3,2)$ が第2列に,第3行 $(1,0,4)$ が第3列になる).
(b) $a_{11}=1,a_{12}=2,a_{13}=-1,a_{21}=0,a_{22}=3,a_{23}=2,a_{31}=1,a_{32}=0,a_{33}=4$ を式\eqref{eq:8-det3-sarrus}に代入する.
$$ \begin{aligned} |A| &= a_{11}a_{22}a_{33} - a_{11}a_{23}a_{32} - a_{12}a_{21}a_{33} + a_{12}a_{23}a_{31} + a_{13}a_{21}a_{32} - a_{13}a_{22}a_{31}\\ &= 1\times3\times4 - 1\times2\times0 - 2\times0\times4 + 2\times2\times1 + (-1)\times0\times0 - (-1)\times3\times1\\ &= 12-0-0+4+0+3\\ &= 19 \end{aligned} $$よって $|A|=19$ である.
(c) 第1行 $(1,2,-1)$ に関して余因子展開する.$A_{11}$ は第1行・第1列を除いた行列式に $(-1)^{1+1}=1$ を,$A_{12}$ は第1行・第2列を除いた行列式に $(-1)^{1+2}=-1$ を,$A_{13}$ は第1行・第3列を除いた行列式に $(-1)^{1+3}=1$ を掛けたものである(定義8.10).
$$ A_{11}=(-1)^{1+1}\begin{vmatrix}3&2\\0&4\end{vmatrix}=3\times4-2\times0=12,\qquad A_{12}=(-1)^{1+2}\begin{vmatrix}0&2\\1&4\end{vmatrix}=-(0\times4-2\times1)=2,\qquad A_{13}=(-1)^{1+3}\begin{vmatrix}0&3\\1&0\end{vmatrix}=0\times0-3\times1=-3 $$したがって,
$$ |A| = a_{11}A_{11}+a_{12}A_{12}+a_{13}A_{13} = 1\times12+2\times2+(-1)\times(-3) = 12+4+3 = 19 $$となり,たしかに(b)の値 $19$ と一致する.
(d) $A$ の3つの列ベクトル $\vec{x}_1=(1,0,1)$,$\vec{x}_2=(2,3,0)$,$\vec{x}_3=(-1,2,4)$ が作る平行六面体の体積は,定理8.5より $|A|$ の絶対値に等しいから,$|19|=19$ である.体積が $0$ でないので,定理8.5より $\vec{x}_1,\vec{x}_2,\vec{x}_3$ は線形従属ではない(線形独立である).((a)〜(d)すべてsympyで検算済み.)
8.7.2 演習問題
演習8.1 転置の性質
$A=\begin{bmatrix}1&-2&3\\0&4&-1\end{bmatrix}$,$B=\begin{bmatrix}2&1&0\\-3&2&5\end{bmatrix}$ に対して,${}^{t}(A+B)={}^{t}A+{}^{t}B$ が成り立つことを確かめよ.また,一般の $(m,n)$ 型行列 $A,B$ についてこの等式が成り立つことを証明せよ.
ヒント:$A+B$(定義8.2)の $(i,j)$ 成分は $a_{ij}+b_{ij}$ である.これを転置してから,${}^{t}A$,${}^{t}B$ の対応する成分の和と比べよ.
演習8.2 少し複雑なべき乗パターンの予想と証明
$M=\begin{bmatrix}2&3\\0&2\end{bmatrix}$ について,$M^{n}$($n$ は自然数)のパターンを,$M^2,M^3$ を計算して予想し,数学的帰納法で証明せよ.
ヒント:$M=2I+N$,$N=\begin{bmatrix}0&3\\0&0\end{bmatrix}$ と分けて考えると見通しがよい.$N^2=O$(零行列)であることに注目せよ.
演習8.3 $n=4$ の置換の偶奇分類
$n=4$ の置換は全部で $4!=24$ 通りある.これらをすべて書き出し,偶置換と奇置換に分類せよ(それぞれ何通りずつあるか).
ヒント:恒等置換 $\begin{pmatrix}1&2&3&4\\1&2&3&4\end{pmatrix}$ から1回の互換で作れる置換(下段の2数を1回入れ替えるだけ)はすべて奇置換である.そこからさらに1回互換すれば偶置換に戻る.8.2.3項の転倒数の数え方(例題8.3)を使うと1つずつ確実に判定できる.
演習8.4 余因子展開による$4$次行列式の計算
$Q=\begin{bmatrix}3&0&1&2\\1&-1&0&4\\0&2&1&0\\2&1&-1&3\end{bmatrix}$ の行列式 $|Q|$ を,余因子展開(定理8.2)を使って求めよ.
ヒント:第3行 $(0,2,1,0)$ に $0$ が2個あるので,第3行に関して展開すると計算量が減る.
演習8.5 行の性質を用いた$4$次行列式の効率的な計算
$R=\begin{bmatrix}1&2&1&0\\2&4&0&1\\1&1&2&2\\0&1&1&3\end{bmatrix}$ の行列式 $|R|$ を,定理8.3の性質(行の基本変形)を使って上三角形に変形してから求めよ.
ヒント:まず $R_2\leftarrow R_2-2R_1$,$R_3\leftarrow R_3-R_1$ とする.そのあと,都合よく非ゼロの成分を対角線上に持ってくるために,行の入れ替え(性質(ii),符号が反転することに注意)が必要になる場面がある.
演習8.6 線形従属性の判定と行列式の幾何学的意味
$\vec p=(1,1,2)$,$\vec q=(2,-1,1)$,$\vec r=(0,3,3)$ は線形従属か.行列式を計算して判定し,定理8.5の意味(体積)と関連付けて説明せよ.
ヒント:線形従属であれば,$\vec r$ を $\vec p,\vec q$ の実数倍の和として具体的に表せるはずである.行列式が $0$ になることと,平行六面体が潰れて体積が $0$ になることの対応を思い出せ.
8.7.3 参考文献
- 望月泰英『数学ノート 線形代数学』(手書き講義ノート).本章の底本.
- 斎藤正彦『線型代数入門』東京大学出版会,1966年.
- 佐武一郎『線型代数学』裳華房,1958年.
関連シミュレーター:行列式(面積・体積としての意味)