gauss-jordanの消去法.pptVIP

* 連立1次方程式の数値解 例 nｘn 連立1次方程式 行列とベクトル記号で表す連立1次方程式 A：係数行列； X：未知数ベクトル； ｂ：定数ベクトル 連立1次方程式 A：係数行列； X：未知数ベクトル； ｂ：定数ベクトル 連立1次方程式 A-1：係数行列Aの逆行列 連立1次方程式 １．直接法 2．反復法 　Gauss-Jordanの消去法 　LU分解法 　ヤコビ（Jacobi）法 　ガウス―ザイデル（Guass-Seidel 　法 　SOR法（逐次緩和法） 　収束条件?終了判定 Gauss-Jordanの消去法 Gauss-Jordanの消去法 連立方程式の性質： 方程式の両辺に０でない数を掛けたり、割ったりしても ある方程式を他の方程式に加えたり、引いたりしても その連立方程式の解が変わらない。 消去法はこれらの性質を利用して、連立方程式の係数を次々に消去して、解の求めやすい形にする。 Gauss-Jordanの消去法の手順 １．Aの拡大行列を作る ２．行列を変形する ① 正規化 ② 方程式の解 ３．単位行列を作って、解を求める ｎ 実用的な消去法 上三角行列にする 方程式の解 応用例：逆行列の計算 消去法 消去法の注意点 演算の途中でakkが０になったとき、適当な行または列同士の入れ換えを行って、０でない要素をピボットに選んで、掃出すことが必要である。 akk ：ピボット（対　　角線上の要素） ガウス消去法の利点 必ず収束で、繰り返し回数≦方程式の元数ｎ。 手順が簡単、計算量が少ない。 多少性質が悪い係数行列にも適用可能。 ガウス消去法の欠点 正規化によって丸め誤差が生じる。 消去の途中の減算によって桁落ちが生じる。 乗算回数が多くなって誤差の累積する可能性がある。 元の係数行列が破壊してしまう。 LU分解法 下三角行列 上三角行列 AのLU分解 方程式の数：ｎｘｎ 未知数の数：ｎｘｎ LUの保存 連立1次方程式?反復法 ヤコビ法（Jacobi Jacobiの反復計算式