clear all
n_max=[3 21 41 81];
N=length(n_max);
t=-1.1:.002:1.1;
omega_0=2*pi;
for k=1:N
n=[];
n=[1:n_max(k)];
b_n=2./(pi*n).*(-1).^(1+n);
x=b_n*sin(omega_0*n'*t);
subplot(N,1,k),plot(t,x,'linewidth',2);
axis([-1.1 1.1 -1.5 1.5]);
line([-1.1 1.1],[0 0],'color','r');
line([0 0],[-1.5 1.5],'color','r');
bt=strcat('最高谐波次数=',num2str(n_max(k)));
title(bt);
end