%内容3.m
N1=16;N2=32,N3=64;Fs=64;T=1/Fs;
Tp1=T*N1;Tp2=T*N2;Tp3=T*N3;n1=0:N1-1;n2=0:N2-1;n3=0:N3-1;
x8t1=cos(8*pi*n1*T)+cos(16*pi*n1*T)+cos(20*pi*n1*T);
x8t2=cos(8*pi*n2*T)+cos(16*pi*n2*T)+cos(20*pi*n2*T);
x8t3=cos(8*pi*n3*T)+cos(16*pi*n3*T)+cos(20*pi*n3*T);
X8K1=fft(x8t1);X8K1=fftshift(X8K1);
X8K2=fft(x8t2);X8K2=fftshift(X8K2);
X8K3=fft(x8t3);X8K3=fftshift(X8K3);
% subplot(3,3,4);stem(X8K1,'.');xlabel('k');ylabel('X8K1(k)');
% subplot(3,3,5);stem(X8K2,'.');xlabel('k');ylabel('X8K2(k)');
% subplot(3,3,6);stem(X8K3,'.');xlabel('k');ylabel('X8K3(k)');
F1=1/Tp1;k=-N1/2:N1/2-1;f1k=k*F1;
subplot(3,1,1);stem(f1k,abs(X8K1),'.');title('16点|DFT[x_6(nT)]|');
xlabel('f(Hz)');ylabel('幅度');axis([-N1*F1/2-1,N1*F1/2-1,0,1.2*max(abs(X8K1))]);
F2=1/Tp2;k=-N2/2:N2/2-1;f2k=k*F2;
subplot(3,1,2);stem(f2k,abs(X8K2),'.');title('32点|DFT[x_6(nT)]|');
xlabel('f(Hz)');ylabel('幅度');axis([-N2*F2/2-1,N2*(F2)/2-1,0,1.2*max(abs(X8K2))]);
F3=1/Tp3;k=-N3/2:N3/2-1;f3k=k*F3;
subplot(3,1,3);stem(f3k,abs(X8K3),'.');title('64点|DFT[x_6(nT)]|');
xlabel('f(Hz)');ylabel('幅度');axis([-N3*F3/2-1,N3*F3/2-1,0,1.2*max(abs(X8K3))]);