时间序列分析论文

all-z; ………………………………由构造序列Рsumw=zeros(1,6);Рsumw1=0;Рfor j=1:50Р sumw1=sumw1+w(j)^2; ..……………………………..计算РendРfor k=0:5Рfor i=1:(b-k)Р sumw(k+1)=sumw(k+1)+w(i)*w(i+k); …………….......计算РendРendРr=sumw/b;Рr0=sumw1/b;Рp=r/r0; ……………………….计算自相关函数Рkk11=p(2); ………………………计算Рa2=[1,p(2);p(2),1]Рa22=inv(a2);Рkk2=a22*p(1,2:3)'; ………………………计算Рkk22=kk2(2,1);Рa3=[1,p(2),p(3);p(2),1,p(3);p(3),p(2),1]Рa33=inv(a3);Рkk3=a33*p(1,2:4)';Рkk33=kk3(3,1);Рa4=[1,p(2),p(3),p(4);p(2),1,p(2),p(3);p(3),p(2),1,p(2);p(4),p(3),p(2),1];Рa44=inv(a4);Рkk4=a44*p(1,2:5)';Рkk44=kk4(4,1);Рa5=[1,p(2),p(3),p(4),p(5);p(2),1,p(2),p(3),p(4);p(3),p(2),1,p(2),p(3);p(4),p(3),p(2),1,p(2);p(5),p(4),p(3),p(2),1]; Рa55=inv(a5);Рkk5=a55*p(1,2:6)';Рkk55=kk5(5,1); ………………..计算Рkk=zeros(1,5);Рkk=[kk11,kk22,kk33,kk44,kk55];РD=r0-kk11*r(2)-kk22*r(3) ………………..计算