1~2万円代のコスパ重視、人気メーカーのサウンドバー、シアターバー【Bluetooth対応】

DeepFaceLabのインストールと実装【ディープフェイク】

【DearMob iPhoneマネージャー】PCからiPhoneへ、iPhoneからPCへ写真転送と管理に最適なソフトのご紹介

PCのバックアップに最適なソフト「EaseUS Todo Backup Home 11.5」を使ってみた感想

【Pythonista3】複数のSceneを切り替える方法

スマホ画面をフッ素コーティングして指紋防止をしてみた

【DearMob iPhoneマネージャー】iPhoneのバックアップやデータ移行に最適なソフトのご紹介

【DearMob iPhoneマネージャー】PCからiPhoneへ、iPhoneからPCへ動画転送と管理に最適なソフトのご紹介

エルゴトロン VESA対応モニターアーム　AmazonベーシックやHPのOEM商品との違い

iCUE制御のゲーミングPCを自作しました

【BlueLounge】MacBookなどの薄型ノートパソコンのおすすめフリップスタンド

【Oculus Quest2】本体、付属品の一覧

【Logicool】G913TKL White おすすめポイントをレビュー！

ポケモンなどの位置情報ゲームに使える位置偽装ツール「iAnyGo」

Pythonista3でアクションゲームをつくってみる【iOSアプリ開発】

【スマホでプログラミング】Pythonでアプリ開発！【ダイアログ表示・センサ取得#2】

Matrix Power Watchが故障したので交換してもらった話

【最新版】人気なワイヤレスゲーミングマウスを比較/Logicool,Corsair,Razer,Steelseries

【Python】ヤコビ法とガウス・ザイデル法で連立一次方程式を解く

Python

2020.04.28 2019.05.30

コンピュータは有限桁の数値しか扱う事はできないので、桁数の多い場合や無限小数の場合は四捨五入され切り捨てられます。なので実際の数値とは多少の誤差が生じますが、これを丸め誤差といいます。

なので、コンピュータを用いて連立方程式等の計算すると多少の誤差が生じますが、誤差が最小になるように高精度で計算をすれば限りなく正しい解を求めることができます。

今回は、連立一次方程式を解くための「ヤコビ法」と「ガウス・ザイデル法」についてプログラムと共に紹介します。

Contents

ヤコビ法 -Jacobi Method-
ガウス・ザイデル法 -Gauss Seidel-
ヤコビ法とガウス・ザイデル法の特徴

ヤコビ法 -Jacobi Method-

ヤコビ法とは、n元の連立一次方程式の連立一次方程式 $A{\vec {x}}={\vec {b}}$

元の方程式である[Ax = b]を[r = b – Ax]と変形し、rの値が小さくなるようにxの値を反復的に変えていき、r = 0を満たすxを見つけることができれば、方程式が解けたことになります。

例えば以下のような3元の連立一次方程式を解くことを考える。

$\left({\begin{array}{ccc}a_{{11}}&a_{{12}}&a_{{13}}\\a_{{21}}&a_{{22}}&a_{{23}}\\a_{{31}}&a_{{32}}&a_{{33}}\end{array}}\right)\left({\begin{array}{c}x_{1}\\x_{2}\\x_{3}\end{array}}\right)=\left({\begin{array}{c}b_{1}\\b_{2}\\b_{3}\end{array}}\right)$

$k$ 回目の反復で得られた $x_1$ の値を $x_{1}^{{(k)}}$ と書く。初期値 ${\vec {x}}^{{(0)}}$ 。

$x_{1}^{{(k+1)}}=(b_{1}-a_{{12}}x_{2}^{{(k)}}-a_{{13}}x_{3}^{{(k)}})/a_{{11}}$

$x_{2}^{(k+1)}=(b_{2}-a_{21}x_{1}^{(k)}-a_{23}x_{3}^{(k)})/a_{22}$

$x_{3}^{(k+1)}=(b_{3}-a_{31}x_{1}^{(k)}-a_{32}x_{2}^{(k)})/a_{33}$

という反復を繰り返していく。また、収束判定として、適当な正数εに対して、

　

また、

が成り立てば終了する。収束判定の値であるεを厳しくすれば精度が高まります。

次のような簡単な3元連立方程式を解いてみる。

3a + 2b + c = 10

a + 4b + c = 12

2a + 2b + 5c = 21

以下にPythonのプログラムです。

実行結果：

ガウス・ザイデル法 -Gauss Seidel-

ガウス・ザイデル法も同様にn元の連立一次方程式を解く反復法ので解く手法の１つである。

以下のように、k回反復した後に、k+1回目の解を求める際にすでに求めた解を次の方程式に利用する。

$x_{{i}}^{{(k+1)}}={\frac {1}{a_{{ii}}}}\left(b_{{i}}-\sum _{{j=1}}^{{i-1}}a_{{ij}}x_{{j}}^{{(k+1)}}-\sum _{{j=i+1}}^{{n}}a_{{ij}}x_{{j}}^{{(k)}}\right)$

このようにすることで、ガウス・ザイデル法はヤコビ法よりも収束の速度が速くなります。またガウス・ザイデル法では反復計算で書き換えた変数を使うだけであるが、ヤコビ法では「1つ前の値」を記憶する必要があるので、プログラム上でもガウス・ザイデル法の方が簡単であることが言える。

同様の3元連立一次方程式を解きます。以下がPythonのプログラムです。

実行結果：

ヤコビ法より、計算回数が非常に速いことがわかります。

ヤコビ法とガウス・ザイデル法の特徴

▶ガウス・ザイデル法は、ヤコビ法よりも収束する速度が速く、計算量が少ない。

▶プログラムもガウス・ザイデル法の方が簡潔である。（反復計算で書き換えた変数をそのまま用いるため）

しかし、ガウス・ザイデル法は書き換えた変数を用いるため逐次計算となる。そのため並列計算には向かない。また、ヤコビ法は各式で独立に計算ができるため、並列計算に向いている。

以上から、一概にガウス・ザイデル法が良いとは言えません。また、ヤコビ法やガウス・ザイデル法を緩和係数を用いて組み合わせた「SOR法」もありますが、今回は割愛します。