Navigation : Top / MATLAB演習問題:データの読み込みと書き出しとフーリエ解析の基礎 問題
*データの読み書き+フーリエ解析 [#h7c1d61a]

-.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 

*データ [#l2234b91]
*データ [#md0bced3]
- [[テキストデータ>Home:files/education/crf_wind.dat]]
- [[Matlabバイナリーデータ>Home:files/education/crf_wind.mat]]

参考
 % 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);