matlab中应用fft求傅立叶变换并画频谱图

kongfeng0824

在matlab中应用fft求傅立叶变换后，如果想画出频谱图，必须用fftshift命令处理变换的结果。例子如下：

clear;

clc;

t=0:0.001:2;

n=2001;

Fs=1000;

Fc=200;

x=cos(2*pi*Fc*t);

y1=fft(x);

y2=fftshift(y1);

f=(0:2000)*Fs/n-Fs/2;

hold on;

plot(f,abs(y1),'r')

plot(f,abs(y2),'b')

结果如下图：

图中红色是没经过fftshift处理的频谱图，蓝色是经过处理之后的。结合程序，显然x的频谱应该位于200Hz处，经过fftshift处理的蓝色频谱是正确的。注意：红色和蓝色的曲线在两边分别关于-250Hz和250Hz对称。这并不是偶然。以下是Matlab的帮助文件中对fftshift的说明：Y = fftshift(X) rearranges the outputs of fft, fft2, and fftn by moving the zero-frequency component to the center of the array. It is useful for visualizing a Fourier transform with the zero-frequency component in the middle of the spectrum. For vectors, fftshift(X) swaps the left and right halves of X.由此可见，fftshift的作用正是让正半轴部分和负半轴部分的图像分别关于各自的中心对称。将信号频率Fc改为100Hz后的频谱如下，蓝色是fftshift处理后的频谱：

棋子灯花

请教一下楼主，程序里面的t,n,fs,fc分别是什么意义呢？

福剑in奋斗

太给力了，非常喜欢，收下了

kongfeng0824

福剑in奋斗 发表于 2013-7-29 12:29

                               
登录/注册后可看大图

太给力了，非常喜欢，收下了

有用就好      

elunxu

非常感谢，学习学习

ghgu

谢谢！非常感谢!

蓝冰草

试了一下，果然。不过这其中有什么物理意义吗？纯数学处理？

navy_yehaijun

学习了，好东西

棋子灯花

其实我就是想知道怎么确定的，采样频率什么的

小虫小虫

{:eb513:}学习了。。。