正方行列に対して定義される「行列式」というスカラー量について,その定義と計算方法について解説します。
置換を用いた行列式(det)の定義
定義(行列式)
n 次正方行列 A=(aij) に対し,
detA=σ∈Sn∑(sgnσ)a1σ(1)a2σ(2)…anσ(n)
を行列式 (determinant) という。A の行列式は,detA や ∣A∣ などとかく。
ここで,Sn は置換全体を集めた対称群であり,∑σ∈Sn は,それら全部を足し合わせることを指します。
また,sgnσ は置換の符号を指します。
行列式の具体例
2次正方行列の行列式
2次正方行列の行列式
∣∣a11a21a12a22∣∣=a11a22−a12a21.
第1項・第2項の置換はそれぞれ (1122),(1221) であり,置換の符号は 1,−1 となります。
3次正方行列の行列式
3次正方行列の行列式
=∣∣a11a21a31a12a22a32a13a23a33∣∣a11a22a33−a11a23a32+a12a23a31−a12a21a33+a13a21a32−a13a22a31.
各項の置換は順に
(112233),(112332),(122331),(122133),(132132),(132231)
に対応しています。3次の行列式は覚えなくてもよいですが,サラスの方法という暗記法があります。以下の記事を参照してください。
三角行列の行列式
三角行列の行列式
上三角行列
A=⎝⎛a110a12a22……⋱a1na2n⋮ann⎠⎞
の行列式は detA=a11a22…ann である。
下三角行列も同じである。
a11a22…ann 以外の項は積の中に必ず 0 が含まれるため,消えるわけです。
行列式の計算手順
行列式の計算は,もちろん定義通り計算しても構いません。しかし,たとえば 4×4 行列の場合,足し算の数は 4!=24 個にもなってしまうため,あまり現実的ではありません。そこで,現実的な計算手順について,順を追って解説します。
行列式の現実的な計算手順
現実的な計算手順は,順に
- 行列の基本変形を用いて0を増やす(特に1行目or1列目)
- 1行目or1列目に注目して展開する
- 得られた小さな行列式に対して手順1-2を繰り返す
- 2次・3次または三角行列の行列式に帰着させ,計算する
です。それぞれ解説していきましょう。
1. 行列の基本変形を用いて0を増やす(特に1行目or1列目)
まずは行列の基本変形を用いて,1行目 or 1列目の成分がある程度 0 になるようにしましょう。行列の基本変形とは,次の3つの変形です。
行列の基本変形
- ある行(列)の c 倍を他の行(列)に加えること
- 2つの行(列)を入れ替えること
- ある行(列)を c=0 倍すること
詳しくは,行列の基本変形についてわかりやすく図解するで解説しています。
基本変形をしたときの行列式は,以下のように変化します。
基本変形 | 行列式 |
---|
ある行(列)の c 倍を他の行(列)に加える | 変化なし |
2つの行(列)を入れ替える | −1 倍 |
ある行(列)を c=0 倍する | c 倍 |
このことを留意して,係数を調整しながら,ある程度計算しやすいところまで基本変形していきましょう。
注意ですが,このステップでは無理に計算しすぎなくて良いです。頑張って計算すれば三角行列にできるため,行列式が求まりますが,たとえば各成分に変数が含まれている場合,複雑な計算をすることになってしまうため,そこまでする必要はありません。ある程度同じ行・列に 0 が集まればOKと思い,次のステップに進みましょう。
2. 1行目or1列目に注目して展開する
行列式には,以下のような性質があります。
定理(行列式に関する等式)
A を n 次正方行列とする。
Aij を A から i 行と j 列を取り除いた n−1 次正方行列の行列式とすると,
detA=a11A11−a12A12+a13A13−⋯+(−1)n−1a1nA1n=a11A11−a21A21+a31A31−⋯+(−1)n−1an1An1.
(−1)i+jAij を余因子と言います。これを図にすると,以下のようになります。
これ自身の証明は,余因子行列の定義と余因子展開~逆行列になる証明~で行っています。これを用いて,行列式を展開しましょう。
このとき,手順1の効果で, a1j=0 (または ai1=0) の項は消えますね。
なお,一般に1行目or1列目以外の行・列に注目して展開することも可能ですが,手順1で「2つの行(列)を入れ替えること」を行って,1行目or1列目に移動させてから展開すると考えればよいでしょう。
3. 得られた小さな行列式に対して手順1-2を繰り返す
手順2を行うと,小さな行列式がたくさん出てくると思います。これについて,再び手順1-2を繰り返しましょう。
4. 2次・3次または三角行列の行列式に帰着させ,計算する
2次・3次・三角行列の行列式にしてしまえば,具体例のところで挙げたとおり,計算できますね。
長いステップですが,実際に手計算するのは3次から5次程度の行列式だと思うので,「何度も手順1-2を繰り返さねばならない」ということはないでしょう。
行列式の計算の例題
実際に,以下の問題に取り組んでみましょう。
例題
∣∣113275012215440−3∣∣ を計算せよ。
計算手順を復習すると,
- 行列の基本変形を用いて0を増やす(特に1行目or1列目)
- 1行目or1列目に注目して展開する
- 得られた小さな行列式に対して手順1-2を繰り返す
- 2次・3次または三角行列の行列式に帰着させ,計算する
でした。これを用いて手順通りに計算すると,解答例は以下のようになります。最短手順とは限らず,別の変形をしても構いません。係数に気を付けましょう。
本記事では,例は一つしか紹介しませんが,もっと多くの例を知りたい場合は,たとえば固有値の定義と求め方をていねいに~計算の手順~の例題も見てみてください。
「固有値」の概念は,現時点で分かる必要はないですが,計算例題は,結局「変数を含む行列式の計算」ですから,そこは理解できるでしょう。
計算機やコンピューターを利用するのもアリ
数学科の学生が行列式を手計算するのは,線形代数のテストや院試のときくらいです(もちろん,計算機で計算できない抽象的な行列を研究する場合はその限りでないです)。
応用数学の研究などでは,計算機を用いて計算するのが一般的です。たとえば,WolframAlphaを用いれば,インターネット上でも計算することができます。本格的に計算するためには,MathematicaやPython, Julia, R, Fortran, Mapleなどのプログラミング言語を利用すればよいでしょう。
行列式の性質
行列式の性質は,以下が有名です。
以下で,複素数 C は実数 R に変えても成立します。
定理(行列式の性質)
1. 行列式は,各列に関して線形性がある。すなわち,ak∈Cn(1≤k≤n),bk∈Cn を n 次元列ベクトル,s,t∈C とすると,
=det(a1,…,ak−1,sak+tbk,ak+1,…,an)sdet(a1,…,ak−1,ak,ak+1,…,an)+tdet(a1,…,ak−1,bk,ak+1,…,an).
また同様に,行列式は各行に関して線形性がある。すなわち,ak∈Cn(1≤k≤n),bk∈Cn を n 次元行ベクトル,s,t∈C とすると,
det⎝⎛a1⋮ak−1sak+tbkak+1⋮an⎠⎞=sdet⎝⎛a1⋮ak−1akak+1⋮an⎠⎞+tdet⎝⎛a1⋮ak−1bkak+1⋮an⎠⎞.
2. 転置行列の行列式は元のそれと同じ。すなわち,detA=detA⊤.
3. (列・行の入れ替え)σ∈Sn を置換とする。このとき,太字を列ベクトルとすると,
det(aσ(1),…,aσ(n))=(sgnσ)det(a1,…,an).
太字を行ベクトルとすると,
det⎝⎛aσ(1)⋮aσ(n)⎠⎞=(sgnσ)det⎝⎛a1⋮an⎠⎞.
特に, 2つの列(または行)を入れ替えると,行列式は −1 倍になる。
4. 2つの列(または行)が等しい行列式の値は 0 である。
5. 行列式と積は可換である。すなわち,A,B を n 次正方行列とするとき,
det(AB)=detAdetB.
6. A が正則のとき,det(A−1)=(detA)−1.
これらの証明については,以下の記事を参照してください。
おわりに
行列式の計算は,線形代数学のテストで頻出です。しっかりと演習を積んで,テストに臨みましょう。
特別な行列式