vxcorr.m

来自「本程序仿真一个从卫星到区域数据中心的下行链路（由一组6个载波构成）」· M 代码 · 共 14 行

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.';
end