vxcorr.m

来自「QPSK系统在无线通信信道中的仿真,通过对有若干条（大于2条）路径的无线多径衰落」· M 代码 · 共 15 行

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