hc251.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第二章例2.5.1程序hc251
% 求卷积
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
x=input('输入x数祖 x=（例如ones(1,10)）  ');	% 输入x和h向量
h=input('输入h数祖 h=（例如exp(-0.1*[0:15])） ');	
Nx=length(x);Nh=length(h);L=Nx+Nh-1;	% 求向量长度
y=conv(x,h);				            % 用conv函数求卷积
subplot(3,1,1),stem(x)                  % 画出输入序列
vx=axis;axis([-2,L+2,vx(3),vx(4)])
title('输入序列'),ylabel('x(n)')
subplot(3,1,2),stem(h)                  % 画出脉冲响应序列
vh=axis;axis([-2,L+2,vh(3),vh(4)])
title('脉冲响应序列'),ylabel('h(n)')
subplot(3,1,3),stem(y)
vy=axis;axis([-2,L+2,vy(3),vy(4)])
title('输出序列'),ylabel('y(n)')        % 画出输出序列
set(gcf,'color','w')                    % 置图形背景色为白