求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
, s: x9 D9 K" k' D% c  y2 lfunction r=cor_self(x,p,q)
%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
5 l  C. `( ^, O; H9 R4 m1 i4 @clc;close all;
n=size(x,1);
, A+ |6 b! o" A5 z: |& pfor k=p:q
# A  h6 S( I! Q- Y  ]5 h    sum=0;
    for i=1:n-k
        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
: [8 M- v. W  B8 a3 C8 M- d    end
, u9 A: E, y! X0 U    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
end
: g1 u. Z) u: \% B9 E  [  _: nk=[p:q];
2 |6 {, _0 @- a3 ^plot(k,r,'-*');
. z( T; }: j7 @( t7 L7 Qrefline(0,0);
end
运行显示
Error using ==> plot
Vectors must be the same lengths.

7 C* v" K, y9 c$ X$ w( ?/ J, |Error in ==> cor_self at 13
4 t0 k& i8 n: Hplot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
2 w7 @4 E* x8 K  F4 u2 m' d! _# B$ i, Xk r 长度不一样   应该是r多了一位   可以进workspace里面看看

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

zjqylcy

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