締切済み

fortran

2007/06/15 22:50

fortranでn次正方行列の対角要素の和を求める関数副プログラムを作ろうと思います。 fortran初心者で、何から始めていいのかさっぱりワカリマセン。ヒント・アドバイスがあればよろしくお願いいたします。

mari0725
お礼率59% (62/104)

その他（プログラミング・開発）
回答数1
ありがとう数1

みんなの回答 （1）
専門家の回答

みんなの回答

Tacosan
ベストアンサー率23% (3656/15482)

2007/06/16 01:58 回答No.1

・行列の表し方・(一般論として) 和の求め方・ループの書き方くらいはおさえときたいねぇ. まあ, 「Fortran の初心者」というだけなら「別の言語で書いてみれば」と言って終わりなんだけど.

関連するQ&A

行列　対角和　トレース

対角和（トレース）について質問させて下さい。対角和（トレース）は、n次正方行列の対角成分の総和を表しますが、この対角和とは一体なにを表すのでしょうか？なんのために対角和を求めるのか素朴な疑問ですが教えて頂けないでしょうか？以上、よろしくお願い致します。
　Fortran90を授業でやっています。

　Fortran90を授業でやっています。入力されたNまでの素数を配列で返す関数をつくって、それを主プログラムで呼び出して素数を出力するという問題です。　素数を関数を用いて配列で返すまでは出来ました。しかし、その関数（　Sosuu(N)としました。　）の要素数を(1000)と設定したため、Sosuu(N)の1番目から100番目までは素数が入っているのですが、それ以降の要素は値がすべて0になってしまいます。　なので、WRITE(*,*)Sosuu(N) とすると、素数が列挙された後さらに0が列挙されてしまいます。 0が出力されない方法は何かありませんか？
行列の固有値とトレースについて

問　正方行列のトレース（対角成分の和）は、その固有値の総和になる。この問題は行列が対角化可能ならば成り立ちますが、対角化不可能の場合でも成り立つのでしょうか？ご指導よろしくお願いします。
FortranにあってCにない関数？について

10年間Fortran90でプログラミングをしてきましたが、 C言語を使う必要があって、プログラムを書き直し中の者です。入門書を数冊読んでいて疑問に思ったのですが、Fortranには、 A=SUM(B) と書けば、配列Bの全ての総和をAに返す、という関数があるのですが、C言語に同じような関数はないのでしょうか？（標準ライブラリに無いという事はないのでしょうか？） Fortranと比べてCは組み込み関数が少ないという印象を持ったのですが実際そうなのでしょうか？また、Fortanでは配列の和をとる時、 do i=0,n a(i)=b(i)+c(i) end do を一行で、 a(0:n)=b(0:n)+c(0:n) と書いても良いですし、もし定義されている全部の配列の要素が対象なら a=b+c と一行で書いても良いのですが、ｃ言語では似たような書き方はできないのでしょうか？つまり、多次元配列を扱う時、Fortranでは単に「a=b+c」と書けば済むところを、Cでは for (i=0,i<=n,i++){ for (j=0,j<=n,ｊ++){ for (k=0,k<=n,ｋ++){ a[i][j][k]=b[i][j][k]+c[i][j][k]; }}} と書かなくてはいけないのか？と悩んでおります。私がCの機能を知らないだけなのか、Cとはそうゆうことが出来ないのか、浅学なため判断できません。ご教授ください、どうぞよろしくお願いいたします。
行列の対角化について

n次正方行列Aがある対角行列と相似の時、行列Aの各特性根αに対する固有空間の次元の和がnになることは分かるのですが、各特性根αの（特性方程式の）重複度と固有空間の次元が一致するのがなぜだか分かりません。どなたか教えてください。お願いします。
Fortran90 初心者です

Fortran初心者です。「3-1=」を文字列データとして入力された式の答えを数値で出力するプログラムをつくりたいのですがよくわかりません。どなたかヒントください。御願いします。
fortran90/95のSUM関数について

fortran90/95のSUM関数について fortran90/95のSUM関数について、SUM関数を使うときは範囲を指定できるのでしょうか？例えば、 integer A(3,2)/1,2,3,4,5,6/ A=1,4 2,5 3,6 の配列があるとします。列の和を求めたいのですが、最後の行（3行目）を足したくありません。結果としては 1+2=3 4+5=9 にしたいのです。このようなときは、やはりdo文にするしかないのでしょうか。実際のプログラムではdo文があまりにも多く、見づらいのでSUM関数を使いたいのですが... どなたかご存知の方、お教えいただけないでしょうか。よろしくお願いします。
C言語プログラム

N次の複素正方行列S,Tの積Uを計算するプログラムを作りたいのですが、実数で正方行列を計算するプログラムと複素数の積、和のプログラム struct complex { double re; double im; }; を作ったのですが、この二つをまとめるとプログラムができるらしいのですがまとめ方が全然わかりません。どのようにしたらいいかヒントなど教えてください。
待ち行列のプログラム（Fortran)

Fortranで待ち行列（レジ待ち）のプログラムを作りたいのですが、途中で行き詰まってしまい、うまくできません！　待ち行列についてまだよく理解できてないのだと思います。指数乱数を発生させるプログラムだけだったらできるのですが、それをヒストグラムで示して、待ち行列がどのように変化するのか観察するプログラムができないのです。　待ち行列のプログラムについて詳しく載っているHPなどがあれば教えてください！！（本でもOKです）
Fortranについて

私は今大学生のものです。今卒業研究でFortranのプログラムを用いています。しかし私はプログラミングには疎く、独学で勉強しようと思うのですがなにか初心者にわかりやすいお勧めの本、サイトなどはあるでしょうか？本屋に行ってもFortranの本は販売されておらず、通販では手にはいるのですがどの本がいいのかわかりません。是非教えていただければと思いここに書き込ませていただきました。また、家でもFortranのプログラムをコンパイル、実行できる環境をつくりたいのですがどうすればよろしいでしょうか？学校ではcygwinをインストールして行っていますが家でもcygwinをインストールすればいいのでしょうか？また今もととなるプログラムはあるのですがこのプログラムはLinux上では正常にコンパイルされないところがあるらしく書き直さなければならないらしいのですがLinux上でのFortranのプログラムは普通のFortranのプログラムとは違うのでしょうか？ちなみに学校のパソコンのOSはWindows Vistaです。初心者の質問ばかりになってしまってすみません。自分でしっかり学びたいと思うのでぜひ教えてください。
fortranのプログラム

fortranのプログラム現在、fortranの勉強をしております。そこで、質問があるのですが、ある関数ｆ(x,y,z)の座標（x,y,z）の値がデータとして与えられているとき、 S=10+f をfortranで計算したいと考えております。ただ、関数ｆは複数（ｆ１、ｆ２、ｆ３）あり、次々とｆに代入してＳを計算したいのですが、どのようにプログラムしたらいいか思いつきません。どなたか、ヒントだけでもいいので、教えてください。ちなみに、私が考えたプログラムは(下のプログラムはポイントだけ書いてあります。endやその他関係ないと思われるところは省いております。) do 100 k=1,3 S=S+fk(x,y,z) continue function f1 f1(s,t,u)=・・・ return end f2(s,t,u)=・・・ return end 使用しているバージョンは、fortran77(本当は90を使っているのですが、77だけで書いています)です。
サラスの公式

サラスの公式（実3次正方行列の行列式の公式）の証明が分からず困っています。「まず、27個の行列式の和に分解し、そのうち0ベクトルの行を持つものを消去し、各行列式を性質（行の入れかえ）によって、対角行列の行列式にし、成分の積に直す。」という方法の証明が知りたいのですが…何かアドバイスやその証明が載っているURLなど、よろしければお願いします！
行列の積の可換条件

線形代数で二つのｎ次正方行列が可換になる条件とはどんなものなのか？特に対角化できない行列Aに対し、交換可能な行列Bはどんなものか？それについて詳しく書いてある本を教えてください。お願いします。
CとFORTRANの計算速度はどちらが速いですか？

マシンスペックは同じとして、OSはUNIXであった場合、科学技術計算（微分方程式を解いたり、行列を対角化したり。。。）の計算速度はどちらが速いですか？普通、科学計算というとFORTRANを使っていますが、これはライブラリの遺産が多いからだと思いますが、実際のところFORTRANとCではどちらが計算速度が速いのでしょうか。
対称行列　対角行列

対角行列と対角化について質問させて頂きます。対角行列は、対角成分以外が0の正方行列です。対称行列は、t^A=Aが成り立つ正方行列Aです。ここで、対称行列の定理で、・対称行列の異なる固有値に属する固有ベクトルは直交する。というものがあるのですが、これは対角行列にも言えるのでしょうか？対角行列は対称行列なので言えると思いますが、テキストに特に記載がなかったので質問させて頂きました。以上、ご回答よろしくお願い致します。
線形・行列の証明がさっぱり。。。

対角成分以外がゼロである正方行列を対角行列という。対角行列の固有値は、対角成分に等しいことを示せ。また、対角成分より左下(右上)の成分がゼロである正方行列を上三角行列(下三角行列)という。上三角行列、下三角行列の固有値が対角成分に等しいことを証明せよ。この証明がさっぱりわかりません。ご指導お願い致します。
ユニタリ行列って？？

ユニタリ行列ってなんですか？ユニタリを満たすと、量子力学や、線形代数学において、どのような意味をもつのでしょうか？行列の要素に簡単な数値を用いて説明してもらえるとうれしいです。次の文章は、自分で調べてみたけど、いまいち意味がわからなかったことです。複素正方行列をUとすると、そのエルミート共役がその逆数に等しいとき、ユニタリと呼ばれるんですか？ U^†＝U^(-1) （１）エルミート行列Aの対角要素は相似（ユニタリ）変換により、要約される。 D=U^†AU ここでUは列が行列Aの直交ベクトルであるユニタリ行列で、実数対角行列で、対角要素は行列Aの固有値であるとありました。
この行列問題を早くとく方法

実３次正方行列ＡＡ＝|(7/2) 0　　-(9/2)　| 　　　| 1　　　1　　　1　　　| 　　　|-(9/2) 0　(7/2)　　| についてＢ^3=A となる実３次正方行列Ｂを一つ求めよという問題ですが、以下の方法より早く解く方法がありましたら教えてください (*は行列同士の積の操作とします) 私はＢは対角化可能の行列とし、そこからＢ^3 =P *(P^(-1)*B*P)^(3)*P^(-1) =A これに左からP^(-1)、右からPをかけ (P^(-1)*B*P)^(3)= P^(-1)*A*P そしてＡの固有値を求め、固有ベクトルからＰとＰ^(-1)を求め、Ａを対角化する（対角要素は固有値）対角化されたＡの要素の三乗根を求め左からＰ、右からP^(-1)をかけＢを求めましたしかし、これよりもっと早く解ける方法があると思うのです。ポイントはＢを一つ挙げよ、というとこだと思うのですが・・どうかおねがいします
線形代数　行列　対角化

対角化について質問させて頂きます。対角化とは、「正方行列を適当な線形変換により、もとの行列と同値な対角行列に帰着させること。」と説明がありました。ここで、同値とは具体的にどのような内容を指すのでしょうか？また、対角化を求める際、正方行列Aに対してP^-1APとなる正則行列Pを求めます。この正則行列Pは正方行列Aより求めた固有値に属する固有ベクトルを並べたものになりますが、これはなぜですか？なぜ、固有ベクトルを並べたものが正則行列Pになるのでしょうか？以上、ご回答よろしくお願い致します。
fortran90

fortran90 現在 fortran90 で計算プログラムを作っております。プログラムはかけたのですが、実行画面にエラーが出ます。解決方法が分かりません。どなたか分かる方いらっしゃいますか。プログラム implicit real*8 (a-h,o-z) parameter(n=900) open (1,file='004.txt',status='old') do 10 i=1,n read (1,*,end=11)xx,yy,zz,ww theta=acos(zz) phi1=acos(xx/sqrt(1-zz*zz)) phi2=asin(yy/sqrt(1-zz*zz)) write(*,*)theta,phi1,phi2 10 continue 11 end プログラム内容あるフォルダーに保存してあるテキストファイルから４×４行列のデータを読み込み、そのデータからθの値とφ１、φ２のアタを求めるプログラムです。実行分のエラー画面には、 acos: DOMAIN error です。

fortran

みんなの回答

関連するQ&A

行列　対角和　トレース

Fortran90を授業でやっています。

行列の固有値とトレースについて

FortranにあってCにない関数？について

行列の対角化について

Fortran90 初心者です

fortran90/95のSUM関数について

C言語プログラム

待ち行列のプログラム（Fortran)

Fortranについて

fortranのプログラム

サラスの公式

行列の積の可換条件

CとFORTRANの計算速度はどちらが速いですか？

対称行列　対角行列

線形・行列の証明がさっぱり。。。

ユニタリ行列って？？

この行列問題を早くとく方法

線形代数　行列　対角化

fortran90

注目のQ&A

カテゴリ
一覧

専門家に質問してみよう
専門家登録

あなたにピッタリな商品が見つかる！ OKWAVE セレクト

fortran

みんなの回答

関連するQ&A

行列 対角和 トレース

Fortran90を授業でやっています。

行列の固有値とトレースについて

FortranにあってCにない関数？について

行列の対角化について

Fortran90 初心者です

fortran90/95のSUM関数について

C言語プログラム

待ち行列のプログラム（Fortran)

Fortranについて

fortranのプログラム

サラスの公式

行列の積の可換条件

CとFORTRANの計算速度はどちらが速いですか？

対称行列 対角行列

線形・行列の証明がさっぱり。。。

ユニタリ行列って？？

この行列問題を早くとく方法

線形代数 行列 対角化

fortran90

注目のQ&A

カテゴリ 一覧

専門家に質問してみよう 専門家登録

あなたにピッタリな商品が見つかる！ OKWAVE セレクト

行列　対角和　トレース

　Fortran90を授業でやっています。

対称行列　対角行列

線形代数　行列　対角化

カテゴリ
一覧

専門家に質問してみよう
専門家登録