参考程序：
t=-pi:0.001:pi;
omega=1;
y=square(t,50);
plot(t,y),grid on
xlabel('t'),ylabel('周期方波信号')
axis([-pi pi -1.5 1.5])
n_max=[1:2:9]; %表示不同的谐波数的叠加%
N=length(n_max);
for k=1:N
n=1:2:n_max(k);
b=4./(pi*n);
x=b*sin(omega*n'*t);
figure;
plot(t,y);
hold on;
plot(t,x);
hold off;
xlabel('t'),ylabel('部分和的波形')
axis([-pi pi -1.5 1.5]),grid on
title(['谐波数=',num2str(n_max(k))])
end
这个是通过谐波合成的过程,其中 n 为谐波数.当n改为单值99时可以看出,已经近乎方波,所以,可知,当n→∞时,可以合成为方波的.
