§D3.4行列式

最終更新

正方行列が正則であるかどうかは、行基本変形を行って簡約階段形が単位行列に一致するかどうかを調べることで判定することができます。行列式は、行変形を実行しなくても、行列の成分だけから計算することができ、その値が00であるかどうかによって正則性を判定することができる量です。nn次正方行列に対して、行列式は各行について線形であること、二つの行が等しいときに値が00になること、および単位行列InI_nでの値が11であることの三つの性質によって特徴づけられ、det⁡A\det Aと書きます。行列式は線形代数における基本的な量の一つであり、たとえば22次ではdet⁡A=ad−bc\det A = ad - bcという成分の式で与えられます。∣det⁡A∣|\det A|は行列AAによる面積・体積の拡大率としても解釈することができ、この解釈は次の記事「数ベクトルと幾何的な意味」で扱います。本記事では、行列式の基本的な性質や求め方について解説します。

1 2次・3次の計算

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

33次はサラスの方法(右下がりの積を足し、右上がりの積を引く)で計算することができます。これらは、以下で与える一般の定義の特別な場合です(例 4.3)。

2 行列式を特徴づける

nn次正方行列AAを、その行a⃗1,…,a⃗n\vec a_1, \dots, \vec a_nを縦に並べたものとみなし、A=(a⃗1;… ;a⃗n)A = (\vec a_1; \dots; \vec a_n)と書きます。行列式は、このnn個の行の関数として特徴づけます。

定義 2.1 (行列式).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とする。関数D:(Kn)n→KD : (\mathbb{K}^n)^n \to \mathbb{K}が次の3条件を満たすとき、DDをK\mathbb{K}上のnn次の行列式 (determinant) という。

  1. 各行について線形である。すなわち各kkについて、他の行を固定すると、任意のα,β∈K\alpha,\beta\in\mathbb{K}とx⃗,y⃗∈Kn\vec x,\vec y\in\mathbb{K}^nに対して D(… ;αx⃗+βy⃗;… )=αD(… ;x⃗;… )+βD(… ;y⃗;… )D(\dots; \alpha\vec x + \beta\vec y; \dots) = \alpha D(\dots; \vec x; \dots) + \beta D(\dots; \vec y; \dots) が成り立つ。
  2. 二つの行が等しければ値は00である。すなわちp≠qp \ne qかつa⃗p=a⃗q\vec a_p = \vec a_qならばD(A)=0D(A) = 0である。
  3. 単位行列での値は11である。すなわちD(In)=1D(I_n) = 1である。

条件 (a)と条件 (b)を満たす関数を交代多重線形形式 (alternating multilinear form) という。

補題 2.2 (行の入れ替えと符号).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とする。F:(Kn)n→KF : (\mathbb{K}^n)^n \to \mathbb{K}が交代多重線形形式であるとき、二つの行を入れ替えるとFFの値の符号が変わる。

証明.p≠qp \ne qを固定し、第pp行と第qq行の双方にx⃗+y⃗\vec x + \vec yを置いた行列を考える。二つの行が等しいので、定義 2.1 条件 (b)より値は00である。第pp行と第qq行のそれぞれについて定義 2.1 条件 (a)で展開すると

0=F(… ;x⃗;… ;x⃗;… )+F(… ;x⃗;… ;y⃗;… )+F(… ;y⃗;… ;x⃗;… )+F(… ;y⃗;… ;y⃗;… )0 = F(\dots; \vec x; \dots; \vec x; \dots) + F(\dots; \vec x; \dots; \vec y; \dots) + F(\dots; \vec y; \dots; \vec x; \dots) + F(\dots; \vec y; \dots; \vec y; \dots)

となる。右辺の第1項と第4項は二つの行が等しいので00である。したがって

F(… ;x⃗;… ;y⃗;… )=−F(… ;y⃗;… ;x⃗;… )F(\dots; \vec x; \dots; \vec y; \dots) = -F(\dots; \vec y; \dots; \vec x; \dots)

であり、第pp行と第qq行を入れ替えると符号が変わる。▨

定義 2.1 条件 (c)を落としても、他の二つの条件から導かれる帰結は変わらず、失われるのは一意性だけです。この記事では、その二つだけを仮定した主張を交代多重線形形式について述べ、行列式そのものへは定義 2.1 条件 (c)を課して適用します。

3 置換とその符号

行列式の値を成分の式として書き下すには、行の番号と列の番号を対応させる置換が必要になります。

定義 3.1 (置換・互換).{1,2,…,n}\{1, 2, \dots, n\}からそれ自身への全単射を置換 (permutation) といい、置換全体の集合をSnS_nと書く。相異なるp,qp, qについて、ppとqqを入れ替え、他をすべて動かさない置換を互換 (transposition) といい(p  q)(p\;q)と書く。恒等写像を恒等置換 (identity permutation) といいid\mathrm{id}と書く。

定義 3.2 (置換の符号).σ∈Sn\sigma \in S_nに対し

sgn⁡(σ)=∏1≤i<j≤nσ(j)−σ(i)j−i\operatorname{sgn}(\sigma) = \prod_{1 \le i < j \le n} \frac{\sigma(j) - \sigma(i)}{j - i}

をσ\sigmaの符号 (sign) という。

命題 3.3 (符号の基本性質).σ,τ∈Sn\sigma, \tau \in S_nについて次が成り立つ。

  1. sgn⁡(σ)∈{1,−1}\operatorname{sgn}(\sigma) \in \{1, -1\}である。
  2. sgn⁡(στ)=sgn⁡(σ)sgn⁡(τ)\operatorname{sgn}(\sigma\tau) = \operatorname{sgn}(\sigma)\operatorname{sgn}(\tau)である。とくにsgn⁡(σ−1)=sgn⁡(σ)\operatorname{sgn}(\sigma^{-1}) = \operatorname{sgn}(\sigma)である。
  3. 互換の符号は−1-1である。
  4. 任意の置換は有限個の互換の積として表される(恒等置換は00個の積とする)。σ\sigmaがrr個の互換の積で表されるならsgn⁡(σ)=(−1)r\operatorname{sgn}(\sigma) = (-1)^rであり、とくにrrの偶奇は表し方によらない。

証明.(1)を示す。σ\sigmaは全単射なので、対応{i,j}↦{σ(i),σ(j)}\{i, j\} \mapsto \{\sigma(i), \sigma(j)\}は{1,…,n}\{1, \dots, n\}の2元部分集合全体の上の全単射である。したがって

∏i<j∣σ(j)−σ(i)∣=∏i<j∣j−i∣\prod_{i < j} |\sigma(j) - \sigma(i)| = \prod_{i < j} |j - i|

であり、両辺は00でないから∣sgn⁡(σ)∣=1|\operatorname{sgn}(\sigma)| = 1である。

(2)を示す。比(σ(l)−σ(k))/(l−k)(\sigma(l) - \sigma(k)) / (l - k)はkkとllを入れ替えても変わらないので、 2元部分集合{k,l}\{k, l\}だけで決まる。この値をRσ({k,l})R_\sigma(\{k, l\})と書くとsgn⁡(σ)=∏{k,l}Rσ({k,l})\operatorname{sgn}(\sigma) = \prod_{\{k, l\}} R_\sigma(\{k, l\})である。ここで

στ(j)−στ(i)j−i=σ(τ(j))−σ(τ(i))τ(j)−τ(i)⋅τ(j)−τ(i)j−i=Rσ({τ(i),τ(j)})⋅Rτ({i,j})\frac{\sigma\tau(j) - \sigma\tau(i)}{j - i} = \frac{\sigma(\tau(j)) - \sigma(\tau(i))}{\tau(j) - \tau(i)} \cdot \frac{\tau(j) - \tau(i)}{j - i} = R_\sigma(\{\tau(i), \tau(j)\}) \cdot R_\tau(\{i, j\})

であり、{i,j}↦{τ(i),τ(j)}\{i, j\} \mapsto \{\tau(i), \tau(j)\}が2元部分集合全体の上の全単射であることから、すべての2元部分集合にわたる積をとるとsgn⁡(στ)=sgn⁡(σ)sgn⁡(τ)\operatorname{sgn}(\sigma\tau) = \operatorname{sgn}(\sigma)\operatorname{sgn}(\tau)を得る。σσ−1=id\sigma\sigma^{-1} = \mathrm{id}とsgn⁡(id)=1\operatorname{sgn}(\mathrm{id}) = 1および(1)から、sgn⁡(σ−1)=sgn⁡(σ)\operatorname{sgn}(\sigma^{-1}) = \operatorname{sgn}(\sigma)が従う。

(3)を示す。まず隣接する二つの数の互換ρ=(k  k+1)\rho = (k\;k+1)を考え、sgn⁡(ρ)=∏{i,j}Rρ({i,j})\operatorname{sgn}(\rho) = \prod_{\{i, j\}} R_\rho(\{i, j\})を2元部分集合の種類ごとに分けて計算する。

  • {k,k+1}\{k, k+1\}:因子は(ρ(k+1)−ρ(k))/((k+1)−k)=(k−(k+1))/1=−1(\rho(k+1) - \rho(k))/((k+1) - k) = (k - (k+1))/1 = -1である。
  • kkとk+1k+1のどちらも含まない2元部分集合:ρ\rhoは両方の数を動かさないので因子は11である。
  • kkとk+1k+1のちょうど一方を含む2元部分集合:i<ki < kについて{i,k}\{i, k\}と{i,k+1}\{i, k+1\}が組になり、因子の積は (k+1)−ik−i⋅k−i(k+1)−i=1\frac{(k+1) - i}{k - i}\cdot\frac{k - i}{(k+1) - i} = 1 である。j>k+1j > k+1について{k,j}\{k, j\}と{k+1,j}\{k+1, j\}が組になり、因子の積は j−(k+1)j−k⋅j−kj−(k+1)=1\frac{j - (k+1)}{j - k}\cdot\frac{j - k}{j - (k+1)} = 1 である。この2種類の組で、一方だけを含む2元部分集合はすべて尽くされる。

したがってsgn⁡((k  k+1))=−1\operatorname{sgn}((k\;k+1)) = -1である。次にp<qp < qについて

(p  q)=(p  p+1)(p+1  p+2)⋯(q−1  q)⋯(p+1  p+2)(p  p+1)(p\;q) = (p\;p{+}1)(p{+}1\;p{+}2)\cdots(q{-}1\;q)\cdots(p{+}1\;p{+}2)(p\;p{+}1)

と、隣接互換2(q−p)−12(q - p) - 1個の積で表す。(2)によりsgn⁡((p  q))=(−1)2(q−p)−1=−1\operatorname{sgn}((p\;q)) = (-1)^{2(q-p)-1} = -1である。

(4)を示す。σ\sigmaが動かす点の個数m=#{k:σ(k)≠k}m = \#\{k : \sigma(k) \ne k\}についての帰納法による。m=0m = 0ならσ=id\sigma = \mathrm{id}で、00個の互換の積である。m>0m > 0とし、σ(k)=l≠k\sigma(k) = l \ne kとなるkkをとる。τ=(k  l)\tau = (k\;l)とおくとτσ\tau\sigmaはkkを固定する。またσ\sigmaが固定する点はτσ\tau\sigmaも固定する(σ(j)=j\sigma(j) = jかつj≠k,lj \ne k, lならτσ(j)=j\tau\sigma(j) = jであり、j=lj = lはσ(k)=l≠l\sigma(k) = l \ne lと全単射性からσ\sigmaが固定しない)。よってτσ\tau\sigmaが動かす点の個数はmmより真に小さく、帰納法の仮定からτσ\tau\sigmaは互換の積である。τ−1=τ\tau^{-1} = \tauなのでσ=τ(τσ)\sigma = \tau(\tau\sigma)も互換の積である。符号については、(2)と(3)からrr個の互換の積の符号が(−1)r(-1)^rになる。sgn⁡(σ)\operatorname{sgn}(\sigma)は表し方によらない値なので、rrの偶奇も表し方によらない。▨

4 存在と一意性、置換にわたる和

定理 4.1 (交代多重線形形式の決定と置換にわたる和).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とする。F:(Kn)n→KF : (\mathbb{K}^n)^n \to \mathbb{K}が交代多重線形形式であるとき、成分がK\mathbb{K}に属する任意のnn次正方行列A=(aij)A = (a_{ij})について

F(A)=(∑σ∈Snsgn⁡(σ) a1σ(1)a2σ(2)⋯anσ(n))F(I).F(A) = \left( \sum_{\sigma \in S_n} \operatorname{sgn}(\sigma)\, a_{1\sigma(1)} a_{2\sigma(2)} \cdots a_{n\sigma(n)} \right) F(I).

証明.e⃗1,…,e⃗n\vec e_1, \dots, \vec e_nをKn\mathbb{K}^nの標準基底とすると、AAの第ii行はa⃗i=∑j=1naije⃗j\vec a_i = \sum_{j=1}^n a_{ij} \vec e_jと書くことができる。定義 2.1 条件 (a)を各行について順に用いて展開すると

F(A)=∑j1=1n⋯∑jn=1na1j1a2j2⋯anjn F(e⃗j1;e⃗j2;… ;e⃗jn)F(A) = \sum_{j_1 = 1}^n \cdots \sum_{j_n = 1}^n a_{1j_1} a_{2j_2} \cdots a_{nj_n}\, F(\vec e_{j_1}; \vec e_{j_2}; \dots; \vec e_{j_n})

となる。j1,…,jnj_1, \dots, j_nのうち二つが等しい項では、対応する二つの行が等しいので定義 2.1 条件 (b)によりF(e⃗j1;… ;e⃗jn)=0F(\vec e_{j_1}; \dots; \vec e_{j_n}) = 0である。したがって残るのは、i↦jii \mapsto j_iが全単射である項、すなわちji=σ(i)j_i = \sigma(i)となる置換σ∈Sn\sigma \in S_nに対応する項だけである。

σ\sigmaを命題 3.3 (4)によりrr個の互換の積として表す。行の入れ替えをrr回行うと(e⃗σ(1);… ;e⃗σ(n))(\vec e_{\sigma(1)}; \dots; \vec e_{\sigma(n)})を(e⃗1;… ;e⃗n)=I(\vec e_1; \dots; \vec e_n) = Iへ移すことができ、補題 2.2により1回ごとに符号が変わるので

F(e⃗σ(1);… ;e⃗σ(n))=(−1)rF(I)=sgn⁡(σ)F(I)F(\vec e_{\sigma(1)}; \dots; \vec e_{\sigma(n)}) = (-1)^r F(I) = \operatorname{sgn}(\sigma) F(I)

である。以上を代入すると主張の等式を得る。▨

定理 4.2 (行列式の存在と一意性).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とする。定義 2.1の3条件を満たす関数D:(Kn)n→KD : (\mathbb{K}^n)^n \to \mathbb{K}はただ一つ存在し、成分がK\mathbb{K}に属する任意のnn次正方行列A=(aij)A=(a_{ij})に対して

det⁡A=∑σ∈Snsgn⁡(σ) a1σ(1)a2σ(2)⋯anσ(n)\det A = \sum_{\sigma \in S_n} \operatorname{sgn}(\sigma)\, a_{1\sigma(1)} a_{2\sigma(2)} \cdots a_{n\sigma(n)}(4.2.1)

で与えられる。したがって、交代多重線形形式FFはすべてF(A)=det⁡(A) F(I)F(A) = \det(A) \, F(I)と書くことができる。

証明. 一意性と最後の主張は定理 4.1から従う。実際、3条件を満たすDDについてF=DF = Dとすると、定義 2.1 条件 (c)によりD(I)=1D(I) = 1であるから、定理 4.1の等式は式 (4.2.1)に一致する。

存在を示すため、式 (4.2.1)の右辺をΔ(A)\Delta(A)とおく。

定義 2.1 条件 (a)を示す。各項a1σ(1)⋯anσ(n)a_{1\sigma(1)} \cdots a_{n\sigma(n)}は、各行から成分をちょうど一つずつ取った積である。したがって第kk行の成分akσ(k)a_{k\sigma(k)}について1次であり、他の行の成分を含む因子は第kk行を動かしても変わらない。よってΔ\Deltaは第kk行について線形である。

定義 2.1 条件 (c)を示す。A=IA = Iのときaiσ(i)a_{i\sigma(i)}はσ(i)=i\sigma(i) = iのとき11、そうでなければ00である。積が00でないのはσ=id\sigma = \mathrm{id}の項だけで、その値はsgn⁡(id)=1\operatorname{sgn}(\mathrm{id}) = 1である。よってΔ(I)=1\Delta(I) = 1である。

定義 2.1 条件 (b)を示す。p≠qp \ne qかつ第pp行と第qq行が等しいとする。すなわちすべてのjjについてapj=aqja_{pj} = a_{qj}である。τ=(p  q)\tau = (p\;q)とおくと、σ↦στ\sigma \mapsto \sigma\tauはSnS_nの全単射であり、στ≠σ\sigma\tau \ne \sigmaである。この対応で移り合う二つの項を比べる。στ\sigma\tauはppをσ(q)\sigma(q)へ、qqをσ(p)\sigma(p)へ移し、他のiiではσ\sigmaと一致するので

∏i=1nai στ(i)=ap σ(q) aq σ(p)∏i≠p,qaiσ(i)=aq σ(q) ap σ(p)∏i≠p,qaiσ(i)=∏i=1naiσ(i)\prod_{i=1}^n a_{i\,\sigma\tau(i)} = a_{p\,\sigma(q)}\, a_{q\,\sigma(p)} \prod_{i \ne p, q} a_{i\sigma(i)} = a_{q\,\sigma(q)}\, a_{p\,\sigma(p)} \prod_{i \ne p, q} a_{i\sigma(i)} = \prod_{i=1}^n a_{i\sigma(i)}

となる。2つめの等号でapj=aqja_{pj} = a_{qj}を用いた。一方sgn⁡(στ)=−sgn⁡(σ)\operatorname{sgn}(\sigma\tau) = -\operatorname{sgn}(\sigma)である。したがってSnS_nを{σ,στ}\{\sigma, \sigma\tau\}という対に分けると、各対の寄与は打ち消し合いΔ(A)=0\Delta(A) = 0を得る。▨

式 (4.2.1)を行列式の置換にわたる和による表示といいます。項の個数はn!n!です。この表示は有限和なので値を求めることができますが、次数が上がると計算する項が急速に増えて非効率になります。次数が大きい行列の計算には、次の節以降で導く行基本変形と余因子展開を使います。

例 4.3 (2次と3次の場合).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とし、行列の成分はK\mathbb{K}に属するとする。n=2n = 2のときS2={id,(1  2)}S_2 = \{\mathrm{id}, (1\;2)\}で、符号はそれぞれ11と−1-1である。式 (4.2.1)は

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

となり、冒頭の式に一致する。n=3n = 3のときS3S_3の6個の置換のうち、恒等置換と2つの巡回置換が符号+1+1、3つの互換が符号−1-1である。式 (4.2.1)は

a11a22a33+a12a23a31+a13a21a32−a13a22a31−a11a23a32−a12a21a33a_{11}a_{22}a_{33} + a_{12}a_{23}a_{31} + a_{13}a_{21}a_{32} - a_{13}a_{22}a_{31} - a_{11}a_{23}a_{32} - a_{12}a_{21}a_{33}

となり、これはサラスの方法が並べる6個の積にほかならない。

5 行基本変形と行列式

命題 5.1 (行基本変形と行列式).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とし、AAを成分がK\mathbb{K}に属するnn次正方行列、ccをK\mathbb{K}の元とする。次が成り立つ。

  1. 二つの行を入れ替えると、行列式は符号が変わる。
  2. 一つの行をcc倍すると、行列式はcc倍になる。
  3. ある行に他の行のcc倍を加えても、行列式は変わらない。

証明.(1)は補題 2.2にほかならない。(2)は、定義 2.1 条件 (a)においてy⃗=x⃗\vec y = \vec x、α=c\alpha = c、β=0\beta = 0とした場合である。(3)を示す。第pp行に第qq行(q≠pq \ne p)のcc倍を加えた行列の行列式は、第pp行についての線形性から

det⁡(… ;a⃗p+ca⃗q;… )=det⁡(… ;a⃗p;… )+cdet⁡(… ;a⃗q;… )\det(\dots; \vec a_p + c\vec a_q; \dots) = \det(\dots; \vec a_p; \dots) + c \det(\dots; \vec a_q; \dots)

である。右辺の第2項は、第pp行と第qq行がともにa⃗q\vec a_qである行列の行列式なので、定義 2.1 条件 (b)により00である。▨

命題 5.2 (三角行列の行列式).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とし、AAを成分がK\mathbb{K}に属するnn次正方行列とする。AAが上三角行列(i>ji > jのときaij=0a_{ij} = 0)または下三角行列(i<ji < jのときaij=0a_{ij} = 0)なら、det⁡A=a11a22⋯ann\det A = a_{11}a_{22}\cdots a_{nn}である。

証明.AAが上三角であるとする。式 (4.2.1)の項a1σ(1)⋯anσ(n)a_{1\sigma(1)} \cdots a_{n\sigma(n)}が00でないためには、すべてのiiについてσ(i)≥i\sigma(i) \ge iが必要である。ところが∑i=1nσ(i)=∑i=1ni\sum_{i=1}^n \sigma(i) = \sum_{i=1}^n iなので、すべてのiiでσ(i)≥i\sigma(i) \ge iが成り立つならすべてのiiで等号、すなわちσ=id\sigma = \mathrm{id}である。よって残る項はa11⋯anna_{11}\cdots a_{nn}のみである。下三角の場合も、σ(i)≤i\sigma(i) \le iとして同じ議論による。▨

例 5.3 (行基本変形で値を求める).

A=(242134375)A = \begin{pmatrix} 2 & 4 & 2 \\ 1 & 3 & 4 \\ 3 & 7 & 5 \end{pmatrix}

の行列式を求める。第1行を22で括ると、命題 5.1 (2)によりdet⁡A=2det⁡B\det A = 2 \det Bである。ここでBBは第1行が(1,2,1)(1, 2, 1)で他はAAと同じ行列である。BBの第2行から第1行を引き、第3行から第1行の33倍を引くと、命題 5.1 (3)により行列式は変わらず

det⁡B=det⁡(121013012)\det B = \det\begin{pmatrix} 1 & 2 & 1 \\ 0 & 1 & 3 \\ 0 & 1 & 2 \end{pmatrix}

である。さらに第3行から第2行を引くと

det⁡B=det⁡(12101300−1)=1⋅1⋅(−1)=−1\det B = \det\begin{pmatrix} 1 & 2 & 1 \\ 0 & 1 & 3 \\ 0 & 0 & -1 \end{pmatrix} = 1 \cdot 1 \cdot (-1) = -1

である(最後の等号は命題 5.2)。したがってdet⁡A=2⋅(−1)=−2\det A = 2 \cdot (-1) = -2である。

6 転置と、列についての性質

命題 6.1 (転置の行列式).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とし、AAを成分がK\mathbb{K}に属するnn次正方行列とする。このときdet⁡A⊤=det⁡A\det A^\top = \det A。

証明.A⊤A^\topの(i,j)(i, j)成分はajia_{ji}なので、式 (4.2.1)から

det⁡A⊤=∑σ∈Snsgn⁡(σ) aσ(1)1aσ(2)2⋯aσ(n)n.\det A^\top = \sum_{\sigma \in S_n} \operatorname{sgn}(\sigma)\, a_{\sigma(1)1} a_{\sigma(2)2} \cdots a_{\sigma(n)n}.

各項でk=σ(i)k = \sigma(i)と置き換えるとi=σ−1(k)i = \sigma^{-1}(k)であり、積は因子の順序によらないので

∏i=1naσ(i) i=∏k=1nak σ−1(k)\prod_{i=1}^n a_{\sigma(i)\,i} = \prod_{k=1}^n a_{k\,\sigma^{-1}(k)}

である。また命題 3.3 (2)からsgn⁡(σ)=sgn⁡(σ−1)\operatorname{sgn}(\sigma) = \operatorname{sgn}(\sigma^{-1})であり、σ↦σ−1\sigma \mapsto \sigma^{-1}はSnS_nの全単射である。したがってτ=σ−1\tau = \sigma^{-1}と書き換えると

det⁡A⊤=∑τ∈Snsgn⁡(τ) a1τ(1)⋯anτ(n)=det⁡A.\det A^\top = \sum_{\tau \in S_n} \operatorname{sgn}(\tau)\, a_{1\tau(1)} \cdots a_{n\tau(n)} = \det A.

▨

系 6.2 (列についての線形性と交代性).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とし、AAを成分がK\mathbb{K}に属するnn次正方行列とする。AAの行列式は各列についても線形であり、二つの列が等しければ値は00である。また、二つの列を入れ替えると符号が変わり、ある列に他の列の定数倍を加えても値は変わらない。

証明. 列に関する操作は、転置した行列の行に関する同じ操作である。命題 6.1により行列式の値は転置で変わらないので、定義 2.1と命題 5.1の各主張がそのまま列についても成り立つ。▨

7 余因子展開

定義 7.1 (小行列式・余因子).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とし、AAを成分がK\mathbb{K}に属するnn次正方行列とする。AAから第ii行と第jj列を除いて得られるn−1n-1次正方行列をMijM_{ij}、その行列式det⁡Mij\det M_{ij}を小行列式 (minor) という。

Aij=(−1)i+jdet⁡MijA_{ij} = (-1)^{i+j} \det M_{ij}

をAAの(i,j)(i, j)余因子 (cofactor) という。

定理 7.2 (余因子展開).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とし、AAを成分がK\mathbb{K}に属するnn次正方行列とする。n≥2n \ge 2ならば、各iiについて

det⁡A=∑j=1naijAij\det A = \sum_{j=1}^n a_{ij} A_{ij}

(第ii行についての展開)であり、各jjについて

det⁡A=∑i=1naijAij\det A = \sum_{i=1}^n a_{ij} A_{ij}

(第jj列についての展開)である。

証明. 第ii行について示す。a⃗i=∑jaije⃗j\vec a_i = \sum_j a_{ij}\vec e_jと書き、第ii行についての線形性を使うと

det⁡A=∑j=1naijDij,Dij=det⁡(a⃗1;… ;a⃗i−1;e⃗j;a⃗i+1;… ;a⃗n)\det A = \sum_{j=1}^n a_{ij} D_{ij}, \qquad D_{ij} = \det(\vec a_1; \dots; \vec a_{i-1}; \vec e_j; \vec a_{i+1}; \dots; \vec a_n)

となる。Dij=(−1)i+jdet⁡MijD_{ij} = (-1)^{i+j}\det M_{ij}を示せばよい。

DijD_{ij}の行列で、第ii行を隣接する行とi−1i-1回入れ替えて第1行へ移し、続いて第jj列を隣接する列とj−1j-1回入れ替えて第1列へ移す。行の入れ替えは補題 2.2、列の入れ替えは系 6.2により、1回ごとに符号が変わる。これらの入れ替えは、除かれる行と列以外の行・列の相対的な順序を変えないので、得られる行列は

C=(10⃗⊤c⃗Mij)C = \begin{pmatrix} 1 & \vec 0^\top \\ \vec c & M_{ij} \end{pmatrix}

の形であり

Dij=(−1)(i−1)+(j−1)det⁡C=(−1)i+jdet⁡CD_{ij} = (-1)^{(i-1)+(j-1)} \det C = (-1)^{i+j}\det C

となる。ここでc⃗∈Kn−1\vec c \in \mathbb{K}^{n-1}は第1列の残りの成分である。CCの第1行は(1,0,…,0)(1, 0, \dots, 0)なので、k≥2k \ge 2について第kk行から第1行のck−1c_{k-1}倍を引くと、命題 5.1 (3)により行列式は変わらず、第1列はc⃗=0⃗\vec c = \vec 0となる。すなわち

det⁡C=det⁡(10⃗⊤0⃗Mij).\det C = \det\begin{pmatrix} 1 & \vec 0^\top \\ \vec 0 & M_{ij} \end{pmatrix}.

最後に、n−1n-1次正方行列XXに対して

G(X)=det⁡(10⃗⊤0⃗X)G(X) = \det\begin{pmatrix} 1 & \vec 0^\top \\ \vec 0 & X \end{pmatrix}

と定めると、GGはXXの各行について線形であり、XXの二つの行が等しければ値は00である。さらにG(In−1)=det⁡In=1G(I_{n-1}) = \det I_n = 1である。よってGGはn−1n-1次の行列式の3条件を満たすので、定理 4.2の一意性からG(X)=det⁡XG(X) = \det Xである。したがってdet⁡C=det⁡Mij\det C = \det M_{ij}、すなわちDij=(−1)i+jdet⁡MijD_{ij} = (-1)^{i+j}\det M_{ij}を得る。

列についての展開は、命題 6.1を用いてA⊤A^\topの行についての展開へ帰着する。A⊤A^\topの(j,i)(j, i)余因子はAAの(i,j)(i, j)余因子に等しいので、主張の形になる。▨

例 7.3 (余因子展開による計算). 行列

A=(102030405)A=\begin{pmatrix}1&0&2\\0&3&0\\4&0&5\end{pmatrix}

の第22行には00が二つある。第22行について定理 7.2を用いると

det⁡A=0⋅A21+3⋅A22+0⋅A23=3(−1)2+2det⁡(1245)=3(5−8)=−9\det A =0\cdot A_{21}+3\cdot A_{22}+0\cdot A_{23} =3(-1)^{2+2}\det\begin{pmatrix}1&2\\4&5\end{pmatrix} =3(5-8)=-9

である。

aij=0a_{ij} = 0である項は和に寄与しないので、00の多い行または列について展開すると、値を求めるために計算する小行列式の個数が減ります。

定理 7.4 (余因子行列).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とし、AAを成分がK\mathbb{K}に属するnn次正方行列とする。n≥2n \ge 2とし、(i,j)(i, j)成分がAjiA_{ji}であるnn次正方行列をA~\widetilde{A}と書き、AAの余因子行列という。このとき

AA~=A~A=(det⁡A)I.A\widetilde{A} = \widetilde{A}A = (\det A) I.

証明.AA~A\widetilde{A}の(i,k)(i, k)成分は∑jaijAkj\sum_{j} a_{ij} A_{kj}である。i=ki = kのときは定理 7.2の第ii行についての展開そのもので、値はdet⁡A\det Aである。

i≠ki \ne kとする。AAの第kk行を第ii行で置き換えた行列をA′A'とおく。A′A'の第kk行以外の行はAAと同じであり、(k,j)(k, j)余因子は第kk行を除いて作るので、A′A'の(k,j)(k, j)余因子はAAの(k,j)(k, j)余因子AkjA_{kj}に等しい。A′A'を第kk行について展開すると

det⁡A′=∑jakj′Akj=∑jaijAkj\det A' = \sum_j a'_{kj} A_{kj} = \sum_j a_{ij} A_{kj}

である。一方A′A'は第ii行と第kk行がともにa⃗i\vec a_iなので、定義 2.1 条件 (b)によりdet⁡A′=0\det A' = 0である。したがってi≠ki \ne kのとき∑jaijAkj=0\sum_j a_{ij}A_{kj} = 0であり、AA~=(det⁡A)IA\widetilde{A} = (\det A) Iを得る。

A~A\widetilde{A}Aの(i,k)(i, k)成分は∑jA~ijajk=∑jajkAji\sum_j \widetilde{A}_{ij} a_{jk} = \sum_j a_{jk} A_{ji}である。i=ki = kのときは定理 7.2の第ii列についての展開そのもので、値はdet⁡A\det Aである。i≠ki \ne kのときは、AAの第ii列を第kk列で置き換えた行列を第ii列について展開した式であり(第ii列を除いて作る余因子は置き換えで変わらない)、その行列は第ii列と第kk列が等しいので、系 6.2により行列式は00である。よってA~A=(det⁡A)I\widetilde{A}A = (\det A) Iも成り立つ。▨

例 7.5 (2次の余因子行列と逆行列).a,b,c,da,b,c,dを実数とし、A=(abcd)A=\begin{pmatrix}a&b\\c&d\end{pmatrix}とする。AAの余因子は

A11=d,A12=−c,A21=−b,A22=aA_{11}=d,\qquad A_{12}=-c,\qquad A_{21}=-b,\qquad A_{22}=a

であるから

A~=(d−b−ca)\widetilde A=\begin{pmatrix}d&-b\\-c&a\end{pmatrix}

である。det⁡A=ad−bc≠0\det A=ad-bc\ne0ならば、定理 7.4の両辺をdet⁡A\det Aで割ることにより

A−1=1det⁡AA~=1ad−bc(d−b−ca)A^{-1}=\frac{1}{\det A}\widetilde A =\frac{1}{ad-bc}\begin{pmatrix}d&-b\\-c&a\end{pmatrix}

を得る。この式は§D3.3 公式 4.1と一致する。

8 基本行列

定義 8.1 (基本行列).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とする。成分がK\mathbb{K}に属する単位行列IIに1回だけ行基本変形を施して得られる行列を基本行列 (elementary matrix) という。次の3種類がある。

  1. 第pp行と第qq行を入れ替えたもの(p≠qp \ne q)をPpqP_{pq}と書く。
  2. 第pp行をcc倍したもの(c∈Kc \in \mathbb{K}、c≠0c \ne 0)をDp(c)D_p(c)と書く。
  3. 第pp行に第qq行のcc倍を加えたもの(c∈Kc \in \mathbb{K}、p≠qp \ne q)をNpq(c)N_{pq}(c)と書く。

命題 8.2 (基本行列の作用と行列式).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とし、EEをK\mathbb{K}上の基本行列、AAを成分がK\mathbb{K}に属するnn次正方行列とする。このときEAEAはAAにEEの作られ方と同じ行基本変形を施した行列である。また

det⁡Ppq=−1,det⁡Dp(c)=c,det⁡Npq(c)=1.\det P_{pq} = -1, \qquad \det D_p(c) = c, \qquad \det N_{pq}(c) = 1.

証明.EAEAの(i,j)(i,j)成分は行列の積の定義により∑k=1neikakj\sum_{k=1}^n e_{ik}a_{kj}であるから、EAEAの第ii行は、EEの第ii行の成分ei1,…,eine_{i1}, \dots, e_{in}を係数とするAAの第11行から第nn行までの一次結合である。3種類それぞれについてEEの行を調べる。

  • PpqP_{pq}の第pp行はe⃗q\vec e_q、第qq行はe⃗p\vec e_p、それ以外の第ii行はe⃗i\vec e_iである。よってPpqAP_{pq}AはAAの第pp行と第qq行を入れ替えた行列である。
  • Dp(c)D_p(c)の第pp行はc e⃗pc\,\vec e_p、それ以外の第ii行はe⃗i\vec e_iである。よってDp(c)AD_p(c)AはAAの第pp行をcc倍した行列である。
  • Npq(c)N_{pq}(c)の第pp行はe⃗p+c e⃗q\vec e_p + c\,\vec e_q、それ以外の第ii行はe⃗i\vec e_iである。よってNpq(c)AN_{pq}(c)AはAAの第pp行に第qq行のcc倍を加えた行列である。

行列式の値は、E=EIE = EIとみて上の作用を単位行列へ適用し、det⁡I=1\det I = 1と命題 5.1 (1)、命題 5.1 (2)、命題 5.1 (3)をそれぞれ用いれば得られる。▨

9 積の行列式

定理 9.1 (積の行列式).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とし、A,BA,Bを成分がK\mathbb{K}に属するnn次正方行列とする。このときdet⁡(AB)=det⁡A⋅det⁡B\det(AB) = \det A \cdot \det Bである。

証明.BBを固定し、F(A)=det⁡(AB)F(A) = \det(AB)と定める。行列の積の定義により、ABABの第ii行はAAの第ii行a⃗i\vec a_iを用いてa⃗iB\vec a_i Bと書くことができる。写像x⃗↦x⃗B\vec x \mapsto \vec x Bは線形なので、FFはAAの各行について線形である。またAAの第pp行と第qq行が等しければ、ABABの第pp行と第qq行も等しいので、定義 2.1 条件 (b)によりF(A)=det⁡(AB)=0F(A) = \det(AB) = 0である。よってFFは交代多重線形形式であり、定理 4.2から

F(A)=det⁡(A) F(I)=det⁡(A)det⁡(IB)=det⁡A⋅det⁡BF(A) = \det(A)\, F(I) = \det(A) \det(IB) = \det A \cdot \det B

を得る。▨

10 行列式と正則性

定理 10.1 (行列式が00でないことと正則であること).K\mathbb{K}をR\mathbb{R}またはC\mathbb{C}とし、AAを成分がK\mathbb{K}に属するnn次正方行列とする。次の3つは同値である。

  1. det⁡A≠0\det A \ne 0である。
  2. AAは正則である。すなわち、AB=BA=IAB = BA = Iを満たす、成分がK\mathbb{K}に属するnn次正方行列BBが存在する。
  3. Ax⃗=0⃗A\vec x = \vec 0を満たすx⃗∈Kn\vec x \in \mathbb{K}^nはx⃗=0⃗\vec x = \vec 0のみである。

さらにAAが正則ならdet⁡(A−1)=(det⁡A)−1\det(A^{-1}) = (\det A)^{-1}であり、n≥2n \ge 2のときはA−1=(det⁡A)−1A~A^{-1} = (\det A)^{-1}\widetilde{A}と余因子行列によって書き表される。K=R\mathbb{K} = \mathbb{R}のときは、上の3つは§D3.3 定理 2.2の各条件とも同値である。

証明.(1)⇒\Rightarrow(2)を示す。n=1n = 1のときA=(a)A = (a)と書くとdet⁡A=a≠0\det A = a \ne 0であり、H=(a−1)H = (a^{-1})がAH=HA=IAH = HA = Iを満たす。n≥2n \ge 2のとき、定理 7.4のAA~=A~A=(det⁡A)IA\widetilde{A} = \widetilde{A}A = (\det A)Iの両辺をdet⁡A≠0\det A \ne 0で割ると、H=(det⁡A)−1A~H = (\det A)^{-1}\widetilde{A}がAH=HA=IAH = HA = Iを満たす。

(2)⇒\Rightarrow(1)を示す。AH=HA=IAH=HA=Iを満たすHHを取る。定理 9.1により

1=det⁡I=det⁡(AH)=det⁡Adet⁡H1=\det I=\det(AH)=\det A\det H

であるから、det⁡A≠0\det A\ne0である。

(2)⇒\Rightarrow(3)を示す。AH=HA=IAH = HA = Iを満たすHHを取る。Ax⃗=0⃗A\vec x = \vec 0ならばx⃗=Ix⃗=H(Ax⃗)=H0⃗=0⃗\vec x = I\vec x = H(A\vec x) = H\vec 0 = \vec 0である。

(3)⇒\Rightarrow(1)を対偶で示す。det⁡A=0\det A=0と仮定し、nnについての帰納法によりAx⃗=0⃗A\vec x=\vec0を満たす0⃗\vec0でないx⃗∈Kn\vec x\in\mathbb{K}^nが存在することを示す。n=1n=1ならばA=(0)A=(0)であるから、x⃗=(1)\vec x=(1)と取ることができる。

n≥2n\ge2とし、n−1n-1次正方行列について主張が成り立つと仮定する。AAの第1列が零ベクトルならばx⃗=e⃗1\vec x=\vec e_1と取ることができる。第1列が零ベクトルでないならば、ap1≠0a_{p1}\ne0を満たすppを取る。p≠1p\ne1ならば第1行と第pp行を入れ替え、得られた第1列の第1成分をα\alphaとして第1行をα−1\alpha^{-1}倍する。続いて、得られた第1列の第ii成分をuiu_iとして、各i=2,…,ni=2,\ldots,nについて第ii行へ第1行の−ui-u_i倍を加える。

これらの行基本変形に対応する基本行列を、操作した順にE1,…,EmE_1,\ldots,E_mとし、E=Em⋯E1E=E_m\cdots E_1とおくと、EAEAの第1列はe⃗1\vec e_1である。各EiE_iに対応する行基本変形には逆向きの行基本変形があり、その基本行列をEi′E_i'とすればEi′Ei=EiEi′=IE_i'E_i=E_iE_i'=Iである。したがってE1′⋯Em′E_1'\cdots E_m'はEEの両側逆行列であり、EEは正則である。

EA=(1r⃗⊤0⃗B)EA=\begin{pmatrix}1&\vec r^{\top}\\\vec0&B\end{pmatrix}

と書くことができる。ここでBBは成分がK\mathbb{K}に属するn−1n-1次正方行列である。定理 9.1によりdet⁡(EA)=det⁡Edet⁡A=0\det(EA)=\det E\det A=0であり、定理 7.2をEAEAの第1列へ適用するとdet⁡(EA)=det⁡B\det(EA)=\det Bである。したがってdet⁡B=0\det B=0である。

帰納法の仮定により、By⃗=0⃗B\vec y=\vec0を満たす0⃗\vec0でないy⃗∈Kn−1\vec y\in\mathbb{K}^{n-1}が存在する。第1成分が−r⃗⊤y⃗-\vec r^{\top}\vec yであり、残りの成分を並べたベクトルがy⃗\vec yであるx⃗∈Kn\vec x\in\mathbb{K}^nを取ると、x⃗≠0⃗\vec x\ne\vec0かつ(EA)x⃗=0⃗(EA)\vec x=\vec0である。EEは正則なので、左からE−1E^{-1}を掛けるとAx⃗=0⃗A\vec x=\vec0を得る。よって対偶が成り立つ。

AAが正則であるとし、AH=HA=IAH = HA = Iを満たすHHを取る。AC=CA=IAC=CA=Iを満たす行列CCがあればH=H(AC)=(HA)C=CH=H(AC)=(HA)C=Cであるから、このようなHHは一意である。これをA−1A^{-1}と書く。定理 9.1によりdet⁡A⋅det⁡(A−1)=det⁡I=1\det A \cdot \det(A^{-1}) = \det I = 1であり、とくにdet⁡A≠0\det A \ne 0であってdet⁡(A−1)=(det⁡A)−1\det(A^{-1}) = (\det A)^{-1}である。n≥2n \ge 2のときは、(1)⇒\Rightarrow(2)の証明で構成した行列(det⁡A)−1A~(\det A)^{-1}\widetilde{A}がAAとの積を両側でIIにするので、一意性によりA−1=(det⁡A)−1A~A^{-1} = (\det A)^{-1}\widetilde{A}である。

K=R\mathbb{K} = \mathbb{R}とする。(2)は§D3.3 定理 2.2 (1)と同じ主張であり、(3)は§D3.3 定理 2.2 (2)と同じ主張である。よって§D3.3 定理 2.2の各条件は、上の3つと同値である。▨

定理 10.2 (正則行列は基本行列の積である).AAを実数を成分とするnn次正則行列とすると、有限個の実数を成分とする基本行列E1,…,EmE_1, \dots, E_mが存在して

A=E1E2⋯EmA = E_1 E_2 \cdots E_m

と書くことができる。

証明.§D3.2 定理 3.2により、AAに有限回の行基本変形を続けて簡約階段形RRへ変形することができる。命題 8.2により、各段の行基本変形は基本行列を左から掛けることに等しいので、基本行列F1,…,FmF_1, \dots, F_mを用いて

Fm⋯F1A=RF_m \cdots F_1 A = R

と書くことができる。

AAは正則であるから、定理 10.1 (2)⇒\Rightarrow(3)によりAx⃗=0⃗A\vec x = \vec 0の解はx⃗=0⃗\vec x = \vec 0のみである。RRはAAに有限回の行基本変形を続けて得られる簡約階段形であるから、§D3.2 系 4.2 (1)⇒\Rightarrow(3)によりR=IR = Iであり

Fm⋯F1A=IF_m \cdots F_1 A = I

が成り立つ。

基本行列は正則であり、その逆行列もまた基本行列である。実際、命題 8.2により基本行列を左から掛けることは対応する行基本変形を施すことであり、第pp行と第qq行の入れ替えを二度続けると元の行列に戻り、第pp行をcc倍したのちc−1c^{-1}倍すると元の行列に戻り、第pp行に第qq行のcc倍を加えたのち−c-c倍を加えると元の行列に戻る。よってEEを基本行列、E′E'を逆向きの操作に対応する基本行列とするとE′E=E′(EI)=IE'E = E'(EI) = IかつEE′=E(E′I)=IEE' = E(E'I) = Iであり、Ppq−1=PpqP_{pq}^{-1} = P_{pq}、Dp(c)−1=Dp(c−1)D_p(c)^{-1} = D_p(c^{-1})(c≠0c \ne 0)、Npq(c)−1=Npq(−c)N_{pq}(c)^{-1} = N_{pq}(-c)であって、いずれも定義 8.1 (1)、定義 8.1 (2)、定義 8.1 (3)のいずれかの形である。したがってFm⋯F1A=IF_m \cdots F_1 A = Iの両辺へ左からF1−1F2−1⋯Fm−1F_1^{-1} F_2^{-1} \cdots F_m^{-1}を掛けると

A=F1−1F2−1⋯Fm−1A = F_1^{-1} F_2^{-1} \cdots F_m^{-1}

を得る。右辺は基本行列の積である。▨

幾何的には、∣det⁡A∣|\det A|はAAによる面積(体積)の拡大率で、det⁡A=0\det A = 0は「面積が00につぶれる(次元が落ちる)」ことを意味します。次項でこの見方を図で確認します。

参考文献

  1. 齋藤正彦『線型代数入門』基礎数学 1, 東京大学出版会, 1966.置換の符号から行列式を定義し、余因子展開と正則性の同値へ進む道筋を参考にしました。
  2. Kenneth Hoffman and Ray Kunze, Linear Algebra, 2nd ed., Prentice Hall, Englewood Cliffs, N.J., 1971.各行についての線形性と交代性によって行列式を特徴づけ、存在と一意性を別々に示す構成を参考にしました。
  3. Serge Lang, Linear Algebra, 3rd ed., Springer, New York, 1987.積の行列式が各因子の行列式の積になることを、一意性から導く議論を参考にしました。

前提記事