1. Introduction
As the core of electricity market operations, accurate electricity price forecasting can help market participators to formulate reasonable competitive strategies and achieve risk minimization as well as benefit maximization [
1]. For power-generation enterprises, according to accurate electricity price prediction decision makers need to assign an appropriate proportion of electricity to participate in day-ahead market transactions, real-time market transactions, auxiliary services market transactions, and bilateral contract transactions. Moreover, they also need to formulate corresponding bidding tactics to obtain maximum return and minimum risks [
2]. For power-supply enterprises and other electricity buyers, based on precise electricity price prediction, they can optimally assign electricity purchasing amounts in the spot market. For electricity regulators, they need to monitor the power market in accordance with accurate electricity price forecasting to assure the stable and ordered development of the power market. Considering various factors affecting electricity prices as well as the uncertain and fluctuant nature of electricity prices, accurate price prediction is becoming an increasingly complicated problem [
3].
The critical role of accurate electricity price prediction for electricity market participants has promoted the appearance of literature on forecasting electricity prices. The techniques used to forecast prices can be divided into two types, which are classical forecasting technologies and intelligent forecasting technologies [
4]. The classical forecasting technologies are made up of some commonly used methods, including auto-regressive integrated moving average (ARIMA) [
5], the dynamic regression (DR) method [
6], as well as the generalized auto-regressive conditional heteroskedastic (GARCH) approach [
7]. The intelligent forecasting technologies employ data-driven methods learning from historical data, containing artificial neural network (ANN) [
8], support vector regression (SVR) based on ARIMA [
9], and fuzzy neural network (FNN) [
10] methods.
Moreover, in view of the nonlinear, fluctuting and inherent uncertaintainty of electricity prices and the complexity of predicting them, a series of hybrid methods has been proposed integrating classical and intelligent prediction techniques, which include the hybrid models combining mutual information (MI), wavelet transform (WT), and FNN on the basis of evolutionary particle swarm optimization (PSO) [
11], MI-WT-SVM (support vector machine) on the basis of the gravitational search algorithm (GSA) [
12], the integrated intelligent approach [
13], and WT-ARIMA [
14]. Among the literature on electricity price forecasting, it has been found that the ARIMA approach is one of the most widely used methods owing to its statistical advantages [
4]. Additionally, Dong et al. [
15] proposed an optimized prediction algorithm separating high fluctuation and daily electricity prices in Australia and employing empirical pattern disintegration method, seasonal arrangement and the ARIMA approach. Mandal et al. [
16] put forward an integrated intelligent method combining a data-screening method according to the WT approach, an optimization algorithm on the basis of the firefly model, and a soft calculating technique using a fuzzy art-map (FA) network to predict day-ahead power priced for the Ontario electricity market. Wu et al. [
17] established a recursive dynamic factor analysis (RDFA) model using an effective subspace tracing approach to track critical components and a Kalman filter approach to trace and predict critical components scores. In [
18], a three-tier front feed neural network model was proposed to predict electricity prices for next-week in the power markets of Spain and California. Abedinia et al. [
19] presented a prediction engine based on a combinatorial neural network (CNN) integrating with a new training method on the basis of a revised chemical reaction optimization method to select optimal weights of CNN for forecasting electricity price data.
Based on the existing research literature in the field of electricity price forecasting, developing more precise prediction approaches for power price is required due to the complex electricity market environment and requirement for electricity price forecasting accuracy. This investigation concentrates on next-day hourly electricity price forecasting that is an electricity price prediction procedure commonly employed by participants in the power market for preparing their bids. The novel dedications of this paper are listed as below:
- (1)
An integrated multi-stage intelligent electricity price forecasting model is proposed combining the Beveridge–Nelson decomposition (B-N-D) model, the least square support vector machine (LSSVM) method, as well as a nature-inspired intelligent optimization approach named the whale optimization algorithm (WOA).
- (2)
In the stage of data processing, the B-N-D approach, for the first time to our knowledge, is used to decompose the original electricity price data series into deterministic component, stochastic item and cyclic trend, taking the fluctuant, random, and uncertain nature of power price data into consideration.
- (3)
In the stage of forecasting, the disintegrated three components will be respectively predicted by employing the LSSVM model, the regularization parameter ‘c’ and the radial basis function (RBF) kernel width ‘σ’ which are optimally selected by the novel optimization method WOA.
The performance of the established integrated multi-stage intelligent electricity price forecasting model will be tested on the day-ahead power market of Pennsylvania–New Jersey–Maryland (PJM) in the United States of America. The forecasting superiority of the proposed model will be verified through comparing it with several models based on B-N-D approach, including single LSSVM, LSSVM optimized by PSO, LSSVM optimized by FOA and ARIMA, and LSSVM optimized by WOA, without disintegrating the initial electricity price data by B-N-D approach.
The other sections of this research are structured as below. Following the introductory section, the basic theories employed in this research are elaborated.
Section 3 depicts the detailed course of the constructed multi-stage intelligent electricity price prediction model.
Section 4 provides the prediction results of the established multi-stage intelligent model.
Section 5 comparatively analyzes the prediction manifestation of the constructed model and another five selected comparison models.
Section 6 obtains the conclusions.
2. Theoretical Framework for Electricity Price Forecasting
The established integrated multi-stage intelligent electricity price forecasting model contains the B-N-D approach, LSSVM method, and an intelligent optimization approach named WOA. This section will elaborate the basic theories of these three approaches.
2.1. Beveridge–Nelson Decomposition (B-N-D) Approach
In order to resolve the first order differential stable data series, Beveridge and Nelson presented the B-N-D approach by which the data series co-integrated at the first order can be disintegrated into transitory component and permanent trend. The permanent trend is composed of a determinacy component and random item. The temporary item has a stable course with zero average data, which is called the cyclic term [
20,
21,
22,
23].
Considering the disintegration qualifications of the B-N-D approach, it is required to demonstrate whether the data series are stable after the first differencing. If they are stable after first order differencing, then the B-N-D approach can be conducted as below.
Suppose the natural logarithmic values of hourly electricity prices to be
in the first order stable form in terms of the Wold theorem:
where
,
Pt represents the power price in time point
t,
μ demonstrates the average value for
in the long-run, ε
t~
(
i.i.d. implies independently and identically distribute),
t indicates time, and λ
i is the coefficient. The expectation data of the above equation can be written as:
where
E is a function employed to calculate the expected data of every variable.
Since ln
Pt represents the natural logarithm data of power price, the first order difference of it means the growth ratio of the power price. Considering Equation (2), the average data of the growth ratio of the power price demonstrates the long-term growth rate. According to the B-N-D theory, the deterministic component
DTt can be disintegrated as:
where
DTt demonstrates the determinacy component at
t time point, and
represents the original logarithm data of the power price.
Moreover, ln
Pt is deemed as the forecasting data in accordance with current information. Morley [
24] held the opinion that such time series can be predicted more precisely via use of a steady uni-variate auto regressive (AR) Equation (1) formula at first order, which can be depicted as:
where
can be calculated through AR Equation (1) formula and
. Based on the Wold form in the AR Equation (1) formula, according to the normality hypothesis, the minimum average squared error (MASE) of period
j ahead forecasting of the first order of Δln
Pt can be written as:
In line with the decomposition identification of B-N-D [
20], the total tendency item of time series represented by
Tt is deemed to be the
MASE forecasting from the long-run level. It is equal to the current level of the sequences and the infinite sum of the
MASE of period
j ahead of the first order prediction:
where
Tt is defined as the total item of data series.
Integrating Equations (5) and (6), based on the infinite sum formula of isometric series (due to
), the total component of ln
Pt after being decomposed by B-N-D under the condition of AR(1) can be calculated by:
The cyclic item of ln
Pt can be calculated as:
where
Ct illustrates the cyclic item at
t time point.
The stochastic item can be computed according to Equations (3), (7), and (8), namely:
where
STt indicates the random term at
t time point.
2.2. Least Square Support Vector Machine (LSSVM) Model
The LSSVM is an improved algorithm of the support vector machine (SVM). Considering that the minimum structural risks by regularization restrict that have built up on the weights of the model, the LSSVM alleviates the problem of convex quadratic programming associated with SVM [
25,
26,
27]. The merits of LSSVM contain setting slack via equality restricts and calculating the regression issues as a series of linear formulas bringing about rapid training as well as more precision and stability [
28]. The basic theoretical framework of LSSVM is elaborated as below.
Set up samples series
, and suppose
as input variables and
as output data of sample
i. Through utilizing the function
, the sample series are mapped into a higher dimensional space from the initial characteristic space, hence it can be written in a linear form:
where
demonstrates the weights vector and
b illustrates the deviation term.
In virginal space, the LSSVM approach with an equality restrict is depicted as:
where
c demonstrates the regularization parameter and
ξi implies the slack component.
The Lagrangian function
L is adopted [
29], and then the qualifications for Karush–Kuhn–Tucker (KKT) for optimality are identified by:
where
ai implies Lagrange multiplier.
Eliminating the series
w and
, the optimization procedure is converted to the linear formula. To be consistent with Mercer’s requirements, the Kernel function is described as:
The LSSVM approach for regression is elaborated as:
Radial basis function (RBF) is selected as the Mercer Kernel function
in this investigation due to the fewer parameters that need to be identified and the advantageous manifestation, which is illustrated as:
In this function, two parameters, namely ‘
c’ and ‘σ’ [
30,
31], are required to be calculated through a complex search procedure. The latest intelligent optimization algorithm WOA is utilized in this paper to select the optimum values of them.
2.3. Whale Optimization Algorithm (WOA)
Prompted by the premium hunting behaviors of humpback whales [
32], the scholar Seyedali Mirjalili established a nature-inspired optimization algorithm called WOA, which is employed to analog the distinctive behaviors of humpback whales [
33]. The hunting behavior explored by Goldbogen et al. discovered that there exist two kinds of plots associated with bubbles, which are ‘upward-spirals’ and ‘double-loops’ [
34]. In the former tactic, the whales dive down several meters and spit bubbles in a twist-form around the hunting objects and move to the surface excerpted from the previous study [
33]. Hence, the mathematical simulation of the twist shape bubble-net hunting tactic is helpful for conducting optimization. The WOA is presented based on such behavior.
The WOA is composed of five steps based on the spiral bubble-net hunting tactic elaborated as:
● Step 1: Parameters installation
The critical parameters of WOA are the number of search individuals SearchAgents_no, the number of investigated variables dim, the maximum amount of iterations Max_iteration, the lower limit lb = [lb1, lb2, ….., lbn−1, lbn], and upper limit ub = [ub1, ub2, ….., ubn−1, ubn] of variables.
● Step 2: Population initialization
Since the WOA is an algorithm based on population, the positions of whales in WOA are set in a matrix as:
where
W is the position matrix of various search individuals,
illustrates the data of
j-th parameter for
i-th search individual,
,
. The parameter
can be obtained through utilizing the random distribution written as the following equation:
where
W(i, j) demonstrates the entry in the
i-th row and
j-th column in matrix
W,
lb(i) and
ub(i) illustrate the bottom and upper limits for
i-th search individual,
rand(i,j) means the random data obtained from the uniform distribution ranging from 0 to 1.
● Step 3: Fitness function identification
During the optimization course, the fitness function
f[*] needs to be identified to evaluate each search individual, and the fitness values of each search individual are stored in the matrix
OW, which are depicted as below:
● Step 4: Iteration procedure
For WOA, the searching individuals regenerate their locations in terms of a randomly chosen searching individual or the best tactic discovered up to now. With the purpose of avoiding local optimum and achieving global exploration, the WOA starts searching with several stochastic tactics, which demonstrates the position of a search individual is renewed with respect to a randomly chosen search individual. The mathematical formula is shown as:
where
demonstrates the stochastic place vector selected from present searching individuals,
means the current position vector,
and
represent coefficient vectors,
t is the generation, | implies calculating the absolute value, and · is the multiplying computing process. The vectors
and
are calculated as:
where
linearly decreases from 2 to 0 in iteration procedure, and
is a stochastic vector in the scope of 0–1.
Values of are decided by the diversification of and the random values of , and the absolute values of determine whether the process is exploration or exploitation. A stochastic individual is selected if , which implies an exploration process. While the best resolution is found if to regenerate the places of the searching individuals, which demonstrates the exploitation process. For , whales swim around the hunting objects in contractible circles and spiral paths.
For the contractible circle path, the searching individuals regenerate their places in accordance with the following formula:
where
implies the place vector storing the best solution discovered up to now. The vector
and
can be obtained in terms of Equations (21) and (22). The contractible circle mechanism can be realized via decreasing the value of
in Equation (21).
For the spiral shape track, it can be imitated as:
where
is the distance of the whale to the hunting object (the best solution found up to now),
b demonstrates a constant employed to identify the path of the logarithmic twist, and
l represents random data generated from the interval [−1, 1].
Considering that the whales move around the optimum solution both in a contractible circle and twist shape path, to simulate such behavior, it is supposed that there exists 50% feasibility to choose the twist path or contractible circle path to regenerate the location of whales in searching optimum solution course. Such modeling process can be mathematically constructed as:
where
p is a random data generated from 0 to 1. If
p < 0.5, the whales move around the optimum solution in a contractible circle. While if
, the whales move around the optimum solution in a twist shape path.
● Step 5: Optimum solution selection.
Through the iteration procedure, the searching individuals regenerate their locations in terms of a random searching individual that is found or the optimum tactic discovered up to now. A stochastic searching individual can be found if , and the searching individuals restores their place based on Equation (20). Otherwise, the optimum solution is selected if , and the searching individuals renew their locations based on Equations (24) and (25) determined by the value of p. Ultimately, if the iteration requirement is met, the iteration of WOA comes to an end.
3. The Course of the Constructed Multi-Stage Intelligent Model for Electricity Price Prediction
In this study, an integrated intelligent algorithm is established to improve the precision of electricity price prediction. Considering the uncertain and fluctuant inherent nature of electricity price data and the interference of stochastic effects on power price data, the B-N-D approach is used to decompose the original hourly electricity price data into determinacy term, cyclic item, and random component. After that, the LSSVM method is applied to forecast three decomposed terms respectively. Finally, through multiplying the prediction results of these three data series, the forecasted electricity price data are obtained. Before forecasting by the LSSVM approach, the optimal values of two parameters ‘σ’ and ‘c’ related to the LSSVM approach need to be selected. This paper employs the latest nature-inspired optimization approach WOA to select the optimal values of those two parameters to enhance the forecasting precision.
The multi-stage forecasting procedures of the electricity price are elaborated as below.
● Step 1: First order stationary examination
In view of the fundamental principles of the B-N-D, it is essential to test whether the first difference form of logarithmic series of the original data is steady. The augmented Dickey–Fuller (ADF) methodology is applied to conduct the first order stability test. Only if the logarithmic form of the initial data is stable after first differencing, the next stage can be stepped.
● Step 2: Original data series disintegration
Since the logarithmic series of original data are stable at first order, the original data can be disintegrated into determinacy trend, cyclic item, and random impact effect in terms of Equations (3), (8) and (9).
● Step 3: Parameters’ setting
In WOA, several parameters need to be installed according to
Section 2.3. Based on the examples in Ref. [
33] and the research object of this paper, we suppose
SearchAgents_no = 50,
dim = 2,
Max_iteration = 100,
lb = 0.01, and
ub = 1000.
● Step 4: Optimization beginning
By utilizing WOA for selecting the optimal values of LSSVM’s two parameters, it is required to determine the fitness function
f[*] firstly. In this investigation, the root mean square error (
RMSE) illustrated as below is employed as
f[*].
where
x(k) demonstrates the historical value of electricity price at time point
k, and
means the calculated value of electricity price at time point
k.
The locations of searching individuals are utilized to indicate the parameters of the LSSVM approach, which are set as and . At the first generation, the fitness values of all searching individuals are calculated employing Equation (27). Then, the population of all searching individuals can be ranked in accordance with the fitness values, and the searching individuals with optimum fitness data are chosen according to the minimum values of RMSE. Since the best searching individual with the minimum RMSE value is discovered, other searching individuals will regenerate their locations related to the best one based on the Equation (26).
● Step 5: Optimization termination
After the first generation, the best searching individual and the optimum fitness value can be generated. Then, the offspring iteration can be conducted, and the searching individuals renew their positions according to the Equations (24)–(26) with respect to the best searching individual selected by the minimum values of RMSE. After 100 iterations, the minimum RMSEs can be discovered, and the optimum data of ‘σ’ and ‘c’ are found. After that, the optimized LSSVM is constructed, and the data of electricity prices are finally calculated and forecasted.
The course of the constructed multi-stage intelligent model of electricity price forecasting is demonstrated as
Figure 1. The pseudo code of the proposed electricity price forecasting model is illustrated in
Figure 2. Through conducting the pseudo code in
Figure 2, we can obtain the optimal parameters of LSSVM and then we can forecast electricity price data via substituting the optimal parameters into LSSVM.
5. Prediction Performance Comparison
With the aim of comparing the prediction performances, five comparison prediction approaches are selected including ARIMA, single LSSVM, LSSVM optimized by fruit-fly optimization algorithm (FOA-LSSVM), and LSSVM optimized by particle swarm optimization (PSO-LSSVM) models based on B-N-D approach, as well as WOA-LSSVM based on the initial electricity price data without B-N-D, testing at three different periods the electricity price data collected from the PJM electricity markets. For the single LSSVM, PSO-LSSVM, and FOA-LSSVM models based on the B-N-D approach, the input data series and output data series are the same as the established WOA-LSSVM based on B-N-D approach, while the diversity among these models is that the mechanism of parameters’ determination. For the ARIMA method, when the deterministic term, periodic term and random item are predicted, only the historical data of component itself are used to predict the corresponding term without considering the data of the last moment and the data at the same moment of yesterday. For the WOA-LSSVM based on the initial electricity price data without the B-N-D, the initial electricity price data are employed to predict future data of hourly electricity price, for which the electricity price data at last moment and at the same moment of yesterday are taken as input data series.
Based on the electricity price data series mentioned in
Section 4, for the single LSSVM approach based on B-N-D (B-N-D-LSSVM), the corresponding values of the determinacy trend, periodic item and stochastic effect at the last moment and at the same moment of yesterday from 1 to 168 h are treated as the input data series to predict three items respectively. Then, through learning the historical data, the optimized values of ‘
σ2‘ and ‘
c’ are found [
36]. Finally, hourly electricity price data in the next 24 h can be predicted.
For the FOA-LSSVM model based on B-N-D (B-N-D-FOA-LSSVM), the parameters ‘σ2‘ and ‘c’ in LSSVM approach are used to respectively predict the deterministic trend, periodic item and stochastic effect, and are determined by FOA through optimization iteration. Similar to WOA, there are also several parameters that need to be initialized in FOA, which are the maximum iteration set as 100, and the amount of searching individuals supposed to be 50, (X_axis, Y_axis) ⊂ [−50, 50], FR ⊂ [−10, 10], X_axis = rands(1, 2), and Y_axis = rands(1, 2). Then, the optimal value of ‘σ2‘ and ‘c’ can be discovered, and future electricity price data can be forecast.
For the PSO-LSSVM method based on B-N-D (B-N-D-PSO-LSSVM), parameters ‘σ2‘ and ‘c’ in the LSSVM approach used to respectively predict determinacy trend, cyclic item and stochastic effect are selected by PSO through iteration. It is also necessary to initialize several parameters in PSO, including max_iteration = 100, swarm size = 30, particle size = 2, the minimum of the particle= [0.01, 100], the maximum of the particle= [0.1, 1000], the lower limit of velocity = −500, the upper limit of velocity = 500, and learning factors c1 = 1.5, c2 = 1.7. Then, the best values of ‘σ2‘ and ‘c’ can be explored by PSO, and the electricity price data can be forecast.
The ARIMA based on B-N-D (B-N-D-ARIMA) has been widely used in previous literature related to electricity price forecasting [
4,
9]. The ARIMA method can be utilized to forecast the non-stationary time series. It regards time series as a stochastic process which can be simulated and depicted by a mathematical model. After the mathematical model is identified, the future data can be predicted by the past and present time series data. This method takes the dynamic and persistent characteristics of time series into account, which can reveal the interaction between the past and present time series and between the present and the future time series. Considering the deterministic term, cyclic component and random trend, the ARIMA method is used to simulate the variation tendency of the past 168 h of these three components, and then the future 24 h of the corresponding terms will be forecast. Therefore, the future 24 h power price data can be forecast by multiplying the forecasted three components.
According to the principles of aforementioned models, the future 24 h electricity price data can be predicted. The parameters values of B-N-D-LSSVM, B-N-D-FOA-LSSVM, B-N-D-PSO-LSSVM, and WOA-LSSVM are illustrated in
Table 4. The forecasting results of B-N-D-LSSVM, B-N-D-FOA-LSSVM, B-N-D-PSO-LSSVM, B-N-D-WOA-LSSVM, B-N-D-ARIMA and WOA-LSSVM are illustrated in
Figure 5. It can clearly be seen that the gap between B-N-D-WOA-LSSVM and actual electricity price data is the smallest, followed by B-N-D-PSO-LSSVM. While the gap between WOA-LSSVM and actual electricity price data is the largest, followed by B-N-D-ARIMA.
To further analyze the forecasting performances of all models,
MAPE and
RMSE are applied according to Equations (27) and (29). The calculated data of
MAPE and
RMSE are demonstrated in
Table 5. It can be seen that the
MAPE and
RMSE values of B-N-D-WOA-LSSVM are the smallest, respectively 0.54% and 0.18, while the
MAPE values of B-N-D-PSO-LSSVM, B-N-D-FOA-LSSVM, and B-N-D-LSSVM are respectively 1.44%, 2.76%, and 3.07%. Therefore, it can be concluded that the optimization algorithm WOA can effectively select the optimal data of two parameters in LSSVM, which can greatly increase the prediction accuracy of LSSVM in electricity price forecasting. By comparing the
MAPE and
RMSE values of B-N-D-LSSVM and B-N-D-ARIMA, it can be found that the LSSVM method is more accurate than ARIMA method in electricity price forecasting. By comparatively analyzing the
MAPE and
RMSE of B-N-D-WOA-LSSVM and WOA-LSSVM, it is apparent that the B-N-D approach can significantly improve the forecasting accuracy by 5.23%. Since the B-N-D method can disintegrate the electricity price data into the determinacy trend, cyclic item and random effect, the uncertainty, stochastic, and fluctuation instincts of electricity price data can comprehensively be considered when these three components are forecast, and hence the forecasting precision can be improved greatly. By comparing the B-N-D-WOA-LSSVM with B-N-D-PSO-LSSVM, B-N-D-FOA-LSSVM, B-N-D-LSSVM, B-N-D-ARIMA and WOA-LSSVM, it is confirmed that the proposed multi-stage intelligent model employing LSSVM optimized by WOA based on the B-N-D approach can effectively enhance the prediction precision of electricity prices.
In order to validate the forecasting performance of B-N-D-WOA-LSSVM, we also selected another 8 days’ hourly electricity price data from 15 January 2018 to 22 January 2018 of the PJM electricity market [
35]. As with the process conducted in the first case, 168 h electricity price data are selected to be the training sample, and 24 h electricity price data are treated as the testing sample. By comparing the prediction performances of B-N-D-WOA-LSSVM, B-N-D-PSO-LSSVM, B-N-D-FOA-LSSVM, B-N-D-LSSVM, B-N-D-ARIMA, and WOA-LSSVM employing the calculated data of
MAPE and
RMSE shown in
Table 6, we can obtain the prices. Although the forecasting accuracy of all models decreased, the proposed B-N-D-WOA-LSSVM method still performs best because its
MAPE value is 0.80%, lower than comparative models, followed by B-N-D-PSO-LSSVM with 1.76% and B-N-D-FOA-LSSVM with 2.52%. The performance of WOA-LSSVM without B-N-D performs worst with 6.72%
MAPE value.
With the purpose of eliminating the effect of the length of selected data on prediction accuracy, another 25 days’ hourly electricity price data from 1 February 2018 to 25 February 2018 of we can obtain PJM electricity markets [
35] are also selected to verify the superiority of the proposed B-N-D-WOA-LSSVM for electricity price forecasting. As with the process conducted above, 480 h electricity price data are treated as the training sample, and 120 h electricity price data are examined as the testing sample. After the prediction performances of B-N-D-WOA-LSSVM, B-N-D-PSO-LSSVM, B-N-D-FOA-LSSVM, B-N-D-LSSVM, B-N-D-ARIMA, and WOA-LSSVM are compared by employing the
MAPE and
RMSE criteria shown in
Table 7, it can be seen that the forecasting precision has generally been improved compared with that of 22 January 2018. The B-N-D-WOA-LSSVM still obtains the best performance among all methods, the
MAPE value of which is 0.71%, much lower than other comparative approaches, followed by B-N-D-PSO-LSSVM and B-N-D-FOA-LSSVM. The prediction precision of WOA-LSSVM without the B-N-D procedure ranks last with 6.41%
MAPE value.
According to the forecasting results of three selected cases, it can be verified that the proposed forecasting model B-N-D-WOA-LSSVM is appropriate for electricity price prediction and can enhance the prediction precision significantly.