求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
function r=cor_self(x,p,q)
+ m( }" u: E" q! t* i%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
$ v7 D& Q  R/ O( ~9 I6 ^( S4 Zclc;close all;
n=size(x,1);
for k=p:q
' a$ C" [! ~; Y2 l  i. S, k% O    sum=0;
0 `5 q* c" ~: O2 z* r0 a    for i=1:n-k
        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    end
$ l* U6 p  g  ?    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);
  o3 Y' l. Y8 y* hend
7 c( i3 U8 T, P' x  T运行显示
Error using ==> plot
- [. s+ ^; r' R8 K1 g5 f& K2 H. w, fVectors must be the same lengths.

Error in ==> cor_self at 13
! G+ W% @6 u! G0 eplot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
( ~7 a5 T8 [; Z7 C4 i4 V8 Ek r 长度不一样   应该是r多了一位   可以进workspace里面看看

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

zjqylcy

我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   
2 p, ^) P- y# E4 r9 I2 `& x  naf一般就是分母为0了