ep341.m

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

%《数字信号处理--原理、实现及应用》第三章例3.4.1 重叠相加法的MATLAB实现程序：ep352.m
% 西安电子科技大学出版社出版  高西全  丁玉美合著 2008年5月
L=41;N=5;M=10;
hn=ones(1,N);hn1=[hn zeros(1,L-N)];	%产生h(n)，补零是为了绘图好看
n=0:L-1;
xn=cos(pi*n/10)+cos(2*pi*n/5);		%产生x(n)的L个样值
yn=fftfilt(hn,xn,M);		%调用fftfilt计算卷积
%===============================================================
%以下为绘图部分
subplot(4,1,1);stem(n,hn1,'.');ylabel('h(n)');axis([0,45,0,1.2])
subplot(4,1,2);stem(n,xn,'.');ylabel('x(n)')
xlabel('n');line([0,L+N-1],[0,0])
subplot(4,1,3);
m=0:length(yn)-1;
stem(m,yn,'.');line([0,L+N-1],[0,0])
xlabel('n');ylabel('y(n)');