Matlab使用fft画出信号频谱图的方法
更新时间：2021年04月23日 11:45:21 作者：非 常 道
这篇文章主要介绍了Matlab使用fft画出信号频谱图的方法,本文通过实例图文相结合给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
（福利推荐：【腾讯云】服务器最新限时优惠活动，云服务器1核2G仅99元/年、2核4G仅768元/3年，立即抢购>>>：9i0i.cn/qcloud）
（福利推荐：你还在原价购买阿里云服务器？现在阿里云0.8折限时抢购活动来啦！4核8G企业云服务器仅2998元/3年，立即抢购>>>：9i0i.cn/aliyun）
做雷原作业，需要对信号进行频谱分析，网上一搜太乱了，很多不是我想要的，特此整理。
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信号频谱内容请搜索程序员之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持程序员之家！
相关文章
这篇文章主要介绍了vscode配置备份的相关知识，本文通过实例代码相结合给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下2021-04-04
这篇文章主要介绍了Scala函数式编程的的相关资料，文中讲解非常细致，帮助大家更好的理解和学习，感兴趣的朋友可以了解下2020-06-06
这篇文章主要介绍了权限控制之粗粒度与细粒度概念及实现简单介绍，具有一定参考价值，需要的朋友可以了解下。2017-10-10
这篇文章主要介绍了作为程序员必须了解的缩写和专业名词，文中讲解非常详细，对想学编程的朋友有所帮助，感兴趣的可以了解下2020-07-07
B站是一个学习网站。一入B站深似海，从此游戏是路人。B站虽然视频资源多，但是内容也是五花八门、参差不齐，本文给大家收集了关于学习计算机编程的视频，这里有我曾经的入门视频，也有跟朋友交流获得的，特此筛选了下面这些视频，分享给大家2021-04-04
本文主要是为了方便线下播放Bilibili的弹幕，而专门制作的一款将弹幕转换为ASS的程序，介绍了程序制作的思路及过程，有需要的朋友可以参考下2014-09-09
向 KernelKiller 致敬，只需要背会以下2012-07-07
这篇文章主要介绍了scala基础语法的的相关资料，文中讲解非常详细，帮助大家更好的理解和学习，感兴趣的朋友可以了解下2020-06-06
对于这个算法，很多书上只是说一下思路，没有具体的实现。我在这里举个例子吧2012-04-04
国外的营销工作者和技术开发者开发了许多能帮助企业创建简单APP应用的工具，通过这些工具，就算是一个对程序一窍不通的普通网民都可以很容易的创建一个企业的APP程序，并可以对程序进行应用更新维护，开展营销等活动2013-11-11
最新评论