もう一度(これでだめな時は更新してください。) : 新しいNと場所で描き直します。
非線形系の反復法たちいろいろな関数ニュートン法の図弓弦法の図正割法の図はさみうち法の図
極値を求める
ニュートン法の図
極値を求める
弓弦法の図
極値を求める
正割法の図
極値を求める
はさみうち法の図
連立方程式関係ガウスの消去法ガウス・ジョルダンの掃き出し法ヤコビ法ガウス・ザイデル法連立方程式の解法まとめて元気よく
LU分解コレスキー法---
フーリエ解析関係周期関数たち周期関数に対する
フーリエ級数展開
非周期関数に対する
フーリエ級数展開
ストークスの波動公式-
ガウス・ジョルダンの掃き出し法(2011年6月5日公開、2011”N06ŒŽ12“ú01:28:27第1回の改訂)
●注意

上はJAVAで作られています。メモリを大量に使ったり、重くなるかもしれません。その時は、ごめんなさい。
実行後に画面をスクロールしたり、アプレット全体が画面に入ってないと、間違った画面になるかもしれないので、気をつけてください。画面の大きさを決めてから”もう一度”をクリックするか、更新(reload)してください。


ガウス・ジョルダンの掃き出し法は、ホームページで数式を書くのは、ちょー面倒なので詳しくは参考文献を読む方がよろしいかとは思いますが、一次の連立方程式Ax=b(Aは行列、bはベクトル、xは答えのベクトル)でxを求める方法の一つです。
ここでは、得られた答えxと、Axの結果、Axとbの差を表示して楽しんでいます。わーい。
なお、Aとbの要素を*.**の形式で表示していますが、実際には桁数をたくさん使っています。表示の都合で*.**にしています。

ここでは、求めたい連立方程式を
a11x1 + a12x2 +・・・+ a1jxj +・・・+ a1nxn = b1
a21x1 + a22x2 +・・・+ a2jxj +・・・+ a2nxn = b2
・・・
ai1x1 + ai2x2 +・・・+ aijxj +・・・+ ainxn = bi
・・・
an1x1 + an2x2 +・・・+ anjxj +・・・+ annxn = bn
と表現します。これを行列とベクトルを使ってあらわすと、ただし、行列とベクトルで使う括弧()はかっこ悪いけど使っていませんが・・・
A =
a11 a12 ・・・ a1j ・・・ a1n
a21 a22 ・・・ a2j ・・・ a2n
・・・
ai1 ai2 ・・・ aij ・・・ ain
・・・
an1 an2 ・・・ anj ・・・ ann
x=
x1
x2
・・・
xi
・・・
xn
b=
b1
b2
・・・
bi
・・・
bn
を使って、Ax=bとなります。

ガウス・ジョルダンの掃き出し法は、
-------------
STEP 1
aij(1) = aij
bi(1) = bi
STEP 2
k = 2,3, ..., n + 1に対して、
aij(k) = aij(k-1) (jk-1のとき)
aij(k) = aij(k-1) (j>=k-1でi = k - 1のとき)
aij(k) = aij(k-1) - ai(k-1)(k-1) a(k-1)j(k-1) / a(k-1)(k-1)(k-1) (jk-1でik-1のとき)
aij(k) = 0 (j = k-1でik-1のとき)
bi(k) = bk-1(k-1) (i = k-1のとき)
bi(k) = bi(k-1) - ai(k-1)(k-1)b(k-1)(k-1) / a(k-1)(k-1)(k-1) (ik-1のとき)
を計算する。
STEP 3
i = 1, 2, ..., nの順に、
xi = bi(n+1) / aii(n)
を計算する。
-------------
というプロセスでxを求める方法です。

●参考文献
有本卓著、コロナ社、数値解析(1)
川上一郎著、岩波書店、数値計算
●プログラムのダウンロード
○Java(gj.java)


ご意見、ご感想、お問い合わせ、お願い等がございましたら、お気軽に、
メール送信フォームからメールを送るか、
●掲示板に書き込むか、
どちらかお好きな方法で、ご連絡お願いいたします。


●大山崇のホームページの利用について
●大山崇のホームページ