%《MATLAB在電子信息課程中的應(yīng)用》第七章例7.15程序q715
% DFT計算
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

%程序q715: 
clear;close all
N=16;
N1=8;
%
%產(chǎn)生序列x1(n),計算DFT[x1(n)]
n=0:N-1;
x1n=exp(j*pi*n/8); %產(chǎn)生x1(n)
X1k=fft(x1n,N);   %計算N點DFT[x1(n)]
Xk1=fft(x1n,N1);  %計算N1點DFT[x1(n)]
%
%產(chǎn)生序列x2(n),計算DFT[x2(n)]
x2n=cos(pi*n/8);
X2k=fft(x2n,N);   %計算N點DFT[x2(n)]
Xk2=fft(x2n,N1);  %計算N1點DFT[x1(n)]
%
%產(chǎn)生序列x3(n),計算DFT[x3(n)]
x3n=sin(pi*n/8);
X3k=fft(x3n,N);   %計算N點DFT[x3(n)]
Xk3=fft(x3n,N1);  %計算N1點DFT[x1(n)]
%
%繪圖
subplot(3,3,1);stem(n,abs(X1k),'.');
title('16點 DFT[x1(n)]');
xlabel('k');ylabel('|X1(k)|')
subplot(3,3,2);stem(n,abs(X2k),'.');
title('16點 DFT[x3(n)]');
xlabel('k');ylabel('|X2(k)|')
subplot(3,3,3);stem(n,abs(X3k),'.');
title('16點 DFT[x3(n)]');
xlabel('k');ylabel('|X3(k)|')
k=0:N1-1;
subplot(3,3,7);stem(k,abs(Xk1),'.');
title('8點 DFT[x1(n)]');
xlabel('k');ylabel('|X1(k)|')
subplot(3,3,8);stem(k,abs(Xk2),'.');
title('8點 DFT[x2(n)]');
xlabel('k');ylabel('|X2(k)|')
subplot(3,3,9);stem(k,abs(Xk3),'.');
title('8點 DFT[x3(n)]');
xlabel('k');ylabel('|X3(k)|')