这个应该如何修改，结果应该是方波

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

max_harmonics=10000;
a=zeros(1,max_harmonics+1);
b=zeros(1,max_harmonics+1);
synth_triang=zeros(size(t));
for n =0:max_harmonics
b(n+1)=(3/pi*(2n-1));
synth_triang=synth_triang+b(n+1)sin(2pi(2*n-1)*t/T);

end
t=linspace(0,10T,1010000);
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')
结果好像应该是1khz的方波
“算不连续处合成波形（与理想波形相比）的超调百

img


分比” 这个应该怎么计算

大玉米同学你好，给你更正如下，更改的地方都已经说明了

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

max_harmonics=10000;
a=zeros(1,max_harmonics+1);
b=zeros(1,max_harmonics+1);
synth_triang=zeros(size(t));
for n =1:max_harmonics %这里要从1开始
b(n+1)=3/(pi*(2*n-1)); % 这里2*n-1要做分母
synth_triang=synth_triang+b(n+1).*sin(2*pi*(2*n-1)*t/T);
end
t=linspace(0,10*T,10*10000);
triang=repmat(triang,1,10);% repmat应该是1行10列，不是10行10列
synth_triang=repmat(synth_triang,1,10); % repmat应该是1行10列
plot(t,triang,'b-',t,synth_triang,'r--')
legend('Triangular wave', 'Synthesised triangular wave')

img

img