求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
function r=cor_self(x,p,q)
%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
clc;close all;
% e2 p. I1 {, V0 g7 k; hn=size(x,1);
% ~2 ?5 |: F+ Q3 ofor k=p:q
! o- Z" D, x$ U8 Q    sum=0;
8 K# c3 C2 q" h3 x, F/ e# M    for i=1:n-k
+ S" o8 a1 ^/ Y8 x' b, O, Z        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    end
8 u6 w$ D/ U9 y# Z+ ?    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
end
k=[p:q];
plot(k,r,'-*');
) Y& d; u: A$ A0 _refline(0,0);
end
运行显示
, |, G6 D2 r) E' b  E3 u  JError using ==> plot
  [. y# q2 j  Z' s( X$ z: |Vectors must be the same lengths.
9 j; N' n! t. n* I8 ^6 U
Error in ==> cor_self at 13
plot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
9 y. p/ b" y) G) j* Nk r 长度不一样   应该是r多了一位   可以进workspace里面看看

" s  @6 x! L% E  A' g这个后来我又试了下，是矩阵的问题，x的输入应该是列向量的形式，不过k的范围是怎么确定的呀，我输入的是从0到12，可是出来的结果第十一和第十二的数显示为naf，帮忙解释下吧，谢谢了~~

zjqylcy

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