q717.m

matlab在电子信息中的应用和实例中的源代码

%《MATLAB及在电子信息课程中的应用》第七章例7.17程序q717
% 时域与频域采样的对偶性验证
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
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,'.')
set(gcf,'color','w')                    % 设置图形背景色为白色
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,'.')
set(gcf,'color','w')                    % 设置图形背景色为白色
title('x1(n)的周期延拓序列')
xlabel('n');ylabel('x(mod(n,32))')