基于改进的灰色NGM%281%2C1%29预测模型研究及应用——以新发传染病为例.pdf

Modeling and Simulation 建模与仿真建模与仿真,2023,12(4),3794-3806 Published Online July 2023 in Hans.https:/www.hanspub.org/journal/mos https:/doi.org/10.12677/mos.2023.124347 文章引用文章引用:李胡蓉,陈燕婷.基于改进的灰色 NGM(1,1)预测模型研究及应用J.建模与仿真,2023,12(4):3794-3806.DOI:10.12677/mos.2023.124347 基于改进的灰色基于改进的灰色NGM(1,1)预测模型研究预测模型研究 及应用及应用 以新发传染病为例以新发传染病为例 李胡蓉，陈燕婷李胡蓉，陈燕婷*上海理工大学管理学院，上海 收稿日期：2023年6月1日；录用日期：2023年7月12日；发布日期：2023年7月18日 摘摘 要要 2020年突发传染病严重威胁人们的生命健康，精准的预测对于传染病的管理有重要意义。基于传染病发年突发传染病严重威胁人们的生命健康，精准的预测对于传染病的管理有重要意义。基于传染病发展的不确定性，建立传统灰色展的不确定性，建立传统灰色NGM(1,1)模型，而模型，而NGM(1,1)模型初始值的选取、背景值的构造以及参数模型初始值的选取、背景值的构造以及参数估计方法是导致该模型精度不稳定的重要原因。首先，本文以三参数的紧邻均值背景值来优化传统估计方法是导致该模型精度不稳定的重要原因。首先，本文以三参数的紧邻均值背景值来优化传统NGM(1,1)模型背景值为基础，以平均相对误差最小为目标，利用粒子群优化算法求最优初始值的取值，模型背景值为基础，以平均相对误差最小为目标，利用粒子群优化算法求最优初始值的取值，然后将三参数的紧邻均值背景值改为加权的三参数背景值，利用粒子群优化算法寻求最优的加权权数，然后将三参数的紧邻均值背景值改为加权的三参数背景值，利用粒子群优化算法寻求最优的加权权数，最后将参数估计方法最小二乘法改为加权最小二乘法求模型的灰参数，构成最后将参数估计方法最小二乘法改为加权最小二乘法求模型的灰参数，构成NGM(1,1)模型初始值的选模型初始值的选取、背景值的构造和参数估计方法的组合优化。通过算例来验证优化后的灰色取、背景值的构造和参数估计方法的组合优化。通过算例来验证优化后的灰色NGM(1,1)模型性能，并将模型性能，并将改进的模型应用于传染病的预测。改进的模型应用于传染病的预测。关键词关键词 NGM(1,1)预测模型预测模型，初始值初始值，背景值背景值，参数估计参数估计，传染病传染病 Research and Application of Improved Gray NGM(1,1)Based Prediction Model Taking Emerging Infectious Diseases as an Example Hurong Li,Yanting Chen*Business School,University of Shanghai for Science and Technology,Shanghai Received:Jun.1st,2023;accepted:Jul.12th,2023;published:Jul.18th,2023 *通讯作者。李胡蓉，陈燕婷 DOI:10.12677/mos.2023.124347 3795 建模与仿真 Abstract The outbreak of infectious diseases in 2020 is a serious threat to peoples life and health,and ac-curate prediction is important for the management of infectious diseases.Based on the uncertain-ty of infectious disease development,the traditional gray NGM(1,1)model is established,and the selection of initial values of NGM(1,1)model,the construction of background values and the pa-rameter estimation methods are important reasons for the unstable accuracy of this model.Firstly,this paper optimizes the background value of the traditional NGM(1,1)model based on the three-parameter tight-neighborhood mean background value,and uses the particle swarm opti-mization algorithm to find the optimal initial value with the objective of minimizing the average relative error,then changes the three-parameter tight-neighborhood mean background value to a weighted three-parameter background value,and uses the particle swarm optimization algorithm to seek the optimal weighting weights,as well as finally changes the parameter estimation method least squares Finally,the least squares method of parameter estimation is changed to the weighted least squares method to find the gray parameters of the model,which constitutes the optimization of the combination of the initial value selection,background value construction and parameter es-timation method of the NGM(1,1)model.The performance of the optimized gray NGM(1,1)model is verified by arithmetic examples,and the improved model is applied to the prediction of infec-tious diseases.Keywords NGM(1,1)Prediction Model,Initial,Background Values,Parameter Estimation,Infectious Diseases Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 灰色预测模型是解决“少数据，贫信息”的小数据的重要方法，灰色预测模型以不需要大量的数据样本，短期预测效果好，运算过程简单的特点，且在已知数据少的情况下仍然可以较好的拟合与预测数据的发展趋势，所以灰色预测广泛应用于各个领域1 2 3。朱晓宵等4以武汉 2020 年新冠肺炎情况为例，首先构建滚动灰色 GM(1,1)模型，用来预测疫情的变化趋势，并与传染病动力学模型预测精度比较，得出滚动 GM(1,1)模型 SEIH1H2RD 传染病动力学模型，这说明对于具有不确定性的传染病的传播趋势，灰色预测模型有很好的预测效果。GM(1,1)模型是目前应用最为广泛的灰色预测模型，但是据学者研究，经典 GM(1,1)有时也会出现较大的误差，导致经典的 GM(1,1)模型的模拟精度不理想，而模拟精度主要与模型参数、模型结构和原始数据序列等有关，针对这一问题，目前学者们主要从模型参数优化和模型结构优化这两个方面研究了 GM(1,1)模型的优化问题。如何提高灰色预测模型来预测事物的发展情况也是本文的研究目标。在模型结构优化方面，考虑到经典的 GM(1,1)预测模型参数估计来源于 GM(1,1)模型的差分形式，时间响应函数来源于微分方程，造成二者的“非同源性”，为了实现模型参数估计和时间响应函数来源相同，谢乃明等提出了离散的 GM(1,1)模型，即 DGM(1,1)模型，并且得到了很好的应用5 6。不管是GM(1,1)模型还是 DGM(1,1)模型都只在近似齐次指数数据序列上可以得到较好的模拟效果，对于近似非齐次指数数据的模拟效果不佳，所以有了 GM(1,1)模型和 DGM(1,1)模型的很多拓展。张和平等7将可以Open AccessOpen Access李胡蓉，陈燕婷 DOI:10.12677/mos.2023.124347 3796 建模与仿真 描述趋势为 S 型过程的 Verhulst 模型应用于网络舆情的预测，得出模型的适用性；邱慧8等将 Verhulst模型应用于人口老龄化的预测；石季雨等9等将对于近似非齐次指数数据模拟效果好的 NGM(1,1)模型应用于中国能源消费情况的预测。在模型参数优化方面，主要包括初始值优化、背景值优化和参数估计方法优化三个方面，这也是目前对于灰色预测模型优化最多的三个方面。经典的灰色预测模型是选择一阶累加生成序列的第一个数据作为模型的初始值，而这个选择并不能突出新信息的重要性。所以考虑到新信息的重要性，Zhao 等10引入滚动机制；也有学者引入新陈代谢模型11 12，用每天更新的新数据来替换最初始数据，给予新数据更高的优先级；丁松等13考虑新信息优先原理，引入权重对生成序列加权求和来优化 GM(1,1)幂模型初始值，然后利用平均相对误差最小来确定权重值；郑坚等14同样考虑了生成序列的综合影响，按生成数据的大小对数据赋权值，并且引入了时间权重函数来反应各个误差平方和的重要性。经典的 GM(1,1)模型的背景值是采用紧邻均值构造法，构造出的梯形面积就是经典的灰色模型的背景值，该值与实际的曲边梯形的面积之差就是经典 GM(1,1)模型背景值构造的误差值。为减小经典背景值选取的误差，卢捷等15将背景值设为变量，然后用平均相对误差最小来确定背景值参数，证明了优化后的模型精度提高，同时也对初始值重新取值对模型进行优化；杨孝良等16考虑到原始序列中极端值的影响，提出了三参数紧邻均值来构造背景值，弱化了原始序列中极端值的影响；刘震等17对 NGM(1,1)的微分方程进行积分，通过分析误差产生原因来重新推导背景值表达式。对于灰色预测模型的参数估计方法，经典 GM(1,1)模型的参数估计方式是最小二乘法，最小二乘法要求数据总体为正态分布，而经典 GM(1,1)模型的数据总体特征不明确，且最小二乘法存在稳健性差的问题。许多研究者对此进行了改进，比如何霞18利用加权最小二乘法来提高模型的精度和稳健性；周德强19利用 LS-SVM 来估计灰色 Verhulst 模型的参数，符合灰色 Verhulst 模型小样本的建模，提高了灰色 Verhulst 预测模型的精度。以上研究都是从灰色预测误差主要来源的一个方面或者两个方面对灰色预测模型进行优化，在一定程度上提高了模型的精度和减少了预测模型误差，但是没有总体解决误差来源问题，所以基于以上研究，本研究首先构建 NGM(1,1)模型，从背景值、初始值和参数估计方法三个方面对 NGM(1,1)灰色预测模型进行组合优化。用优化后的灰色 NGM(1,1)预测模型进行算例分析，计算出模型的相对平均误差、方差比(C)和小误差概率(p)三个指标值，通过与其他模型对比三个指标值的大小，比较模型的性能。最后将本文构建的组合优化的模型应用于武汉市新型冠状病毒肺炎新增感染人数模拟，分别采用不同模型建模并与本文模型进行对比，结果表明所提出的优化模型能够获得较高精度。2.NGM(1,1)模型优化模型优化 NGM(1,1)模型是基于经典 GM(1,1)模型构建的一种能够同时模拟近似齐次指数数据和近似非齐次指数数据的灰色预测模型，而 GM(1,1)模型是针对近似齐次指数数据有较为理想的预测效果，由于现实生活中大多数数据都是近似非齐次指数的，近似齐次指数数据只是一种理想状态，所以基于现实数据，NGM(1,1)模型具有更好的预测效果。本节基于文献16依次从初始值取值、背景的构造和参数估计方法三个方面来优化 NGM(1,1)模型。2.1.三参数背景值的三参数背景值的 NGM(1,1)模型模型 NGM(1,1)模型20原始数列为：()()()()(0)(0)(0)(0)1,2,Xxxxn=，累加生成()()()()(1)(1)(1)(1)1,2,Xxxxn=，其中()()(1)(0)1kixkxk=，其中1,2,kn=，用一阶微分方程拟合(1)X，得到 NGM(1,1)的白化微分方程形式：李胡蓉，陈燕婷 DOI:10.12677/mos.2023.124347 3797 建模与仿真 (1)(1)ddxaxbtct+=+(1)生成紧邻均值生成序列()()()()(1)(1)(1)(1)2,3,Zzzzn=，其中()()()()(1)(1)(1)11,2,3,2zkxkxkkn=+=，作为模型的背景值，已知数据是离散的而不是连续的，将()()()()()(1)(1)(1)(1)1(0)(0)(0)11()1dd1kkiixkxkxxxixixktttt=，2,3,kn=，和背景值带入式(1)得到 NGM(1,1)的基本形式：()()(0)(1)xkazkbkc+=+(2)式(2)a、b 和 c 为待估灰参数。根据杨孝良等16为了提高 NGM(1,1)模型背景值的平滑效果，减少模型中极端值对模型精度的影响，提出三参数紧邻均值作为 NGM(1,1)模型的背景值，即：()()()()()(1)(1)(1)(1)1123zkxkxkxk=+(3)其中3,4,kn=，代替传统 NGM(1,1)背景值()()()()(1)(1)(1)112zkxkxk=+，其中2,3,kn=，将三 参数背景值带入式(2)，利用最小二乘法求解参数 a、b 和 c。将()()()(0)(1)(1)1xkxkxk=和式(3)带入式(2)得：()()()(1)(1)(1)1331211113333aabcxkxkxkkaaaa=+(4)令11313aaa=+，2313aaa=+，313baa=+，413caa=+，得到离散的 NGM(1,1)模型：()()()(1)(1)(1)123412xka xka xka ka=+(5)由于式(5)是三阶差分方程求解很困难，所以直接利用递推算法进行求解()(1)xk，然后利用()(0)xk与()(1)xk和()(1)1xk 关系求()(0)xk的预测值()(0)xk。2.2.优化三参数背景值优化三参数背景值 NGM(1,1)模型的初始值模型的初始值 杨孝良等16中是将具有三参数背景值的 NGM(1,1)模型的初始值()(1)1x和()(1)2x看为已知数据，即令()()(1)(1)11xx=，()()(1)(1)22xx=，直接给预测值赋原值，虽然序列前两个值的没有预测误差，但是初始值不是参与模型运算得出可能会导致后面的预测值出现较大的误差，所以这里将初始值()(1)1x和()(1)2x设为未知，令：()(1)11x=(6)()(1)22x=(7)利用式(5)和()(0)xk与()(1)xk和()(1)1xk 关系求预测值()(0)xk，然后以平均相对误差式(8)最小为目标，利用粒子群优化算法确定()(1)1x和()(1)2x的最优值1和2。()()()(0)(0)(0)111nkxkxkfitnessnxk=(8)李胡蓉，陈燕婷 DOI:10.12677/mos.2023.124347 3798 建模与仿真 2.3.优化三参数背景值优化三参数背景值 NGM(1,1)模型的背景值模型的背景值 杨孝良等16提出的三参数背景值与传统的 NGM(1,1)模型的背景值构造一样是取的紧邻生成序列的均值，为了提高模型预测精度，将紧邻均值生成序列的均值改为未知的加权数，即将式(3)改为：()()()()(1)(1)(1)(1)12zkxkxkxk=+(9)并满足1+=。将式(9)和()(0)xk与()(1)xk和()(1)1xk 的关系式代入式(2)得：()()()(1)(1)(1)1121111bcxkxkxkkaaaa=+(10)将寻求的最优初始值1和2带入模型中，再次以平均相对误差最小为目标(式(8)，通过粒子群优化 算法寻求最优的、和值，当13=，背景值变为杨孝良等16提出的紧邻三参数均值背景值 构造。2.4.优化三参数背景值优化三参数背景值 NGM(1,1)模型的参数估计方法模型的参数估计方法 传统的 NGM(1,1)模型及其他的灰色预测模型都是利用最小二乘法求解模型的灰参数，即利用残差平方和最小为目标来求灰参数，即：()()()2(0)(0)3nksxkxk=(11)式(11)没有考虑到新旧信息数据对于模型的贡献不同，对于序列长度为 n 的原始数据，根据数据值发生的时间的不同，给予不同重要程度，当数据值发生时间距离当前时刻越远，则其在建模时所占重要程度就越小，应该占有更小的比重；当数据值发生时间离当前时刻越近的数据，则在建模时的重要程度就越高，应该占有更大的比重。本文参考何霞等18用赋予残差平方和(式(11)权值 w 来反应建模序列值的重要程度，令：()()()()2(0)(0)3nksw kxkxk=(12)其中()()()2(1)1w kzk=。将式(2)改为()()(0)(1)xkazkbkc=+(13)代入式(11)得：()()()()2(0)(1)3nksw kxkazkbkc=+(14)令()()()()()()(0)(0)(0)3344 wxwxyw nxn=，()()()()()()()()()()()()(1)(1)(1)3333344444wzwwwzwwBw nznnw nw n=，T,da b c=使 s 最小的灰参数 a、b 和 c 应该满足：李胡蓉，陈燕婷 DOI:10.12677/mos.2023.124347 3799 建模与仿真 ()()()()()()()()()()()()()()()()()()(0)(1)(1)3(0)(1)3(0)(1)32020210nknknksw kxkazkbkcw kzkasw kxkazkbkckw kbsw kxkazkbkcw kc=+=+=+=(15)根据式(15)解得()()()()()22222T222222222 2,22CECEZF HZEIFC EZECZ FFCEZF GAFCHZCAE Ida b cC EZECZ FFCEAEZEGZCAE HZIFFC EZECZ FF+=+其中：()()()22(1)3NkAwkzk=，()()2(1)3nkZwk zk k=，()()2(1)3nkCwk zk=，()23nkEwk k=，()23nkFwk=，()()()2(1)(0)3nkGwk zk xk=，()()2(0)3nkHwk xk=，()()2(0)3nkIwk xk=()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()T(1)(1)(1)(1)T(1)(1)22(1)2(1)2(1)33333333334444444444 wzwwwzwwwzwwwzwwB Bw n znnw nw nw n znnw nw nwkzkwk zk kwk z =()()()()()()()()()3332(1)22223332(1)22333 nnnkkknnnkkknnnkkkkAZCwk zk kwk kwk kEZEFCEFwk zkwk kwk=()()()()()()()()()()()()()()()()()()()()()()()()()2(1)(0)T(1)(0)3(1)(0)T2(0)3(1)(0)2(0)333333334444444 nknknkwk zk xkwzwwwxwzwwwxB Ywk xk kw n znnw nw nw n xnwk xk=GHI=李胡蓉，陈燕婷 DOI:10.12677/mos.2023.124347 3800 建模与仿真 22T22C EB BZECZ FF=由克拉默法则求得：()221T2220CEZECEZFFCEZFAFCZCAEB BCEAEZEZZCAEFF=则：()()()()()()221TT222222220 CEZECEZFFGCEZFAFCZCAEB BB YHICEAEZEZZCAEFFCECEZF HZEIFCEZF GAFCHZCAE ICEAEZEGZCAE HZIFF=+=+()T,a b c=所以利用加权最小二乘法求解模型参数依然满足：()1TTT,da b cB BB Y=(16)3.NGM(1,1)模型应用模型应用 首先，本节利用的数据实例来验证本文构建的组合优化模型的可行性，然后，通过与其他模型进行对比得出本文构建的组合优化模型的优化性，最后将本文构建的组合优化模型应用于武汉市 2020 年 2 月6 日2 月 15 日的新增感染人数的拟合与预测，来检验本文构建模型的适用性。3.1.算例分析算例分析 通过文献16中的算例分析数据，分别建立 GM(1,1)模型(模型 1)、DGM(1,1)模型(模型 2)、文献14改进初始值模型(模型 3)及本文建立背景值、初始值和参数估计方式组合优化模型(模型 4)，并且计算模型的相对误差()、平均相误差(avg_)、方差比(C)和小误差概率(P)来比较模型的拟合性能和精度。李胡蓉，陈燕婷 DOI:10.12677/mos.2023.124347 3801 建模与仿真 Table 1.Raw data 表表 1.原始数据 年份 2002 2003 2004 2005 2006 2007 数据 21.62 22.83 23.74 23.89 24.50 26.88 年份 2008 2009 2010 2011 2012 数据 27.66 26.31 32.44 37.67 39.50 根据表 1 中的原始数据，计算灰色预测模型模拟值、相对误差和平均相对误差如表 2 所示。其中相对误差()和平均相误差(avg_)的计算公式分别为：()()()()(0)(0)(0)100%xkxkkxk=(17)()1_1nkkavgn=(18)通过寻优得到最优初始值为121.5204=，144.3504=，代入式(10)，通过寻优得到的最优的三参数背景值加权权数为0.2365=，0.7635=，100.1744 10=，最后通过加权最小二乘法(式(16)计算得到灰参数0.3904a=，8.8525b=，30.7150c=，从而通过计算得参数11.4302a=，1227.0521 10a=，39.7532a=，433.8402a=，根据式(10)可知，表 1 数据序列的 NGM(1,1)模型为：()()()(1)(1)12(1)1.430217.5021 1029.753233.8402xkxkxkk=+(19)其中3,4,kn=。同时计算 GM(1,1)模型(模型 1)、DGM(1,1)模型(模型 2)、文献14改进初始值模型(模型 3)及本文建立背景值、初始值和参数估计方式组合优化模型(模型 4)得模拟值、相对误差和平均相对误差来比较模型的模拟性能，如表 2 所示。Table 2.Average error and average relative error of simulated values 表表 2.模拟值的平均误差和平均相对误差 原始数据 模型 1 模型 2 模型 3 模型 4 模拟值 相对误差(%)模拟值 相对误差(%)模拟值 相对误差(%)模拟值 相对误差(%)21.62 21.62 0.00 21.62 0.00 21.62 0.00 21.62 0.00 22.83 20.73 9.21 20.76 9.09 20.84 8.72 22.83 0.00 23.74 22.15 6.70 22.18 6.57 22.27 6.19 23.74 0.02 23.89 23.67 0.91 23.70 0.77 23.80 0.38 23.98 0.37 24.5 25.30 3.26 25.33 3.40 25.44 3.82 24.73 0.95 26.88 27.04 0.58 27.07 0.72 27.18 1.12 25.64 4.60 27.66 28.89 4.46 28.93 4.60 29.05 5.02 27.00 2.38 26.31 30.88 17.36 30.92 17.52 31.04 18.00 28.90 9.86 32.44 33.00 1.72 33.04 1.86 33.18 2.27 31.62 2.53 37.67 35.27 6.38 35.31 6.25 35.46 5.88 35.47 5.83 39.5 37.69 4.59 37.74 4.45 37.89 4.07 40.95 3.67 平均相对 误差(%)5.517 5.525 5.547 3.021 李胡蓉，陈燕婷 DOI:10.12677/mos.2023.124347 3802 建模与仿真 为了进一步比较模型的模拟精度，计算出四个模型的方差比(C)和小误差概率(P)，首先计算出四个模型得出的模拟值与原始序列(表 1)的残差值，残差计算公式：()()()(0)(0)ikxkxk=(20)计算得到残差序列如表 3 所示：Table 3.Residual 表表 3.残差 模型 1 残差 1 模型 2 残差 2 模型 3 残差 3 模型 4 残差 4 0.00 2.10 1.59 0.22 0.80 0.16 1.23 4.57 0.56 2.40 0.00 1.99 1.47 0.09 0.94 0.30 1.39 4.73 0.74 2.21 0.00 1.99 1.47 0.09 0.94 0.30 1.39 4.73 0.75 2.36 0.00 0.00 0.00 0.09 0.23 1.24 0.66 2.59 0.82 2.20 通过表 3 计算得到的残差序列值，进一步计算原始序列的均值(0)x)和标准差(1s)以及三个模型与原始序列的残差序列均值(i)与标准差(1is)、方差比(C)和小误差概率(P)，计算公式如下：()(0)(0)11nkxxkn=(21)()()(0)(0)11nkxkxsn=(22)411iiin=(23)4111iiisn=(24)11iisCs=(25)()()10.6745iiiPpts=(26)式(20)(21)计算表 1 原始数据的均值为 27.91，标准差为 5.76，通过式(22)(26)计算四个模型的平均相对残差(q)、方差比(C)和小误差概率(P)如表 4 所示。小误差概率和方差比对于模型预测精度的评价准则如下：李胡蓉，陈燕婷 DOI:10.12677/mos.2023.124347 3803 建模与仿真 Table 4.Accuracy evaluation criteria 表表 4.精度评价准则 C 和 P 的取值 精度等级 0.35C,0.95P 一级 0.350.50C,0.800.9P 二级 0.500.65C,0.60.8P,0.6P 四级 Table 5.Mean relative error,varianceratio,and small error probability test 表表 5.平均相对残差、方差比和小误差概率检验 指标 模型 1 模型 2 模型 3 模型 4 残差均值 0.10 0.04 0.23 0.20 标准差 1.90(s11)1.91(s12)1.91(s13)1.18(s14)C 0.550 0.329 0.329 0.212 P 0.818 1 0.909 1 精度等级 2 1 1 1 通过表 2、表 5 和表 4 的评价准则比较 4 个模型相关指标得出每个模型的性能，从表 2 可以看出本文所构建初始值、背景值和参数估计方法组合优化的模型(模型 4)的平均相对误差为 3.021%，都小于GM(1,1)模型(模型 1)、将作为初始值的改进 GM(1,1)模型(模型 2)和文献14改进初始值模型(模型 3)，通过表 4 的预测精度评价准则以及表 5 计算得出的指标值比较 4 个模型的方差比(C)和小误差概率(P)，得出本文所构建模型(模型 4)的模拟精度等级、将作为初始值的改进 GM(1,1)模型(模型 2)和文献14改进初始值模型(模型 3)的模拟精度等级最高，且本文所构建模型(模型 4)的模拟精度最高，GM(1,1)模型(模型 1)的模拟精度等级为 2 最差。通过该算例得出本文将初始值通过文献16将初始值设为未知变量1和2，以平均相对误差最小为目标利用粒子群优化算法求最优初始值的取值，将三参数的紧邻均值背景值改为加权的三参数背景值，利用粒子群优化算法寻求最优的加权权数，最后将参数估计方法最小二乘法改为加权最小二乘法求模型的灰参数，构成 NGM(1,1)模型初始值的选取、背景值的构造和参数估计方法的组合优化一定程度上使得模型精度提高，因此本文所构建组合优化模型优于其他三种模型。3.2.模型应用模型应用 下面通过将模型应用于武汉市新型冠状病毒肺炎新增感染人数的预测来分析本文的组合优化模型的效果，将武汉市 2020 年 2 月 6 日2020 年 2 月 15 日(封城阶段，假设不存在人口流动)作为原始数据，比较传统的 GM(1,1)模型、文献16模型、文献14模型和本文组合优化模型平均相对误差来比较模型的预测精度，如表 6 所示。通过寻优得到最优初始值为11686=，23671=，代入式(10)，通过寻优得到的最优的三参数背景值加权权数为0.4910=，146.1000 10=，0.5090=，最后通过加权最小二乘法(式(16)计算得到灰参数0.8661a=，31.3818 10b=，2816.5495c=，从而通过计算得参数11.7399a=，20.7670a=，332.4042 10a=，344.9004 10a=，根据式(10)可知，表 1 数据序列的 NGM(1,1)模型为：()()()(1)(1)(1)331.739910.767022.4042 104.9004 10 xkxkxkk=+(27)其中3,4,kn=，通过式(27)计算出本文组合优化模型的()(1)xk的模拟值，之后根据()(1)xk、()(1)1xk 李胡蓉，陈燕婷 DOI:10.12677/mos.2023.124347 3804 建模与仿真 与()(0)xk的关系，计算出原始数据序列的模拟值，根据模拟值与原始数据本文组合优化模型的相对误差和平均相对误差。同时计算 GM(1,1)模型、文献16模型、文献14模型的模拟值、相对误差和平均相对误差。从表 6 可以看出本文组合优化模型的平均相对误差小于 GM(1,1)模型，且远小于文献16和文献14模型的平均相对误差，且本文基于的文献16的三参数紧邻均值背景值 NGM(1,1)模型平均相对误差最大，这也说明本文在此基础上对初始值、背景值和参数估计的组合优化的有效性。为了更加直观的比较GM(1,1)模型、文献16、文献14和本文组合优化模型对于武汉市 2020 年 2 月 6 日2020 年 2 月 15 日的新增感染肺炎人数的模拟情况，本文绘制 GM(1,1)模型、文献16、文献14和本文组合优化模型的相对误差折线图，如图 1 所示：Table 6.Average error and average relative error 表表 6.平均误差和平均相对误差 GM(1,1)模型 文献16模型 文献14模型 本文组合优化模型 模拟值 相对误差 模拟值 相对误差 模拟值 相对误差 模拟值 相对误差 1501 1501.00 0.00 1501.00 0.00 1501.00 0.00 1686.00 0.12 1985 1935.67 0.02 1985.00 0.00 2526.53 0.27 1985.00 0.00 1379 1832.85 0.33 693.01 0.50 2669.51 0.94 1233.19 0.11 1921 1735.48 0.10 2363.53 0.23 2820.58 0.47 1727.56 0.10 1552 1643.29 0.06 3671.80 1.37 2980.19 0.92 1547.32 0.00 1104 1555.99 0.41 4222.18 2.82 3148.85 1.85 1612.88 0.46 13436 1473.33 0.89 4346.67 0.68 3327.04 0.75 1588.73 0.88 3910 1395.07 0.64 4329.85 0.11 3515.32 0.10 1596.98 0.59 1923 1320.96 0.31 4296.97 1.23 3714.26 0.93 1592.81 0.17 1548 1250.78 0.19 4278.69 1.76 3924.45 1.54 1591.88 0.03 平均相对误差 0.330 0.967 0.863 0.274 Figure 1.Relative error line chart 图图 1.相对误差折线图 12345678910 x轴2020年2月6 日-2020年2月15日00.511.522.53y轴相对误差传统GM(1,1)模型相对误差文献16模型相对误差文献14模型相对误差本文组合优化模型相对误差李胡蓉，陈燕婷 DOI:10.12677/mos.2023.124347 3805 建模与仿真 从图 1 可以更加清楚地看出本文构建的组合优化模型优于 GM(1,1)模型，且远优于文献16和文献14模型。这一结果表明本文构建的组合优化模型具有合理性和有效性。4.结论结论 本文构建对于近似非齐次指数数据模拟效果好的 NGM(1,1)预测模型，考虑到灰色模型的背景值容易受到极端值的影响，所以为了减少原始数据中极端值对于灰色预测模型模拟性能的影响，达到提高灰色预测模型模拟性能稳定性的目的，本文参考文献16构建 NGM(1,1)模型的三参数紧邻均值背景值，基于此，依次对三参数紧邻均值背景值的 NGM(1,1)模型的初始值、背景值和参数估计方式进行优化。首先，以平均相对误差最小为目标利用粒子群优化算法求最优初始值的取值，然后将三参数的紧邻均值背景值改为加权的三参数背景值，利用粒子群优化算法寻求最优的加权权数，最后将参数估计方法最小二乘法改为加权最小二乘法求模型的灰参数，实现 NGM(1,1)模型初始值、背景值和参数估计方式的协同优化。通过算例分析，验证了本文构建的初始值、背景值和参数估计方式组合优化模型模拟性能优于经