matlab卷积

参考自：http://sepwww.stanford.edu/data/media/public/sep//prof/index.html
free book Earth Sounding Analysis Processing Vers.pdf
y0=x0b0;
y1=x1b0+x0b1;
y2=x2b0+x1b1+x0b2;
y3=x3b0+x2b1+x1*b2;
…
卷积定义：

function y = my_conv( x,h )
%
%   function y = my_conv( x,h )用来计算y(n) = h(n)*x(n)的卷积
nx = length(x); nh = length(h);
y = zeros(1,nx+nh-1); 

% for index = 1:nx
%     indexSum = x(index)*h;
%     y(1,index:index+nh-1) = y(1,index:index+nh-1)+indexSum; 
% end
for ib=1:nh
    for index = 1:nx
%         indexSum = x(index)*h;
%         y(1,index:index+nh-1) = y(1,index:index+nh-1)+indexSum; 
    y(index+ib-1)=y(index+ib-1)+x(index)*h(ib);
    end
end