Short-Term Electricity Price Forecasting Model Using Interval-Valued Autoregressive Process

The uncertainty that dominates in the functioning of the electricity market is of great significance and arises, generally, because of the time imbalance in electricity consumption rates and power plants’ production capacity, as well as the influence of many other factors (weather conditions, fuel costs, power plant operating costs, regulations, etc.). In this paper we try to incorporate this uncertainty in the electricity price forecasting model by applying interval numbers to express the price of electricity, with no intention of exploring influencing factors. This paper represents a hybrid model based on fuzzy C-mean clustering and the interval-valued autoregressive process for forecasting the short-term electricity price. A fuzzy C-mean algorithm was used to create interval time series to be forecasted by the interval autoregressive process. In this way, the efficiency of forecasting is improved because we predict the interval, not the crisp value where the price will be. This approach increases the flexibility of the forecasting model.


Introduction
Over the last few decades, many models for forecasting electricity prices have been developed using different mathematical methods.Generally, these models are based on simulation techniques; time series analysis; autoregressive processes such as autoregressive, autoregressive integrated, autoregressive moving average and seasonal autoregressive methods; artificial neural network; fuzzy logic; grey system theory; game theory; and combinations of two or more methods.Weron explained the complexity of available solutions, their strengths and weaknesses, and the opportunities and threats that the forecasting tools offer or that may be encountered in the process of electricity price forecasting [1].A very useful tutorial review of probabilistic electricity price forecasting that presented much-needed guidelines for the rigorous use of methods, measures and tests, in line with the paradigm of 'maximizing sharpness subject to reliability,' was created by Nowotarski and Weron [2].With the introduction of deregulation to the power industry, the price of electricity has been the key to all activities in the power market.Accurately and efficiently forecasting electricity prices becomes more and more important.Accordingly, Shahidehpour et al. mainly described short-term electricity price forecasting models based on simulations, neural networks, conditional probability distribution, and risk analysis [3].Electricity markets suffer from a severe case of over-dimensionality due to the lack of economic storage of the commodity and the scarcity and hedging formulations can help users to form dynamic replicating portfolios.By using this approach they can reduce the complexity of possible future price combinations [4].Forecasting is an integral part of revenue management.Designing dynamic prices requires forecasts of future demand, and scheduling consumption requires forecasts of future prices.Dutta and Mitra described research on prices in the electricity sector.They showed how different methods, such as the dynamic regression model, the transfer function model, the autoregressive integrated moving average model, the artificial neural network technique, the kernel-based method, principal component analysis, exponential smoothing, Energies 2018, 11, x complexity of possible future pri management.Designing dynamic consumption requires forecasts of electricity sector.They showed ho transfer function model, the autor network technique, the kernel-base ɛ-insensitive loss function support parametric additive model, can be in electricity price forecasting have techniques have been covered: stoc models [6].
In a data-sparse electricity ma to forecast the mid-term electricity electricity market using autoregres al. to forecast today the 24 market Brownian motion has been propos modeling of spot power prices [9 forecast the electricity price withi electricity price forecasting model Partanena was based on the autor wavelet transformation-based neur price profile in a deregulated elec autoregressive integrated moving a to forecast short-term electricity p developed an adaptive short-term the pre-processing of data and a m spike detection were used by Cerja forecasting [15].A method using m for forecasting short-term electricit moving average models, comple heteroscedasticity models and th forecasting of Singapore's weekly extensive empirical study of shor question of whether the optimal multivariate.They provided evide performance, the multivariate mod one across all 12 considered data outperformed by the latter [18].Go attempts to generalize the standard structure of the proposed model functional variables [19].A new a methodology is divided into two integrated moving average mode moving average models per hour -insensitive loss function support vector regression, trigonometric gray prediction, and the semi-parametric additive model, can be used for forecasting the price [5].The main methodologies used in electricity price forecasting have been reviewed by Aggarwal et al.The following price-forecasting techniques have been covered: stochastic time series, causal models, and artificial intelligence-based models [6].
In a data-sparse electricity market, Cheng et al. applied a novel interval grey prediction model to forecast the mid-term electricity market clearing price [7].Day-ahead price forecasting of a daily electricity market using autoregressive integrated moving average models was used by Contreras et al. to forecast today the 24 market clearing prices of tomorrow [8].A modified version of geometric Brownian motion has been proposed by Barlow (2002) as a jump diffusion model for the stochastic modeling of spot power prices [9].Fu and Li applied simulation of power system equipment to forecast the electricity price within the context of a competitive electricity market [10].A hybrid electricity price forecasting model for the Finnish electricity spot market proposed by Voronin and Partanena was based on the autoregressive moving average process and a neural network [11].A wavelet transformation-based neural network model has been applied by Aggarwal et al. to forecast price profile in a deregulated electricity market [12].Ghodsi and Zakerinia made a hybrid of the autoregressive integrated moving average process, an artificial neural network, and fuzzy regression to forecast short-term electricity prices [13].Yamin et al. used an artificial neural network and developed an adaptive short-term electricity price forecasting model [14].Statistical techniques for the pre-processing of data and a multi-layer neural network for forecasting electricity price and price spike detection were used by Cerjan et al. to develop a hybrid model for short-term electricity price forecasting [15].A method using modular feed-forward neural networks and fuzzy inference system for forecasting short-term electricity price was developed by Esfahani [16].Autoregressive integrated moving average models, complemented by the use of generalized autoregressive conditional heteroscedasticity models and their variants to account for volatility, were used for univariate forecasting of Singapore's weekly wholesale electricity prices [17].Ziel and Weron conducted an extensive empirical study of short-term electricity price forecasting to address the long-standing question of whether the optimal model structure for electricity price forecasting is univariate or multivariate.They provided evidence that, despite having the edge in terms of overall predictive performance, the multivariate modeling framework does not uniformly outperform the univariate one across all 12 considered datasets, seasons of the year, or hours of the day, and at times is outperformed by the latter [18].Gonzalez et al. proposed a new functional forecasting method that attempts to generalize the standard seasonal ARMAX time series model to the L 2 Hilbert space.The structure of the proposed model is a linear regression where functional parameters operate on functional variables [19].A new approach to forecast day-ahead electricity market prices, whose methodology is divided into two parts (forecasting of the electricity price through autoregressive integrated moving average models and construction of a portfolio of autoregressive integrated moving average models per hour using stochastic programming) was created by Nieta et al. [20].Neupane et al. employed an ensemble prediction model in which a group of different algorithms participates in forecasting one hour ahead the price for each hour of a day.They proposed two different strategies, namely, the fixed weight method and the varying weight method, for selecting each hour's expert algorithm from the set of participating algorithms [21].
Alvarez et al. applied a partitioning-clustering technique to the electricity price time series to discover the behavior patterns of a series [22].A combination of a fuzzy C-mean clustering algorithm and a teaching-learning-based optimization algorithm was used to forecast the next day's electricity price [23].Sokhanvar et al. proposed a method that involves input-output decomposition and a simple clustering algorithm to classify the data points.The prediction is then generated by using the weighted mean of forecasted outputs of some clusters with the highest probabilities [24].A hybrid model of fuzzy C-elliptotypes of fuzzy clustering and an advanced general radial basis function network was developed to forecast electricity prices in the market [25].
This paper represents a hybrid model based on the fuzzy C-mean clustering and interval-valued autoregressive process for forecasting the short-term electricity price.A fuzzy C-mean algorithm is used to separate a set of monitored data into an adequate number of fuzzy states.Accordingly, we can create chronological fuzzy time series by inspection of what fuzzy state the monitored value belongs to.In many cases, creating fuzzy relations among fuzzy time series is very difficult.To avoid it, we represent a fuzzy time series by an interval time series.Every interval corresponds to an adequate fuzzy state, but is expressed by closed bounded set of real numbers.The lower bound of an interval is equal to the electricity price with a minimum value in the fuzzy state while the upper bound is the price with the maximum value.Generally, a fuzzy C-mean clustering algorithm can be treated as a data preprocessing phase, i.e., the preparation of data to be forecasted by the interval-valued autoregressive process.The main target of the data preprocessing is to decrease the volatility of the monitored states, because sometimes it is very hard to describe it adequately.Calibration of the forecasting model considers two main parameters.The first parameter represents the number of fuzzy states for a monitored electricity price series, while the second gives the order of the autoregressive process.To define the number of fuzzy states we developed an approach based on the autocorrelation function of the interval series created with respect to a given number of fuzzy states.The order of the autoregressive process is defined by Akaike and the Schwartz-Bayesian information criterion.The model was tested on the real monitored electricity price series and the accuracy of the model is very high.

Short-Term Forecasting Model
This paper investigates time invariant fuzzy time series of daily electricity prices.The investigation includes two main components.The first component refers to the fuzzification of monitored data using the fuzzy C-mean clustering algorithm, while the second refers to the creation of a forecasting model based on the interval-valued autoregressive process.

Fuzzification of Monitored Data
The fuzzy time series was originally defined by Song and Chissom [26,27].The basic definitions of the fuzzy time series, and time-variant and time-invariant fuzzy time series, are given as follows.
Suppose that X(t) ∈ R 1 , t = 0, 1, 2, . . . is the universe of monitored data that should be fuzzified by fuzzy sets f i (t), i = 1, 2, . . .and let F(t) be a set of f i (t).Then, F(t) is called a fuzzy time series on X(t).If F(t) is caused by F(t − 1), this dependency is defined as , where R(t, t − 1) is a fuzzy relation of the first order between F(t) and F(t − 1) and * is the fuzzy operator.If F(t) is caused by F(t − 1) and R(t, t − 1) is constant for each t, then F(t) is a time-invariant fuzzy series.
The intention of fuzzification is to transform a series of monitored data into a series of fuzzy states.Using this approach, each value within a series is transferred into an adequate fuzzy state and we obtain the order of states instead of the order of discrete values.For that purpose, we apply the fuzzy C-mean clustering algorithm [28][29][30][31] over the set of daily monitored electricity prices, A fuzzy C-mean algorithm is a method based on the minimization of a generalized least-squared error function within groups.Let P = {p i } R N×1 be the univariate series of monitored data and C m = {c m }, m = 1, 2, . . ., M be the electricity price vector composed of M cluster centers.Every cluster is a fuzzy set defined by the relative closeness to space P = {p i } R N×1 .The membership degree u im ∈ [0, 1] indicates to what degree the relative closeness P N belongs to the cluster center vector C m , which results in a fuzzy partition matrix U = |u im | N×M .Let u im be the membership, C m the center of the cluster, N the number of observed data points, and M the number of clusters.Then, the objective function, which should be minimized in fuzzy clustering, is defined as follows: where d mi is the Euclidean distance between the observation and the center of the cluster; . ., N, and exponent ω is the fuzzy index (ω > 1).In this paper we set up ω = 2.The objective function J is minimized subject to the following constraints: An iterative algorithm is used to find the minimum of the objective function.In j-th iteration, the values of u ω mi and c m are updated as follows: The iteration process stops when J (j+1) − J (j) < δ , where δ represents the minimum amount of improvement or a given threshold value.
The number of clusters represents the number of fuzzy sets or states.According to the obtained fuzzy partition matrix (U = |u im | N×M ) we can define the fuzzy state matrix S(i) of the observed data as follows: If we take into consideration that the fuzzy C-mean algorithm produces a sequence of M cluster centers, then, on the basis of the maximum value of the membership function we can determine the fuzzy state to which the i-th monitored value belongs.Finally, the sequence of fuzzy states S(i) = {s(1), s(2), . . ., s(N)} represents a fuzzy time series on The creation of fuzzy relations among obtained fuzzy states can be very complex.To overcome this problem, we transform the sequence of fuzzy states into an adequate stationary interval series and apply an interval-valued autoregressive process to create a short-term forecasting model.

Interval-Valued Forecasting Model
The process of transformation of a fuzzy time series into an interval time series is based on the fact that each fuzzy state S m (i), i ∈ {1, 2, . . ., N}, m ∈ {1, 2, . . ., M} can be represented by a triplet; , where c min m (i) is equal to the element of the cluster with minimum value, c max m (i) is equal to the element with the maximum value, and c m (i) is the center of the cluster.If we Energies 2018, 11, 1911 5 of 17 take into consideration that we monitored the daily electricity price then transformation produces the following interval time series: where, p l 1 (m) is the lowest and p u 1 (m) is the highest daily electricity price for the m-th fuzzy state.Accordingly, we can say that we created chronological sequence of interval-valued data, i.e., we have a time series of interval data.
An interval variable X is defined as a closed, bounded set of real numbers in the form of X = x l , x u , where x l is the lower bound, x u is the upper bound of interval, and x l , x u ∈ .An interval time series is a sequence of interval variables monitored in successive time points and expressed as a two-dimensional vector In this paper we use the basic interval arithmetic operations introduced by Moore et al. [32]: The attraction of interval arithmetic is that it would not be necessary to analyze whether the conventional point-wise floating-point computations are safe.As the interval results contain all possible values, a narrow interval indicates success.A wide interval does not prove that a conventionally computed result is wrong, but it does indicate a risk.Interval arithmetic has the property of correctness; result intervals are guaranteed to contain the real number that is the value of the expression [33].
Every monitored series contains information about the process that generated it.This process can be realized by modeling the current value as a function of its previous values.An autoregressive process defines the current value of the variable as a weighted linear sum of its previous values and is defined as the following stochastic difference equation: where, β 0 is a constant, β 1 , β 2 , . . ., β k are coefficients of the model, ε t is white noise ( ε t ∼ iidN(0, σ 2 )), and k is the order of the autoregressive process (lag).In the context of interval variables, the autoregressive process becomes an interval autoregressive process: Another representation of the interval can be done in terms of the so-called midpoint and radius, with The midpoint is the center of an interval (the location), whereas the radius is the half-width of an interval (a measure of imprecision), with X R ≥ 0 [34].Let Y(i) = [y C i , y R i ] be the interval response variable and X(j) = [x C j , x R j ] be the set of interval explanatory variables monitored on N units.Accordingly, Equation ( 11) can now be expressed by the following linear relationships: where: the center and radius of the i-th interval response variable, x C i−j , x R i−j -the center and radius of the j-th interval explanatory variables, β C 0 , β R 0 -the center and radius of the constant, β C j , β R j -the center and radius of the regression coefficients of j-th interval explanatory variables, ε C i , ε R i -the center and radius of the residual of the i-th interval response variable.The optimal coefficients are obtained according to the least-squares approach by minimizing the sum of squared residuals of y C i and y R i subject to non-negativity of β R j .The criterion of minimization is defined as: To define the values of the coefficients we apply the center and range method developed by Neto and de Carvalho [35].The values of the coefficients β C 0 , β C 1 , . . ., β C j and β R 0 , β R 1 , . . ., β R j are the solutions of the following system of 2(k + 1) normal equations written in matrix form: where A is a 2(k + 1) × 2(k + 1) matrix and b is a 2(k + 1) × 1 vector described as: .
If the interval time series is stationary, then d = 0; otherwise d > 0, where d denotes the order of differencing.If the interval series is non-stationary, then it is necessary to transform it into a stationary series by d-differencing.First-order differencing is performed in the following manner: After differencing, the interval time series becomes stationary and we can apply the k-th interval autoregressive process (Equation ( 11)) on differenced series.Values of the interval differenced series will be forecasted as follows: Energies 2018, 11, 1911 7 of 17 Equation ( 19) represents values for forecasting the first-order differenced series, but we want to forecast monitored electricity price series.Therefore, we must transform these values into the Ŷ(i) = ŷl i , ŷu i , i = 2, . . ., N form as follows: If we take into consideration that the monitored series is represented by the series composed of M fuzzy states, then it is necessary to find out to which fuzzy state of M values the forecasted state corresponds.We propose the following procedure for expressing the sequence of forecasted fuzzy states by the set of M fuzzy states (Algorithm 1):

Algorithm 1:
Input: forecasted series by AR(k); P(i), i = 1, . . ., N Output: forecasted fuzzy states expressed by M fuzzy states Step 1. create series of forecasted centers; pC i = Step 2. create a fuzzy partition matrix; Step 3. create forecasted fuzzy state matrix; Ŝ(i) Step 4. create the sequence of forecasted fuzzy states; Ŝ(i) = {sˆ(1), sˆ(2), . . ., sˆ(N)} Step 5. according to Ŝ(i) create interval forecasted series; P(i, m) = [ pl i (m), pu i (m)] Algorithm 1 represents the forecasting model.This algorithm is applied to forecast the future electricity price for N + h, h = 1, 2, . . ., H. Note that this is a short-term forecasting model and H should be set to a few days.
Calibration of the forecasting model is related to the determination of two main parameters.The first parameter considers the number of fuzzy states (M-clusters) for a monitored electricity price series, while the second considers the order of the autoregressive process (k).
We developed the following iterative algorithm to define the minimum number of possible fuzzy states (clusters), (Algorithm 2):   ) .( 21) The autocorrelation function (ACF) of any series gives the correlation between y t and y t−k for k = 1, 2, 3, . . . .Theoretically, the autocorrelation between y t and y t−k is calculated as: Since the minimum number of clusters is defined, we continue the iteration process until the condition (Step 7) is met.Suppose the condition is not met in the j-th iteration, in that case, we stop the iterative process and take the number of clusters from j − 1 iteration to be the maximum number of clusters.The algorithm that considers the maximum number of fuzzy states is as follows (Algorithm 3): Since the minimum number of clusters is defined, we continue the iteration process until the condition (Step 7) is met.Suppose the condition is not met in the j-th iteration, in that case, we stop the iterative process and take the number of clusters from j − 1 iteration to be the maximum number of clusters.The algorithm that considers the maximum number of fuzzy states is as follows (Algorithm 3): Figure 2 shows a case when the condition defined by Step 7 is not met.Algorithm 2 produces the minimum number of fuzzy states that can be used to transform a monitored series into an interval one, which can be further forecasted by the interval autoregressive process alone, while Algorithm 3 produces the maximum number of fuzzy states.In this way we can increase the flexibility of the forecasted model because any number of fuzzy states in that range can be taken for the purposes of transformation.Finally, the main aim of the proposed algorithms is to prepare a monitored series to be forecasted by the autoregressive model only.
The order of the autoregressive model (lag length) is defined according to information criteria that measure the quantity of information about the dependent variable contained in a set of process alone, while Algorithm 3 produces the maximum number of fuzzy states.In this way we can increase the flexibility of the forecasted model because any number of fuzzy states in that range can be taken for the purposes of transformation.Finally, the main aim of the proposed algorithms is to prepare a monitored series to be forecasted by the autoregressive model only.
The order of the autoregressive model (lag length) is defined according to information criteria that measure the quantity of information about the dependent variable contained in a set of regressors.In this paper we will use the Akaike information criterion (AIC) [36] and the Schwartz-Bayesian information criterion (SBIC) [37].They are calculated by the following equations: where N is the sample length; k is the total number of estimated coefficients; and ûi are the residuals.
The sum of the squared residuals represents the main part of both criteria and we want to minimize it.Thus, we select the model with the smallest value of ACI or SBIC.
The absolute percentage error of the model for the i-th electricity price is calculated as follows: To demonstrate the effectiveness of the model, the forecasted results are analyzed by the mean absolute percentage error (MAPE):

Numerical Example
In this section, we focus on the practical performance of the proposed approach.To test the efficiency of the model we used the data of daily monitored electricity price recorded by Nord Pool, N2EX Day Ahead Auction Prices [38] (Figure 3).
The set of data was divided into a training subset and validation subset.The training subset was composed of data from 1 March 2017 to 31 March 2018.We used this subset to check the confidence of the developed model.The validity subset was composed of data from 1 April 2018 to 10 April 2018, and has been used to check the validity of the model.
The results of Algorithms 2 and 3 are represented in Figure 4.
Energies 2018, 11, x 9 of 17 The sum of the squared residuals represents the main part of both criteria and we want to minimize it.Thus, we select the model with the smallest value of ACI or SBIC.
The absolute percentage error of the model for the i-th electricity price is calculated as follows: To demonstrate the effectiveness of the model, the forecasted results are analyzed by the mean absolute percentage error (MAPE):

Numerical Example
In this section, we focus on the practical performance of the proposed approach.To test the efficiency of the model we used the data of daily monitored electricity price recorded by Nord Pool, N2EX Day Ahead Auction Prices [38] (Figure 3).
The set of data was divided into a training subset and validation subset.The training subset was composed of data from 1 March 2017 to 31 March 2018.We used this subset to check the confidence of the developed model.The validity subset was composed of data from 1 April 2018 to 10 April 2018, and has been used to check the validity of the model.The results of Algorithms 2 and 3 are represented in Figure 4.    Obviously, six fuzzy states have met the conditions of Algorithm 2 and seven fuzzy states have not met the conditions of Algorithm 3. Accordingly, the minimum number of fuzzy states is six and the maximum number of fuzzy states is six as well.Hence, the expected number of fuzzy states, i.e., clusters, is six.7) we obtained the fuzzy time series on daily monitored electricity price (see Table 1).The transformation of six fuzzy states into six intervals is represented by Table 2.For example, fuzzy state s1 appears 58 times in the monitored sequence.The minimum value of the electricity price for the state s1 is 33.88 EUR/MWh, which corresponds to 16 June 2017, while the maximum value is 44.78 EUR/MWh, which corresponds to 27 May 2017.In this way we defined the first interval that corresponds to fuzzy state s1.The remaining intervals are defined in the same way.The original time series was first transformed into an ordered fuzzy time series whereby each fuzzy state was expressed by an adequate center of cluster.After that, the ordered fuzzy time series was further transferred into an ordered interval-valued time series, where the lower bound of the interval was equal to the electricity price of the fuzzy state with a minimum value and the upper bound was equal to the electricity price of the fuzzy state with maximum value.The process of transformation is expressed by Figure 5.The original time series was first transformed into an ordered fuzzy time series whereby each fuzzy state was expressed by an adequate center of cluster.After that, the ordered fuzzy time series was further transferred into an ordered interval-valued time series, where the lower bound of the interval was equal to the electricity price of the fuzzy state with a minimum value and the upper bound was equal to the electricity price of the fuzzy state with maximum value.The process of transformation is expressed by Figure 5.The stationary nature of the time series is the main prerequisite for the application of an autoregressive process.Figure 6 shows the presence of a trend in the fuzzy time series, i.e., the series is non-stationary.After first-order differencing of fuzzy time series and interval-valued time series, we can see that these two series became stationary series (Figure 7).The stationary nature of the time series is the main prerequisite for the application of an autoregressive process.Figure 6 shows the presence of a trend in the fuzzy time series, i.e., the series is non-stationary.The original time series was first transformed into an ordered fuzzy time series whereby each fuzzy state was expressed by an adequate center of cluster.After that, the ordered fuzzy time series was further transferred into an ordered interval-valued time series, where the lower bound of the interval was equal to the electricity price of the fuzzy state with a minimum value and the upper bound was equal to the electricity price of the fuzzy state with maximum value.The process of transformation is expressed by Figure 5.The stationary nature of the time series is the main prerequisite for the application of an autoregressive process.Figure 6 shows the presence of a trend in the fuzzy time series, i.e., the series is non-stationary.After first-order differencing of fuzzy time series and interval-valued time series, we can see that these two series became stationary series (Figure 7).After first-order differencing of fuzzy time series and interval-valued time series, we can see that these two series became stationary series (Figure 7).After differencing, data were prepared to be modeled by an interval autoregressive process.Hereinafter, we present the interval autoregressive process for k = 1 in detail.According to the fuzzy time series (see  7).The obtained stationary series was then expressed in the form of a midpoint and radius (see Table 3).After differencing, data were prepared to be modeled by an interval autoregressive process.Hereinafter, we present the interval autoregressive process for k = 1 in detail.According to the fuzzy time series (see  5).Since the interval time series was non-stationary, we have done first-order differencing and obtained the following stationary series:  7).The obtained stationary series was then expressed in the form of a midpoint and radius (see Table 3).According to the center and range method, the following matrices ∆A and ∆b are defined as follows:  In the two cases, the model missed the original fuzzy state within a range of ±1 state.The mean absolute percentage error for a monitored electricity price series of 10 days and forecasted fuzzy states is 1.41%.In the two cases, the model missed the original fuzzy state within a range of ±1 state.The mean absolute percentage error for a monitored electricity price series of 10 days and forecasted fuzzy states is 1.41%.In the two cases, the model missed the original fuzzy state within a range of ±1 state.The mean absolute percentage error for a monitored electricity price series of 10 days and forecasted fuzzy states is 1.41%.

Algorithm 2 :Figure 1
Figure1shows a case when the condition defined by Step 7 is met.

Step 5 .Figure 1
Figure1shows a case when the condition defined by Step 7 is met.

Figure 1 .
Figure 1.Autocorrelation function that satisfies the conditions.The autocorrelation function (ACF) of any series gives the correlation between yt and yt−k for k = 1, 2, 3,….Theoretically, the autocorrelation between yt and yt−k is calculated as: =

Figure 1 .
Figure 1.Autocorrelation function that satisfies the conditions.

Algorithm 3 : 2 .
Input: monitored series Y(i), i = 1, 2, . . ., N Output: maximum number of clusters m max Step 1. set up the initial value of M to m = m min and k = 20 Step 2. do the fuzzy C-mean Step 3. create interval series [y l i (m), y u i (m)], i = 1, 2, . . ., N Step 4. if interval series is non-stationary do d-differencing Step 5. create center series y C i , i = 1, . . ., N Step 6. do autocorrelation function for 20 lags ACF(y C i , 20) Step 7. if there is one and only one lag such that ACF(y C i , 20) > ±2 × standard error, set m+1 and go to Step 2. Otherwise stop and get m = m max from the previous iteration Accordingly, the number of possible fuzzy states is defined as M = [m min , m max ].The expected number of clusters equals the center of the interval; E( M) = m min +m max Figure 2 shows a case when the condition defined by Step 7 is not met.Energies 2018, 11, x 8 of 17

Algorithm 3 :
Input: monitored series ( ), = 1,2, … , Output: maximum number of clusters mmax Step 1. set up the initial value of M to m = mmin and k = 20 Step 2. do the fuzzy C-mean Step 3. create interval series ( ), ( ) , = 1,2, … , Step 4. if interval series is non-stationary do d-differencing Step 5. create center series , = 1, … , Step 6. do autocorrelation function for 20 lags ( , 20) Step 7. if there is one and only one lag such that ( , 20) > ±2 × , set m+1 and go to Step 2. Otherwise stop and get = from the previous iteration Accordingly, the number of possible fuzzy states is defined as = [ , ].The expected number of clusters equals the center of the interval; = .

Figure 2 .
Figure 2. Autocorrelation function that does not satisfy the conditions.

Figure 2 .
Figure 2. Autocorrelation function that does not satisfy the conditions.

Figure 3 .
Figure 3. Daily monitored series of electricity price.

Figure 3 .
Figure 3. Daily monitored series of electricity price.

Figure 3 .
Figure 3. Daily monitored series of electricity price.

Figure 5 .
Figure 5. From monitored electricity price time series to interval-valued time series: (a) fuzzy time series; (b) interval-valued time series.

Figure 6 .
Figure 6.Trend of fuzzy time series monitored electricity price time series to interval-valued time series.

Figure 5 .
Figure 5. From monitored electricity price time series to interval-valued time series: (a) fuzzy time series; (b) interval-valued time series.

Figure 5 .
Figure 5. From monitored electricity price time series to interval-valued time series: (a) fuzzy time series; (b) interval-valued time series.

Figure 6 .
Figure 6.Trend of fuzzy time series monitored electricity price time series to interval-valued time series.

Figure 6 .
Figure 6.Trend of fuzzy time series monitored electricity price time series to interval-valued time series.

Figure 11 .
Figure 11.Monitored and forecasted values of electricity price for h = 10 days.

Figure 11 .
Figure 11.Monitored and forecasted values of electricity price for h = 10 days.

Figure 11 .
Figure 11.Monitored and forecasted values of electricity price for h = 10 days.

Table 1 .
Fuzzy time series of monitored series (six fuzzy states).

Table 2 .
From fuzzy state to interval.