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的方波
“算不连续处合成波形(与理想波形相比)的超调百
分比” 这个应该怎么计算