Kode untuk Plot Respon Frekuensi

Saya gunakan kode berikut untuk membuat plot respon frekuensi dari data respon waktu. Plot dibuat dengan fungsi fft (salah satu algoritma transofrmasi fourier diskrit). Diasumsikan, data berupa matriks dua kolom; kolom 1 memuat t (waktu), dan kolom 2 memuat nilai y(t), dan periode sampling diketahui.

y=data(:,2); % input time-series signal
L=length(y);
Fs=1/Ts;     % Frequency sampling
NFFT = 2^nextpow2(L); % Next power of 2 from length of y

Y = fft(y,NFFT);
f = Fs/NFFT*(0:NFFT/2);
y = 2*abs(Y(1:NFFT/2));

plot(f(1:NFFT/2),y(1:NFFT/2));
title('Frequency Response');
xlabel('Frequency (Hz)');
ylabel(‘|Y(f)|’);

Ref: MATLAB R2008a Product Help