求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
function r=cor_self(x,p,q)
%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
- m; c6 Z. Z  f- V$ ^" K! b" |clc;close all;
% j0 u4 `. j. S8 K% B" Wn=size(x,1);
6 |/ l0 e! o) a; H  V* yfor k=p:q
    sum=0;
    for i=1:n-k
" k& @: q  L2 d. }4 ]6 @9 P! p( v        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    end
. c' F1 l' W- @2 a0 _    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
& b9 j) e) r6 R8 B; ~- K* `end
k=[p:q];
plot(k,r,'-*');
9 \4 c4 `5 j9 }$ M. ^' I5 Xrefline(0,0);
. z3 `1 H/ z; \3 W) B/ J! \9 n7 Mend
运行显示
  Z  F- J, Y, x: e! a* qError using ==> plot
& X6 g" d( _" F2 q# _0 RVectors must be the same lengths.

Error in ==> cor_self at 13
plot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
; n7 x7 v8 t# G' Ck r 长度不一样   应该是r多了一位   可以进workspace里面看看

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

zjqylcy

我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   
5 Y! _: r; Z$ T; a0 ~  naf一般就是分母为0了