Navigation :
Top/MATLAB演習問題:データの読み込みと書き出しとフーリエ解析 問題
MATLAB演習問題:データの読み込みと書き出しとフーリエ解析 問題
データの読み書き+フーリエ解析 †
- .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
データ †
参考
% 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);