求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
- i0 b0 U2 \  n7 ?function r=cor_self(x,p,q)
%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
+ Y7 F) v( q7 e" \  j5 s+ bclc;close all;
n=size(x,1);
  g0 k/ C6 E) a# @7 efor k=p:q
+ N6 ^+ {1 Y4 V1 Y" p8 d' }! @    sum=0;
    for i=1:n-k
        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    end
, E( Q' I7 m& o5 d# ^( b    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
* [' T% V& _9 g( |- x8 Cend
k=[p:q];
: f5 Z7 Q* n8 D/ N8 Lplot(k,r,'-*');
refline(0,0);
end
* _* I! ~1 ~# b! z9 {运行显示
7 w: k& J+ K- l) `2 G$ NError using ==> plot
Vectors must be the same lengths.
* P3 W7 Q0 F+ E6 t: F" S- u  r# p
+ D% p2 H+ [; z, M9 o4 PError in ==> cor_self at 13
plot(k,r,'-*');

zjqylcy

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

梓爱

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

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

zjqylcy

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