求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
function r=cor_self(x,p,q)
%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
$ T+ Z& R0 K9 a. Z8 V: e+ Xclc;close all;
- o! \5 b$ Q. L8 ~8 yn=size(x,1);
for k=p:q
    sum=0;
" |, a. X/ u# D9 R& @    for i=1:n-k
: o4 R: v4 F& v6 C7 ]        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
4 L$ [2 Y8 W2 ?$ d- v    end
    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
9 c5 _0 Q) ]+ Wend
k=[p:q];
/ F% @$ G8 j; z5 @9 {/ k9 U+ Hplot(k,r,'-*');
4 `5 N  l% L; w* Erefline(0,0);
1 w* o8 K& R, }9 H8 fend
运行显示
Error using ==> plot
& |- _# \! `7 f* tVectors must be the same lengths.
6 j) P# U: @' r. [& k& s+ q9 l
2 A0 i! @. B  m0 D7 Q) kError in ==> cor_self at 13
" w" A5 i6 d( y; wplot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
; ]* A( Q3 ]& n: |k r 长度不一样   应该是r多了一位   可以进workspace里面看看

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

zjqylcy

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