covariance.m

来自「频谱感知 多用户合作 能量检测 伪随机序列」· M 代码 · 共 19 行

function coresult=covariance(x,N,tal)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   x   采样信号样本
%   N   周期所含样本数
%   tal 延迟
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
T=length(x);
S=floor(T/N);
y=[x,zeros(1,T)]
Rc(N,N)=0;
for m=1:1:N
    for n=1:1:N
        for s=0:1:S-1
        Rc(m,n)=Rc(m,n)+y(n+s*N)*y(n+s*N+tal)*y(m+s*N)*y(m+s*N+tal);
        end
    end
end
coresult=Rc/S
coresult=inv(coresult)