求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
function r=cor_self(x,p,q)
; U2 ]" P5 l$ |' w/ Q% n%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
: S; b! l+ s0 C  U* W( \+ Aclc;close all;
. N, @2 M4 Y, e. }n=size(x,1);
for k=p:q
; E4 S1 G! |9 q: b; d1 b. Z1 ?    sum=0;
- |* g/ J! ^7 j4 f! M1 d  H2 k8 A    for i=1:n-k
2 b' |; B$ F$ n3 Q8 L* n  x        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);
) [. ^! y0 J4 [# nend
4 Q7 U5 q1 Z: {2 [& fk=[p:q];
plot(k,r,'-*');
3 z$ Q; t0 G' L( y! trefline(0,0);
end
1 I: w- L1 ]9 i# l4 N3 O运行显示
Error using ==> plot
; X0 `3 q2 J! l5 G# Y2 TVectors must be the same lengths.
/ U  J* n! n4 C- J# X
Error in ==> cor_self at 13
plot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
k r 长度不一样   应该是r多了一位   可以进workspace里面看看

3 {5 Y# k- \& g6 o9 ^( v+ ^+ F( `这个后来我又试了下，是矩阵的问题，x的输入应该是列向量的形式，不过k的范围是怎么确定的呀，我输入的是从0到12，可是出来的结果第十一和第十二的数显示为naf，帮忙解释下吧，谢谢了~~

zjqylcy

我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   
  naf一般就是分母为0了