conv_tp.m

来自「matlab例子」· M 代码 · 共 19 行

function [y,H]=conv_tp(h,x)
%--------------
%[y,H]=conv_tp(h,x)
%-----------------
nh=1:length(h);
n=length(h)+length(x)-1;
H=zeros(n,n);
[h1,nh1]=sigfold(h,nh);
for i=length(h)+1:n 
   h1=[h1,0];
end
h2=zeros(1,length(h1));
h2(1)=h1(1);
H=toeplitz(h2,h1);
H1=H(:,n+1-length(x));
for i=n+2-length(x):n 
   H1=[H1,H(:,i)];
end
y=(H1*x')';