yanzhengfft.m

来自「呵呵,第一次来很多东西不懂 还请多多指教啊! 几个编的简单的FFT和滤波器设计 」· M 代码 · 共 27 行

%%%对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信号')