如何使用matlab进行频域分析.docVIP

如何使用matlab进行频域分析 Matlab可以说是一个非常有用且功能齐全的工具，在通信、自控、金融等方面有广泛的应用。 本文讨论使用Matlab对信号进行频域分析的方法。 说到频域，不可避免的会提到傅里叶变换，傅里叶变换提供了一个将信号从时域转变到频域的方法。之所以要有信号的频域分析，是因为很多信号在时域不明显的特征可以在频域下得到很好的展现，可以更加容易的进行分析和处理。 FFT Matlab提供的傅里叶变换的函数是FFT，中文名叫做快速傅里叶变换。快速傅里叶变换的提出是伟大的，使得处理器处理数字信号的能力大大提升，也使我们生活向数字化迈了一大步。 接下来就谈谈如何使用这个函数。 fft使用很简单，但是一般信号都有x和y两个向量，而fft只会处理y向量，所以想让频域分析变得有意义，那么就需要用户自己处理x向量 一个简单的例子 从一个简单正弦信号开始吧，正弦信号定义为： y(t)=2sin?(2πf0t) 我们现在通过以下代码在Matlab中画出这个正弦曲线 1? fo = 4;?? %frequency of the sine wave2? Fs = 100; %sampling rate3? Ts = 1/Fs; %sampling time interval4? t = 0:Ts:1-Ts; %sampling period5? n = length(t); %number of samples6? y = 2*sin(2*pi*fo*t); %the sine curve78? %plot the cosine curve in the time domain9? sinePlot = figure;10? plot(t,y)11? xlabel(time (seconds))12? ylabel(y(t))13? title(Sample Sine Wave)14? grid 这就是我们得到的： 当我们对这条曲线fft时，我们希望在频域得到以下频谱（基于傅里叶变换理论，我们希望看见一个幅值为1的峰值在-4Hz处，另一个在+4Hz处） 使用FFT命令 我们知道目标是什么了，那么现在使用Matlab的内建的FFT函数来重新生成频谱 1? %plot the frequency spectrum using the MATLAB fft command2? matlabFFT = figure;? %create a new figure3? YfreqDomain = fft(y); %take the fft of our sin wave, y(t)45? stem(abs(YfreqDomain));? %use abs command to get the magnitude6? %similary, we would use angle command to get the phase plot!7? %well discuss phase in another post though!89? xlabel(Sample Number)10? ylabel(Amplitude)11? title(Using the Matlab fft command)12? grid13? axis([0,100,0,120]) 效果如下： 但是注意一下，这并不是我们真正想要的，有一些信息是缺失的 · x轴本来应该给我们提供频率信息，但是你能读出频率吗？· 幅度都是100· 没有让频谱中心为0 为FFT定义一个函数来获取双边频谱 以下代码可以简化获取双边频谱的过程，复制并保存到你的.m文件中 1? function [X,freq]=centeredFFT(x,Fs)2? %this is a custom function that helps in plotting the two-sided spectrum3? %x is the signal that is to be transformed4? %Fs is the sampling rate56? N=length(x);78? %this part of the code generates that frequency axis9? if mod(N,2)==010??? k=-N/2:N/2-1; % N even11? else12??? k=-(N-1)/2:(N-1)/2; % N odd13? end14? T=N/Fs;15? freq=k/T;? %the frequency axis1617? %takes the fft of the sign