求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
: ]' @% [" j  r- e$ A6 [function r=cor_self(x,p,q)
+ c! O. o. u$ |, Y* b3 C+ \. J%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
clc;close all;
" N$ a2 X: K4 Z9 p- H( L0 Z* v7 k  i, un=size(x,1);
2 j# |  H' T! O: ?; Cfor k=p:q
6 Z" c0 B4 r, M    sum=0;
    for i=1:n-k
        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
5 Z* F: e" h, c6 l7 ^0 v    end
( Z6 ]4 c3 `3 N8 m. v4 o% i    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
end
/ P/ ^; o% o6 _k=[p:q];
plot(k,r,'-*');
refline(0,0);
; ?4 u2 P9 W' G2 W  N; J, V- Rend
运行显示
Error using ==> plot
Vectors must be the same lengths.
; `' a1 |4 e3 T1 ~! D1 Q, d
Error in ==> cor_self at 13
4 r# y  q0 }* yplot(k,r,'-*');

zjqylcy

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

梓爱

zjqylcy 发表于 2012-7-13 16:31
/ z7 V$ d! J3 B+ B' kk r 长度不一样   应该是r多了一位   可以进workspace里面看看

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

zjqylcy

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