q717.m

MATLAB及在电子信息课程中的应用例题索引与程序

%《MATLAB在电子信息课程中的应用》第七章例7.17程序q717
% 时域与频域采样的对偶性验证
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

clear;close all
M=40;N=64;n=0:M;
%产生M长三角波序列x(n)
xa=0:floor(M/2);  xb= ceil(M/2)-1:-1:0; xn=[xa,xb];
Xk=fft(xn,64);%64点FFT[x(n)]
X1k=Xk(1:2:N);%隔点抽取Xk得到X1(K)
x1n=ifft(X1k,N/2);%32点IFFT[X1(k)]得到x1(n)
nc=0:3*N/2;		% 取97点为观察区
xc=x1n(mod(nc,N/2)+1);%x1(n)的周期延拓序列
subplot(3,2,1);stem(n,xn,'.')
title('40点三角波序列x(n)');xlabel('n');ylabel('x(n)')
k=0:N-1;
subplot(3,2,5);stem(k,abs(Xk),'.')
title('64点DFT[x(n)]');xlabel('k');ylabel('|X(k)|')
k=0:N/2-1;
subplot(3,2,6);stem(k,abs(X1k),'.')
title('X(k)的隔点抽取');xlabel('k');ylabel('|X1(k)|')
n1=0:N/2-1;
subplot(3,2,2);stem(n1,x1n,'.')
title('32点IDFT[X1(k)]');xlabel('n');ylabel('x1(n)')
figure(2)
subplot(3,1,1);stem(nc,xc,'.')
title('x1(n)的周期延拓序列');xlabel('n');ylabel('x(mod(n,32))')