k=8;n1=[0: 1: 19];
xa1=sin(2*pi*n1/k);
subplot(2,2,1)
plot(n1,xa1);
xlabel('t/T' );
xk1=fft(xa1);xk1=abs(xk1);
subplot(2,2,2)
stem(n1,xk1);
xlabel ('k') ;
n2=[0:1:15] ;
xa2=sin(2*pi*n2/k);
subplot(2,2,3)
plot(n2,xa2);
xlabel ('t/T' ) ;
xk2=fft(xa2);
xk2=abs(xk2);
subplot(2,2,4)
stem(n2,xk2)
xlabel('k');
ylabel('X(K)');
N=64;
n=[0:1:N-1];
n1=2*n;
n2=2*n+1;
k=[0: 1:N-1];
xn1=cos(2* pi/N*7*n1)+1/2*cos(2* pi/N*19*n1);
xn2=cos(2*pi/N*7*n2)+1/2*cos(2* pi/N*19*n2) ;
XK1=fft(xn1);
XK2=fft(xn2);
X1=XK1+exp(-j*pi*k/N).*XK2;
X2=XK1-exp(-j*pi*k/N) .*XK2;
X1=[X1 zeros(1,N)];
X2=[zeros(1,N) X2] ;
XK=X1+X2;
k=[0: 1:2*N-1] ;
XK=abs(XK);
stem(k,XK);
xlabel('k');
%ylabel (lX(k' )I ' );
title('X(k)=DFT[x(n) ]2N' );
N=64;k=[0:1:63];
xk=1./(1-0.8* exp(-j*2*pi*k/N));
xn=ifft(xk,64);
stem(k,xn);
xlabel('k') ;
ylabel('x(n)' );
disp('xn序列');
disp(xn);