Abstract

Traditional discrete grey forecasting model can effectively predict the development trend of the stabilizing system. However, when the system has disturbance information, the prediction result will have larger error, and there will appear significant downward trend in the stability of the model. In the presence of disturbance information, this paper presents a fractional-order linear time-varying parameters discrete grey forecasting model to deal with the system that contains both linear trend and nonlinear trend. The modeling process of the model and calculation method are given. The perturbation bounds of the new model are analyzed by using the least-squares method of perturbation theory. And it is compared with that of the first-order linear time-varying parameters discrete grey forecasting model. Finally, two real cases are given to verify the effectiveness and practicality of the proposed method.

1. Introduction

The grey system theory was proposed by Deng in 1982. It only needs a small amount of data to establish the model for analysis [1, 2]. Grey prediction model is the main part of grey system theory, and its theoretical basis is the grey of accumulation generation. The prediction model is established by using the characteristic of grey exponential. Since the grey prediction model was proposed, it has been widely concerned by scholars. The existing studies mainly focus on the theoretical development from the aspects of background value optimization [3, 4], parameter optimization [57], model expansion [810], and the application of natural gas [11, 12], electric power [1316], environment [1720], economy [21, 22], and transportation [23, 24].

Dang gave the new method to select initial value [25]. Cui proposed the NGM(1,1,k) model and described the modeling mechanism and modeling process of the model [26]. Zhou proposed the generalized GM(1,1) model and used the new model to simulate and predict China's fuel output from 2003 to 2010 [27]. Combined with the concept of Bernoulli differential equation, Chen proposed the NGBM(1,1) model on the basis of GM(1,1) model [28]. Li proposed 3spGM(1, 1) model and applied the new model to the failure data sets of electric product manufacturing systems [29]. These methods further improve the modeling effect of grey prediction model. However, the transformation of difference equation and differential equation is required in the solving process. There is still error in the grey prediction model for the sequences that conform to the exponential features.

Xie proposed the discrete grey model (DGM(1,1) model). The relationship between DGM(1,1) and GM(1,1) was studied deeply, and the reason for the instability of GM(1,1) was found. It only needs to use difference equation to solve the equation and does not need to convert the difference equation to differential equation. Therefore, the modeling accuracy is effectively improved [30]. Wu proposed the discrete grey prediction model based on fractional-order accumulation, discussed the properties of the model, and gave the calculation method of the model [3133]. Liu proposed the fractional-order reverse accumulation discrete grey prediction model and discussed the properties of the model [34].

The research mentioned above has positive significance for improving the accuracy of grey prediction model. However, for complex systems with disturbance, the model’s robustness is insufficient; how to deal with the system disturbance information is particularly important. To solve this problem, this paper presents a fractional-order accumulation linear time-varying parameters discrete grey prediction model (FTDGM(1,1)model). The modeling process and parametric calculation method of the model are given. It is proven that the model has good stability using the theory of matrix perturbation analysis. Finally, two real cases are given. And the calculation results show that FTDGM(1,1) model can effectively reduce the disturbance caused by disturbance information. The model's robustness and prediction accuracy are improved, and the validity and practicality of the model are further verified.

2. The Fractional-Order Accumulative Linear Time-Varying Parameters Discrete Grey Model

Definition 1 (see [30]). Assume that the nonnegative sequence is . is the first-order accumulative sequence of .
Among them, The equation is called discrete grey prediction model (DGM(1,1) model).

Theorem 2 (see [30]). The parameters of the DGM(1,1) model can be solved by using the following least-squares estimation:and, among them,

Definition 3 (see [35]). Assume that the nonnegative sequence is . is the first-order accumulative sequence of .
Among them,The equationis called linear time-varying parameters discrete grey model (TDGM).

Theorem 4. The parameters of the TDGM(1,1) model can be solved by using the following least-squares estimation:and, among them,

Definition 5 (see [18]). Assume that the nonnegative sequence is .
is called the fractional-order accumulative sequence of
Among them,

Definition 6. Assume the nonnegative sequence ; is defined as Definition 5. The equationis called fractional-order accumulative linear time-varying parameters discrete grey model (FTDGM).

Theorem 7. The parameters of the FTDGM(1,1) model can be solved by using the following least-squares estimation:and, among them,The predicted value of FTDGM(1,1)model is as follows:According to the calculation formula of fractional-order accumulation, it is not difficult to calculate the reduced value of the predicted sequence as follows:

3. Disturbance Analysis of TDGM(1,1) Model and FTDGM(1,1) Model

Theorem 8 (see [36, 37]). Suppose that is the generalized inverse matrix of . When the column vector of has linear independence, the function has a unique solution.

Theorem 9 (see [36, 37]). Suppose that is the generalized inverse matrix of . . Suppose that the solutions of function and are and , respectively. When and , we have the following result:Among them,

3.1. Disturbance Analysis of TDGM(1,1) Model

The perturbation bounds of TDGM(1,1) model will be analyzed in the following section.

Theorem 10. The solution of TDGM model can be given as the following function: . Suppose that the solution of the TDGM(1,1) model is , and . Among them, is the disturbance information.
Then

Proof. Assume that the solution of the new model is , and the disturbance is . The equation has a unique solution due to the linear independence of column vectors of .
Sincewe haveThen, the following result can be obtained according to Theorem 9:

Theorem 11. Assume that the conditions of Theorem 9 remain unchanged, and .
Then the perturbation bound of the solution is as follows:

Proof. Assume that the solution of the new model is , and the disturbance is . The equation has a unique solution due to the linear independence of column vectors of .
Sincewe haveThen, the following result can be obtained according to Theorem 9:

3.2. Disturbance Analysis of TDGM(1,1) Model and FTDGM(1,1) Model

Theorem 12. The solution of TDGM(1,1) model can be given as the following function: . Suppose that the solution of the TDGM(1,1) model is , and . Among them, is the disturbance information.
Then

Proof. Assume that the solution of the new model is , and the disturbance is . The equation has a unique solution due to the linear independence of column vectors of .
Sincewe haveThen, the following result can be obtained according to Theorem 9:Sincewe haveIt is not hard to get .

Theorem 13. The solution of TDGM(1,1) model can be given as the following function: . Suppose that the solution of the TDGM(1,1) model is , and . Among them, is the disturbance information.
Then

Proof. We haveThen, the following result can be obtained according to Theorem 9.We haveThen

It can be seen from the conclusion of Theorems 12 and 13 that when , the disturbance bound of FTDGM(1,1) model is smaller than that of TDGM(1,1) model. Generally speaking, compared with the TDGM(1,1) model, the FTDGM(1,1) model has better robustness. The proposed FTDGM(1,1) model can effectively reduce the prediction error caused by system disturbance and improve the prediction accuracy of the grey forecasting model. And the solution of optimal parameter can be given by genetic algorithms.

4. Numerical Illustrations

In order to test the modeling effect of the model, two real cases will be given in the following section.

Case 1. The GDP of Guangdong province in 2001-2009 is used to build different grey prediction model. The data is shown in Table 1. Different grey prediction models were established based on the given data. And the advantages of the model were tested by comparing the prediction accuracies of different models.
The highest prediction accuracy in [38] is time-varying parameters grey model (TVGM(1,1) model). This Paper established DGM(1,1) model, TVGM(1,1) model, and FTDGM(1,1) model, respectively. The calculation results are shown in Table 2 and Figure 1. The parameters of FTDGM(1,1) model are as follows: It can be seen from Table 2 that, because of the nature of exponential function, large errors will appear after the second prediction of the DGM(1,1) model. TVGM(1,1) model can effectively increase the prediction accuracy along with adjusting and optimizing the time-varying parameters. However, the prediction error is still relatively large. The proposed FTDGM(1,1) model in this paper shows good robustness, and the prediction error does not increase with the passage of time. It indicated that the FTDGM(1,1) model has better anti-interference and long-term memory, which can be used to predict medium-term goals.

Case 2. As one of the important symbols of China's transportation modernization, highway reflects the degree and level of a country's modernization. Compared with railway and air or water transportation, highway transportation is the more used mode in passenger and cargo transportation. Highway transportation is point-to-point direct, flexible, and convenient. It is very important to accurately predict the length of highway transportation route. The length of Chinese highway transportation route in 2010-2017 is used to build different grey prediction model. The unit of highway mileage is ten thousand kilometers. The calculation results of different models are shown in Table 3 and Figure 2. The parameters of FTDGM(1,1) model are as follows: ,

It can be seen from the calculation results that the TDGM(1,1) model is better than the DGM(1,1) model in describing the internal evolution. However, neither the DGM(1,1) model nor the TDGM(1,1) model can accurately describe the development trend of the system. In particular, the prediction error of the second step is relatively large, indicating that the memory of the integer-order model is insufficient for the FTDGM(1,1) model. The prediction error of the FTDGM(1,1) model is only 0.64%, which shows that the model has strong extrapolation ability and good memory.

5. Concluding Remarks

On the basis of the traditional discrete grey prediction model, this paper proposed the FTDGM(1,1) model. The parametric solution method of the model was given. By using the matrix disturbance theory, the disturbance boundary of the model was analyzed, and it was proved that the FTDGM(1,1) model had better robustness than the TDGM(1,1) model.

Two real cases were used to test the effect of the proposed FTDGM(1,1) model. It was found that the prediction accuracy is higher than that of the existing models, which further verifies the superiority of the proposed FTDGM(1,1) model. At the same time, the FTDGM(1,1) model had better memory than other discrete grey prediction models. Research results of this paper further expand the application scope of the grey prediction model. And the reasons for the existence of short-term memory model deserve further discussion.

Data Availability

The data are from China Statistical Yearbook, website of China’s national bureau of statistics, http://www.stats.gov.cn/, Guangdong Statistical Yearbook, and Guangdong Statistical Information Network, http://www.gdstats.gov.cn/tjsj/gmjjhs/.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (71871084), 2016 Shanghai Maritime University National Key Project Development Project (A20201161107X), and 2015 Shanghai Maritime University Research Start-Up Funding Project (A15101154505Z). At the same time, the authors would like to acknowledge the support of the project of Chinese postdoctoral science foundation (2018M632777) and the training program of youth backbone teacher in institutions and universities in Henan province (2018GGJS115).