




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、建立GM (1, 1)模型对产品销售额预测祁诗阳冯晓凯申静某大型企业1999年至2004年的产品销售额如下表,试建立GM(1,1)预测模型,并预测 2005年的产品销售额。年份199920002001200220032004销售额 (亿元)2.673.133.253.363.563.72有题目知X(0)= (X(0)(1),. X(0)(6) = (2.67,3.13,3.25,3.36,3.56,3.72)构造累加生成序列X (1) = (X (1),. X (6) = (2.67,5.8,9.05,12.41,15.97,19.69)对X (1作紧邻均值生成_ 1,、z (k) = _ (
2、z (k) + z (k -1)2k = 2,.6编程如下:x=2.67 5.8 9.05 12.41 15.97 19.69;z(1)=x(1);for i=2:6z(i)=0.5*(x(i)+x(i-1);endformat long gz结果如下:z =Columns 1 through 42.67 4.235 7.425 10.73Columns 5 through 614.19 17.83因此z (1) = (z (1)(1),z (5) = (4.235,7.425,10.73,14.19,17.83)于是构造b矩阵和y矩阵如下:-4.235 1、(3.13)-7.425 13.2
3、5B =-10.73 1y =3.36-14.19 13.5617.83 L、3.72 /对参数&进行最小二乘估计,采用matlab编程完成解答如下:B=-4.235 -7.425 -10.73 -14.19 -17.83】,ones(5,1);Y=3.13 3.25 3.36 3.56 3.72;format long ga=inv(B*B)*B*Y结果如下: a =-0.04396098154759662.92561659879905即 a =-0.044, u=2.96- =-66.55d则GM(1,1)白化方程为dxcr0.044x = 2.96 dt预测模型为:X(k +1) = 6
4、9.22e0.044*k - 66.551、关联度检验法:米用matlab编程得到模拟序列for i=1:6X(i)=69.22*exp(0.044*(i-1)-66.55;endformat long gx(1)=X(1);for i=2:6x(i)=X(i)-X(i-1);endX结果如下:x =Columns 1 through 42.673.11367860537808 3.25373920141375 3.40010005288617Columns 5 through 63.553044560121343.71286887145915因此模拟序列为(0)= (0)(1),. (0)(
5、6) = (2.67,3.113,3.253,3.40,3.553,3.712)求模拟序列和原始序列的相关度 (0) = ( (0)(1),. (0)(6) = (2.67,3.13,3.25,3.36,3.56,3.72)初始化,即将该序列所有数据分别除以第一个数据。原始序列变为 x1=(1,1.172,1.217,1.258,1.333,1.393)模拟序列变为 x2=(1,1.166,1.218,1.273,1.331,1.390)序列差=(0,0.006, 0.001,0.002, 0.003)两级差 M=maxmax A =0.006m=minmin A =0计算关联系数取P = 0
6、.50.003门(k)= 一:A(k )| + 0.003n (1) = 1 n (2) = 0.3333) = 0.75 4) = 0.65) = 0.5计算关联度Y = 1 寸 n (k) = 0.6366 0.65k=1因此此模型符合,预测出来的2005年的产品销售额也可信。2、残差检验模拟序列为 (0) = ( (0)(1),. (0)(6) = (2.67,3.113,3.253,3.40,3.553,3.712)原始序列为 (0) = ( (0)(1),. (0)(6) = (2.67,3.13,3.25,3.36,3.56,3.72)由 e(k )= |x(。)(k )一文(。)
7、(k 得,残差序列为:e = (e(1),e(6) = (0,0.017,0.003,0.04,0.007,0.008)所以,相对相对误差为:rel (k) = M)= (0,0.005,0.001,0.012,0.002,0.002) x(0)(k)平均相对误差为:。=1 寸 rel(k) = 0.0037k=1从上述可得到平均精度为99.63%,所以模型符合,预测结果可信。3、后验差检验模拟序列:x (0) = (x (0)(1),. x (0)(6) = (2.67,3.113,3.253,3.40,3.553,3.712) 原始序列:x (0) = (x (0)(1),. x (0)(
8、6) = (2.67,3.13,3.25,3.36,3.56,3.72)残差序列e = (e(1),e(6) = (0,0.017,0.003,0.04,0.007,0.008)x = J_ X(0)(k) = 3.28 k=1采用VC编程完成均方差比值C的解答程序:#include#includevoid main() ( int i;double x6=2.67, 3.13, 3.25, 3.36, 3.56, 3.72;/x 为初始序列double y6=2.67, 3.113, 3.253, 3.4, 3.553, 3.712;/y 为模拟序列double b6;double a=0.
9、00,s,c=0.00,d,e=0.00,f,w;/f 为 S2,s 为 s1,w 为均方差比值 Cfor(i=0;i6;i+)a+=(xi-3.28)*(xi-3.28);s=sqrt(a/6);printf(a=%f,s=%fn”,a,s);for(i=0;i6;i+)bi=xi-yi;printf(b%d=%fn”,i,bi);c+=bi;d=c/6;printf(c=%f,d=%fn”,c,d);for(i=0;i6;i+)e+=(bi-d)*(bi-d);f=sqrt(e/6);w=f/s;printf(f=%f,w=%fn”,f,w);S = :1 寸 X(0)(k)-X=s = 0.8204 1*5k=1S = ;1(k) J = f = 0.03282 6k=1c = % = w = 0.0399(0.5精度为2级,合格1小误差概率:0.5533598)= 1精度为1级好p = p|(k)-e 0.6745S=
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 应试策略:保安证试题及答案
- 达州市大竹县2025届三下数学期末联考模拟试题含解析
- 中国美术学院《Labview》2023-2024学年第二学期期末试卷
- 中国人民警察大学《高等隧道工程》2023-2024学年第一学期期末试卷
- 江西省上饶市“山江湖”协作体统招班2025年高三下学期第三次月考试卷生物试题试卷含解析
- 迎接挑战保安证考试试题及答案
- 大连工业大学艺术与信息工程学院《中华才艺专题一》2023-2024学年第二学期期末试卷
- 西安职业技术学院《硬笔书法》2023-2024学年第二学期期末试卷
- 商洛职业技术学院《景观生态规划》2023-2024学年第二学期期末试卷
- 2025届湖北省七市教科研协作体高三年级第二学期2月周测试英语试题卷含解析
- DLT 560 电力安全工作规程(高压试验室部分)
- 2024年广东广州黄埔区长岭街道森林消防护林员招聘笔试冲刺题(带答案解析)
- ① 小数的意义(课件) 2023-2024学年四年级下册数学人教版
- 巴巴爸爸搬大树
- 2024年云南呈贡区城市投资集团有限公司招聘笔试参考题库含答案解析
- 第2课中华文化的世界意义教学设计-高中历史选择性必修3文化交流与传播
- 儿童康复的健康宣教课件
- 2024年苏州健雄职业技术学院高职单招(英语/数学/语文)笔试历年参考题库含答案解析
- 新生儿围生期感染课件
- 《大数据技术在财务中的应用》 课件 项目1-3 Python概述
- 2024-2025常州新课结束考试化学试卷与答案
评论
0/150
提交评论