求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
function r=cor_self(x,p,q)
5 W4 H4 p! l  A; {$ }( Z9 _%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
clc;close all;
/ j" _2 I4 g5 en=size(x,1);
for k=p:q
    sum=0;
+ W' a* ?3 {6 u: y7 _    for i=1:n-k
, `4 @" B! z: Q: Q2 d! r. v* g7 c        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    end
8 H% Z2 ^! F/ z9 J- U/ }! d+ I    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
' s& E$ M& s8 k/ w5 M6 Vend
# L) ~: U: n9 c2 D& l% I0 c. q/ qk=[p:q];
plot(k,r,'-*');
) z+ o( N" q) o; j% q% l9 K+ ~refline(0,0);
end
0 v9 i5 `: r& r% ^% E1 }& g运行显示
Error using ==> plot
Vectors must be the same lengths.
0 b5 A% M0 X8 O' \
! H) D, z- |8 U& t3 F/ f- tError in ==> cor_self at 13
: ^/ `; ~- U7 G3 \plot(k,r,'-*');

zjqylcy

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

梓爱

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

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

zjqylcy

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