r.m

利用一般最小二乘法和SVD-TLS方法 进行AR参数估计和正弦波频率估计

%计算自相关
function r=R(N,x,T)       %N is the data number,T is the time delay
r=sum(x((T+1):N).*x(1:(N-T)))/N;