矩形波的幅频图怎么画啊

yuanping09

t=0:0.001:1;
x = square(2*pi*10*t,50);
plot(t,x)

大家帮帮忙啊，由这几句程序生成的矩形波，它的幅频图怎么画啊，程序怎么编写啊？怎么也做不出来了:@L ，请大侠们帮帮忙吧，谢谢了:handshake ！

ChaChing

幅频图的绘制与波形是什么有关吗?
建议动手搜下

yuanping09

tau_T=1/4;                                            % 占空比1/4
omega_0=2*pi;                                         % 基波频率
n0=0;n1=15;
n=n0:n1;
F_n=[(sin(tau_T*pi*n))/n*pi]*exp(-j*tau_T*n.*pi);
F_n=2*F_n;
Fn_max=max(abs(F_n));
Fn_min=min(abs(F_n));
subplot(2,1,1),stem(n,abs(F_n),'.');            
axis([n0 n1 Fn_min-0.1 Fn_max+0.1]);
line([n0 n1],[0 0],'color','r');                 
title('单边幅度频谱');   
请问一下，我这个程序怎么画出来幅度是一个值啊，哪里错了呢？

yuanping09

t = 0:1/1000:0.05;
y = square(2*pi*100*t,20);%f=100hz
Y=fft(y,512);%1024为采样点数
df = 1000/1024;%频率分辨率
f=df*(0:511)
n=1:512
plot(f,abs(Y))
ylable('magnitude')
xlabel('frequency (Hz)')
周期函数的频谱应该是离散的，我这个画出来的怎么不是离散的呢，还有幅值好像不对啊，到底哪里错了啊，初学，请帮帮忙吧！我上网查了可是还是解决不了。

ChaChing

个人水平专业有限, 建议楼主看看这里有没有你要问的问题－－信号处理专栏话题索引
http://forum.vibunion.com/forum/thread-49413-1-1.html