rxx.m

arma模型的功率谱估计代码。还包括系统的极零点分析等完整软件包。

function rxx = Rxx(x, k)
%估计X(n)的自相关Rxx(k)
k = abs(k);
N = length(x);
rxx = 0;
if k>=N
    rxx = 0;
    return;
end
for n=1:N-k
    rxx = rxx+x(n)*x(n+k);
end
%rxx = rxx/N;    %it's not right, just for test
rxx = rxx/(N-k);