Fractional-Order Model Predictive Frequency Control of an Islanded Microgrid

Optimal frequency control of an islanded microgrid has been a challenging issue in the research field of microgrids. Recently, fractional-order calculus theory and some related control methods have attempted to handle this issue. In this paper, a novel fractional-order model predictive control (FOMPC) method is proposed to achieve the optimal frequency control of an islanded microgrid by introducing a fractional-order integral cost function into model predictive control (MPC) algorithm. Firstly, a discrete state-space model is derived for the optimal frequency control problem of an islanded microgrid. Afterward, a fractional-order integral cost function is designed to guide the FOMPC algorithm to obtain optimal control law by borrowing the Grünwald-Letnikov (GL) definition of fractional order calculus. Six simulation studies have been carried out to illustrate the superiority of FOMPC to conventional MPC under dynamical load disturbances, perturbed system parameters and random dynamical power fluctuation of wind turbines.


Introduction
In the past decade, microgrids have been widely studied because of their potential environmental and economic benefits in regard to the integration distributed energy resources, especially renewable energy systems [1][2][3].Due to the different types of distributed energy resources, dynamical load characteristics, power quality constraints, and perturbed system parameters, how to optimize control of a microgrid, for it to operate well under grid-connected and islanded modes, has been a challenging issue [4][5][6].When a microgrid is islanded, an optimal controller should be designed to minimize the frequency deviations.In recent years, some research works have been reported related to the optimal frequency control issue of an islanded microgrid [7][8][9][10], but most works have focused on proportional-integral-derivative (PID) control method and its different kinds of improved versions based on intelligent optimization algorithms.
As one of the most widely-applied control methods in industry, model predictive control (MPC) ranks second only after PID control algorithm [11].Recently, MPC has been extended to solve the optimal load frequency control (LFC) problems of power systems [12][13][14][15][16] and operation optimization issues of microgrids [17][18][19][20].Especially, it should be noted that there are some existing works regarding MPC-based frequency control of microgrids [12,[19][20][21][22][23].Reference [12] developed a MPC-based coordinated load frequency control method for a microgrid under the condition regarding varied number of electric vehicles.By introducing a rule-based fuzzy controller to dynamically tune the weighting parameter used in cost function, a fuzzy MPC is proposed for the optimal LFC issue of an isolated microgrid [21].In Reference [22], a MPC method is developed to implement the torque compensation of each wind turbine generator for the optimal frequency control of an isolated grid.Recently-reported work [23] focuses on a conventional MPC-based optimal frequency control method in terms of a traditional integer-order cost function for a standalone microgrid.The simulation results have illustrated the superiority of MPC to traditional proportional-integral (PI) controller.
On the other hand, fractional-order control theories and applications of traditional integer-order and fractional-order systems have attracted increasing attention over the past two decades [24].The fractional-order PI/PID control algorithm has been considered as one of the most widely and successfully used methods in complex systems, such as servo systems [25,26], a nonlinear vertical tank system [27], electric balance vehicle system [28], an induction motor drive system [29], permanent magnetic synchronous generator [30], automatic voltage regulator systems [31], single-area delayed power systems [32], and multi-area interconnected power systems [33].One particular concern is the fractional-order frequency control of an isolated microgrid.Multi-objective extremal optimization algorithm has been adopted to optimize the fractional-order PID frequency controller for islanded microgrids [34].Additionally, some MPC or fractional-order MPC methods have been developed for traditional integer-order or fractional-order systems [35][36][37][38][39][40][41][42].However, to the best knowledge of the authors, few works regarding the fractional-order model predictive frequency control of islanded microgrids have been reported so far.
As one of the most important applications concerning fractional-order calculus theory in single processing and optimal control areas, a novel fractional-order integral performance index has been introduced as an improved cost function to improve the system performance.As a seminal research work, Romero et al. [43] proposed a new type of arbitrary real-order integral cost function for infinite impulse response filter design.The work in Reference [44] used a real-coded genetic algorithm to optimize the weighting matrices used in Linear Quadratic Regulator (LQR) systems by designing a fractional-order integral cost function.Furthermore, Romero et al. [36] also introduced a fractional-order definite integration operator into the cost function and then presented a fractional-order generalized predictive control method to deal with the optimal speed control issue for gasoline-propelled cars.Moreover, Bigdeli [37] designed a genetic algorithm-based fractional-order predictive functional controller (FOPFC) for fractional-order systems by adopting a fractional-order integral cost function.In Reference [45], Grünwald-Letnikov (GL) definition of fractional-order calculus theory has also been introduced in the cost function of a FOPFC to enhance the performance of fractional-order industrial processes.In addition, Zhang el al. [40] proposed a fractional-order MPC method for the optimal temperature control issue of industrial heating furnaces by using a fractional-order cost function based on the GL definition.These above research works have illustrated the superiority of fractional-order integral cost functions to conventional integer-order ones in the applications of signal processing, optimal LQR design and optimal predictive functional control.
Inspired by the aforementioned descriptions, this work proposes a fractional-order MPC, termed the FOMPC method, to achieve optimal frequency control performance for an islanded microgrid by introducing a fractional-order integral cost function into a model predictive control (MPC) algorithm.Firstly, an extended discrete state-space model is derived for the optimal frequency control problem of an islanded microgrid by considering some disturbance variables.Then, a fractional-order integral cost function is designed to guide the model predictive control algorithm to obtain optimal control law by borrowing the GL definition of fractional-order calculus theory.Six simulation studies on a typical microgrid have been carried out to illustrate the superiority of FOMPC over a conventional MPC under dynamical load disturbances, perturbed system parameters and the dynamic fluctuation of wind.
The major contribution of this paper can be summarized by the following three points: (1) A fractional-order model predictive control (FOMPC) method is proposed for the optimal frequency control problem of an islanded microgrid for the first time.
(2) The proposed FOMPC method extends the traditional quadratic cost function and introduces a fractional-order integral cost function based on the GL fractional-order integral definition.(3) In comparison with the conventional integer-order MPC algorithm, the proposed FOMPC can improve the performance of an islanded microgrid especially under several conditions of load disturbance, parameter uncertainty and random dynamical fluctuations of wind turbines.
The rest of this work is structured as follows.The discrete-time space model of an islanded microgrid is derived in Section 2. Section 3 proposes the FOMPC-based frequency control method.In Section 4, the simulation results for six cases and the effects of two fractional-order parameters on frequency control performance will be provided.Finally, the conclusions and some open issues are presented in Section 5.

Discrete-Time Space Model
The block diagram of the frequency control issue for an islanded microgrid with a FOMPC controller is shown in Figure 1.For convenient, fair comparison with the traditional integer-order MPC method [23], the system model of the microgrid tested in this paper is the same as that in Reference [23].The transfer functions and related system parameters of the distributed generations in an islanded microgrid under study are described in Table 1.The distributed subsystems in this microgrid under study include a wind turbine generator (WTG), diesel-engine generator (DEG), and energy storage system (ESS).As mentioned in Chapter 11 of Reference [46], EES can be considered as a fuel cell (FC)/battery, flywheel energy storage system (FESS), superconducting magnetic energy storage (MESS) or a combination of these.Here, we adopt FESS as the ESS in the microgrid under study.
Energies 2018 3 of 21 (1) A fractional-order model predictive control (FOMPC) method is proposed for the optimal frequency control problem of an islanded microgrid for the first time.
(2) The proposed FOMPC method extends the traditional quadratic cost function and introduces a fractional-order integral cost function based on the GL fractional-order integral definition.(3) In comparison with the conventional integer-order MPC algorithm, the proposed FOMPC can improve the performance of an islanded microgrid especially under several conditions of load disturbance, parameter uncertainty and random dynamical fluctuations of wind turbines.
The rest of this work is structured as follows.The discrete-time space model of an islanded microgrid is derived in Section 2. Section 3 proposes the FOMPC-based frequency control method.In Section 4, the simulation results for six cases and the effects of two fractional-order parameters on frequency control performance will be provided.Finally, the conclusions and some open issues are presented in Section 5.

Discrete-Time Space Model
The block diagram of the frequency control issue for an islanded microgrid with a FOMPC controller is shown in Figure 1.For convenient, fair comparison with the traditional integer-order MPC method [23], the system model of the microgrid tested in this paper is the same as that in Reference [23].The transfer functions and related system parameters of the distributed generations in an islanded microgrid under study are described in Table 1.The distributed subsystems in this microgrid under study include a wind turbine generator (WTG), diesel-engine generator (DEG), and energy storage system (ESS).As mentioned in Chapter 11 of Reference [46], EES can be considered as a fuel cell (FC)/battery, flywheel energy storage system (FESS), superconducting magnetic energy storage (MESS) or a combination of these.Here, we adopt FESS as the ESS in the microgrid under study., and ∆P m (t) represent the power deviations of ESS, WTG, governor, and turbine, respectively.∆u(t) represents the incremental form of control signal, ∆w(t) represents the input power fluctuation of WTG, ∆P L (t) is the load fluctuation, and ∆P e (t) represents the power deviation of the whole microgrid system, i.e., ∆P e (t) = ∆P w (t) − ∆P E (t) − ∆P m (t) − ∆P L (t).The continuous-time state space model for the frequency control of an islanded microgrid is defined by Equations ( 1) and (2).

Subsystem Transfer Function Parameters
.
where A c , B c1 , B c2 , and C are defined as follows: The discrete-time state space model of an islanded microgrid is obtained by the following equations: x where It should be noted, the discrete-time space model can also be derived by a similar method for a more complex microgrid system, with more types of distributed generations, e.g., hydric, thermal, or photovoltaic generations.The only difference is that x(t), u(t), w(t), A c , B c1 , B c2 , and C in Equations ( 1) and ( 2) should be redefined according to the topology structure and block diagram of a new microgrid system.

FOMPC-Based Optimal Frequency Control Method
In this section, we propose a novel FOMPC method for the optimal frequency control issue of an islanded microgrid.
Firstly, let us define the predictive output vector as , the predictive state vector as , the predictive disturbance , the predictive control vector as , where N x , N y , N u and N w are the dimensions of x(k), y(k), u(k) and w(k), respectively.Then, the predictive output vector Y P (k) is computed according to Equations ( 9)-( 13): . . .
Energies 2019, 12, 84 6 of 21 where P and M are the prediction horizon and the control horizon, respectively.Reference trajectory y r (k+p|k) used in this paper is defined as follows [11]: where η represents a soften factor, and c(k) represents the set value of system output.In terms of vector-based form, Equation ( 14) can be also reformulated using the following equation: From the perspective of traditional MPC, the optimal frequency control issue of an islanded microgrid is formulated as follows: where Q is a weight vector to adjust the performance of predicted errors, and R is another weight vector for the future control vector.u min , ∆u min , and y min are defined as the lower limits of u(k), ∆u(k), and y(k), respectively, and u max , ∆u max , and y max are the corresponding upper limits.
For an islanded microgrid, Figure 2 represents a schematic diagram of FOMPC-based optimal frequency control method.More specifically, by introducing the GL definition of fractional order calculus [43] into Equation ( 16), a fractional-order integral cost function J FO (k) is approximated by Equations ( 18)-( 21).The detailed derivation process from Formulation (16) to Formulation (18) can be found in Reference [43].
where η represents a soften factor, and c(k) represents the set value of system output.In terms of vector-based form, Equation ( 14) can be also reformulated using the following equation: From the perspective of traditional MPC, the optimal frequency control issue of an islanded microgrid is formulated as follows: where Q is a weight vector to adjust the performance of predicted errors, and R is another weight vector for the future control vector.umin, ∆umin, and ymin are defined as the lower limits of u(k), ∆u(k), and y(k), respectively, and umax, ∆umax, and ymax are the corresponding upper limits.
Based on the gradient descent method, we can obtain the following optimal control vector u(k): where I N u represents an identity matrix with N u rows and N u columns.
The flowchart of the proposed FOMPC-based optimal frequency control method for an islanded microgrid is shown in Figure 3.
where β and γ are the fractional-order parameters.
Based on the gradient descent method, we can obtain the following optimal control vector u(k): ( ) ( ) where u N I represents an identity matrix with Nu rows and Nu columns.
The flowchart of the proposed FOMPC-based optimal frequency control method for an islanded microgrid is shown in Figure 3.In conclusion, the detailed steps of the proposed FOMPC-based optimal frequency control method for an islanded microgrid are summarized as follows: Step 1: Input a discrete-time state space model with sampling period Ts for the frequency control issue of an islanded microgrid system described by Equations ( 7) and (8).In conclusion, the detailed steps of the proposed FOMPC-based optimal frequency control method for an islanded microgrid are summarized as follows: Step 1: Input a discrete-time state space model with sampling period T s for the frequency control issue of an islanded microgrid system described by Equations ( 7) and (8).
Energies 2019, 12, 84 8 of 21 Step 2: Initialize the algorithm parameters of FOMPC including maximum number of sampling times N max , P, M, Q, R, β and γ, and set k = 1.
Step 6: Update the optimal y(k) and x(k) under u(k).
Step 7: Set k = k + 1, and repeat the loop from Step 4 to Step 7 until k = N max .

Simulation Results
The simulation results in Reference [23] have illustrated the superiority of MPC to traditional PI controller for the optimal frequency control issue of an islanded microgrid.To make a fair comparison between the proposed FOMPC and traditional MPC method [23], simulation studies are carried out for the same microgrid system model, as shown in Figure 1 [23].The parameters of system model are set as in Table 1.Six experiments are designed as Table 2 to demonstrate the superiority of FOMPC to the traditional integer-order MPC method under dynamical load disturbances, perturbed system parameters and random dynamical power fluctuation of wind turbines.The parameter settings for FOMPC and MPC are shown in Table 3.All the following simulation experiments have been performed by the aid of MATLAB 2014b software on a 2.50 GHz PC with i7-6500U processor and 8GB RAM.

Experiment Condition
Case 1 Normal system parameters, ∆P L = 0.02 p.u. after 1 s.Case 2 Normal system parameters, dynamical fluctuations of ∆P L .Case 3 Normal system parameters, random dynamical fluctuations of ∆w.
Table 3.The parameters settings for FOMPC and MPC.

Case 1: Normal System Parameters and Step Load Fluctuations
In case 1, we consider the condition of normal system parameters of an islanded microgrid and ∆P L = 0.02 p.u. after 1 s. Figure 4 presents the comparative results of frequency deviation ∆f, and power deviation ∆P e obtained by FOMPC and MPC for an islanded microgrid under case 1. Obviously, the transient-state and steady-state responses of ∆f and ∆P e obtained by FOMPC are better than those by MPC [23].The performance of FOMPC and MPC methods is evaluated by the indices including the integral of absolute value of the error (IAE) for ∆f and ∆P e denoted as IAE 1 and IAE 2 , respectively, the overshoot of ∆f and ∆P e denoted as M p1 and M p2 , respectively, the rising time of ∆f and ∆P e denoted as t u1 and t u2 , respectively, the settling time of ∆f and ∆P e denoted as t s1 and t s2 , respectively, Energies 2019, 12, 84 9 of 21 and the steady-state error of ∆f and ∆P e denoted as E ss1 and E ss2 , respectively.The IAE 1 and IAE 2 are defined as follows: Table 4 shows the performance comparison of FOEMPC and MPC for an islanded microgrid under case 1.Clearly, all the indices except Ess 1 obtained by FOMPC are better than those by MPC [23].Table 4 shows the performance comparison of FOEMPC and MPC for an islanded microgrid under case 1.Clearly, all the indices except Ess1 obtained by FOMPC are better than those by MPC [23].

Case 2: Normal System Parameters and Dynamical Load Fluctuations
For case 2, the dynamical fluctuations of ∆PL (e.g., Figure 5a) and normal system parameters are considered.Figure 5b and Figure 5c compare ∆f and ∆Pe obtained by FOMPC and MPC for an islanded microgrid, respectively.From the perspective of quantitative comparison, the control performance of ∆f and ∆Pe obtained by FOMPC and MPC is evaluated by the indices including IAE1, IAE2, ITAE1, ITAE2, ISE1, ISE2, ITSE1, and ITSE2.To be more specific, the ITAE1, ITAE2, ISE1, ISE2, ITSE1, and ITSE2 are defined as follows:

Case 3: Normal System Parameters and Random Dynamical Fluctuations of Δw
In case 3, we designed an experiment under normal system parameters and random dynamical fluctuations of Δw shown in Figure 6a.It should be noted that the random dynamical fluctuations of Δw are modeled based on Reference [34]. Figure 6b,c compares ∆f and ∆Pe of FOMPC and MPC, respectively.Intuitively, the fluctuations of ∆f and ∆Pe obtained by FOMPC are also less than those by MPC [23].From a quantitative perspective, Table 6 compares the performance indices, such as IAE1, IAE2, ITAE1, ITAE2, ISE1, ISE2, ITSE1 and ITSE2 obtained by FOMPC and MPC.It is evident that FOMPC outperforms MPC under case 3.In case 3, we designed an experiment under normal system parameters and random dynamical fluctuations of ∆w shown in Figure 6a.It should be noted that the random dynamical fluctuations of ∆w are modeled based on Reference [34]. Figure 6b,c compares ∆f and ∆P e of FOMPC and MPC, respectively.Intuitively, the fluctuations of ∆f and ∆P e obtained by FOMPC are also less than those by MPC [23].From a quantitative perspective, Table 6 6a.It should be noted that the random dynamical fluctuations of Δw are modeled based on Reference [34]. Figure 6b,c compares ∆f and ∆Pe of FOMPC and MPC, respectively.Intuitively, the fluctuations of ∆f and ∆Pe obtained by FOMPC are also less than those by MPC [23].From a quantitative perspective, Table 6

Case 4: Perturbed System Parameters and Dynamical ∆PL
In order to test the robustness of FOMPC and MPC against perturbed system parameters and dynamical ∆PL, we design an experiment under the following condition: Tt = 0.6 (increase 50%), R = 1.5 (decrease 50%), H = 0.0833 (decrease 50%), and dynamical fluctuations of ∆PL (e.g., Figure 7a).Figure 7b and Figure 7c present the responses of ∆f and ∆Pe by FOMPC and MPC, respectively.It is apparent that the responses of ∆f and ∆Pe obtained by FOMPC are better than those by MPC because of its more smooth flucutations and faster stabilization.Furthermore, the corresponding comparison of IAE1, IAE2, ITAE1, ITAE2, ISE1, ISE2, ITSE1, and ITSE2 is shown in Table 7.It is obvious  In order to test the robustness of FOMPC and MPC against perturbed system parameters and dynamical ∆P L , we design an experiment under the following condition: T t = 0.6 (increase 50%), R = 1.5 (decrease 50%), H = 0.0833 (decrease 50%), and dynamical fluctuations of ∆P L (e.g., Figure 7a).Figure 7b,c present the responses of ∆f and ∆P e by FOMPC and MPC, respectively.It is apparent that the responses of ∆f and ∆P e obtained by FOMPC are better than those by MPC because of its more smooth flucutations and faster stabilization.Furthermore, the corresponding comparison of IAE 1 , IAE 2 , ITAE 1 , ITAE 2 , ISE 1 , ISE 2 , ITSE 1 , and ITSE 2 is shown in Table 7.It is obvious that the proposed FOMPC method is superior to the traditional MPC [23] under case 4.  For case 5, we consider the condition of perturbed system parameters, e.g., T t = 0.6 (increase 50%), R = 1.5 (decrease 50%), H = 0.0833 (decrease 50%), and random dynamical fluctuations of ∆w shown in Figure 8a. Figure 8b,c compare ∆f and ∆P e obtained by FOMPC and MPC under case 5, respectively.It is obvious that these fluctuations of ∆f and ∆P e obtained by FOMPC are less than MPC.Additionally, Table 8 also shows the comparison of the aforementioned eight indices.Clearly, the eight peformance indices of FOMPC are all better than those by MPC [23].In other words, FOMPC is more robust than MPC under case 5.For case 5, we consider the condition of perturbed system parameters, e.g., Tt = 0.6 (increase 50%), R = 1.5 (decrease 50%), H = 0.0833 (decrease 50%), and random dynamical fluctuations of ∆w shown in Figure 8a. Figure 8b and Figure 8c compare ∆f and ∆Pe obtained by FOMPC and MPC under case 5, respectively.It is obvious that these fluctuations of ∆f and ∆Pe obtained by FOMPC are less than MPC.Additionally, Table 8 also shows the comparison of the aforementioned eight indices.Clearly, the eight peformance indices of FOMPC are all better than those by MPC [23].In other words, FOMPC is more robust than MPC under case 5.

Case 6: The Worst Case
In case 6, we consider the worst case, i.e., under the following combined condition: Tt = 0.6 (increase 50%), R = 1.5 (decrease 50%), H = 0.0833 (decrease 50%), dynamical fluctuations of ∆PL shown in Figure 9a, and random dynamical fluctuations of Δw shown in Figure 9b.For case 6, Figure 9c and Figure 9d show the reponses of ∆f and ∆Pe by FOMPC and MPC, respectively.Obviously, It is apparent that the fluctuations of ∆f and ∆Pe by FOMPC are also more smooth than those by MPC even under the worst case.Moreover, Table 9 compares the corresponding performance indices obtained by FOMPC and MPC.Clearly, FOMPC performs better than MPC in terms of all the indices, which indicates that FOMPC has better robustness than the traditional MPC [23] under the worst case.In case 6, we consider the worst case, i.e., under the following combined condition: T t = 0.6 (increase 50%), R = 1.5 (decrease 50%), H = 0.0833 (decrease 50%), dynamical fluctuations of ∆P L shown in Figure 9a, and random dynamical fluctuations of ∆w shown in Figure 9b.For case 6, Figure 9c,d show the reponses of ∆f and ∆P e by FOMPC and MPC, respectively.Obviously, It is apparent that the fluctuations of ∆f and ∆P e by FOMPC are also more smooth than those by MPC even under the worst case.Moreover, Table 9 compares the corresponding performance indices obtained by FOMPC and MPC.Clearly, FOMPC performs better than MPC in terms of all the indices, which indicates that FOMPC has better robustness than the traditional MPC [23] under the worst case.

Case 6: The Worst Case
In case 6, we consider the worst case, i.e., under the following combined condition: Tt = 0.6 (increase 50%), R = 1.5 (decrease 50%), H = 0.0833 (decrease 50%), dynamical fluctuations of ∆PL shown in Figure 9a, and random dynamical fluctuations of Δw shown in Figure 9b.For case 6, Figure 9c and Figure 9d show the reponses of ∆f and ∆Pe by FOMPC and MPC, respectively.Obviously, It is apparent that the fluctuations of ∆f and ∆Pe by FOMPC are also more smooth than those by MPC even under the worst case.Moreover, Table 9 compares the corresponding performance indices obtained by FOMPC and MPC.Clearly, FOMPC performs better than MPC in terms of all the indices, which indicates that FOMPC has better robustness than the traditional MPC [23] under the worst case.As the aforementioned description of the FOMPC method in Section 3, fractional-order parameters β and γ play important roles in adjusting the frequency control performance.To illustrate the effects of β and γ on the frequency control performance, Figure 10 presents the frequency deviation ∆f under a set of different values for β and γ.Intuitively, the dynamical responses become faster, and the steady-state performance becomes better as the values of β and γ increase.

Conclusions and Open Issues
In this work, a novel FOMPC-based optimal frequency control method is designed for an islanded microgrid.The proposed FOMPC method has extended the traditional quadratic cost function and introduced a fractional-order integral cost function based on the GL fractional-order integral definition.Simulation results have demonstrated the superiority of FOMPC to conventional MPC [23] with integer-order cost function under dynamical load disturbances, perturbed system parameters and random dynamical power fluctuation of wind turbines.Moreover, this paper discussed the effects of two important fractional-order parameters used in the FOMPC method on frequency control performance.In conclusion, the FOMPC method can be considered as a promising solution for the optimal frequency control issue of an islanded microgrid.Nevertheless, some of the following open issues are still rather remarkable for future research.Firstly, the fractional-order calculus theory can also be introduced into the gradient descent method used in FOMPC.Secondly, it is still challenging to tune a set of optimal algorithm parameters used in FOMPC from the multi-objective optimization perspective.Furthermore, the extension of the proposed FOMPC to more complex power systems is another interesting research issue.Foundation of China (Grant Nos.LY16F030011 and LZ16E050002).

Conclusions and Open Issues
In this work, a novel FOMPC-based optimal frequency control method is designed for an islanded microgrid.The proposed FOMPC method has extended the traditional quadratic cost function and introduced a fractional-order integral cost function based on the GL fractional-order integral definition.Simulation results have demonstrated the superiority of FOMPC to conventional MPC [23] with integer-order cost function under dynamical load disturbances, perturbed system parameters and random dynamical power fluctuation of wind turbines.Moreover, this paper discussed the effects of two important fractional-order parameters used in the FOMPC method on frequency control performance.In conclusion, the FOMPC method can be considered as a promising solution for the optimal frequency control issue of an islanded microgrid.Nevertheless, some of the following open issues are still rather remarkable for future research.Firstly, the fractional-order calculus theory can also be introduced into the gradient descent method used in FOMPC.Secondly, it is still challenging to tune a set of optimal algorithm parameters used in FOMPC from the multi-objective optimization perspective.Furthermore, the extension of the proposed FOMPC to more complex power systems is another interesting research issue.

Figure 1 .
Figure 1.The block diagram of an islanded microgrid under study with a FOMPC-based frequency controller.

Figure 1 .
Figure 1.The block diagram of an islanded microgrid under study with a FOMPC-based frequency controller.

Figure 2 .
Figure 2. The schematic diagram of FOMPC-based optimal frequency control method.

Figure 2 .
Figure 2. The schematic diagram of FOMPC-based optimal frequency control method.

Figure 3 .
Figure 3.The flowchart of FOMPC-based optimal frequency control method for an islanded microgrid.

Figure 3 .
Figure 3.The flowchart of FOMPC-based optimal frequency control method for an islanded microgrid.

Figure 4 .
Figure 4. Responses of ∆f and ∆Pe obtained by FOMPC and MPC under case 1.

Figure 4 .
Figure 4. Responses of ∆f and ∆P e obtained by FOMPC and MPC under case 1.

Figure 5 .
Figure 5. Responses of ∆f and ∆Pe by FOEMPC and MPC under case 2.

Figure 5 .
Figure 5. Responses of ∆f and ∆P e by FOEMPC and MPC under case 2.

Figure 6 .
Figure 6.Responses of ∆f and ∆Pe by FOEMPC and MPC under case 3.

Figure 6 .
Figure 6.Responses of ∆f and ∆P e by FOEMPC and MPC under case 3.

Figure 7 .
Figure 7. Reponses of ∆f and ∆Pe by FOMPC and MPC under case 4.

Figure 8 .
Figure 8. Responses of ∆f and ∆P e by FOMPC and MPC under case 5.

Figure 8 .
Figure 8. Responses of ∆f and ∆Pe by FOMPC and MPC under case 5.

Figure 10 .
Figure 10.Effects of fractional-order parameters on the frequency control performance of an islanded microgrid.

Figure 10 .
Figure 10.Effects of fractional-order parameters on the frequency control performance of an islanded microgrid.

Table 1 .
The transfer functions and the related parameters used in each subsystem.

Table 1 .
The transfer functions and the related parameters used in each subsystem.∆P E (t), ∆P w (t), ∆P g (t)

Table 2 .
The conditions of six experiments for an islanded microgrid.

Table 4 .
Control performance indices obtained by FOMPC and MPC for an islanded microgrid under case 1.The best values are marked in bold.

Table 4 .
Control performance indices obtained by FOMPC and MPC for an islanded microgrid under case 1.The best values are marked in bold.For case 2, the dynamical fluctuations of ∆P L (e.g., Figure5a) and normal system parameters are considered.Figure5b,ccompare ∆f and ∆P e obtained by FOMPC and MPC for an islanded microgrid, respectively.From the perspective of quantitative comparison, the control performance of ∆f and ∆P e obtained by FOMPC and MPC is evaluated by the indices including IAE 1 , IAE 2 , ITAE 1 , ITAE 2 , ISE 1 , ISE 2 , ITSE 1 , and ITSE 2 .To be more specific, the ITAE 1 , ITAE 2 , ISE 1 , ISE 2 , ITSE 1 , and ITSE 2 are defined as follows: )

Table 5
also shows the control performance indices obtained by FOMPC and MPC for an islanded microgrid under case 2. Obviously, all the indices obtained by FOMPC are better than those byMPC 23.

Table 5 .
Performance comparison of FOMPC and MPC for an islanded microgrid under case 2. The best values are marked in bold.

Table 5
also shows the control performance indices obtained by FOMPC and MPC for an islanded microgrid under case 2. Obviously, all the indices obtained by FOMPC are better than those byMPC 23.

Table 5 .
Performance comparison of FOMPC and MPC for an islanded microgrid under case 2. The best values are marked in bold.

Table 6 .
Performance indices obtained by FOMPC and MPC for an islanded microgrid under case 3. The best values are marked in bold.

Table 6 .
Performance indices obtained by FOMPC and MPC for an islanded microgrid under case 3. The best values are marked in bold.

Table 7 .
Performance indices obtained by FOMPC and MPC for an islanded microgrid under case 4. The best values are marked in bold.

Table 8 .
Performance indices obtained by FOMPC and MPC for an islanded microgrid under case 5.The best values are marked in bold.

Table 8 .
Performance indices obtained by FOMPC and MPC for an islanded microgrid under case 5.The best values are marked in bold.

Table 8 .
Performance indices obtained by FOMPC and MPC for an islanded microgrid under case 5.The best values are marked in bold.

Table 9 .
Performance indices obtained by FOMPC and MPC for an islanded microgrid under case 6.The best values are marked in bold.Dynamical responses of ∆f and ∆P e obtained by FOMPC and MPC under case 6.

Table 9 .
Performance indices obtained by FOMPC and MPC for an islanded microgrid under case 6.The best values are marked in bold.As the aforementioned description of the FOMPC method in Section 3, fractional-order parameters β and γ play important roles in adjusting the frequency control performance.To illustrate the effects of β and γ on the frequency control performance, Figure10presents the frequency deviation ∆f under a set of different values for β and γ.Intuitively, the dynamical responses become faster, and the steady-state performance becomes better as the values of β and γ increase.