高级计算器在灰色预测GM(1,1)模型精度检验中的使用

1、高级计算器在灰色预测GM（1, 1）模型精度检验中的使用摘要：目的 灰色预测GM（1, 1）模型精度检验计算工作繁琐，介绍一种不使用ENT键，进行程序写入和运算的方法。方法 使用CASIO 区X3600 P计算器二利用X、K寄存器输入程序变量，与SD机能组合运算。结果P1区程序写入17步，计算出了 R值。P2区程序写入21步， 计算出了， 3,转换在SD状态，读取n、云、s2, So结论 灰色预测GM（1，1）模型精度检验的程序正好38步完成了计算工作，释放 了计算器的潜在功能，提高了工作效率，计算器在精度检验工作中得到了正确的使用。关键词：计算器；灰色预测；精度检验；程序；使用Applica

2、tion of Advanced Calculator in the Precision Test of Grey Prediction GM(1,1) ModelABSTRACT: Objective The calculation of GM (1,1) model precision test for grey prediction is tedious. A method of program writing and operation without using ENT key is introduced. Methods CASIO 区 -3600P calculator was

3、used, X and K registers were used to input program variables, and combined with SD function. Results The program of P1 area was written in 17 steps and the value of R was calculated. The P2 program was written in 21 steps, calculating 中， k、后 +, transforming in the SD state, reading n、玉、 s2, S. Concl

4、usion The program of accuracy test for grey prediction GM (1K1) model has completed the calculation work in 38 steps, released the latent function of the calculator, improved the working efficiency, and has been used correctly in the precision test work. KEY WORDS: Calculator; Grey prediction; Accur

5、acy test; Program; Use0引言灰色预测GM（1,1）模型应用广泛，为保证建模方法的准确 性、可靠性和可行性，需要对已知数据进行精度检验。CASIOX 3600 P计算器，置有2组最大38步程序存储器和SD机能， 利用X、K寄存器输入变量，编程运算，满足了精度检验的计算 要求。现介绍如下。（限于篇幅公式从略）1写入程序1.1准备MODE O（LRN状态）INV PCL（程序存储器清零）MODE -（RUN状态）INV KAC（K寄存器清零）令：1 Kin1 2 Kin 4Kin 5 （输入数值给需要使用的寄存器）MODE O（LRN状态）1.2程序P1（进入程序P1区）Kin

6、 21Kin + 3Kout 4 Kin + 2+ Kout 1 =一 Kout 2= Kin +5 INV HLT Kout 5 一 Kout 3 =(17步显示2.75 )INV P2（进入程序P2区）Kout 5 Kout 4 Kin 5 INV %INVZ2 INV 丁INV HLTKin + 6Kout 5 INVc2Kin + 1 INV 丁INV HLT Kin + 21 Kin + 3 Kout 6 一 Kout 3 =(21 步显示18.75 )1.3退出与运算MODE -表1灰色预测精度检验等级参照表精度等级相对误差关联度均方差比值小误差概率一级0.010.900.350.

7、95二级0.050.800.500.80三级0.100.700.650.70四级0.200.600.800.602运算实例呼和浩特市19922001年城市人口的统计资料（见表2 ）， 这些数据符合建立灰色预测模型的要求，进行精度检验，对2002年呼和浩特市城市人口做出预测112.1使用P2区程序，计算，：，值。操作显示说明MODE -RUN状态INV KACK寄存器清零MODE 7 2设定小数令：67.10 Kin 4输入x266.13 Kin 5输入X 2INV P21.452表2呼和浩特市城市人口规模（万人）精度检验计算成果表序号12345678910年份19921993199419951

8、99619971998199920002001nXS2xQ66.4067.1068.3069.3070.2072.6075.5078.4079.3081.701072.885.2230 x66.4066.1367.8869.6871.5373.4375.3777.3779.4281.531.450.610.551.891.140.171.310.150.210.970.420.381.330.830.131.030.120.17100.53801.451.030.871.131.130.971.020.910.830.83/100=0.0083M0.5380.4320.1 180.1580.7

9、920.2920.4080.4920.4180.368R1.0000.4070.6130.6360.3330.4450.8360.3920.8470.796100.631RUN0.97ARUN1.45568.30Kin 4输入X367.88Kin 5输入X 3INV P20.615 3RUN0.42A3RUN1.035MODE -RUN状态MODE 7 3设定小数令：0.538 K0.538输入A0.97 =K0.432M20.42 =K-0.118M30.17 =K-0.368M10说明说明RUN状态K寄存器清零设定小数（省略）输入max与p的积R181.70 Kin 4输入Xm输入X 10

10、5 10A。581.53 Kin 5INV P2RUN0RUN00.21L17l832.2转换在SD状态，读取2n， A，s，So操作显示说明MODE 3SD状态MODE 7 4设定小数令：1 Kin + 3n = 10INV X0.5380AINV Xbn0.4412S2INV X crn-10.4650S2.3在SD状态，输入x读取n， xS2, So操作显示说明MODE 3SD状态INV KACK寄存器清零令：66.40 DATA输入X167.10 DATA输入X281.70 DATA输入X10Kout 310nINV x72.8800XINV Xbn5.2230S2INV X crn-

11、15.5055S2.4使用Pl区程序，计算R值。min配）=0,0.97,0.42, 0.38, 1.33, 0.83, 0.13, 1.03, 0.12, 0.17）=0max函）=0,0.97,0.42, 0.38, 1.33, 0.83, 0.13, 1.03, 0.12, 0.17=1.33p =0.5操作显示MODE -INV KACMODE 7 3令：（0 Kin 1 ）1.33 x 0.5 = 0.665 Kin 40 P1 1.0003精度检验及预测3.1相对误差检验本例：后=0.00830.0083 0.6的检验准则。.0.631 0.6精度等级判定：一级3.3后验差检验本例

12、：C = 0.4412 一 5.2230 = 0.0845C = 0.4650 一 5.5055 = 0.0845.0.0845 0.35 精度等级判定：一级3.4小误差概率检验本例：p = P | At - 0.538 | 0.6745 x 5.5055 = 3.7135 = 1.P二0.95精度等级判定：一级根据以上精度检验，表示预测等级良好，模型可用，外推预 测2002年呼和浩特市城市人口为83.69万人。4讨论4.1选择使用文献12程序，输入表2原始数据x。，快速得到下A一行估计值X。，给精度检验提供了基础数据，简单，方便。P1区写入程序时，第1步将变量存入匕寄存器备用，第4 步读取K

13、6寄存器数值后，第5步才与K2寄存器数值相加，再 进行下一步得出结果，平行处理，简化了程序。第2步、第3 步将1累加到K3寄存器为计算均数提供了样本含量。P1区程序运算时，需要提前将残差序列的最小值存入K1 寄存器，最大值与p值的积存入匕寄存器。当p值变动时， 可以在程序运算前调整输入，有利于写入程序和运算。P2区写入程序时，首先将公式（实际值-估计值）/实际 值x 100%修改为（估计值-实际值）/实际值x 100%，因为 公式分子的差是绝对值，所以将被减数和减数位置互换，利用 该类计算器%率键的功能，只用4步就完成了计算，然后进行 平方和开方，显示相对误差的绝对值。P2区写入程序时，根据计算器SD