灰色预测模型论文.docVIP

《统计决策与预测》论文 灰色预测模型 学 院：理学院 专 业：统计学 姓 名: 张 峰 学 号：0914142516 指导教师： 张 军 论文提交日期：二Ο一二年十二月 作业二 成绩: 姓名: 张峰 学号:091412516 问题描述： 设有时间序列数据，如下表所示： 时间 2002 2003 2004 2005 2006 2007 k 1 2 3 4 5 6 15 10 14 7 10 15 试预测2009年的数据。 第一步，构造累加生成列 生成列为： 第二步，构造矩阵B和数据向量 第三步，计算 即： 第四步，得出预测模型 第五步，残差检验 （1）计算 k 0 1 2 3 4 5 15.0000 24.8944 35.3965 46.5437 58.3757 70.9344 （2）累减生成序列 k 1 2 3 4 5 6 15.00000 9.8944 10.5021 11.1472 11.8319 12.5587 （3）计算绝对误差序列和相对误差序列 绝对误差序列为： 相对误差序列为： 第六步，进行关联度检验 min{}={}=0 max{}={}= （1）计算关联系数 由于只有两个序列，故不再寻找第二级最小差及最大差： i 1 2 3 4 5 6 1.0 0.9515 0.3721 0.3333 0.5309 0.4592 （2）计算关联度 满足时的检验准则r0.6。 第七步，后验差检验 计算原始序列的标准差 计算残差的标准差 （3）计算C （4）计算小误差概率 所有的e都小于,故P=1,C=0.06530.35 第八步，模型应用于预测 预测公式为： 当k=7时，可预测出第8期的预测值为： 附录 randint(1,6,[6,15]) ans = 15 10 14 7 10 15 x0=[ 15,10,14,7,10,15] n=length(x0); x1(1)=x0(1); for i=2:n x1(i)=x1(i-1)+x0(i); end x1 B=[(-0.5*[x1(1:5)+x1(2:6)]),ones(5,1)] Y=(x0(2:6)) C=(inv(B*B))*B*Y; a=C(1) b=C(2) for k=1:n x11(k)=[x0(1)-b/a]*exp(-a*(k-1))+b/a; end x11 x00(1)=x11(1); for i=2:n x00(i)=x11(i)-x11(i-1); end x00 wc=abs(x00-x0) format long xdwc=wc./x0 a1=min(wc) a2=max(wc) s=(a1+0.5*a2)./(wc+0.5*a2) r=mean(s) ux0=mean(x0) s1=sqrt([power(x0(1)-ux0,2)+power(x0(2)-ux0,2)+power(x0(3)-ux0,2)+power(x0(4)-ux0,2)+power(x0(5)-ux0,2)+power(x0(6)-ux0,2)]/5) uxdwc=mean(xdwc) s2=sqrt([power(xdwc(1)-uxdwc,2)+power(xdwc(2)-uxdwc,2)+power(xdwc(3)-uxdwc,2)+power(xdwc(4)-uxdwc,2)+power(xdwc(5)-uxdwc,2)+power(xdwc(6)-uxdwc,2)]/5) CC=s2/s1 s0=0.6745*s1 e=abs(xdwc-uxdwc*ones(1,6)) for k=1:n+2 x111(k)=[x0(1)-b/a]*exp(-a*(k-1))+b/a; end x111 x009(1)=x111(1); for i=2:n+2 x009(i)=x111(i)-x111(i-1); end x009