求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
) F% l. F8 D# _) C  O  Xfunction r=cor_self(x,p,q)
. o# A& V& u7 U  _1 z9 `% t0 M%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
& S' @( v5 e+ M7 k0 Q# l8 X6 sclc;close all;
$ \9 F1 z  k) En=size(x,1);
for k=p:q
  I: r5 q5 ~" h& J# e, o: L    sum=0;
, S/ v$ R: {5 Y) s4 D, O. ^    for i=1:n-k
0 v2 D& @. E' `2 _; [        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
7 E/ w3 T2 e" u6 o- e3 U$ E    end
    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
, Q) v  a6 x* Y! nend
) @+ p9 n! L5 A, [k=[p:q];
plot(k,r,'-*');
refline(0,0);
end
+ E# H# x0 u& v  Z) d" d2 e/ f2 l运行显示
Error using ==> plot
Vectors must be the same lengths.
# T! _" T9 ^: y9 ~
( N) P8 I1 P; oError in ==> cor_self at 13
# ]8 M/ z; r4 [; Y7 V6 Q& zplot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
k r 长度不一样   应该是r多了一位   可以进workspace里面看看

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

zjqylcy

我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   
& e2 G  ]( a6 Q: v4 }  naf一般就是分母为0了