求助：绘制频谱图 有点急

努卡

如何绘制一周期为T,振幅为F0的矩形波的幅值频谱图和相位频谱图.
表达式为                 F0            0<t<1/2T
                     F(t)=   -F0          -1/2T<t<T  
希望哪位大哥给出程序 自己好琢磨

yangzj

照理说应该先滤波才能做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');

[ 本帖最后由 yangzj 于 2006-10-30 12:46 编辑 ]

努卡

不是很清楚每句语句的意思,能给出解释么?麻烦了

noble_herb

还没有找到答案，继续

suqiao

clear all;
fs=256; %采样频率
N=256; %采样点数
T=fs*N;%信号周期

F0=1;
s=[F0*ones(1,N/2) -F0*ones(1,N/2)]%前办个周期个数据为F0,后半个周期数据位-F0
sf=fft(s);%傅里叶变换
sf=sf(1:N/2)/N*2;%转换为实际的幅值

df=fs/N;%FFT频率间隔
f=(0:N/2-1)*df;%对应的实际频率

figure();%作图
subplot(211);
stem(f,abs(sf),'marker','none');%幅度谱

subplot(212);
stem(f,angle(sf),'marker','none');%相位谱

都是基本的，楼主多读些程序就好了~