zxg.m

用Matlab编写的功率谱估计

function R=ZXG(x)
N=length(x);
X=zeros(N);
for m=1:N
    for n=1:N-m
        X(m)=X(m)+x(n)*x(n+m-1);
    end
    R(m)=(1/(N-m))*X(m);
end