两种不同方法计算自相关的结果疑问

re-us

我用下面的公式求滞后相关. 用autocorr和自己编的程序结果有差异. 请高手指点.
Tjan代表气温,样本是26.  我自己编的程序就是按照下面的公式.

clc; clear

Tjan=[2.7 5.9 3.4 4.7 3.8 5.3 5.0 4.3 5.7 3.6 3.1 3.9 3.0 4.9 5.7 4.8 5.6 6.4 5.6 4.2 4.9 4.1 4.2 3.3 3.7 7.6];

n=length(Tjan); s=cov(Tjan)*(n-1)/n; (这个就是为了和书本上一样才用(n-1)/n的)

ave=mean(Tjan); m=6;

for i=0:1:m, r(i+1)=sum((Tjan(1:n-i)-ave).*(Tjan(1+i:n)-ave))/s/(n-i); end

[ACF, Lags, Bounds] = autocorr(Tjan,6);

算出来的结果是: r = [1.0000, -0.0003, 0.0758, -0.0176, -0.2125, -0.2291, -0.4031]
用autocorr算出来的结果是: ACF' = [1.0000, -0.0003, 0.0700, -0.0155, -0.1798, -0.1851, -0.3101]

[ 本帖最后由 ChaChing 于 2010-2-3 09:39 编辑 ]

sigma665

精度问题吧,看上去差个0.1