Matlab使用fft画出信号频谱图的方法
更新时间：2021年04月23日 11:45:21 作者：非 常 道
这篇文章主要介绍了Matlab使用fft画出信号频谱图的方法,本文通过实例图文相结合给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
做雷原作业，需要对信号进行频谱分析，网上一搜太乱了，很多不是我想要的，特此整理。
DSP还没学到fft，就不对fft原理进行详细解释了，直接上代码。
fs=500;%采样率
f1=5;%信号频率
f2=10;%信号频率
T=1;%时宽1s
n=round(T*fs);%采样点个数
t=linspace(0,T,n);%时域横坐标
x = 3+cos(2*pi*f1*t) + 2.*cos(2*pi*f2*t);%形成三频信号,注意第二个频率信号幅度为2，直流幅度为3
figure(1);
plot(t,x);%画时域图
xlabel("t/s")
grid on
X = fftshift(fft(x./(n))); %用fft得出离散傅里叶变换
f=linspace(-fs/2,fs/2-1,n);%频域横坐标，注意奈奎斯特采样定理，最大原信号最大频率不超过采样频率的一半
figure(2)
plot(f,abs(X));%画双侧频谱幅度图
xlabel("f/Hz")
ylabel("幅度")
grid on
显然，该信号有三个频率分量，直流（0频），5Hz和10Hz，对应的幅度分别为3、1、2，其时域波形图如下:
转换为频域的关键函数是X = fftshift(fft(x./(n)));一定注意需要除以总样本数n，然后用fftshift将曲线挪一下位置。
做完这步操作后，得出的是双边频谱，频率范围从-fs/2到fs/2，这是因为奈奎斯特采样定理，给定采样频率为fs，那么原信号的最大频率不超过fs/2。横坐标的点数和时域信号的采样点数相同，这由离散傅里叶变换的性质给出。
绘制频谱幅度图如下：
可以看到，直流分量的幅度是3，这与时域的幅度相同，5Hz和10Hz的幅度分别为0.5和1，这比时域的幅度减小了一半，这是由于这个频谱图是双边频谱。沿着x=0线“对折”过去加起来，就变成单边频谱，频域幅度就和时域幅度对应上了。
到此这篇关于Matlab使用fft画出信号频谱的文章就介绍到这了,更多相关Matlab画fft信号频谱内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！
相关文章
这个网站的最大优点，就是省去了安装数学绘图软件或计算软件的麻烦，只要打开浏览器就能使用了。看了介绍之后，可别忘了把这个好网站加到书签2021-08-08
这篇文章主要介绍了在VScode上使用latex编辑pdf的相关知识，本文通过实例代码给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下2021-04-04
此错误仅适用于奔腾 处理器、奔腾处理器与MMX技术、奔腾OverDrive处理器 和奔腾OverDrive处理器带有MMX技术。 奔腾Pro处理器, 奔腾II处理器和i486及更早版本处理器不受影响2012-07-07
下面小编就为大家分享一篇Hbuilder连远程接服务器上传代码的图文教程，具有很好的参考价值，一起跟随小编过来看看吧，希望对大家有所帮助2017-11-11
FFmpeg是一套可以用来记录、转换数字音频、视频，并能将其转化为流的开源计算机程序。这篇文章主要介绍了ffmpeg播放器实现详解视频显示,需要的朋友可以参考下2020-07-07
本文主要初步带你了解代理和负载均衡的知识，文中对正向、反向代理以及反向代理与负载均衡的关系等做了详细讲解，感兴趣的朋友可以参考一下这篇文章2021-09-09
这篇文章主要介绍了密码哈希函数 Bcrypt的最大密码长度限制详解的相关资料,需要的朋友可以参考下2017-03-03
DLL(Dynamic Linkable Library) 详解说明
DLL文件(Dynamic Linkable Library 即动态链接库文件)，是一种不能单独运行的文件，它允许程序共享执行特殊任务所必需的代码和其他资源2008-12-12
这篇文章主要介绍了作为程序员必须了解的缩写和专业名词，文中讲解非常详细，对想学编程的朋友有所帮助，感兴趣的可以了解下2020-07-07
这篇文章主要介绍了Eclipse 误删文件的恢复，代码的恢复的相关资料,需要的朋友可以参考下2016-09-09
最新评论