求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
function r=cor_self(x,p,q)
' L/ ?: t" o* Z%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
clc;close all;
  W' o; U9 w8 ]n=size(x,1);
for k=p:q
    sum=0;
, s$ ]) k, F* `. v( Z    for i=1:n-k
        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
9 k* r3 w% x3 ]1 B8 G& J4 ]    end
2 o. c6 o' |! d! j6 _2 B, b    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
/ r+ E0 W& R6 y1 [! |  p# iend
; P- r0 p+ t% G! F6 X+ S" `k=[p:q];
. I( Z4 y2 q7 h% v; kplot(k,r,'-*');
refline(0,0);
: T* ?2 ?: g1 ^% N5 V9 Z- aend
1 h" f: P+ k5 T# t9 B运行显示
8 k$ G( v3 b- s! s3 B+ b' Q" HError using ==> plot
Vectors must be the same lengths.

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

zjqylcy

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

梓爱

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

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

zjqylcy

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