求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
% D, {( A* {1 ufunction r=cor_self(x,p,q)
%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
" o( ~' H& \6 k  Gclc;close all;
n=size(x,1);
for k=p:q
    sum=0;
    for i=1:n-k
. ]4 S# G0 Q# |        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    end
    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
' C9 Z9 m' b2 }9 ^end
- _5 w9 y- _9 l1 }k=[p:q];
: ^3 O( Z0 F- j( Qplot(k,r,'-*');
4 P* B& R/ ~) I1 f" w8 M  Vrefline(0,0);
* B- J  P2 }. \5 tend
/ j( m; C) [$ U3 a运行显示
8 ?) l  I2 q6 g4 y9 V( K! EError using ==> plot
% T6 W2 l& `' {  A+ b- j& @% FVectors must be the same lengths.

4 b' E: X! E/ y( BError in ==> cor_self at 13
8 |0 T4 U2 L7 e2 i2 t/ g5 D! Gplot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
+ \1 }( f1 w5 vk r 长度不一样   应该是r多了一位   可以进workspace里面看看

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

zjqylcy

我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   
# W* d$ s. l2 b6 \0 b3 }  naf一般就是分母为0了