ylabel('Magnitude');
title('N=1024 With Noise') grid
x=sin(2*pi*f1*t)+sin(2*pi*f2*t); y=fft(x,N); mag=abs(y);
f=(0:length(y)-1)'*fs/length(y); subplot(212)
plot(f(1:N/2),mag(1:N/2)); xlabel('Frequency(Hz)'); ylabel('Magnitude');
title('N=1024 Without Noise') grid
其频谱图如下：
N=1024 With Noise600Magnitude4002000050100150200250300350Frequency(Hz)N=1024 Without Noise400450500600400Magnitude2000050100150200250300Frequency(Hz)350400450500
四、实验环境
MATLAB 7.0
五、实验感想
通过了解FFT变换，及其库函数fft及ifft，更加明确了
运用快速傅里叶变换时点数N必须为2的整数次幂，另外通过实验截图可以看出随着取样点数不足N点的进行补零，超过N点的截取N点，并且取样点数越多，快速傅里叶变换后的幅值越大。另外当信号中存在造噪声时，可以通过对FFT后频谱的分析可以将周期信号从噪声中提取出来。