unctionerr=fs(n) %n=100; %n为傅里叶级数 E=1;%E为方波信号的高度， T1=0.5;%T1为方波的周期; Fn=2000;%Fn为方波的采样频率 t0=T1/2;%t0为脉冲的宽度 T=2;%T为方波的截取的时间长度 t=-T/2:1/Fn:T/2;%t为横轴的时间坐标 squ=E*(1+square(2*pi/T1*t))/2;%产生方波脉冲信号，周期是T1，高度是E squ=cat(2,squ(T1/4*Fn:end),squ(1:T1/4*Fn-1));%把方波信号移位1/4周期，形成偶对称的脉冲