MathLabs

代数学

行列と行列式

数を長方形状に並べたものと、それが体積をどれだけ拡大縮小するかを測るスカラー量。

直観直感:空間を変形させる装置としての行列

行列は数を長方形に並べたものであるが、その最大の有用性は、ある一定の線形規則によって一つのベクトルを別のベクトルに変える働きにある:ここでは伸ばし、あそこでは回転させ、時には反転させる。行列式は正方行列に結び付けられた一つの数であり、非常に具体的な問いに答える:この行列は面積(平面内)または体積(空間内)をどの倍率で拡大縮小するか、そしてその過程で向き(左手系が右手系になる、など)を反転させるかどうか。行列式が 00 であるとは、行列が空間を完全に押しつぶし、一つの次元を丸ごと失わせることを意味する——この一つの事実が、その行列が「元に戻せる」かどうかを含め、他のほとんどすべての性質を左右する。

線形変換を示す、4 次元超立方体の回転する対話的な射影図。
2×22\times 2 行列 A=(a11a12a21a22)A = \begin{pmatrix}a_{11}&a_{12}\\a_{21}&a_{22}\end{pmatrix} による平面の変形:単位正方形は符号付き面積 det⁡A=a11a22−a12a21\det A = a_{11}a_{22} - a_{12}a_{21} の平行四辺形に移る。

中高行列の乗法と行列式の公式

定義: 行列の乗法

AA が m×nm\times n、BB が n×pn\times p のとき、その積 C=ABC=AB は成分が Cij=∑k=1nAikBkjC_{ij} = \sum_{k=1}^{n} A_{ik}B_{kj} で与えられる m×pm\times p 行列である:CC の成分 (i,j)(i,j) は AA の第 ii 行と BB の第 jj 列の内積である。

Cij=∑k=1nAikBkjC_{ij} = \sum_{k=1}^{n} A_{ik}B_{kj}

ここで AikA_{ik} は AA の第 ii 行第 kk 列の成分、BkjB_{kj} は BB の第 kk 行第 jj 列の成分である。和は共通の内側添字 kk にわたって取られ、これが積が定義されるために AA の列数と BB の行数が一致しなければならない理由である。

定義: 行列式(ライプニッツの公式)

n×nn\times n の正方行列 AA について、行列式は det⁡(A)=∑σ∈Snsgn⁡(σ)∏i=1nai,σ(i)\det(A)=\sum_{\sigma\in S_n}\operatorname{sgn}(\sigma)\prod_{i=1}^n a_{i,\sigma(i)} によって定義される。これは {1,…,n}\{1,\dots,n\} のすべての置換 σ\sigma にわたる和であり、σ\sigma が偶置換か奇置換かに応じて sgn⁡(σ)=+1\operatorname{sgn}(\sigma)=+1 または −1-1 となる。n=2,3n=2,3 の場合、これは以下のおなじみの余因子展開の公式に展開される。

det⁡(abcd)=ad−bc\det\begin{pmatrix}a&b\\ c&d\end{pmatrix}=ad-bc

2×22\times 2 行列について、det⁡(abcd)=ad−bc\det\begin{pmatrix}a&b\\ c&d\end{pmatrix}=ad-bc はまさに行ベクトル (a,b)(a,b) と (c,d)(c,d) が張る平行四辺形の符号付き面積である:{1,2}\{1,2\} の 2 つの置換とは、恒等置換(+ad+ad を与える)と交換(−bc-bc を与える)である。

det⁡(a11a12a13a21a22a23a31a32a33)=a11(a22a33−a23a32)−a12(a21a33−a23a31)+a13(a21a32−a22a31)\det\begin{pmatrix}a_{11}&a_{12}&a_{13}\\ a_{21}&a_{22}&a_{23}\\ a_{31}&a_{32}&a_{33}\end{pmatrix}=a_{11}(a_{22}a_{33}-a_{23}a_{32})-a_{12}(a_{21}a_{33}-a_{23}a_{31})+a_{13}(a_{21}a_{32}-a_{22}a_{31})

3×33\times 3 行列について、det⁡(a11a12a13a21a22a23a31a32a33)=a11(a22a33−a23a32)−a12(a21a33−a23a31)+a13(a21a32−a22a31)\det\begin{pmatrix}a_{11}&a_{12}&a_{13}\\ a_{21}&a_{22}&a_{23}\\ a_{31}&a_{32}&a_{33}\end{pmatrix}=a_{11}(a_{22}a_{33}-a_{23}a_{32})-a_{12}(a_{21}a_{33}-a_{23}a_{31})+a_{13}(a_{21}a_{32}-a_{22}a_{31}) は第 1 行に沿った余因子展開である:第 11 行に沿って展開し、各成分 a1ja_{1j} に、行 11 と列 jj を除いた後に残る 2×22\times 2 行列の行列式を掛け、j=1,2,3j=1,2,3 にわたって符号を +,−,++,-,+ と交互に変える。

基本行操作が行列式に与える影響
行操作det⁡\det への影響
2 つの行を入れ替えるdet⁡\det に −1-1 を掛ける
1 つの行を kk 倍するdet⁡\det に kk を掛ける
ある行の kk 倍を別の行に加えるdet⁡\det は変化しない
2 つの行が等しいか比例するdet⁡=0\det=0

大学定理:行列式の乗法性と正則性の判定条件

任意の n×nn\times n 行列 A,BA,B について、det⁡(AB)=det⁡(A)det⁡(B)\det(AB)=\det(A)\det(B) が成り立つ。

なぜ正しいのか?

行列式は線形写像が体積をどれだけ拡大縮小するかを測るものなので、まず写像 BB を適用し、その後写像 AA を適用すると、体積は BB の倍率で、続いて AA の倍率でスケールされるはずである、すなわち 2 つの倍率の積になるはずである——乗法性とは、線形写像を合成すると、それらの体積スケーリング係数も合成されるという事実の代数的な表現にすぎない。

証明

まず、3 つの基本行操作が行列式に与える影響を振り返る(これはライプニッツの和 det⁡(A)=∑σ∈Snsgn⁡(σ)∏i=1nai,σ(i)\det(A)=\sum_{\sigma\in S_n}\operatorname{sgn}(\sigma)\prod_{i=1}^n a_{i,\sigma(i)} から直接確立される。この和の各項は各行について別々に線形だからである):2 つの行を入れ替えると det⁡\det に −1-1 を掛け、1 つの行を kk 倍すると det⁡\det に kk を掛け、1 つの行の何倍かを別の行に加えても det⁡\det は変化しない。

AA に対する各基本行操作は、AA に対応する基本行列 EE を左から掛けることと同じである(単位行列 II に同じ操作を適用して得られる)。前段落と比較すると、det⁡(E)\det(E) はちょうどその操作のスケーリング係数(−1-1、kk、または 11)に等しいので、すべての基本行列 EE について det⁡(EA)=det⁡(E)det⁡(A)\det(EA)=\det(E)\det(A) が成り立つ。

AA が正則な場合:ガウスの消去法は、有限個の基本行操作の列を用いて、任意の正則行列を単位行列に変形する。すなわち、基本行列 E1,…,EmE_1,\dots,E_m に対して Em⋯E1A=IE_m\cdots E_1A=I となるので、A=E1−1⋯Em−1A=E_1^{-1}\cdots E_m^{-1} となり、これ自体が基本行列の積である(基本行列の逆行列も同じ種類の基本行列である)。前段落の等式を繰り返し適用すると det⁡(A)=det⁡(E1−1)⋯det⁡(Em−1)\det(A)=\det(E_1^{-1})\cdots\det(E_m^{-1}) となり、これを再び AB=E1−1⋯Em−1BAB=E_1^{-1}\cdots E_m^{-1}B に適用すると det⁡(AB)=det⁡(E1−1)⋯det⁡(Em−1)det⁡(B)=det⁡(A)det⁡(B)\det(AB)=\det(E_1^{-1})\cdots\det(E_m^{-1})\det(B)=\det(A)\det(B) が得られる。

AA が特異な場合:このとき det⁡(A)=0\det(A)=0 である(特異行列は単位行列まで完全には変形できず、上記の行操作の規則から、到達可能な行階段形には常に零行が残ることが分かり、その行に沿った余因子展開により det⁡=0\det=0 が強制される)。特異であることは rank⁡(A)<n\operatorname{rank}(A)<n を意味し、任意の BB に対して rank⁡(AB)≤rank⁡(A)<n\operatorname{rank}(AB)\leq\operatorname{rank}(A)<n であるから、積 ABAB も特異であり、したがって det⁡(AB)=0\det(AB)=0 となる。よって det⁡(AB)=0=0⋅det⁡(B)=det⁡(A)det⁡(B)\det(AB)=0=0\cdot\det(B)=\det(A)\det(B) となり、この場合も等式は成り立ち、すべての AA について証明が完成する。

正方行列 AA が正則であるのは det⁡(A)≠0\det(A)\neq 0 のとき、かつそのときに限る。正則であるとき、A−1=1det⁡(A)adj⁡(A)A^{-1}=\dfrac{1}{\det(A)}\operatorname{adj}(A) が成り立つ。ここで adj⁡(A)\operatorname{adj}(A) は余因子行列(余因子からなる行列の転置)である。

なぜ正しいのか?

線形写像が「元に戻せる」のは、空間をより低い次元に押しつぶさない場合に限られる——行列式はまさに「体積をゼロに押しつぶすかどうか」という問いを測るものなので、正則性と det⁡≠0\det\neq 0 は、実は 2 つの角度から見た同じ主張である:一つは代数的な角度(逆行列が存在するか)、もう一つは幾何学的な角度(体積スケーリング係数がゼロでないか)。

証明

(⇒\Rightarrow)AA が正則であると仮定すると、AA−1=IAA^{-1}=I を満たす行列 A−1A^{-1} が存在する。上で証明した乗法性の定理を適用すると、det⁡(A)det⁡(A−1)=det⁡(AA−1)=det⁡(I)=1\det(A)\det(A^{-1})=\det(AA^{-1})=\det(I)=1 となる。det⁡(A)\det(A) と det⁡(A−1)\det(A^{-1}) という 2 つの数の積が 11 に等しいので、どちらも 00 にはなり得ない——よって det⁡(A)≠0\det(A)\neq 0 である。

(⇐\Leftarrow)、2×22\times 2 で具体的に示す:A=(abcd)A=\begin{pmatrix}a&b\\ c&d\end{pmatrix} が det⁡(A)=ad−bc≠0\det(A)=ad-bc\neq 0 を満たすとする。adj⁡(A)=(d−b−ca)\operatorname{adj}(A)=\begin{pmatrix}d&-b\\ -c&a\end{pmatrix} と定義する。直接計算すると、A⋅adj⁡(A)=(abcd)(d−b−ca)=(ad−bc−ab+abcd−cd−cb+ad)=(ad−bc00ad−bc)=det⁡(A) IA\cdot\operatorname{adj}(A)=\begin{pmatrix}a&b\\ c&d\end{pmatrix}\begin{pmatrix}d&-b\\ -c&a\end{pmatrix}=\begin{pmatrix}ad-bc&-ab+ab\\ cd-cd&-cb+ad\end{pmatrix}=\begin{pmatrix}ad-bc&0\\ 0&ad-bc\end{pmatrix}=\det(A)\,I となり、これはちょうど等式 A⋅adj⁡(A)=det⁡(A) IA\cdot\operatorname{adj}(A)=\det(A)\,I である。逆の順序で同じ計算をすると adj⁡(A)⋅A=det⁡(A) I\operatorname{adj}(A)\cdot A=\det(A)\,I も得られる。

det⁡(A)≠0\det(A)\neq 0 なので、A⋅adj⁡(A)=det⁡(A)IA\cdot\operatorname{adj}(A)=\det(A)I の両辺をスカラー det⁡(A)\det(A) で割ると A⋅adj⁡(A)det⁡(A)=IA\cdot\dfrac{\operatorname{adj}(A)}{\det(A)}=I となり、同様に逆の順序から adj⁡(A)det⁡(A)⋅A=I\dfrac{\operatorname{adj}(A)}{\det(A)}\cdot A=I が得られる。これにより明示的な両側逆行列 A−1=1det⁡(A)adj⁡(A)A^{-1}=\dfrac{1}{\det(A)}\operatorname{adj}(A) が示され、AA は正則である。

一般の nn の場合(n=3n=3 で例示)、同じ議論が成り立つ。ラプラスの余因子展開により、∑jaijCij=det⁡(A)\sum_j a_{ij}C_{ij}=\det(A)(第 ii 行に沿った展開)であり、i≠ki\neq k のとき ∑jaijCkj=0\sum_j a_{ij}C_{kj}=0 となる(この和は、第 ii 行を第 kk 行にコピーした行列の余因子展開であり、2 つの行が等しいため上の表により行列式が 00 になる)。これらの行ごとの等式を行列の形にまとめると、任意の nn について、ちょうど A⋅adj⁡(A)=det⁡(A)IA\cdot\operatorname{adj}(A)=\det(A)I が得られる。したがって det⁡(A)≠0\det(A)\neq 0 である限り、2×22\times 2 の場合と同じ除算の議論により A−1=adj⁡(A)/det⁡(A)A^{-1}=\operatorname{adj}(A)/\det(A) が得られる。逆に det⁡(A)=0\det(A)=0 の場合、第一段落の対偏を用いれば正則性は完全に排除される。なぜなら det⁡(A)det⁡(A−1)=1\det(A)\det(A^{-1})=1 は det⁡(A)≠0\det(A)\neq 0 を要求するからである。これにより判定条件が両方向で証明される。

大学実世界での応用と具体例

行列式を用いると、エンジニアは消去法の代わりに明示的な公式(クラメルの公式)で小規模な線形連立方程式を解くことができ、コンピュータグラフィックスのプログラムは、変換が向きを保つか反転させるか、そして面積をどれだけ拡大縮小するかを、たった一つの数で確認できる——どちらも det⁡≠0\det\neq 0 が「正則」を意味し、∣det⁡∣|\det| が「拡大縮小率」を意味することの直接的な応用である。

例: 工学:クラメルの公式で 2 ループ回路を解く

2 ループ回路にキルヒホッフの法則を適用すると、2 つのループ電流 x,yx,y(単位はアンペア)について線形連立方程式 3x+y=113x+y=11、x+2y=8x+2y=8 が得られる。行列式を用いて xx と yy を求めよ。

解答

連立方程式を行列方程式として書く:(3112)(xy)=(118)\begin{pmatrix}3&1\\1&2\end{pmatrix}\begin{pmatrix}x\\y\end{pmatrix}=\begin{pmatrix}11\\8\end{pmatrix}。係数行列は det⁡=3×2−1×1=6−1=5≠0\det=3\times2-1\times1=6-1=5\neq 0 を満たすので、正則性の判定条件により連立方程式は一意な解を持ち、クラメルの公式が適用できる。

クラメルの公式によれば x=det⁡(11182)det⁡(3112)x=\dfrac{\det\begin{pmatrix}11&1\\8&2\end{pmatrix}}{\det\begin{pmatrix}3&1\\1&2\end{pmatrix}} である:係数行列の第 1 列を右辺で置き換え、その行列式を取り、元の行列式で割る。分子は 11×2−1×8=22−8=1411\times2-1\times8=22-8=14 なので、x=14/5=2.8x=14/5=2.8 A である。

同様に y=det⁡(31118)5y=\dfrac{\det\begin{pmatrix}3&11\\1&8\end{pmatrix}}{5} である:今度は第 2 列を置き換える。分子は 3×8−11×1=24−11=133\times8-11\times1=24-11=13 なので、y=13/5=2.6y=13/5=2.6 A である。

検算:3(2.8)+2.6=8.4+2.6=113(2.8)+2.6=8.4+2.6=11 ✓ および 2.8+2(2.6)=2.8+5.2=82.8+2(2.6)=2.8+5.2=8 ✓ となり、x=2.8x=2.8 A、y=2.6y=2.6 A が回路方程式を正確に満たすことが確認できる。

例: コンピュータグラフィックス:変換は図形を反転させるか、また面積を何倍にスケールするか

ある 2D グラフィックスエンジンが、面積 55 平方単位の図形のすべての点に変換行列 M=(2143)M=\begin{pmatrix}2&1\\4&3\end{pmatrix} を適用する。変換後の図形が元の図形に対して鏡映反転しているかどうかを判定し、変換後の図形の面積を求めよ。

解答

MM の行列式を計算する:det⁡(M)=2×3−1×4=6−4=2\det(M)=2\times3-1\times4=6-4=2。

det⁡(M)\det(M) の符号は向きについて教えてくれる:det⁡(M)=2>0\det(M)=2>0 なので、この変換は向きを保つため、図形は鏡映反転していない(行列式が負であれば反転を示す。これは行列に奇数回の反射が組み込まれていることを意味するからである)。

絶対値 ∣det⁡(M)∣|\det(M)| は面積のスケーリング係数を教えてくれる:変換後、どの領域の面積もその形に関わらず ∣det⁡(M)∣=2|\det(M)|=2 倍される。これはまさに、先に確立した行列式の幾何学的意味そのものである。

したがって、変換後の図形の面積は 5×∣det⁡(M)∣=5×2=105\times|\det(M)|=5\times2=10 平方単位であり、元の図形と同じ向き(反転していない)を保つ。

det⁡(2513)\det\begin{pmatrix}2&5\\1&3\end{pmatrix} を求めよ。

2 つの n×nn\times n 行列について det⁡(A)=3\det(A)=3、det⁡(B)=4\det(B)=4 のとき、det⁡(AB)\det(AB) はいくつか。

逆行列が存在しないのはどの行列か。

あるグラフィックス変換行列が det⁡=−3\det=-3 を持つ。これは何を意味するか。

参考文献

  1. Eric W. Weisstein (MathWorld) (2024). Determinant
  2. Gilbert Strang (2016). Introduction to Linear Algebra