努卡 发表于 2006-10-29 16:26

求助：绘制频谱图 有点急

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

yangzj 发表于 2006-10-29 19:21

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

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

F0=1;
s=;
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 编辑 ]

努卡 发表于 2006-10-30 09:09

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

noble_herb 发表于 2010-9-17 11:37

还没有找到答案，继续

suqiao 发表于 2010-9-27 11:14

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

F0=1;
s=%前办个周期个数据为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');%相位谱

都是基本的，楼主多读些程序就好了~
页: [1]
查看完整版本: 求助：绘制频谱图 有点急