Frontiers in DataDriven Methods for Understanding, Prediction, and Control of Complex Systems
View this Special IssueReview Article  Open Access
Lifeng Wu, Xiaohui Gao, Yan Chen, "Memory Property of Grey Accumulation Generation Sequence", Complexity, vol. 2019, Article ID 9459038, 10 pages, 2019. https://doi.org/10.1155/2019/9459038
Memory Property of Grey Accumulation Generation Sequence
Abstract
The influence of grey accumulation generation operator on the memory of time sequence is discussed. Generalized Hurst exponent (GHE) approach is used to calculate the memory in three cases, namely, the M3Competition data, the air quality index in Beijing, Xingtai, and Handan, and the power generating capacity, and car production index in Hebei province. The result indicates that one order accumulation generation operator (1AGO) can weaken the volatility and strengthen the memory of time sequence. It also explains the reason that oneorder accumulation can be used in grey prediction.
1. Introduction
Since the introduction of grey prediction theory, it is has been widely applied in many fields, such as Beijing’s tertiary industry [1], natural gas demand in China [2], return flow of endoflife vehicles [3], etc. Typically, a novel multivariable grey forecasting model was proposed to forecast the number of motor vehicles in Beijing [4]. By using a selfadapting intelligent grey model, the forecasting results were substantially improved [5]. A grey model, based on a data grouping approach, was used to predict quarterly hydropower production in China [6]. The nonlinear grey multivariable model was supposed to be potentially suitable for forecasting carbon emissions due to fossil energy consumption [7]. To explore the factors causing smog weather, the grey dynamic trend incidence model was developed [8]. The grey prediction model with a nonlinear optimized time response method was presented to forecast electricity consumption in China [9].
The AGO can change the process from grey to white. It also occupies an extremely important position in the grey system theory [10]. It can reveal the patterns (laws) in chaotic original data and smooth the irregularities in the raw data. Therefore, the AGO is an indispensable part of grey forecasting [11]. To enhance the prediction accuracy, many scholars have made significance contribution to the problem. Using the local sequence grey accumulation, the forecasting accuracy can be improved to some degree [12]. The seasonal rolling grey forecasting model with a cycle truncation AGO is a complement to grey prediction [13]. The fractionalorder accumulation can reduce errors from the inverse AGO of the grey model [14]. However, no study has been conducted to determine whether the grey AGO can change the memory of time sequence. Because the Hurst exponent can represent the memory of time sequence, the GHE is utilized to estimate the memory in this paper. The influence of the AGO on the memory of time sequence is analyzed through a series of experimental cases.
The rest of this paper is organized as follows. The proposed calculation method is described in Section 2. The three experimental cases are detailed in Section 3. The comparison of firstorder single variable grey forecasting model and other models is shown in Section 4. The conclusion is presented in Section 5.
2. Method
The GHE was presented to calculate the memory [15]. It is defined aswhere is the length of time sequence and is the statistic. Equation (1) can yield a series of values by changing the value of . The relationship between and can be described as follows:Thus, we can estimate the Hurst exponent of time sequence in (1)(2) for . When , the time sequence reveals its antipersistence. In addition, the closer H to 0, the greater intensity of this antipersistence behavior [16]. When , the time sequence follows a random walk [17–22]. When , the time sequence implies strong persistence. When , the tendency of future can be predicted by current state [23].
The process of estimation the Hurst exponent is shown as follows.
Step 1. For positive time sequence , the oneorder accumulation of the time sequence can be obtained byEquation (3) is called the 1AGO.
Step 2. Calculate a series of by using the following:
Step 3. According to the formula, we have The leastsquare method is utilized to estimate the Hurst exponent .
3. Experiments and Analysis
All the experiments in this study were conducted by using the software R i386 3.4.1. The Hurst exponent values and the corresponding analyses of three cases are described in detail below.
Case 1. The data are from the M3Competition (N2831, N2832, and N2833) and are downloaded from the Internet site: https://cran.rproject.org/web/packages/Mcomp/. The oneorder accumulation sequence is presented in Figure 1. The Hurst exponent value of three groups of data is given in Table 1. The memory of three groups of data is shown in Figure 2.
In Figure 1, the original data from N2831 to N2833 have different distribution and strong randomness. However, oneorder accumulation sequence from N2831 to N2833 nearly distributes in a straight line; it implied that oneorder accumulation sequence has substantial developing tendency and stronger regularity.
As shown in Table 1, the values of the original sequence are less than 0.5; this indicated that the original time sequence reveals its antipersistence. In contrast, the values of accumulation sequence are all larger than 0.90; it implied that the memory of time sequence can be strengthened and the randomness of the time sequence can be weakened by AGO.
It can be noted in Figure 2 that the accumulation sequence of time sequence is regularly distributed in a straight line. The value of accumulation sequence is larger than that of the original time sequence.

(a) The original value of N2831
(b) The 1AGO of N2831
(c) The original value of N2832
(d) The 1AGO of N2832
(e) The original value of N2833
(f) The 1AGO of N2832
(a) The memory of original N2831 value
(b) The memory of 1AGO N2831 value
(c) The memory of original N2832 value
(d) The memory of 1AGO N2832 value
(e) The memory of original N2833 value
(f) The memory of 1AGO N2833 value
Case 2. Beijing, Handan, and Xingtai are the three cities selected for comparison, where the air is seriously polluted. The air quality index (AQI) data are collected from http://www.zhb.gov.cn/. Since 2013, the air pollution has been noticed by citizens. The government has took urgent measures to control this situation. The result of in the three cities is in Table 2. The AQI of Beijing, Handan, and Xingtai from 2014 to 2016 are given in Figure 3. The memory results of the three group data are plotted respectively in Figure 4.
It can be seen in Figure 3 that the AQI values in winter are higher than that of in the other seasons. In this situation, the AQI time sequence in Beijing, Handan, and Xingtai has fluctuated randomly. However, oneorder accumulation of AQI has stronger regularity. It is possible to establish a suitable model.
As can be seen in Table 2, the values of the original sequence are all close to 0; it implied that the original AQI time sequence have quite strong antipersistence and mutability. However, the values of the accumulation sequence are all close to 1.0; it indicated that the accumulation AQI time sequence has higher regularity. Consequently, it can be more extensively applied for prediction in various fields.
It is observed in Figure 4 that the original AQI sequence is disperse randomly while the accumulation sequence is regularly distributed in a straight line; it revealed that the value of the accumulation sequence is larger than that of the original time sequence.

(a) The AQI in Beijing
(b) The 1AGO of Beijing AQI
(c) The AQI in Handan
(d) The 1AGO of Handan AQI
(e) The AQI in Xingtai
(f) The 1AGO of Xingtai AQI
(a) The memory of original Beijing AQI
(b) The memory of 1AGO Beijing AQI
(c) The memory of original Handan AQI
(d) The memory of 1AGO Handan AQI
(e) The memory of original Xingtai AQI
(f) The memory of 1AGO Xingtai AQI
Case 3. Hebei province is located in the North China Plain. Its power generating capacity and car production index are considered to analyze the memory. These data are collected from National Bureau of Statistics of China. The values of the two industry indexes in Hebei province are listed in Table 3. The original value and oneorder accumulation value of two indexes are presented in Figure 5. The memory results of two indexes in Hebei province are plotted in Figure 6.
It is implied that the tendency of both indexes are increasing slightly in Figure 5. Furthermore, oneorder accumulation sequence of both indexes has strong regularity. It is clear that the trend of on order accumulation sequence is more obvious than the original sequence.
As seen in Table 3, the values of the accumulation sequence are larger than those of the original time sequence for both indexes.
In Figure 6, it is noted that the points of the original value are scattered whereas those of the accumulation sequence are evenly distributed in a straight line. Therefore, the accumulation sequence yields superior persistence.
Actually, for the big size cases, the memory of sequence is enhanced by the 1AGO. For the small size cases, the memory of sequence is also enhanced by the 1AGO.

(a) The original car production
(b) The 1AGO of car production
(c) The original power capacity
(d) The 1AGO of power capacity
(a) The memory of original car production
(b) The memory of lAGO car production
(c) The memory of power capacity
(d) The memory of 1AGO power capacity
4. Comparison of GM(1,1) Model and the Other Models
As the basis of the grey forecasting theory, GM(1,1) is the firstorder single variable grey forecasting model, and it has been widely investigated and studied by many scholars in theoretical and practical fields [24]. The steps of GM(1,1) are shown in the following.
Step 1. An original positive sequence is . Oneorder accumulation sequence is .
Step 2. The GM(1,1) model is , , .
Its whitenization equation is . and are coefficients. is the developing coefficient, and is the grey input. and are obtained by the leastsquares method as follows:where
Step 3. The solution of is . According to , is obtained as .
Step 4. To obtain the predicted value of the data at time , the forecasting sequence is and the inverse accumulating generation operation is applied as follows: .
GM(1,1) is one of the most common methods used for studying data sets characterized by small size. The above cases have big size; GM(1,1) is not the better model for above cases. Therefore, the series N1 (sample size is 14) from the M3 forecasting competition is selected as the evaluation data. The data is downloaded from the Internet site: http://www.insead.fr/facultyresearch/forecasting.
The memory of accumulation sequence is evaluated by comparing the traditional GM(1,1) model with the other models in this section. Therefore, mean absolute percentage error () compares the actual values with the forecasting values. The last 6 observations are used to evaluate their accuracy. The forecasting results of the series N1 are listed in Table 4 and are plotted in Figure 7.

Winter model, ARAMRA, and AUTOANN are the most common methods used for studying data sets characterized by small size. In Table 4, it indicated that the MAPE of GM(1,1) are smaller than that of other models. We can say that the GM(1,1) model is better than the others. The original time sequence always shows strong antipersistent behavior. It is difficult to dig out the information in the sequence. The AGO can strengthen the memory of the time sequence. Therefore, it plays a pivotal role to establish model. With the application of AGO in GM(1,1) model, the case imply that GM(1,1) model has higher forecasting performance.
5. Conclusion
For a time sequence, the current state is the result of a series of interrelated events. The development state has a longterm (shortterm) memory. In general, the Hurst exponent of original time sequence is less than or close to 0.5. In this situation, the time sequence always shows stronger antipersistent behavior; it means that the time sequence follows a random walk with mutability or volatility. As can be seen from the experimental cases, the Hurst exponent of oneorder accumulation sequence was very close to 1, the randomness of the time sequence was significantly weakened and the regularity was strengthened. By the 1AGO, the memory of sequence is enhanced and the accumulation sequence yields superior persistence. The future trend can be predicted based on the present state.
It is worthwhile to use 1AGO in grey prediction. On one hand, the randomness of the original sequence is weakened by 1AGO. On the other hand, grey accumulation generation sequence can provide additional information for modeling and enable the data to construction a more appropriate model. Because the Hurst exponent of one order accumulation sequence was very close to the maximum (i.e., 1), for the grey forecasting model, the accumulationorder number cannot be more than 1 in practical applications.
Three cases are discussed to demonstrate the memory of grey AGO. In the future, more cases are needed to be discussed to demonstrate the memory. All values in this paper are positive, so the accumulation frequency is always increasing. In this case, the Hurst Exponent is close to 1. Of course, this property may change when the original time sequence contains negative values. In the future, more cases with negative values are needed to be discussed.
Conflicts of Interest
All authors declare that they have no conflicts of interest.
Acknowledgments
The relevant researches in this paper are supported by the National Natural Science Foundation of China (Nos. 71871084, 71401051), the Excellent Young Scientist Foundation of Hebei Education Department (No. SLRC2019001), and the project of highlevel talent in Hebei province.
References
 Q. Wang, L. Liu, S. Wang, J.Z. Wang, and M. Liu, “Predicting Beijing's tertiary industry with an improved grey model,” Applied Soft Computing, vol. 57, pp. 482–494, 2017. View at: Publisher Site  Google Scholar
 F. Shaikh, Q. Ji, P. H. Shaikh, N. H. Mirjat, and M. A. Uqaili, “Forecasting China’s natural gas demand based on optimised nonlinear grey models,” Energy, vol. 140, pp. 941–951, 2017. View at: Publisher Site  Google Scholar
 S. Ene and N. Öztürk, “Grey modelling based forecasting system for return flow of endoflife vehicles,” Technological Forecasting & Social Change, vol. 115, pp. 155–166, 2017. View at: Publisher Site  Google Scholar
 B. Zeng, C. Luo, C. Li, and S. Liu, “A novel multivariable grey forecasting model and its application in forecasting the amount of motor vehicles in Beijing,” Computers & Industrial Engineering, vol. 101, pp. 479–489, 2016. View at: Publisher Site  Google Scholar
 B. Zeng and C. Li, “Forecasting the natural gas demand in China using a selfadapting intelligent grey model,” Energy, vol. 112, pp. 810–825, 2016. View at: Publisher Site  Google Scholar
 Z.X. Wang, Q. Li, and L.L. Pei, “Grey forecasting method of quarterly hydropower production in China based on a data grouping approach,” Applied Mathematical Modelling, vol. 51, pp. 302–316, 2017. View at: Publisher Site  Google Scholar
 Z.X. Wang and D.J. Ye, “Forecasting Chinese carbon emissions from fossil energy consumption using nonlinear grey multivariable models,” Journal of Cleaner Production, vol. 142, pp. 600–612, 2017. View at: Publisher Site  Google Scholar
 J. Wang, K. W. Hipel, and Y. Dang, “An improved grey dynamic trend incidence model with application to factors causing smog weather,” Expert Systems with Applications, vol. 87, pp. 240–251, 2017. View at: Publisher Site  Google Scholar
 N. Xu, Y. Dang, and Y. Gong, “Novel grey prediction model with nonlinear optimized time response method for forecasting of electricity consumption in China,” Energy, vol. 118, pp. 473–480, 2017. View at: Publisher Site  Google Scholar
 L. Wu, S. Liu, D. Chen, L. Yao, and W. Cui, “Using gray model with fractional order accumulation to predict gas emission,” Natural Hazards, vol. 71, no. 3, pp. 2231–2236, 2014. View at: Publisher Site  Google Scholar
 C.J. Chang, J.Y. Lin, and P. Jin, “A grey modeling procedure based on the data smoothing index for shortterm manufacturing demand forecast,” Computational and Mathematical Organization Theory, vol. 23, no. 3, pp. 409–422, 2017. View at: Publisher Site  Google Scholar
 K. L. Wen, Grey prediction theory and application, Open Tech Company, Taipei, Taiwan, 2002.
 X. Xiao, J. Yang, S. Mao, and J. Wen, “An improved seasonal rolling grey forecasting model using a cycle truncation accumulated generating operation for traffic flow,” Applied Mathematical Modelling, vol. 51, pp. 386–404, 2017. View at: Publisher Site  Google Scholar
 L. Wu, S. Liu, L. Yao, R. Xu, and X. Lei, “Using fractional order accumulation to reduce errors from inverse accumulated generating operator of grey model,” Soft Computing, vol. 19, no. 2, pp. 483–488, 2014. View at: Publisher Site  Google Scholar
 T. Di Matteo, T. Aste, and M. M. Dacorogna, “Scaling behaviors in differently developed markets,” Physica A: Statistical Mechanics and its Applications, vol. 324, no. 12, pp. 183–188, 2003. View at: Publisher Site  Google Scholar
 J. P. RamosRequena, J. E. TrinidadSegovia, and M. A. SánchezGranero, “Introducing Hurst exponent in pair trading,” Physica A: Statistical Mechanics and its Applications, vol. 488, pp. 39–45, 2017. View at: Publisher Site  Google Scholar
 L. Salim and S. Amir, “Accurate classification of seizure and seizurefree intervals of intracranial EEG signals from epileptic patients,” IEEE Transactions on Instrumentation, vol. 68, no. 3, pp. 791–796, 2019. View at: Google Scholar
 S. Lahmiri, “An accurate system to distinguish between normal and abnormal electroencephalogram records with epileptic seizure free intervals,” Biomedical Signal Processing and Control, vol. 40, pp. 312–317, 2018. View at: Publisher Site  Google Scholar
 L. Salim, “Glioma detection based on multifractal features of segmented brain MRI by particle swarm optimization techniques,” Biomedical Signal Processing and Control, vol. 31, pp. 148–155, 2017. View at: Publisher Site  Google Scholar
 T. Blachowicz, A. Ehrmann, and K. Domino, “Statistical analysis of digital images of periodic fibrous structures using generalized Hurst exponent distributions,” Physica A: Statistical Mechanics and its Applications, vol. 452, pp. 167–177, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 M. Ausloos, “Generalized Hurst exponent and multifractal function of original and translated texts mapped into frequency and length time series,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 86, no. 3, Article ID 031108, 2012. View at: Google Scholar
 C. Dou, X. Wei, and J. Lin, “Fault diagnosis of gearboxes using nonlinearity and determinism by generalized hurst exponents of shuffle and surrogate data,” Entropy, vol. 20, no. 5, p. 364, 2018. View at: Publisher Site  Google Scholar
 W. Wang, Long Memory Characteristics And Prediction of Financial Time Series, Nankai University Press, 2018.
 B. Zeng, H. Duan, and Y. Zhou, “A new multivariable grey prediction model with structure compatibility,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 75, pp. 385–397, 2019. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2019 Lifeng Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.