# Prediction of Air Quality Based on Hybrid Grey Double Exponential Smoothing Model

Authors: Zhicun Xu; Meng Dun; Lifeng Wu
Journal: Complexity
 (2020)

---

## Abstract

To predict the concentration of air pollutants accurately, conformable fractional accumulation grey double exponential smoothing and adjacent accumulation grey double exponential smoothing are proposed, respectively. A hybrid model with conformable fractional accumulation and adjacent accumulation is proposed to further improve the prediction accuracy. The hybrid model is applied to the prediction of AQI, PM2.5, and PM10. Compared with the traditional double exponential smoothing model, the new model has a higher prediction accuracy.


---

## Body

## 1. Introduction

Air quality pollution can induce a variety of respiratory diseases and seriously endanger human health. Therefore, environmental pollution is highly valued. A series of policies to deal with environmental issues have also been proposed, such as the “Air Pollution Prevention and Control Action Plan” and “Notice of Winning the Three-Year Plan of Action to Defend the Blue Sky.” With the vigorous control of the country, the air pollution problem has gradually improved. In order to have a better air level, it is necessary to continue to study air quality issues. Due to the promulgation of the new air quality standards, there are less available data on air quality under the new standards. Grey system theory provides an opportunity for accurate prediction of air quality due to its characteristics of “less data” and “poor information.” Since Professor Deng published the first grey system paper “The control problems of grey systems” in “Systems & Control Letters” in 1982, this paper symbolizes the advent of the discipline of grey systems theory [1]. Since the emergence of the grey system theory, it has been widely used in various fields, such as construction demand [2], evaluating the interaction between cryptocurrencies [3], traffic-related emissions [4], energy consumption [5], environmental protection investment [6], medicine and health [7], and missile development cost [8]. With the development of grey system theory, different scholars have improved and expanded the grey system theory in terms of cumulative generation methods, parameters, and background values. For example, a grey multivariable convolution model with new information priority accumulation is proposed by Wu and Zhang [9]. Conformable fractional accumulation is proposed by Ma et al. [10]. Nonhomogenous discrete grey model with fractional-order accumulation is proposed by Wu et al. [11]. The error feedback is used to optimize the fractional accumulation grey model to predict per capita power generation by Yang and Xue [12]. Original condition, power exponent, and background value are optimized by Lu et al. The optimized nonlinear grey Bernoulli model is applied for traffic flow prediction [13]. Luo et al. used particle swarm optimization to optimize parameters and proposed a discrete grey polynomial model DGPM (1, 1, N) with new information priority accumulation, which effectively improved the prediction accuracy [14]. Xiong et al. proposed a new GM (1, 1) model to predict China’s energy consumption and production based on the optimized original condition [15]. In order to enhance the forecasting accuracy of the traditional grey model, the grey model with the fractional Hausdorff derivative is proposed by Chen et al. [16]. All the above scholars have improved and optimized the grey model, but they are all a single model, and there is no hybrid prediction model.The hybrid prediction method was first proposed by Bates and Granger [17]. Due to its applicability and high accuracy, it has received much attention. In predictive competitions, M-series competitions are very authoritative. The M4 competition is based on the previous three sessions of M1, M2, and M3. An important conclusion drawn from the M4 competition is that hybrid prediction can improve the calculation accuracy [18]. A hybrid approach for obtaining optimal cutting conditions during turning of commercially pure titanium (CP-TI) grade 2 was proposed by Khan and Maity [19]. This approach can generate adequate dimensional accuracy. The purpose of a hybrid model is to make up for the shortcomings of individual errors. For example, Zhu et al. believed that the existing prediction methods could not fully reflect the information of the pollution index sequence, so they proposed two hybrid models and conducted an empirical study taking Xingtai City as an example [20]. Gulia et al. combined the deterministic-based models with suitable statistical distribution models to predict the entire range of pollutant concentration distribution [21]. A hybrid model was established by Lai et al. from the perspective of air quality pollution sources. As emissions from pollution sources decrease, PM2.5 would also decrease [22]. Smyl proposed the combination of exponential smoothing and neural networks to predict time series [23]. Sharma et al. combined online sequential extreme learning machines with empirical pattern decomposition algorithms to predict hourly air quality [24]. Wu and Lin considered air pollution factors, such as PM2.5, PM10, SO2, CO, NO2, and O3. Wavelet decomposition, variational mode decomposition improved by sample entropy, and long-term memory neural network were used to predict air quality index [25]. The above scholars used various hybrid models to evaluate and predict air quality but did not involve the field of grey system. This paper introduces the grey accumulation operator into the statistical model. Then, the grey double exponential smooth hybrid forecasting model is proposed. The conformable fractional accumulation proposed by Ma and the adjacent accumulation will be introduced to grey double exponential smoothing, respectively, in this paper. Finally, the two grey double exponential smoothing prediction results are hybridized for prediction using the coefficient of variability. The validity of the model was verified by applying the hybrid model for the air quality prediction in Chongqing and Foshan.This article is divided into four parts: Section2 introduces the conformable fractional accumulation grey double exponential smoothing (CFES), adjacent accumulation grey double exponential smoothing (AAES), and hybrid model. Section 3 uses Nanchong as an example to predict AQI through CFES and AAES, respectively. Then, CFES and AAES are compared with the traditional grey double exponential smoothing (ES). In order to obtain better prediction results, the grey double exponential smoothing model of the two accumulation methods is used to make a hybrid prediction with the variation coefficient. Finally, this paper compares the hybrid model with ES. Section 4 is a summary of the full paper.


## 2. Grey Double Exponential Smoothing Model

CFES and AAES differ from ES in that grey exponential smoothing contains accumulation operators. The following model will explain in detail.Assume thatX0=x01,x02,x03,…,x0n−1,x0n is a nonnegative time series. xrk=∑i=1kx0i/ir−ri=1,2,…,k is the conformable fractional accumulation generating operator of order r, r∈0,1. In special cases, r takes 1 as the traditional grey first-order accumulation, and the accumulation operator is x1k=∑i=1kx0i [26]. So Xr=xr1,xr2,xr3,…,xrn−1,xrn can be deduced that the r-order accumulation sequence of X0 is(1)xr1,xr2,xr3,⋯,xrk−1,xrk=x01,x02,x03,⋯,x0k−1,x0k111⋯11012r−r12r−r⋯12r−r12r−r0013r−r⋯13r−r13r−r⋮⋮⋮⋮⋮000⋯1k−1r−r1k−1r−r000⋯01kr−r.Equations (2)(5) are derived from [27, 28].CFES based on theXr sequence is(2)S1k=αxrk+1−αS1k−1,(3)S2k=αS1k+1−αS2k−1.S1k is the grey single exponential smoothing value of the k period, and S2k is the grey double exponential smoothing value of the k period.(4)ak=2S1k−S2k,(5)bk=α1−αS1k−S2k.The fitting formula is(6)x^rk=ak+bk.The prediction formula is(7)x^rk=ak+mbk,where m is the out-of-sample size.Finally, we transform thex^rk back to the x^0k:(8)x^0k=k1−rx^rk−x^rk−1.The original sequence ofXr is x^0i=k1−rx^ri−x^ri−1i=2,3,…,k.Furthermore, the original sequence can be also expressed in the form of a matrix as follows:(9)x^01,x^02,x^03,…,x^0k−1,x^0k=x^r1,x^r2,x^r3,…,x^rk−1,x^rk1−21−r⋯00021−r⋯00⋮⋮⋮⋮0000⋯⋯−k−11−rk−11−r0−k1−r00⋯0k1−r=a1+b1,a2+b2,…,ak+bk1−21−r⋯00021−r⋯00⋮⋮⋮⋮0000⋯⋯−k−11−rk−11−r0−k1−r00⋯0k1−r.The prediction reduction formula is(10)x^01,x^02,…,x^0k,…x^0k+m=a1+b1,a2+b2,…,ak+bk,ak+2bk,…,ak+mbk1−21−r⋯0021−r⋯0⋮⋮⋮00⋯−k+m1−r00⋯k+m1−r.The above is the calculation method of the CFES. The difference between the AAES and the CFES is the accumulated method.Assume thatX0=x01,x02,x03,…,x0n−1,x0n is a nonnegative time series:(11)x11=x01,x1i=λx0i−1+x0i,i=2,3,…,n,where λ∈−1,1 is a parameter accumulated nearby, and its purpose is to adjust the weight between old and new information. After accumulation, we can obtain X1=x11,x12,x13,…,x1n−1,x1n, and the accumulation sequence of X0 can be expressed as follows:(12)x11,x12,x13,…,x1n−1,x1n=x01,x02,x03,…,x0n−1,x0n1λ0⋯0001λ⋯00001⋯00⋮⋮⋮⋮⋮000⋯λ0000⋯1λ000⋯01.Then, we useX1 for traditional double exponential smoothing. Finally, equation (13) is used to transform the prediction value to original sequence:(13)x^0k=x^1k−λx^0k−1.So the original sequence ofX0 can be also expressed as follows:(14)x^01,x^02,x^03,…,x^0k=x^r1,x^r2,x^r3,…,x^rk−1,x^rk1−λλ2…−1n−1λn−101−λ…−1n−1λn−2001…−1n−1λn−3⋮⋮⋮⋮000…λ2000…−λ000…1=a1+b1,a2+b2,⋯,an+bn1−λλ2…−1n−1λn−101−λ…−1n−1λn−2001…−1n−1λn−3⋮⋮⋮⋮000…λ2000…−λ000…1.The prediction values can be shown in the following formula.When(15)k=1,x^01=x^11,k=2,x^02=x^12−λx^01=x^12−λx^11,k=3,x^03=x^13−λx^02=x^13−λx^12−λx^11=x^13−λx^12+λ2x^11.Without loss of generality,(16)k=4,x^04=x^14−λx^02=x^14−λx^13−λx^12+λ2x^11=x^14−λx^13+−λ2x^12−λ3x^11.Ifn is an even number, then(17)x^0n−1=x^1n−1−λx^0n−2=x^1n−1−λx^0n−2+⋯+λn−3x^12−λn−2x^11,x^0n=x^1n−λx^1n−1+⋯+λn−2x^12−λn−1x^11.Ifn is an odd number, then(18)x^0n−1=x^1n−1−λx^1n−2=x^1n−1−λx^1n−2+⋯+λn−3x^12−λn−2x^11,x^0n=x^1n−λx^1n−1+⋯−λn−2x^12−λn−1x^11.When the prediction of a single model is completed, CFES and AAES are hybridized with the variation coefficient to determine the weight. The hybrid formula is shown in equations (19)(24). Suppose there are m single prediction models, and each original sequence has n data; then, the matrix formed by the sequence values xij0 of a single prediction model is Xij0:(19)Xij0=x110⋯x1n0⋮⋱⋮xm10⋯xmn0,i=1,2,…,m;j=1,2,…,n.The mean of the individual sequence prediction is(20)x¯i0=1n∑i=1nxij0.The standard deviation of the single prediction series is(21)Si=1n∑j=1nxij0−x¯i02.The variation coefficient of the single method is(22)Vi=Six¯i0,i=1,2,…,m.We can obtain the weight of the hybrid model is(23)wi=Vi∑i=1mVi.We can hybridize the individual forecasting sequences from the obtained weights:(24)Y=∑i=1mwiyj.In this paper, CFES and AAES are used to compare with ES, respectively. Then, the two models are hybridized with the variation coefficient to compare with ES again to prove that the cumulative grey double exponential smoothing has a better prediction effect. This article uses MAPE as the prediction criterion, as shown in following equation:(25)MAPE=1n∑k=1nx0k−x^0kx0k.


## 3. Application of Grey Double Exponential Smoothing in Air Quality Prediction

This part will use data to verify the validity of CFES, AAES, and hybrid model, respectively. The data in this article come fromhttps://www.aqistudy.cn/historydata/.


### 3.1. Application of CFES in Air Quality Prediction

The whale algorithm is proposed by Mirjalili and Lewis [29]. The algorithm mimics the social behavior of humpback whales, and it is inspired by the bubble-net hunting strategy. The whale algorithm is used to search for the smallest fitting error in this paper. The original sequence is X0=51,78,90,83,72,80,84,94,133,109. The conformable fractional order is r=0.5563, and the smoothing coefficient α=0.5032. Then, X0.5563=51,108.35,163.63,208.50,243.75,279.88,315.31,352.67,402.85,442.09. For convenience of calculation, the first value of the original sequence is used as the original value of the single exponential smoothing value and the double exponential smoothing value, respectively. Equations (2) and (3) are used to calculate S1k and S2k. The results are shown as follows:(26)S1k=51,79.86,122.02,165.54,204.90,242.64,279.21,316.18,359.79,401.21,S2k=51,65.52,93.95,129.98,167.68,205.40,242.55,279.60,319.96,360.85.We can obtainak and bk by equations (4) and (5):(27)ak=51,94.20,150.08,201.10,242.12,279.87,315.87,352.76,399.63,441.57,bk=0,14.52,28.43,36.03,37.70,37.72,37.14,37.06,40.36,40.89.Soan=441.57andbn=40.89x0.5563k=51,108.72,178.51,237.12,279.82,317.59,353.01,389.81,439.99,482.46,523.35,564.24,605.13 The data of December 1–10, 2019 were used to predict the AQI concentration value of December 11–13. The final prediction is(28)x0k=51,78.51,113.62,108.42,87.20,83.62,84.00,92.58,133.00,117.97,118.48,123.14,127.60.The prediction results between ES and CFES are shown in Table1. In Table 1, in order to facilitate the comparison of the prediction accuracy of the two models, the smoothing coefficient between ES and CFES takes the same value α=0.5032 . As can be seen in Table 1, the prediction error of ES is 18.82%, and the prediction error of CFES is 8.08%. It can be seen that CFES makes the prediction error significantly lower.Table 1
Comparison of CFES and ES prediction results.
DateActual valueESCFES11-Dec102128.95118.4812-Dec118135.39123.1413-Dec123141.83127.60MAPE18.82%8.08%


### 3.2. Application of AAES in Air Quality Prediction

The original sequence isX0=51,78,90,83,72,80,84,94,133,109. The whale algorithm is used to search for the smallest fitting error. The adjacent accumulation is λ=−0.3863, and the smoothing coefficient α=0.4995. The sequence X1=51,58.30,59.87,48.23,39.94,52.19,53.10,61.55,96.69,57.62 is obtained by adjacent accumulation. We can get S1k and S2k by equations (2) and (3):(29)S1=51,54.65,57.25,52.75,46.35,49.26,51.18,56.36,76.50,67.07,S2=51,52.82,55.04,53.89,50.12,49.70,50.44,53.39,64.94,66.00,where ak and bk can be obtained by equations (4) and (5):(30)ak=51,56.47,59.47,51.60,42.57,48.83,51.92,59.32,88.07,68.14,bk=0,1.82,2.21,−1.14,−3.77,−0.43,0.74,2.96,11.54,1.07.Soan=68.14 and bn=1.07.(31)x1k=51,58.29,61.69,50.46,38.81,48.40,52.66,62.28,99.61,69.21,70.27,71.34,72.41.The final prediction is(32)x^0k=51,77.99,91.82,85.93,72.00,76.22,82.10,94.00,135.92,121.71,117.29,116.65,117.47.As can be seen in Figure1, the fitted value and predicted value of the AAES are closer to the actual value. The fitting error and forecasting error of AAES are 2.64% and 6.88%, respectively. The fitting error and forecasting error of ES are 4.94% and 18.93%. It can be seen that AAES has a higher prediction accuracy than ES.Figure 1
Comparison of AAES and ES forecasting results.Forecasting AQI in Nanchong demonstrated that the grey double exponential smoothing with conformable fractional accumulation and adjacent accumulation can improve the prediction effect. It is proved that the grey double exponential smoothing with accumulation mode is meaningful. In order to obtain a more accurate prediction effect, a hybrid prediction model is proposed.


### 3.3. Application of Hybrid Model in Air Quality Prediction

Case 1.
Forecasting AQI, PM2.5, and PM10 in Chongqing.
Compared with ES, the grey exponential smoothing with the accumulation operator has improved the prediction accuracy to a certain extent, but individual prediction values sometimes have large errors. The conformable fractional accumulation is the principle of old information first, and adjacent accumulation is the principle of new information first. In order to make the model more applicable, a hybrid model is used to balance new information with old information to make up for the shortcomings of a single model.
Table2 shows the original data of Chongqing AQI, PM2.5, and PM10 in December 2019. Chongqing is foggy and is known as “fog Chongqing.” The annual average foggy day in Chongqing is 104 days. The average annual foggy day in London, England, which is known as the world ’s foggy city, is only 94 days. Therefore, it is necessary to study Chongqing air quality.Table 2
The original data of Chongqing AQI, PM2.5, and PM10 in December 2019.
DateAQIPM2.5PM101-Dec4632462-Dec6746753-Dec7555864-Dec94701055-Dec8462886-Dec7555847-Dec7756848-Dec107801219-Dec1299813810-Dec1279613711-Dec13510314212-Dec14310914613-Dec145111156This section uses AQI as an example to elaborate the model. In the ES model, the fitting error is minimal whenα=0.4.We can obtain(33)S1=46,54.40,62.64,75.18,78.71,77.23,77.14,89.08,105.05,113.83,S2=46,49.36,54.67,62.88,69.21,72.42,74.30,80.22,90.15,99.62,by(34)S1k=0.4×X0k+0.6×S1k−1,S2k=0.4×S1k+0.6×S2k−1.We can obtain(35)ak=46,59.44,70.61,87.49,88.21,82.04,79.97,97.95,119.95,128.04,bk=0,3.36,5.31,8.20,6.33,3.21,1.89,5.91,9.93,9.47,by(36)ak=2S1k−S2k,bk=0.67×S1−S2.Soan=128.04 and bn=9.47.The forecasting result of ES can be obtained as follows:(37)x^0k=46,62.80,75.92,95.70,94.54,85.24,81.85,103.86,129.88,137.51,146.98,156.45,165.93.When the conformable fractional grey double exponential smoothing takesr=0.5385andα=0.5,(38)x^0.5385k=46,67,95.20,115.41,103.66,84.44,77.00,103.36,131.79,135.14,124.60,129.70,134.58.When the adjacent accumulation grey double exponential smoothing takesλ=−0.4540andα=0.5122,(39)x^1k=46,67,74.99,94.00,87.20,75.04,73.96,105.14,135.54,138.98,144.96,152.08,159.72,where x^0.5385k obtained by CFES and x^1k obtained by AAES are combined using the variation coefficient. A matrix Xij0 composed of the sequence values xij0 of a single prediction model obtained by two accumulation methods is(40)Xij0=466774.9994.0087.2075.0473.96105.14135.54138.98144.96152.08159.72466795.20115.41103.6684.4477.00103.36131.79135.14124.60129.70134.58.We can getx¯ij0=104.20103.68,Si=32.8026.15, Vi=0.310.25, and wi=0.560.44 by equations (11)(14).So the final fitting values and forecasting values can be obtained bywi:(41)Y=466783.98103.5394.5279.2275.31104.35133.87137.27135.90142.12148.83.As can be seen from Figure2, the fitting error of CFES is large and the deviation from the true value is higher, but the fitting of the AAES model to the actual value is higher. After combining the two models of CFES and AAES, the curve fitting effect of the hybrid model is significantly better than the CFES model. Judging from the prediction results, the predicted value polyline of the combined model is closer to the actual value polyline, so it can explain the superiority of the hybrid model compared to the single model.Figure 2
Prediction results of AQI concentration in four models.Table3 is the PM2.5 concentration values predicted by ES, CFES, AAES, and hybrid model. The smoothing coefficients and parameters of CFES and AAES are obtained by using the whale algorithm in Matlab R2014b. The prediction error obtained by the hybrid model in Table 3 is 0.97%, which is far lower than the traditional double exponential smoothing prediction error of 10.98%, so the hybrid model has better prediction performance.Table 3
The forecasting results of PM2.5 concentration.
DateCFESAAESHybrid modelESParameterα=0.5,r=0.5216α=0.5041,λ=−0.4905wc=0.45,wA=0.55α=0.41-Dec323232322-Dec46.0045.9945.9943.203-Dec68.9654.5861.0354.884-Dec85.7870.0077.0771.105-Dec77.0164.7870.2770.376-Dec62.2655.2758.4063.167-Dec56.0053.6454.7060.068-Dec77.0677.9377.5470.509-Dec100.27102.52101.5198.5710-Dec102.71105.38104.18104.47MAPE9.83%2.69%5.81%5.60%11-Dec93.85109.91102.71112.0112-Dec97.84115.27107.45119.5613-Dec101.66121.03112.34127.10MAPE9.18%7.16%0.97%10.98%Table4 shows the predicted concentrations of PM10 under the four models. In PM2.5 concentration prediction, the single prediction error of CFES is relatively large. In PM10 concentration prediction, the single prediction error of AAES is relatively large. The hybrid model makes up for the lack of single errors, integrates ES new information and old information first, and improves prediction accuracy. Comparing the hybrid model with ES, the prediction error of ES is significantly higher than the prediction error of the hybrid model, so the hybrid model has more advantages.Table 4
The forecasting results of PM10 concentration.
DateCFESAAESHybrid modelESParameterα=0.5000,r=0.5488α=0.4588,λ=−0.5153wc=0.46,wA=0.54α=0.51-Dec464646462-Dec75.0074.5674.7675.003-Dec108.5187.1896.9093.254-Dec129.48105.00116.15115.005-Dec109.9493.05100.74100.946-Dec92.8484.0088.0390.197-Dec84.0080.9382.3385.958-Dec116.88114.41115.53121.419-Dec142.03141.09141.52147.1710-Dec145.09147.28146.28150.32MAPE9.72%2.65%5.87%5.90%11-Dec135.48154.40145.78161.0912-Dec140.90162.00152.39171.8713-Dec146.08169.84159.02182.64MAPE4.82%9.52%2.99%16.08%The above content uses Chongqing as an example to explain the grey double exponential smoothing. The following content further verifies the model with Foshan, and the calculation process is the same, so it will not be repeated.Case 2.
Forecasting AQI, PM2.5, and PM10 in Foshan
Foshan is located in the central and southern part of Guangdong Province, China. It is an important part of the Pearl River Delta economic circle. While vigorously developing the economy, the country is also actively addressing air pollution. Compared with the air quality in December 2013, air quality has improved significantly in 2019. With the increasing pursuit in a better life, a high-quality atmospheric environment has received more attention. This section uses Foshan as an example to predict the three air quality indicators of AQI, PM2.5, and PM10.
Table5 shows the prediction results of the AQI in Foshan City. The optimal parameters of CFES α=0.5174 and r=0.4383 and the optimal parameters of AAES α=0.4605 and λ=0.1814 are obtained from the whale algorithm search. The two sets of prediction results obtained based on the optimal parameters were used to make a fixed weight combination prediction using the variation coefficient. The CFES weight was 0.3943, and the AAES weight was 0.6057. It can be seen from Table 3 that the prediction error of the single model has exceeded 10% and is no longer suitable for single prediction, but the prediction error of CFES is only more than 0.03% of the MAPE standard, so a hybrid prediction can be considered. The hybrid fitting error is 4.38%, and the prediction error is 8.99%. Therefore, compared with the single model, the hybrid model can improve the prediction accuracy.Table 5
The forecasting results of AQI concentration.
DateActual valueCFESAAESHybrid modelParameter—α=0.5174,r=0.4383α=0.4605,λ=0.1814wc=0.3943,wA=0.60571-Dec636363632-Dec4445.5344.6044.973-Dec5670.7553.6660.404-Dec5874.1458.0064.365-Dec5464.6055.6959.216-Dec5861.5358.3359.597-Dec5757.0558.1157.698-Dec6764.7167.0266.119-Dec100100.00100.00100.0010-Dec108115.69116.97116.46MAPE—9.40%1.95%4.38%11-Dec11197.46126.93115.3112-Dec112102.35138.16124.0413-Dec118107.05149.16132.56MAPE—10.03%21.37%8.99%When using ES to predict AQI, the smoothing coefficient is determined to be 0.5 based on the minimum fitting error. The comparison between ES and the hybrid prediction results is shown in Figure3. The fitting error of ES is 2.66%, the fitting error of the hybrid model is 4.38%, and the fitting error of the two models is far below 10%, while the prediction error of the two models is quite different, and the fitting error of ES is 25.92%, more than 10%, so it can be said that the combination model is meaningful.Figure 3
Compared forecasting results between hybrid model and ES.It can be seen in Table6 that the prediction error of CFES is 19.05%, the prediction error of AAES is 12.34%, the hybrid prediction error is 4.04%, and the prediction accuracy has been greatly improved. Compare the hybrid prediction result with the ES. As shown in Figure 4, the predicted polyline of the hybrid model is closer to the actual value of PM2.5.Table 6
The forecasting results of PM2.5 concentration.
DateActual valueCFESAAESHybrid modelParameter—α=0.5029,r=0.3760α=0.4573,λ=0.4571wc=0.4057,wA=0.59431-Dec353535352-Dec1919.1119.0019.043-Dec2935.2128.7731.384-Dec3240.6732.0035.525-Dec2632.9728.5830.366-Dec3032.0330.1830.937-Dec3434.3334.5134.448-Dec3333.0034.6533.989-Dec5149.9450.8650.4910-Dec6669.0369.3669.23MAPE—9.03%2.32%4.95%11-Dec6653.8070.8663.9412-Dec6756.8080.1370.6713-Dec7859.7185.8575.25MAPE—19.05%12.34%4.04%Figure 4
Compared forecasting results between hybrid model and ES in PM2.5.Table7 is similar to Table 6, and Figure 5 is similar to Figure 4. The hybrid model is superior to the single model. In general, when the error of the single model is greater than or close to 10%, to obtain a higher fitting accuracy, a hybrid model can be used. In the above charts, the superiority of the hybrid model can be verified. The new information-first principle of adjacent accumulation and the old information-first principle of conformable fractional-order accumulation are fused in this paper to obtain a higher prediction accuracy. Therefore, the proposed hybrid model has certain practical significance.Table 7
The forecasting results of PM10 concentration.
DateActual valueCFESAAESHybrid modelParameter—α=0.5000,r=0.5837α=0.5037,λ=0.2386wc=0.3632,wA=0.63681-Dec757575752-Dec4444.0043.9043.943-Dec6174.0257.7763.674-Dec6683.1966.0072.245-Dec5872.0860.2264.526-Dec6571.2165.2267.407-Dec6366.2264.4065.068-Dec6059.9960.8160.519-Dec9794.3697.0096.0410-Dec120126.12129.08128.01MAPE—9.42%2.08%4.07%11-Dec12295.41138.30122.7312-Dec12698.93152.98133.3513-Dec154102.28166.35143.08MAPE—25.62%14.26%4.51%Figure 5
Compared forecasting results between hybrid model and ES in PM10.


## 3.1. Application of CFES in Air Quality Prediction

The whale algorithm is proposed by Mirjalili and Lewis [29]. The algorithm mimics the social behavior of humpback whales, and it is inspired by the bubble-net hunting strategy. The whale algorithm is used to search for the smallest fitting error in this paper. The original sequence is X0=51,78,90,83,72,80,84,94,133,109. The conformable fractional order is r=0.5563, and the smoothing coefficient α=0.5032. Then, X0.5563=51,108.35,163.63,208.50,243.75,279.88,315.31,352.67,402.85,442.09. For convenience of calculation, the first value of the original sequence is used as the original value of the single exponential smoothing value and the double exponential smoothing value, respectively. Equations (2) and (3) are used to calculate S1k and S2k. The results are shown as follows:(26)S1k=51,79.86,122.02,165.54,204.90,242.64,279.21,316.18,359.79,401.21,S2k=51,65.52,93.95,129.98,167.68,205.40,242.55,279.60,319.96,360.85.We can obtainak and bk by equations (4) and (5):(27)ak=51,94.20,150.08,201.10,242.12,279.87,315.87,352.76,399.63,441.57,bk=0,14.52,28.43,36.03,37.70,37.72,37.14,37.06,40.36,40.89.Soan=441.57andbn=40.89x0.5563k=51,108.72,178.51,237.12,279.82,317.59,353.01,389.81,439.99,482.46,523.35,564.24,605.13 The data of December 1–10, 2019 were used to predict the AQI concentration value of December 11–13. The final prediction is(28)x0k=51,78.51,113.62,108.42,87.20,83.62,84.00,92.58,133.00,117.97,118.48,123.14,127.60.The prediction results between ES and CFES are shown in Table1. In Table 1, in order to facilitate the comparison of the prediction accuracy of the two models, the smoothing coefficient between ES and CFES takes the same value α=0.5032 . As can be seen in Table 1, the prediction error of ES is 18.82%, and the prediction error of CFES is 8.08%. It can be seen that CFES makes the prediction error significantly lower.Table 1
Comparison of CFES and ES prediction results.
DateActual valueESCFES11-Dec102128.95118.4812-Dec118135.39123.1413-Dec123141.83127.60MAPE18.82%8.08%


## 3.2. Application of AAES in Air Quality Prediction

The original sequence isX0=51,78,90,83,72,80,84,94,133,109. The whale algorithm is used to search for the smallest fitting error. The adjacent accumulation is λ=−0.3863, and the smoothing coefficient α=0.4995. The sequence X1=51,58.30,59.87,48.23,39.94,52.19,53.10,61.55,96.69,57.62 is obtained by adjacent accumulation. We can get S1k and S2k by equations (2) and (3):(29)S1=51,54.65,57.25,52.75,46.35,49.26,51.18,56.36,76.50,67.07,S2=51,52.82,55.04,53.89,50.12,49.70,50.44,53.39,64.94,66.00,where ak and bk can be obtained by equations (4) and (5):(30)ak=51,56.47,59.47,51.60,42.57,48.83,51.92,59.32,88.07,68.14,bk=0,1.82,2.21,−1.14,−3.77,−0.43,0.74,2.96,11.54,1.07.Soan=68.14 and bn=1.07.(31)x1k=51,58.29,61.69,50.46,38.81,48.40,52.66,62.28,99.61,69.21,70.27,71.34,72.41.The final prediction is(32)x^0k=51,77.99,91.82,85.93,72.00,76.22,82.10,94.00,135.92,121.71,117.29,116.65,117.47.As can be seen in Figure1, the fitted value and predicted value of the AAES are closer to the actual value. The fitting error and forecasting error of AAES are 2.64% and 6.88%, respectively. The fitting error and forecasting error of ES are 4.94% and 18.93%. It can be seen that AAES has a higher prediction accuracy than ES.Figure 1
Comparison of AAES and ES forecasting results.Forecasting AQI in Nanchong demonstrated that the grey double exponential smoothing with conformable fractional accumulation and adjacent accumulation can improve the prediction effect. It is proved that the grey double exponential smoothing with accumulation mode is meaningful. In order to obtain a more accurate prediction effect, a hybrid prediction model is proposed.


## 3.3. Application of Hybrid Model in Air Quality Prediction

Case 1.
Forecasting AQI, PM2.5, and PM10 in Chongqing.
Compared with ES, the grey exponential smoothing with the accumulation operator has improved the prediction accuracy to a certain extent, but individual prediction values sometimes have large errors. The conformable fractional accumulation is the principle of old information first, and adjacent accumulation is the principle of new information first. In order to make the model more applicable, a hybrid model is used to balance new information with old information to make up for the shortcomings of a single model.
Table2 shows the original data of Chongqing AQI, PM2.5, and PM10 in December 2019. Chongqing is foggy and is known as “fog Chongqing.” The annual average foggy day in Chongqing is 104 days. The average annual foggy day in London, England, which is known as the world ’s foggy city, is only 94 days. Therefore, it is necessary to study Chongqing air quality.Table 2
The original data of Chongqing AQI, PM2.5, and PM10 in December 2019.
DateAQIPM2.5PM101-Dec4632462-Dec6746753-Dec7555864-Dec94701055-Dec8462886-Dec7555847-Dec7756848-Dec107801219-Dec1299813810-Dec1279613711-Dec13510314212-Dec14310914613-Dec145111156This section uses AQI as an example to elaborate the model. In the ES model, the fitting error is minimal whenα=0.4.We can obtain(33)S1=46,54.40,62.64,75.18,78.71,77.23,77.14,89.08,105.05,113.83,S2=46,49.36,54.67,62.88,69.21,72.42,74.30,80.22,90.15,99.62,by(34)S1k=0.4×X0k+0.6×S1k−1,S2k=0.4×S1k+0.6×S2k−1.We can obtain(35)ak=46,59.44,70.61,87.49,88.21,82.04,79.97,97.95,119.95,128.04,bk=0,3.36,5.31,8.20,6.33,3.21,1.89,5.91,9.93,9.47,by(36)ak=2S1k−S2k,bk=0.67×S1−S2.Soan=128.04 and bn=9.47.The forecasting result of ES can be obtained as follows:(37)x^0k=46,62.80,75.92,95.70,94.54,85.24,81.85,103.86,129.88,137.51,146.98,156.45,165.93.When the conformable fractional grey double exponential smoothing takesr=0.5385andα=0.5,(38)x^0.5385k=46,67,95.20,115.41,103.66,84.44,77.00,103.36,131.79,135.14,124.60,129.70,134.58.When the adjacent accumulation grey double exponential smoothing takesλ=−0.4540andα=0.5122,(39)x^1k=46,67,74.99,94.00,87.20,75.04,73.96,105.14,135.54,138.98,144.96,152.08,159.72,where x^0.5385k obtained by CFES and x^1k obtained by AAES are combined using the variation coefficient. A matrix Xij0 composed of the sequence values xij0 of a single prediction model obtained by two accumulation methods is(40)Xij0=466774.9994.0087.2075.0473.96105.14135.54138.98144.96152.08159.72466795.20115.41103.6684.4477.00103.36131.79135.14124.60129.70134.58.We can getx¯ij0=104.20103.68,Si=32.8026.15, Vi=0.310.25, and wi=0.560.44 by equations (11)(14).So the final fitting values and forecasting values can be obtained bywi:(41)Y=466783.98103.5394.5279.2275.31104.35133.87137.27135.90142.12148.83.As can be seen from Figure2, the fitting error of CFES is large and the deviation from the true value is higher, but the fitting of the AAES model to the actual value is higher. After combining the two models of CFES and AAES, the curve fitting effect of the hybrid model is significantly better than the CFES model. Judging from the prediction results, the predicted value polyline of the combined model is closer to the actual value polyline, so it can explain the superiority of the hybrid model compared to the single model.Figure 2
Prediction results of AQI concentration in four models.Table3 is the PM2.5 concentration values predicted by ES, CFES, AAES, and hybrid model. The smoothing coefficients and parameters of CFES and AAES are obtained by using the whale algorithm in Matlab R2014b. The prediction error obtained by the hybrid model in Table 3 is 0.97%, which is far lower than the traditional double exponential smoothing prediction error of 10.98%, so the hybrid model has better prediction performance.Table 3
The forecasting results of PM2.5 concentration.
DateCFESAAESHybrid modelESParameterα=0.5,r=0.5216α=0.5041,λ=−0.4905wc=0.45,wA=0.55α=0.41-Dec323232322-Dec46.0045.9945.9943.203-Dec68.9654.5861.0354.884-Dec85.7870.0077.0771.105-Dec77.0164.7870.2770.376-Dec62.2655.2758.4063.167-Dec56.0053.6454.7060.068-Dec77.0677.9377.5470.509-Dec100.27102.52101.5198.5710-Dec102.71105.38104.18104.47MAPE9.83%2.69%5.81%5.60%11-Dec93.85109.91102.71112.0112-Dec97.84115.27107.45119.5613-Dec101.66121.03112.34127.10MAPE9.18%7.16%0.97%10.98%Table4 shows the predicted concentrations of PM10 under the four models. In PM2.5 concentration prediction, the single prediction error of CFES is relatively large. In PM10 concentration prediction, the single prediction error of AAES is relatively large. The hybrid model makes up for the lack of single errors, integrates ES new information and old information first, and improves prediction accuracy. Comparing the hybrid model with ES, the prediction error of ES is significantly higher than the prediction error of the hybrid model, so the hybrid model has more advantages.Table 4
The forecasting results of PM10 concentration.
DateCFESAAESHybrid modelESParameterα=0.5000,r=0.5488α=0.4588,λ=−0.5153wc=0.46,wA=0.54α=0.51-Dec464646462-Dec75.0074.5674.7675.003-Dec108.5187.1896.9093.254-Dec129.48105.00116.15115.005-Dec109.9493.05100.74100.946-Dec92.8484.0088.0390.197-Dec84.0080.9382.3385.958-Dec116.88114.41115.53121.419-Dec142.03141.09141.52147.1710-Dec145.09147.28146.28150.32MAPE9.72%2.65%5.87%5.90%11-Dec135.48154.40145.78161.0912-Dec140.90162.00152.39171.8713-Dec146.08169.84159.02182.64MAPE4.82%9.52%2.99%16.08%The above content uses Chongqing as an example to explain the grey double exponential smoothing. The following content further verifies the model with Foshan, and the calculation process is the same, so it will not be repeated.Case 2.
Forecasting AQI, PM2.5, and PM10 in Foshan
Foshan is located in the central and southern part of Guangdong Province, China. It is an important part of the Pearl River Delta economic circle. While vigorously developing the economy, the country is also actively addressing air pollution. Compared with the air quality in December 2013, air quality has improved significantly in 2019. With the increasing pursuit in a better life, a high-quality atmospheric environment has received more attention. This section uses Foshan as an example to predict the three air quality indicators of AQI, PM2.5, and PM10.
Table5 shows the prediction results of the AQI in Foshan City. The optimal parameters of CFES α=0.5174 and r=0.4383 and the optimal parameters of AAES α=0.4605 and λ=0.1814 are obtained from the whale algorithm search. The two sets of prediction results obtained based on the optimal parameters were used to make a fixed weight combination prediction using the variation coefficient. The CFES weight was 0.3943, and the AAES weight was 0.6057. It can be seen from Table 3 that the prediction error of the single model has exceeded 10% and is no longer suitable for single prediction, but the prediction error of CFES is only more than 0.03% of the MAPE standard, so a hybrid prediction can be considered. The hybrid fitting error is 4.38%, and the prediction error is 8.99%. Therefore, compared with the single model, the hybrid model can improve the prediction accuracy.Table 5
The forecasting results of AQI concentration.
DateActual valueCFESAAESHybrid modelParameter—α=0.5174,r=0.4383α=0.4605,λ=0.1814wc=0.3943,wA=0.60571-Dec636363632-Dec4445.5344.6044.973-Dec5670.7553.6660.404-Dec5874.1458.0064.365-Dec5464.6055.6959.216-Dec5861.5358.3359.597-Dec5757.0558.1157.698-Dec6764.7167.0266.119-Dec100100.00100.00100.0010-Dec108115.69116.97116.46MAPE—9.40%1.95%4.38%11-Dec11197.46126.93115.3112-Dec112102.35138.16124.0413-Dec118107.05149.16132.56MAPE—10.03%21.37%8.99%When using ES to predict AQI, the smoothing coefficient is determined to be 0.5 based on the minimum fitting error. The comparison between ES and the hybrid prediction results is shown in Figure3. The fitting error of ES is 2.66%, the fitting error of the hybrid model is 4.38%, and the fitting error of the two models is far below 10%, while the prediction error of the two models is quite different, and the fitting error of ES is 25.92%, more than 10%, so it can be said that the combination model is meaningful.Figure 3
Compared forecasting results between hybrid model and ES.It can be seen in Table6 that the prediction error of CFES is 19.05%, the prediction error of AAES is 12.34%, the hybrid prediction error is 4.04%, and the prediction accuracy has been greatly improved. Compare the hybrid prediction result with the ES. As shown in Figure 4, the predicted polyline of the hybrid model is closer to the actual value of PM2.5.Table 6
The forecasting results of PM2.5 concentration.
DateActual valueCFESAAESHybrid modelParameter—α=0.5029,r=0.3760α=0.4573,λ=0.4571wc=0.4057,wA=0.59431-Dec353535352-Dec1919.1119.0019.043-Dec2935.2128.7731.384-Dec3240.6732.0035.525-Dec2632.9728.5830.366-Dec3032.0330.1830.937-Dec3434.3334.5134.448-Dec3333.0034.6533.989-Dec5149.9450.8650.4910-Dec6669.0369.3669.23MAPE—9.03%2.32%4.95%11-Dec6653.8070.8663.9412-Dec6756.8080.1370.6713-Dec7859.7185.8575.25MAPE—19.05%12.34%4.04%Figure 4
Compared forecasting results between hybrid model and ES in PM2.5.Table7 is similar to Table 6, and Figure 5 is similar to Figure 4. The hybrid model is superior to the single model. In general, when the error of the single model is greater than or close to 10%, to obtain a higher fitting accuracy, a hybrid model can be used. In the above charts, the superiority of the hybrid model can be verified. The new information-first principle of adjacent accumulation and the old information-first principle of conformable fractional-order accumulation are fused in this paper to obtain a higher prediction accuracy. Therefore, the proposed hybrid model has certain practical significance.Table 7
The forecasting results of PM10 concentration.
DateActual valueCFESAAESHybrid modelParameter—α=0.5000,r=0.5837α=0.5037,λ=0.2386wc=0.3632,wA=0.63681-Dec757575752-Dec4444.0043.9043.943-Dec6174.0257.7763.674-Dec6683.1966.0072.245-Dec5872.0860.2264.526-Dec6571.2165.2267.407-Dec6366.2264.4065.068-Dec6059.9960.8160.519-Dec9794.3697.0096.0410-Dec120126.12129.08128.01MAPE—9.42%2.08%4.07%11-Dec12295.41138.30122.7312-Dec12698.93152.98133.3513-Dec154102.28166.35143.08MAPE—25.62%14.26%4.51%Figure 5
Compared forecasting results between hybrid model and ES in PM10.


## 4. Conclusion

The conformable fractional accumulation proposed by Ma is the old information first principle, and the adjacent accumulation is the new information-first principle. In this paper, the conformable fractional order and adjacent are added to the double exponential smoothing model, respectively. Taking the AQI of Nanchong in December 2019 as an example, two sets of grey exponential smoothing models are used to obtain two sets of sequences. In order to facilitate the comparison of model priorities, the conformable fractional smoothing coefficient is the same as the traditional double exponential smoothing coefficient, proving that CFES is better than ES. Under the condition that adjacent cumulative grey double exponential smoothing and traditional double exponential smoothing have the same smoothing coefficient, AAES is also better than ES. Therefore, it is proved that the exponential smoothing model with accumulation is better than the ES.CFES and AAES models have their advantages and disadvantages in general. Generally speaking, the prediction accuracy of CFES is better than the fitting accuracy, while the fitting accuracy of AAES is better than its own prediction accuracy. Considering the two situations, in order to get more good results, this paper uses a combination of the variation coefficient of the CFES and AAES to obtain a fixed weighted combination of predicted values. From the case of this article, it can be seen that the hybrid prediction result is significantly better than the prediction result using the single grey double exponential smoothing model. Comparing the hybrid results with the prediction results of ES, it is also verified that the hybrid model is better than ES.In general, from the perspective of predicting air quality, it is of practical significance to put two new grey accumulation methods of conformable fractional accumulation and adjacent accumulation into the exponential smoothing model and their hybrid model, respectively. In future work, we can combine these two new accumulation methods with Holt–Winters model for further research.



---

*Source: 9427102-2020-07-26.xml*