这段程序写的对吗，求指证

clear
T=1
t=linspace(0,T,1000);
triang=zeros(size(t));

max_harmonics=10;
a=zeros(1,max_harmonics+1);
b=zeros(1,max_harmonics+1);
synth_triang=zeros(size(t));
for n =0:max_harmonics
if mod (n,10)
b(n+1)=1/n;
synth_triang=synth_triang+b(n+1)sin(2pin(1/T)t);
end
end
t=linspace(0,10T,10*1000);
triang=repmat(triang,10,10);
synth_triang=repmat(synth_triang,10,10);
plot(t,triang,'b-',t,synth_triang,'r--')
legend('Triangular wave', 'Synthesised triangular wave')

f（t）=sinwt + （1/3）sin3wt +(1/5）sin5wt+（1/7）sin7wt +（1/9）sin9wt......
这个函数的傅里叶合成

然后想加一段显示频谱图的代码