频谱图的画法

1、本人原创：在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处，经过ffts

2、hift处理的蓝色频谱是正确的。注意：红色和蓝色的曲线在两边分别关于-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 t

3、he middle of the spectrum. For vectors, fftshift(X) swaps the left and right halves of X. 由此可见，fftshift的作用正是让正半轴部分和负半轴部分的图像分别关于各自的中心对称。将信号频率Fc改为100Hz后的频谱如下，蓝色是fftshift处理后的频谱：如何画一个信号的频谱今天终于搞明白了，这么简单的东西今天才明白如何快速的用matlab画出，真是惭愧。光啃书本知识，会做题目真的是不够的，以后我要多多注意如何实践！%比如clear;clc;t=0:0.001:2;n=2001;Fs=1000;Fc=2

4、00;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')  这里需要在做了2001点FFT后将0到2000映射到真正的频谱-fs/2到+fs/2上才能成为真正的频谱又如%  .正弦波定义f1=50;  %  频率1f2=100; %  频率2fs=2*(f1+f2);  %&

5、#160; 采样频率Ts=1/fs;  %  采样间隔N=120;    %  采样点数n=1:N;y=sin(2*pi*f1*n*Ts)+sin(2*pi*f2*n*Ts);  %  正弦波混合figure(1)subplot(2,1,1)plot(y);title('Signal')subplot(2,1,2)f=(1:120)*fs/200-fs/2;   %这里相当于做了120点FFT，需要把1：120映射到-fs/2到+

6、fs/2F=abs(fftshift(fft(y);stem(f,F);title('Amplitude Spectrum')%再比如：clc;clear all;fs=500;t=0:1/fs:1.5;f1=40; f2=20;x=sin(2*pi*t*f1)+sin(2*pi*t*f2); %input signalfigure(1);plot(t,x);title('input signal xn'); %input signal wavefigure(3);f=(0:750)*fs/751-fs/2;plot(f,abs(fftshift(fft(x); % plot input signal&#