1. Introduction
Let an original non-negative and uniformly spaced sequence be
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
. Then, establish a first-order grad forecasting differential equation on the sequence
. 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
is given as follows
where
From Equation (3), we can see that the 1-AGO sequence
has the property of monotonicity-increasing. Suppose that
meets the following first-order grad forecasting differential equation
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
is as follows
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
, then we can get
that is
or
Let background value be
. To calculate the background value
, we need to integrate
, 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
. We use the piecewise linear polynomial interpolation
to approximate
in the classical GM(1,1) model, see [
1,
2], then we get the estimated background value
as follows
For each interval
,
, by substituting the estimated background value
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
where
Finally, we get the following estimated solution to the differential Equation (4) with the initial condition
as follows
We thus get the following grey prediction equation
From (7), it can be found that the classical GM(1,1) model uses average of adjacent values to estimate the background value
. 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
, 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 (
) 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
. In [
6], Li and Dai reconstructed
by a high-order Newton interpolation polynomial. In addition, they estimated the background value
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
by using the piecewise quadratic interpolation polynomial to reconstruct
on the interval
. 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
and thus got the estimated background value
. 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
to be reconstructed. If the resulting reconstructed curve
loses the monotonicity-increasing feature, there will also result in a large error on the background value
. Therefore, an accurate approximation of the reconstructed curve
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
and thereby give a new scheme to estimate the background value
. The rest of this paper is structured as follows.
Section 2 represents the construction of the
monotonicity-preserving cubic interpolation splines. In
Section 3, based on
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. Monotonicity-Preserving Piecewise Cubic Interpolation Spline
According to Equation (3), the 1-AGO sequence has the property of monotonicity-increasing, that is , . The fitting exponential curve to the 1-AGO sequence is also monotonicity-increasing and has infinite smoothness. Therefore, we develop a monotonic-preserving cubic interpolation spline to interpolate the 1-AGO sequence, to reconstruct the curve .
For the discrete data
, we denote
as the derivative value at node
. For
, a cubic interpolation spline with local parameter
is constructed as follows
where
. From Equation (10), after some computations, we have
which indicates that
. This means that the cubic interpolation spline defined by Equation (10) is
continuous for arbitrary nonzero local parameter. Here,
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
, 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
We shall derive sufficient conditions for the
cubic interpolation spline preserving monotonicity. For the monotonicity-increasing 1-AGO sequence
, it is obvious that the derivative value determined by Equation (12) is non-negative, which means
. Without loss of generality, for
, direct computation gives that
Thus, we can see that the following conditions are sufficient to guarantee
Summarizing the above discussion, we can obtain the following sufficient conditions for
preserving monotonicity
where
serves as free parameter. In practice, we recommend all
. If there is no special explanation below, we will take all
.
3. Establish New GM(1,1) Model
For the original non-negative sequence
, we first calculate its 1-AGO sequence
. Then for the 1-AGO sequence
, we use the sufficient conditions (14) with all the free parameters
to determine a
monotonicity-preserving piecewise cubic interpolation spline
to reconstruct the exponential curve
. For the interval
, we estimate the background value
by the following method
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)
where
Finally, we get the following estimated solution to the differential Equation (4) with the initial condition
as follows
We thus get the following grey prediction equation
We shall give several examples to show that the new GM(1,1) model based on
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
Example 1. In this example, we use the exponential function 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 and the reconstruct exponential curve . From the right Figure 2, we can see that the interpolation spline has 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.