基于MATLAB分析,基于Matlab对信号进行频域分析的方法
描述
Matlab可以说是一个非常有用且功能齐全的工具,在通信、自控、金融等方面有广泛的应用。
本文讨论使用Matlab对信号进行频域分析的方法。
说到频域,不可避免的会提到傅里叶变换,傅里叶变换提供了一个将信号从时域转变到频域的方法。之所以要有信号的频域分析,是因为很多信号在时域不明显的特征可以在频域下得到很好的展现,可以更加容易的进行分析和处理。
FFT
Matlab提供的傅里叶变换的函数是FFT,中文名叫做快速傅里叶变换。快速傅里叶变换的提出是伟大的,使得处理器处理数字信号的能力大大提升,也使我们生活向数字化迈了一大步。
接下来就谈谈如何使用这个函数。
fft使用很简单,但是一般信号都有x和y两个向量,而fft只会处理y向量,所以想让频域分析变得有意义,那么就需要用户自己处理x向量
一个简单的例子
从一个简单正弦信号开始吧,正弦信号定义为:
我们现在通过以下代码在Matlab中画出这个正弦曲线
fo = 4; %frequency of the sine wave
Fs = 100; %sampling rate
Ts = 1/Fs; %sampling time interval
t = 0:Ts:1-Ts; %sampling period
n = length(t); %number of samples
y = 2*sin(2*pi*fo*t); %the sine curve
%plot the cosine curve in the time domain
sinePlot = figure;
plot(t,y)
xlabel(‘time (seconds)’)
ylabel(‘y(t)’)
title(‘Sample Sine Wave’)
grid
这就是我们得到的:
当我们对这条曲线fft时,我们希望在频域得到以下频谱(基于傅里叶变换理论,我们希望看见一个幅值为1的峰值在-4Hz处,另一个在+4Hz处)
使用FFT命令
我们知道目标是什么了,那么现在使用Matlab的内建的FFT函数来重新生成频谱
%plot the frequency spectrum using the MATLAB fft command
matlabFFT = figure; %create a new figure
YfreqDomain = fft(y); %take the fft of our sin wave, y(t)
stem(abs(YfreqDomain)); %use abs command to get the magnitude
%similary, we would use angle command to get the phase plot!
%we‘ll discuss phase in another post though!
xlabel(’Sample Number‘)
ylabel(’Amplitude‘)
title(’Using the Matlab fft command‘)
grid
axis([0,100,0,120])
效果如下:
但是注意一下,这并不是我们真正想要的,有一些信息是缺失的
x轴本来应该给我们提供频率信息,但是你能读出频率吗?
幅度都是100
没有让频谱中心为
为FFT定义一个函数来获取双边频谱
以下代码可以简化获取双边频谱的过程,复制并保存到你的.m文件中
function [X,freq]=centeredFFT(x,Fs)
%this is a custom function that helps in plotting the two-sided spectrum
%x is the signal that is to be transformed
%Fs is the sampling rate
N=length(x);
%this part of the code generates that frequency axis
if mod(N,2)==0
k=-N/2:N/2-1; % N even
else
k=-(N-1)/2:(N-1)/2; % N odd
end
T=N/Fs;
freq=k/T; %the frequency axis
%takes the fft of the signal, and adjusts the amplitude accordingly
X=fft(x)/N; % normalize the data
X=fftshift(X); %shifts the fft data so that it is centered
这个函数输出正确的频域范围和变换后的信号,它需要输入需要变换的信号和采样率。
接下来使用前文的正弦信号做一个简单的示例,注意你的示例.m文件要和centeredFFT.m文件在一个目录下
[YfreqDomain,frequencyRange] = centeredFFT(y,Fs);
centeredFFT = figure;
%remember to take the abs of YfreqDomain to get the magnitude!
stem(frequencyRange,abs(YfreqDomain));
xlabel(’Freq (Hz)‘)
ylabel(’Amplitude‘)
title(’Using the centeredFFT function‘)
grid
axis([-6,6,0,1.5])
效果如下:
这张图就满足了我们的需求,我们得到了在+4和-4处的峰值,而且幅值为1.
为FFT定义一个函数来获取右边频谱
从上图可以看出,FFT变换得到的频谱是左右对称的,因此,我们只需要其中一边就能获得信号的所有信息,我们一般保留正频率一侧。
以下的函数对上面的自定义函数做了一些修改,让它可以帮助我们只画出信号的正频率一侧
function [X,freq]=positiveFFT(x,Fs)
N=length(x); %get the number of points
k=0:N-1; %create a vector from 0 to N-1
T=N/Fs; %get the frequency interval
freq=k/T; %create the frequency range
X=fft(x)/N; % normalize the data
%only want the first half of the FFT, since it is redundant
cutOff = ceil(N/2);
%take only the first half of the spectrum
X = X(1:cutOff);
freq = freq(1:cutOff);
和前面一样,使用正弦信号做一个示例,下面是示例代码
[YfreqDomain,frequencyRange] = positiveFFT(y,Fs);
positiveFFT = figure;
stem(frequencyRange,abs(YfreqDomain));
set(positiveFFT,’Position‘,[500,500,500,300])
xlabel(’Freq (Hz)‘)
ylabel(’Amplitude‘)
title(’Using the positiveFFT function‘)
grid
axis([0,20,0,1.5])
效果如下:
打开APP精彩内容
点击阅读全文
基于MATLAB分析,基于Matlab对信号进行频域分析的方法相关推荐
- matlab的局放仿真,基于Matlab/Simulink的电力电缆局部放电信号传播的仿真方法
基于Matlab/Simulink的电力电缆局部放电信号传播的仿真方法 [技术领域] [0001] 本发明涉及一种电力电缆局部放电信号传播的仿真方法,尤其是涉及一种基于 Matlab/Simulink ...
- 2dpsk调制解调实验matlab_ila抓取数据,matlab分析,调试AD9361信号通路
关注.星标公众号,直达精彩内容 公众号:ZYNQ 作者:watchman AD9361功能介绍 上次我们介绍了AD9361的配置,代码,采样,脚本转换,软件使用等等,今天我们来调试并使用ila抓取实测 ...
- 扫描电镜图片的Matlab分析,基于Matlab的图像增强与复原技术在SEM图像中的应
0引 言 根据国内外的相关文献,研究和发展图像处理工具,改善图像质量是当今研究的热点.图像增强与复原是一种基本的图像处理技术.其按照特定的需要突出一幅图像中的某些信息或强化某些感兴趣的特征,将原来不 ...
- 使用matlab分析一段声音信号的均方根
修改记录:2022年1月15日17:37:18. 感谢me141414,W.Scofield,m0_57478585等几位好友,帮我发现提供的源代码中会提示"未定义变量youxiao&quo ...
- 期权组合策略用matlab分析,通过matlab 运用蒙特卡罗模型运算一个期权投资组合的VaR....
通过matlab 运用蒙特卡罗模型运算一个期权投资组合的VaR. 2016/03/01 18:22 通过matlab 运用蒙特卡罗模型运算一个期权投资组合的VaR. 题目: a bank has wr ...
- 基于Matlab App Designer的语音信号分析与处理(二):IIR和FIR滤波器的设计,语音信号的滤波
接上文:https://blog.csdn.net/weixin_53877178/article/details/122470759 目录 一.课题的任务 二.内容.步骤和要求 (1)语音信号的采集 ...
- 基于MATLAB的声信号的采集与分析,基于Matlab的声音信号采集与分析处理
2005年第6期 文章编号:100622475(2005)0620091202 计算机与现代化 JISUANJI YU XIANDAIHUA 总第118期 基于Matlab的声音信号采集与分析处理 陈 ...
- matlab随机信号分析,基于MATLAB的随机信号分析方法.ppt
<基于MATLAB的随机信号分析方法.ppt>由会员分享,可在线阅读,更多相关<基于MATLAB的随机信号分析方法.ppt(31页珍藏版)>请在人人文库网上搜索. 1.基于MA ...
- 连续函数matlab采样,基于 MATLAB 的时域信号采样及频谱分析(转)
一.课程设计题目: 基于 MATLAB 的时域信号采样及频谱分析 二.基本要求: ① 掌握数字信号处理的基本概念.基本理论和基本方法: ② 学会 MATLAB 的使用,掌握 MATLAB 的程序设计方 ...
- 基于matlab的语音信号去噪毕业论文,MATLAB的FIR数字滤波器语音信号的去噪研究和仿真...
摘要:语音信号作为人类语言交流的重要手段,语音处理的质量直接影响人们的正常通信.本文基于MATLAB对语音信号及加噪信号进行时域和频域分析,设计了FIR数字滤波器,完成对加噪信号滤波的处理.结果表明设 ...
最新文章
- 关于Jvm知识看这一篇就够了
- 静态html文件js读取url参数
- 如何构建分布式SFU/MCU媒体服务器?
- vue本地没事放到服务器上无限循环,解决vue中的无限循环问题
- 【计算机网络】Quiz集合
- python分支语句_Python中的分支语句和循环语句及案例
- POJ 1472 Instant Complexity
- 拓端tecdat|R语言如何做马尔可夫转换模型markov switching model
- 计算机设计大赛答辩提问,电子设计大赛答辩常见问题合集
- TextWatcher基本用法
- 计算机ping使用的端口,mac系统ping端口命令怎么使用
- codesmith oracle 64,CodeSmith OracleModel写法
- MAC-MAC-MAC-MAC
- 解压报错 你需要从上一压缩卷启动解压命令以便解压 解决方案及WinRAR怎么分卷压缩详解
- 必应每日一图php,自建必应每日一图API接口
- H.264 高度压缩数字视频编解码器标准
- centos代理服务器搭建(linux搭建http代理服务器)
- Java学习: wait和notify方法、ThreadLocal
- 社会工程学软件_信息工程学院大学生心理健康教育讲座
- concourse(concourse ci)