Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Reliability-Aware Cooperative Routing with Adaptive Amplification for Underwater Acoustic Wireless Sensor Networks
Previous Article in Journal
Performance Analysis of Preemptive Priority Retrial Queueing System with Disaster under Working Breakdown Services
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New GM(1,1) Model Based on Cubic Monotonicity-Preserving Interpolation Spline

1
School of Mathematics, South China University of Technology, Guangzhou 510640, China
2
School of Business Administration, South China University of Technology, Guangzhou 510640, China
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(3), 420; https://doi.org/10.3390/sym11030420
Submission received: 28 February 2019 / Revised: 18 March 2019 / Accepted: 19 March 2019 / Published: 21 March 2019

Abstract

:
In the classical GM(1,1) model, an accumulated generating operation is made on the original non-negative sequence to obtain a monotone increasing 1-AGO sequence, and the forecasting model is established based on the 1-AGO sequence. A great number of scholars have improved the accuracy of grey model prediction through better developed background value and the equation for the time response. In this work, we reconstruct the background value based on a new developed monotonicity-preserving piecewise cubic interpolations spline, and thereby establish a new GM(1,1) model. Numerical examples show that the new GM(1,1) model has better prediction quality of data than the original GM(1,1) model and improves the precision of prediction in practice.

1. Introduction

Let an original non-negative and uniformly spaced sequence be
X ( 0 ) = x ( 0 ) ( 1 ) , x ( 0 ) ( 2 ) , , x ( 0 ) ( n ) .
The main idea of the classical grey forecasting GM(1,1) model proposed by Deng [1,2] is to make an accumulated generating operation on the original sequence, so as to reduce the randomization of the original data and obtain an obviously monotone increasing 1-AGO sequence X ( 1 ) . Then, establish a first-order grad forecasting differential equation on the sequence X ( 1 ) . In addition, further use the least-square method to numerically solve the differential equation to estimate the parameters. Finally, the original data simulation and prediction are carried out by using the inverse accumulated generating operation.
The 1-AGO sequence X ( 1 ) is given as follows
X ( 1 ) = x ( 1 ) ( 1 ) , x ( 1 ) ( 2 ) , , x ( 1 ) ( n ) ,
where
x ( 1 ) ( k ) = i = 0 k x ( 0 ) ( k ) = x ( 1 ) ( k 1 ) + x ( 0 ) ( k ) , k = 1 , 2 , , n .
From Equation (3), we can see that the 1-AGO sequence X ( 1 ) has the property of monotonicity-increasing. Suppose that x ( 1 ) ( t ) meets the following first-order grad forecasting differential equation
d x ( 1 ) ( t ) d t + a x ( 1 ) ( t ) = b ,
where the grey developmental coefficient a and the grey control parameter b are the parameters in the model to be estimated.
The solution of the Equation (4) with the initial condition X ˜ ( 1 ) ( 1 ) = X ( 1 ) ( 1 ) is as follows
x ( 1 ) ( t ) = x ( 1 ) ( 1 ) b a e a ( t 1 ) + b a .
Therefore, to obtain the prediction model of the raw data sequence, we need to identify the effect of the grey development coefficient a and the grey control parameter b in Equation (4). For this purpose, we do the integral accumulation on both sides of Equation (4) for k , k + 1 , k = 1 , 2 , , n 1 , then we can get
k k + 1 d x ( 1 ) ( t ) d t d t + a k k + 1 x ( 1 ) ( t ) d t = b ,
that is
x ( 1 ) ( k + 1 ) x ( 1 ) ( k ) + a k k + 1 x ( 1 ) ( t ) d t = b ,
or
x ( 0 ) ( k + 1 ) + a k k + 1 x ( 1 ) ( t ) d t = b .
Let background value be z ( 1 ) ( k + 1 ) : = k k + 1 x ( 1 ) ( t ) d t . To calculate the background value z ( 1 ) ( k + 1 ) , we need to integrate x ( 1 ) ( t ) , which requires the values of a and b to be given in advance. However, the values of a and b need to be determined from the Equation (6). Consequently, to estimate the values of a and b, we must use some methods to estimate the background value z ( 1 ) ( k + 1 ) . We use the piecewise linear polynomial interpolation L ( t ) : = k + 1 t x ( 1 ) ( k ) + t k x ( 1 ) ( k + 1 ) to approximate x ( 1 ) ( t ) in the classical GM(1,1) model, see [1,2], then we get the estimated background value z ( 1 ) ( k + 1 ) as follows
z ( 1 ) ( k + 1 ) = k k + 1 x ( 1 ) ( t ) d t k k + 1 L ( t ) d t = 1 2 x ( 1 ) ( k ) + x ( 1 ) ( k + 1 ) .
For each interval [ k , k + 1 ] , k = 1 , 2 , , n 1 , by substituting the estimated background value z ( 1 ) ( k + 1 ) into Equation (6) and further applying the least-square method, we estimate the values of the parameters a and b by the formula as follows
a b = G T G 1 G T X ,
where
X = x 0 ( 2 ) x 0 ( 3 ) x 0 ( n ) , G = z ( 1 ) ( 2 ) 1 z ( 1 ) ( 3 ) 1 z ( 1 ) ( n ) 1 .
Finally, we get the following estimated solution to the differential Equation (4) with the initial condition X ˜ ( 1 ) ( 1 ) = X ( 1 ) ( 1 ) as follows
x ˜ ( 1 ) ( t ) = x ( 1 ) ( 1 ) b a e a t 1 + b a , k = 1 , 2 , .
We thus get the following grey prediction equation
x ˜ ( 0 ) ( k + 1 ) = x ˜ ( 1 ) ( k + 1 ) x ˜ ( 1 ) ( k ) , k = 1 , 2 , .
From (7), it can be found that the classical GM(1,1) model uses average of adjacent values to estimate the background value z ( 1 ) ( k + 1 ) . Its geometric meaning is to use the area of straight ladder to replace the trapezoidal area which is based on the edge of exponential curve x ( 1 ) ( t ) , as shown in Figure 1. This method has a shortcoming that when the 1-AGO data sequence varies greatly, the result of prediction will have large error ( Δ S ) with the exponential increasing. Thus, if we apply the classical GM(1,1) model in practice, there exist some restrictions. As pointed out in [3,4,5], the accuracy of prediction in GM(1,1) model depends on the estimation of the background value z ( 1 ) ( k + 1 ) . In [6], Li and Dai reconstructed x ( 1 ) ( t ) by a high-order Newton interpolation polynomial. In addition, they estimated the background value z ( 1 ) ( k + 1 ) based on the Newton-Cores integral. However, as shown in [6], when there is a mass of data, the high-order Newton interpolation polynomial may have the Runge phenomenon and thus the truncation error may be very large. Furthermore, the numerical stability is not guaranteed when calculating the approximate value of Newton-Cores integral. In [7], Tang and Xiang estimated the background value z ( 1 ) ( k + 1 ) by using the piecewise quadratic interpolation polynomial to reconstruct x ( 1 ) ( t ) on the interval [ k , k + 1 ] . The advantage is that it has the characteristics of less computation and good numerical stability. In [8], Wang et al. used piecewise cubic interpolation spline to reconstruct x ( 1 ) ( t ) and thus got the estimated background value z ( 1 ) ( k + 1 ) . The advantages of the given method are that it can avoid the Runge phenomenon resulted from high-order polynomial. However, the methods we mentioned above all ignore the important monotonicity-increasing feature of the curve x ( 1 ) ( t ) to be reconstructed. If the resulting reconstructed curve x ( 1 ) ( t ) loses the monotonicity-increasing feature, there will also result in a large error on the background value z ( 1 ) ( k + 1 ) . Therefore, an accurate approximation of the reconstructed curve x ( 1 ) ( t ) is the key to enhance the estimation of the background value.
Recently, lots of scholars proposed many methods to improve the classical GM(1,1) model. Liu and his colleagues have proposed some new methods of grey development clustering based on the multiattribute difference, which enriches the grey fixed-weight cluster theory, see [9,10,11] for example. In [12,13,14], Xie and his colleagues proposed some new grey models, including Verhulst model and the segmental corrected new information discrete GM(1,1) model. In [15], Yang and his colleagues presented a dynamic GM(1,1) model based on grey system theory and cubic spline function interpolation principle. Many scholars have put many efforts on the establishments of new grey prediction models, see for example [16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] and the references therein.
In this paper, we shall propose a monotonicity-preserving piecewise cubic interpolation spline to reconstruct the curve x ( 1 ) ( t ) and thereby give a new scheme to estimate the background value z ( 1 ) ( k + 1 ) . The rest of this paper is structured as follows. Section 2 represents the construction of the C 1 monotonicity-preserving cubic interpolation splines. In Section 3, based on C 1 monotonicity-preserving cubic interpolation spline, a new GM(1,1) model is constructed in detail. Several numerical examples are also given. In addition, the conclusion is given in Section 4.

2. C 1 Monotonicity-Preserving Piecewise Cubic Interpolation Spline

According to Equation (3), the 1-AGO sequence X ( 1 ) has the property of monotonicity-increasing, that is x ( 1 ) ( k ) x ( 1 ) ( k + 1 ) , k . The fitting exponential curve x ( 1 ) ( t ) to the 1-AGO sequence X ( 1 ) is also monotonicity-increasing and has infinite smoothness. Therefore, we develop a C 1 monotonic-preserving cubic interpolation spline to interpolate the 1-AGO sequence, to reconstruct the curve x ( 1 ) ( t ) .
For the discrete data k , x ( 1 ) ( k ) , k = 1 , 2 , , n , we denote d ( 1 ) k as the derivative value at node t = k . For t k , k + 1 , a cubic interpolation spline with local parameter α k is constructed as follows
B ( t ) = 1 s 3 x ( 1 ) ( k ) + 3 1 s 2 s x ( 1 ) ( k ) + d ( 1 ) k α k + 3 1 s s 2 x ( 1 ) ( k + 1 ) d ( 1 ) k + 1 α k + 1 + s 3 x ( 1 ) ( k + 1 ) ,
where s = t k 0 , 1 . From Equation (10), after some computations, we have
B ( k ) = x ( 1 ) k , B ( k + 1 ) = x ( 1 ) k + 1 , B ( k ) = 3 d ( 1 ) k α k , B ( k + 1 ) = 3 d ( 1 ) k + 1 α k + 1 ,
which indicates that B k = B k + , B k = B k + . This means that the cubic interpolation spline defined by Equation (10) is C 1 continuous for arbitrary nonzero local parameter. Here, C 1 continuity means that a function together with its first-order derivative function is continuous. In addition, it is of interest to note that for all α k = 3 , the cubic interpolation spline given in Equation (10) returns into the classic cubic Hermite interpolation spline.
In practical application, we should estimate the derivative values of the cubic interpolation spline at the nodes at first. In this paper, we calculate the derivative values by the following method
d ( 1 ) 1 = x ( 1 ) 2 x ( 1 ) 1 , d ( 1 ) k = 1 2 x ( 1 ) k + 1 x ( 1 ) k 1 , k = 2 , 3 , , n 1 , d ( 1 ) n = x ( 1 ) n x ( 1 ) n 1 .
We shall derive sufficient conditions for the C 1 cubic interpolation spline preserving monotonicity. For the monotonicity-increasing 1-AGO sequence X ( 1 ) , it is obvious that the derivative value determined by Equation (12) is non-negative, which means d ( 1 ) k 0 , k . Without loss of generality, for t k , k + 1 , direct computation gives that
B ( t ) = 3 1 s 2 d ( 1 ) k α k + s 2 d ( 1 ) k + 1 α k + 1 + 6 1 s s x ( 1 ) k + 1 x ( 1 ) k d ( 1 ) k α k + d ( 1 ) k + 1 α k + 1 .
Thus, we can see that the following conditions are sufficient to guarantee R ( t ) 0
α k > 0 , α k + 1 > 0 , x ( 1 ) k + 1 x ( 1 ) k d ( 1 ) k α k + d ( 1 ) k + 1 α k + 1 = x ( 1 ) k + 1 x ( 1 ) k 2 d ( 1 ) k α k + x ( 1 ) k + 1 x ( 1 ) k 2 d ( 1 ) k + 1 α k + 1 0 .
Summarizing the above discussion, we can obtain the following sufficient conditions for B ( t ) preserving monotonicity
α 1 = max 0 , 2 d ( 1 ) 1 x ( 1 ) 2 x ( 1 ) 1 + a 1 , α k = max 0 , 2 d ( 1 ) k x ( 1 ) k + 1 x ( 1 ) k , 2 d ( 1 ) k x ( 1 ) k x ( 1 ) k 1 + a k , k = 2 , , n 1 , α n = max 0 , 2 d ( 1 ) n x ( 1 ) n x ( 1 ) n 1 + a n ,
where a i > 0 serves as free parameter. In practice, we recommend all a k = 3 . If there is no special explanation below, we will take all a k = 3 .

3. Establish New GM(1,1) Model

For the original non-negative sequence X ( 0 ) = x ( 0 ) ( 1 ) , x ( 0 ) ( 2 ) , , x ( 0 ) ( n ) , we first calculate its 1-AGO sequence X ( 1 ) = x ( 1 ) ( 1 ) , x ( 1 ) ( 2 ) , , x ( 1 ) ( n ) . Then for the 1-AGO sequence X ( 1 ) , we use the sufficient conditions (14) with all the free parameters a k = 3 to determine a C 1 monotonicity-preserving piecewise cubic interpolation spline B ( t ) to reconstruct the exponential curve x ( 1 ) ( t ) . For the interval k , k + 1 , we estimate the background value z ( 1 ) k + 1 = k k + 1 x ( 1 ) t d t by the following method
z ( 1 ) k + 1 = k k + 1 x ( 1 ) t d t k k + 1 B t d t = 1 4 x ( 1 ) ( k ) + 1 4 x ( 1 ) ( k ) + d ( 1 ) k α k + 1 4 x ( 1 ) ( k + 1 ) d ( 1 ) k + 1 α k + 1 + 1 4 x ( 1 ) ( k + 1 ) .
Then by substituting the estimated background value into the grey differential Equation (6), we further apply the following least-square method to solve Equation (6)
a b = G T G 1 G T X ,
where
X = x 0 ( 2 ) x 0 ( 3 ) x 0 ( n ) , G = z ( 1 ) ( 2 ) 1 z ( 1 ) ( 3 ) 1 z ( 1 ) ( n ) 1 .
Finally, we get the following estimated solution to the differential Equation (4) with the initial condition X ˜ ( 1 ) ( 1 ) = X ( 1 ) ( 1 ) as follows
x ˜ ( 1 ) ( t ) = x ( 1 ) ( 1 ) b a e a t 1 + b a .
We thus get the following grey prediction equation
x ˜ ( 0 ) ( k + 1 ) = x ˜ ( 1 ) ( k + 1 ) x ˜ ( 1 ) ( k ) = 1 e a x ( 1 ) ( 1 ) b a e a k , k = 1 , 2 , .
We shall give several examples to show that the new GM(1,1) model based on C 1 monotonicity-preserving piecewise cubic interpolation spline has better predict accuracy than the classical GM(1,1) model. In the following examples, the relative error is computed by
ε = x ¯ ( 0 ) ( k ) x ( 0 ) ( k ) x ( 0 ) ( k ) .
Example 1.
In this example, we use the exponential function f ( t ) = η e λ t in [20] to generate the original data with convexity. By letting the parameters g and k take a fixed value respectively, the original data sequence can be obtained when the parameter t takes different values. In addition, we compare the results predicted by our new GM(1,1) model with the GM(1,1) model and the method proposed in [20]. The results show that our new model achieves the best results among the three prediction models and it performs very well in predicting data with the exponential growth trend, see Table 1. On the left of Figure 2, the figure shows the 1-AGO data of Table 1 and the curves of piecewise linear interpolant, monotonic-preserving cubic interpolation spline B ( t ) and the reconstruct exponential curve X ( 1 ) ( t ) . From the right Figure 2, we can see that the interpolation spline B ( t ) has C 1 continuity.
Example 2.
In this example, we directly use the numerical example provided in [21]. We compare our new GM(1,1) model with the classical GM(1,1) model and the prediction models presented in [21]. The results turn out that our new GM(1,1) model prediction accuracy is significantly higher than the classical GM(1,1)model. While the method proposed in [21] performs the best in this example. The reason we suppose is that the method given in [21] simulates and predicts the sequence of exponential distribution by optimizing the background value of the grey differential equations, so that the model has better simulation and prediction accuracy. Table 2 and Figure 3 give the numerical results.
Example 3.
In this example, we use the example in [8] to test the new proposed GM(1,1) model. To verify the applicability of the model to predict electricity consumption in the smart grid, we use the electricity consumption data of Jiangsu province in 2008 into the smart grid as the raw data (in KWh). The results show that the new GM(1,1) has improved prediction accuracy compared to the methods proposed in [7] and the classical GM(1,1) model. Table 3 and Figure 4 give the numerical results.
Example 4.
In this example, we select the grey prediction data from [9]. This set of data is based on Chinese health statistics of people with syphilis (in millions) from 2000 to 2010. We compare the new GM(1,1) model with the classical GM(1,1) model and the model in [9]. From Table 4, we can see that our new GM(1,1) model prediction accuracy is improved compared to the classical GM(1,1) model and the model in [9]. Figure 5 gives the graphic results of this example.
Example 5.
In [11], the data selected by the authors to show a downward trend is very representative, so we use this data sequence in this example. The results show that our new GM(1,1) model is still applicable to the descending sequence, and the accuracy is still slightly improved compared to the classical GM(1,1) model. While the method proposed in [11] performs the best in this example. The reason we suppose is that the method given in [11] simulates and predicts the sequence of exponential distribution by optimizing the background value of the grey differential equations, so that the model has better simulation and prediction accuracy. Table 5 and Figure 6 give the numerical results.
Example 6.
In this example, we selected the relevant data from the 2012 China Statistical Yearbook in [24] to analysis of whether our GM(1,1) model can effectively predict general statistics. Thus, we only compare the prediction results of the classical GM(1,1) model with our new GM(1,1) model. The results show that our new GM(1,1) model applies to general sociological statistics, and our model performs better accuracy of prediction compared to classical GM(1,1) model. This shows that our new GM(1,1) model has the ability for further promotion and development. Table 6 and Figure 7 give the numerical results.
According to the results of the above numerical examples 1–6, the prediction accuracy of our new GM(1,1) model is improved for all the numerical examples compared to classical GM(1,1) model. In addition, our new model performs better than the methods proposed in [7,9,20]. There are different degrees of improvement for different data features. Based on the above data characteristics, we make the conclusion that the data applicable to our new GM(1, 1) model are generally with the continuously increasing feature over time. In particular, the exponential growth data can show better prediction.

4. Conclusions

By using a new developed C 1 monotonicity-preserving piecewise cubic interpolation spline to reconstruct the background value, we have established a new GM(1,1) model. Numerical examples show that the new GM(1,1) model can improve the forecasting quality, especially in prediction reliability and this model performs better when the original data are presented with convexity in time series. Future work will concentrate on exploring more applications of the new GM(1,1) model, such as scientific decision-making in electricity production and manufactures.

Author Contributions

Conceptualization, Y.Z. and Z.J.; Methodology, Y.Z. and Z.J.; Software, Y.Z. and Z.J.; Validation, Y.Z.; Formal Analysis, Y.Z.; Investigation, W.C.; Resources, Y.Z. and W.C.; Data Curation, Y.D. and J.F.; Writing—Original Draft Preparation, Y.Z. and Z.J.; Writing—Review & Editing, Y.Z., Z.J. and Y.D.; Visualization, W.C. and J.F.; Supervision, Y.Z., Z.J. and Y.D.

Funding

This research was funded by the National Natural Science Foundation of China (No. 61802129), the Postdoctoral Science Foundation of China (No. 2015M571931), the Fundamental Research Funds for the Central Universities (No. 2017MS121), the Natural Science Foundation Guangdong Province, China (No. 2018A030310381) and the National Training Program of Innovation and Entrepreneurship for Undergraduates 2018.

Acknowledgments

We thank to the anonymous reviewers for their insightful suggestions and recommendations, which led to the improvements of presentation and content of the paper.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Deng, J. Control problems of Grey system. Syst. Control Lett. 1982, 5, 288–294. [Google Scholar]
  2. Deng, J. Grey system basic method. In The Basis of Grey Theory; Julong, D., Ed.; Press of Huazhong University of Science and Technology: Wuhan, China, 2002; pp. 210–313. [Google Scholar]
  3. Tan, G. The structure method and application of background value in grey system GM(1,1) model (I). Syst. Eng. Theory Pract. 2010, 4, 98–103. [Google Scholar]
  4. Tan, G. The structure method and application of background value in grey system GM(1,1) model (II). Syst. Eng. Theory Pract. 2010, 5, 125–132. [Google Scholar]
  5. Tan, G. The structure method and application of background value in grey system GM(1,1) model (III). Syst. Eng. Theory Pract. 2010, 6, 70–74. [Google Scholar]
  6. Li, J.; Dai, W. A new approach of background value building and its application based on data interpolation and Newton-Cotes formula. Syst. Eng. Theory Pract. 2014, 4, 22–126. [Google Scholar]
  7. Tang, W.; Xiang, C. The improvements of forecasting method in GM(1,1) model based on quadratic interpolation. Chin. J. Manag. Sci. 2006, 14, 109–112. [Google Scholar]
  8. Wang, X.J.; Yang, S.L.; Ding, J.; Wang, H.J. Dynamic GM(1,1) model based on cubic spline for electricity consumption prediction in smart grid. Chin. Commun. 2010, 7, 83–88. [Google Scholar]
  9. Wu, L.; Liu, S.; Yao, L.; Yan, S. The effect of sample size on the grey system model. Appl. Math. Model. 2013, 37, 6577–6583. [Google Scholar] [CrossRef]
  10. Lin, Y.; Ting, Z.; Bingting, Q.; Liu, S. Object multi-attribute differences based grey dynamic clustering method and its application. Oper. Res. Manag. Sci. 2018, 27, 57–63. [Google Scholar]
  11. Wu, L.; Liu, S.; Fang, Z.; Xu, H. Properties of the GM(1,1) with fractional order accumulation. Appl. Math. Comput. 2015, 252, 287–293. [Google Scholar] [CrossRef]
  12. Cui, J.; Liu, S.; Zeng, B.; Xie, N. Parameters characteristics of grey Verhulst prediction model under multiple transformation. Control Decis. 2013, 28, 605–608. [Google Scholar]
  13. Yao, T.X.; Liu, S.F.; Xie, N.M. Study on the properties of new information discrete GM(1,1) model. J. Syst. Eng. 2010, 25, 164–170. [Google Scholar]
  14. Wang, B.; Xie, N. Unified representation and properties of generalized grey relational analysis models. Syst. Eng. Theory Pract. 2019, 39, 226–235. [Google Scholar]
  15. Hou, L.; Yang, S.; Wang, X. Mid-term load forecasting based on buffer operator and modified grey model. J. Syst. Simul. 2013, 25, 1–5. [Google Scholar]
  16. Fan, J.; Gijbels, I. Local Polynomial Modeling and Its Applications; Chapman and Hall: London, UK, 1996. [Google Scholar]
  17. Bohlin, T. Practical Grey-Box Process Identification: Theory Applications; Springer: Berlin, Germany, 2006. [Google Scholar]
  18. Wang, Z.; Dang, Y.; Liu, S. An optimal GM(1,1) based on the discrete function with exponential law. Syst. Eng. Theory Pract. 2008, 2, 61–67. [Google Scholar] [CrossRef]
  19. Tien, T. A new grey prediction model FGM(1,1). Math. Comput. Model. 2009, 49, 1416–1426. [Google Scholar] [CrossRef]
  20. Wang, Y.; Dang, Y.; Li, Y.; Liu, S. An approach to increase prediction precision of GM(1,1) model based on optimization of the initial condition. Expert Syst. Appl. 2010, 37, 5640–5644. [Google Scholar] [CrossRef]
  21. Wang, Y.; Liu, Q.; Tang, J.; Cao, W.; Li, X. Optimization approach of background value and initial item for improving prediction precision of GM(1,1) model. J. Syst. Eng. Electron. 2014, 25, 77–82. [Google Scholar] [CrossRef]
  22. Liu, J.; Xiao, X.; Guo, J.; Mao, S. Error and its upper bound estimation between the solutions of GM(1,1) grey forecasting models. Appl. Math. Comput. 2014, 246, 648–660. [Google Scholar] [CrossRef]
  23. Chen, P.; Yu, H. Foundation settlement prediction based on a novel NGM model. Math. Probl. Eng. 2014, 2, 1–8. [Google Scholar] [CrossRef]
  24. He, Z.; Shen, Y.; Li, J.; Wang, Y. Regularized multivariable grey model for stable grey coefficients estimation. Expert Syst. Appl. 2015, 42, 1806–1815. [Google Scholar] [CrossRef]
  25. Sidorov, D. Integral Dynamical Models: Singularities, Signals and Control; Vol. 87 of World Scientific Series on Nonlinear Science Series A; World Scientific Publ. Pte Ltd.: Singapore, 2015. [Google Scholar]
  26. Duan, H.; Xiao, X.; Pei, L. Forecasting the short-term traffic flow in the intelligent transportation system based on an inertia nonhomogenous discrete gray model. Complexity 2017, 2017, 1–16. [Google Scholar] [CrossRef]
  27. Ma, X.; Liu, Z. Application of a novel time-delayed polynomial grey model to predict the natural gas consumption in China. J. Comput. Appl. Math. 2017, 324, 17–24. [Google Scholar] [CrossRef]
  28. Wer, B.; Xie, N.; Hu, A. Optimal solution for novel grey polynomial prediction model. Appl. Math. Model. 2018, 62, 717–727. [Google Scholar]
  29. Wang, Z.; Li, Q.; Pei, L. A seasonal GM(1,1) model for forecasting the electricity consumption of the primary economic sectors. Energy 2018, 154, 522–534. [Google Scholar] [CrossRef]
  30. Zeng, B.; Duan, H.; Bai, Y.; Meng, W. Forecasting the output of shale gas in China using an unbiased grey model and weakening buffer operator. Energy 2018, 151, 238–249. [Google Scholar] [CrossRef]
  31. Wang, J.; Du, P.; Lu, H.; Yang, W.; Niu, T. An improved grey model optimized by multi-objective ant lion optimization algorithm for annual electricity consumption forecasting. Appl. Soft Comput. 2018, 72, 321–337. [Google Scholar] [CrossRef]
Figure 1. Prediction error source diagram of conventional GM(1,1) model.
Figure 1. Prediction error source diagram of conventional GM(1,1) model.
Symmetry 11 00420 g001
Figure 2. Graphic results for Example 1.
Figure 2. Graphic results for Example 1.
Symmetry 11 00420 g002
Figure 3. Graphic results for Example 2.
Figure 3. Graphic results for Example 2.
Symmetry 11 00420 g003
Figure 4. Graphic results for Example 3.
Figure 4. Graphic results for Example 3.
Symmetry 11 00420 g004
Figure 5. Graphic results for Example 4.
Figure 5. Graphic results for Example 4.
Symmetry 11 00420 g005
Figure 6. Graphic results for Example 5.
Figure 6. Graphic results for Example 5.
Symmetry 11 00420 g006
Figure 7. Graphic results for Example 6.
Figure 7. Graphic results for Example 6.
Symmetry 11 00420 g007
Table 1. Numerical results for Example 1.
Table 1. Numerical results for Example 1.
x ( 0 ) Classical GM(1,1)New GM(1,1)The Model in [20]
Prediction Data ε (%)Prediction Data ε (%)Prediction Data ε (%)
2.98362.983602.983602.98360
4.45114.38041.58164.35312.20214.45610.1123
6.64026.50062.09036.52221.78356.61320.4066
9.90619.64692.59949.77201.35699.81460.9237
14.778114.31623.103914.64130.934414.56571.4373
22.046421.24543.606921.93680.501321.61681.9486
32.889331.52854.106532.86750.079332.08122.4570
ε ¯ (%)2.84811.14291.2143
Table 2. Numerical results for Example 2.
Table 2. Numerical results for Example 2.
x ( 0 ) Classical GM(1,1)New GM(1,1)The Model in [21]
Prediction Data ε (%)Prediction Data ε (%)Prediction Data ε (%)
21.121.1021.1021.10000
26.624.41668.208324.07799.481623.360612.1782
36.135.71981.053132.56489.792835.78580.8704
52.352.25570.084750.31163.801854.81984.8180
80.176.44664.561177.73002.958883.97774.8411
126.8111.836111.8012120.09065.2913128.64431.4545
196.3163.608716.6537185.53655.4832197.06840.3914
ε ¯ (%)7.06046.13494.0923
Table 3. Numerical results for Example 3.
Table 3. Numerical results for Example 3.
x ( 0 ) Classical GM(1,1)New GM(1,1)The Model in [7]
Prediction Data ε (%)Prediction Data ε (%)Prediction Data ε (%)
110,852110,8520110,8520110,8520
135,175117,98012.72130,0783.77127,8215.41
153,647119,11722.47128,57416.32126,66417.66
120,296128,2646.62127,0875.64125,8304.68
96,362121,42226.27125,61730.36124,38029.23
90,798122,59235.01124,16436.75123,25335.70
102,591123,77320.65122,72819.63122,13719.11
150,534124,96516.99121,30819.41121,03119.63
175,123126,16827.95119,90531.53119,93431.52
127,148113,38310.83118,5186.79114,8489.76
102,085128,61025.98117,14714.75117,77215.47
97,103116,70520.19115,79219.25116,70520.21
ε ¯ (%)21.7318.5618.91
Table 4. Numerical results for Example 4.
Table 4. Numerical results for Example 4.
x ( 0 ) Classical GM(1,1)New GM(1,1)The Model in [9]
Prediction Data ε (%)Prediction Data ε (%)Prediction Data ε (%)
5.085.0805.0805.080
4.803.5526.044.2411.513.6523.9
4.674.523.215.2512.574.611.28
4.505.7628.005.5122.445.8429.78
7.127.343.098.0513.097.424.21
9.679.343.419.973.059.452.28
12.8011.907.0312.333.6412.045.94
15.8815.154.5915.263.8815.373.21
19.4919.300.9718.893.0719.650.82
23.0724.576.5023.381.3424.578.92
26.8631.2916.4928.937.7231.2919.73
ε ¯ (%)9.037.488.29
Table 5. Numerical results for Example 5.
Table 5. Numerical results for Example 5.
x ( 0 ) Classical GM(1,1)New GM(1,1)The Model in [11]
Prediction Data ε (%)Prediction Data ε (%)Prediction Data ε (%)
1610.711610.7101610.7101610.710
1245.281363.919.5241365.229.6321342.767.833
1347.711274.955.4021275.755.3401280.325.004
1382.451191.7913.7951192.1413.7661204.5812.872
1018.451114.069.3811114.019.3831122.6410.203
1014.961041.392.5971041.002.5661040.182.476
949.46973.472.519972.782.456960.341.147
ε ¯ (%)6.1746.1635.648
Table 6. Numerical results for Example 6.
Table 6. Numerical results for Example 6.
x ( 0 ) Classical GM(1,1)New GM(1,1)
Prediction Data ε (%)Prediction Data ε (%)
64,832.0564,832.05064,832.050
71,847.0957,476.7720.00157,308.5020.235
78,646.3067,165.2114.59867,022.8314.779
86,293.1078,486.769.04678,383.839.166
93,887.9591,716.702.31291,670.622.362
105,557.09107,176.711.534107,209.651.566
125,761.85125,242.710.413125,382.700.301
143,143.63146,353.962.243146,636.242.440
168,850.20171,023.781.287171,492.461.565
198,739.27199,852.010.560200,562.040.917
245,352.80233,539.604.815234,559.184.400
278,541.09272,905.572.023274,319.161.516
334,839.41318,907.394.758320,818.834.187
386,086.72372,663.283.477375,200.632.820
ε ¯ (%)4.7914.732

Share and Cite

MDPI and ACS Style

Zhu, Y.; Jian, Z.; Du, Y.; Chen, W.; Fang, J. A New GM(1,1) Model Based on Cubic Monotonicity-Preserving Interpolation Spline. Symmetry 2019, 11, 420. https://doi.org/10.3390/sym11030420

AMA Style

Zhu Y, Jian Z, Du Y, Chen W, Fang J. A New GM(1,1) Model Based on Cubic Monotonicity-Preserving Interpolation Spline. Symmetry. 2019; 11(3):420. https://doi.org/10.3390/sym11030420

Chicago/Turabian Style

Zhu, Yuanpeng, Zehua Jian, Yurui Du, Wenqing Chen, and Jiwei Fang. 2019. "A New GM(1,1) Model Based on Cubic Monotonicity-Preserving Interpolation Spline" Symmetry 11, no. 3: 420. https://doi.org/10.3390/sym11030420

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop