この問題、どのように解けばよいですか？

質問

中心極限定理のシュミレーションとして出された問題なのですが、乱数表(00から99まで)から無作為に数を抽出し、その数が00から49なら１、50から79なら２．80から99なら３に変換して、数の抽出を25回行い、変換した数25個の平均値を出すという作業を50回繰り返したものをグラフにするという問題なのですが、どのように解けばよいでしょうか？プログラムを使って作るようにいわれたのですが、プログラムに関してはほとんど知識はないので、(FORTRAN77ならかろうじて学校でやっているので少しはわかるのですが、もしできればFORTRAN77を使ったプログラムでお願いします)できるだけ丁寧に教えてください

回答No.1

wolv · ベストアンサー

プログラムの概要は
------------------------------------------------------------
for i=1,50
＃２５回の抽出作業
＃平均値の計算
＃平均値を出力
next i
------------------------------------------------------------
のようになります．＃部分は数行ずつ必要になります．
まず，「＃２５回の抽出作業」の部分は，
------------------------------------------------------------
sum=0
for j=1,25
＃数値の00から99の数値を無作為に抽出→xとする．
＃xを変換しyとする
sum=sum+y
next j
------------------------------------------------------------
となります．「＃平均値の計算」と「＃平均値の出力」は，
------------------------------------------------------------
ave=sum/25
write "F16.8",ave
------------------------------------------------------------
のようになります．write文の書式は記憶違いしているかもしれないので，
furenntirappuさんの知識をもとにして，必要なら修正してください．
さて，
＃数値の00から99の数値を無作為に抽出→xとする．
の部分ですが，
0から1の乱数を取り出す関数ってありませんか？忘れてしまいましたが．
その値を百倍し，整数部分だけとりだせば，00から99のうちの一つになります．

あとは，
if文で0<=x かつx<=49ならy=1
if文で50<=x かつx<=79ならy=2
if文で80<=x かつx<=99ならy=3
をそのまま書いてyに変換すれはよいです．
以上まとめると，
------------------------------------------------------------
for i=1,50

sum=0
for j=1,25

0から1の乱数を取り出し百倍して整数化する．
if文で0<=x かつx<=49ならy=1
if文で50<=x かつx<=79ならy=2
if文で80<=x かつx<=99ならy=3
sum=sum+y

next j
ave=sum/25
write "F16.8",ave

next i
------------------------------------------------------------
fortranの文法をよく覚えていないので，適当に修正してください．

グラフ化もfortranでやるのなら，「write "F16.8",ave」
の部分に少し工夫が必要です．

お礼

質問してから早い時間にありがとうございました。おかげでできました、さすがにグラフまではできませんでしたが。ちなみに少し修正させていただくとwrite(*,*)でforはdoになります。

回答No.2

wolv

グラフの表示は，例えば
　　　０　　　　１　　　　２　　　　３
　１：　　　　　　　＊
　２：　　　　　＊
　３：　　　　　　　　　＊
　　　　　：
２５： 　　　　　　　　＊

というように，折れ線グラフの線のないものを縦にしたようなグラフを書くのが
一番簡単です．
------------------------------------------------------------
write "F16.8",ave
------------------------------------------------------------
のかわりに，
------------------------------------------------------------
文字列のクリア（string=""など）
for k=1,int(ave*10)-1
文字列＝文字列＋”　”
next k
文字列＝文字列＋”＊”
write 文字列
------------------------------------------------------------
のようなことをやればグラフになります．

for文を使って，ave（の10倍）の長さに相当する文字列を作り，
それを表示してます．

特定の長さの文字列を生成する関数があったかもしれませんが，
自信がないので，あえて自分で作ってみました．

------------------------------------------------------------
この他に，
配列などを利用して度数分布を計算し，50回の作業が終わったところで
一気にグラフとして出力することも考えられます．