一般の連立一次方程式の解法(2)

一般の連立一次方程式の解法を示します。

一般の連立一次方程式(Ax=bA \bm{x} = \bm{b})は、拡大係数行列を行標準形に変形することで解を得ることができます。

一般の連立一次方程式の解法

概要

斉次連立一次方程式との違い(解の有無判定)

斉次連立一次方程式 Ax=0A \bm{x} = \bm{0} の場合と異なり、一般の連立一次方程式 Ax=bA \bm{x} = \bm{b} は必ずしも解を持つとは限りません。

このため、一般の連立一次方程式 Ax=bA \bm{x} = \bm{b} を解く際は、まず、Ax=bA \bm{x} = \bm{b} が解を持つか否かを判定する必要があります。解の有無を判定する方法は前項の通りです。

斉次連立一次方程式との違い(解法の手順)

解法の手順も、斉次連立一次方程式の解法と一般の連立一次方程式の解法は若干異なります。

連立一次方程式が斉次(Ax=0A \bm{x} = \bm{0})の場合、係数行列 AA を行標準形に変形することで解が得られました(斉次連立一次方程式の解法)。一方で、一般の連立一次方程式(Ax=bA \bm{x} = \bm{b})の場合、係数拡大行列 (Ab)(\, A \mid \bm{b} \,) を行標準形に変形することで解が得られます。

クラメルの公式による解法との違い

我々は、一般の連立一次方程式の解法として、クラメルの公式を既に得ています。当然ながら、クラメルの公式は斉次連立一次方程式にも適用できます。

しかしながら、クラメルの公式は行列式を用いた解法であり、特に式の数(mm)と変数の数(nn)が等しい場合にのみ適用できる解法です。また、変数の数に応じて計算しなければならない行列式の数が増えていきます( (n+1)(n + 1) 個の行列式を計算する必要がある)。

このような理由から、クラメルの公式はあまり実用的でない場合がほとんどです。

行列の基本変形による解法

以下に整理する解法は行列の基本変形による解法であり、本節のはじめ(基本的な考え方)に示した掃き出し法の考え方に則っています。

解を得るために必要な手続きが明確であり、計算量も(比較的)少ないため、連立一次方程式の解法として実用的といえます。また、この解法はコンピュータによる計算にも適してます。

手順

AA(m,n)(m, n) 型行列とすると、次の手順(11\sim33)により連立一次方程式 Ax=bA \bm{x} = \bm{b} の解が得られます。

(1)解の有無の判定

連立一次方程式の解の有無の判定方法にしたがって、連立一次方程式 Ax=bA \bm{x} = \bm{b} が解を持つか否か判定します。詳細は、前項を参照ください。

(1-1)階段行列への変形
  • 拡大係数行列 (Ab)(\, A \mid \bm{b} \,) を階段行列に変形します。
  • 変形の結果得られる行列は、次のような形になります。
係数拡大行列の階段行列へ変形(一般の連立一次方連立の解の有無判定法)

(1-2)解の有無判定
  • br+1b_{r+1} の値により解の有無を判定します。

  • br+1=0b_{r+1} = 0 であれば、Ax=bA \bm{x} = \bm{b} は解を持ちます。

  • br+10b_{r+1} \neq 0 であれば、Ax=bA \bm{x} = \bm{b} は解を持ちません。

  • 連立一次方程式 Ax=bA \bm{x} = \bm{b} が解を持つ場合、次の手順22)行標準形への変形に進みます。

(2)行標準形への変形

  • 拡大係数行列 (Ab)(\, A \mid \bm{b} \,) を行標準形に変形します。
    • 11)解の有無の判定において階段行列まで変形が完了しています。
    • この階段行列に対して更に行基本変形と列の入れ替えの操作を行い、行標準形まで変形を進めます。
  • 変形の結果得られる行列 (Ab)(\, A^{\prime} \mid \bm{b}^{\prime} \,) は、次のような行列となります。
    • ここで、r=rankA=rank(Ab)r = \text{rank} \, A = \text{rank} \, (\, A \mid \bm{b} \,) が成り立ちます(定理 5.14(行標準形))。
(Ab)=(11a1r+1a1narr+1arnb1brOO) (\, A^{\prime} \mid \bm{b}^{\prime} \,) = \left( \begin{array} {c:c|c} \begin{matrix} 1 & & \\ & \ddots & \\ & & 1 \end{matrix} & \begin{matrix} a_{1 \, r+1} & \cdots & a_{1n} \\ \vdots & & \vdots \\ a_{r \, r+1} & \cdots & a_{r n} \end{matrix} & \begin{matrix} b_{1} \\ \vdots \\ b_{r} \end{matrix} \\ \hline \begin{matrix} & & \\ & \Large{O} & \\ & & \end{matrix} & \begin{matrix} & & \\ & & \Large{O} \\ & & \end{matrix} & \\ \end{array} \right)

(3)解の整理

  • 行標準形への変形により、簡単になった斉次連立一次方程式 Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} を解きます。
    • (Ab)(\, A \mid \bm{b} \,) から (Ab)(\, A^{\prime} \mid \bm{b}^{\prime} \,) への変形は可逆的であり(定理 5.7(基本変形の可逆性))、対応する連立一次方程式の変形も可逆的です(定理 5.16(基本変形の可逆性))。
    • したがって、Ax=bA \bm{x} = \bm{b}Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} は同じ解を持ちます。
  • Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} は次のような連立一次方程式となります。
{x1+  a1r+1xr+1++a1nxn=b1  x2+  a2r+1xr+1++a2nxn=b2            xr+  arr+1xr+1++arnxn=br    {  x1=b1(a1r+1xr+1++a1nxn)x2=b2(a2r+1xr+1++a2nxn)    xr=br(arr+1xr+1++arnxn) \begin{gather*} \left\{ \begin{array} {ll} x_{1} & + \; a_{1 \, r+1} x_{r+1} + \cdots + a_{1 n} x_{n} = b_{1} \\ \quad \; x_{2} & + \; a_{2 \, r+1} x_{r+1} + \cdots + a_{2 n} x_{n} = b_{2} \\ \quad \; \quad \; \ddots & \; \vdots \quad \quad \quad \quad \quad \, \vdots \\ \quad \; \quad \; \quad \; \quad x_{r} & + \; a_{r \, r+1} x_{r+1} + \cdots + a_{r n} x_{n} = b_{r} \\ \end{array} \right. \tag{5.4.3} \\ \\ \iff \left\{ \; \begin{split} x_{1} &= b_{1} - (\, a_{1 \, r+1} x_{r+1} + \cdots + a_{1 n} x_{n} \,) \\ x_{2} &= b_{2} - (\, a_{2 \, r+1} x_{r+1} + \cdots + a_{2 n} x_{n} \,) \\ & \; \; \vdots \\ x_{r} &= b_{r} - (\, a_{r \, r+1} x_{r+1} + \cdots + a_{r n} x_{n} \,) \\ \end{split} \right. \end{gather*}
  • ここで、dr+1,,dnKd_{r+1}, \cdots, d_{n} \in K を任意のスカラーとして、(nr)(n - r) 個の変数 xr+1,,xnx_{r+1}, \cdots, x_{n} について xr+1=dr+1,,xn=dnx_{r+1} = d_{r+1}, \cdots, x_{n} = d_{n} とすると、Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} の解は次のように表すことができます。
{  x1=b1(a1r+1dr+1++a1ndn)x2=b2(a2r+1dr+1++a2ndn)    xr=br(arr+1dr+1++arndn)xr+1=dr+1    xn=dn(5.4.4) \left\{ \; \begin{split} x_{1} &= b_{1} - (\, a_{1 \, r+1} d_{r+1} + \cdots + a_{1 n} d_{n} \,) \\ x_{2} &= b_{2} - (\, a_{2 \, r+1} d_{r+1} + \cdots + a_{2 n} d_{n} \,) \\ & \; \; \vdots \\ x_{r} &= b_{r} - (\, a_{r \, r+1} d_{r+1} + \cdots + a_{r n} d_{n} \,) \\ x_{r+1} &= d_{r+1} \\ & \; \; \vdots \\ x_{n} &= d_{n} \\ \end{split} \right. \tag{5.4.4}
  • 基本変形の可逆性から、これは Ax=0A \bm{x} = \bm{0} の解に他なりません。

以上から、斉次連立一次方程式 Ax=bA \bm{x} = \bm{b} の解を求めることができました。解法の要所は拡大係数行列 (Ab)(\, A \mid \bm{b} \,) を行標準形 (Ab)(\, A^{\prime} \mid \bm{b}^{\prime} \,) に変形するところであり、基本変形により得られた (Ab)(\, A^{\prime} \mid \bm{b}^{\prime} \,) の形により解が定まります。

考察

一般解と基本解

5.4.45.4.4)式より、連立一次方程式の解は、次のように (nr)+1(n - r) + 1 個の nn 項列ベクトルの線型結合として表すことができます。

(x1x2xrxr+1xr+2xn)=(b1b2br000)+dr+1(a1r+1a2r+1arr+1100)++dn(a1na2narn001) \begin{pmatrix} x_{1} \\ x_{2} \\ \vdots \\ x_{r} \\ x_{r+1} \\ x_{r+2} \\ \vdots \\ x_{n} \end{pmatrix} = \begin{pmatrix} \, b_{1} \, \\ \, b_{2} \, \\ \, \vdots \, \\ \, b_{r} \, \\ \, 0 \, \\ \, 0 \, \\ \, \vdots \, \\ \, 0 \, \end{pmatrix} + d_{r+1} \begin{pmatrix} -a_{1 \, r+1} \\ -a_{2 \, r+1} \\ \vdots \\ -a_{r \, r+1} \\ 1 \\ 0 \\ \vdots \\ 0 \end{pmatrix} + \cdots + d_{n} \begin{pmatrix} -a_{1 n} \\ -a_{2 \, n} \\ \vdots \\ -a_{r n} \\ 0 \\ 0 \\ \vdots \\ 1 \end{pmatrix}

それぞれの nn 項列ベクトルを左から x0,xr+1,,xn\bm{x}_{0}, \bm{x}_{r+1}, \cdots, \bm{x}_{n} とすると、Ax=bA \bm{x} = \bm{b} の解は次のように表すことができます。

x=x0+dr+1xr+1++dnxn \begin{equation*} \bm{x} = \bm{x}_{0} + d_{r + 1} \bm{x}_{r + 1} + \cdots + d_{n} \bm{x}_{n} \end{equation*}

これは、定理 5.6(連立一次方程式の解の形)に示した連立一次方程式の一般解の形に他なりません。

すなわち、Ax=bA \bm{x} = \bm{b} の解は、11 つの特殊解 x0\bm{x}_{0} と、Ax=bA \bm{x} = \bm{b} と同じ係数行列を持つ斉次連立一次方程式 Ax=0A \bm{x} = \bm{0} の基本解 xr+1,,xn\bm{x}_{r + 1}, \cdots, \bm{x}_{n} の線型結合との和として表すことができるということです。

斉次連立一次方程式の解との関係

拡大係数行列 (Ab)(\, A \mid \bm{b} \,) を行標準形に変形する操作は b\bm{b} によりません。(Ab)(\, A \mid \bm{b} \,) を変形して得られる行標準形 (Ab)(\, A^{\prime} \mid \bm{b}^{\prime} \,) の左側のブロック AA^{\prime} は、係数行列 AA を行標準形に変形して得られる行列 AA^{\prime} と等しくなります。

したがって、斉次連立一次方程式 Ax=0A \bm{x} = \bm{0} の任意の解を y\bm{y} とすれば、y\bm{y} は次のように表すことができます。

y=dr+1xr+1++dnxn \begin{equation*} \bm{y} = d_{r + 1} \bm{x}_{r + 1} + \cdots + d_{n} \bm{x}_{n} \end{equation*}

以上から、連立一次方程式 Ax=bA \bm{x} = \bm{b} の解 x\bm{x} と、同じ係数行列を持つ斉次連立一次方程式 Ax=0A \bm{x} = \bm{0} の解 y\bm{y} の間には次の関係が成り立ちます。

x=x0+y \bm{x} = \bm{x}_{0} + \bm{y}

この考察もまた、定理 5.5(連立一次方程式の解の集合)の主張に整合するものです。

解法の注意点

基本変形による解法の妥当性

行列の基本変形による解法において、もとの連立一次方程式 Ax=bA \bm{x} = \bm{b} と、変形後の Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} が同じ解を持つということは重要な点です。

連立一次方程式の変形が可逆的であり 22 つの連立一次方程式が同じ解を持つということが、基本変形によるこの解法に妥当性を与えているからです。

基本変形に対応する操作

行列の基本変形のすべての操作が可逆的であることは定理 5.7(基本変形の可逆性)により担保されます。

一方で、連立一次方程式の基本変形に関して、可逆的であることが証明されているのは行基本変形に対応する操作についてのみです(定理 5.16(基本変形の可逆性))。つまり、列基本変形に対応する連立一次方程式の変形に関しては、必ずしも可逆的であるとはいえません。

上記の解法では、行標準形を得るために列の入れ替えの操作を行いますが、これは問題ないのでしょうか?係数行列における列の入れ替えの操作は、連立一次方程式に対するどのような操作に対応しているのでしょうか?

列の入れ替えに対応する操作

結論からいえば、係数拡大行列における列の入れ替えは、対応する連立一次方程式における変数の和の順番の入れ替えに相当し、行っても問題ありません。

例えば、係数拡大行列において、次のように第 ii 列と第 jj 列を入れ替えることを考えます。

(a1ia2iamia1ja2jamjb1b2bm)(a1ja2jamja1ia2iamib1b2bm) \begin{gather*} & \left( \begin{array} {ccccc|c} \cdots & \fbox{a1ia2iami\begin{matrix} a_{1i} \\ a_{2i} \\ \vdots \\ a_{mi} \end{matrix}} & \cdots & \fbox{a1ja2jamj\begin{matrix} a_{1j} \\ a_{2j} \\ \vdots \\ a_{mj} \end{matrix}} & \cdots & \begin{matrix} b_{1} \\ b_{2} \\ \vdots \\ b_{m} \end{matrix} \end{array} \right) \\ \\ \longrightarrow & \left( \begin{array} {ccccc|c} \cdots & \fbox{a1ja2jamj\begin{matrix} a_{1j} \\ a_{2j} \\ \vdots \\ a_{mj} \end{matrix}} & \cdots & \fbox{a1ia2iami\begin{matrix} a_{1i} \\ a_{2i} \\ \vdots \\ a_{mi} \end{matrix}} & \cdots & \begin{matrix} b_{1} \\ b_{2} \\ \vdots \\ b_{m} \end{matrix} \end{array} \right) \end{gather*}

これは、対応する連立一次方程式において第 ii 項と第 jj 項を入れ替えることに相当します。

{+++a1ix1a2ixiamixi++++++a1jxja2jxjamjxj+=b1+=b2+=bm    {+++a1jxja2jxjamjxj++++++a1ix1a2ixiamixi+=b1+=b2+=bm \begin{gather*} & \left\{ \begin{array} {ccccc} \begin{matrix} \cdots + \\ \cdots + \\ \phantom{\vdots} \\ \cdots + \end{matrix} & \fbox{a1ix1a2ixiamixi\begin{matrix} a_{1i} x_{1} \\ a_{2i} x_{i} \\ \vdots \\ a_{mi} x_{i} \end{matrix}} & \begin{matrix} + \cdots + \\ + \cdots + \\ \phantom{\vdots} \\ + \cdots + \end{matrix} & \fbox{a1jxja2jxjamjxj\begin{matrix} a_{1j} x_{j} \\ a_{2j} x_{j} \\ \vdots \\ a_{mj} x_{j} \end{matrix}} & \begin{matrix} + \cdots = b_{1} \\ + \cdots = b_{2} \\ \phantom{\vdots} \\ + \cdots = b_{m} \end{matrix} \end{array} \right. \tag{\ast} \\ \\ & \iff \left\{ \begin{array} {ccccc} \begin{matrix} \cdots + \\ \cdots + \\ \phantom{\vdots} \\ \cdots + \end{matrix} & \fbox{a1jxja2jxjamjxj\begin{matrix} a_{1j} x_{j} \\ a_{2j} x_{j} \\ \vdots \\ a_{mj} x_{j} \end{matrix}} & \begin{matrix} + \cdots + \\ + \cdots + \\ \phantom{\vdots} \\ + \cdots + \end{matrix} & \fbox{a1ix1a2ixiamixi\begin{matrix} a_{1i} x_{1} \\ a_{2i} x_{i} \\ \vdots \\ a_{mi} x_{i} \end{matrix}} & \begin{matrix} + \cdots = b_{1} \\ + \cdots = b_{2} \\ \phantom{\vdots} \\ + \cdots = b_{m} \end{matrix} \end{array} \tag{\ast^{\prime}} \right. \end{gather*}

それぞれの方程式において、各項の和の順序は任意です。したがって(\ast)と(\ast^{\prime})はまったく同じ連立一次方程式といえます。

つまり、上記の解法の手順おいて(行基本変形に加えて)列の入れ替えを行っても問題なく、Ax=bA \bm{x} = \bm{b}Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} は同じ解を持つ連立一次方程式であるといえます。

定数項に対応する最右列(b\bm{b})は入れ替えてはいけない

ただし、拡大係数行列の最右列とその他の列を入れ替えてはいけません。

拡大係数行列 (Ab)(\, A \mid \bm{b} \,) の最右列にある b\bm{b} は、連立一次方程式の各方程式の右辺にある定数項(b1,,bmb_{1}, \cdots, b_{m})に対応しています。これを他の項と入れ替えてしまうと、拡大係数行列と連立一次方程式の対応が崩れてしまいます。

したがって、拡大係数行列を変形して行標準形を得る際は、行基本変形に加えて、最右列を除く列の入れ替えは行っても問題ない、といえます。

解の変数の順序

列の入れ替え操作を行った場合、注意すべき点があります。それは、得られた解における変数の順序です。

上記の\ast\ast^{\prime}をそれぞれ Ax=bA \bm{x} = \bm{b}Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} と表したとき、係数拡大行列 (Ab)(\, A \mid \bm{b} \,) から (Ab)(\, A^{\prime} \mid \bm{b}^{\prime} \,) への変形に伴って(解にあたる)変数ベクトル x\bm{x} がどのように変化するか考えてみましょう。

注意深くみれば、もとの\astの解に対応する x\bm{x} に対して、変形後の\ast^{\prime}の解に対応する x\bm{x}^{\prime} では第 ii 項と第 jj 項が入れ替わっていることがわかります。

x=(xixj),x=(xjxi) \begin{array} {ccc} \bm{x} = \begin{pmatrix} \, \vdots \, \\ \, x_{i} \, \\ \, \vdots \, \\ \, x_{j} \, \\ \, \vdots \, \end{pmatrix}, & & \bm{x}^{\prime} = \begin{pmatrix} \, \vdots \, \\ \, x_{j} \, \\ \, \vdots \, \\ \, x_{i} \, \\ \, \vdots \, \end{pmatrix} & \end{array}

すなわち、係数拡大行列を行標準形に変形する際に列の入れ替えを行った場合、得られた解において、対応する変数の順序も(列の順序に合わせて)入れ替わっていなければなりません。

基本変形による解法を用いて連立一次方程式を解く際は、この点を忘れないよう充分注意する必要があります。


解法の例

行列の基本変形による解法の例として、次の 22 つの連立一次方程式を行列の基本変形を用いて解いてみます。


例題1(行基本変形のみの場合)

次の連立一次方程式を解け。

{  xy+3z=1y+z=32xy+5z=1 \begin{align*} \left\{ \; \begin{alignat*} {5} & & & x & & \phantom{y} + {} & 3 & z = {} & {} & 1 \\ & & & & & y + {} & & z = {} & {} & 3 \\ & & 2 & x - {} & & y + {} & 5 & z = {} & - {} & 1 \\ \end{alignat*} \right. \end{align*}



解答(例題1)

与えられた連立一次方程式を Ax=bA \bm{x} = \bm{b} のように表すと、拡大係数行列 (Ab)(\, A \mid \bm{b} \,) は、行基本変形により次のような階段行列に変形することができる。

(Ab)=(103101132151)(i)  (103101130113)(ii)  (103101130000) \begin{split} (\, A \mid \bm{b} \,) &= \left( \begin{array} {ccc|c} 1 & 0 & 3 & 1 \\ 0 & 1 & 1 & 3 \\ 2 & -1 & 5 & -1 \\ \end{array} \right) \\ &\overset{(\text{i})}{\longrightarrow} \; \left( \begin{array} {ccc|c} 1 & 0 & 3 & 1 \\ 0 & 1 & 1 & 3 \\ 0 & -1 & -1 & -3 \\ \end{array} \right) \\ &\overset{(\text{ii})}{\longrightarrow} \; \left( \begin{array} {ccc|c} 1 & 0 & 3 & 1 \\ 0 & 1 & 1 & 3 \\ 0 & 0 & 0 & 0 \\ \end{array} \right) \\ \end{split}

したがって、係数行列 AA の階数と拡大係数行列 (Ab)(\, A \mid \bm{b} \,) の階数はともに 22 に等しく、Ax=bA \bm{x} = \bm{b} は解を持つ。上記の変形により得られた拡大係数行列は既に行標準形であり、対応する斉次連立一次方程式は次の通り。また、変形は可逆的であるからこれは Ax=bA \bm{x} = \bm{b} と同じ解をもつ。

{  xy+3z=1x+y+z=3 \begin{align*} \left\{ \; \begin{alignat*} {5} & & & x & & \phantom{y +} & 3 & z = {} & {} & 1 \\ & & & \phantom{x +} & & y + {} & & z = {} & {} & 3 \\ \end{alignat*} \right. \end{align*}

したがって、dd を任意のスカラーとすると Ax=bA \bm{x} = \bm{b} の解は次のとおり。

{    x=13dy=3dz=d \begin{gather*} \left\{ \; \; \begin{split} x &= 1 - 3 d \\ y &= 3 - d \\ z &= d \\ \end{split} \right. \end{gather*}

これをベクトルとして表せば次のようになる。

x=(130)+d(311) \bm{x} = \begin{pmatrix} \, 1 \, \\ \, 3 \, \\ \, 0 \, \end{pmatrix} + d \begin{pmatrix} \, -3 \, \\ \, -1 \, \\ \, 1 \, \end{pmatrix}



解答の考え方(例題1)

基本変形による解法の手順に従います。(11)拡大係数行列 (Ab)(\, A \mid \bm{b} \,) を階段行列に変形して解の有無を判定した後に(22)更に行標準形 (Ab)(\, A^{\prime} \mid \bm{b}^{\prime} \,) まで標準化を進め(33)簡単になった連立一次方程式 Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} を解きます。

(1)解の有無の判定

階段行列への変形は次のような手順で行います。

(Ab)=(103101132151)(i)  (103101130113)(ii)  (103101130000) \begin{split} (\, A \mid \bm{b} \,) &= \left( \begin{array} {ccc|c} 1 & 0 & 3 & 1 \\ 0 & 1 & 1 & 3 \\ 2 & -1 & 5 & -1 \\ \end{array} \right) \\ &\overset{(\text{i})}{\longrightarrow} \; \left( \begin{array} {ccc|c} 1 & 0 & 3 & 1 \\ 0 & 1 & 1 & 3 \\ 0 & -1 & -1 & -3 \\ \end{array} \right) \\ &\overset{(\text{ii})}{\longrightarrow} \; \left( \begin{array} {ccc|c} 1 & 0 & 3 & 1 \\ 0 & 1 & 1 & 3 \\ 0 & 0 & 0 & 0 \\ \end{array} \right) \\ \end{split}

変形の結果得られる階段行列は次のようになります。

(103101130000) \left( \begin{array} {ccc|c} 1 & 0 & 3 & 1 \\ 0 & 1 & 1 & 3 \\ \hline 0 & 0 & 0 & 0 \\ \end{array} \right)

  • 左側のブロックは係数行列 AA を変形して得られる階段行列に等しく、AA の階数が 22 であることがわかります(定理 5.11(階段行列))。
  • 右下のブロックをみると、(3,4)(3, 4) 成分が 00 に等しいことがわかります。すなわち、拡大係数行列 (Ab)(\, A \mid \bm{b} \,) の階数も 22 であり、係数行列の階数に等しくなります。

よって、rankA=rank(Ab)\text{rank} \, A = \text{rank} \, (\, A \mid \bm{b} \,) が成り立ち、Ax=bA \bm{x} = \bm{b} が解を持つための条件(定理 5.1(連立一次方程式が解を持つための条件))を満たすことが確かめられました。

(2)行標準形への変形

11)解の有無の判定で得られた階段行列は、既に行標準形となっています。この場合、更なる変形の操作は不要です。

(3)解の整理

変形により得られた連立一次方程式 Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} は次の通りです。

{  xy+3z=1x+y+z=3    {  x=13zy=3z \begin{gather*} \left\{ \; \begin{alignat*} {5} & & & x & & \phantom{y +} & 3 & z = {} & {} & 1 \\ & & & \phantom{x +} & & y + {} & & z = {} & {} & 3 \\ \end{alignat*} \right. \\ \iff \left\{ \; \begin{split} x &= 1 - 3 z \\ y &= 3 - z \\ \end{split} \right. \end{gather*}

ここで、dd を任意定数として z=dz = d   (dK)\; (\, d \in K \,) とすることで解が得られます。

{    x=13zy=3zz=d \begin{gather*} \left\{ \; \; \begin{split} x &= 1 - 3 z \\ y &= 3 - z \\ z &= d \\ \end{split} \right. \end{gather*}



例題2(列の入れ替えを含む場合)

次の連立一次方程式を解け。

{  3x2+3x32x4=4x1+x2+2x3+3x4=2x1+2x2+3x3+2x4=1x1+3x2+4x3+2x4=1 \begin{align*} \left\{ \; \begin{alignat*} {7} & & & & 3 & x_{2} + {} & 3 & x_{3} - {} & 2 & x_{4} = {} & - {} & 4 \\ & & & x_{1} + {} & & x_{2} + {} & 2 & x_{3} + {} & 3 & x_{4} = {} & & 2 \\ & & & x_{1} + {} & 2 & x_{2} + {} & 3 & x_{3} + {} & 2 & x_{4} = {} & & 1 \\ & & & x_{1} + {} & 3 & x_{2} + {} & 4 & x_{3} + {} & 2 & x_{4} = {} & - {} & 1 \\ \end{alignat*} \right. \end{align*}



解答(例題2)

与えられた連立一次方程式を Ax=bA \bm{x} = \bm{b} のように表すと、拡大係数行列 (Ab)(\, A \mid \bm{b} \,) は、行基本変形により次のような階段行列に変形することができる。

(Ab)=(03324112321232113421)(i)  (11232033240111102213)(ii)  (11232011110001100011)(iii)  (11232011110001100000) \begin{split} (\, A \mid \bm{b} \,) &= \left( \begin{array} {cccc|c} 0 & 3 & 3 & -2 & -4 \\ 1 & 1 & 2 & 3 & 2 \\ 1 & 2 & 3 & 2 & 1 \\ 1 & 3 & 4 & 2 & -1 \\ \end{array} \right) \\ &\overset{(\text{i})}{\longrightarrow} \; \left( \begin{array} {cccc|c} 1 & 1 & 2 & 3 & 2 \\ 0 & 3 & 3 & -2 & -4 \\ 0 & 1 & 1 & -1 & -1 \\ 0 & 2 & 2 & -1 & -3 \\ \end{array} \right) \\ &\overset{(\text{ii})}{\longrightarrow} \; \left( \begin{array} {cccc|c} 1 & 1 & 2 & 3 & 2 \\ 0 & 1 & 1 & -1 & -1 \\ 0 & 0 & 0 & 1 & -1 \\ 0 & 0 & 0 & 1 & -1 \\ \end{array} \right) \\ &\overset{(\text{iii})}{\longrightarrow} \; \left( \begin{array} {cccc|c} 1 & 1 & 2 & 3 & 2 \\ 0 & 1 & 1 & -1 & -1 \\ 0 & 0 & 0 & 1 & -1 \\ 0 & 0 & 0 & 0 & 0 \\ \end{array} \right) \\ \end{split}

したがって、係数行列 AA の階数と拡大係数行列 (Ab)(\, A \mid \bm{b} \,) の階数はともに 33 に等しく、Ax=bA \bm{x} = \bm{b} は解を持つ。また、(Ab)(\, A \mid \bm{b} \,) は、更に行基本変形と列の入れ替えを行うことで、行標準形に変形することができる。

(Ab)(iv)  (10143011110001100000)(v)  (10107011020001100000)(vi)  (10017010120010100000) \begin{split} \phantom{(\, A \mid \bm{b} \,)} & \\ &\overset{(\text{iv})}{\longrightarrow} \; \left( \begin{array} {cccc|c} 1 & 0 & 1 & 4 & 3 \\ 0 & 1 & 1 & -1 & -1 \\ 0 & 0 & 0 & 1 & -1 \\ 0 & 0 & 0 & 0 & 0 \\ \end{array} \right) \\ &\overset{(\text{v})}{\longrightarrow} \; \left( \begin{array} {cccc|c} 1 & 0 & 1 & 0 & 7 \\ 0 & 1 & 1 & 0 & -2 \\ 0 & 0 & 0 & 1 & -1 \\ 0 & 0 & 0 & 0 & 0 \\ \end{array} \right) \\ &\overset{(\text{vi})}{\longrightarrow} \; \left( \begin{array} {cccc|c} 1 & 0 & 0 & 1 & 7 \\ 0 & 1 & 0 & 1 & -2 \\ 0 & 0 & 1 & 0 & -1 \\ 0 & 0 & 0 & 0 & 0 \\ \end{array} \right) \\ \end{split}

vi\text{vi})で第 33 列と第 44 列を入れ替えたので、変形により得られた斉次連立一次方程式は次のようになり、これは Ax=bA \bm{x} = \bm{b} と同じ解をもつ。

{  x1x2+x4++x3=7x1+x2x4++x3=2x1+x2+x4x3=1 \begin{align*} \left\{ \; \begin{alignat*} {7} & & & x_{1} & & \phantom{x_{2} +} & & \phantom{x_{4} +} + {} & & x_{3} = {} & {} & 7 \\ & & & \phantom{x_{1} +} & & x_{2} & & \phantom{x_{4} +} + {} & & x_{3} = {} & - {} & 2 \\ & & & \phantom{x_{1} +} & & \phantom{x_{2} +} & & x_{4} & & \phantom{x_{3}} = {} & - {} & 1 \\ \end{alignat*} \right. \end{align*}

したがって、dd を任意のスカラーとすると Ax=bA \bm{x} = \bm{b} の解は次のとおり。

{    x1=7dx2=2dx3=dx4=1 \begin{gather*} \left\{ \; \; \begin{split} x_{1} &= 7 - d \\ x_{2} &= -2 - d \\ x_{3} &= d \\ x_{4} &= -1 \\ \end{split} \right. \end{gather*}

これをベクトルとして表せば次のようになる。 x=(7201)+d(1110) \bm{x} = \begin{pmatrix} \, 7 \, \\ \, -2 \, \\ \, 0 \, \\ \, -1 \, \end{pmatrix} + d \begin{pmatrix} \, -1 \, \\ \, -1 \, \\ \, 1 \, \\ \, 0 \, \end{pmatrix}



解答の考え方(例題2)

基本変形による解法の手順に従います。(11)拡大係数行列 (Ab)(\, A \mid \bm{b} \,) を階段行列に変形して解の有無を判定した後に(22)更に行標準形 (Ab)(\, A^{\prime} \mid \bm{b}^{\prime} \,) まで標準化を進め(33)簡単になった連立一次方程式 Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} を解きます。

(1)解の有無の判定

階段行列への変形は次のような手順で行います。なお、解の有無判定までは前項の例題1とまったく同じです。

(Ab)=(03324112321232113421)(i)  (11232033240111102213)(ii)  (11232011110001100011)(iii)  (11232011110001100000) \begin{split} (\, A \mid \bm{b} \,) &= \left( \begin{array} {cccc|c} 0 & 3 & 3 & -2 & -4 \\ 1 & 1 & 2 & 3 & 2 \\ 1 & 2 & 3 & 2 & 1 \\ 1 & 3 & 4 & 2 & -1 \\ \end{array} \right) \\ &\overset{(\text{i})}{\longrightarrow} \; \left( \begin{array} {cccc|c} 1 & 1 & 2 & 3 & 2 \\ 0 & 3 & 3 & -2 & -4 \\ 0 & 1 & 1 & -1 & -1 \\ 0 & 2 & 2 & -1 & -3 \\ \end{array} \right) \\ &\overset{(\text{ii})}{\longrightarrow} \; \left( \begin{array} {cccc|c} 1 & 1 & 2 & 3 & 2 \\ 0 & 1 & 1 & -1 & -1 \\ 0 & 0 & 0 & 1 & -1 \\ 0 & 0 & 0 & 1 & -1 \\ \end{array} \right) \\ &\overset{(\text{iii})}{\longrightarrow} \; \left( \begin{array} {cccc|c} 1 & 1 & 2 & 3 & 2 \\ 0 & 1 & 1 & -1 & -1 \\ 0 & 0 & 0 & 1 & -1 \\ 0 & 0 & 0 & 0 & 0 \\ \end{array} \right) \\ \end{split}

変形の結果得られる階段行列は次のようになります。

(11232011110001100000) \left( \begin{array} {cccc|c} 1 & 1 & 2 & 3 & 2 \\ 0 & 1 & 1 & -1 & -1 \\ 0 & 0 & 0 & 1 & -1 \\ \hline 0 & 0 & 0 & 0 & 0 \\ \end{array} \right)

  • 左側のブロックは係数行列 AA を変形して得られる階段行列に等しく、AA の階数が 33 であることがわかります(定理 5.11(階段行列))。
  • 右下のブロックをみると、(4,5)(4, 5) 成分が 00 に等しいことがわかります。すなわち、拡大係数行列 (Ab)(\, A \mid \bm{b} \,)33 であり、係数行列 AA の階数に等しくなります。

よって、rankA=rank(Ab)\text{rank} \, A = \text{rank} \, (\, A \mid \bm{b} \,) が成り立ち、Ax=bA \bm{x} = \bm{b} が解を持つための条件(定理 5.1(連立一次方程式が解を持つための条件))を満たすことが確かめられました。

(2)行標準形への変形

上記の階段行列に対して、更に行基本変形と列の入れ替えを行い、行標準形まで標準化を進めます。行標準形への変形は次のような手順で行います。

(Ab)  (11232011110001100000)(iv)  (10143011110001100000)(v)  (10107011020001100000)(vi)  (10017010120010100000) \begin{split} \phantom{(\, A \mid \bm{b} \,)} & \\ & \; \left( \begin{array} {cccc|c} 1 & 1 & 2 & 3 & 2 \\ 0 & 1 & 1 & -1 & -1 \\ 0 & 0 & 0 & 1 & -1 \\ 0 & 0 & 0 & 0 & 0 \\ \end{array} \right) \\ &\overset{(\text{iv})}{\longrightarrow} \; \left( \begin{array} {cccc|c} 1 & 0 & 1 & 4 & 3 \\ 0 & 1 & 1 & -1 & -1 \\ 0 & 0 & 0 & 1 & -1 \\ 0 & 0 & 0 & 0 & 0 \\ \end{array} \right) \\ &\overset{(\text{v})}{\longrightarrow} \; \left( \begin{array} {cccc|c} 1 & 0 & 1 & 0 & 7 \\ 0 & 1 & 1 & 0 & -2 \\ 0 & 0 & 0 & 1 & -1 \\ 0 & 0 & 0 & 0 & 0 \\ \end{array} \right) \\ &\overset{(\text{vi})}{\longrightarrow} \; \left( \begin{array} {cccc|c} 1 & 0 & 0 & 1 & 7 \\ 0 & 1 & 0 & 1 & -2 \\ 0 & 0 & 1 & 0 & -1 \\ 0 & 0 & 0 & 0 & 0 \\ \end{array} \right) \\ \end{split}

変形の結果得られる行標準形 (Ab)(\, A^{\prime} \mid \bm{b}^{\prime} \,) は次のようになります。

A=(10017010120010100000) A^{\prime} = \left( \begin{array} {cccc|c} 1 & 0 & 0 & 1 & 7 \\ 0 & 1 & 0 & 1 & -2 \\ 0 & 0 & 1 & 0 & -1 \\ \hline 0 & 0 & 0 & 0 & 0 \\ \end{array} \right)

(3)解の整理

変形により得られた連立一次方程式 Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} は次の通りです。

{  x1x2+x4++x3=7x1+x2x4++x3=2x1+x2+x4x3=1 \begin{gather*} \left\{ \; \begin{alignat*} {7} & & & x_{1} & & \phantom{x_{2} +} & & \phantom{x_{4} +} + {} & & x_{3} = {} & {} & 7 \\ & & & \phantom{x_{1} +} & & x_{2} & & \phantom{x_{4} +} + {} & & x_{3} = {} & - {} & 2 \\ & & & \phantom{x_{1} +} & & \phantom{x_{2} +} & & x_{4} & & \phantom{x_{3}} = {} & - {} & 1 \\ \end{alignat*} \right. \end{gather*}

  • 例題1と異なり、行基本変形に加えて、列の入れ替えの操作を行うことで行標準形への変形が完了します。
  • 変形(vi\text{vi})における第 33 列と第 44 列を入れ替えに伴って、変数の順序(x3x_{3}x4x_{4} の順序)も入れ替わっている点に注意が必要です。

Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} を解くと、次のようになります。

{  x1x2+x4++x3=7x1+x2x4++x3=2x1+x2+x4x3=1    {  x1=7x3x2=2x3x4=1 \begin{gather*} \left\{ \; \begin{alignat*} {7} & & & x_{1} & & \phantom{x_{2} +} & & \phantom{x_{4} +} + {} & & x_{3} = {} & {} & 7 \\ & & & \phantom{x_{1} +} & & x_{2} & & \phantom{x_{4} +} + {} & & x_{3} = {} & - {} & 2 \\ & & & \phantom{x_{1} +} & & \phantom{x_{2} +} & & x_{4} & & \phantom{x_{3}} = {} & - {} & 1 \\ \end{alignat*} \right. \\ \iff \left\{ \; \begin{split} x_{1} &= 7 - x_{3} \\ x_{2} &= -2 - x_{3} \\ x_{4} &= - 1 \\ \end{split} \right. \end{gather*}

ここで、dd を任意定数として x3=dx_{3} = d   (dK)\; (\, d \in K \,) とすることで解が得られます。

{    x1=7dx2=2dx3=dx4=1 \begin{gather*} \left\{ \; \; \begin{split} x_{1} &= 7 - d \\ x_{2} &= -2 - d \\ x_{3} &= d \\ x_{4} &= -1 \\ \end{split} \right. \end{gather*}


まとめ

  • 一般の連立一次方程式(Ax=bA \bm{x} = \bm{b})は、拡大係数行列を行標準形に変形することで解を得ることができる。

  • 解法の手順は次の通り。

    11)係数拡大行列 (Ab)(\, A \mid \bm{b} \,) を階段行列に変形し、解の有無を判定する。
    22)係数拡大行列を行標準形 (Ab)(\, A^{\prime} \mid \bm{b}^{\prime} \,) まで変形する。
    33)簡単になった連立一次方程式 Ax=bA^{\prime} \bm{x}^{\prime} = \bm{b}^{\prime} を解いて、解を整理する。

  • 手順(22)行標準形への変形において列の入れ替えの操作を行った場合、連立方程式の対応する変数の順序を入れ替える必要がある。

  • 手順(22)行標準形への変形において 、連立一次方程式の定数項に対応する最右列(b\bm{b})と他の列を入れ替えてはいけない。


参考文献

[1] 齋藤正彦. 線型代数入門. 東京大学出版会. 1966.
[2] 永田雅宣 他. 理系のための線型代数の基礎. 紀伊國屋書店. 1986.
[3] 川久保勝夫. 線形代数学 [新装版]. 日本評論社. 2010.
[4] 松坂和夫. 線型代数入門 [新装版]. 岩波書店. 2018.
[5] 三宅敏恒. 線形代数学 初歩からジョルダン標準形へ. 培風館. 2008.
[6] S. Lang. Linear Algebra Third Edition. Springer. 1987.
[7] T. Miyake. Linear Algebra From the Beginnings to the Jordan Normal. Springer. 2022.
[8] 雪江明彦. 代数学 11 群論入門. 日本評論社. 2010.
[9] 雪江明彦. 代数学 22 環と体とガロア理論. 日本評論社. 2010.
[10] 桂利行. 代数学 I\text{I} 群と環. 東京大学出版会. 2004.
[11] 松坂和夫. 代数系入門. 岩波書店. 1976.
[12] 高木貞治. 代数学講義 [改訂新版]. 共立出版. 1965.
[13] S. Lang. Algebra Revised Third Edition. Springer. 2002.
[14] M. Artin. Algebra Second Edition. Pearson Education Limited. 2014.
[15] 青本和彦 他. 数学入門辞典. 岩波書店. 2005.


初版:2023-07-29   |   改訂:2024-12-04