求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
5 P+ Y! ?- ~3 N2 W$ _function r=cor_self(x,p,q)
%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
clc;close all;
6 z! f' b4 {3 n6 g4 \" \5 qn=size(x,1);
for k=p:q
    sum=0;
4 K$ l: `5 v8 d0 I    for i=1:n-k
$ U4 f: ~5 i5 A% T1 O- C- \' {        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
9 M: d: B; V% v7 }    end
    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
2 Q3 ^3 }  y6 D) s' Z2 \end
9 J" l8 ^# s) a  }. O! i$ z$ G+ sk=[p:q];
4 J0 b, g- e2 x; C$ @. pplot(k,r,'-*');
8 P5 r/ Q) o4 r- vrefline(0,0);
7 ~: g/ [8 A5 w; ]% Iend
运行显示
Error using ==> plot
Vectors must be the same lengths.
1 o! ]; R- x! q' F6 p$ Z+ ]
Error in ==> cor_self at 13
, x8 U6 n0 P8 ^" _- i8 d, Splot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
  K# I  S; c* i( ek r 长度不一样   应该是r多了一位   可以进workspace里面看看

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

zjqylcy

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