% 2021-08-13 Ver 1. % フーリエ変換するデーターを読込み、プロット F0 = importdata('FFTTest-01.txt'); plot (F0) % F0の1から300番目を取り出しプロット F0s = F0(1:300,:); plot (F0s) % フーリエ変換し、プロット t = 0.01; % サンプリング周期 l = 1024; % データー点数 F1 = fft(F0); % フーリエ変換 F2 = abs(F1/l); F3 = F2(1:l/2+1); ts = 1/t; tt = ts*(0:(l/2))/l; plot (tt,F3) % 横軸を拡大 xlim([0,20])