求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
2 D# X9 }" C6 afunction r=cor_self(x,p,q)
%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
clc;close all;
n=size(x,1);
for k=p:q
% L  t1 B1 g0 E( N2 \9 R    sum=0;
    for i=1:n-k
7 `+ x; @2 @% z6 u. w0 u+ Q        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
    end
    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
8 v5 m1 ?5 @% X% Y) Oend
k=[p:q];
8 `! b( _$ ]: R7 C$ p7 M- Iplot(k,r,'-*');
refline(0,0);
2 u2 _4 T1 B7 T' p3 t$ T- [end
: j" U! e5 i" T" _' U% T- e. q% b% ^3 Y运行显示
Error using ==> plot
# _/ U. P  i8 S* I& E. CVectors must be the same lengths.
' L  K" s0 {, e6 Z6 N' N* s, o
Error in ==> cor_self at 13
plot(k,r,'-*');

zjqylcy

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

梓爱

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

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

zjqylcy

我看不太明白的是   你前面K是一个变 后面K又成了一个矩阵   
1 P9 r* i6 T8 e  l  naf一般就是分母为0了