求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
* b' x4 w, T0 {; d5 Z9 z/ }6 s8 D1 ]function r=cor_self(x,p,q)
" h6 Q: d( D3 X* \. R%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
- A& V. Y- H( {  t' }clc;close all;
n=size(x,1);
for k=p:q
* J' D/ K4 Y& ?$ a  I    sum=0;
) }) A: R. H4 M& S% O    for i=1:n-k
; q. U9 A' E6 q- `# V" I" T        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
" g/ v+ K. G1 r8 ]; |$ S) y2 p    end
8 X, h- P  o* A    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
9 @0 x2 V3 w# v( t3 k0 ~end
k=[p:q];
; a6 @# E, ^( J  M7 `plot(k,r,'-*');
; E5 Y5 M) G- C3 \# s+ trefline(0,0);
end
运行显示
* K3 K4 }: j. M& O! e+ p9 G$ jError using ==> plot
Vectors must be the same lengths.
8 B+ l8 y0 {. n$ Y3 a3 e
3 M6 h( h) }' |0 N/ {. oError in ==> cor_self at 13
) P; P7 U. q7 z/ r% c2 _plot(k,r,'-*');

zjqylcy

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

梓爱

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

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

zjqylcy

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