求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
function r=cor_self(x,p,q)
1 o6 I+ T: P) U0 e%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
; m# E) j0 S( p3 p8 P8 i9 Uclc;close all;
n=size(x,1);
for k=p:q
/ J3 ^) z3 r, a, X5 [4 Y3 V; G+ b( s    sum=0;
  l& z. x8 C  b' i+ |    for i=1:n-k
        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
9 r" s+ N' r  G' r( E5 C    end
4 ^0 U( {% r6 p! E" t* F8 P    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
! T& {2 y! T1 B. @, g! M! t- V9 W  Rend
; Z+ R, r3 h& N7 J" Y0 Lk=[p:q];
0 b' m7 R* e' S8 B7 W9 `plot(k,r,'-*');
refline(0,0);
- ^; B, O* ]9 R7 y- o; X' Vend
运行显示
Error using ==> plot
Vectors must be the same lengths.

Error in ==> cor_self at 13
- q( E) v3 x' {8 Yplot(k,r,'-*');

zjqylcy

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

梓爱

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

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

zjqylcy

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