%《數字信號處理教程——MATLAB釋義與實現》第二章例2.5.2程序hc252
% 用向量－toeplitz矩陣相乘求卷積
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
x=input('x=');h=input('h=');		% 輸入x和h向量
Nx=length(x);Nh=length(h);L=Nx+Nh-1;	% 求向量長度
% 生成toeplitz 矩陣H
H=toeplitz([h(1),zeros(1,Nx-1)],[h,zeros(1,Nx-1)]);
y=x*H				% 用向量－toeplitz矩陣相乘求卷積