求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
, A2 Z: V; e  z  z; B1 S2 mfunction r=cor_self(x,p,q)
%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
  j8 V0 M! n7 P( N! Q! l/ w* Gclc;close all;
8 S& m% f& \- w& `2 Gn=size(x,1);
: H# k/ z( Q, `: gfor k=p:q
    sum=0;
    for i=1:n-k
& l& w/ j+ \! m- T        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    end
6 k3 I1 M" `3 N    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
end
' }" R8 o  D! x- A& L' zk=[p:q];
plot(k,r,'-*');
" `2 m3 f9 b* n1 U9 A3 q" trefline(0,0);
end
# {) e1 T6 F1 y+ ?6 \# r+ v运行显示
8 q3 t3 j  V+ jError using ==> plot
& n, D% C) P- K: R* s# ?  B% j3 i' ^1 ^Vectors must be the same lengths.
* q1 x. O4 H/ X8 S9 z
Error in ==> cor_self at 13
- c. s$ G2 o! g* E8 N/ bplot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
: ]( W! a7 _# Y) ]( E: O. P5 ]k r 长度不一样   应该是r多了一位   可以进workspace里面看看

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

zjqylcy

我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   
( O1 X8 `. [9 ^8 @- p0 Z6 g  naf一般就是分母为0了