Electrical Load Prediction Using Interval Type-2 Atanassov Intuitionist Fuzzy System: Gravitational Search Algorithm Tuning Approach

: Establishing accurate electrical load prediction is vital for pricing and power system management. However, the unpredictable behavior of private and industrial users results in uncertainty in these power systems. Furthermore, the utilization of renewable energy sources, which are often variable in their production rates, also increases the complexity making predictions even more difﬁcult. In this paper an interval type-2 intuitionist fuzzy logic system whose parameters are trained in a hybrid fashion using gravitational search algorithms with the ridge least square algorithm is presented for short-term prediction of electrical loading. Simulation results are provided to compare the performance of the proposed approach with that of state-of-the-art electrical load prediction algorithms for Poland, and ﬁve regions of Australia. The simulation results demonstrate the superior performance of the proposed approach over seven different current state-of-the-art prediction algorithms in the literature, namely: SVR, ANN, ELM, EEMD-ELM-GOA, EEMD-ELM-DA, EEMD-ELM-PSO and EEMD-ELM-GWO. proposed system parameter tuning approach on previous identiﬁcation benchmark problems. The comparisons with several other methods including Type-1 TSK FNS, Type-2 TSK FNS, Feedorward Type-2 FNN, SIT2FNN, SEIT2 FNN, TSCIT2FNN, IT2 FNN-GD, IT2 FNN-SMC, IT2 FNNPSO+ SMC, IT2 IFLS -DEKF+GD, and IT2FLS with Modiﬁed SVR support the idea that the proposed approach is an efﬁcient approach in system identiﬁcation problems. The proposed approach is then investigated on electrical load prediction for ﬁve regions of Australia and Poland in the presence of noise and uncertainty which inherently exist in these datasets. Statistical properties of these datasets are presented that show they can be either stationary, trend stationary or difference stationary. In the case of the Poland electrical load dataset, the inputs to the IT2AIFLS are considered to be current values as well as three consecutive time delays of data. In the case of the ﬁve regions in Australia, current data values as well as four consecutive lags are used as the inputs to the fuzzy logic system. In both cases, one-step ahead prediction is considered. For the Poland dataset, the obtained prediction results are compared with several other algorithms including IT2FLS DEKF+GD, IFLS DEKF+GD and IT2IFLS DEKF+GD. The comparisons made in this paper show that the proposed algorithm results in superior performance when it is compared to these methods. For the case of the ﬁve regions of Australia, it is observed that the proposed prediction can perform with a much higher performance as compared with SVR, ANN, ELM, EEMD-ELM-GOA, EEMD-ELM-DA, EEMD-ELM-PSO and EEMD-ELM-GWO methods of prediction. Hence, the proposed approach is an automated solution to develop the antecedent part parameters of an IT2AIFLS, which successfully outperforms current state-of-the-art approaches in literature when applied to an electrical load prediction problem. A possible


Introduction
Electrical power is vital for our life, as it has illuminated our living areas, workplaces, and planet. Electrical load consumption in a region varies with several parameters such as population density, wealth, social factors, climate and distribution of use between home and industry [1]. It is required to have an accurate prediction of electrical load to continue to generate electrical power to meet demand and decrease the cost of electrical power delivery. In the UK, a study has shown that one percent load prediction error is equivalent to 10 million GBP in operating costs per year for the UK power system [2,3]. Electricity generation companies use load forecasting data from a couple of hours to a week ahead to produce energy at what is perceived to be the correct volume and to plan maintenance [4]. Traditional flat methods of pricing electrical load are not efficient as they result in high peaks corresponding to people's social behavior and the habits of industry [4]. To avoid such peaks, manage the network and reduce bills, newer methods such as variable peak pricing tariff and real time pricing are proposed which heavily depend upon realtime load prediction. The output of the prediction algorithms is further utilized in real time pricing, decision making processes and control algorithms. Therefore, for power system management purposes, improving the performance of predictors is of high interest. Renewable energy sources such as solar energy [5], wind power [6] and offshore energy [7] further increase the complexity and make predictions even more difficult due to their natural power fluctuations. Soft sensors technology [8] can estimate solar irradiation [5] which may contribute to decreasing the uncertainty introduced by this source of energy to the overall power system. However, along with reducing the uncertainty associated with renewable energy sources, usage of more powerful prediction approaches to deal with these uncertainties is highly desirable.
A variety of machine learning techniques are already used for electrical load prediction. Among classical approaches, the autoregressive integrated moving average method is most commonly used in the literature [9,10]. However, more recent approaches to electricity load prediction have explored the use of fuzzy logic due to its power to describe the real world in terms of IF-THEN statements as well as its learning capabilities. Among fuzzy logic approaches, interval type-2 fuzzy logic systems used in [1] and interval type-2 Atanassov fuzzy logic systems (IT2AIFLS) used in [11] are previous approaches which have been used to predict electrical load data in different regions of Australia and Poland, respectively. The IT2AIFLSs benefit from more degrees of freedom compared to interval type-2 fuzzy logic systems, as in addition to membership grades they benefit from non-membership grades. Based on previous studies, the greater degrees of freedom in this type of fuzzy system may improve its overall performance [11,12]. However, the success of such fuzzy logic systems is highly dependent on the methods used to solve the complex task of training them. This motivated us to challenge the use of new optimization methods to increase the performance of IT2AIFLSs.
Different methods currently exist to train IT2AIFLSs including gradient descent approaches as well as hybrid computational methods [11,12]. A full gradient descent method for training this structure is investigated in [13], where statistical analysis supported the null hypothesis (α = 0.05) that the IT2AIFLS outperforms type-1 Attanassov intuitionist fuzzy logic systems. Although gradient descent methods provide a computational approach to train IT2AIFLSs, they suffer from falling in a local minima and instability problems.
Hybrid computational methods instead use different training methods in the antecedent and consequent parts of an IT2AIFLS to improve performance. The hybrid training algorithms investigated in [11] use gradient descent for the antecedent part of the IT2AIFLS and Kalman filter for its consequent part. This prediction model is then used for prediction of the Mackey-Glass chaotic time series, the Santa Fe time series and the Box-Jenkins time series, in addition to the electrical load consumption of Poland in the 1990s [11]. This approach has also been applied for the prediction of energy, stock market and house price datasets [12]. Hybrid training approaches utilizing intelligent optimization algorithms previously implemented on IT2FLSs include particle swarm optimization plus Kalman filter [14], particle swarm optimization plus recursive least square [14] and particle swarm optimization plus gradient descent [15]. However, to the best of the authors' knowledge a hybrid intelligent optimization approach for the antecedent part of an IT2AIFLS and a computational method for its conclusion part parameters has not been investigated in the literature. Hybrid training methods offer better performance than solely using a computational method as they do not have local minima and instability issues. Since intelligent optimization approaches are more appropriate options for the antecedent part parameters which appear nonlinearly in the output of an IT2AIFLS, they are more appropriate options for training them.
In this paper, a new hybrid training approach is proposed to improve the performance of the hybrid training methods. The gravitational search algorithm (GSA) is chosen for optimizing the antecedent part parameters of the IT2FLS. The GSA is a physics inspired optimization algorithm [16] which defines each solution in terms of an object with mass, position, velocity and acceleration. The mass assigned to each object is proportional to its cost function and gravitational force among various objects absorbs masses with worst cost function to better object while scanning the space between them to find the optimum solution to the problem [17,18]. The GSA was chosen due to its high performance, benefiting from multiple solutions and less probability to fall in a local minimum. This algorithm has already outperformed several optimization algorithms such as particle swarm optimization, the real genetic algorithm, the differential evolution algorithm and central force optimization in benchmark optimization problems [18][19][20]. Additionally, the ridge least square (R-LS) algorithm is chosen for the consequent part parameter which solves an l 2 cost function as a summation over identification error square plus the l 2 -norm of the weights of the prediction system [21]. The R-LS algorithm is chosen as it is widely known to have superior generalization performance than the simple least square algorithm [21]. Comparisons are made between the proposed approach and some of the existing approaches in the field including hybrid gradient descent Kalman filter training for an IT2AIFLS [11] for electrical load prediction in Poland. Another comparison is provided between the proposed estimation method and seven different estimation algorithms for electrical load prediction in five regions of Australia, namely: SVR, ANN, ELM, EEMD-ELM-GOA, EEMD-ELM-DA, EEMD-ELM-PSO and EEMD-ELM-GWO. The comparison results show the superior performance of the proposed approach over state-of-the-art approaches in the literature. In particular, it outperforms the approach investigated in [11] which uses the IT2AIFLS but hybrid gradient descent Kalman filter as the training algorithm.
This paper is organized as follows. Section 2 introduces the general structure of interval type-2. The methodology of this paper is presented in Section 3. Simulation results are presented in Section 4. Finally concluding marks are presented in Section 5.

Attanassov Intuitionist Fuzzy System
Atanassov intuitionist fuzzy logic systems are a newer variant of the fuzzy logic family which have been successfully applied for prediction purposes [11] as well as pattern recognition [22]. Non-membership grade (ν) for an input to an ordinary fuzzy MF (µ) is simply calculated as the complement of its membership grade as (ν = 1 − µ). However, for an intuitionist fuzzy set (ν + µ) does not necessarily need to be equal to one [23]. This introduces a degree of hesitation or intuition for the fuzzy set. Atanassov in 1986 defined π ∈ [0, 1] as the degree of hesitation which complements the membership and non-membership grades of an input such that π + ν + µ = 1.

Attanassov Intutionist Fuzzy Set
A fuzzy set is an ordinary fuzzy set if the degrees of membership and non-membership for every single input value add up to one. However, if the degrees of membership and non-membership do not add up to one for some input values, the fuzzy set is an intuitionist one. Let X be the universe of discourse and x be an individual value selected from it. The intuitionist fuzzy set A * is defined as presented in the following equation: where µ A * (x) : X → [0, 1] is the membership grade and ν A * (x) : X → [0, 1] is the nonmembership grade for the input value x ∈ X, and we have 0 ≤ µ A * (x) + ν A * (x) ≤ 1 [24]. In the special case, when ν A * (x) = 1 − µ A * (x) the intutionist fuzzy set reduces to an ordinary fuzzy set. However, if this equality does not hold the intuition index of X in A is represented by π A * (x) and defined by:

Structure of Attanassov Intutionist Fuzzy System
Let T be the total number of inputs for the IT2AIFLS, with each sample containing n−dimensional input values x ∈ R n and m−dimensional output values y ∈ R m . The membership functions considered for this structure are interval type-2 Gaussian MFs with uncertain σ values. The inference mechanism to calculate the output of the IT2AIFLS is demonstrated in Figure 1 and is explained as follows: Layer 1: The input layer is the first layer of this system which consists of n nodes passing input values to the fuzzification layer.
Layer 2: The fuzzification layer is the second layer which uses interval type-2 Atanassov membership functions. The inputs to this layer are the outputs of the previous layer and its outputs are the degrees of membership and non-membership which are themselves interval values. The π c,ik (x i ) represents the IF-index or hesitation of center and π var,ik (x i ) is the IF-index of variance and are defined by: such that 0 ≤ π_c (x)≤ 1, 0 ≤ _π_var (x) ≤ 1 and 0 ≤ π_var (x) ≤ 1 and the degrees of membership as well as non-membership are defined as follows: The parameters σ 2i,k , σ 1i,k , π c,ik , π var,ik and c are premise part parameters associated with interval type-2 intuitionistic fuzzy MFs. Furthermore, n is the number of inputs to the system and M is the total number of rules in the fuzzy system. Layer 3: The rule layer is the third layer of the IT2AIFLS which calculates the firing values of the rules of the fuzzy system are presented by [12]: Layer 4: The output layer is the last layer of the system which performs the defuzzification + type reduction and calculates the output of the fuzzy system as follows.
where the parameter β ∈ [0, 1] is the coefficient which determines the weight of its corresponding terms in the output. Moreover, the firing of the rules corresponding to the membership functions are defined as follows: and the ones corresponding to non-membership values are defined as follows: Furthermore, F µ k and F ν k are defined as follows: The output of the IT2AIFLS can be rewritten as follows [25]: The Equation (10) can be further compacted to the following form.
To be able to apply ridge least square algorithm to update the parameters of the IT2AIFLS, we need to write it as follows: where the parameters φ µ and θ µ are as follows: and the parameters φ ν and θ ν are as follows:

Methodology
The parameters of the IT2AIFLS were estimated using the GSA and the R-LS algorithm. The GSA was used to estimate the antecedent part parameters of the IT2AIFLS which appears nonlinearly in the output of the system. The R-LS algorithm was used to estimate the consequent part parameters of the IT2AIFLS as they appeared linearly in the output. Because the GSA is initiated from multiple start points, the probability of falling into a local minimum for this algorithm is less than computational approaches, such as gradient descent, that are initiated from a single point.

GSA Optimization of Antecdent Part Parameters
The GSA was employed to optimize the antecedent part parameters which include the interval associated with sigma value σ 1,ik σ 2,ik and the crisp center (c ik ) of intutionist Gaussian functions in (4). Using the GSA, the antecedent part parameters were encoded in terms of the positions of particles in the GSA, with their velocity vector and acceleration terms updated using the GSA. Each solution in the GSA is defined as the particle positions in a d−dimensional search space representing the antecedent part parameters. The position vector associated with the GSA is as follows: where σ l 1,ik , σ l 2,ik and c ik refer to antecedent part parameters as appeared in (4) which correspond to the k-th rule for i-th input and l represents the l-th solution for the particle, σ l 2,ik refers to σ 2 value as appeared in (4) corresponding to the k-th rule for i-th input and l is the solution counter. The total number of antecedent part parameters to be estimated in this case is equal to 3 × n × M. To evolve the antecedent part parameters according to the GSA, a mass value was assigned to each particle according to its merit, with a higher mass representing an antecedent part with a lower mean squared error and a lower mass representing a higher mean squared error. This makes the particles with worse performance move towards better particles. The mass of particles was updated and normalized at the t-th iteration as [5]: where m i (t) is a non-normalized mass value corresponding to the i-th particle at iteration number t. The values of m i (t) represent the quality of a solution and are defined as [5]: where f p l is the mean squared error value corresponding to p l after estimating the consequent part parameters corresponding to p l parameters using the R-LS algorithm. The R-LS is summarized in Section 3.3, f worst (t) and f best (t) are updated at every iteration as follows: In each iteration, the acceleration term and velocity value must be calculated to update the position vector p l . The position term in each step was updated as follows: where v l represents the velocity vector and was updated using the acceleration term a l as follows: where v l (t) ∈ R d represents the d−dimensional velocity of the particles at t th iteration and r i ∈ [0, 1] is a uniform random number. The acceleration term, a l (t), was then updated as follows: where k b represents the number of best solutions selected, . stands for the Euclidean norm, ε is a small value added to prevent division by zero, r p is the power considered for the Euclidean distance between two particles, G(t) is the gravitational constant and r j ∈ [0, 1] is a uniform random value. The gravitational constant was then updated at each iteration using the following equation: where G 0 has a constant real value and t max is the maximum value of algorithm iteration.

Ridge Least Square Estimation of Consequent Part Parameters
The consequent part parameters of the IT2AIFLS were trained using the R-LS algorithm, also known as Tikhonov regularization. To be able to apply the R-LS algorithm to the consequent part parameters we needed to write them in a matrix form as in (23). Then the R-LS algorithm of (27) could be applied to estimate the consequent part parameters. The R-LS algorithm, also known as Tikhonov regularization, is the solution to the following l 2 -norm cost function [21]. min where . 2 denotes the l 2 -norm and λ is a regulation parameter. Furthermore, the parameters Φ, θ and Y are as follows: where φ µ and φ ν for each sample are defined in (13) and (14), S is the total number of samples, θ µ and θ ν are defined in (13) and (14), y s , s = 1, . . . , S are the measured values. This solution to this problem prevents over-fitness and the analytical solution to the problem was obtained by: where I is the identity matrix. It was then necessary to select the parameter λ for the R-LS algorithm. Large values for λ needed to be avoided to maintain the prediction accuracy of the system. An appropriate value for this parameter was found to be 0.01 in this paper through trial and error.

Performance Measurement
Root mean squared error (RMSE), as the most common performance evaluation method, was used to illustrate the identification performance with its mathematical formula being as follows.
where N is the number of test samples, y actual (t) is the actual value of the system target and y IT2AIFLS (t) is the output of the IT2AIFLS. The flow chart of the proposed tuning method for the IT2AIFLS using GSA for the antecedent part parameters and the R-LS algorithm for the consequent part parameters is illustrated in Figure 2. The first step was to create the random initial population, then the consequent part parameters were estimated using R-LS according to (28). The output of the IT2AIFLS can then be calculated using Equations (2)- (14). The cost function associated with each member of the population was calculated using the RMSE Equation of (29). The GSA algorithm was then iterated for a single iteration to generate the next positions according to (16)- (22). If the stop conditions were not met, we had to iterate the algorithm until the stop condition was satisfied.

Benchmark Identification Problem
Although the proposed prediction method is mainly designed for electricity load prediction, to test its efficacy, it is implemented on a benchmark second order nonlinear dynamic system with time-varying parameters [26]. This nonlinear dynamic system has previously been tested in several papers with other prediction methods to show their efficacy. The output of this dynamic system is a nonlinear time-varying function of inputs, with time delays of input and output as follows [26].
where the nonlinear function f (.) is defined as follows: and parameters a, b and c are time-varying parameters defined by: with T, the total number of samples, taken to be equal to 1000. The input signal to the system u(t) is taken as follows: where: The first 80% of the generated data is used for training and the last 20% is chosen for testing purposes. The comparison results with several other methods including Type-1 TSK FNS [27], Type-2 TSK FNS [27], Feedorward Type-2 FNN [11], SIT2FNN [28], SEIT2 FNN [29], TSCIT2FNN [30], IT2 FNN-GD [26], IT2 FNN-SMC [26], IT2 FNNPSO+ SMC [26], IT2 IFLS -DEKF+GD [11], and IT2FLS with Modified SVR [31] are presented in Table 1,where results support the idea that the proposed approach is effective at system identification by outperforming the other tested algorithms. The behavior of the proposed identification system for the training and tesing data is presented in Figures 3 and 4, respectively. As these figures show, the error between the real data and the output of the IT2AIFLS is very low.
Motivated by the fact that the proposed system identification method can successfully outperform several state-of-the-art system identification approaches in the literature, this method was then used for electrical load prediction in Sections 4.2 and 4.3, where the real-time electrical load of Poland and five different regions in Australia are considered, respectively. Comparisons with state-of-the-art prediction models are then presented to illustrate the performance of the proposed prediction method.

Electrical Load Prediction for Poland
The dataset selected in this part shows the performance of the proposed approach to deal with electrical load prediction using the Poland electricity load dataset available online [32] which presents electricity load values of Poland in the 1990s on a daily basis [33]. A statistical analysis is conducted for the Poland dataset to examine if this time series is stationary or not. To evaluate this property of the Poland electrical load dataset, Augmented Dickey Fuller (ADF) and Kwiatkowski-Phillips-Schmidt-Shin (KPSS) tests are utilized. These two tests are available under the statsmodels Python package. The test statistic value under these two tests as well as their critical vales are given in Table 2 where it can be seen that the test value for ADF is greater than the critical value with a 5% confidence interval, and the test value for KPSS is less than the critical value with a 5% confidence interval. This means that this type of dataset is a trend stationary one. Although the stationary property of the signal is required in some statistical approaches, it is expected that the IT2AIFLS can handle this trend stationary signal due to its nonlinear nature. As separated on the website, 1400 sample data are selected for training and 201 data samples are selected for testing. In this approach the one-step ahead prediction problem is investigated. The inputs taken for this prediction are the current value of electricity load and its time delays as [y(t), y(t − 1) y(t − 2) y(t − 3)]. Here, 14 rules (number of r µ k s and r ν k s) are considered for the fuzzy system and the number of rules was obtained by trial and error to maximize performance. The performance comparison between the proposed approach and previous approaches in the literature are presented in Table 3. As can be seen from the table, IT2AIFLS with GSA-R-LS outperforms IT2FLS DEFKF+GD, IFLS DEKF+GD, IT2 IFLS DEKF+GD and IT2 AIFLS DEKF+GD (previously studied in [11] on the same dataset) by at least 5.7%. The prediction performance of the proposed algorithm for the training and testing data is presented in Figures 5 and 6, respectively. As can be seen from these figures, the prediction output of the IT2AIFLS replicates the real data with high performance. The enlarged portion of the plot given in Figure 5 shows that the prediction output of the IT2AIFLS closely follows not only the low frequency part of the measured data but also the high frequency oscillations of the measured data present in the train dataset. Table 3. Performance comparison of the proposed hybrid algorithm for the IT2AIFLS on the Poland electrical load dataset in terms of RMSEs. The best results are highlighted in bold.

Model Train/Test RMSE Train RMSE Test
IT2FLS DEKF+GD [11] 1395/196 0.0564 0.0595 IFLS DEKF+GD [11] 1395/196 0.0589 0.0599 IT2 IFLS DEKF+GD [11] 1395/196 0.   Figure 7 shows the regression analysis for the IT2AIFLS on the Poland electrical load dataset for test data. It can be seen that the predicted value-target value graph is close to the ideal graph of (Y = T). Furthermore, the R value for this prediction is equal to 0.94 which is close to one. Hence, the prediction is performed with high quality.

Electrical Load Prediction for Five regions in Australia
In this section, the proposed hybrid training method is used for the prediction of  Table 4 presents the statistical parameters associated with data samples which show that this dataset has high levels of variation with a large standard deviation. Furthermore, there exist some cases in Table 4 for which the Skewness value is larger than 0.5 which means that data are moderately right skewed in these cases. Example moderately rightskewed cases are the train data for QLD, the train data for SA, and the train and test data for VIC. Furthermore, the skewness for the test data in the case of the NSW region is −0.69 which is less than −0.5. This means that test data in the NSW region is moderately left-skewed. In the remaining cases, the absolute value of skewness belongs to the interval of [−0.5,0.5] which means that those datasets are fairly symmetrical. Table 5 presents the results of the ADF and KPSS tests for the electrical load datasets of five regions of Australia. As Table 5 shows, the datasets associated with the electrical load for the NWS, SA, TAS and VIC regions of Australia are stationary; however, the dataset associated with the QLD region is difference stationary. Usually, it is easier to deal with stationary signals. Despite this, it is expected that the IT2AIFLS can handle difference stationary signals, where the difference of the signal is stationary, due to its nonlinear nature. The comparison results of the proposed training algorithm for the IT2AIFLS with seven different approaches using the RMSE performance index on the test data are presented in Table 6. As Table 6 shows, the IT2AIFLS with the proposed training method outperforms seven other previously studied algorithms in [34], namely: SVR, ANN, ELM, EEMD-ELM-GOA, EEMD-PSO-ELM, and EEMD-ELM-GWO. The percentage improvement over these algorithms is significant and is at least 67% of improvement. Figure 8 presents plots of the prediction output of the proposed algorithm against the measured data for the electrical load data in each of the five regions of Australia. These plots show that the prediction output of the proposed algorithm closely follows the measured data which indicated that the performance of the proposed algorithm is satisfactory.  To provide greater insight on the obtained results, regression analysis is performed on the results. Figure 9 provides regression analysis for all five regions of Australia for the test data. As Figure 9 shows, the R value in all cases is very close to 1 and the bias term is reasonable. This indicates that the electrical load forecasting for five regions of Australia is performed with high quality. Moreover, although the electrical load dataset of the QLD region is a difference stationary signal, the R value obtained for it is 0.9967 which is very close to one and supports the idea that an IT2AIFLS, because of its nonlinear nature, can handle this class of data as well.

Conclusions
In this paper, a IT2AIFLS is trained using a hybrid method containing the GSA for the antecedent part parameters and the R-LS algorithm for consequent part parameters. A benchmark system identification problem is studied to investigate the efficacy of the proposed system parameter tuning approach on previous identification benchmark problems . The comparisons with several other methods including Type-1 TSK FNS, Type-2  TSK FNS, Feedorward Type-2 FNN, SIT2FNN, SEIT2 FNN, TSCIT2FNN, IT2 FNN-GD,  IT2 FNN-SMC, IT2 FNNPSO+ SMC, IT2 IFLS -DEKF+GD, and IT2FLS with Modified SVR support the idea that the proposed approach is an efficient approach in system identification problems. The proposed approach is then investigated on electrical load prediction for five regions of Australia and Poland in the presence of noise and uncertainty which inherently exist in these datasets. Statistical properties of these datasets are presented that show they can be either stationary, trend stationary or difference stationary. In the case of the Poland electrical load dataset, the inputs to the IT2AIFLS are considered to be current values as well as three consecutive time delays of data. In the case of the five regions in Australia, current data values as well as four consecutive lags are used as the inputs to the fuzzy logic system. In both cases, one-step ahead prediction is considered. For the Poland dataset, the obtained prediction results are compared with several other algorithms including IT2FLS DEKF+GD, IFLS DEKF+GD and IT2IFLS DEKF+GD. The comparisons made in this paper show that the proposed algorithm results in superior performance when it is compared to these methods. For the case of the five regions of Australia, it is observed that the proposed prediction can perform with a much higher performance as compared with SVR, ANN, ELM, EEMD-ELM-GOA, EEMD-ELM-DA, EEMD-ELM-PSO and EEMD-ELM-GWO methods of prediction. Hence, the proposed approach is an automated solution to develop the antecedent part parameters of an IT2AIFLS, which successfully outperforms current state-of-the-art approaches in literature when applied to an electrical load prediction problem.
A possible direction of future work is the development of an interpretable IT2AIFLS which would lead to IF-THEN rules to describe the overall behavior of the electrical load system. Such IF-THEN rules are easy to communicate and may help experts to make more efficient decisions using their knowledge and experience. An additional direction of future work would be to apply the model to other time series prediction or system identification problems to test its applicability to enhance prediction in other industries.