fortran gauss消去法のhilbelt行列について

今fortran90をやっています。で、ついこないだhilbelt行列をgauss消去法で解いたのですが、どうやっても４元までのhilbelt行列じゃないと答えが狂ってきます。どうしてでしょうか

No.1 回答者： bob 回答日時： 2001/07/09 02:06

これだけの情報量ではなんとも言いがたいのですが、とりあえず５次以上の対角行列とか簡単な問題では正解が得られているんですよね？

Hilbert行列は定義から分かるように次数が上がる程対角成分が小さくなるので直接法でそのまま解くと誤差がたまりやすくなります。しかし、５次で解けなくなる程かと言われると、やってみないと分からないといった感じです...

その辺りを考慮して演算の順番を誤差がたまりにくいように考慮するとか、なんらかの前処理を加えてみるなどの工夫をしてみてはどうでしょう？

どうしても駄目だったら分解法なり反復法なりを試してみるしかないですね。

この回答へのお礼

ありがとうございます。かなり理解しました。前処理を施して、もう一度プログラミングしてみます。

お礼日時：2001/07/10 02:07