如何利用matlab合成周期方波信号

如何利用matlab合成周期方波信号

参考程序：
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→∞时，可以合成为方波的。