The purpose of this paper is to explore modeling mechanism of a nonhomogeneous multivariable grey prediction NMGM(1, m, kα) model and its application. Although multi-variable grey prediction MGM(1, m) model has been employed in many fields, its prediction results are not always satisfactory. Traditional MGM(1, m) model is constructed on the hypothesis that original data sequences are in accord with homogeneous index trend; however, the nonhomogeneous index data sequences are the most common data existing in all systems, and how to handle multivariable nonhomogeneous index data sequences is an urgent problem. This paper proposes a novel nonhomogeneous multivariable grey prediction model termed NMGM(1, m, kα) to deal with those data sequences that are not in accord with homogeneous index trend. Based on grey prediction theory, by least square method and solutions of differential equations, the modeling mechanism and time response function of the proposed model are expounded. A case study demonstrates that the novel model provides preferable prediction performance compared with traditional MGM(1, m) model. This work is an extension of the multivariable grey prediction model and enriches the study of grey prediction theory.
National Natural Science Foundation of China714732331. Introduction
Grey system theory (GST), with the superiority of dealing with uncertain problems that have partially unknown parameters, has been developed greatly since it was applied to system theory [1]. Traditional system analysis theory such as the qualitative theory of dynamic systems [2–18] is built on the assumption of acknowledging the system structure. However, with the rapid development of science and technology, it is impossible to completely master system structure as the system is becoming more complex and uncertain. Grey system theory has been adopted to various aspects of fields including systems analysis, forecasting, and decision-making due to its advantages in tackling semicomplex uncertainty problems.
Forecasting a future development is always of significant importance in energy [19], science and technology [20], and some other fields. A large number of studies on forecasting models and applications have been reported, such as garch-types models needing plenty of data, and sample size would limit the predictive accuracy of those methods [21]. Nevertheless, searching enough effective data is almost impossible either in physical system or in generalized system. Grey prediction models show excellent ability in dealing with small data problems since they were proposed by Professor Deng [22].
Grey forecasting models can be divided into two categories according to the number of variables, single-variable grey forecasting models where GM(1, 1) is the core [23–26] and multivariable grey forecasting models, represented by GM(1, n) [27–29] and MGM(1, m) [30]. Zeng et al. [27] pointed out that there were some structure deficiencies in GM(1, n), which may lead to greater errors. The MGM(1, m) model can reflect mutual relationship among systematic variables and performed better prediction accuracy compared with single variable prediction models. Compared to signal grey prediction models, multivariable grey prediction model MGM(1, m) is a distinct grey model being adequate for considering the mutual interactions of multiple variables in a system, and it is of vital importance in simulating and forecasting the multivariable data sequences.
The nonhomogeneous index data sequences are the most common data existing in all systems, Cui et al. [31] put forward a novel NGM(1, 1, k) model in order to solve the nonhomogeneous index function and established the foundation of our study. This grey model (NGM) is a novel tool to tackle the nonhomogeneous data sequence, which attracted considerable interest of research [32]. How to handle multivariable nonhomogeneous data sequences is an urgent problem; this paper expounds a novel nonhomogeneous multivariable grey prediction model termed NMGM(1, m, kα) to tackle the multivariable nonhomogeneous index data sequences.
In recent years, great attention has been devoted to optimizing and expanding applications of MGM(1, m) model [33–38]. Xiong et al. [33] provided background values to improve the original MGM model, which can be used to eliminate the random fluctuations and errors of the observational data. Guo et al. [37] constructed SMGM through coupling self-memory principle of dynamic system to MGM; the example showed that SMGM had superior predictive performance over other traditional grey prediction models. Karaaslan and O¨zden [38] analyzed Turkey’s ratings credit by multivariable grey prediction model and expanded the application scope of multivariable grey prediction model. Those studies illustrated that multivariable grey prediction model is widely used in many fields.
Through analyzing the existing research on MGM(1, m) model, we find that most of scholars optimized the model from the view of modeling parameters to better fit data sequences with index law rather than optimizing the model from modeling structure of MGM(1, m) model. Traditional MGM(1, m) model was constructed on the hypothesis that data sequences are suitable for grey exponential law resulted from accumulated generation operation. The simulated original data sequences are usually in the form of (1)X^0k=eBkA,k≥2.However, in multivariable system analysis, there are only a few data with characteristic of this hypothesis but more for other hypotheses that original data sequences are in accord with nonhomogeneous index trend. Therefore, it is necessary to propose a novel multivariable grey prediction model that is proper for data sequences that are not in accord with homogeneous index trend. This paper presents the novel model NMGM(1, m, kα) to handle the nonhomogeneous index data sequences in the form of (2)X^0k=eBkA+C,k≥2.By least square method and differential equations, we obtain parameters identification values and time response function of the novel model. The prediction accuracy is theoretically analyzed and a case study is presented to illustrate the effectiveness of the proposed model. The remainder of this paper is organized as follows. A novel nonhomogeneous multivariable grey prediction model and its modeling mechanism are presented in Section 2. The precision analysis and a case study are adopted to demonstrate the effectiveness and practicality of the novel model in Section 3. Our conclusions and future work are given in Section 4.
2. Grey NMGM(1,m,kα) Model
Multivariable grey prediction model is one of the frequently used grey forecasting models. In this section, we present the modeling mechanism and time response function of the novel model NMGM(1, m, kα). The constructing process of NMGM(1, m, kα) is presented below.
Definition 1.
Assume that X(0)=(X1(0),X2(0),…,Xm(0))T is a nonnegative original data matrix and the original nonnegative data vector Xj(0) is (3)Xj0=xj01,xj02,…,xj0n,j=1,2,…,m.The data matrix X(1)=(X1(1),X2(1),…,Xm(1))T is called the first-order accumulated generation vector, where (4)Xj1=xj11,xj12,…,xj1n,j=1,2,…,mand (5)xj1k=∑i=1kxj0i,k=1,2,…,n.The adjacent neighbour average sequence is Z(1)=(Z1(1),Z2(1),…,Zm(1))T and (6)Zj1=zj11,zj12,…,zj1n,j=1,2,…,m,where (7)zj1k=0.5xj1k+xj1k-1,k=2,3,…,n.
The original form of nonhomogeneous multivariable grey prediction model abbreviated NMGM(1, m, kα) is defined as follows:(8)dx11tdt=a11x11t+a12x21t+⋯+a1mxm1t+b1tα,dx21tdt=a21x11t+a22x21t+⋯+a2mxm1t+b2tα,⋮dxm1tdt=am1x11t+am2x21t+⋯+ammxm1t+bmtα,where α≥0. We denote the notation (9)A=a11a12⋯a1ma21a22⋯a2m⋮⋮⋱⋮am1am2⋯amm andB=b1⋮bm.For the convenience of the reader, (8) can be written in matrix, which is (10)dX1tdt=AX1t+Btα.
Definition 2.
Assume that X(0) is nonnegative original data matrix, X(1) is the first-order accumulated generation sequences, and Z(1) is adjacent neighbour average sequences. The differential equation (11)dX1tdt=AZ1t+Btαis said to be basic form of nonhomogeneous multivariable grey prediction model abbreviated as NMGM(1, m, kα). The differential equation (12)dX1tdt=AX1t+Btαis said to be whitening differential equation of grey NMGM(1, m, kα) model. The discrete form of NMGM(1, m, kα) model is (13)xj0k=∑l=1majlzl1k+bjkα,j=1,2,…,m,k=1,2,…,n.
The novel model NMGM(1, m, kα) contains a nonlinear correction term Bkα, and Bkα is said to be time term of the model. The restored values of original data sequences can be adjusted through their coefficients Bkα, which is more suitable for time series prediction. It is easy to see that MGM(1, m) model is a special case of NMGM(1, m, kα) model when α=0.
In the following, we illustrate the modeling mechanism of NMGM(1, m, kα) model.
Theorem 3.
Assume that X(0) is nonnegative original data matrix, X(1) is first-order accumulated generation sequences, and Z(1) is adjacent neighbour average sequences. If aj=(aj1,…,ajm,bj)T, then (14)a^1,a^2,…,a^m=PTP-1PTQ1,Q2,…,Qm,where (15)P=z112z212…zm122αz113z213…zm133α⋮⋮⋱⋮⋮z11nz21n…zm1nnαand (16)Qj=xj02,xj03,…,xj0nT,j=1,2,…,m.
Proof.
Substituting all data values into the discrete form of the model, we obtain (17)x10k=a11z11k+a12z21k+⋯+a1mzm1k+b1kα,x20k=a21z11k+a22z21k+⋯+a2mzm1k+b2kα,⋮xm0k=am1z11k+am2z21k+⋯+ammzm1k+bmkα.For one of the fixed equations j, setting k=2,3,…,n, we have (18)xj02=aj1z112+aj2z212+⋯+ajmzm12+2αbj,xj03=aj1z113+aj2z213+⋯+ajmzm13+3αbj,⋮xj0n=aj1z11n+aj2z21n+⋯+ajmzm1n+nαbj.The matrix form of the above equations is Qj=Paj,j=1,2,…,m.
In order to get parameters vector a^j=(a^j1,…,a^jm,b^j)T(j=1,2,…,m), substituting ∑l=1majlzl(1)(k)+bjkα with ∑l=1ma^jlzl(1)(k)+b^jkα(k=2,3,…,n,j=1,2,…,m), we obtain the error sequence εj=Qj-Pa^j. Let (19)Sj=εjTεj=Qj-Pa^jTQj-Pa^j=∑k=2nxj0k-∑l=1ma^jlzl1k-b^jkα2.Then the parameters vector a^j=(a^j1,…,a^jm,b^j)T making Sj minimum should satisfy (20)∂Sj∂a^j1=∑k=2nb^jkα+∑l=1ma^jlzl1k-xj0kz11k=0,∂Sj∂a^j2=∑k=2nb^jkα+∑l=1ma^jlzl1k-xj0kz21k=0,⋮∂Sj∂b^j=∑k=2nb^jkα+∑l=1ma^jlzl1k-xj0kkα=0.Hence, we deduce that (21)∑k=2nz11k2⋯∑k=2nkαz11k⋮⋱⋮∑k=2nkαz11k⋯∑k=2nk2αa^j1⋮a^jmb^j=z112⋯z11n⋮⋱⋮zm12⋯zm1n2α⋯nαxj02⋮xj0n.Thus, we obtain the desired solution (22)a^j=PTP-1PTQj,where P and Qj are defined in Theorem 3. We can get the results of A^ and B^ by letting j=1,2,…,m.
The NMGM(1, m, kα) model reduces to MGM(1, m) when α=0, and NMGM(1, m, kα) model becomes NMGM(1, m, k) model when α=1. In the following, we give the time response function of MGM(1, m), NMGM(1, m, k), and NMGM(1, m, kα), respectively.
Theorem 4.
Suppose that X(0) is nonnegative original data matrix, X(1) is first-order accumulated generation sequences and Z(1) is adjacent neighbour average sequences. The parameters A and B are obtained by Theorem 3. Then
The time response function of MGM(1, m) model is (23)X^1k=eAk-1X11+A-1B-A-1B,k≥2.
The inverse accumulated generation is (24)X^0k=X^1k-X^1k-1,k≥2.
Theorem 5.
Assume that X(0), X(1), and Z(1) are defined as Theorem 4. The parameters A and B are obtained by Theorem 3. Then
The time response sequence of discrete NMGM(1, m, k) model is (25)X^1k=x11k,x21k,…,xm1kT=eAk-1X01+A-1B+A-12B-A-1Bk+A-12B,k≥2.
The inverse accumulated generation is (26)X^0k=X^1k-X^1k-1,k≥2.
Proof.
(1) Multiplying the whitening differential equation of NMGM(1, m, k) model by e-At we have (27)e-AtdX1tdt-e-AtAX1t=e-AtBt,which is equivalent to (28)de-AtX1tdt=e-AtBt.Integrating the above equation from t0 to t, we get (29)e-AtX1t-e-At0X1t0=-e-AtA-1Bt-e-AtA-12B+e-At0A-1Bt0+e-At0A-12B.Multiplying the equation by eAt, we obtain (30)X1t=eAt-t0X0t0+A-1Bt0+A-12B-A-1Bt+A-12B.Set t0=1 and t=k. We get time response function of discrete NMGM(1, m, k) model (31)X^1k=x11k,x21k,…,xm1kT=eAk-1X01+A-1B+A-12B-A-1Bk+A-12B,k≥2.(2) The restored data is easy to obtain by the definition of first-order accumulated data sequences and so is omitted.
Theorem 6.
Assume that X(0), X(1), and Z(1) are defined as Theorem 4. The parameters A and B are obtained by Theorem 3. The time response function of NMGM(1, m, kα) model is (32)X^1t=eAt∫tαe-AtdtB.
3. Precision Analysis
In this part, we compare the precision of MGM(1, m) and NMGM(1, m, k) model to illustrate the practicality of our results.
3.1. Theoretical Analysis
The predictive performance is of great importance in multivariable grey prediction when constructing a new model. Though many scholars conducted much research to improve the precision of MGM(1, m) model and obtained some progress, their findings are not always satisfactory. Most researchers investigated MGM(1, m) model based on the hypothesis that original data sequences are in accord with homogeneous sequences. The time response function of MGM(1, m) model indicates that the restored data of MGM(1, m) model is (33)X^0k=X^1k-X^1k-1=eAk-2eA-IX11+A-1B,k≥2,where I is unit matrix. Similarly, we obtain the restored values of NMGM (1, m, k) model (34)X^0k=X^1k-X^1k-1=eAk-2eA-IX01+A-1B+A-12B-A-1B,k≥2.Compared with restored data functions of MGM and NMGM model, we find that the restored function of MGM(1, m) model is a data sequence with pure exponential growth law. The simulation and prediction function of NMGM(1, m, k) model is in accordance with nonhomogeneous grey exponential law.
To demonstrate the practicability and maneuverability of nonhomogeneous multivariable grey prediction NMGM model, a case study is employed to compare the predictive performance of MGM and NMGM model.
3.2. Case Study
Per-capita net income of rural households and per-capita disposable income of urban residents are considered to be two important indicators that can reflect people’s living standard and economic level. Therefore, it is necessary and helpful to control the variation trend of per-capita income. However, making such a prediction is challenging because per-capita net income of rural households and per-capita disposable income of urban residents are influenced by many factors. Grey prediction model is fairly appropriate for this problem and shows excellent ability in solving such problems [39, 40]. Zhao et al. [39] forecasted per-capita annual net income of rural households in China by optimized grey GM(1, 1) model, which demonstrated that grey prediction model can be used effectively. In the following, the per-capita net income of rural households and per-capita disposable income of urban residents of Jiangsu province are forecasted by multivariable grey prediction model.
Jiangsu is a representative province of urban development in China, which enters into a new phase of economic development and undergoes a rapid development. Per-capita net income of rural households and per-capita disposable income of urban residents from 2001 to 2013 in Jiangsu province are chosen [41], which are obtained from Jiangsu Province Statistical Yearbook. The average annual growth rate of per-capita net income of rural households is 13.7%, and the average annual growth rate of per-capita disposable income of urban residents is only 11.25%. The per-capita net income of rural households and per-capita disposable income of urban residents are all affected by some certain factors, and there exist interaction and interrelation between them. The data set 2001–2013 is stable, which is suitable for constructing grey prediction model.
The data is divided into two groups: the data set 2001–2007 are used as original data while those 2008–2013 as test. Assume that 2001 is k=1 and so on as shown in Table 1. The input data sequence is the per-capita net income of rural households and per-capita disposable income of urban residents from 2001 to 2007 of Jiangsu province; then parameters of MGM and NMGM model can be determined. On the basis of this, the forecasting results of 2008–2013 can be calculated by Theorems 4 and 5.
The original data of X1 and X2 (unit: yuan).
Time
Variable
k=1
k=2
k=3
k=4
k=5
k=6
k=7
X1(0)
7375
8178
9262
10482
12319
14084
16378
X2(0)
3785
3996
4239
4754
5276
5813
6561
In the following, we construct MGM(1, m) model and NMGM(1, m, k) model to compare the prediction accuracy. Let X1 be the per-capita disposable income of urban residents and X2 be the per-capita net income of rural households.
By traditional MGM(1, m) model, we construct MGM(1, 2) model (35)dx11kdk=0.2732x11k-0.3044x21k+6772.164,dx21kdk=0.08618x11k-0.096x21k+3526.2495.By NMGM(1, m, k) proposed in this paper, we construct NMGM(1, 2, k) model (36)dx11kdk=3.7136x11k-11.2009x21k+15103.6306k,dx21kdk=1.8827x11k-5.7843x21k+7878.2759k.Actual values and prediction values of MGM(1, m) and NMGM(1, m, k) model are presented in Table 2. As can be seen from Table 2, the prediction values of NMGM(1, 2, k) are more close to actual values compared with MGM(1, 2) prediction values.
The actual and prediction values of X1 and X2 (unit: yuan).
Actual values
MGM
NMGM
Time
X1(0)
X2(0)
X^1(0)
X^2(0)
X^1(0)
X^2(0)
k=8
18680
7357
19011.8
7385.7
10547.9
4025.6
k=9
20965
8108
22235.2
8402.1
20863.7
7974.6
k=10
23217
8980
26083.5
9615.7
24427.2
9100
k=11
26341
10805
30677.9
11064.5
28672.5
10441.5
k=12
29677
12202
36163.4
12794.3
33730
12039
k=13
32538
13598
42712.4
14859.5
39755.1
13942.7
To analyze reliability of the model, the accuracy of models should be tested. Two criteria are employed to measure the model, including relative error (RE) and mean relative error (MRE), by which we can investigate the effectiveness of the proposed model, where (37)RE=xik-x^ikxik×100%,MRE=1n∑k=1nxik-x^ikxik×100%.
The relative errors of per-capita disposable income of urban residents and per-capita net income of rural households predicted by MGM(1, 2) and NMGM(1, 2, k) model are listed in Table 3. The mean relative error (MRE) indicates that NMGM(1, 2, k) model is more accurate than traditional MGM(1, 2) model. From the comparative analysis, we know that NMGM(1, m, k) model can better follow the tendency of per-capita disposable income of urban residents and per-capita net income of rural households. The case study shows that the novel nonhomogeneous NMGM(1, m, kα) model has certain advantages compared with traditional multiple grey prediction model.
The relative errors of MGM and NMGM (unit: %).
X1
X2
Time
MGM
NMGM
MGM
NMGM
k=8
1.78
4.32
0.39
4.44
k=9
6.06
0.48
3.62
1.65
k=10
12.34
5.2
7.07
1.34
k=11
16.46
8.8
2.4
3.36
k=12
21.86
12.01
4.85
1.33
k=13
31.27
22.18
9.27
2.53
MRE
14.96
8.83
4.6
2.44
4. Conclusions and Future Research
Forecasting a future development is always an important issue in system analysis; however, because of the limitation of information and knowledge, only part of system structure could be fully realized. Grey forecasting models demonstrated its superiority in dealing with problems that partial information known and partial information unknown.
This paper put forward a novel model named NMGM(1, m, kα) to tackle the nonhomogeneous multivariable data sequences, and the novel model makes up the deficiency of traditional MGM(1, m) model. By least square method and differential equations, we present the modeling mechanism and time response function of the proposed model. The case study shows that the improved model is more accurately than traditional MGM(1, m) model.
The proposed NMGM(1, m, kα) model is fairly appropriate for the systems that are affected by other relative factors and their characteristic values are not in accordance with exponential law completely, and the novel NMGM(1, m, kα) model can improve the adaptability of the traditional grey model.
There are also some problems needed to be solved in our future work, such as the optimization of background value, properties of novel model, finding the best α for some specific cases, and integrating other kinds of optimization techniques with the novel model in order to further improve the prediction accuracy. This work is an extension of the multivariable grey prediction MGM(1, m) model and a great contribution to the development of multivariable grey prediction theory.
Data Availability
The data used to support the findings of this study are from Jiangsu Province Statistical Yearbook and the data are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research is supported by National Natural Science Foundation of China (Grant no. 71473233).
LiuS. F.ForrestJ.YangY. J.Advances in grey systems research2013252118WangP.LiuX.Rapid convergence for telegraph systems with periodic boundary conditions2017201710198256810.1155/2017/1982568MR3729693LiT.RogovchenkoY. V.Oscillation criteria for second-order superlinear Emden-Fowler neutral differential equations2017184348950010.1007/s00605-017-1039-9MR3709308LiT.RogovchenkoY. V.On asymptotic behavior of solutions to higher-order sublinear Emden-Fowler delay differential equations201767535910.1016/j.aml.2016.11.007MR3600247Zbl1362.34106LiT.RogovchenkoY. V.Oscillation criteria for even-order neutral differential equations201661354110.1016/j.aml.2016.04.012MR3518445Zbl1347.34106BohnerM.HassanT. S.LiT.Fite-Hille-Wintner-type oscillation criteria for second-order half-linear dynamic equations with deviating arguments201829254856010.1016/j.indag.2017.10.0062-s2.0-85035246029ShahR.LiT.The thermal and laminar boundary layer flow over prolate and oblate spheroids20181216076192-s2.0-8504166817710.1016/j.ijheatmasstransfer.2017.12.130FengQ.MengF.Explicit solutions for space-time fractional partial differential equations in mathematical physics by a new generalized fractional Jacobi elliptic equation-based sub-equation method201612719745074582-s2.0-8497976920410.1016/j.ijleo.2016.05.147FengQ.MengF.Traveling wave solutions for fractional partial differential equations arising in mathematical physics by an improved fractional Jacobi elliptic equation method201740103676368610.1002/mma.4254MR3664602Zbl1370.350732-s2.0-85007315065ShaoJ.ZhengZ.MengF.Oscillation criteria for fractional differential equations with mixed nonlinearities2013201311910.1186/1687-1847-2013-323MR3213901XuR.Oscillation criteria for nonlinear fractional differential equations20131691797135710.1155/2013/971357MR3142583MengF.HuangY.Interval oscillation criteria for a forced second-order nonlinear differential equations with damping201121851857186110.1016/j.amc.2011.06.069Zbl1235.341042-s2.0-80052264108ChatzarakisG. E.LiT. X.Oscillation criteria for delay and advanced differential equations with nonmonotone arguments2018201818823763410.1155/2018/8237634XuR.MengF.New Kamenev-type oscillation criteria for second order neutral nonlinear differential equations200718821364137010.1016/j.amc.2006.11.004MR2335635Zbl1142.343602-s2.0-34248175708MengF.XuR.Oscillation criteria for certain even order quasi-linear neutral differential equations with deviating arguments2007190145846410.1016/j.amc.2007.01.040MR2335466Zbl1131.343192-s2.0-34249817112LiuL. H.BaiY. Z.New oscillation criteria for second-order nonlinear neutral delay differential equations2009231265766310.1016/j.cam.2009.04.009MR2549731Zbl1175.340872-s2.0-67949115607LiuH.MengF.LiuP.Oscillation and asymptotic analysis on a new generalized Emden-Fowler equation201221952739274810.1016/j.amc.2012.08.106MR2988149Zbl1308.340842-s2.0-84868213011XuR.MengF.Oscillation criteria for second order quasi-linear neutral delay differential equations2007192121622210.1016/j.amc.2007.01.108MR2385585Zbl1193.341372-s2.0-34548444240XuN.DangY.GongY.Novel grey prediction model with nonlinear optimized time response method for forecasting of electricity consumption in China20171184734802-s2.0-8500598097410.1016/j.energy.2016.10.003YinK. D.LiS. Y.LiX. M.Research on China’s strategic emerging marine industry based on a new GRA model2016284114KimH. Y.WonC. H.Forecasting the volatility of stock price index: A hybrid model integrating LSTM with multiple GARCH-type models201810325372-s2.0-8504345158010.1016/j.eswa.2018.03.002DengJ. L.Introduction to grey system theory198911124MR1025897Zbl0701.90057WangC.-N.NguyenH.-K.LiaoR.-Y.Partner selection in supply chain of vietnam’s textile and apparel industry: the application of a hybrid DEA and GM (1,1) approach2017201716782684010.1155/2017/7826840RenF.GuL.Study on transition of primary energy structure and carbon emission reduction targets in China based on Markov chain model and GM (1, 1)20162016849129352-s2.0-85008712481YuanC.LiuS.FangZ.Comparison of China's primary energy consumption forecasting by using ARIMA (the autoregressive integrated moving average) model and GM(1,1) model201610038439010.1016/j.energy.2016.02.0012-s2.0-84958998908YaoT. X.LiuS. F.XieN. M.On the properties of small sample of GM(1, 1) model20093341894190310.1016/j.apm.2008.03.017Zbl1205.60085ZengB.LuoC.LiuS.BaiY.LiC.Development of an optimization method for the GM(1,N) model2016553533622-s2.0-8498994548210.1016/j.engappai.2016.08.007TienT.-L.A research on the grey prediction model GM(1,n)201221894903491610.1016/j.amc.2011.10.055MR2870015WangZ.-X.A GM(1,N)-based economic cybernetics model for the high-tech industries in China20144356726852-s2.0-8490251192610.1108/K-10-2013-0227ZhaiJ.ShengJ. M.FengY. J.The grey model MGM(1, n) and its application1997175109113CuiJ.LiuS.-F.ZengB.XieN.-M.A novel grey forecasting model and its optimization20133764399440610.1016/j.apm.2012.09.052MR3020581Zbl1272.930142-s2.0-84872610875MaX.HuY.-S.LiuZ.-B.A novel kernel regularized nonhomogeneous grey model and its applications201748516210.1016/j.cnsns.2016.12.017MR3607357XiongP. P.DangY. G.WuX. H.LiX. M.Combined model based on optimized multi-variable grey model and multiple linear regression201122461562010.3969/j.issn.1004-4132.2011.04.0102-s2.0-80051915794HanY.XuS.-G.YuC.-W.Multi-variable grey model (MGM (1, n, q)) based on genetic algorithm and its application in urban water consumption simulation20082017453345362-s2.0-52349105168LiS.YangX.LiR.Forecasting China’s coal power installed capacity: a comparison of MGM, ARIMA, GM-ARIMA, and NMGM models201810250610.3390/su10020506WangJ.LiuS.ShaoJ.LongM.WangJ.TangY.Study on dual pre-warning of transmission line icing based on improved residual MGM-Markov theory20181345615692-s2.0-8504123839210.1002/tee.22601GuoX. J.LiuS. F.WuL. F.GaoY. B.YangY. J.A multi-variable grey model with a self-memory component and its application on engineering prediction201542829310.1016/j.engappai.2015.03.0142-s2.0-84928882510KaraaslanA.ÖzdenK. Ö.Forecasting turkey’s credit ratings with multivariate grey model and grey relational analysis201715358361010.1007/s40953-016-0064-1ZhaoZ.WangJ.ZhaoJ.SuZ.Using a Grey model optimized by Differential Evolution algorithm to forecast the per capita annual net income of rural households in China201240552553210.1016/j.omega.2011.10.0032-s2.0-80755125913MengZ. Y.XuC. Y.LanH. X.Fuzzy GM(1, 1) model based per capital income predicted of farmers in the world natural and cultural heritage area: take Leshan city for an exampleProceedings of the Eleventh International Conference on Management Science and Engineering Management201710071018XiongP. P.2012Nanjing, ChinaNanjing University of Aeronautics and Astronautics