求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
: O9 T* a5 H! b4 N" b* jfunction r=cor_self(x,p,q)
* |4 ^' M  A# M8 m8 K%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
5 O3 J, I6 C" Y- Aclc;close all;
n=size(x,1);
$ c/ O7 a% y3 a: g& xfor k=p:q
    sum=0;
$ l- Q# X/ r( p) i$ M9 y    for i=1:n-k
        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    end
# F5 E' e& R9 `) @    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
end
  b0 B+ I- ~2 N7 ~; P' u7 Fk=[p:q];
: R9 }" G6 i; q$ z* V8 Eplot(k,r,'-*');
refline(0,0);
4 d# o" W8 u' K( L2 f4 a- Eend
; f2 N+ v% I0 ]+ b5 n运行显示
Error using ==> plot
Vectors 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
) z5 @& [$ r6 v: z& ^, O7 s% gk r 长度不一样   应该是r多了一位   可以进workspace里面看看

1 T, r" n% M# I4 |2 r2 O7 h) J这个后来我又试了下，是矩阵的问题，x的输入应该是列向量的形式，不过k的范围是怎么确定的呀，我输入的是从0到12，可是出来的结果第十一和第十二的数显示为naf，帮忙解释下吧，谢谢了~~

zjqylcy

我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   
* {) B3 l8 Z. z# i0 u  naf一般就是分母为0了