MATLAB 快速傅里叶变换绘制频谱【fft函数】

fs=1000;
t = 0:1/fs:(length(data)-1)/fs;
MixSignal = data;
AmpFreRes = 20*log10(abs(fft(MixSignal))); % The amplitude frequency response(Power)

f = 0:fs/2/floor(length(AmpFreRes)/2):fs/2; % conver axis to Hz
p = AmpFreRes(1:length(f));
subplot(2,1,1)
plot(t,MixSignal);
xlabel('时间(s)');
ylabel('幅度(v)');
title('信号时域波形');
subplot(2,1,2)
plot(f,p);
xlabel('频率(f)');
ylabel('功率(dBW)');
title('信号频谱图');

感谢博主【通信牛肉干】的分享。

---

参考文献：

• https://blog.csdn.net/flypassion/article/details/106882812