用EXCEL实现灰色数列模型GM(1，1)的预测.docVIP

用EXCEL实现灰色数列模型GM(1，1)的预测 1 1.1 资料来源 1998年6月份第二次国家卫生服务总调查的浙江省11个县(市区)、90个乡镇(街道)的卫生机构卫生人 力资源的数据资料。 1.2 灰色数列模型GM(1，1)预测方法 设原始数列排成时间数列X(t=0,1,2,…,n)，其中X表示第t时刻的原始数列。 tt 1.2.1 累加生成 通常原始数据呈现离乱现象，灰色理论将无规律的原始数据按(1)式累加生成，使其变为较有规律的生 成数列Y。 t (1) 1.2.2 移动平均数生成 对累加生成数据Y按式(2)作移动平均数生成Z。 tt (2) 1.2.3 建立模型GM(1,1) Y的一阶线性微分方程为： (3) t (3)式即为GM(1，1)模型，其中α、μ为特定系数。 按微分方程的求解方法得到： (4) 其中X为初始时刻的原始数据。根据最小二乘法估计参数得： 0 (5) (6) (7) 其中 1.2.4 求预测值 (8) 因灰色数列模型实际上是生成数列模型，对累加生成数据必须经过逆生成—累减还原后才能使用，即 GM(1,1)模型计算所得结果是预测值的累加和，故预测值由上式计算可得。 1.2.5 计算后验差比值C ? 计算残差δ及δ的样本总体的标准偏差S 1 ? 计算实际数X(1,2,…,n)的样本总体的标准偏差S t2 ? 计算C值 C=S/S, C值是进行预测精确度检验的统计量。 12 1.2.6 外推预测 按式(4)、(8)进行。 1.3 建立电子表格 1.3.1 调整电子表格的结构，如表1所示。 1.3.2 输入原始数据 即1994～1997年的浙江省卫生技术人员的实际数和年份。 1.3.3 编程 2(t=0,1,2,…,n)、移动平均数生成Z以及Z和ZX(t=1,2,…,n) ? 计算累加生成Yttttt 方法：先输入D3=C3,D4=D3+C4,E4=(1)/(2)(D3+D4),F4=POEWR(E4,2),G4=E4×C4，然后同时选中D4、E4、F4、G4，将之内容往下复制两行(即已知年份对应的行)。 ? 计算 方法：输入C7=SUM(C4:C6),E7=SUM(E4:E6),F7=SUM(F4:F6),G7=SUM(G4:G6)。 ? 按公式(7)、(5)、(6)分别计算D、α、μ 方法：输入D10=ROUND((B6×F7-POWER(E7,2)),4),E10=ROUND((C7×E7-B6×G7)/D10,4),F10=ROUND ((F7×C7-E7×G7)/D10,4)。 ? 按公式(4)、(8)计算Y、 tt 方法：先输入H4=ROUND((C3-F10/E10)×EXP(-E10×B4)+F10/E10,0),I4=H4-H3,然后选中H4、I4，再将之内容往下复制两行(注：H3内容可以由H4复制而得)。 ? 计算实际值与预测值的残差δ和残差均数以及残差样本总体标准差S 1 方法：输入J4=C4-I4,J8=ROUND (AVERAGE (J4:J6),4),J10=ROUND (STDEVP (J4:J6),4),并且将J4内容往下复制两个格子。 ? 计算已知年份中非出始年份的实际数的样本总体标准差S和后验差比值C 2 方法：输入C12=ROUND (STDEVP (C4:C6),4),J13=ROUND (J10/C12,2)。 ? 计算|δ(t)-|和0.6745S 2 方法：输入J15=ROUND (ABS(J4-J8),2),J19-ROUND(0.6745×C12,2),再将J15的内容往下复制两格。经过以上几步，就可得到一张GM(1，1)的电子表格。 2 只要输入1994～1997年的实际数，就可以得到目标年份的预测值和预测精确度的检验结果，见表1。 表1 浙江省卫生技术人员的预测 A B C D E F G H I J 2t X年份 Y Z Z ZX y 后验差检验 ttttttt t1994 0 33404 33404 33404 残差(δ) 1995 1 33946 67350 50377 253784212 171009764 67312 33908 38 1996 2 35169 102519 84934.5 7213869290 2987061431 102543 35231 -62 1997 3 35543 139162 120840.5 14602426440 4427958442 139150 36607 36 1998 4 105758 256152 243541378 912511751 177186 38036 1999 5 216708 39522 4 2000 6 D α