原创 绘制频谱图

http://www.chinavib.com/forum/thread-29926-1-121.html

如何绘制一周期为T,振幅为F0的矩形波的幅值频谱图和相位频谱图.
表达式为                 F0            0<t<1/2T
                     F(t)=   -F0          -1/2T<t<T 
-------------------------------------

照理说应该先滤波才能做FFT,但只要谱线数足够大的话频率反折现象也可忽略了.

clear all;
fs=256;
N=256;
T=fs*N;

F0=1;
s=[F0*ones(1,N/2) -F0*ones(1,N/2)];
sf=fft(s);
sf=sf(1:N/2)/N*2;

df=fs/N;
f=(0:N/2-1)*df;

figure();
subplot(211);
stem(f,abs(sf),'marker','none');

subplot(212);
stem(f,angle(sf),'marker','none');