求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
$ ~5 v& z# n0 nfunction r=cor_self(x,p,q)
%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
1 R3 c$ s& j& ^' d) s) eclc;close all;
, B* y9 i7 \2 t+ v! F" r: Tn=size(x,1);
for k=p:q
. S2 K, J9 c  s1 c, k2 N+ z    sum=0;
    for i=1:n-k
4 b+ i1 |6 _; r+ N! x5 `        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
1 v& X6 {% ~$ \+ j4 J! O9 G% R$ d    end
2 w( ?, o* H/ W2 u; C& }# ^    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
end
k=[p:q];
plot(k,r,'-*');
refline(0,0);
) U: D, D8 {+ c# j2 @5 \% Bend
运行显示
Error using ==> plot
Vectors must be the same lengths.
2 `1 L; g# K# q% s
Error in ==> cor_self at 13
plot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
) M! ?/ y% A6 A+ F  s7 Bk r 长度不一样   应该是r多了一位   可以进workspace里面看看

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

zjqylcy

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