求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
, }7 v' Q, Y5 i7 [function r=cor_self(x,p,q)
%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
) x# l+ R6 _/ J/ pclc;close all;
# W8 ]7 L; g3 S: K, |n=size(x,1);
( E  b' Y6 e; d7 [' _for k=p:q
    sum=0;
    for i=1:n-k
" {" Q+ }( V5 d1 @0 S# G$ X5 W1 o& r        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
4 Y& G) f" ^) v3 f9 Y' B    end
    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
end
* b; o/ U) r( `) B2 m- ?( jk=[p:q];
6 y/ |' q& P& s' L2 splot(k,r,'-*');
; F1 Y+ R5 X; _9 Jrefline(0,0);
end
运行显示
, j, [! `- Q2 y$ H7 zError using ==> plot
Vectors must be the same lengths.

Error in ==> cor_self at 13
& J/ W1 d( T+ x. c  Oplot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
' c: V+ x" l7 ^; v4 W; Ok r 长度不一样   应该是r多了一位   可以进workspace里面看看

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

zjqylcy

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