線形代数I/連立一次方程式 のバックアップの現在との差分(No.15)

更新


  • 追加された行はこの色です。
  • 削除された行はこの色です。
#contents

培風館「教養の線形代数(五訂版)」に沿って行っている授業の授業ノート(の一部)です。

* 2.1 連立一次方程式と行列 [#q942c267]
** 連立一次方程式とは [#i953febc]

どういう意味か?

:連立|複数ある
:一次方程式|全ての項の次数が最大でも一次である方程式~
-連立 = 複数ある
-一次方程式 = 全ての項の次数が最大でも一次である方程式

「項」とは?

:項|+やーで繋がれたもの~
&math(3x-2ay+(c+d)); なら &math(3x); や &math(-2ay); の他、&math((c+d)); を項と呼ぶことも
:因子|積で繋がれたもの~
&math(ax); なら &math(a); や &math(x);
:項| +やーで繋がれたもの~
 &math(3x-2ay+(c+d)); なら &math(3x); や &math(-2ay); の他、&math((c+d)); を項と呼ぶことも
:因子| 積で繋がれたもの~
 &math(ax); なら &math(a); や &math(x);

「次数」とは?

   項を構成する因子に「変数が」何個含まれているか~
   変数 &math(x,y,z,w,\dots);, 定数 &math(-3,\sqrt{2},a,b,c,\dots);

「一次」とは?

:0次| &math(3,\ -\sqrt{5},\ a,\ 3b,\ \cdots);
:一次| &math(x,\ -2y, 5z,\ w/\sqrt{2},\ \cdots ); ← 変数が1つ含まれたこういうやつ
:二次| &math(-2x^2,\ ay^2,\ 5xy\, \cdots);
:三次| &math(2x^3,\ xy^2,\ \pi xyz,\ \cdots);
:0次| &math(3,\ -\sqrt{5},\ a,\ 3b,\ \cdots);
:一次| &math(x,\ -2y, 5z,\ w/\sqrt{2},\ \cdots ); ← 変数が1つ含まれたこういうやつ
:二次| &math(-2x^2,\ ay^2,\ 5xy\, \cdots);
:三次| &math(2x^3,\ xy^2,\ \pi xyz,\ \cdots);
:……|

∴「項の次数が0次か一次の方程式が複数ある」=「連立一次方程式」

** 連立一次方程式の例 [#qa2db057]

たとえば、
&math(\left\{\begin{array}{l} x+2y=-1 \\ 3x+y=2 \end{array}\right .); 
を解くと、
&math(\left\{\begin{array}{l} x=1 \\ y=-1 \end{array}\right .);

** うまく解けない場合もある [#q01d398b]

- 解がたくさん出てくる
- 解が1つもない

*** 解がたくさん出てくる [#oa5d8bb5]

&math(\left\{\begin{array}{l} x+2y=-1 \\ 2x+4y=-2 \end{array}\right .); 
を解くと、
&math(\left\{\begin{array}{l} x=-2a-1 \\ y=a \end{array}\right .);~
任意の &math(a); について解になる。

(下の式が上の式のちょうど2倍になっているため、
実効的には式が1つしかないのと同じ事になっている)

*** 解が1つもない [#xb5593b1]

&math(\left\{\begin{array}{l} x+2y=-1 \\ 2x+4y=-3 \end{array}\right .); 

下の式から上の式の2倍を引くと次式を得る。

&math(0=-1);~

どんな &math(x,y); を持ってきてもこの式を満たすことはできない。

したがって「解無し」

** 連立一次方程式の一般形 [#nb322acb]

変数の数も、式の数も、いくつでも良い。

&math(\left\{\begin{array}{c@{}c@{}c@{}c@{}c@{}c@{}c} a_{11}x_1&+&a_{12}x_2&+\cdots+&a_{1n}x_n&=&b_1 \\ a_{21}x_1&+&a_{22}x_2&+\cdots+&a_{2n}x_n&=&b_2 \\ \vdots&&\vdots&&\vdots&&\vdots\\ a_{m1}x_1&+&a_{m2}x_2&+\cdots+&a_{mn}x_n&=&b_m \\ \end{array}\right .);

&math(x_1,x_2,\cdots,x_n); が変数~
&math(a_{11},\cdots,a_{mn}); が定数(係数)~
&math(b_1,b_2,\cdots,b_m); も定数

変数に使う文字:~
変数が3つまでであれば、&math(x,y,z); が使われ、~
変数が4つまでであれば、&math(x,y,z,w); が使われる事が多い。

それ以上の時は、&math(x_i); の形で表すとよい。

** ベクトル方程式として表す [#n192d560]

係数行列、変数ベクトル、定数ベクトル、をそれぞれ

&math(A=\begin{bmatrix}a_{11}&a_{12}&\cdots&a_{1n}\\a_{21}&a_{22}&\cdots&a_{2n}\\\vdots&\vdots&\ddots&\vdots\\a_{m1}&a_{m2}&\cdots&a_{mn}\end{bmatrix});、
&math(\bm x = \begin{bmatrix}x_1\\x_2\\\vdots\\x_n\end{bmatrix});、
&math(\bm b = \begin{bmatrix}b_1\\b_2\\\vdots\\b_m\end{bmatrix});~

と置くと、連立一次方程式は次のベクトル方程式の形に書ける。

&math(A\bm x=\bm b);

** 方程式の分類 [#f5014358]

任意の連立一次方程式は、先の例と同様に、

- 1つだけに決まる
- たくさんある
- 1つもない
- 解が1つだけに決まる
- 解がたくさんある
- 解が1つもない

のいずれかに分類できる。

&math(A\bm x=\bm b); なら &math(\bm x=A^{-1}\bm b); 
だから、解は必ず1つだけ求まるんじゃないの?~
&math(A\bm x=\bm b); なら両辺に左から &math(A^{-1}); をかければ &math(\bm x=A^{-1}\bm b); 
となるから、解は必ず1つだけ求まるんじゃないの?~
というのはちょっと甘い。

- &math(A); が正則でない場合もあるし
- そもそも &math(A); が正方行列にならない場合も考えないと

** 本章の目的 [#x0a1ce96]
* ガウスの消去法(定理2.2 変形定理) [#e5502340]

連立一次方程式はガウスの消去法を用いることで必ず解くことができる。

** 本項目の目的 [#x0a1ce96]

物理学やその他のシミュレーションでは、
何万、何十万の変数を持つ連立一次方程式を解かなければならない。

- 大きな次数の計算は「行き当たりばったり」では困るので、
「必ず解ける解法」を知りたい
- 方程式の性質を行列を使って理解したい

というのが、以降で扱う内容。
というのが、以降で扱う内容になる。

(計算機で必要になる「本当に効率の良い計算方法」はやらないけど)
(計算機で必要になる「本当に効率の良い計算方法」はやらない)

必ずしも「計算を楽にするため」に勉強しているのではなく、
その後に続くものなので、「別に他のやり方でも解けるし」
などと考えることなく、ここで学ぶ解法をしっかりマスターして欲しい。
「ガウスの消去法」は一次連立方程式を解くための1つの方法ではあるが、
これは必ずしも「計算を楽にするため」に勉強しているのではなく、
「後の話を理解するために必要だから」勉強しているものなので、
「別に他のやり方でも解けるし」などと考えることなく、
この解法をしっかりマスターして欲しい。

* ガウスの消去法(定理2.2 変形定理) [#e5502340]

** 連立一次方程式を解く(普通に) [#l59c1a6b]

&math(\left\{\!\begin{array}{r@{\,}c@{\,}r@{\,}c@{\,}r@{\,}c@{\,}r@{\ \ \ \ \ \ \ }c}x&-&2y&+&3z&=&1&\MARU{1}\\3x&+&y&-&5z&=&-4&\MARU{2} \\-2x&+&6y&-&9z&=&-2&\MARU{3}\end{array}\right .);

を解こう。

&math(\begin{array}{r@{\,}c@{\,}r@{\,}c@{\,}r@{\,}c@{\,}r@{\ \ \ \ \ \ \ }ll}\phantom{-2x}&\phantom{-}&7y&-&14z&=&-7&\MARU{2}-3\times\MARU{1}&=\MARU{4}\\&&2y&-&3z&=&0&\MARU{3}+2\times\MARU{1}&=\MARU{5} \\&&y&-&2z&=&-1&\MARU{4}/7&=\MARU{6}\\x&&&-&z&=&-1&\MARU{1}+2\times\MARU{6}&=\MARU{7}\\&&&&z&=&2&\MARU{5}-2\times\MARU{6}&=\MARU{8}\\x&&&&&=&1&\MARU{7}+\MARU{8}&=\MARU{9}\\&&y&&&=&3&\MARU{6}+2\times\MARU{8}&=\MARU{10}\\\end{array});

∴&math(\MARU{9}\MARU{10}\MARU{8});より &math(x=1,\ y=3,\ z=2);

あるいは、

&math(\begin{bmatrix}x\\y\\z\end{bmatrix}=\begin{bmatrix}1\\3\\2\end{bmatrix});

が解となる。

** どうやったか? [#e099220e]

連立方程式の解を変化させずに以下の式変形が可能:
連立方程式の解を変化させずに以下の式変形(同値変形)が可能:
- ある式を定数倍(&math(\neq 0);)する
- ある式を定数倍して別の式に加える
- (ある式を別の式を入れ替える) ←必要なら
- (ある式と別の式との順番を入れ替える) ←必要なら

与えられた式に上記の操作を繰り返して、最終的に

&math(\left\{\begin{array}{c@{\ }c@{\ }c@{\ }c@{\ }c}x&&&=&?\\&y&&=&?\\&&z&=&?\end{array}\right .);

の形にできれば、解けたことになる。
の形が得られれば、解けたことになる。

** 係数だけ取り出して考える [#gd31968b]

&math(\left\{\!\begin{array}{r@{\,}c@{\,}r@{\,}c@{\,}r@{\,}c@{\,}r}x&-&2y&+&3z&=&1\\3x&+&y&-&5z&=&-4 \\-2x&+&6y&-&9z&=&-2\end{array}\right .);

のような連立一次方程式に対して、係数行列と定数ベクトルを並べた

&math(A^*=[\,A\ \bm b\,]=\begin{bmatrix}1&-2&3&1\\3&1&-5&-4\\-2&6&-9&-2\end{bmatrix});~
を「拡大係数行列」と呼ぶ。

これに、
- ある行を定数倍(&math(\neq 0);)する
- ある行を定数倍して別の行に加える
- ある行を別の行と入れ替える

を行っても、対応する方程式の解は変わらない。~
→ これら3つは「(行に対する)基本変形」と呼ばれる

したがって、拡大係数行列に基本変形を繰り返すことで、

&math(\left[\begin{array}{ccccl}1&0&\cdots&0&b_1^\prime\\0&1&&\vdots&\vdots\\\vdots&&\ddots&0&b_{m-1}^\prime\\0&\cdots&0&1&b_m^\prime\end{array}\right]=[\,I\ \bm b^\prime\,]);

のように 「単位行列と(始めと異なる)定数ベクトルが並ぶ形」 にできれば、
これは、

&math(\left\{\begin{array}{c@{\ }c@{\ }c@{\ }c@{\ }c}x&&&=&b_1^\prime\\&y&&=&b_2^\prime\\&&z&=&b_3^\prime\end{array}\right .);

のような式に対応し、

&math(\bm x=\bm b^\prime);

が解となる。

(必ずこううまく行くとは限らないけど、うまく行った=解けた、である)

*** 注意 [#b6685d64]
繰り返しになるが、「別に係数だけを取り出した行列を作らなくても計算できるし、その方が分かりやすい!」
と思っても、ここで説明される方法の通りの手順で解けるようになるよう練習して欲しい。

なぜならここで学ぶ「解き方」は必ずしも「楽に解く」ためのものではなく、
連立一次方程式の持つ性質を理解する助けとするためのものだから。

** 実際にやってみる → ガウスの消去法(掃出し法) [#i2e196a4]

先の問題を、拡大係数行列を変形する方法で解いてみる。

ここで採用する方法は「ガウスの消去法」あるいは「掃出し法」と呼ばれる。
ここで採用する方法は「ガウスの消去法」あるいは「&ruby(はきだ){掃出};し法」と呼ばれる。

&math(\left\{\!\begin{array}{r@{\,}c@{\,}r@{\,}c@{\,}r@{\,}c@{\,}r}x&-&2y&+&3z&=&1\\3x&+&y&-&5z&=&-4 \\-2x&+&6y&-&9z&=&-2\end{array}\right .);

より、拡大係数行列は

&math([\,A\ \bm b\,]=\begin{bmatrix}1&-2&3&1\\3&1&-5&-4\\-2&6&-9&-2\end{bmatrix});

まず、1列目を目的の形(単位行列)と一致させたい。

そこで、(1, 1) 要素を軸にして1列目を掃き出す
(この「軸にする」「掃き出す」は専門用語)。

具体的には、(1, 1) 要素が1なので、1行目の3倍を2行目から引き、1行目の2倍を3行目に足せば、(2, 1)要素と(3, 1)要素をゼロにできる。

&math(\sim \begin{bmatrix}1&-2&3&1\\0&7&-14&-7\\0&2&-3&0\end{bmatrix}\begin{array}{l}\bm r_1\\\bm r_2-3\bm r_1\\\bm r_3+2\bm r_1\\\end{array});

行列の右側に出てくる &math(\bm r_1); などは、教科書と同じ書き方で、直前の行列の &math(i); 行目の行ベクトルを &math(\bm r_i); と書いている。あくまでメモ代わりだが、書いておくと計算間違いなどを防ぎやすい。

これで1列目を目的の形にできた。次は2列目。

(2, 2) 要素を軸に掃き出したいが、他の行が7で割り切れないので、まずは2行目に 1/7 を掛けて (2, 2) 要素を1にしておく。

&math(\sim \begin{bmatrix}1&-2&3&1\\0&1&-2&-1\\0&2&-3&0\end{bmatrix}\begin{array}{l}\bm r_1\\\bm r_2/7\\\bm r_3\\\end{array});

そして、(2, 2) 要素を軸に2列目を掃き出す。ここでは1行目に2行目の2倍を、3行目に2行目の−2倍を加えればよい。

&math(\sim \begin{bmatrix}1&0&-1&-1\\0&1&-2&-1\\0&0&1&2\end{bmatrix}\begin{array}{l}\bm r_1+2\bm r_2\\\bm r_2\\\bm r_3-2\bm r_2\\\end{array});

(3, 3) 要素は1なので、そのまま軸にして3列目を掃き出す。1行目に3行目をそのまま、2行目には2倍して加えればよい。

&math(\sim \begin{bmatrix}1&0&0&1\\0&1&0&3\\0&0&1&2\end{bmatrix}\begin{array}{l}\bm r_1+\bm r_3\\\bm r_2+2\bm r_3\\\bm r_3\\\end{array});

この形が求める物であり、

&math(\left\{\begin{array}{c@{\ }c@{\ }c@{\ }c@{\ }c}x&&&=&1\\&y&&=&3\\&&z&=&2\end{array}\right .);

という方程式に対応する。

「式の解を変えない変形」でこの結果を得たのであるから、これは元の式の解でもある。

すなわち、&math(x=1,\ y=3,\ z=2); あるいは、
&math(\begin{bmatrix}x\\y\\z\end{bmatrix}=\begin{bmatrix}1\\3\\2\end{bmatrix}); が与式の解となる。

*** 解答の書き方 [#r3f50ad8]

問:

&math(\left\{\!\begin{array}{r@{\,}c@{\,}r@{\,}c@{\,}r@{\,}c@{\,}r}x&-&2y&+&3z&=&1\\3x&+&y&-&5z&=&-4 \\-2x&+&6y&-&9z&=&-2\end{array}\right .);

を解け。

----

解答:

拡大係数行列に行に対する基本変形を適用する。

&math([\,A\ \bm b\,]=\begin{bmatrix}1&-2&3&1\\3&1&-5&-4\\-2&6&-9&-2\end{bmatrix});

&math(\sim \begin{bmatrix}1&-2&3&1\\0&7&-14&-7\\0&2&-3&0\end{bmatrix}\begin{array}{l}\bm r_1\\\bm r_2-3\bm r_1\\\bm r_3+2\bm r_1\\\end{array});

&math(\sim \begin{bmatrix}1&-2&3&1\\0&1&-2&-1\\0&2&-3&0\end{bmatrix}\begin{array}{l}\bm r_1\\\bm r_2/7\\\bm r_3\\\end{array});

&math(\sim \begin{bmatrix}1&0&-1&-1\\0&1&-2&-1\\0&0&1&2\end{bmatrix}\begin{array}{l}\bm r_1+2\bm r_2\\\bm r_2\\\bm r_3-2\bm r_2\\\end{array});

&math(\sim \begin{bmatrix}1&0&0&1\\0&1&0&3\\0&0&1&2\end{bmatrix}\begin{array}{l}\bm r_1+\bm r_3\\\bm r_2+2\bm r_3\\\bm r_3\\\end{array});

したがって、

&math(\left\{\begin{array}{c@{\ }c@{\ }c@{\ }c@{\ }c}x&&&=&1\\&y&&=&3\\&&z&=&2\end{array}\right .);

より、解は

&math(\begin{bmatrix}x\\y\\z\end{bmatrix}=\begin{bmatrix}1\\3\\2\end{bmatrix});

*** 注意点 [#p7584dcf]

上記では基本変形を表すのに "〜" を使っている。

ここでの "〜" は、同じ解を持つ連立方程式を表す行列であることを表す「同値関係」(これは専門用語なので余力があれば各自調べること)を表わす。
ここでの "〜" は同じ解を持つ連立方程式を表す行列であるという「同値関係」(これは専門用語なので余力があれば各自調べること)を表わす。

基本変形の前後で行列同士が等しいわけではないため "=" で繋いではならない。

** ガウスの消去法:復習 [#h0fcb226]

上で使った方針は、

- 左側の列から順に、対角成分を軸にして掃き出す
- 最後が &math([\,I\ \bm b^\prime\,]); の形になれば解けたことになる

この方法は、任意の大きさの連立一次方程式に適用可能。

しかし、これだけではうまく行かないケースがある。

** うまく行かないケース(1):行の入れ替えが必要 [#kef3e707]

基本変形による計算途中で、以下のように軸に取るべき要素がゼロになってしまう場合がある。

&math([\,A\ \bm b\,] = \begin{bmatrix}2&1&1&2\\4&2&3&1\\-2&-2&0&-1\end{bmatrix});

&math(\sim \begin{bmatrix}2&1&1&2\\0&0&1&-3\\0&-1&1&1\end{bmatrix}\begin{array}{l}\bm r_1\\\bm r_2+2\bm r_1\\\bm r_3+\bm r_1\\\end{array});

次に (2, 2) を軸に掃出しをしたいが、(2, 2) がゼロなので掃出しができない。

ここでは、他の行と入れ替えることで掃出しが行える。(2行目と3行目を入れ替える)

&math(\sim \begin{bmatrix}2&1&1&2\\0&-1&1&1\\0&0&1&-3\end{bmatrix}\begin{array}{l}\bm r_1\\\bm r_3\\\bm r_2\\\end{array});

&math(\sim \begin{bmatrix}2&1&1&2\\0&1&-1&-1\\0&0&1&-3\end{bmatrix}\begin{array}{l}\bm r_1\\-\bm r_2\\\bm r_3\\\end{array});

&math(\sim \begin{bmatrix}2&0&2&3\\0&1&-1&-1\\0&0&1&-3\end{bmatrix}\begin{array}{l}\bm r_1-\bm r_2\\\bm r_2\\\bm r_3\\\end{array});

&math(\sim \begin{bmatrix}2&0&0&9\\0&1&0&-4\\0&0&1&-3\end{bmatrix}\begin{array}{l}\bm r_1-2\bm r_3\\\bm r_2+\bm r_3\\\bm r_3\\\end{array});

&math(\sim \begin{bmatrix}1&0&0&9/2\\0&1&0&-4\\0&0&1&-3\end{bmatrix}\begin{array}{l}\bm r_1/2\\\bm r_2\\\bm r_3\\\end{array});

&math(\therefore \begin{bmatrix}x\\y\\z\end{bmatrix}=\begin{bmatrix}9/2\\-4\\-3\end{bmatrix});

** うまく行かないケース(2):軸に取るべき要素がない [#rce7a4e0]

入れ替えようとしてもできない場合がある。

たとえば、計算途中で次のようになった場合、2列目を掃き出せない~
→ (2, 2) 要素がゼロだし、それより下にも要素がない

&math(\begin{bmatrix}1&2&-2&-3\\0&0&2&2\\0&0&3&3\end{bmatrix});

その場合には、仕方がないのでその列は飛ばして、次列以降の掃き出しを続ける。

&math(\sim \begin{bmatrix}1&2&-2&-3\\0&0&1&1\\0&0&3&3\end{bmatrix}\begin{array}{l}\bm r_1\\\bm r_2/2\\\bm r_3\\\end{array});

&math(\sim \begin{bmatrix}1&2&0&-1\\0&0&1&1\\0&0&0&0\end{bmatrix}\begin{array}{l}\bm r_1+2\bm r_2\\\bm r_2\\\bm r_3-3\bm r_2\\\end{array});

&math(\sim \begin{bmatrix}1&2&0&-1\\0&0&1&1\\0&0&0&0\end{bmatrix}\begin{array}{l}\bm r_1+2\bm r_2\\\bm r_2\\\bm r_3-3\bm r_2\\\end{array});

4列目も掃き出せないのでこれで終了。

- 1,3列目:掃き出せた = ある行に1を残して残りは全て0になっている
- 2,4列目:掃き出せなかった = ごちゃごちゃなまま

これに対応する方程式は、

&math(\left\{ \begin{array}{l}x+2y\phantom{+z}=-1\\\phantom{x+2y+}z=\phantom{-}1\\\phantom{x+2y+}0=\phantom{-}0\end{array} \right .);

3列目の &math(0=0); は、任意の &math(x,y,z); に対して成立するため、
あってもなくても同じであるが、ここでは確認のため明示している。

このような場合、''「掃出しのできなかった列に対応する変数をパラメータに置く」''のが良い~
→ ここでは2列目 = &math(y);

&math(y=a); ただし &math(a); は任意の実数とすると、

&math(\left\{ \begin{array}{l}x=-2a-1\\y=a\\z=1\end{array} \right .);

とすれば全ての式を満たす。すなわち、

&math(\begin{bmatrix}x\\y\\z\end{bmatrix}=a\begin{bmatrix}-2\\1\\0\end{bmatrix}+\begin{bmatrix}-1\\0\\1\end{bmatrix});

が一般解となる。

「一般解」とは?
:一般解|すべての解がこの形に表せることを保証できる形
:特殊解|その値自体は解になるが、他にも解があることを否定しない形
:一般解| すべての解がこの形に表せることを保証できる形
:特殊解| その値自体は解になるが、他にも解があることを否定しない形

*** 行数が足りない場合も同様 [#s9f2a158]

&math(\left\{\begin{array}{l}x+2y-z+3w=-3\\x+y+z+w=1\end{array}\right.);

に対して、1列目、2列目を掃き出した結果は次のようになる。~
→ 掃き出せた列は軸とした列に1があり、残りは0になっている

&math(\sim \begin{bmatrix}1&0&3&-1&5\\0&1&-2&2&-4\end{bmatrix});

3列目、4列目はこれ以上掃き出せない。対応する方程式は、

&math(\left\{ \begin{array}{l}x\phantom{+y}+3z-\phantom{2}w=\phantom{-}5\\\phantom{x+}y-2z+2w=-4\end{array}\right .);

掃出しのできなかった列に対応する変数を &math(z=a,\ w=b); と置けば、

&math(\left\{ \begin{array}{l}x=-3a+b+5\\y=2a-2b-4\\z=a\\w=b\end{array} \right .);

すなわち、

&math(\begin{bmatrix}x\\y\\z\\w\end{bmatrix}=a\begin{bmatrix}-3\\2\\1\\0\end{bmatrix}+b\begin{bmatrix}1\\-2\\0\\1\end{bmatrix}+\begin{bmatrix}5\\-4\\0\\0\end{bmatrix});

が一般解となる。

** うまく行かないケース(3):定数項まで掃き出せてしまう [#rce7a4e0]

&math(\begin{bmatrix}1&2&-2&-3\\0&0&2&2\\0&0&3&4\end{bmatrix});

について掃出しを進めると、

&math(\sim \begin{bmatrix}1&2&-2&-3\\0&0&1&1\\0&0&3&4\end{bmatrix}\begin{array}{l}\bm r_1\\\bm r_2/2\\\bm r_3\end{array});

&math(\sim \begin{bmatrix}1&2&0&-1\\0&0&1&1\\0&0&0&1\end{bmatrix}\begin{array}{l}\bm r_1+2\bm r_2\\\bm r_2\\\bm r_3-3\bm r_2\end{array});

ここまででおかしな事になっているのだけれど、形式的には4列目も掃き出せてしまう。

&math(\sim \begin{bmatrix}1&2&0&0\\0&0&1&0\\0&0&0&1\end{bmatrix}\begin{array}{l}\bm r_1+\bm r_3\\\bm r_2-\bm r_3\\\bm r_3\end{array});

対応する方程式は、

&math(\left\{\begin{array}{l}x+2y\phantom{+z}=0\\\phantom{x+2y+}z=0\\\phantom{x+2y+}0=1\end{array} \right .);

この連立方程式(特に3番目の式)は、どんな &math(x,\ y,\ z); に対しても成り立たないため、
「解無し」が結論になる。

本来なら、係数部分がすべてゼロで、定数部分が非ゼロの行ベクトルが現れた時点で「解無し」を結論できる。

** 連立一次方程式の解法 [#h5f54f81]

+ 掃出しにより &math([\,A\ \bm b\,]\sim[\,I\ \bm b^\prime\,]); とできれば、解は &math(\bm x=\bm b^\prime);
+ 掃き出せない列があれば、飛ばして計算を続け、最後にその列をパラメータに置く
+ 定数項(&math(\ne 0);)を残して他がゼロの行が現れたら「解無し」

により、どんな連立一次方程式も解ける(任意の次数でOK)

注)「解ける」=「一般解が求まる」 or 「解無しを結論する」

注2)繰り返しになるが、この方法は必ずしも「効率的な方法」を求めた物ではないが、
注2)口を酸っぱく繰り返しになるが、この方法は必ずしも「効率的な方法」を求めた物ではないが、
「どんな連立方程式もこの方法で必ず解ける」という事実が後で重要になる。
後のためにもガウスの消去法は完璧にマスターしておくように。

* 逆行列を求める (§2.5 を一部先取り) [#hb26fec8]

** 準備 [#reb713f0]

上記のように &math(A\bm x=\bm b); を &math(\bm x); について解けるようになれば、
与えられた正方行列 &math(A); について逆行列 &math(A^{-1}); を求める事もできる。
任意の正方行列 &math(A); について逆行列 &math(A^{-1}); を求める事もできる。

というのも、

&math(
A\bm x_1=\bm e_1=\begin{bmatrix}1\\0\\\vdots\\0\end{bmatrix},\ \ 
A\bm x_2=\bm e_2=\begin{bmatrix}0\\1\\\vdots\\0\end{bmatrix},\ \ 
\dots,\ \ 
A\bm x_n=\bm e_n=\begin{bmatrix}0\\0\\\vdots\\1\end{bmatrix}
);

となるような &math(\bm x_1,\bm x_2,\dots,\bm x_n); を求め、
それを並べた行列

&math(X=\begin{bmatrix}\bm x_1&\bm x_2&\dots&\bm x_n\end{bmatrix});

を考えると、

&math(AX&=\Big[\ \ A\ \ \Big] \begin{bmatrix}\bm x_1&\bm x_2&\dots&\bm x_n\end{bmatrix}\\
&=\begin{bmatrix}A\bm x_1&A\bm x_2&\dots&A\bm x_n\end{bmatrix}\\
&=\begin{bmatrix}\bm e_1&\bm e_2&\dots&\bm e_n\end{bmatrix}\\
&=I);

となる。すなわち、&math(AX=I); より

&math(X=A^{-1});

である(かもしれない)。

→ 上記 &math(n); 個の連立方程式の解を並べて書けば、それが逆行列である!(かもしれない)

***「かもしれない」の意味は・・・ [#y6fef9ea]

逆行列の定義は &math(AX=XA=I); であったのに、&math(XA=I); を確かめていないから、そこが問題。

&math(A); が正方行列であれば、&math(AX=I); から &math(XA=I); を導けることは後に証明する。

** 求め方 [#td64e8ed]

&math(A\bm x_k=\bm e_k); (&math(k=1\dots n);) という &math(n); 個の連立方程式を全て解くのは大変だが、例えば &math(A\bm x_1=\bm e_1); を解く手順と &math(A\bm x_2=\bm e_2); を解く手順を並べて書けば

&math(
&\begin{bmatrix}
&1\ \\
&0\ \\
\hspace{1cm}A\hspace{1cm}&0\ \\
&\vdots\ \\
&0\ \\
\end{bmatrix}\\
&\sim
\begin{bmatrix}
&\phantom{1\ }\\
&\\
\hspace{1cm}\phantom{A}\hspace{1cm}&\phantom{0}\ \\
&\\
&\\
\end{bmatrix}\\
&\sim
\begin{bmatrix}
&\phantom{1\ }\\
&\\
\hspace{1cm}\phantom{A}\hspace{1cm}&\phantom{0}\ \\
&\\
&\\
\end{bmatrix}\\
&\hspace{1cm}\vdots
);   
&math(
&\begin{bmatrix}
&0\ \\
&1\ \\
\hspace{1cm}A\hspace{1cm}&0\ \\
&\vdots\ \\
&0\ \\
\end{bmatrix}\\
&\sim
\begin{bmatrix}
&\phantom{1\ }\\
&\\
\hspace{1cm}\phantom{A}\hspace{1cm}&\phantom{0}\ \\
&\\
&\\
\end{bmatrix}\\
&\sim
\begin{bmatrix}
&\phantom{1\ }\\
&\\
\hspace{1cm}\phantom{A}\hspace{1cm}&\phantom{0}\ \\
&\\
&\\
\end{bmatrix}\\
&\hspace{1cm}\vdots
);

となるが、「係数行列の等しい方程式の掃き出し手順は定数項によらず等しくなる」ことに気付く。~
← 1列目を掃き出す手順は &math(A); の1列目で決まる~
← 2列目を掃き出す手順は &math(A); の2列目で決まる~
…~
← &math(n); 列目(定数列の直前)を掃き出す手順は &math(A); の &math(n); 列目で決まる~
← ← &math(\bm e_k); の値に依らない!

これを利用して、次のように &math(n); 個の方程式を一度に解くことができる。

例 = 教科書P47 2.2(2):

&math(A=\begin{bmatrix}1&2&3\\2&4&5\\3&5&6\end{bmatrix}); の逆行列を求める。

&math(A); の右に単位行列 &math(I); を並べた行列に掃き出し法を適用する。

&math(\Bigg[\ A\ \ \ I\ \ \Bigg]=&\begin{bmatrix}
1&2&3&1&0&0\\
2&4&5&0&1&0\\
3&5&6&0&0&1\\
\end{bmatrix}\begin{array}{lll}\bm r_1\\\bm r_2\\\bm r_3\end{array}\\
\sim&\begin{bmatrix}
1&2&3&1&0&0\\
0&0&-1&-2&1&0\\
0&-1&-3&-3&0&1\\
\end{bmatrix}\begin{array}{lll}\bm r_1\\\bm r_2-2\bm r_1=\bm r_4\\\bm r_3-3\bm r_1=\bm r_5\end{array}\\
\sim&\begin{bmatrix}
1&2&3&1&0&0\\
0&-1&-3&-3&0&1\\
0&0&-1&-2&1&0\\
\end{bmatrix}\begin{array}{lll}\bm r_1\\\bm r_5\\\bm r_4\end{array}\\
\sim&\begin{bmatrix}
1&0&-3&-5&0&2\\
0&-1&-3&-3&0&1\\
0&0&-1&-2&1&0\\
\end{bmatrix}\begin{array}{lll}\bm r_1+2\bm r_5=\bm r_6\\\bm r_5\\\bm r_4\end{array}\\
\sim&\begin{bmatrix}
1&0&0&1&-3&2\\
0&-1&0&3&-3&1\\
0&0&-1&-2&1&0\\
\end{bmatrix}\begin{array}{lll}\bm r_6-3\bm r_4=\bm r_7\\\bm r_5-3\bm r_4=\bm r_8\\\bm r_4\end{array}\\
\sim&\begin{bmatrix}
1&0&0&1&-3&2\\
0&1&0&-3&3&-1\\
0&0&1&2&-1&0\\
\end{bmatrix}\begin{array}{lll}\bm r_7\\-\bm r_8\\-\bm r_4\end{array}\\
&\hspace{2cm}\underbrace{\smash{\hspace{2cm}}}_{\displaystyle A^{-1}}
);

&math(
\therefore A^{-1}=\begin{bmatrix}
1&-3&2\\
-3&3&-1\\
2&-1&0\\
\end{bmatrix}
);

上記手順は、例えば 5, 6 列目を隠してみると &math(\Big[\ A\ \ \bm e_1\ \Big]); 
を掃き出す手順と全く同じであることが分かる。

&math(\Big[\ A\ \ \bm e_1\ \Big]); を掃き出したとき、右端に残るのが &math(\bm x_1); 
であること

&math(
\Big[\ A\ \ \bm e_1\ \Big]\sim
\Big[\ I\ \ \bm x_1\ \Big]
);

と比べれば、

&math(
\Bigg[\ A\ \ \bm e_1\ \bm e_2\ \dots\ \bm e_n\ \Bigg]
); を掃き出した結果は

&math(
\sim\Bigg[\ I\ \ \bm x_1\ \bm x_2\ \dots\ \bm x_n\ \Bigg]=\Big[\ I\ \ X\ \Big]=\Big[\ I\ \ A^{-1}\ \Big]
);

であることが分かる。

例 = 教科書P47 2.2(1):

&math(
A=\begin{bmatrix}
7&8\\
8&9
\end{bmatrix}
);
の逆行列を求める。

&math(
\Big[A\ I\,\Big]=&\begin{bmatrix}
7&8&1&0\\
8&9&0&1
\end{bmatrix}
\begin{array}{l}\bm r_1\\\bm r_2\end{array}\\
\sim&\begin{bmatrix}
-1&-1&1&-1\\
8&9&0&1
\end{bmatrix}
\begin{array}{l}\bm r_1-\bm r_2=\bm r_3\\\bm r_2\end{array}\\
\sim&\begin{bmatrix}
-1&-1&1&-1\\
0&1&8&-7
\end{bmatrix}
\begin{array}{l}\bm r_3\\\bm r_2+8\bm r_3=\bm r_4\end{array}\\
\sim&\begin{bmatrix}
-1&0&9&-8\\
0&1&8&-7
\end{bmatrix}
\begin{array}{l}\bm r_3+\bm r_4=\bm r_5\\\bm r_4\end{array}\\
\sim&\begin{bmatrix}
1&0&-9&8\\
0&1&8&-7
\end{bmatrix}
\begin{array}{l}-\bm r_5\\\bm r_4\end{array}\\
);

高校で習った公式を用いて検算してみる。

&math(
A^{-1}=&\frac{1}{ad-bc}\begin{bmatrix}d&-b\\-c&a\end{bmatrix}\\
A^{-1}=\begin{bmatrix}a&b\\c&d\end{bmatrix}^{-1}=&\frac{1}{ad-bc}\begin{bmatrix}d&-b\\-c&a\end{bmatrix}\\
=&\frac{1}{7\times 9-8\times 8}\begin{bmatrix}9&-8\\-8&7\end{bmatrix}\\
=&-1\times \begin{bmatrix}9&-8\\-8&7\end{bmatrix}\\
=&\begin{bmatrix}-9&8\\8&-7\end{bmatrix}\\
);

** まとめ [#q0e3e0e5]

逆行列の求め方:

+ &math(\Big[\,A\ I\ \Big]); を作って左から掃き出していく
+ &math(\Big[\,I\ X\ \Big]); の形になれば、&math(X=A^{-1}); である。
+ 途中で掃き出しができなくなったら、&math(A); は正則でない(逆行列は存在しない)

いろいろ証明せずに進んでいる所があるが、次回以降「行列の積」を用いて上記方法を理解し、
証明できるようになる。

* コメント [#q5588fba]

#article_kcaptcha
**無題 [#xe577618]
>[[Nice]] (&timetag(2015-07-23T02:42:01+09:00, 2015-07-23 (木) 11:42:01);)~
~
Nicework~

//

#comment_kcaptcha


Counter: 366995 (from 2010/06/03), today: 13, yesterday: 0