Takuya Miyashita
This site
Web
Start:
*Julia演習問題:フーリエ解析の基礎 [#ydc4cd76]
***問題 [#zfd69032]
+matファイル "crf_wind.mat" を読み込み,一番初めの行成分 ...
u は風速データであり単位は [m/s],2次元目は時系列を表し,...
#codeprettify(lang-julia){{
matdata = MAT.matopen() # v1.0.0 以降,正常にimportできな...
var_name = read(matdata,"var_name")
}}
+text データ "crf_wind.dat" を読み込み,一番初めの行成分 ...
#codeprettify(lang-julia){{
f = open(filename,"r")
readline(f) # or
readlines(f)
parse(Float64, stringdata[1])
}}
+u[1,:] のスペクトルを計算し,対数軸でプロットしなさい.~
#codeprettify(lang-julia){{
F0 = FFTW.fft(dataorg,1)
P = abs.(F0)
plt = Plots.plot(..., xscale=:log10, yscale=:log10)
}}
+u[1,:] のフーリエ解析結果から高周波成分を取り除きなさい...
#codeprettify(lang-julia){{
abs.(freq) .> freq_cutoff
FFTW.ifft(F0)
}}
***データ [#mff0554e]
-mat 形式ファイル~
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i...
-txt 形式ファイル
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i...
txt ファイルのフォーマットは以下の通り.~
1行目: nt, nz
2行目: dt
3行目: z(1:nz)
4行目以降:
u(1,1) u(1,2) ... u(1,nt)
u(2,1) u(2,2) ... u(2,nt)
...
u(nz,1) ... u(nz,nt)
***解答例(参考) [#m22bec5d]
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i...
↓↓↓の図,縦軸は powerではなくてフーリエ振幅 abs(F) です....
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i...
End:
*Julia演習問題:フーリエ解析の基礎 [#ydc4cd76]
***問題 [#zfd69032]
+matファイル "crf_wind.mat" を読み込み,一番初めの行成分 ...
u は風速データであり単位は [m/s],2次元目は時系列を表し,...
#codeprettify(lang-julia){{
matdata = MAT.matopen() # v1.0.0 以降,正常にimportできな...
var_name = read(matdata,"var_name")
}}
+text データ "crf_wind.dat" を読み込み,一番初めの行成分 ...
#codeprettify(lang-julia){{
f = open(filename,"r")
readline(f) # or
readlines(f)
parse(Float64, stringdata[1])
}}
+u[1,:] のスペクトルを計算し,対数軸でプロットしなさい.~
#codeprettify(lang-julia){{
F0 = FFTW.fft(dataorg,1)
P = abs.(F0)
plt = Plots.plot(..., xscale=:log10, yscale=:log10)
}}
+u[1,:] のフーリエ解析結果から高周波成分を取り除きなさい...
#codeprettify(lang-julia){{
abs.(freq) .> freq_cutoff
FFTW.ifft(F0)
}}
***データ [#mff0554e]
-mat 形式ファイル~
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i...
-txt 形式ファイル
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i...
txt ファイルのフォーマットは以下の通り.~
1行目: nt, nz
2行目: dt
3行目: z(1:nz)
4行目以降:
u(1,1) u(1,2) ... u(1,nt)
u(2,1) u(2,2) ... u(2,nt)
...
u(nz,1) ... u(nz,nt)
***解答例(参考) [#m22bec5d]
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i...
↓↓↓の図,縦軸は powerではなくてフーリエ振幅 abs(F) です....
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i...
Page:
Edit with a page name which already exists