求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
function r=cor_self(x,p,q)
; t, m. O& C+ u/ Y( k% [%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
2 h6 {! y$ m! ]; ?! Lclc;close all;
n=size(x,1);
4 A6 _! d3 G2 b( b5 nfor k=p:q
    sum=0;
" C0 w% X( l6 h! E  L6 @4 D' u    for i=1:n-k
        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    end
9 X) @7 E0 a4 N    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
end
5 U! ~9 n; A% B* ik=[p:q];
plot(k,r,'-*');
. S; G4 Y/ m! H# Frefline(0,0);
0 I1 r' Y9 P/ S" n" _* m* zend
4 t$ v% K4 j' Z- w% y运行显示
Error using ==> plot
) x* z  h) X3 t/ ~' m; yVectors must be the same lengths.

4 Q6 ?! |. D- _7 H, w0 _2 {) L4 MError in ==> cor_self at 13
: o) Z& S" `0 Lplot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
  ^( v8 M" M8 C. B" @$ t2 pk r 长度不一样   应该是r多了一位   可以进workspace里面看看

5 M' X4 Y& ?7 [% p+ \, {* n# A这个后来我又试了下，是矩阵的问题，x的输入应该是列向量的形式，不过k的范围是怎么确定的呀，我输入的是从0到12，可是出来的结果第十一和第十二的数显示为naf，帮忙解释下吧，谢谢了~~

zjqylcy

我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   
" ^; I6 n* W4 D# p9 I3 U% r  naf一般就是分母为0了