スペクトル解析実習

練習 4

答えのわかっているスペクトルを計算してみましょう (その 4).

ファイルから乱数を読んでスペクトルを計算してみましょう.

用意したデータファイルはこちら. (円周率の数字を 4 桁ずつ区切って生成した 0 から 10 までの実数の数列.)

データファイルは, wget では下のようにダウンロードできます.

$ wget http://itpass.scitec.kobe-u.ac.jp/~yot/itpass/2022_spectrum/data/pi.txt

円周率の値はこちらから取得しました.

上記ファイルをプログラムで読みましょう.

まず, 番号を読むための変数を定義します.

integer :: Num

その上で, 下のようにすればファイルから NData 個のデータを読むことができるでしょう.

open( 50, file='pi.txt', status='unknown' )
do j = 0, NData-1
  read( 50, * ) Num, x_Y(j)
end do
close( 50 )

NData は...

コンパイルして実行してみましょう.

どのようなスペクトルが計算されたか確認してみましょう.