求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
* H8 ~$ m' w9 P0 N( n# u* R6 E" {5 u+ d- ofunction r=cor_self(x,p,q)
4 ~( q+ }8 E: R# x( G! ^%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
* A4 C* e5 {( L# Dclc;close all;
  e# x" F9 i7 A, w& e8 r* z4 I  Yn=size(x,1);
for k=p:q
6 a, t; q9 k7 y8 k    sum=0;
    for i=1:n-k
        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    end
0 f5 M+ O" M. K' h8 l6 t2 `    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
end
k=[p:q];
$ o, y- e* j' [7 yplot(k,r,'-*');
refline(0,0);
3 O  _/ L; N" e6 f7 Bend
  D3 G3 J& M, J3 D$ E运行显示
+ I- q% k: E9 WError using ==> plot
% S9 R9 j* f. e* n# z4 F; QVectors must be the same lengths.
6 N1 c" p$ |1 a3 g' e% G! l/ q! b$ r
. g( k5 x' @. H6 uError in ==> cor_self at 13
plot(k,r,'-*');

梓爱

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

  _# n7 D1 j5 s+ q这个后来我又试了下，是矩阵的问题，x的输入应该是列向量的形式，不过k的范围是怎么确定的呀，我输入的是从0到12，可是出来的结果第十一和第十二的数显示为naf，帮忙解释下吧，谢谢了~~