An Improved Prediction Model Combining Inverse Exponential Smoothing and Markov Chain

. On the basis of the triple exponential smoothing prediction model, this paper introduces the reverse prediction idea and es-tablishes the reverse triple exponential smoothing model by setting parameters such as threshold value and iteration times and reasonably correcting its initial value. This method can eﬀectively reduce the error of early prediction value. At the same time, aiming at the problem that the predicting advantage of the reverse triple exponential smoothing model weakens in the later period, Markov theory is introduced to correct its error value, and an improved prediction model combining inverse exponential smoothing and Markov chain is further established. The improved model combines the advantages of index model trend prediction and Markov ﬂuctuation prediction, and the prediction accuracy and stability of the model are signiﬁcantly improved through case tests.


Introduction
National defense expenditure is an important component of national financial expenditure. It is the source of funds and important support for national defense and military construction. It reflects the economic level of a country's investment in national defense construction and embodies national defense policy and national defense strategy [1]. Predicting a country's national defense expenditure is not only helpful to analyze the trend of the country's national defense and military construction but also helpful to analyze the relationship between its national defense expenditure and economic growth. erefore, it is of far-reaching practical significance to make a reasonable prediction of national defense expenditure. In this paper, a new model is established on the basis of exponential smoothing model to effectively predict defense expenditure.
Exponential smoothing method is a time series analysis and prediction method.
is method predicts the future trend according to the current situation and data by calculating the smoothing value of the index and combining with a reasonable time series prediction model [2][3][4]. e exponential smoothing method can be divided into single exponential smoothing method, double exponential smoothing method, and triple exponential smoothing method according to exponential times. Among them, the triple exponential smoothing method is often used to fit and predict nonlinear time series and has achieved good prediction effect especially in short-term and medium-term prediction of nonlinear time series, with small error fluctuation range and strong credibility. At present, it has been widely used in the fields of public transportation passenger volume prediction [5,6], economic output value prediction [7,8], spare parts prediction [9,10], wind speed prediction of wind farms [11,12], building displacement prediction [13,14], GPS PWV prediction [15], Docker container resource load prediction [16], etc. Qi and Huo proposed a single exponential smoothing model based on self-adaptation. By introducing the approximate dynamic programming method and combining with the actual traffic flow data, the exponential smoothing coefficient was optimized to make it update automatically with the prediction process, thus ensuring the real-time accuracy of the prediction [17]. Wang et al. proposed an adaptive dynamic cubic exponential smoothing prediction method. In this method, the carpet search method is used, the best smoothing coefficient is obtained according to the principle of minimum sum of squares error, and the prediction effect of the model is verified by an example of wind speed data [18]. Mi et al. proposed a short-term power load forecasting method based on improved exponential smoothing gray model. is method combined the exponential smoothing model and gray model and used the 0.618 method to search for the best smoothing coefficient, which achieved good prediction effect [19]. Liu et al. proposed a new short-medium satellite clock error prediction algorithm based on the modified exponential smoothing method, improved the weighted parameters in ES, and proposed the dynamic weighted parameters based on the sliding window. e gray scale model (GM) is introduced to learn the prediction error of DES, which improves the prediction accuracy of the algorithm [20].
In this paper, the triple exponential smoothing predicting model is applied to the field of defense expenditure prediction. On the basis of the triple exponential smoothing model, the error trend and fluctuation of the initial data are fully considered, the reverse predicting idea and Markov state transition matrix are introduced to correct its data fluctuation, and a reverse triple exponential smoothing model based on Markov correction is established. After example verification, the new model has higher prediction accuracy in the field of national defense expenditure than the traditional triple exponential smoothing model.

Triple Exponential Smoothing Model
Exponential smoothing model is a weighted average model that uses dynamic weight coefficients to weigh the original data. And the biggest characteristic of this method is that it focuses on the influence of recent data on the prediction model [21]. In other words, the more recent the data, the greater the weight coefficient and the smaller the weight coefficient of the earlier data. e triple exponential smoothing method is to add another exponential smoothing on the basis of the first exponential smoothing and the second quadratic exponential smoothing. By estimating the parameters of the quadratic curve model, the nonlinear time series can be adjusted to eliminate irregular disturbances and random errors. It is suitable for numerical prediction of quadratic curve trend of original data.

Traditional Triple Exponential Smoothing Model. Let the time series be
where x i is the time series data at time i, x 1 is the first group of time series data, x 2 is the second group of time series data, x 3 is the third group of time series data, and x n is the nth group of time series data. Single exponential smoothing series: Double exponential smoothing series: Triple exponential smoothing series: Among them, S 1 t , S 2 t , and S 3 t are exponential smoothing values (t � 1, 2, . . . , n); S 1 0 , S 2 0 , and S 3 0 are exponential smoothing initial values, which generally take the first original value or the average of the previous original values. According to the initial value of exponential smoothing and the original time series, the exponential smoothing value at the following time is determined. e recurrence formulas are as follows: Among them, α is called smoothing coefficient (0 < α < 1). Existing literatures usually use MSE, MAE, AARE minimum principle, or artificial subjective test to determine the reasonable value of α. e larger α, the higher the emphasis on new data in the prediction, the greater the role of new data, the higher the sensitivity of prediction results, and the better the ability to adapt to new levels. e smaller α, the higher the emphasis on old data in the prediction, the more conservative the prediction results, and the slower the response to the changes of the actual data, and lag is easy to occur [22,23]. In this paper, the AARE minimum principle is used to determine the value of α. Parameters a t , b t , and c t are usually calculated by using exponential smoothing values. e parameter estimation formulas are as follows: e quadratic parabola model is established by using parameters a t , b t , and c t . T is used to represent the number of predictive lead periods to predict the future value at t + T: Generally speaking, when the number of time series of original data is large (i.e., for data with >25 items), the triple exponential smoothing method takes S 1 0 � S 2 0 � S 3 0 � x 1 . When the original data is 25 items or less, the average value of previous periods of data is often taken as the initial value of exponential smoothing. However, in general, the selection of initial values has great subjectivity, which will cause certain error influence on the prediction trend of the later model [24]. Although this error often has little influence in the medium-and long-term prediction composed of a large amount of data, it cannot be ignored in the short-term prediction.
erefore, this paper introduced the reverse prediction. Due to the unreasonable selection of initial values, the error of future predicted values may be large. erefore, a reverse cubic exponential smoothing model is established to solve this problem.

Reverse Triple Exponential Smoothing Model.
e inverse triple exponential smoothing model is based on the triple exponential smoothing model and uses the idea of reverse prediction to correct the initial value of exponential smoothing. Firstly, the traditional triple exponential smoothing model is used to predict and obtain n + 1 to n + j predicted values (j is selected according to needs; usually j ≤ n, and this paper only takes the n + 3 value as an example). And at this time, the initial value selected by the model is the first value of the original data. en, the obtained n + 1, n + 2, and n + 3 predicted values and n − 1 actual values are used for reverse prediction to obtain the first 3 reverse predicted values of the original data. According to the quadratic parabola trend of triple exponential smoothing model, if this set of data fitting effect is good, the initial value obtained by reverse prediction is closely related to the first three items of actual data. And it has similar numerical values. e modified initial value can be obtained by weighting the predicted initial value and the initial value of the original data. e specific steps are as follows.
be that original sequence of single reverse prediction: (i) Step 1. e traditional triple exponential smoothing model is used to fit the original data containing n items, and the prediction can be obtained: (ii) Step 2. Establish a reverse prediction: Step 3. Establish a triple exponential smoothing model for the newly established reverse single prediction original sequence X (0) 1 to obtain (iv) Step 4. e first three initial values are weighted and combined to obtain an initial value C 1 (k) after single reverse prediction correction: where η is the adjustment factor (usually 0.5). (v) Step 5. Set the threshold index.
According to the actual demand, the corresponding threshold value is set, and the fitting accuracy of the new model is judged by using different accuracy test indexes. If it meets the requirements, the model can be used for prediction. If it does not meet the requirements, the next reverse prediction is carried out until the initial value C n (k) meets the accuracy requirements after n iterations, and then the corresponding prediction is carried out by using the model. (vi) Step 6. Output the predicted value. e improved triple exponential smoothing model is established by using the modified initial values C n (k), and the predicted values can be obtained:

Improved Prediction Model Combining Inverse Exponential Smoothing and
Markov Chain e Markov model is a random time series analysis method, which predicts the future state of things by studying different states and state transition probability matrices of things. It has high scientific accuracy and adaptability [14,[25][26][27].
is method requires less historical data and only needs the recent data and information of the predicted object to make prediction. It has better error correction effect for data with large random fluctuation. Markov theory is a branch of stochastic process, which is a method to predict future system development according to the transition probability between states. e Markov prediction model is a random and variable mathematical process; the core of modeling is to master the law of system state transition. e basic idea of the Markov probabilistic prediction model is to analyze the current situation of the system and use Markov chain to solve the probability of a particular state to which the system may change in the future.
Assuming time series X(t), t ∈ T { } and the observable data in time series are discrete, we divide the range of error values into r intervals. So, this time series has r states E � E 1 , E 2 , . . . , E r . Suppose that the probability of the sequence X(t) in state E i is a i (t) � P x i � i . When the probability of sequence X(t) to be transferred from the first state E i to another state E j is p ij , then p ij � P(x t+i � j | x t � i); when the time series value at time t + 1 is only related to the transition probability and time series value at time t, then this time series becomes a time series with Markov property.
Based on the inverse triple exponential smoothing model, Markov's state transition probability matrix is used to predict the error fluctuation, thus correcting the error and further improving the prediction accuracy of the model. e specific steps are as follows: (i) Step 1. According to the nature of the state transition and full probability formula, deduce the equation of Markov chain.
where a i (t) � P x i � i is the probability of the sequence X(t) in state E i and p ij is the probability of sequence X(t) to be transferred from the first state According to the prediction results of the inverse triple exponential smoothing model, the error between it and the actual value is calculated, the corresponding interval threshold is set, and the error is divided into several intervals. e error selection here usually adopts relative error or actual residual error [28]. (iii) Step 3. Calculate the initial probability.
Assuming that the definition domain of time series . , E r according to certain requirements. For this time series }, all we know is the transition state observed previously, and the transition state of the last term is unknown. Calculating the initial probability requires the number of E i state data in the state of the previous n − 1 data. Suppose that M i data are in the state; the occurrence frequency of state E i is (iv) Step 4. Construct Markov state transition matrix and calculate the transition probability. In this paper, the transition probability of error from one state to another is defined as the state transition probability matrix. For example, the number of transitions from state E i to state E j is m ij , and the total number of transitions from state E i is M i ; then, the probability of transitions from state E i to state E j is p ij � (m ij /M i ), and a state transition probability matrix P is constructed [29,30]: After predicting backward n times from the current state, the transition probability matrix correspondingly becomes P n . (v) Step 5. Error correction.
According to the actual situation of the observed data falling into the state E i , it is also clear that the ith row element of P is p ij � max p i1 , p i2 , . . . , p in because the probability of transition to E j may be greater than other states, so the final transition state can be predicted. And according to Step 4, after the state transition probability matrix is determined, the next state that is most likely to occur can be predicted, and the error interval of the predicted state can be grasped so that the error of this state can be reasonably estimated. Let us define the error to be Δy t+T . In this paper, the median value of the range of state values is used to correct Markov error. (vi) Step 6. Output the corrected predicted value.
e overall flowchart of building the model is shown in Figure 1. First is the selection stage of the research object. e trend analysis, seasonal analysis, periodic analysis, and other methods are adopted to identify the time series and determine the applicable model. In this paper, the exponential smoothing model is selected according to the data instance. Secondly, the traditional exponential smoothing model is studied. Different principles are adopted to determine the relevant parameters, which are substituted into the recursive formula to calculate the coefficient of the exponential smoothing model. e third is the improvement stage of exponential smoothing model. e inverse prediction method is used to determine the number of iterations and other parameters, and the initial value is corrected so as to output the inverse triple exponential smoothing initial value correction model. In the fourth part, an improved inverse exponential smoothing and Markov combination prediction model is established. On the basis of inverse triple exponential smoothing model, Markov theory is introduced to correct the fluctuation error by dividing the state interval and establishing the probability transition matrix. Finally, in the stage of verification and analysis, the accuracy of the traditional exponential smoothing model, inverse triple exponential smoothing model, and improved prediction model combining inverse exponential smoothing and Markov chain was tested, respectively, and the prediction effects were analyzed to output the optimal prediction model.

Data Selection.
is paper selects India's defense expenditure data from 1990 to 2017 (Table 1) to verify the prediction accuracy of the model. e self-fitting of MATLAB data shows that this group of data shows an obvious quadratic parabola trend (Figure 2), and its growth trend is obvious and relatively stable, which meets the data requirements of the triple exponential smoothing model. Among them, India's defense expenditure data from 1990 to 2012 are taken as the original time series, and India's defense expenditure data from 2013 to 2017 are taken as the predicted test data.

Data Inspection Method.
At present, the commonly used data inspection methods mainly include the following (y j represents the actual value of the j th time series data, and y j represents the predicted value of the j th time series data).
e accuracy division range of AARE is shown in Table 2.
AARE is usually expressed as a percentage. e smaller the value, the higher the prediction accuracy.
e value of IC is between 0 and 1. e closer to 1, the worse the prediction accuracy, and the closer to 0, the higher the prediction accuracy.

Prediction Results and Analysis of Results.
e data of India's defense expenditure from 1990 to 2012 are substituted into the traditional triple exponential smoothing prediction model, and the initial value is set as the average of the first three data:   Table 3 and Figure 3, and the predicted values are shown in Table 4(model I refers to the   triple exponential smoothing model and model II refers to  the reverse triple exponential smoothing model). From Table 5, it can be seen that the MAE and AARE indexes of model II are lower and the IC indexes are higher than those of model I. According to the fitting data of the two models from 1991 to 2012, it can be found that the errors of the first three data are too large, which seriously affects the evaluation of indexes. is is because the reverse prediction model carries out reverse prediction correction for the first three initial values, while for the original data less than or equal to three items, it does not meet the reverse prediction conditions, resulting in excessive error of the first three fitting values. erefore, the first three fitting data are removed to make a precision comparison (see Table 6).
As can be seen from Table 6, after removing the first three fitting values, the three indexes of model II are better than those of model I. It should be emphasized that the cubic exponential smoothing model in this paper is established based on dynamic data. erefore, when predicting the leading period number T �1, 2, 3, 4, and 5, respectively, the model comprehensively considers the error results of the fitting value and the predicted value and calculates the optimal smoothing coefficient α so as to generate the corresponding parameter values a 2012 , b 2012 , and c 2012 , respectively (predicting the benchmark year is 2012), as shown in Table 7.
According to models I and II, India's defense spending from 2013 to 2017 is predicted, respectively, and compared with actual data (Table 4), and the prediction precision of model I and model II is compared (Table 8).
It can be seen from Tables 4 and 8 that the prediction effect of model II was better than that of model I from 2012 to 2017, with the mean absolute error (MAE) reduced by 12.49%, the average absolute relative error (AARE) reduced by 14.89%, and the inequality coefficient (IC) reduced by 6.67%. Overall, the prediction accuracy is improved compared with the traditional exponential smoothing prediction model. However, for specific data, it can be seen that the absolute error of reverse triple exponential smoothing model II is obviously reduced on the predicted values from 2013 to 2015, but for the predicted data from 2016 to 2017, the error has basically not changed. erefore, on the basis of model II, Markov theory is introduced to further correct the error. According to the calculation results of model II, its relative error is divided into several states, and according to Markov theory, the concentration degree of error range is divided so that each interval meets the objective law of state change [31][32][33][34]. e standard of state interval division is shown in Table 9. According to the relative error of the fitting values obtained in model II, we can divide them into different states, as shown in Table 10.
Since the state transition is random, probability must be used to describe the possibility of the state transition, that is, the state transition probability [35,36]. According to the relevant probability calculation formula in Section 3 (Improved Prediction Model Combining Inverse Exponential Smoothing and Markov Chain), the single state transition probability matrix is constructed as follows: As can be seen from Table 10, the data in 2012 are in state E3, and according to the state transition matrix, it can be seen that the predicted data in 2013 are most likely to be in state E3. e median value of the state value range is taken for Markov error correction, and combined with the prediction value of reverse triple exponential smoothing model, the final prediction value for 2013 is 2786.298 billion rupees.
Similarly, the n times state transition probability matrix P n can be obtained according to the one-time state transition probability matrix P. After the transition matrix is determined, the possible states in 2013-2017 can be predicted according to the state in 2012, and then the final prediction values of an improved reverse exponential smoothing and Markov combination prediction model (set as model III) for 2013-2017 can be obtained, respectively, as shown in Table 11.      According to the three indexes of MAE, AARE, and IC, the prediction accuracy of the three models are compared, as shown in Table 12 and Figures 4-6.
From the analysis of Table 12 and Figures 4-6, it can be seen that the prediction accuracy of model III is significantly improved compared with models I and II, and its MAE, AARE, and IC are greatly reduced. e MAE index of model III is 43.30% lower than model I, the AARE index of model III is 40.76% lower than model I, and the IC index of model III is 52.03% lower than model I. For specific data, the prediction error of model III increased from 2013 to 2014, but its prediction error from 2015 to 2016 decreased significantly compared with that of models I and II, and the overall curve was closer to the actual value.

Conclusion
In this paper, triple exponential smoothing model I, inverse triple exponential smoothing model II, and improved prediction model III combining inverse exponential smoothing and Markov chain are established, respectively. Taking India's defense expenditure data from 1990 to 2012 as the original time series and the data from 2013 to 2017 as the unknown test values, this paper draws the following conclusions through comparative analysis: (1) e quadratic curve fitting trend of India's defense expenditure data from 1990 to 2017 is good. Verified by models I, II, and III, the average relative error of the fitting values is below 10%, and the average relative error of the predicted values is below 4%, with high accuracy, which shows that it is more reasonable to use the triple exponential smoothing model to predict India's defense expenditure. (2) Compared with model I, the overall prediction accuracy of model II is improved. Its MAE is relatively reduced by 12.49%, AARE is relatively reduced by 14.89%, and IC is relatively reduced by 6.67%, but its prediction advantage is weakened in the later period.     Mathematical Problems in Engineering respectively. On the basis of making full use of time series to predict the trend, the improved prediction model combining inverse exponential smoothing and Markov chain uses Markov theory to predict the fluctuation, thus obviously reducing the error fluctuation range and making the predicted value closer to the actual value and the prediction effect more stable.
is paper only predicts the national defense expenditure from the perspective of time series, but in fact, national defense expenditure will inevitably be affected by multidimensional factors such as economic development, military construction, and national defense policies [37,38]. In the following research, if the time series and multidimensional influencing factors can be integrated to establish a prediction model, it will more effectively support the prediction decision.
Data Availability e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.