掲示板:数学月間(第17回)

10_フーリエ変換

投稿日時: 08/13 システム管理者

% 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])