对方程(6)可以采用分离变量法,求解得
, (7)
为了利用简单的线性最小二乘法估计这个模型的参数 和 ,将方程(6)表为
, . (8)
利用1978年至2010年的河南省企业职工历年平均工资数据,用MATLAB软件拟合得 , ,将其代入(7)式得到阻滞模型为
. (9)
由上式,可以预测出2013-2035年的河南省职工年平均工资(见表6)
表6 阻滞模型预测河南省年平均工资
年份 2013 2014 2015 2016 2017 2018 2019
年平均工资 45562 51133 57327 64199 71805 80201 89442
年份 2020 2021 2022 2023 2024 2025 2026
年平均工资 99581 110667 122741 135837 149977 165170 181409
年份 2027 2028 2029 2030 2031 2032 2033
年平均工资 198668 216903 236048 256016 276700 297972 319689
年份 2034 2035
年平均工资 341695 363826
该模型考虑了几个阻滞因素,对未来企业职工工资的预测相对较小.
可考虑灰色模型,这是最近几年较为流行的预测模型,它对于影响因素要求不高,不必考虑多种因素对经济发展影响的情况下可以保证较高的求解精度.接下来建立灰色预测模型GM(1,1)对2013—2035年的河南省企业职工年平均工资进行预测,用于对阻滞模型的对比分析.
记 其中 表示第 年河南省职工年平均工资值.
第一步:令 为GM(1,1)建模序列,表示灰导数
第二步:令 为 的AGO序列,对 作累加生成,即得到新的序列 ,
第三步:令 为 的均值序列,表示白化背景值
则得到GM(1,1)的灰微分方程模型为 (11)
式中: 为待估计参数,分别称为发展灰度和内生控制灰度.
其中,
经变换后得到
. (12)
在(12)两端同时乘以 得, ,
即 .
将代入上式中,可得 , 企业退休职工养老金收支问题研究+MATLAB数学建模(6):http://www.youerw.com/shuxue/lunwen_49.html