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

A 回答 (1件)

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



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

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

どうしても駄目だったら分解法なり反復法なりを試してみるしかないですね。
    • good
    • 0
この回答へのお礼

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

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

お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!

今、見られている記事はコレ!

おしトピ編集部からのゆる~い質問を出題中

お題をもっとみる


このカテゴリの人気Q&Aランキング

おすすめ情報

カテゴリ