clc; clear all; pack;

　　Fs = 10e9;

　　Nsamp = 2e4;

　　t = [0:Nsamp-1].*(1/Fs);

　　f1 = 1e6;

　　x0 = 3.3/2;

　　x1 = x0 + 1.65*sin(2*pi*f1*t);

　　x3 = x0;

　　for n=1:2:3

　　x3 = x3 + 3.3*2/(pi*n) * sin(2*pi*n*f1*t);

　　end

　　x5 = x0;

　　for n=1:2:5

　　x5 = x5 + 3.3*2/(pi*n) * sin(2*pi*n*f1*t);

　　end

　　x7 = x0;

　　for n=1:2:7

　　x7 = x7 + 3.3*2/(pi*n) * sin(2*pi*n*f1*t);

　　end

　　figure

　　subplot(221)

　　plot(x1)

　　subplot(222)

　　plot(x3)

　　subplot(223)

　　plot(x5)

　　subplot(224)

　　plot(x7)

　　x217 = x0;

　　for n=1:2:217

　　x217 = x217 + 3.3*2/(pi*n) * sin(2*pi*n*f1*t);

　　end

　　figure

　　plot(x217)

　　figure

　　plot(x217,'k')

　　hold on

　　plot(x1,'b')

　　plot(x3,'g')

　　plot(x7,'r')

　　hold off

　　axis([8000 12000 -0.5 4])