求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
function r=cor_self(x,p,q)
  n' b& f5 L8 {, o% t1 z) h%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
; v9 P; T% e/ l3 r+ ]* oclc;close all;
n=size(x,1);
) m4 C3 R% k2 L8 n  u/ n4 Vfor k=p:q
    sum=0;
& D, Q; E( N; |% u( G. }    for i=1:n-k
+ d; T% K- ^$ H6 Y' g( B        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
1 e% T  k: z7 s; L    end
0 H, Q0 P6 z9 m! `. W    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
end
& B, p  p* @: S- Y. jk=[p:q];
plot(k,r,'-*');
refline(0,0);
. v/ E, M- J, G6 S; C+ D+ cend
运行显示
, r% ]. G* V8 V* T6 C9 bError using ==> plot
Vectors must be the same lengths.
/ K4 K4 N8 A0 N, h& r! Y8 B5 y; |- i
7 e+ R  s: m# F$ eError in ==> cor_self at 13
plot(k,r,'-*');

zjqylcy

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

梓爱

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

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

zjqylcy

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