%《數字信號處理教程——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')                    % 置圖形背景色為白