Multistep-Ahead Stock Price Forecasting Based on Secondary Decomposition Technique and Extreme Learning Machine Optimized by the Differential Evolution Algorithm

.e prediction research of the stock market prices is of great significance. Based on the secondary decomposition techniques of variational mode decomposition (VMD) and ensemble empirical mode decomposition (EEMD), this paper constructs a new hybrid prediction model by combining with extreme learning machine (ELM) optimized by the differential evolution (DE) algorithm. .e hybrid model applies VMD technology to the original stock index price sequence to obtain different modal components and the residual item, then applies EEMD technology to the residual item, and then superimposes the prediction results of the DE-ELM model for each modal component and the residual item to obtain the final prediction results. In order to verify the validity of the model, this paper constructs a series of benchmark models and, respectively, tests the samples of the S&P 500 index and the HS300 index by one-step, three-step, and five-step forward forecasting. .e empirical results show that the hybrid model proposed in this paper achieves the best prediction performance in all prediction scenarios, which indicates that the modeling idea focusing on the residual term effectively improves the prediction performance of the model. In addition, the prediction effect of the model combined with the decomposition technology is superior to the single DE-ELM model, where the secondary decomposition technique has a significant decomposition advantage compared to the single decomposition technique.


Introduction
e stock index is an important indicator reflecting the development of the stock market. e stock index price forecast can not only provide investors with a reasonable investment basis, but also help relevant policy makers to monitor systemic risks and asset price bubbles in the stock market in a timely manner. However, the price fluctuations of the stock index are affected by complex factors such as the development of the overall financial market and investor sentiment [1][2][3]. Compared to a single stock, the price of a stock index has more uncertainty. It has a series of characteristics such as nonstationarity, nonlinearity, and long memory.
erefore, the majority of scholars have been committed to continuously improving the accuracy of stock index prediction [4][5][6].
According to the different variables involved in predicting stock index price models, it can be divided into two categories: multivariable prediction and univariate prediction. Multivariable forecasting refers to that the forecasting model includes not only historical stock price data, but also technical analysis tools, various macroeconomic variables, and investor sentiment that affect stock index prices [7][8][9]. Among them, technical analysis tools have been widely used by stockbrokers. In the research of stock price prediction, technical analysis tools are often combined with other prediction models, especially artificial intelligence models [10,11]. e univariate prediction refers to that the input variables in the prediction model only include historical stock index price data [12][13][14][15]. is article belongs to the research scope of univariate prediction.
Scholars make in-depth research on the prediction of univariate prediction. ere are three main types of prediction models: traditional econometric models, artificial intelligence models, and hybrid models. Traditional econometric models include the autoregressive conditional heteroscedasticity (ARCH) model [16], generalized autoregressive conditional heteroscedasticity (GARCH) [17,18], autoregressive moving average (ARMA) model [19,20], autoregressive integrated moving average (ARIMA) model [21], and vector autoregressive model [22]. However, the traditional econometric models often need to follow a series of strict assumptions when dealing with data, which are usually only suitable for dealing with stable linear data and do not have good prediction ability for the stock index. erefore, the artificial intelligence model with good prediction ability for nonlinear time series has been widely used in stock market price prediction in recent years. e artificial intelligence model improves the prediction accuracy of the model by training historical price data. Common representative models include the artificial neural network (ANN) [23][24][25], feedforward neural network (FFNN) [26], support vector machine (SVM) [27,28], and long short-term memory neural network (LSTM) [29,30]. Unfortunately, stock index prices are highly volatile, and there is a lot of noise in the data series. However, the traditional feedforward network has difficulty in adjusting model parameters, and the model does not have good antinoise ability. A large number of studies have shown that the prediction effect of hybrid models is often better than that of single models [31][32][33]. e typical method in the hybrid model is the TEI@I complex system research methodology [34]. e method firstly uses the decomposition algorithm to process the original sequence with a large amount of noise and then models each decomposed component with the prediction model. is strategy significantly improves the prediction accuracy of time series and has been widely used in various fields of prediction in recent years [35][36][37][38][39].
Common decomposition algorithms include wavelet decomposition (WD) [40] and empirical mode decomposition (EMD) [41]. Compared with WD, EMD is more suitable for processing nonlinear complex sequences. Wei [42] proposed a mixed time series model based on empirical mode decomposition to predict the stock prices of the Taiwan Stock Exchange Capitalization Weighted Stock Index (TAIEX) and Hang Seng Stock Index (HSI). e empirical results prove the effectiveness of the EMD algorithm in forecasting time series. Cheng and Wei [43] combined empirical mode decomposition and the support vector regression model to build a mixed model for forecasting time series, taking the stock price of Taiwan Stock Exchange (TAIEX) as an empirical target. e empirical results show that the mixed model is superior to the autoregressive model and the support vector regression model in forecasting performance.
e EMD algorithm has the disadvantage of easily superimposing various modes in the result of decomposing the original sequence. erefore, the ensemble empirical mode decomposition (EEMD) [44] was further proposed and applied. Li and Feng [45] applied ensemble empirical mode decomposition (EEMD) to the processing of nonlinear and nonstationary financial time series. Fang [46] used ensemble empirical mode decomposition (EEMD) to explore the prediction ability of investor sentiment in different periods of stock market returns. Although EEMD eliminates the phenomenon of mode aliasing to some extent, the EEMD algorithm cannot completely separate components with similar frequencies. In addition, the EEMD algorithm has the defect of overenveloping or incomplete enveloping the data, and the decomposition of the original data is still insufficient.
In recent years, as an improvement to the above decomposition algorithm, variational mode decomposition (VMD) [47] has performed excellent in processing signals interfered by noise. VMD can effectively separate components with similar frequencies and improve the accuracy of original sequence decomposition. It has been widely combined with various prediction models by scholars in empirical research and applied in the prediction research of the energy price market or the financial market [36,48]. Lahmiri [49] combined variational mode decomposition (VMD) and the backpropagation neural network (BPNN) to construct the VMD-PSO-BPNN model to predict stock prices. e empirical results show that compared with the benchmark model PSO-BPNN, the VMD-PSO-BPNN model has advantages in prediction performance. Bisoi et al. [50] combined the variational model decomposition (VMD) and the optimized limit learning organization to build a DE-VMD-RKELM model to predict the BSE S&P 500 Index (BSE), Hang Seng Index (HSI), and Financial Times Stock Exchange 100 Index (FTSE). Wang et al. [51], based on fast ensemble empirical mode decomposition (FEEMD), VMD, and backpropagation (BP) neural network, established a mixed model of two-level decomposition to predict the electricity price.
As a method to form the prediction part of the TEI@I complex system research methodology, any regression method is adequate. However, there are differences in the prediction accuracy and the completion time of different models. e artificial intelligence models mentioned above (SVR, ANN, and LSTM) are better choices than traditional econometric models. Unfortunately, the prediction ability of the traditional feedforward network is limited by parameter adjustment, which has the defects of slow calculation speed and unstable generalization ability. erefore, extreme learning machine (ELM) [52] has been widely used in various prediction fields due to its advantages in learning convergence speed and parameter setting [53,54].
In summary, in the existing stock index price prediction research, the TEI@I complex system research methodology represented by "decomposition-prediction" has gained recognized advantages in the prediction field. However, the existing research still has the following deficiencies: (1) In the existing research using VMD technology, the residual term is discarded or treated as an ordinary component. None of these studies focused on the residuals that contain rich information.
(2) Existing studies using ELM algorithms only reduce the reliability of prediction results and reduce the stability of model prediction performance. (1) e secondary decomposition technique combining VMD and EEMD is applied to the price sequence, in which EEMD is used to process the residual term of VMD. e secondary decomposition technology can better capture the overall features of the data and improve the overall prediction accuracy of the hybrid model. (2) e DE algorithm, which has been proven effective [44,56,57], is introduced to optimize the input weight and hidden layer thresholds of ELM during the training process, thus improving the stability and reliability of ELM prediction results.

Introduction of Methodology
is paper builds a combination model (VMD-RES.-EEMD-DE-ELM) based on the secondary decomposition technique and the machine learning method to predict the price of grain futures. Before building the model, it is necessary to briefly introduce the components of the model combination: VMD technology, EEMD technology, extreme learning machine, differential evolution algorithm, and VMD-RES.-EEMD-DE-ELM hybrid model constructed in this paper.

VMD.
VMD is used to decompose the original input signal into k band-limited intrinsic mode functions with certain sparsity. In this paper, it is simply called VMF, i.e., the set of u k , in which the signal evolution of each mode is expanded around the center frequency w k . By minimizing the sum of the bandwidth estimates of all modes, the corresponding mode component signal u k and related parameters can be obtained. e VMD signal decomposition process is also the solution process of the variational constraint problem, which is shown as follows: where u k : � u 1 , . . . , u K are each modal component VMF after VMD; w k : � w 1 , . . . , w K are, respectively, the center frequencies corresponding to each VMF; z t represents the partial derivative to t; δ(t) is the impact function; * is the convolution symbol; f is the original input signal. e analytic signal of correlation u k (t) is obtained by the Hilbert transform to obtain its unilateral spectrum. e exponential term e − jw k t is used to adjust the estimated value of each w k and integrate the spectrum of u k into the basic frequency band. In order to solve the above constrained variational problem, the constrained variational problem needs to be converted into an unconstrained problem. By introducing the quadratic penalty factor α and the Lagrange operator λ(t), the following Lagrange expression describing the unconstrained variational problem is obtained: where the penalty factor α is used to ensure the accuracy of signal reconstruction; Lagrange operator λ(t) can maintain strict constraint conditions. Further, an alternative direction method of multiplications (ADMM) is adopted to carry out iterative search to find the saddle point of the Lagrange function, thus obtaining the optimal solution of the unconstrained variational problem of equation (2). e expressions of VMF u k and center frequency w k are, Mathematical Problems in Engineering 3 respectively: e specific implementation steps of the VMD decomposition method are as follows: Step 1. Initialize the values of parameters such as VMF and center frequency, that is, u 1 k w 1 k , n � 0, and select the appropriate number of VMFs, that is, the value of K.
Step 2. Update the values of u k and w k according to equations (3) and (4), respectively.
Step 3. Update the value of λ: Step 4. Given the determination accuracy ε > 0, if the following condition is met: en, stop the iteration; otherwise, return to Step 2. In the above equation, u n k (w), f(w), and λ n (w) are the Fourier transforms corresponding to u n k , f(t), and λ n .

EEMD.
Compared with empirical mode decomposition (EMD), EEMD has stronger adaptability and local variation characteristics [44] and can effectively identify the nonlinearity and nonstationarity of time series. For the unparsed original signal, the brief implementation of the EEMD algorithm is as follows [56]: Step 1. After adding white noise sequence Step 2. x ′ (t) is decomposed into IMFs and residuals.
Step 3. Repeat steps 1 and 2 p times, using different white noises each time (p is the integration number).
Step 4. e influence of white noise on the original components is eliminated by calculating the average of all IMF components and residual terms. So far, the complete time series can be expressed as the sum of IMF and residual terms: where IMF i (t) is the inherent modal function, r m (t) represents the final residual component, m is the total number of IMF, and i is the modal index. e IMF should meet the following two conditions [41,57]: first, the number of extreme values and the number of zero crossings must be equal or the difference between the two must be at most 1 over the entire length of the IMF; secondly, in any case, the average value of the upper envelope defined by the local maximum and the lower envelope defined by the local minimum is zero. erefore, IMF 1 represents the maximum amplitude and the highest frequency, while the other IMF subsequently has lower amplitude and frequency, and the residual term represents the slow change process around the long-term average.

DE-ELM
2.3.1. ELM. Extreme learning machine (ELM) was originally proposed by Huang Guangbin. As a special kind of the single-hidden layer feedforward neural network (SLFN), it contains only one hidden layer, and its network structure is shown in Figure 1.
In the ELM network training process, the connection weights of the input layer and the hidden layer and the threshold value of the hidden layer can be randomly given and need not be adjusted. Assuming there are N arbitrary training samples ( given the hidden layer number of ELM is L, the randomly generated connection weight matrix of the input layer and the hidden layer is W i � [w i1 , w i2 , . . . , w in ], and the offset of the hidden layer is b i i � 1, . . . , L; the network output can be expressed as where T is the connection weight between the hidden layer and the output layer, g(x) is the corresponding activation function, and W i .X j is the inner product of W i and X j . e goal of the ELM network training is to make the error function reach the minimum, that is: Given the input weight W and the hidden layer offset b (equation (9)), the solution process of the network is transformed into searching for the optimal output weight β e matrix expression is where H (equation (12) is the output matrix of hidden layer nodes: An expression further derived from formula (13) is as follows: where H † (equation (14)) is the left inverse matrix of H. erefore, once the input weight W and the hidden layer offset b are randomly given in the ELM training process, the output matrix H of the hidden layer can be uniquely determined without adjusting the parameters in the training process.

Differential Evolution.
DE algorithm is an evolutionary computation method based on population difference and has strong global optimization ability. Firstly, a mutation operation is performed on the difference of parent individuals to generate mutation individuals. en, the parent individual and the variant individual perform crossover operation according to a certain probability to generate mixed individuals. Finally, the selection operation of survival of the fittest is carried out between the parent and the mixed individuals according to the fitness to achieve the evolution of the population. e basic process is shown in Figure 2. e basic steps are as follows: Step 1. Initialize the related parameters: for example, the population size is set to NP; the mutation factor is F; the cross factor is CR; the spatial dimension is D; the evolutionary algebra t � 0.
Initialize the parent population: Step 3. Calculate the fitness value of each individual, that is, the objective function of the model.
Step 4. Perform the mutation operation: the mutant T are generated by the difference between the parental individuals x t i . e generation process is shown in the following equation: where x t r 1 , x t r 2 , and x t r 3 (equation (15)) are, respectively, the three parent individuals randomly selected from the population and r 1 ≠ r 2 ≠ r 3 ≠ i.
Step 5. Crossover operation: a crossover operation is performed between the parent individual x t i and the mutant individual v t i to generate a mixed individual u t i � (u t i1 , u t i2 , . . . , u t i D ) T . e generation process is as follows:

Mathematical Problems in Engineering
where rand[0, 1] is a random number valued in the range [0, 1] and j rand is a random integer valued at [1, D].
Step 6. Selection operation: the differential evolution algorithm keeps excellent individuals and eliminates inferior individuals through continuous evolution. It selects the best one with the best fitness value from the parent individual x t i and the mixed individual u t i as the next generation individual x t+1 i , and the selection process is shown in the following formula: where fitness(·) is the fitness function, which is the objective function that needs to be optimized, and the purpose is to find its minimum value.
Step 7. Iteration termination test: if the error requirement is met or the maximum number of iterations is reached in the generated next-generation population NP , the iteration is stopped and the optimal individual is output. Otherwise, mutation, crossover, and selection operations continue until the iteration stop condition is met.

Construction of DE-ELM Model.
Since the input weight matrix and hidden layer threshold of the model are randomly set during the network training of the basic ELM, the training stability of the model is weakened. In order to improve the training stability of the ELM, the DE algorithm is adopted to optimize the input weight matrix and hidden layer threshold of the ELM model. It is worth mentioning that in the ELM model, the optimization of the input weight matrix and the hidden layer threshold needs to first determine the number of hidden layer neurons. Based on experience, this paper sets the number of hidden layer neutrinos in the ELM network at 35. On this basis, the DE-ELM model is proposed in this paper. Its specific steps are as follows: Step 1. e random input weight and the hidden layer threshold of the ELM model are encoded to initialize the population.
Step 2. Initialize relevant parameters of the DE algorithm. e mutation factor of the DE algorithm is 0.5, the crossover factor is 0.9, the maximum number of iterations is 20, the dimension of each individual is 245, and the population size is 10. is paper is based on the same settings as above for all hybrid models.
e fitness value of each individual in the population is calculated, which is the root mean square error value of the predicted output of the ELM model.
Step 4. Mutation, crossover, and selection operations are performed in sequence, and iterative termination check is performed. If the termination condition is satisfied, the optimal input weight and the hidden layer threshold are output; otherwise, the iteration is continued until the optimal output is obtained.
Step 5. e input weight and hidden layer threshold of the ELM model are set as the optimal individuals optimized by the DE algorithm, thus obtaining the optimized and improved DE-ELM model.
In the above steps of optimizing the ELM model, the specific decision variables are the input weight matrix and the hidden layer threshold, and the goal of optimization is to minimize the value of the root mean square error (RMSE) index as a function of fitness. e specific formula of RMSE can be found in equation (19).
In addition, the input layer dimension of the ELM network in the DE-ELM model is set to 6; that is, the data of the 7th day are predicted based on the data of the previous 6 days, and by analogy, rolling prediction is carried out. at is, , where x t is the predicted value of the sequence at time t and x t−1 , . . ., x t−6 are the values of the previous 6 days of the sequence at time t, respectively.

Construction of VMD-RES.-EEMD-DE-ELM Hybrid
Model. As we all know that the price sequence of the stock index shows nonstationary and nonlinear complex characteristics. e VMD technology can effectively remove noise from the original price series and extract the main features hidden in the original price series. In addition, unlike EEMD technology, the residual term remaining after VMD technology processes the original sequence contains rich information. However, previous studies directly discarded the residual term or treated the residual term as a common component, which failed to capture the complete characteristics of the time series and reduced the overall prediction accuracy of the combined model. erefore, this paper proposes a secondary decomposition technology based on the combination of VMD and EEMD technology. After VMD processes the original sequence, EEMD technology is used to further process the residual term. Further, this paper combines the forecast advantages of the DE-ELM model to build the VMD-RES.-EEMD-DE-ELM hybrid model. e step diagram of visualization is shown in Figure 3. e detailed modeling steps of the model are as follows: Step 1. VMD technology is applied to the original price sequence of the stock index to obtain various modal components (denoted as VMF) and the residual item containing complex information.
Step 2. e different modal components are normalized. After selecting an appropriate number of training sample data and test sample data, the DE-ELM model is applied to train and test the training samples to obtain the respective prediction results of each modal component.
Step 3. e EEMD technique is applied to the residual term to obtain different modal components (denoted as IMF). e content in step 2 is repeated, and the DE-ELM model is used to predict each component. After the superposition, the predicted value of the residual term can be obtained.
Step 4. e predicted values of each modal component obtained after applying the VMD technology are superimposed, and finally the predicted values of the original price series are obtained.

Sources and Processing of Data.
is paper selects the S&P 500 index and the HS300 index of mainland China as representatives of the stock market [58][59][60] to test the prediction accuracy of the VMD-RES.-EEMD-DE-ELM hybrid model constructed in this paper. e S&P 500 index covers a wide range of stocks in the market, while the HS300 index consists of 300 large-scale and highly liquid-constituent stocks on China's Shanghai and Shenzhen Stock Exchanges, both of which are representatives of the securities market.
e sample of the S&P 500 index selected in this paper covers a period of 1,511 transaction data from January 2, 2014, to December 27, 2019. e sample of the HS300 index covers 1,465 transaction data from January 2, 2014, to January 2, 2020. All models in this paper are run by using MATLAB 2019b software package. All the data come from the Choice Financial Terminal of Oriental Fortune. e samples in this paper are composed of the training set and the test set, respectively. In the sample of the S&P 500 index, the first 1361 data are training sets, and the next 150 data are set as test sets. In the sample of the HS300 index, the first 1315 data are training sets, and the next 150 data are set as test sets. Figure 4 shows the trend of the S&P 500 index and the HS300 index samples. Table 1 lists the information of descriptive statistical analysis of the S&P 500 index and the HS300 index, including the maximum (Max), minimum (Min), average (Ave), median (Med), and standard deviation (Std).
In addition, in order to further test the robustness of the model, all models in this paper are predicted by one-step, three-step, and five-step-ahead forecasting; i.e., the data of the previous 6 trading days are used to predict the data of the 7th, 9th, and 11th trading days, respectively. at is, in all the prediction models constructed in this paper, the input variables of the model are the data of the previous 6 trading days, and in the one-step-ahead forecasting, the output variable is the forecast data of the 7th trading day; in the three-step-ahead forecasting, the output variable is the forecast data of the 9th trading day; in the five-step-ahead forecasting, the output variable is the forecast data of the 11th trading day.
When the original sequence is processed by VMD technology, the number of preset modal components must be given in advance. Assuming that the number of   decomposed modal components is n, this paper chooses n as 10 according to the principle that the n-power of 2 does not exceed the sample data. Furthermore, in the process of applying the DE-ELM model to train the data, in order to ensure the training effect of the DE-ELM model, this paper adopts the normalization process of "min-max dispersion standardization" to the sequence data of each modal component generated after decomposition. e specific mathematical expression is as follows: where x ′ is the normalized subsequence data, x is its original value, and x max and x min are its maximum and minimum values, respectively.

Evaluation Index of Model Prediction Results.
In order to visually compare the prediction effects of different prediction models. is article introduces three evaluation indicators of prediction performance, namely, root mean square error (RMSE), mean absolute error (MAE), and mean absolute percentage error (MAPE). e specific calculation formula is as follows: where y i and y i are the real and predicted values of the S&P 500 index and the HS300 index, respectively, n is the scale data of the test sample, i is the serial number of the test sample points, and e RMSE , e MAE , and e MAPE indicate the prediction accuracy of the model; the smaller the value of the index, the higher the prediction accuracy of the model.

Comparative Analysis of Model Prediction Effects.
In order to prove the superiority of the VMD-RES.-EEMD-DE-ELM model, this paper constructs seven other different benchmark models for comparison. Figure 5 shows the multistep forward prediction results of each model for the S&P 500 index and the HS300 index, respectively. In addition, Tables 2-4, respectively, record the specific values of evaluation indexes generated by the eight prediction models for the S&P 500 index and the HS300 index prediction. e empirical results show that the VMD-RES.-EEMD-DE-ELM hybrid model proposed in this paper has achieved the best performance in all cases compared with all other benchmark models in the prediction results of the S&P 500 index and the HS300 index. Detailed analysis of the results is described below.  Table 2 and the prediction model combined  with decomposition technology in Table 3, it can be found that the prediction accuracy of the combined model has an overwhelming advantage over the noncombined model. e prediction accuracy of the single prediction model is poor,  In addition, by comparing the prediction results of the EEMD-DE-ELM model and the VMD-DE-ELM model, it can be found that the prediction accuracy of the VMD-DE-ELM model is better than that of the EEMD-DE-ELM model in the same decomposition level. Although the EEMD-VMD-DE-ELM combination model with secondary decomposition technology has achieved better prediction results than the EEMD-DE-ELM model, its performance is still inferior to the VMD-DE-ELM model. e results show that the performance of VMD technology in extracting sequence data features and processing complex signals is better than EMMD technology. Table 4 is a combined model that uses the DE-ELM model to predict each VMF component and residual term after applying VMD technology to the original sequence. e VMD-RES.-EEMD-DE-ELM model is the hybrid model proposed in this paper.

Comparative Analysis of Mixed Models considering Residual Decomposition. VMD-RES.-DE-ELM in
By comparing the prediction results in Tables 3 and 4, it can be known that in the results of the multistep forward    DE-ELM model, it can be seen that the residual term generated by applying VMD technology to the original sequence does contain important information, and its complexity is relatively high. When the DE-ELM model is used to directly predict the residual term, its prediction effect is very limited. erefore, in this study, the VMD-RES.-EEMD-DE-ELM model is the best choice, and the empirical results reasonably confirm that the proposed secondary decomposition technology can combine the advantages of VMD and EEMD to generate smoother subsequences, thus further obtaining more accurate prediction results.

Conclusion
Accurate stock index price predictions are not only of great significance to regulators, but also crucial to the stability of the financial system. Based on the advantages of VMD, EEMD technology, and extreme learning machine, this paper proposes a combined model VMD-RES.-EEMD-DE-ELM, which combines secondary decomposition technology and the artificial intelligence algorithm. Further, this article takes the S&P500 index and the HS300 index as experimental samples and e RMSE , e MAE , and e MAPE as evaluation indexes to test the performance of the VMD-RES.-EEMD-DE-ELM model compared with other seven benchmark models in the results of one-step, three-step, and five-step forward prediction. e conclusions of empirical analysis are as follows: (1) e VMD-RES.-EEMD-DE-ELM hybrid model proposed in this paper not only takes full advantage of the secondary decomposition, but also solves the problem that the residual term is not fully considered in the traditional time series prediction model based on VMD technology. In the prediction research of the price of the S&P 500 index and the HS300 index, the VMD-RES.-EEMD-DE-EELM model has achieved the best performance in all cases, indicating that the hybrid model proposed in this paper can fully capture the characteristics of the original sequence and has good prediction performance. (2) In the prediction research of the S&P 500 index and HS300 index prices, the VMD technology has higher decomposition accuracy of the original sequence than the EEMD technology. In the empirical prediction test combined with the DE-ELM algorithm, the VMD technology can better improve the model prediction effect than the EEMD technology. (3) In the prediction research of the S&P 500 index and HS300 index prices, the secondary decomposition technology has demonstrated its superiority compared to the single decomposition technology. e two-layer decomposition technology can better handle complex original sequences and effectively improve the prediction accuracy of the hybrid model combined with the prediction model.
However, the trend of the stock market's price series has complex influencing factors, such as the economic development of different countries, geopolitical conflicts, and the effects of linkages between international markets, the scale of industrial development and business conditions of different industries. erefore, multidimensional and complex influencing factors can be included in the research scope in future research to further improve the overall prediction effect.

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.