Demo FFT dengan MATLAB

Membangun sinyal Sinusoidal murni sebesar : 5 Sin (2T)

(Amplitudo=5, periode (T)=1 sekon, Frekuensi = 1 Hz)

%

k=0:99;

T=0.02;

Tx=k*T;

N=100;

y=5*sin(2*pi*Tx);

plot(Tx,y);grid on;

%

Hasil plot :

Kemudian, sinyal ini kita transformasi ke domain frekuensi dengan FFT,

%

k=0:99;

T=0.02;

Tx=k*T;

N=100;

y=5*sin(2*pi*Tx);

subplot(2,1,1),plot(Tx,y);grid on;

hertz = k*(1/(N*T));

X=fft(y);

magX=abs(X)*2./N;

subplot(2,1,2),plot(hertz,magX);grid;

%

Hasil Plot:

Dari gambar 2, terlihat bahwa FFT melakukan perhitungan pada sinusoidal murni dengan menunjukkan sifat Simetri dan periodik.

Sinyal Sinusoidal yang digabung (dalam domain waktu) dengan beda waktu sebesar :

%

k=0:99;

T=0.02;

Tx=k*T;

N=100;

y=5*sin(2*pi*Tx)+(2*sin(2*pi*Tx)+pi);

subplot(2,1,1),plot(Tx,y);grid on;

hertz = k*(1/(N*T));

X=fft(y);

magX=abs(X)*2./N;

subplot(2,1,2),plot(hertz,magX);grid;

%

Hasil plot FFT akan terlihat bahwa disana terdapat 2 buah frekuensi dengan amplitudo tertentu, frekuensi tersebut adalah frekuensi dari sinyal sinus dengan beda rentang sebesar , ditambahkan pada sinyal pertama, hasil plot sebagai berikut :

Atau, jika sinyal awal dicampur dengan noise (sinyal acak):

%

k=0:99;

T=0.02;

Tx=k*T;

N=100;

y=5*sin(2*pi*Tx)+rand(1,length(k));

subplot(2,1,1),plot(Tx,y);grid on;

hertz = k*(1/(N*T));

X=fft(y);

magX=abs(X)*2./N;

subplot(2,1,2),plot(hertz,magX);grid;

%

Maka akan terlihat frekuensi lain muncul disekitar frekuensi sinyal awal.

Untuk hasil yang lebih baik, gunakan N yang lebih besar.

Selamt belajar :)