求：用MATLAB求自相关系数

梓爱

下面这个程序为什么运行出错呢
1 D$ g8 w9 z/ s0 t$ R! t- ~/ ffunction r=cor_self(x,p,q)
%连续计算一个序列(列向量）的自相关系数，其中滞后时间k由p变化到q
clc;close all;
8 L! l$ \7 D( z9 `n=size(x,1);
for k=p:q
    sum=0;
    for i=1:n-k
7 e" N- _* c7 f4 y. E        sum=sum+(x(i)-mean(x(1:n-k)))*(x(i+k)-mean(x(1+k:n)));
# S- B! ^+ c+ R& O8 p& k    end
    r(k+1)=sum/std(x(1:n-k))/std(x(1+k:n))/(n-k);
2 }' u5 ?8 S- \5 cend
6 Q) d; ?3 J* i' \5 U# Kk=[p:q];
4 J, c: E$ \3 N/ T) Bplot(k,r,'-*');
# E: F$ x6 r* I$ u& B8 @8 nrefline(0,0);
" Q1 l+ V8 S& N: pend
* }! n; @0 x1 D* ^) [' l( m& b; D运行显示
Error using ==> plot
Vectors must be the same lengths.

  F1 B& I+ U" @2 CError in ==> cor_self at 13
/ s% k5 W" M$ \7 qplot(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又成了一个矩阵   
( K* R* R2 f/ S: f  naf一般就是分母为0了