ep331.m

高西全2009年最新版<数字信号处理>源码程序

%《数字信号处理(第三版）》第三章例3.3.1程序ep331.m
% 频域采样理论验证
M=26;N=32;n=0:M;
%产生M长三角波序列x(n)
xa=0:M/2;  xb= ceil(M/2)-1:-1:0; xn=[xa,xb];
Xk=fft(xn,512);%512点FFT[x(n)]
X32k=fft(xn,32);%32点FFT[x(n)]
x32n=ifft(X32k);%32点IFFT[X32(k)]得到x32(n)
X16k=X32k(1:2:N);%隔点抽取X32k得到X16(K)
x16n=ifft(X16k,N/2);%16点IFFT[X16(k)]得到x16(n)
subplot(3,2,2);stem(n,xn,'.')
title('(b) 三角波序列x(n)');xlabel('n');ylabel('x(n)');axis([0,32,0,15])
k=0:511;wk=2*k/512;
subplot(3,2,1);plot(wk,abs(Xk));title('(a)FT[x(n)]');
xlabel('\omega/\pi');ylabel('|X(e^j^\omega)|');axis([0,1,0,200])
k=0:N/2-1;
subplot(3,2,3);stem(k,abs(X16k),'.')
title('(c) 16点频域采样');xlabel('k');ylabel('|X_1_6(k)|');axis([0,8,0,200])
n1=0:N/2-1;
subplot(3,2,4);stem(n1,x16n,'.');
title('(d) 16点IDFT[X_1_6(k)]');xlabel('n');ylabel('x_1_6(n)');axis([0,32,0,15])
k=0:N-1;
subplot(3,2,5);stem(k,abs(X32k),'.')
title('(e) 32点频域采样');xlabel('k');ylabel('|X_3_2(k)|');axis([0,16,0,200])
n1=0:N-1;
subplot(3,2,6);stem(n1,x32n,'.');
title('(f) 32点IDFT[X_3_2(k)]');xlabel('n');ylabel('x_3_2(n)');axis([0,32,0,15])