スペクトル解析実習

練習 2

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

サンプルプログラムの入力波形を変更してスペクトルを計算してみましょう.

入力波形の変更

入力に波数 1, 7 の波を追加してみましょう.

do j = 0, NData-1
  x_Y(j) = 3.0d0 * cos( 2.0d0*PI/((DelX*dble(NData))/4.0d0)*x_X(j) - PI/3.0d0) &
    + 1.0d0 * cos( 2.0d0*PI/((DelX*dble(NData))/1.0d0)*x_X(j) ) &
    - 5.0d0 * cos( 2.0d0*PI/((DelX*dble(NData))/7.0d0)*x_X(j) - PI/2.0d0 )
end do

コンパイルと実行

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

結果の確認

予定通りのスペクトルが計算されたことを確認してみましょう.