Multistep Wind Power Prediction Using Time-Varying Filtered Empirical Modal Decomposition and Improved Adaptive Sparrow Search Algorithm-Optimized Phase Space Reconstruction–Echo State Network

: Accurate wind power prediction is vital for improving grid stability. In order to improve the accuracy of wind power prediction, in this study, a hybrid prediction model combining time-varying ﬁltered empirical modal decomposition (TVFEMD), improved adaptive sparrow search algorithm (IASSA)-optimized phase space reconstruction (PSR) and echo state network (ESN) methods was proposed. First, the wind power data were decomposed into a set of subsequences by using TVFEMD. Next, PSR was used to construct the corresponding phase space matrix for sequences, which were then divided into training sets, validation sets, and testing sets. Then, ESN was used for subsequence prediction. Finally, the predicted values of all the subseries were used to determine the ﬁnal predicted power. To enhance the model performance, the sparrow search algorithm was improved in terms of the discoverer position update strategy, the follower position update strategy, and the population structure. IASSA was employed to synchronously optimize multiple parameters of PSR-ESN. The results revealed that the proposed model has higher applicability and prediction accuracy than existing models.


Introduction
Due to the exacerbation of the global energy crisis, renewable energy sources, such as wind, solar, and hydropower, have garnered increased attention. While wind energy represents an environmentally friendly source of energy, its unpredictable nature renders the integration of wind power into the power grid a cause for concern, frequently resulting in a decline in the grid's power quality [1].
Therefore, accurate prediction of wind power can effectively improve the stability of the power system. With the development of new power grids, the prediction of wind power forecasting plays an increasingly important role [2,3]. In recent years, power forecasting has become increasingly important in electric system planning studies. The progress in wind power research is of mutual scientific value to the development of other renewable energy sources, such as photovoltaic power generation, hydroelectric power, and nuclear energy. With the increasing penetration of smart grids, it is obvious that power forecasts play an essential role in the intelligence of the power grid [4][5][6].
Existing methods for improving wind power prediction mainly include physical models, statistical models, and machine learning (ML) techniques [7][8][9]. Physical models combine physical factors such as weather, temperature, and wind direction to estimate future wind power [10][11][12]. However, the development of such models is limited by data resources. Statistical models predict future wind power by mining information from his-torical data, mainly by using autoregressive, moving autoregressive, and multi-class autoregressive moving average (ARMA) [13,14]. However, solving the nonlinearity problem in time series by using these models is difficult, thereby hindering high-accuracy prediction. Compared with physical models and statistical models, prediction models based on ML methods, such as extreme learning machine (ELM), backpropagation neural network (BPNN), recurrent neural network (RNN), convolutional Neural Networks (CNN), and long short-term memory (LSTM), can better analyze nonlinear time series and have thus been favored by numerous researchers [15][16][17][18][19]. Ding et al. [20] employed the numerical weather prediction wind speed, trend, and detail terms as the inputs of the weighted time series and used the two-way gated recursive unit neural network to correct the wind speed error of the weather forecast and used the modified data to predict the final wind speed. Due to the volatile nature of wind power, a single model cannot provide accurate predic-tions for wind power generation. Ruiz-Aguilar et al. [21] proposed a hybrid prediction model to decompose the clustering preprocessing method by using empirical modal de-composition (EMD) and PE, that inputting the clustered components into the artificial neural network for prediction. However, due to the modal aliasing phenomenon encoun-tered in EMD, the final prediction performance has limited room for improvement. Zhang et al. [22] used complete ensemble empirical mode decomposition adaptive noise (CEEMDAN) to process data, which effectively addresses issues related to modal decom-position following signal partitioning, and effectively resolves the transference of white noise from high to low frequency, thereby elevating the quality of subsequent analyses. Liu and Zhao [23] proposed a prediction model for convolutional networks and long short-term memory (LSTM) networks based on clustering EMD and experimentally demonstrated that this method yields high-accuracy prediction results [24]. Wind power, being a time series signal, is essentially associated with data processing. Consequently, employing intelligent data processing techniques can significantly enhance its efficiency [25,26].
The echo state network (ESN) is a novel form of RNN and can effectively predict nonlinear, chaotic time series by using the echo characteristics of its internal memory pool [27]. In the ESN, the input layer weight, savings pool weight, and connection of ESN are randomly set; thus, it offers the advantages of fast calculation speed and strong generalization ability and is widely used in the field of wind power prediction [28]. However, in ML-based prediction models, setting hyperparameters is challenging. Setting appropriate hyperparameter can greatly affect the prediction results; thus, various optimization algorithms have been proposed to obtain high-accuracy prediction results [29]. Li et al. [30] used the crow search algorithm to improve the ELM, improved the input weight of the ELM, solved the problem of unstable bias of the hidden layer, and improved the accuracy of the model. Samadianfard et al. [31] combined multilayer perceptrons with the whale optimization algorithm (WOA) to enhance the prediction capability. In the present study, an improved adaptive Sparrow Search Algorithm (IASSA) was developed to enhance the search effectiveness of the original SSA and improve the prediction performance of Echo State Network (ESN) by synchronously optimizing Phase Space Reconstruction (PSR) parameters. The application of IASSA has demonstrated significant effectiveness in enhancing the accuracy of wind power prediction.
For wind power forecasting, Du et al. [32] utilized long short-term memory (LSTM) networks. LSTM, as a specialized type of recurrent neural network, has the capability to learn and retain time-series patterns from historical data, thus enabling accurate forecasts. However, when confronted with the non-linearity and non-stationarity inherent in wind power data, LSTM may encounter challenges in achieving high predictive accuracy. Tian et al. [16] employed the Backpropagation Neural Network (BPNN) for wind power forecasting. BPNN is trained using a backpropagation algorithm, which can address the non-linearity of wind power forecasting to a certain extent. However, BPNN has inherent limitations. Firstly, when dealing with the non-stationarity of wind power data, its predictive performance is often restricted. Additionally, BPNN may get stuck in local optima, thereby compromising predictive precision. In contrast, we proposed a short-term wind power prediction model by combining time-varying filtered EMD (TVFEMD), PSR, ESN, and an IASSA-based synchronous optimization strategy. First, TVFEMD was used to preprocess the original wind power time series to reduce the non-stationarity of the wind power time series. Next, the decomposed subsequences were converted into a phase space matrix using PSR. The training, validation, and test sets were divided, and each set of subsequences was predicted using the ESN. At this stage, IASSA was used to synchronize the optimization of parameters in PSR and ESN. Finally, the prediction results of all subcomponents were accumulated to obtain the final predicted value of wind power.
In conclusion, our method exhibits superior performance in tackling the complexity and challenges of wind power forecasting compared to LSTM, BPNN, and other methods such as SSA-ESN and EMD-SSA-PSR-ESN. Our experimental results further confirm that our method significantly outperforms the aforementioned techniques, thereby affirming the efficacy and superiority of our approach. The novel contributions of this paper primarily encompass the following aspects: 1. The proposed model decomposes the wind power of the original time series through TVFEMD, which effectively solves the problem of it being difficult to predict non-stationary signals. A synchronous optimization strategy based on phase space reconstruction (PSR) and echo state network (ESN) is proposed.
2. In order to enhance the performance of the prediction model, the sparrow search algorithm (SSA) is improved in three aspects, namely, finder position update strategy, follower location update strategy and population structure. Additionally, the proposed IASSA is verified in the test function.
3. It is proposed that IASSA synchronously optimize PSR-ESN to achieve optimal prediction.

Time-Varying Filtering-Based Empirical Mode Decomposition
TVFEMD is an improved EMD method based on time-varying filters proposed by Li et al. [33]. TVFEMD provides an alternative scheme for single-component signals, namely local narrowband signals, which greatly improves the problem of modal aliasing encountered in EMD. TVFEMD performs time-varying filtering of the signal in order to determine the local cutoff frequency. The main steps involved are as follows: (1) The original time series x(t) is transformed using the Hilbert technique to obtain the instantaneous amplitude A(t) and instantaneous frequency ϕ (t).
(2) The local maximum value A({t max }) and minimum value A({t min }) of A(t) are interpolated to obtain β 1 (t) and β 2 (t), respectively, and the instantaneous mean a 1 (t) and instantaneous envelope a 2 (t) are calculated: (3) ϕ ({t min })A2({t min }) and ϕ ({t max })A2({t max }) are interpolated to obtain η 1 (t) and η 2 (t), respectively, to calculate the instantaneous frequency: (4) The cutoff frequency ϕ bis (t) is calculated:  (5) To solve the intermittent problem, the cutoff frequency ϕ bis (t) is rearranged. (6) The extreme value of the signal h(t) = cos ϕ bis (t)dt is taken as the node, x(t) is approximated using B-spline interpolation, and the approximation result m(t) is obtained.
(7) It is determined whether the cutoff criterion is met, and the bandwidth threshold ξ is set. If θ(t) ≤ ξ, x(t) is determined as a subcomponent; otherwise, x 1 (t) = x(t) − m(t), and steps 1-7 are repeated: where ϕ avg (t) is the weighted mean instantaneous frequency, and B Loughlin (t) is the instantaneous Loughlin bandwidth.

Phase Space Reconstruction
PSR is a mathematical method for analyzing chaotic time series [34] and can be used to convert discrete one-dimensional time series x(t) into a d-dimensional vector with delay time τ. In this study, the historical load data were reconstructed into the input and output vectors of the prediction model by using PSR, as follows: where {x i | i = 1, 2, . . . , n} represents the original load sequence, n is the length of the sequence, L = n − (d − 1)·τ − h, and h is the prediction step.
In this experiment, we adopted a combined PSR-ESN model; PSR was used to construct the corresponding phase space matrix for sequences, which were then divided into training sets, validation sets, and testing sets. Then, ESN was used for the subsequence prediction.

Echo State Network
The ESN is a variant of the RNN. The network structure of the ESN includes three parts, namely an input layer, a reserve pool, and an output layer, as shown in Figure 1 [35]. Because the input layer weight matrix and the reserve pool internal connection matrix are randomly generated and fixed, the computational amount of model training is reduced. The ESN solves the fitting regression problem of time series connected hidden layer with a sparsely connected reserve pool; th network along with the output equation is as follows: Input Output The ESN solves the fitting regression problem of time series by replacing the fully connected hidden layer with a sparsely connected reserve pool; the updated state of the network along with the output equation is as follows: where tanh(·) represents the activation function and is used to obtain the network echo characteristics, a is the leakage rate used to control the update weight of the ESN network, W in represents a matrix of input weights randomly generated in the range [−1, 1], R is the connection matrix with a sparse structure inside the reserve pool, u(t) represents the input at time t, x(t) represents the t-moment state of the reserve pool, and y(t) is the output at time t. The output matrix W out of the ESN can be solved via ridge regression, with the following optimization objectives: where λ is the regularization coefficient used to prevent the phenomenon of overfitting in the ESN training set, and I is the identity matrix. The prediction data are substituted into Equations (7) and (8) to obtain the final prediction result. The ESN structure is simple and practical; however, its prediction performance is affected by the parameter settings, such as the reserve pool connection matrix scaling parameter R h , number of reserve pool network nodes N, input data scaling coefficient I S , reserve pool sparsity degree S, and leakage value a. Using the appropriate parameter settings can effectively enhance the prediction ability of the ESN.

Sparrow Search Algorithm
SSA is a meta-heuristic optimization algorithm [36] inspired by the foraging and antipredation behavior of sparrows [36,37]. In SSA, individuals are classified as discoverers, followers, and watchers. The natural behavior and mathematical expression of the sparrow can be described as follows: Discoverers: Each generation of discoverers represents a point in the population that is closer to food, and its main function is to provide directions for the entire population to explore food. The mathematical expression of the location update of the discoverer is as follows: where x t i,j represents the position of the i-th sparrow in the j-dimensional when the number of iterations is t, iter max is the maximum number of iterations, R 2 and α are randomly generated values between 0 and 1, ST is the safety threshold, Q is a random number that follows a normal distribution, and L is a vector with all 1s of the elements. When R 2 < ST, the discoverer can search extensively for food; in contrast, R 2 ≥ ST indicates that there is danger and that the discoverer needs to withdraw from the danger zone.
Followers: The role of the follower is to follow the discoverer for food; their update strategy is as follows: where x t worst and x t+1 best , respectively, represent individuals with the worst current fitness and the best individuals, and A is a vector with a value of 1 or −1 randomly generated for each element. i > n/2 indicates that the current follower is in a poor position in the population and that the food is widely searched for by diverging from the worst individual; otherwise, the search for food is performed by competing with the best individuals.
Vigilantes: Vigilantes are a random proportion of individuals in the population, whose main role is to alert the foraging area. Their mathematical expression is as follows: where β is the update step with a mean of 0 and a variance of 1; k is a random number between −1 and 1; f i , f g , and f w are fitness functions of the current alert, optimal individual, and worst individual, respectively; and ε is used to avoid constants with a denominator of 0.
f i = f g indicates that the sparrow is at the edge of the population and thus approaches the optimal individual; and f i = f g indicates that the sparrow is at the center of the population, feels threatened and moves closer to its own species, reducing the risk of predation.

Improved Adaptive Sparrow Search Algorithm
Although SSA offers the advantages of a comprehensive search mechanism and fast convergence speed, it easily falls into local optimum tendencies. In addition, the population ratio of discoverers to followers in SSA must be set empirically; this affects the optimization performance of SSA. In order to improve the comprehensive optimization capability of SSA, we proposed IASSA in this paper. Improvement 1: In the discoverer position update strategy of the original SSA, the use of the mathematical model as a method to control the flight distance by the normally distributed random number Q is not obvious because the random number generated by the standard normal distribution reduces the probability of the discoverer withdrawing from the danger zone at a distance; this affects the global search ability of the algorithm. Therefore, generating Q via the Cauchy distribution allows the discoverer to quickly fly to a distant location. Moreover, varying the safety threshold RT with the iteration time allows the discoverer to search extensively for food in the early stage of the iteration, and as such the discoverer has a greater probability of escaping the current position in the late iteration to ensure population diversity in the later stage of the iteration. The formulas for generating the random number Q and the safety threshold ST are as follows: Improvement 2: In the follower's renewal strategy, the follower searches for food by competing with the best individual. The process of the follower approaching the optimal individual is single, and food can be searched on only one side of the optimal individual close to the direction of the follower. Inspired by the whale optimization algorithm (WOA) [38], in this study we realized the purpose of searching for food around the optimal individual by updating the position by encircling the optimal individual with a spiral of the better followers. The specific process is as follows: where b is a constant that defines the shape of the logarithmic spiral, and l is a random number in [−1, 1]. Improvement 3: In the original SSA, the proportional population selection of discoverers and followers was a complex process. This is because the discoverer mainly provides the ability to search globally, while the follower represents the local search capability in the algorithm. The selection of discoverer and follower ratios also affects the search efficiency of the algorithm; thus, in this paper, we proposed an adaptive population structure strategy and determined the proportion of discoverers in the population by using the following formula: where P is the proportion of t-generation discoverers in the population, and P max and P min are, respectively, the upper and lower limits of the proportion of discoverers. P decays from P max to P max iteratively.
In conclusion, these three improvements are interconnected and collectively optimize the sparrow search algorithm. The first improvement enhances global search capability, the second improvement increases local search and exploration capabilities, and the third improvement further enhances search efficiency and adaptability through the adaptive population structure strategy. Together, these improvements enable the improved sparrow search algorithm to achieve better performance and effectiveness in problem-solving scenarios.

Algorithm Evaluation
To verify the performance of the proposed algorithm, the algorithm was evaluated using multiple benchmark functions, as shown in Table 1. In addition, the proposed algorithm was compared with other commonly used optimization algorithms to better demonstrate the improvement effect of the IASSA.  The optimization results of different algorithms for the benchmark function are listed in Table 2, the convergence effect of different algorithms is illustrated in Figure 2 and the convergence time of different algorithms for various functions are outlined in Table 3. The comparative analysis revealed that the IASSA has better convergence efficiency in unimodal functions F1-F3 and can jump out of local optimal solutions in multimodal functions F4-F6, and it also has a certain advantage in time. To ensure impartiality, the population size and iteration count were set to 30 and 200 across all algorithms. Each aforementioned algorithm was independently executed 30 times, with their respective mean values (AVE) and standard deviations (STD) recorded each time. Herein, the mean value represents the optimal performance of the algorithm, whereas the standard deviation signifies its stability.  The maximum number of iterations was set to 200 based on previous resea preliminary experimentation. It was observed that, when the iteration count excee the convergence of various algorithm components significantly improved or app stability. Considering the constraints of computational resources and time, the i count was set to 200.  The maximum number of iterations was set to 200 based on previous research and preliminary experimentation. It was observed that, when the iteration count exceeded 200, the convergence of various algorithm components significantly improved or approached stability. Considering the constraints of computational resources and time, the iteration count was set to 200.

Wind Power Prediction Model Design
We developed a wind power hybrid prediction model by using TVFEMD, PSR, and ESN. To improve the prediction accuracy, the IASSA optimization algorithm was used to determine the optimal parameters of PSR and ESN simultaneously. The parameters τ and d in PSR and Rh, N, IS, and S in the ESN were encoded together in the IASSA. In addition, the root mean-square error (RMSE) was used as the fitness function to solve the parameters. The steps involved in the hybrid model are shown in Figure 3 and are described as follows:  (1) Raw wind power data were collected.
(2) The original wind power time series was divided into a series of subsequences by using TVFEMD. (3) All subsequences were converted into phase space matrices by using PSR and divided into the training set, validation set, and test set. (4) The ESN model was trained using the training set to verify that the RMSE of the set was the minimum, and all parameters in PSR and ESN were optimized synchronously using IASSA. (5) The optimized optimal parameters were introduced into PSR-ESN, and the test set was ran in order to obtain the prediction results of the subseries. (6) The prediction results of all subcomponents were accumulated to obtain the final wind power prediction results.

Case Analysis
The wind power generation data of the Irish energy system were used for experi- (1) Raw wind power data were collected.
(2) The original wind power time series was divided into a series of subsequences by using TVFEMD. (3) All subsequences were converted into phase space matrices by using PSR and divided into the training set, validation set, and test set. (4) The ESN model was trained using the training set to verify that the RMSE of the set was the minimum, and all parameters in PSR and ESN were optimized synchronously using IASSA. (5) The optimized optimal parameters were introduced into PSR-ESN, and the test set was ran in order to obtain the prediction results of the subseries. (6) The prediction results of all subcomponents were accumulated to obtain the final wind power prediction results.

Case Analysis
The wind power generation data of the Irish energy system were used for experimental analysis, the wind power generation time was collected from 0:00 on 1 January 2020 to 0:00 on 16 January 2020, and the sampling interval of experimental data was 15 min, yielding a total of 1152 pieces of data.

Description of the Experiment
In order to verify the prediction performance of the proposed TVFEMD-ISSA-PSR-ESN model, the single models BPNN and LSTM, and the combined models SSA-ESN, EMD-SSA-PSR-ESN, and TVFEMD-SSA-PSR-ESN were compared. For the performance evaluation of the prediction models, three indicators were used: root mean-square error (RMSE), mean absolute error (MAE), and mean absolute proportional error (MAPE): For the model involved in the experiment, the B-spline interpolation n and bandwidth threshold ξ for TVFEMD decomposition were set as 26 and 0.3, respectively. The parameters, namely delay time, embedding dimension d, network update weight a, reserve pool connection matrix scaling parameter Rh, reserve pool node number N, input data scaling coefficient W in , and reserve pool sparsity S, were determined using the proposed IASSA and SSA optimization algorithms. Other parameters are presented in Table 4. In the experiment, the training set, validation set, and test set accounted for 60%, 15%, and 25% of all subsequences, respectively.

Data Decomposition and Normalization
After the TVFEMD decomposition of the original wind power data, seven subsequences were obtained and were recorded as IMF1-IMF7, as shown in Figure 4. In order to better train the law of subsequences, all subsequences were normalized separately. The normalization formula is: where y is the normalized data; x is the original wind power data; x max and x min are, respectively, the maximum and minimum values of wind power; y max and y min are, respectively, the maximum and minimum values of the mapping space; and the mapped space is [−1, 1].
Sustainability 2023, 15, x FOR PEER REVIEW Figure 4. Decomposition results of wind power data.

Prediction Results and Analysis
In this subsection, six experimental models are discussed in detail, and the p results are evaluated to demonstrate the validity of the proposed combined m experimental model prediction performance indicators, that is, RMSE, MAE, an including 1-, 3-, and 5-step-ahead predictions, are presented in Table 4, and the er chart is shown in Figure 5.

Prediction Results and Analysis
In this subsection, six experimental models are discussed in detail, and the prediction results are evaluated to demonstrate the validity of the proposed combined model. All experimental model prediction performance indicators, that is, RMSE, MAE, and MASE, including 1-, 3-, and 5-step-ahead predictions, are presented in Table 4, and the error radar chart is shown in Figure 5.
In this subsection, six experimental models are discussed in detail, and the predict results are evaluated to demonstrate the validity of the proposed combined model. experimental model prediction performance indicators, that is, RMSE, MAE, and MA including 1-, 3-, and 5-step-ahead predictions, are presented in Table 4, and the error ra chart is shown in Figure 5. As can be seen from Table 5 and Figure5, in the case of 1, 3, and 5 steps ahead, proposed method predicted the wind power better, and the prediction error gradua increased with the increase in the number of prediction steps. A comparison of the eva ation indicators of different prediction models revealed the following points:  As can be seen from Table 5 and Figure 5, in the case of 1, 3, and 5 steps ahead, the proposed method predicted the wind power better, and the prediction error gradually increased with the increase in the number of prediction steps. A comparison of the evaluation indicators of different prediction models revealed the following points: Thus, the optimization performance of the IASSA not only has a certain effect on the classical test function but also exhibits good capacity to solve the actual parameter optimization problem.

Model 1 Step Ahead 3 Step Ahead 5 Step Ahead
In conclusion, the data analyses presented here highlight the effectiveness of the SSA-ESN model in terms of predictive performance when compared to BPNN and LSTM models. Furthermore, the inclusion of signal decomposition methods, PSR, and the proposed IASSA algorithm demonstrates significant improvements in the models' predictive capabilities. These findings underscore the importance of selecting appropriate techniques and algorithms to optimize data analysis models for the achievement of enhanced performances and better results in real-world applications.
In order to verify the performance of the proposed method more intuitively, the experimental prediction effect was further visualized and analyzed. The error curves, fitting curves, and error distributions of the multistep prediction of the six models are shown in Figure 6, Figure 7, and Figure 8, respectively. To demonstrate the advantages of the proposed model more intuitively, the distribution value of the error was studied; the abscissa of the proposed method was found to be smaller, and its MAE and standard deviation were the smallest-3.2511 and 4.1478 (1 step), 7.5864 and 9.8941 (3 steps), and 11.2086 and 14.6149 (5 steps), respectively. This demonstrated the superiority of the proposed method.

Conclusions
In this paper, we proposed a hybrid model based on TVFEMD, PSR, ESN, and IASS to predict short-term multistep wind power generation. In the proposed model, first, th wind power of the original time series is decomposed into multiple subcomponents b using TVFEMD, which effectively solves the problem of the difficult prediction of no stationary signals. Next, the subcomponents are constructed by PSR to construct the ou put vectors of the prediction model, and the ESN is then used as the final prediction mod to predict the data. Because PSR and ESN involve the adjustment of multiple parameter the IASSA algorithm is used to optimize PSR and ESN parameters synchronousl

Conclusions
In this paper, we proposed a hybrid model based on TVFEMD, PSR, ESN, and IASSA to predict short-term multistep wind power generation. In the proposed model, first, the wind power of the original time series is decomposed into multiple subcomponents by using TVFEMD, which effectively solves the problem of the difficult prediction of nonstationary signals. Next, the subcomponents are constructed by PSR to construct the output vectors of the prediction model, and the ESN is then used as the final prediction model to predict the data. Because PSR and ESN involve the adjustment of multiple parameters, the IASSA algorithm is used to optimize PSR and ESN parameters synchronously. Through a comprehensive experimental analysis of real wind power, the effectiveness of the proposed model is verified.
Based on the experimental results, the following conclusions can be drawn: (1) A comparison of the performance of EMD and TVFEMD decomposition of the original wind power time series revealed that TVFEMD, combined with the prediction model, yielded an improved prediction effect: the difference was expressed via number of steps, namely, 1 step, 3 steps, and 5 steps, and the effect was more obvious. (2) Compared with the original SSA, the proposed IASSA exhibited a better optimization performance in the classical test function and the actual prediction application. The accuracy was effectively improved in the three performance evaluations of root meansquare error (RMSE), mean absolute error (MAE) and mean absolute proportional error (MAPE). (3) Compared with six existing models, the experimental analysis of multistep prediction revealed that the proposed prediction model combining TVFEMD, PSR, ESN, and IASSA synchronous optimization strategies were able to effectively improve the prediction accuracy.
Based on the above three points, the combined prediction model proposed in this paper improved the accuracy of wind power research and provided a new feasible strategy for the multistep prediction of wind power.

Data Availability Statement:
The data used during the study are available from the first author by request.