參考程式：

t=-pi:0.001:pi;

omega=1;

y=square(t,50);

plot(t,y),gridon

xlabel("t"),ylabel("週期方波訊號")

axis([-pipi-1.51.5])

n_max=[1:2:9];%表示不同的諧波數的疊加%

n=length(n_max);

fork=1:n

n=1:2:n_max(k);

b=4./(pi*n);

x=b*sin(omega*n"*t);

figure;

plot(t,y);

holdon;

plot(t,x);

holdoff;

xlabel("t"),ylabel("部分和的波形")

axis([-pipi-1.51.5]),gridon

title(["諧波數=",num2str(n_max(k))])

end

這個是透過諧波合成的過程，其中n為諧波數。當n改為單值99時可以看出，已經近乎方波，所以，可知，當n→∞時，可以合成為方波的。