求傅里叶级数系数出错，请高手帮忙看看

jww100

周期为0.1s，幅值为正负1的周期方波信号，
求它的傅里叶级数的系数
下面是我写的程序，。

clear all;
fs=1000;
N=1024;
delt=1/fs;
t=0:delt:(N-1)*delt;
T=0.1;
delf=1/T;
y=square(2*pi*t*delf);
f=(0:N-1)*fs/N;

syms z;
n=0:N-1;
an=int(y.*cos(2*pi*n*z/T),z,-T/2,T/2)/(T/2);
bn=int(y.*sin(2*pi*n*z/T),z,-T/2,T/2)/(T/2);
figure(3);
subplot(2,1,1);
stem(f,an);
subplot(2,1,2);
stem(f,bn);

[ 本帖最后由 jww100 于 2008-10-12 11:31 编辑 ]

ChaChing

1.楼主的方波y并非symbolic, 而是数值, int的结果应不是楼主要的
2.int积分後是symbolic, 而不是数值, 绘图前需先转换如stem(f,double(an))
方波如何用symbolic, 对不起, 目前我还懒得去学/试

jww100

谢谢chaching大哥指点