灰色轨迹 发表于 2006-12-5 21:14

请教如果画出一个时间序列的自相关函数和偏自相关函数图？？

高手帮忙啊，就像下图中的那样。

lxq 发表于 2006-12-5 22:55

不会~~
不过你不是有图了么？

zhengdan66 发表于 2006-12-6 00:14

他给的图是截屏来的，不是他画的

你应该去看看 自相关函数和偏自相关函数的表达式，然后写程序就行了

happy 发表于 2006-12-6 10:26

xcorr
然后画图就行

灰色轨迹 发表于 2006-12-6 10:35

原帖由 zhengdan66 于 2006-12-6 00:14 发表
他给的图是截屏来的，不是他画的

你应该去看看 自相关函数和偏自相关函数的表达式，然后写程序就行了
我是新手啊，不会画啊，谁能不能给点代码呀:@Q

happy 发表于 2006-12-6 10:39

http://forum.vibunion.com/forum/viewthread.php?tid=25454

[ 本帖最后由 eight 于 2007-4-20 16:38 编辑 ]

灰色轨迹 发表于 2006-12-11 15:12

原帖由 灰色轨迹 于 2006-12-6 10:35 发表

我是新手啊，不会画啊，谁能不能给点代码呀:@Q
看了你给的链接，可是画不了我给的那种图啊:'(

灰色轨迹 发表于 2006-12-14 20:47

ding

xj1983 发表于 2007-4-20 16:15

是不是必须要是零均值序列才画得出，我用零均值就是这个效果，但是如不是零均值则不是这样的图．

[ 本帖最后由 eight 于 2007-4-20 16:37 编辑 ]

xj1983 发表于 2007-4-26 12:14

clear;
%function y=zxg(N)
N=50;
X=;
%Y=X-70.9094
%X=;
                  
%X=;
%X=;
x=0;
for i=1:N
    x=x+X(i);
end
x=x/N;
Y=X-x;
e=0;
b=0;
r=ones(1,N);
p=zeros(1,N);
for c=0:N-1
    sum=0;
    e=e+1;
    for d=e:N
      sum=sum+Y(d)*Y(d-b);
    end
      b=b+1;
      r(c+1)=sum;
      p(c+1)=r(c+1)/r(1);
end
subplot(211);plot(p);
n=N-2;
B=0;
C=0;
A=zeros(n);
%for e=1;N-1
%p(e)=p(e+1);
%end
A(1,1)=p(2);
for k=1:n;
    for i=1:k
      for j=1:k
            B=B+p(k+2-j)*A(k,j);
            C=C+p(j+1)*A(k,j);
      end
      C=1/(1-C);
      A(k+1,k+1)=(p(k+2)-B)*C;
      A(k+1,k+1-i)=A(k,k+1-i)-A(k+1,k+1)*A(k,i);
      B=0;
      C=0;
    end
end
L=ones(1,N);
for i=1:n
    L(i)=A(i,i);
end
subplot(212);plot(L);
其中Ｘ是你的序列，Ｎ是序列长度，这样可以画出一个时间序列的自相关函数和偏自相关函数图！

aski 发表于 2008-4-10 15:27

谢了

200810100659 发表于 2010-11-30 15:18

{:{39}:}

aliu 发表于 2011-5-4 11:50

支持下，，呵呵，也在学习中，，

chayemianbao 发表于 2014-5-12 21:30

xj1983 发表于 2007-4-26 12:14
clear;
%function y=zxg(N)
N=50;


如果N=1000，怎么办？总不能把那1000个数据一个个的敲到里面吧？！

ChaChing 发表于 2014-5-19 15:26

chayemianbao 发表于 2014-5-12 21:30
如果N=1000，怎么办？总不能把那1000个数据一个个的敲到里面吧？！

help load
页: [1]
查看完整版本: 请教如果画出一个时间序列的自相关函数和偏自相关函数图？？