クラメルの公式
連立一次方程式を行列(係数行列)を用いて表す方法を導入します。また、行列式を用いた連立一次方程式の解法であるクラメルの公式を示します。
クラメルの公式は、係数行列が正則であるような、特定の形の連立一次方程式においてのみ使える解の公式です。
係数行列の導入
まず、連立一次方程式を行列(係数行列)を用いて表す方法を導入します。
連立一次方程式の行列による表現
個の変数 についての 個の式からなる、次のような連立一次方程式について考えます。
いま、次のように行列 と列ベクトル を次のように置くと、
連立一次方程式(3.6.9)式は、次のように表すことができます。
係数行列とは
ここで、 は 行列、 は 項列ベクトル、 は 項列ベクトルです。 を特に係数行列()といいます。
係数行列を用いることで、一般の連立方一次方程式を のように簡単に表すことができます。このような表し方は汎用的であり、基本的にどのような連立一次方程式も(3.6.10)式のように表すことができます。
また、連立一次方程式を係数行列で表す方法は便利であるため、基底と次元の準備や、行列の基本変形による連立一次方程式の解法、固有値と固有ベクトルなど、線形代数学の様々な領域において用いられます。
クラメルの公式
次に、行列式を用いた連立一次方程式の解法であるクラメルの公式を示します。
定理 3.23(クラメルの公式)
を 次の正方行列とする。 が正則であるとき、 を係数行列として持つ連立一次方程式 の解はただ つに定まり、次の式により与えられる。
ここで、 であり、 の分母は係数行列の行列式、分子は係数行列の第 列を に置き換えた行列の行列式である。
解説
行列式による連立一次方程式の解法
定理 3.23(クラメルの公式)は、行列式による連立一次方程式の解法を与えるものであり、クラメルの公式()と呼ばれます。
すなわち、係数行列 が正則であるような連立一次方程式 の解 の成分 は、(3.6.11)式により与えられます。そのような意味で、クラメルの公式は、特定の形の連立一次方程式ついてのみ成り立つ、解の公式であるといえます。
クラメルの公式が使える条件
クラメルの公式が適用できるのは、「係数行列 が正則である(または、 である)」 場合に限られます。
- 係数行列 が正則である
- である
より具体的には、定理 3.23(クラメルの公式)は「連立一次方程式の変数の個数 と式の個数 が等しく、かつ係数行列 が正則である」ときのみ成り立つといえます。
変数の個数 と式の個数 が等しければ、係数行列は正方行列となります。また、係数行列が正方行列であることは、係数行列が正則であることの必要条件です(正則行列の定義)。したがって、クラメルの公式が成り立つ条件は、単に「係数行列が正則である」ことと表すことができます。
また、「係数行列が正則である」という条件は、前項の定理 3.22(逆行列を持つための条件)より、「 である」ことと同値です。
連立一次方程式の解法としての実用性
具体的に与えられた連立一次方程式を解くにあたって、クラメルの公式を用いることは必ずしも効率的ではありません。
クラメルの公式が実用的でない理由
(3.6.11)式からわかるとおり、解を得るためには 次行列式を 個(分母となる係数行列の行列式 個と、分子の行列式 個)計算する必要があります。これは非常に面倒で手間のかかる作業です。
行列の基本変形による解法(より一般的な解法)
連立一次方程式の解法としては、一般的には、行列の基本変形による方法が実用的です。基本変形による連立一次方程式の解法の手順は、次の通りです。
クラメルの公式が有効な場合
クラメルの公式が有効なのは、下記の例題 2 のような場合です。すなわち、係数が文字であり、かつ規則的に並んでいるような連立一次方程式の場合などです。
証明
係数行列 は正則であるから、 の逆行列 が存在する。逆行列は一意的であるから、 の解は により与えられ、ただ つに定まる。また、定理 3.22(逆行列を持つための条件)より、 が正則であるとき、 であるから、次が成り立つ。
ここで、右辺の列ベクトルの第 成分 は、係数行列 の第 列を で置き換えた行列を として、 の行列式を第 列に関して展開したものに他ならない。すなわち、次が成り立つ。
したがって、次が成り立つ。
証明の考え方
解の唯一性の確認
- 仮定より、係数行列
は正則であり、逆行列A A を持ちます。A − 1 A^{-1} - 定理 2.4(逆行列の一意性)より、
に対して逆行列A A は一意に定まります。A − 1 A^{-1}
- 定理 2.4(逆行列の一意性)より、
に左からA x = b A \bm{x} = \bm{b} を掛けることで、A − 1 A^{-1} を得ます。x = A − 1 b \bm{x} = A^{-1} \bm{b} - 当然ながら、これは連立一次方程式
の解となります。A x = b A \bm{x} = \bm{b}
- 当然ながら、これは連立一次方程式
- また、逆行列の一意性から、
の解は、A x = b A \bm{x} = \bm{b} のみであることがわかります。x = A − 1 b \bm{x} = A^{-1} \bm{b}
行列式による解の導出
の解である、A x = b A \bm{x} = \bm{b} をより計算しやすい形に変形することを考えます。x = A − 1 b \bm{x} = A^{-1} \bm{b} いま、
が正則であることが分かっていますので、定理 3.22(逆行列を持つための条件)より、次が成り立ちます。A A x = ( i ) A − 1 b = ( ii ) 1 det A A ~ b = ( iii ) 1 det A ( a ~ 11 a ~ 21 ⋯ a ~ n 1 a ~ 12 a ~ 22 ⋯ a ~ n 2 ⋮ ⋮ ⋱ ⋮ a ~ 1 n a ~ 2 n ⋯ a ~ n n ) ( b 1 b 2 ⋮ b n ) = ( iv ) 1 det A ( b 1 a ~ 11 + b 2 a ~ 21 + ⋯ + b n a ~ n 1 b 1 a ~ 12 + b 2 a ~ 22 + ⋯ + b n a ~ n 2 ⋮ b 1 a ~ 1 n + b 2 a ~ 2 n + ⋯ + b n a ~ n n ) = ( v ) 1 det A ( ∑ i b i a ~ i 1 ∑ i b i a ~ i 2 ⋮ ∑ i b i a ~ i n ) \begin{split} \bm{x} &\overset{(\text{i})}{=} A^{-1} \bm{b} \\ &\overset{(\text{ii})}{=} \frac{\, 1 \,}{\, \det A \,} \tilde{A} \, \bm{b} \\ &\overset{(\text{iii})}{=} \frac{\, 1 \,}{\, \det A \,} \begin{pmatrix} \, \tilde{a}_{11} & \tilde{a}_{21} & \cdots & \tilde{a}_{n1} \, \\ \, \tilde{a}_{12} & \tilde{a}_{22} & \cdots & \tilde{a}_{n2} \, \\ \, \vdots & \vdots & \ddots & \vdots \, \\ \, \tilde{a}_{1n} & \tilde{a}_{2n} & \cdots & \tilde{a}_{nn} \, \\ \end{pmatrix} \begin{pmatrix} \, b_{1} \\ b_{2} \\ \vdots \\ b_{n} \, \end{pmatrix} \\ &\overset{(\text{iv})}{=} \frac{\, 1 \,}{\, \det A \,} \begin{pmatrix} \, b_1 \, \tilde{a}_{11} + b_2 \, \tilde{a}_{21} + \cdots + b_n \, \tilde{a}_{n1} \, \\ \, b_1 \, \tilde{a}_{12} + b_2 \, \tilde{a}_{22} + \cdots + b_n \, \tilde{a}_{n2} \, \\ \, \vdots \, \\ \, b_1 \, \tilde{a}_{1n} + b_2 \, \tilde{a}_{2n} + \cdots + b_n \, \tilde{a}_{nn} \, \\ \end{pmatrix} \\ &\overset{(\text{v})}{=} \frac{\, 1 \,}{\, \det A \,} \begin{pmatrix} \, \displaystyle \sum_{i} \; b_i \, \tilde{a}_{i1} \, \\ \, \displaystyle \sum_{i} \; b_i \, \tilde{a}_{i2} \, \\ \, \vdots \, \\ \, \displaystyle \sum_{i} \; b_i \, \tilde{a}_{in} \, \\ \end{pmatrix} \\ \end{split} (
)定理 3.22(逆行列を持つための条件)より、ii \text{ii} が成り立ちます。A − 1 = 1 det A A ~ A^{-1} = \dfrac{1}{\det A} \tilde{A} (
)において、右辺の列ベクトルの第v \text{v} 成分(第j j 行)は次の通りです。j j b 1 a ~ 1 j + b 2 a ~ 2 j + ⋯ + b n a ~ n j = ∑ i b i a ~ i j \begin{align*} b_1 \, \tilde{a}_{1j} + b_2 \, \tilde{a}_{2j} + \cdots + b_n \, \tilde{a}_{nj} = \sum_{i} \; b_i \, \tilde{a}_{ij} \end{align*} これは、定理 3.19(行列式の展開 1)における列に関する展開と同じ形となっています。
a 1 j a ~ 1 j + a 2 j a ~ 2 j + ⋯ + a n j a ~ n j = ∑ i a i j a ~ i j \begin{align*} a_{1j} \, \tilde{a}_{1j} + a_{2j} \, \tilde{a}_{2j} + \cdots + a_{nj} \, \tilde{a}_{nj} = \sum_{i} \; a_{ij} \, \tilde{a}_{ij} \end{align*} つの式を比較してみると、2 2 の第A A 列の成分j j が、a 1 j , a 2 j , ⋯ , a n j a_{1j}, a_{2j}, \cdots, a_{nj} の列ベクトルの成分b \bm{b} に置き換わっていることが分かります。b 1 , b 2 , ⋯ , b n b_1, b_2, \cdots, b_n すなわち、
は、係数行列∑ i b i a ~ i j \displaystyle \sum_{i} \; b_i \, \tilde{a}_{ij} の第A A 列をj j で置き換えた行列をb \bm{b} として、A j A_j の行列式を第A j A_j 列に関して展開したものと等しいということがわかります。j j ∑ i b i a ~ i j = ∣ a 11 ⋯ b 1 ⋯ a 1 n a 21 ⋯ b 2 ⋯ a 2 n ⋮ ⋮ ⋮ a n 1 ⋯ b n ⋯ a n n ∣ = det A j \begin{align*} \sum_{i} \; b_i \, \tilde{a}_{ij} &= \begin{vmatrix} \; a_{11} & \cdots & b_{1} & \cdots & a_{1n} \; \\ \; a_{21} & \cdots & b_{2} & \cdots & a_{2n} \; \\ \; \vdots & & \vdots & & \vdots \; \\ \; a_{n1} & \cdots & b_{n} & \cdots & a_{nn} \; \\ \end{vmatrix} \\ \\ &= \det A_j \end{align*}
したがって、
は次のようになります。これはすべてのx j x_j について成り立ちます。j ( = 1 , 2 , ⋯ , n ) j \; (= 1, 2, \cdots, n) x j = ∣ a 11 ⋯ b 1 ⋯ a 1 n ⋮ ⋮ ⋮ a n 1 ⋯ b n ⋯ a n n ∣ ∣ a 11 ⋯ a 1 j ⋯ a 1 n ⋮ ⋮ ⋮ a n 1 ⋯ a n j ⋯ a n n ∣ \begin{align*} x_j = \frac{\; \begin{vmatrix} \; a_{11} \; \cdots \; \, b_{1} \, \; \cdots \; a_{1n} \; \\ \; \vdots \quad \qquad \vdots \qquad \quad \vdots \; \; \\ \; a_{n1} \; \cdots \; \, b_{n} \, \; \cdots \; a_{nn} \; \\ \end{vmatrix} \; }{ \; \begin{vmatrix} \; a_{11} \; \cdots \; a_{1j} \; \cdots \; a_{1n} \; \\ \; \vdots \quad \qquad \vdots \qquad \quad \vdots \; \; \\ \; a_{n1} \; \cdots \; a_{nj} \; \cdots \; a_{nn} \; \\ \end{vmatrix} \; } \end{align*} 以上から(3.6.11)式が成り立つことが示されました。
クラメルの公式による解法の例
クラメルの公式による解法の例として、次の
例題 1(変数と式の数が 3 3 の場合)
次の連立一次方程式の解を求めよ。
解答(例題 1)
いま、
解答の考え方(例題 1)
クラメルの公式が使えること(係数行列が正則であること)を確認した上で、公式を適用します。
- 例題 1 では、係数行列の行列式(分母)と
の分子の行列式を求めるため、x , y , z x, y, z 次の行列式を合わせて3 3 つ計算しなければなりません。これは、かなり面倒です。4 4 - この場合、クラメルの公式を適用するよりも、行列の基本変形による解法の方が適していると考えられます(連立一次方程式の解法としての実用性を参照)。
例題 2(係数行列が規則的である場合)
次の連立一次方程式の解を求めよ。ただし、
解答(例題 2)
また、このとき、
いま、
解答の考え方(例題 2)
クラメルの公式が使えること(係数行列が正則であること)を確認した上で、公式を適用します。
- 例題 2 では、係数行列は規則的な配列であり、その行列式は
と、きれいな形になります。( a − b ) ( b − c ) ( c − a ) (a-b) (b-c) (c-a) - クラメルの公式を適用する際、例えば
を求めるとき、x x を含む第a a 列を1 1 を含む列に置き換えればよいため、分子の行列式の値はd d をa a に置き換えたd d になることが直ちにわかります。( d − b ) ( b − c ) ( c − d ) (d-b) (b-c) (c-d) についても同様です。y , z y, z - このように、係数が文字で規則的に並んでいるような連立一次方程式の場合、クラメルの公式を用いることで効率的に解を得ることができます。
まとめ
クラメルの公式:係数行列
が正則であるとき、連立一次方程式A A の解はただA x = b A \bm{x} = \bm{b} つに定まり、次の式により与えられる。1 1 x j = ∣ a 11 ⋯ b 1 ⋯ a 1 n ⋮ ⋮ ⋮ a n 1 ⋯ b n ⋯ a n n ∣ ∣ a 11 ⋯ a 1 j ⋯ a 1 n ⋮ ⋮ ⋮ a n 1 ⋯ a n j ⋯ a n n ∣ ( j = 1 , ⋯ , n ) \begin{equation*} \begin{array} {ccc} x_{j} = \frac{ \; \begin{vmatrix} \; a_{11} \; \cdots \; \, b_{1} \, \; \cdots \; a_{1n} \; \\ \; \vdots \quad \qquad \vdots \qquad \quad \vdots \; \; \\ \; a_{n1} \; \cdots \; \, b_{n} \, \; \cdots \; a_{nn} \; \\ \end{vmatrix} \; }{ \; \begin{vmatrix} \; a_{11} \; \cdots \; a_{1j} \; \cdots \; a_{1n} \; \\ \; \vdots \quad \qquad \vdots \qquad \quad \vdots \; \; \\ \; a_{n1} \; \cdots \; a_{nj} \; \cdots \; a_{nn} \; \\ \end{vmatrix} \; } & (\; j = 1, \cdots, n\; ) \end{array} \end{equation*} - クラメルの公式は、係数行列が正則である場合にのみ成り立つ。
- 文字で規則的に並んでいる場合などに有効であるが、一般に効率的な解法とはいえない。
参考文献
[1] 齋藤正彦. 線型代数入門. 東京大学出版会. 1966.[2] 永田雅宣 他. 理系のための線型代数の基礎. 紀伊國屋書店. 1986.
[3] 川久保勝夫. 線形代数学 [新装版]. 日本評論社. 2010.
[4] 松坂和夫. 線型代数入門 [新装版]. 岩波書店. 2018.
[5] S. Lang. Linear Algebra Third Edition. Springer. 1987.
[6] 雪江明彦. 代数学
[7] 雪江明彦. 代数学
[8] 桂利行. 代数学
[9] 松坂和夫. 代数系入門. 岩波書店. 1976.
[10] 高木貞治. 代数学講義 [改訂新版]. 共立出版. 1965.
[11] S. Lang. Algebra Revised Third Edition. Springer. 2005.
[12] M. Artin. Algebra Second Edition. Pearson Education Limited. 2014.
[13] 青本和彦 他. 数学入門辞典. 岩波書店. 2005.