rx_est.m

自己编的基于最抖下降法

function [Rx]=Rx_est(X,M)
% [Rx]=Rx_est(X,M)
% 函数估计随机过程的自相关系数
% 随机过程通过X输入，M是要求的自相关系数的个数
% 注意：本函数只会求取R(0),R(1),……R(M)
N=length(X);
Rx=zeros(1,M+1);    
for m=1:M+1, 
for n=1:N-m+1,
Rx(m)=Rx(m)+X(n)*X(n+m-1); 
end;
Rx(m)=Rx(m)/(N-m+1);      
end;