求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
0 D  h6 d+ y. v7 \' z5 `, G3 hfunction r=cor_self(x,p,q)
) B& g& E7 i7 b; W%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
3 d+ Y; ~4 ~( I2 ~. jclc;close all;
n=size(x,1);
for k=p:q
    sum=0;
& `. `& i- y" c; U0 Q    for i=1:n-k
        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    end
! E. s5 q3 q6 Q) f, ]    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
end
8 v4 U2 g; Z: nk=[p:q];
8 S) U# e7 K, V% fplot(k,r,'-*');
8 S# o2 [% V1 M' y+ Orefline(0,0);
end
运行显示
' V0 X5 D9 H+ wError using ==> plot
Vectors must be the same lengths.

Error in ==> cor_self at 13
" v, s+ F, e( ^" Z  Zplot(k,r,'-*');

zjqylcy

k r 长度不一样   应该是r多了一位   可以进workspace里面看看

梓爱

zjqylcy 发表于 2012-7-13 16:31
0 w, q, K* l/ F8 ?: f! @9 V1 nk r 长度不一样   应该是r多了一位   可以进workspace里面看看

/ h0 V# b. B! @7 `4 ^/ w3 z这个后来我又试了下，是矩阵的问题，x的输入应该是列向量的形式，不过k的范围是怎么确定的呀，我输入的是从0到12，可是出来的结果第十一和第十二的数显示为naf，帮忙解释下吧，谢谢了~~

zjqylcy

我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   
  naf一般就是分母为0了