求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
6 u. e/ m" x. c0 |  c. ifunction r=cor_self(x,p,q)
2 b, q8 w( O2 r  H  r" w6 \%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
clc;close all;
n=size(x,1);
! B7 \" R7 I/ ffor k=p:q
5 B* e+ v, i% M% ?, Q7 S5 P& q, Q: }+ _    sum=0;
    for i=1:n-k
6 u# i4 ?' D' {8 _- a        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
0 ?2 ~9 l, \' L/ C" c    end
1 ?9 u/ p. i4 p& K# i' S    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
end
  o/ ?( k4 y$ Y) j* U. U8 M; uk=[p:q];
plot(k,r,'-*');
refline(0,0);
end
运行显示
6 v7 `( U. c# q& h" K( TError using ==> plot
Vectors must be the same lengths.

Error in ==> cor_self at 13
plot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
1 J* q. l1 J, y/ |k r 长度不一样   应该是r多了一位   可以进workspace里面看看

' A+ {0 F0 B. F+ ], B4 ^& a% L这个后来我又试了下，是矩阵的问题，x的输入应该是列向量的形式，不过k的范围是怎么确定的呀，我输入的是从0到12，可是出来的结果第十一和第十二的数显示为naf，帮忙解释下吧，谢谢了~~

zjqylcy

我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   
* M& Y5 ^" D9 ^" [" s  B  o0 {  naf一般就是分母为0了