function [X,eT]=convolve(A,B)
la=length(A);
lb=length(B);
N=la+lb-1;
X =zeros(1,la+lb-1);   % lunghezza del segnale risultante
B =fliplr(B);% ribalta il segnale B
XX=[zeros(1,N),A,zeros(1,N)];
HH=[B,zeros(1,length(XX)-length(B))];
HH=HH';
tic
for k=1:(la+lb-1)
    H=circshift(HH,(k-1));
    X(k)=XX*H;
end
eT=toc;
end
