連立方程式の解法〜行列でやると
連立方程式の問題は、[br]小学生なら消去算として、[br]中高では連立方程式として、[br]教科書や参考書でさんざんやってきているでしょう。[br][br]この回では、[br][b][color=#0000ff][size=150]連立方程式の解法を[br]一般化すること[/size][/color][/b]にチャレンジしよう。[br]いろいろなやり方がある。[br][br]連立方程式を解くときの[color=#0000ff][b]係数に着目して[/b][/color][br][b][color=#0000ff]行列を使って[/color][/b]解の定式化してみよう。[br][br]未知数が5個、[color=#0000ff][b][size=150]列ベクトル[/size][/b][/color]にして[br][b]x[/b]=[math]\left(\begin{matrix}x1\\x2\\x3\\x4\\x5\\\end{matrix}\right)[/math]に対して、[br]係数が5✕5の要素の行列([color=#0000ff][b]5個の行ベクトル[/b][/color][b]a1,a2,a3,a4,a5[/b])[br][b][i]A=[/i][/b][math]\left(\begin{matrix}a1\\a2\\a3\\a4\\a5\\\end{matrix}\right)[/math][br]係数行列[b][i]A[/i][/b]と列ベクトル[b]x[/b]の[color=#0000ff][b]積を列ベクトル[/b][/color][br][b]b[/b]=[math]\left(\begin{matrix}\begin{matrix}b1\\b2\\b3\\b4\\b5\\\end{matrix}\end{matrix}\right)[/math]とする。[br][br][size=200][b][color=#0000ff][i]つまり、連立方程式は[br]Ax=b[br][/i][math]\left(\begin{matrix}a1\\a2\\a3\\a4\\a5\end{matrix}\right)[/math]x=[math]\left(\begin{matrix}a1x\\a2x\\a3x\\a4x\\a5x\end{matrix}\right)[/math]=[math]\left(\begin{matrix}\begin{matrix}b1\\b2\\b3\\b4\\b5\\\end{matrix}\end{matrix}\right)[/math][/color][/b][/size]とも、かけるね。[br][br][size=150][b]Aのi行目だけを取り出すと、[size=200][color=#0000ff][i]a[/i][sub]i[/sub]x=b[sub]i[/sub][/color][/size]となる。[br][/b][size=100]この内積の式を全部で5行並べて書き出すと、普通の連立方程式になるね。[br][/size][/size][br][size=150][b]<余因子行列式法>[br][/b][size=100]Aの逆行列が出せると、xが出せるはずだ。[br][/size][b]Ax=bの両辺に左からA[/b][sup]-1[/sup][b]をかけると、x=A[/b][sup]-1[/sup][b]bとなる。[br][/b][size=100]だから、連立方程式を解くことは、Aの逆行列の求め方の問題に置き換えることもできるね。[br][/size][/size][i][size=150][b][color=#0000ff]x=A[sup]-1[/sup]b j列[br]A[sup]-1[/sup]=1/det(A){… … …[br] i行 … (Aji) …[br] … … …}[br]=1/det(A ) [/color][color=#ff0000]t [/color][color=#0000ff]{{A11,A12,A13,A14,A15},[br] {A21,A22,A23,A24,A25},[br] {A31,A32,A33,A34,A35},[br] {A41,A42,A43,A44,A45},[br] {A51,A52,A53,A54,A55}}[br]A[/color][sup]-1[/sup][color=#0000ff]の[u]i行j列[/u]= [/color][i][size=150][color=#0000ff]1/det(A) (A[u]ji[/u]) ([/color][color=#ff0000]★転置 [ transpose][/color][color=#0000ff]していることに注意しよう)[br][/color][/size][/i][/b][/size][/i][size=150]Aij[/size]は、Aのi行j列の成分aijに対する[b][size=150]余因子行列式[/size][/b](Aのi行とj列をぬいた残り要素の行列式に符号つけ)[br][color=#0000ff][br]たとえば、[/color]3次のときなら、[br][b]A[/b]={{a11, a12, a13},[br] {a21, a22, a23},[br] {a31, a32, a33}}[br][b]b[/b]={b1,b2,b3}[br]det([b]A[/b])=[size=150] a[size=100]11a22a33+a12a23a31+a13a21a32 -( a11a23a32+a12a21a33+a13a22a31) [br][/size][/size][b]A[sup]-1[/sup][/b]=[math]\frac{1}{det\left(A\right)}^t\left(\begin{matrix}A11,A12,A13\\A21,A22,A23\\A31,A32,A33\end{matrix}\right)[/math] [br]=1/det([b]A[/b]) [color=#0000ff][b]t[/b][/color]{{+(a22a33-a23a32), -(a21a22-a23a31), +(a21a32-a22a31)},[br] { -(a12a33-a13a32), +(a11a33-a13a31), -(a11a32-a12a31)},[br] {+(a12a23-a13a22), -(a11a23-a13a21), +(a11a22-a12a21)}}[br][size=150][color=#0000ff](例)[br][/color][size=100][color=#0000ff]A={{0, −1, 1},[br] {1, 2, −1},[br] {2, 1, 2}}の逆行列は?[br][/color]det(A)=[size=150] (0+2+1) -(4-2+0)=1 [br][/size][b]A[sup]-1[/sup][/b]=[math]\frac{1}{det\left(A\right)}^t\left(\begin{matrix}A11,A12,A13\\A21,A22,A23\\A31,A32,A33\end{matrix}\right)[/math] [br] t {{+(2・2-1・(-1)), -(1・2-2・(-1)), +(1・2-2・2)}[br] { -(−1・2-1・1), +(0−1・2), -(0ー2・(-1))},[br] {+((-1)(-1)-2・1), -(0-1・1) , +(0ー1(-1))}}[br]=[math]t\left(\begin{matrix}5,-4,-3\\3,-2,-2\\-1,1,1\end{matrix}\right)[/math][br][br]=[math]\left(\begin{matrix}5,3,-1\\-4,-2,1\\-3,-2,1\end{matrix}\right)[/math] [/size][/size]
[size=150][b]<列ごとおきかえ行列式法>[br][/b]逆行列を出さずに、直接、解xを計算する解法もある。[b][br]クラメールの公式[/b]とも言われているね。[br][/size][color=#0000ff][b][size=150]xi =|Ai|/det(A) ( |Ai|はAのi列を列ベクトルbにおきかえた行列式。)[br][/size][/b][/color][color=#0000ff]・[/color]3次のとき、[br][b]A[/b]={{a11, a12, a13},[br] {a21, a22, a23},[br] {a31, a32, a33}}[br][b]b[/b]=[math]\left(\begin{matrix}b1\\b2\\b3\end{matrix}\right)[/math][br]det([b]A[/b]) =[size=150] a[size=100]11a22a33+a12a23a31+a13a21a32 -( a11a23a32+a12a21a33+a13a22a31) [br]det([b]A1[/b])=[size=150] b1[size=100]a22a33 +a12a23b3 +a13b2a32 -( b1a23a32 +a12b2a33 +a13a22b3) [br][/size][/size]det([b]A2[/b])=[size=150] a[size=100]11b2a33 +b1a23a31 +a13a21b3 -( a11a23b3 +b1a21a33 +a13 b2a31) [br][/size][/size]det([b]A3[/b])=[size=150] a[size=100]11a22b3 +a12b2a31 +b1a21a32 -( a11b2a32 +a12a21b3 +b1a22a31) [br][/size][/size][b]x [b]=[math]\left(\begin{matrix}x1\\x2\\x3\\\end{matrix}\right)[/math][/b][/b]= [math]\frac{1}{det\left(A\right)}\left(\begin{matrix}det\left(A1\right)\\det\left(A2\right)\\det\left(A3\right)\end{matrix}\right)[/math][br][/size][/size][color=#0000ff](例)[/color][br][color=#0000ff] x -2y +z =-6[br]2x+y +3z =-5[br] x -y -2z =5[br]をクラメールの方式で解くと?[br][/color][b]A[/b]=[br]{{1, −2, 1},[br] {2, 1, 3},[br] {1, −1, −2}}[br][b]b[/b]=[math]\left(\begin{matrix}-6\\-5\\5\end{matrix}\right)[/math][br]det([b]A[/b]) =[size=150] (-2)[size=100]+(-6)+(-2) - (1+8+(-3)) = -10-6=-16[br]det([b]A1[/b])=[size=150] (12)+(-30)+(5) - ( 5+(-20)+18)= -13-3=-16[br][/size]det([b]A2[/b])=[size=150] (10+(-18)+10) -(-5+24+15)=2-34=-32[br][/size]det([b]A3[/b])=[size=150] (5+10+12)- (-6+(-20)+5)=27-(-21)=48[br][/size][/size][/size][math]\left(\begin{matrix}x\\y\\z\end{matrix}\right)[/math]=[math]\frac{1}{-16}\left(\begin{matrix}-16\\-32\\48\end{matrix}\right)[/math]=[b][math]\left(\begin{matrix}1\\2\\-3\end{matrix}\right)[/math] [/b]
[b]余因子やクラメールは2次・3次くらいなら扱いやすい。[br]しかし、次数が多くなると扱いにくいし、解が不定のときは使えない。[br]それに対して、行列と対話しながら、行列を変形を繰り返すことで、[br]段階的に、解や逆行列を求める方法がある。[br][color=#0000ff][size=150]行数、列数が多くても実行できるという利点、[br][/size]次数よりも情報が少ない場合でも、解を表現する手がかりが得られる利点がある。[br][/color][br]Ax[/b]=[b]b[/b]の[br]係数行列[b]A[/b]と値列ベクトル[b]b[/b]を横に並べて合体した[[b]A|b[/b]]を[b]拡大係数行列[/b]と呼ぼう。[br][size=200][color=#0000ff][A|b]が[E|u]に変形できたときに、x=uが解になる。[br]つまり、[br][/color][/size][color=#0000ff][size=150][size=200][A|E]が[E|A[sup]-1[/sup]]に変形できたときに, 逆行列A[sup]-1[/sup]が出る。[br][/size][/size][/color][br][b][size=150][size=200][color=#9900ff]・Aの行基本変形は3種ある。[br] i行とj行入れ替え (A[sub]i↔j[/sub])[br] i行の定数c倍 (A[sub]i[/sub](c))[br] j行の定数c倍を, j行にたす ( A[sub]i←j[/sub](c))[br][br][/color][/size][/size][/b]Aの部分の[color=#0000ff][b][size=150]対角成分を1に[/size][/b][/color]それ以外を[b][color=#0000ff][size=150]0に掃き出す[/size][/color][/b]イメージで[br][b][size=150][color=#0000ff]掃き出し法[Sweeping method]。[br][/color][/size][color=#0000ff](例)[br] x -2y +z =-6[br]2x+y +3z =-5[br] x -y -2z =5[br]掃き出し法で解くと?[br][/color][A|b][/b]=[br]{{1, -2, 1 , - 6},[br] {2,1, 3 , - 5},[br] {1, -1, -2, 5}}[br]1行目の-2倍を2行目にたす、1行目を3行目からひく。[br]{{1, -2, 1 , - 6},[br] {0,5, 1 , 7},[br] {0, 1, -3, 11}}[br]3行目の2倍を1行目にたし、2行目のー0.2倍を3行目をたす。[br]{{1, 0, -5 , 16}}[br] {0,5, 1 , 7},[br] {0, 0, -3.2, 9.6}}[br]3行目の1/(-3.2)倍。[br]{{1,0, -5 , 16},[br] {0,5, 1 , 7},[br] {0,0,1,-3}}[br]3行目の5倍を1行目にたし、3行目の−1倍を2行目からひく。[br]{{1,0, 0, 1 },[br] {0,5,0 , 10 },[br] {0,0,1,-3 }}[br]2行目を0.2倍する。[br]{{1,0, 0, 1 },[br] {0,1,0 , 2 },[br] {0,0,1,-3 }}[br][b]=[E|x][/b][br][b](x,y,z)=(1,2,3)[br][/b][color=#0000ff](例)[br][b]{{-3, 2, 2},[br] {-5,4, 3},[br] {1, -1, -1}}の逆行列は?[br][/b][/color][A|E]=[br]{{-3, 2, 2, 1, 0, 0},[br] {-5,4, 3, 0, 1, 0},[br] {1, -1, -1, 0, 0, 1}}[br]1行目の-5/3倍を2行目にたす、1行目の1/3倍を3行目にたす。[br]1行目を-3で割る。[br]{{1, -2/3, -2/3, -1/3, 0, 0},[br] {0,2/3, -1/3, -5/3, 1, 0},[br] {0, -1/3,-1/3, 1/3, 0, 1}}[br]2行目を1行目にたし、2行目の1/2倍を3行目をたす。[br]2行目を3/2倍する。[br]{{1, 0, -1, -2, 1, 0},[br]{0,1, -1/2, -5/2,3/2,0},[br]{0, 0,-1/2, -1/2, 1/2, 1}}[br]3行目の-2倍を1行目にたし、3行目の−1倍を2行目にたす。[br]3行目を−2倍する。[br]{{1, 0, 0, -1,0,-2},[br] {0,1, 0, -2, 1,-1},[br] {0, 0,1,1, -1, -2}}[br]=[E|A[sup]-1[/sup]][br][b]Aの逆行列は、[br]{{ -1, 0, -2},[br] { -2, 1, -1},[br] { 1,-1, -2}}[br][/b]
[size=150][size=100]連立方程式を解くためには、拡大係数行列のAの部分が単位行列になるまで[br]行の基本変形を続ける必要があったね。[br]解を求めるためではなく、解があるかどうか、あるとしたら1つか不定かを知りたいときがあるね。[br]その場合は、上の行から順に0でない係数が1列目、2列目、3列目と階段のようにならぶ[br]形になるように変形をしていけばよいね。[br][br][/size][/size][size=150][b]行の基本変形[/b]をして、0になる行を作っていくとき、0以外のある行数を[b][color=#0000ff]階数[rank][/color][/b]という。[br]言い換えると、もとの行数から、0だけの行ベクトルの行数を引いたものが階数だね。[br][/size][size=150]そのような変形で簡約化した行列のことを[color=#0000ff][b]階段行列[step matrix][/b][/color]と呼ぶことがもあるね。[br][/size][br]そのときの各行で0でない最初の列の数を[b]主成分[/b]といい、[b]主成分が1[/b]になるまで基本変形した場合、[br]その行列を[b]簡約な行列[/b]と呼ぶこともある。[br][br][size=150][size=200]n次の正方行列Aで、rankA=nならAは正則。[br]rankがn未満では正則でない。[br][color=#0000ff][b]n - rankAを自由度[degree of freedom][/b][/color]という。[br][/size][/size] [br][size=150][b]<同次連立方程式>[br][/b]拡大係数行列の最後の列が0ベクトルになる連立方程式を[b]同次とか斉次[/b]とか呼んだりする[br][/size][b][size=200]Ax=0 のこと。[br][/size][color=#0000ff][size=150]3D空間では、3つの行は、3つの原点通過平面に対応しているね。[br][/size][/color][br][/b]自由度=0だと解は原点[b]x=(x,y,z)=(0,0,0)[/b]の自明な、つまらないものだけになる。[br]自由度が正なら自明でない解がある。自由度1は直線、自由度2は平面に対応するね。[br][br][b][color=#0000ff](例)[/color]自由度0[br]A=[br][/b][color=#0000ff]{{1,−1,1,0},[br] {3,1,1,0},[br] {2,2,3,0}}[br][/color]⇒[br][b]{{1,−1,1,0},[br] {0,4,−2,0},[br] {0,0, 3,0}}[br][/b]⇒[br][b]n=3, rankA=3[br]自由度 n - rankA = 0[br]0パラメータ[/b][br][b]x=[/b](0,0,0)のみ。[br]3D空間では[b]解[/b]は[b]原点Oのみ[/b]。[br]だから、3行の3平面を表し、3平面が原点で交わっている関係だね。[br][br][b][color=#0000ff](例)[/color]自由度1[br]A=[br][/b][color=#0000ff]{{1,1,-1,0},[br] {4,3,-1,0},[br] {6,5,-3,0}}[br][/color]⇒[br][b]{{1,1,-1,0},[br] {0,-1,3,0},[br] {0,0,0,0}}[br][/b]⇒[br][b]n=3, rankA=2[br]自由度 n - rankA = 1[br]1パラメータp[br]x=[/b](-2p,3p,p)のみ。[br]3D空間では解は[b]原点通過の直線[/b]。[br]だから、3行の3平面を表し、3平面が原点を通る直線で交わっている関係だね。[br][br][b][color=#0000ff](例)[/color]自由度2[br]A=[br][/b][color=#0000ff]{{3,2,-1,0},[br] {6,4,-2,0},[br] {-12,-8,4,0}}[br][/color]⇒[br][b]{{3,2,-1,0},[br] {0,0,0,0},[br] {0,0,0,0}}[br][/b]⇒[br][b]n=3, rankA=1[br]自由度 n - rankA = 2[br]2パラメータp,q[br]x=[/b](p,q,3p+2q)のみ。[br]3D空間では解は[b]原点通過の平面[/b]。[br]だから、3行が見かけ上は3平面を表し、実際は1つの原点を通る平面になる関係だね。[br]自由度が2だと、自由すぎてつまらない。3平面が1つに重なっているからね。[br]
[size=150][b]<非同次連立方程式>[br][/b][size=100]拡大係数行列の最後の列が0ベクトルでない連立方程式を非同次とか非斉次とか呼んだりする。[br][/size][/size][b][size=200]Ax=b[/size][br][/b]3D空間では、各行は原点を通るとは限らない平面に対応している。[br][b]係数行列Aのランクr=rank(A)[/b]は平面の異なる[b]法線ベクトルの数[/b]に対応し、[br][b][b]拡大係数行列A~=[A|b][/b]のランクr~=rank(A~)[/b]は平面の法線と通過点でみた[b]事実上の平面数[/b]に対応する。[br][color=#0000ff][b][size=150]だから、係数行列のランクrより、拡大係数行列のランクr~が大きいと[br]法線ベクトルが同じで通過点が異なる平面、つまり、[u]平行面が混入[/u]する。[br]だから、共有空間はないので、解なし。[br][/size][/b][/color]拡大したものと[color=#0000ff][b]ランクが同じなら[u]平行面はないので、解がある[/u][/b][/color]。[br]さらに、[br]自由度=n - r=0なら、ただ1組の解をもつ。[br]自由度が正なら不定解。[br]自由度は同次連立方程式と同じで、解空間の次元(直線、平面、。。。)を表すね。[br][br][b][color=#0000ff](例)[/color]自由度0[br][/b]A=[br][color=#0000ff]{{1,2,-2,2},[br] {2,5,-1,8},[br] {-2,-1,11,18}}[br][/color]⇒[br][b]{{1,2,-2,2},[br] {0,1,3,4},[br] {0,0, 1,-5}}[br][/b]⇒[br][b]拡大前Aも後A~もランクは3。[br]n=3, rankA=3[br]自由度 n - rankA = 0[br]0パラメータ[/b][br][b]x=[/b](-46,19,-5)のみ。[br]3D空間では[b]点[/b]。[br][br][b][color=#0000ff](例)[/color]自由度1[br][/b]A=[br][color=#0000ff]{{1,3,-2,1},[br] {3,10,-3,4}}[br][/color]⇒[br][b]{{1,3,-2,1},[br] {0,1,3,1}}[br][/b]⇒[br][b]拡大前Aも後A~もランクは2。[br]n=3, r=rankA=2[br]自由度 n - r= 1[br]1パラメータp[br]x=[/b](11p-2,-3p+1,p)のみ。[br]3D空間では[b]原点を通らない直線。[br][br][b][color=#0000ff](例)ランクの不一致[/color][br]A=[br][/b][color=#0000ff]{{1,2,-3,2},[br] {3,6,-9,10}}[br][/color]⇒[br][b]{{1,2,-3,2},[br] {0,0,0,4}}[br][/b]⇒[br][b]係数のランクr=1で、拡大ランクr~=2[br]解なし。[br][/b]3D空間では、2つの平面の法線ベクトルが同じで交点なし。[br]平行な2平面は交わらない。[br][br][/b]
・解が1つに決まる連立方程式をできるだけいろいろな方法で解こう。[br][color=#0000ff] 3x +y +2z= 9[br] x -2y +z = 8[br] -2x - y +z= -3[br][/color][br]A={{3, 1, 2 }, b={{9},[br] {1, -2, 1}, {8},[br] {-2, -1,1}} {-3}}[br][b](掃き出し法)[br][/b]拡大係数行列は[A|b][br]A={{3, 1, 2, 9 },[br] {1, -2, 1, 8},[br] {-2,-1, 1, -3}}[br]行交換A[sub]2↔1[/sub][br]A={{1, -2, 1, 8 },[br] {3, 1, 2, 9 },[br] {-2,-1, 1, -3}}[br]A[sub]2←1(-3)[/sub], A[sub]3←1(2)[/sub][br]A={{1, -2, 1, 8 },[br] {0, 7, -1, -15 },[br] {0,-5, 3, 13}}[br]A[sub]2←3[/sub],[br]A[sub]2(1/2)[br][/sub]A={{1, -2, 1, 8 },[br] {0, 1, 1, -1 },[br] {0,-5, 3, 13}}[br]A[sub]3[/sub]←2(5), [br]A[sub]3(1/8)[br][/sub]A={{1, -2, 1, 8 },[br] {0, 1, 1, -1 },[br] {0, 0 1, 1}}[br]A[sub]1←2(2),[br]A1←3(-3)[br][/sub]A={{1, 0, 0, 3 },[br] {0, 1, 1, -1 },[br] {0, 0 1, 1}}[br]A2←3(-1)[br]A={{1, 0, 0, 3 },[br] {0, 1, 0, -2 },[br] {0, 0 1, 1}}[br](x,y,z)=(3,-2,1)[br]掃き出し法のいいところは、公式を覚えなくてもよいこと。ただし、基本変形を大量にするため、[br]そのことによる計算ミスをしやすい。[br][b](余因子)[br]A={{3, 1, 2 }, [br] {1, -2, 1}, [br] {-2, -1,1}} [br] [/b]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[br]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,[br]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,[br]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[br]A-1=1/-16tr{{-1, -3, -5},[br] {-3, 7, 1},[br] {5, -1, -7}}[br]-1/16{{-1, -3, 5},[br] {-3, 7, -1},[br] {-5, 1, -7}}[br](x,y,z)=[br]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)}[br]=-1/16{ -9-24-15, -27+56+3, -45+8+21}[br]=-1/16{-48, 32, -16}[br] =(3, -2, 1) [br]余因子法のいいところは、行列式の個数が増えても、次数が下がること。[br]ー符号を1つおきに交互につけるのを忘れたり、転置するのを忘れるとアウト。[br][b](クラメール)[br][/b]det(A) =[u]3[/u]・(-2)・1+ [u]1[/u]・(-1)・2+[u](-2)[/u]・1・1 -{[u](-2)[/u]・(-2)・2+[u]1[/u]・1・1+[u]3[/u]・(-1)・1}=-6-2-2-{8+1-3}=-10-6=-16[br]det(A1)=[u]9[/u]・(-2)・1+ [u]8[/u]・(-1)・2+[u](-3)[/u]・1・1 -{[u](-3)[/u]・(-2)・2+[u]8[/u]・1・1+[u]9[/u]・(-1)・1}=-18-16-3-{12+8-9}=-37-11=-48[br]det(A) =3・[u](-2)[/u]・1+ 1・[u](-1)[/u]・2+(-2)・[u]1[/u]・1 -{(-2)・[u](-2)[/u]・2+1・[u]1[/u]・1+3・[u](-1)[/u]・1}[br]det(A2) =3・[u]8[/u]・1+ 1・[u](-3)[/u]・2+(-2)・[u]9[/u]・1 -{(-2)・[u]8[/u]・2+1・[u]9[/u]・1+3・[u](-3)[/u]・1}=24-6-18-{-32+9-9}=0+32=32[br]det(A) =3・(-2)・[u]1[/u]+ 1・(-1)・[u]2[/u]+(-2)・1・[u]1[/u] -{(-2)・(-2)・[u]2[/u]+1・1・[u]1[/u]+3・(-1)・[u]1[/u]}=[br]det(A3) =3・(-2)・([u]-3)[/u]+ 1・(-1)・[u]9[/u]+(-2)・1・[u]8[/u] -{(-2)・(-2)・[u]9[/u]+1・1・([u]-3)[/u]+3・(-1)・[u]8[/u]}=18-9-16-{36-3-24}=-7-9=-16[br](x,y,z)=1/-16 { -48, 32, -16}= (3, -2, 1)[br]クラメールのいいところは、行列式の個数が未知数の数だけしかないこと。ただし、bのAへの代入で[br]行列そのものを変えてしまったもので行列式を計算するためもとの行列と変わったところの確認が大切。