function [c,lags]=vxcorr(a,b)
a=a(:);
b=b(:);
M=length(a);
maxlag=M-1;
lags=[-maxlag:maxlag]';
A=fft(a,2^nextpow2(2*M-1));
B=fft(b,2^nextpow2(2*M-1));
c=ifft(A.*conj(B));
c=[c(end-maxlag+1:end,1);c(1:maxlag+1,1)];
[nr nc]=size(a);
if(nr>nc)
    c=c.';
    lags=lags.';
end