ep321.m

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

%《数字信号处理(第三版）》第三章例3.2.1程序ep321.m
% 西安电子科技大学出版社出版  高西全  丁玉美 合著 2008年
% 循环卷积计算（为了编程简单，采用DFT计算）
hn=[1 1 1 1]; xn=[1 2 3 4];
X4k=fft(xn);        %4点FFT[x(n)]
X8k=fft(xn,8);      %8点FFT[x(n)]
H4k=fft(hn);        %4点FFT[h(n)]
H8k=fft(hn,8);      %8点FFT[h(n)]
Y4k=X4k.*H4k;y4n=ifft(Y4k);
Y8k=X8k.*H8k;y8n=ifft(Y8k);
n=0:3;
subplot(3,2,1);stem(n,xn,'.')
title('(a) x(n)');xlabel('n');ylabel('x(n)');axis([0,8,0,5])
subplot(3,2,2);stem(n,hn,'.')
title('(b) h(n)');xlabel('n');ylabel('h(n)');axis([0,8,0,1.5])
subplot(3,2,5);stem(n,y4n,'.')
title('(c) 4点循环卷积');xlabel('n');ylabel('yc(n)');axis([0,8,0,12])
n=0:7;
subplot(3,2,6);stem(n,y8n,'.')
title('(d) 8点循环卷积');xlabel('n');ylabel('yc(n)');axis([0,8,0,12])