%%%對fft和ifft的驗證
clear all;
%產生一個正弦信號
N=256;
f1=2;
fs=20;
w=2*pi/fs;
n=0:N-1;
y=5*sin(w*f1*n);
%求fft
subplot(411)
stem(y(1:N/4));
ylabel('原始信號');
f=-0.5:1/N:0.5-1/N;
Y=fft(y);
t=fftshift(abs(Y));
subplot(412)
stem(f,t);
ylabel('FFT 頻譜');
subplot(413)
stem(abs(Y));
%求ifft
z=ifft(Y);
subplot(414)
z1=real(z);
stem(z1(1:N/4));
ylabel('ifft信號')