q715.m

MATLAB及在电子信息课程中的应用,电子工业出版社出版 陈怀琛 吴大正 高西全合著

%《MATLAB及在电子信息课程中的应用》第七章例7.15程序q715
% DFT计算
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
%程序q715: 
clear;close all
N=16;
N1=8;
%
%产生序列x1(n),计算DFT[x1(n)]
n=0:N-1;
x1n=exp(j*pi*n/8); %产生x1(n)
X1k=fft(x1n,N);                 % 计算N点DFT[x1(n)]
Xk1=fft(x1n,N1);                % 计算N1点DFT[x1(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)]
%
%产生序列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),'.');
set(gcf,'color','w')            % 设置图形背景色为白色
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)|')