Navigation :
Top/MATLAB演習問題:CSVデータの読み込みと書き出しとトレンド解析の基礎 問題
- バックアップ一覧
- 差分 を表示
- ソース を表示
- バックアップ を表示
- MATLAB演習問題:CSVデータの読み込みと書き出しとトレンド解析の基礎 問題 へ行く。
- 1 (2016-07-14 (木) 11:18:37)
- 2 (2016-07-15 (金) 09:36:37)
- 追加された行はこの色です。
- 削除された行はこの色です。
(2017年までに作成) *データの読み書き+フーリエ解析 [#p842e191] -.matデータ'crf_wind.mat'を読み込み,一番始めの行成分[u(1,:)]をプロット しなさい. uは風速データであり次元は[m/s],横軸は時間で増分はdt,次 元は[s]である. # load -textデータ'crf_wind.dat'を読み込み,一番始めの行成分[u(1,:)]をプロット し,.matデータのu(1,:)と比較しなさい.なお'crf_wind.dat'のデータ フォーマットは付録1である. # fopen, fscanf, fclose -u(1,:)のスペクトルを計算し,対数軸でプロットしなさい. # loglog -u(1,:)のスペクトルからべき乗則を求め,Kolmogorovスペクトルと比較しなさい. # loglog, polyfit + csv(コンマ区切り)データ'dat_climate.csv'を読み込んで下さい. -- 1行目は日時,2行めがデータです. # textscan, formatは'%s %f' + 読み込んだ時間データ(文字列)を数値に変換して下さい. # datenum + データを1日毎に補間して下さい. # interp1 + 補間データを元に,線形トレンドを求めて図にして下さい. # smooth, polyfit, polyval &ref(example_climatedata_1.jpg); - データをスペクトルを計算し,固有日数(1/f)を求めて下さい. # pwelch, semilogx &ref(example_climatedata_2.jpg); *データ [#md102557] *データ [#db9c1858] - [[テキストデータ>Home:files/education/crf_wind.dat]] - [[Matlabバイナリーデータ>Home:files/education/crf_wind.mat]] - CSVデータ:&ref(dat_climate.csv); 参考 % crf_wind.datの出力フォーマット filename_tmp = strcat('filename_w','.dat') fid = fopen( 'crf_wind.dat', 'w' ); fprintf( fid, '%5d %5d\n', nt, nz ); fprintf( fid, '%12.5e\n', dt ); fprintf( fid, '%12.5e', z(1:nz) ); fprintf( fid, '\n' ); for iz=1:nz fprintf( fid, '%12.5e', u(iz,:) ); fprintf( fid, '\n' ); end fclose(fid);