- 締切済み
fortran gauss消去法のhilbelt行列について
今fortran90をやっています。で、ついこないだhilbelt行列をgauss消去法で解いたのですが、どうやっても４元までのhilbelt行列じゃないと答えが狂ってきます。どうしてでしょうか
- yusuke641
- お礼率37% (29/77)
- その他（プログラミング・開発）
- 回答数1
- ありがとう数2
- みんなの回答 （1）
- 専門家の回答
みんなの回答
- bob
- ベストアンサー率50% (52/103)
これだけの情報量ではなんとも言いがたいのですが、とりあえず５次以上の対角行列とか簡単な問題では正解が得られているんですよね？ Hilbert行列は定義から分かるように次数が上がる程対角成分が小さくなるので直接法でそのまま解くと誤差がたまりやすくなります。しかし、５次で解けなくなる程かと言われると、やってみないと分からないといった感じです... その辺りを考慮して演算の順番を誤差がたまりにくいように考慮するとか、なんらかの前処理を加えてみるなどの工夫をしてみてはどうでしょう？ どうしても駄目だったら分解法なり反復法なりを試してみるしかないですね。
関連するQ&A
- 待ち行列のプログラム（Fortran)
Fortranで待ち行列（レジ待ち）のプログラムを作りたいのですが、途中で行き詰まってしまい、うまくできません！ 待ち行列についてまだよく理解できてないのだと思います。指数乱数を発生させるプログラムだけだったらできるのですが、それをヒストグラムで示して、待ち行列がどのように変化するのか観察するプログラムができないのです。 待ち行列のプログラムについて詳しく載っているHPなどがあれば教えてください！！（本でもOKです）
- ベストアンサー
- その他（プログラミング・開発）
- fortran 掃き出し法で逆行列
大学の授業でfortranを使って掃き出し法による３×３の逆行列を求める問題がでましたが求め方がわかりません。 どなたか教えてください！ ちなみに行列は 1 2 1 2 1 -2 3 5 1 です。
- 締切済み
- その他（プログラミング・開発）
- C言語とfortranでのデータの受け渡しについて
はじめましてc_f_gaussと申します。 現在fortranで動的に配列を得るためにC言語を使用し 配列を確保しそのアドレスをfortranで定義してある 配列のアドレスと交換したいのですがどのようにしたら よいでしょうか？ 今の段階ではC言語側でアドレスを交換してもfortran 側に戻ると元のアドレスに戻ってしまいます。 fortranとC言語でデータ受け渡すには何か必要な手続き があるのでしょうか？
- ベストアンサー
- C・C++・C#
- fortran
fortranでn次正方行列の対角要素の和を求める関数副プログラムを作ろうと思います。 fortran初心者で、何から始めていいのかさっぱりワカリマセン。 ヒント・アドバイスがあればよろしくお願いいたします。
- 締切済み
- その他（プログラミング・開発）
- 逆行列求めたい
今逆行列計算を行いたいと思い、LAPACK(Fortran)でがんばっておりますが、簡単な行列だとうまくいきますが、状態の悪い行列ではうまくいかないみたいです。計算上このような行列を計算しなければならなく困ってます。 ZGETRI,ZGETRFをつかってますが,いろいろと調べるとZLANGE,ZGECONを使っているとも見つけましたが、使い方がいまいちわかりません。 どなたかご教授いただけないでしょうか？
- 締切済み
- その他（プログラミング・開発）
- FORTRANについて
今、大学の研究室でFORTRANの勉強をしています。そこで、院生から問題を出されました。それは、 テキストファイルから数字を読み込み、それらをFORTRANの中で計算し、その答えをまたテキストファイルに出力する、というプログラムです。 計算するところまではできたのですが、どうしてもテキストファイルに出力するプログラムがわかりません。 出力の仕方だけでも教えてください。教科書にもなく、図書館でも古い時代のものばかりで。しかも、僕は学校の関係でFORTRANを全く知りません。
- 締切済み
- その他（プログラミング・開発）
- FORTRAN
初めまして。 今大学でFORTRANについて議論が交わされてるのですが FORTRANについて知ってるいる方がいたら是非教えてください。 今のところ求めている情報としては １．FORTRANの強み ２．FORTRANの弱み ３．FORTRANのこれから ４．FORTRANが抱えている問題 URLなどでも結構ですので是非情報をお持ちでしたら教えてください。
- ベストアンサー
- その他（プログラミング・開発）
- 行列の消去法のコツなど教えてください。
只今、学校にて行列を習っているわけですが、最近行列を使った消去法を習い始めました。 たとえば ３ １ －７ ０ ４ －１ －１ ５ １ －１ ２ ２ このような行列があったとします。 習った方法は、 (1)一つの行に０でない数をかける。 (2)一つの行にある数をかけたものを他の行に加える。 (3)二つの行を交換する。 １ ０ ０ ３ ０ １ ０ ５ ０ ０ １ ２ このような式に変形してx=3,y=5,z=2みたいな感じにするということでしたが、 今回教えていただきたいことは、 →1度に前述の3つの式を何回も使っていいのか。 →うまく変形するコツ。 の二つです。 やり方自体はなんとなくわかるのですが、単位行列に持っていくまでの手順がイマイチ難しくわからないので、よろしければご教授願います。 2月頭辺りからテストなのでズバリを突いて欲しいと思います。 よろしくお願いします。
- ベストアンサー
- 数学・算数
- FORTRANについての質問です。
最近課題でFORTRANをやっているのです。 サブルーチンのMINVS/D/Q/C/Bというのを使って３×３の行列の逆行列を出力するプログラミングを作っているのですが、うまくできません。 サブルーチンを使うのはほとんどはじめてなのでどう使えばいいのかさえ良くわかっていません。 なにかアドバイスがありましたらお願いします。
- ベストアンサー
- その他（プログラミング・開発）
お礼
ありがとうございます。かなり理解しました。前処理を施して、もう一度プログラミングしてみます。