連立方程式の解法〜行列でやると

余因子を使う逆行列

１．余因子と逆行列

連立方程式の問題は、小学生なら消去算として、中高では連立方程式として、教科書や参考書でさんざんやってきているでしょう。この回では、 連立方程式の解法を一般化することにチャレンジしよう。いろいろなやり方がある。連立方程式を解くときの係数に着目して 行列を使って解の定式化してみよう。未知数が５個、列ベクトルにしてｘ=

に対して、係数が5✕5の要素の行列（５個の行ベクトルa1,a2,a3,a4,a5） A=

係数行列Aと列ベクトルｘの積を列ベクトル b=

とする。 つまり、連立方程式は Ax=b ｘ＝＝とも、かけるね。 Aのi行目だけを取り出すと、a_iｘ＝b_iとなる。 この内積の式を全部で5行並べて書き出すと、普通の連立方程式になるね。 ＜余因子行列式法＞ Aの逆行列が出せると、ｘが出せるはずだ。 Ax=bの両辺に左からA^-1をかけると、ｘ＝A^-1ｂとなる。だから、連立方程式を解くことは、Aの逆行列の求め方の問題に置き換えることもできるね。 x=A^-1b　　　　　j列 A^-1=1/det(A)｛…　…　… 　　　　　i行 … (Aji) … 　　　　　　 …　…　…｝＝1/det(A ) t {{A11,A12,A13,A14,A15}, {A21,A22,A23,A24,A25}, {A31,A32,A33,A34,A35}, {A41,A42,A43,A44,A45}, {A51,A52,A53,A54,A55}} A^-1のi行j列＝ 1/det(A) (Aji) （★転置 [ transpose]していることに注意しよう） Aijは、Aのi行j列の成分aijに対する余因子行列式（Aのi行とj列をぬいた残り要素の行列式に符号つけ）たとえば、3次のときなら、 A={{a11, a12, a13}, {a21, a22, a23}, {a31, a32, a33}} b={b1,b2,b3} det(A)= a11a22a33+a12a23a31+a13a21a32 -( a11a23a32+a12a21a33+a13a22a31) A^-1=

　 =1/det(A) ｔ{{+(a22a33-a23a32), -(a21a22-a23a31), +(a21a32-a22a31)}, { -(a12a33-a13a32), +(a11a33-a13a31), -(a11a32-a12a31)}, {+(a12a23-a13a22), -(a11a23-a13a21), +(a11a22-a12a21)}} （例） A={{０, −１, １}, {１, ２, −１}, {２, １, ２}}の逆行列は？ det(A)= (0+2+1) -(4-2+0)=1 A^-1=

　　　t {{+(2・2-1・(-1)), -(1・2-2・(-1)), +(1・2-2・2)} { -(−1・2-1・1), +(0−1・2), -(0ー2・(-1))}, {+((-1)(-1)-2・1), -(0-1・1) , +(0ー１(-1))}} =

クラメールの公式

２．クラメールの公式

＜列ごとおきかえ行列式法＞ 逆行列を出さずに、直接、解ｘを計算する解法もある。 クラメールの公式とも言われているね。 xi =|Ai|/det(A)　　　（ |Ai|はAのi列を列ベクトルｂにおきかえた行列式。） ・3次のとき、 A={{a11, a12, a13}, {a21, a22, a23}, {a31, a32, a33}} b=

det(A) = a11a22a33+a12a23a31+a13a21a32 -( a11a23a32+a12a21a33+a13a22a31) det(A1)= b1a22a33 +a12a23b3 +a13b2a32 -( b1a23a32 +a12b2a33 +a13a22b3) det(A2)= a11b2a33 +b1a23a31 +a13a21b3 -( a11a23b3 +b1a21a33 +a13 b2a31) det(A3)= a11a22b3 +a12b2a31 +b1a21a32 -( a11b2a32 +a12a21b3 +b1a22a31) x ==

（例） x -2y +z =-6 2x+y +3z =-5 x -y -2z =5 をクラメールの方式で解くと？ A= {{１, −２, １}, {２, １, ３}, {１, −１, −２}} b=

det(A) = (-2)+(-6)+(-2) - (1+8+(-3)) = -10-6=-16 det(A1)= (12)+(-30)+(5) - ( 5+(-20)+18)= -13-3=-16 det(A2)= (10+(-18)+10) -(-5+24+15)=2-34=-32 det(A3)= (5+10+12)- (-6+(-20)+5)=27-(-21)=48

３．掃き出し法

余因子やクラメールは２次・３次くらいなら扱いやすい。しかし、次数が多くなると扱いにくいし、解が不定のときは使えない。それに対して、行列と対話しながら、行列を変形を繰り返すことで、段階的に、解や逆行列を求める方法がある。行数、列数が多くても実行できるという利点、次数よりも情報が少ない場合でも、解を表現する手がかりが得られる利点がある。 Ax=bの係数行列Aと値列ベクトルbを横に並べて合体した[A|b]を拡大係数行列と呼ぼう。 [A|b]が[E|u]に変形できたときに、ｘ＝uが解になる。つまり、 [A|E]が[E|A^-1]に変形できたときに,　逆行列A^-1が出る。 ・Aの行基本変形は３種ある。　i行とj行入れ替え (A_i↔j) 　i行の定数c倍 (A_i(c)) 　j行の定数c倍を, j行にたす ( A_i←j(c)) Aの部分の対角成分を１にそれ以外を０に掃き出すイメージで 掃き出し法[Sweeping method]。（例） x -2y +z =-6 2x+y +3z =-5 x -y -2z =5 掃き出し法で解くと？ [A|b]= {{１, -２, １ , - ６}, {２,１, ３ , - ５}, {１, -１, -２, ５}} １行目の-２倍を２行目にたす、１行目を３行目からひく。 {{１, -２, １ , - ６}, {０,５, １ , ７}, {０, １, -３, １１}} ３行目の２倍を１行目にたし、2行目のー０.2倍を３行目をたす。 {{１, ０, -5 , １６}} {０,５, １ , ７}, {０, ０, -3.2, 9.6}} ３行目の1/(-3.2)倍。 {{１,０, -5 , 16}, {０,5, １ , 7}, {０,０,１,-3}} ３行目の５倍を１行目にたし、３行目の−１倍を２行目からひく。 {{１,０, ０, １ }, {０,５,０ , 10 }, {０,０,１,-３ }} ２行目を0.2倍する。 {{１,０, ０, １ }, {０,１,０ , ２ }, {０,０,１,-３ }} =[E|x] （ｘ，ｙ，ｚ）＝（１，２，３） （例） {{-３, ２, ２}, {-５,４, ３}, {１, -１, -1}}の逆行列は？ [A|E]= {{-３, ２, ２, 1, 0, 0}, {-５,４, ３, 0, 1, 0}, {１, -１, -1, 0, 0, 1}} １行目の-5/3倍を２行目にたす、１行目の1/3倍を３行目にたす。 1行目を-3で割る。 {{1, -2/3, -2/3, -1/3, 0, 0}, {０,2/3, -1/3, -5/3, 1, 0}, {０, -1/3,-1/3, 1/3, 0, 1}} 2行目を１行目にたし、2行目の1/2倍を３行目をたす。２行目を3/2倍する。 {{1, ０, -１, -２, １, ０}, {０,1, -1/2, -5/2,3/2,０}, {０, 0,-1/2, -1/2, 1/2, １}} ３行目の-２倍を１行目にたし、３行目の−１倍を２行目にたす。３行目を−２倍する。 {{１, ０, ０, -１,０,-２}, {０,１, ０, -２, １,-１}, {０, ０,１,１, -１, -２}} =[E|A^-1] Aの逆行列は、 {{ -１, ０, -２}, { -２, １, -１}, { １,-１, -２}}

原点通過と自由度

４．原点を通る３平面の交わり

連立方程式を解くためには、拡大係数行列のAの部分が単位行列になるまで行の基本変形を続ける必要があったね。解を求めるためではなく、解があるかどうか、あるとしたら１つか不定かを知りたいときがあるね。その場合は、上の行から順に０でない係数が１列目、２列目、３列目と階段のようにならぶ形になるように変形をしていけばよいね。 行の基本変形をして、０になる行を作っていくとき、０以外のある行数を階数[rank]という。言い換えると、もとの行数から、０だけの行ベクトルの行数を引いたものが階数だね。そのような変形で簡約化した行列のことを階段行列[step matrix]と呼ぶことがもあるね。そのときの各行で０でない最初の列の数を主成分といい、主成分が１になるまで基本変形した場合、その行列を簡約な行列と呼ぶこともある。ｎ次の正方行列Aで、rankA=nならAは正則。 rankがｎ未満では正則でない。 n - rankAを自由度[degree of freedom]という。　 ＜同次連立方程式＞ 拡大係数行列の最後の列が０ベクトルになる連立方程式を同次とか斉次とか呼んだりする Ax=0　のこと。３D空間では、３つの行は、３つの原点通過平面に対応しているね。 自由度＝０だと解は原点ｘ＝(x,y,z)=(0,0,0)の自明な、つまらないものだけになる。自由度が正なら自明でない解がある。自由度１は直線、自由度２は平面に対応するね。 （例）自由度０ A= {{１,−１,１,０}, {３,１,１,０}, {２,２,３,０}} ⇒ {{１,−１,１,０}, {０,４,−２,０}, {０,０, ３,０}} ⇒ n=3, rankA=3 自由度 n - rankA = 0 ０パラメータ ｘ＝（０，０，０）のみ。 3D空間では解は原点Oのみ。だから、３行の３平面を表し、３平面が原点で交わっている関係だね。 （例）自由度１ A= {{１,１,-１,０}, {４,３,-１,０}, {６,５,-３,０}} ⇒ {{１,１,-１,０}, {０,-１,３,０}, {０,０,０,０}} ⇒ n=3, rankA=2 自由度 n - rankA = 1 １パラメータp ｘ＝（-2p，3p，p）のみ。３D空間では解は原点通過の直線。だから、３行の３平面を表し、３平面が原点を通る直線で交わっている関係だね。 （例）自由度２ A= {{３,２,-１,０}, {６,４,-２,０}, {-12,-８,４,０}} ⇒ {{３,２,-１,０}, {０,０,０,０}, {０,０,０,０}} ⇒ n=3, rankA=1 自由度 n - rankA = 2 2パラメータp,q ｘ＝（p，q，3p+2q）のみ。３D空間では解は原点通過の平面。だから、３行が見かけ上は３平面を表し、実際は１つの原点を通る平面になる関係だね。自由度が２だと、自由すぎてつまらない。３平面が１つに重なっているからね。

一般の平面の交わり

５．一般の平面の交わり

＜非同次連立方程式＞ 拡大係数行列の最後の列が０ベクトルでない連立方程式を非同次とか非斉次とか呼んだりする。 Ax=b ３D空間では、各行は原点を通るとは限らない平面に対応している。 係数行列Aのランクr=rank(A)は平面の異なる法線ベクトルの数に対応し、 拡大係数行列A~=[A|b]のランクr~=rank(A~)は平面の法線と通過点でみた事実上の平面数に対応する。だから、係数行列のランクrより、拡大係数行列のランクr~が大きいと法線ベクトルが同じで通過点が異なる平面、つまり、平行面が混入する。だから、共有空間はないので、解なし。拡大したものとランクが同じなら平行面はないので、解がある。さらに、自由度＝n - r=０なら、ただ１組の解をもつ。自由度が正なら不定解。自由度は同次連立方程式と同じで、解空間の次元（直線、平面、。。。）を表すね。 （例）自由度０ A= {{１,２,-２,２}, {２,５,-１,８}, {-２,-１,11,18}} ⇒ {{１,２,-２,２}, {０,１,３,４}, {０,０, １,-５}} ⇒ 拡大前Aも後A~もランクは３。 n=3, rankA=3 自由度 n - rankA = 0 ０パラメータ ｘ＝（-46，19，-5）のみ。 3D空間では点。 （例）自由度１ A= {{１,３,-２,１}, {３,10,-３,４}} ⇒ {{１,３,-２,１}, {０,１,３,１}} ⇒ 拡大前Aも後A~もランクは２。 n=3, r=rankA=2 自由度 n - r= 1 1パラメータp ｘ＝（11p-2，-3p+1，p）のみ。３D空間では原点を通らない直線。 （例）ランクの不一致 A= {{１,２,-３,２}, {３,６,-９,１０}} ⇒ {{１,２,-３,２}, {０,０,０,４}} ⇒ 係数のランクr=１で、拡大ランクr~=２解なし。 ３D空間では、２つの平面の法線ベクトルが同じで交点なし。平行な２平面は交わらない。

一般の直線の交わり

６．演習

・解が１つに決まる連立方程式をできるだけいろいろな方法で解こう。　3x +y +2z= 9 x -2y +z = 8 -2x - y +z= -3 A={{3, 1, 2 }, b={{9}, {1, -2, 1}, {8}, {-2, -1,1}} {-3}} （掃き出し法） 拡大係数行列は[A|b] A={{3, 1, 2, 9 }, {1, -2, 1, 8}, {-2,-1, 1, -3}} 行交換A_2↔1 A={{1, -2, 1, 8 }, 　{3, 1, 2, 9 }, 　 {-2,-1, 1, -3}} A_2←1(-3), A_3←1(2) A={{1, -2, 1, 8 }, 　{0, 7, -1, -15 }, 　 {0,-5, 3, 13}} A_2←3, A_2(1/2)A={{1, -2, 1, 8 }, 　{0, 1, 1, -1 }, 　 {0,-5, 3, 13}} A₃←2(5), A_3(1/8)A={{1, -2, 1, 8 }, 　{0, 1, 1, -1 }, 　 {0, 0 1, 1}} A_{1←2(2),
A1←3(-3)}A={{1, 0, 0, 3 }, 　{0, 1, 1, -1 }, 　 {0, 0 1, 1}} A2←3(-1) A={{1, 0, 0, 3 }, 　{0, 1, 0, -2 }, 　 {0, 0 1, 1}} (x,y,z)=(3,-2,1) 掃き出し法のいいところは、公式を覚えなくてもよいこと。ただし、基本変形を大量にするため、そのことによる計算ミスをしやすい。 （余因子） A={{3, 1, 2 }, {1, -2, 1}, {-2, -1,1}} det(A)=3･(-2)･1+ 1･(-1)･2+(-2)･1･1 -{(-2)･(-2)･2+1･1･1+3･(-1)･1}=-6-2-2-{8+1-3}=-10-6=-16 A11=det({-2,-1},{-1,1}}=(-2+1)=-1。A12=-det({1,1},{-2,1})=-(1+2)=-3, A13=det({1,-2},{-2,-1})=(-1-4)=-5, A21=-det({1,2},{-1,1})=-(1+2)=-3, A22=det({3,2},{-2,1})=(3+4)=7, A23=-det({3,1},{-2,-1})=-(-3+2)=1, A31=det({1,2},{-2,1})=(1+4)=5, A32=-det({3,2},{1,1})=-(3-2)=-1, A33=det({3,1},{1,-2})= -6-1=-7 A-1=1/-16tr{{-1, -3, -5}, {-3, 7, 1}, {5, -1, -7}} -1/16{{-1, -3, 5}, {-3, 7, -1}, {-5, 1, -7}} (x,y,z)= A-1{{9},{8},{-3}}=-1/16{ (-1)･9+(-3)･8+5･(-3), (-3)･9+7･8+(-1)･(-3), (-5)･9+1･8+(-7)･(-3)} =-1/16{ -9-24-15, -27+56+3, -45+8+21} =-1/16{-48, 32, -16} 　＝(3, -2, 1) 余因子法のいいところは、行列式の個数が増えても、次数が下がること。ー符号を１つおきに交互につけるのを忘れたり、転置するのを忘れるとアウト。 （クラメール） det(A) =3･(-2)･1+ 1･(-1)･2+(-2)･1･1 -{(-2)･(-2)･2+1･1･1+3･(-1)･1}=-6-2-2-{8+1-3}=-10-6=-16 det(A1)=9･(-2)･1+ 8･(-1)･2+(-3)･1･1 -{(-3)･(-2)･2+8･1･1+9･(-1)･1}=-18-16-3-{12+8-9}=-37-11=-48 det(A) =3･(-2)･1+ 1･(-1)･2+(-2)･1･1 -{(-2)･(-2)･2+1･1･1+3･(-1)･1} det(A2) =3･8･1+ 1･(-3)･2+(-2)･9･1 -{(-2)･8･2+1･9･1+3･(-3)･1}=24-6-18-{-32+9-9}=0+32=32 det(A) =3･(-2)･1+ 1･(-1)･2+(-2)･1･1 -{(-2)･(-2)･2+1･1･1+3･(-1)･1}= det(A3) =3･(-2)･(-3)+ 1･(-1)･9+(-2)･1･8 -{(-2)･(-2)･9+1･1･(-3)+3･(-1)･8}=18-9-16-{36-3-24}=-7-9=-16 (x,y,z)=1/-16 { -48, 32, -16}= (3, -2, 1) クラメールのいいところは、行列式の個数が未知数の数だけしかないこと。ただし、ｂのAへの代入で行列そのものを変えてしまったもので行列式を計算するためもとの行列と変わったところの確認が大切。