WKT_AC < test.dat

#---------------------------------------------------------------------------

# ウィーナー・ヒンチンの定理

#---------------------------------------------------------------------------

# パワースペクトルを逆フーリエ変換し、自己相関関数を計算

#  WKT ACF: パワースペクトルから求めた自己相関関数（ウィーナー・ヒンチンの定理）

# DATA ACF: データから直接求めた自己相関関数

#      gap: WKT ACF - DATA ACF

#---------------------------------------------------------------------------

#---------------------------------------------------------------------------

# Wiener-Khinchin theorem (using non-cyclic data)

#---------------------------------------------------------------------------

# inverse Fourier transform of the power spectrum and calculation of the auto-correlation function

#  WKT ACF: auto-correlation function obtained from the power spectrum (Wiener-Khinchin theorem)

# DATA ACF: auto-correlation function obtained directly from the data

#      gap: WKT ACF - DATA ACF

#---------------------------------------------------------------------------

# lag   WKT ACF   DATA ACF       gap

   0   4.756758   4.756758   -7.10543e-015

   1   2.892687   2.892687   -2.66454e-015

   2   3.289990   3.289990   -3.10862e-015

   3   2.400743   2.400743   -3.9968e-015

   4   2.649403   2.649403   -3.9968e-015

   5   2.556849   2.556849   -4.88498e-015

   6   2.677648   2.677648   -4.88498e-015

   7   2.185213   2.185213   -4.44089e-015

   8   2.058504   2.058504   -2.66454e-015

   9   1.086492   1.086492   2.22045e-016

   10   1.216006   1.216006   2.22045e-016

   11   0.697044   0.697044   -3.33067e-015

   12   1.246480   1.246480   -3.10862e-015

   13   0.526165   0.526165   5.21805e-015

   14   0.680647   0.680647   5.55112e-016

   15   0.133933   0.133933   3.05311e-016

#---------------------------------------------------------------------------

# 自己相関関数をフーリエ変換し、パワースペクトルを計算

# DATA PS: データから直接求めたパワースペクトル

#  WKT PS: 自己相関関数から求めたパワースペクトル（ウィーナー・ヒンチンの定理）
#     gap: DATA PS - WKT PS

#---------------------------------------------------------------------------

# Fourier transform the autocorrelation function and calculate the power spectrum

# DATA PS: power spectrum obtained directly from the data

#  WKT PS: power spectrum obtained from the auto-correlation function (Wiener-Khinchin theorem)

#     gap: DATA PS - WKT PS

#---------------------------------------------------------------------------

# freq     DATA PS   WKT PS       gap

   0     57.352366   57.352366   -2.13163e-014

   1     23.704925   23.704925   7.10543e-015

   2     0.045497   0.045497   3.44863e-015

   3     1.056312   1.056312   6.43929e-015

   4     1.248207   1.248207   1.68754e-014

   5     5.417398   5.417398   4.44089e-015

   6     0.947489   0.947489   2.04281e-014

   7     0.731677   0.731677   -1.88738e-015

   8     0.936949   0.936949   2.68674e-014

   9     0.869644   0.869644   -4.44089e-015

   10     0.114270   0.114270   1.02279e-014

   11     1.535367   1.535367   -3.33067e-015

   12     0.915796   0.915796   3.0087e-014

   13     3.081823   3.081823   -1.77636e-014

   14     1.451739   1.451739   -6.88338e-015

   15     1.656916   1.656916   -1.06581e-014