Next Article in Journal
Study on Abundant Dust-Ion-Acoustic Solitary Wave Solutions of a (3+1)-Dimensional Extended Zakharov–Kuznetsov Dynamical Model in a Magnetized Plasma and Its Linear Stability
Previous Article in Journal
Advances in Ostrowski-Mercer Like Inequalities within Fractal Space
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Marine-Predator-Algorithm-Based Optimum FOPID Controller for Enhancing the Stability and Transient Response of Automatic Voltage Regulators

by
Abdullah M. Noman
1,*,
Sulaiman Z. Almutairi
1,
Mokhtar Aly
2,
Mohammed H. Alqahtani
1,
Ali S. Aljumah
1 and
Emad A. Mohamed
1,3
1
Department of Electrical Engineering, College of Engineering, Prince Sattam bin Abdulaziz University, Al Kharj 16278, Saudi Arabia
2
Facultad de Ingeniería, Arquitectura y Diseño, Universidad San Sebastián, Bellavista 7, Santiago 8420524, Chile
3
Department of Electrical Engineering, Faculty of Engineering, Aswan University, Aswan 81542, Egypt
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(9), 690; https://doi.org/10.3390/fractalfract7090690
Submission received: 9 June 2023 / Revised: 29 August 2023 / Accepted: 30 August 2023 / Published: 16 September 2023

Abstract

:
An improved design optimization method for fractional-order-based proportional integral derivative (FOPID) controllers is proposed in this paper to enhance the stability and transient response of automatic voltage-regulator systems. The FOPID represents a higher degree-of-freedom controller through having five tunable parameters, compared with only three parameters in the integer-order PID controller. In the literature, the performance of the FOPID is highly determined through the design method and its parameter-determination process. Recently, optimum design of the FOPID has found wide employment in several engineering applications through using metaheuristic optimization algorithms. In this paper, an improved method for the FOPID’s parameter optimization is proposed for AVR applications using the marine predator optimization algorithm (MPA). The proposed MPA–FOPID controller is verified through comparing its performance with other featured and newly developed optimization algorithms. The proposed MPA–FOPID comparative analysis has been proven to have better stability, frequency response, robustness, faster response, and disturbance-rejection ability over the other studied methods in this paper.

1. Introduction

1.1. Overview

Recently, wide continuous variations exist in electrical power systems due to the changes in the type of energy source, grid modeling, and load behaviors. Recently, renewable-based generation started to denominate the new sources’ installation, followed by the corresponding grid proprieties changing. Based on that, achieving stable and constant voltage magnitudes and frequencies have become essential objectives for the design of control systems [1,2]. Fluctuating frequency and voltage magnitudes may lead to the degraded performance of their connected loads, especially their reliability and useful operating lifetime. The power system’s active power, reactive power, and associated power losses are directly linked with the variations in their frequency and voltage magnitudes. In cases of existing slight deviations in voltage magnitudes, sufficient changes exist in the reactive power’s magnitude. When voltage magnitude deviates out of the predefined ±5% limit around its rated voltage, the useful lifetime and operating efficiency of connected power system components/appliances are accordingly highly reduced [3].
Fortunately, voltage deviations can be properly controlled, and their effects consequently mitigated. The voltage can be regulated at various existing levels in power systems on the generation side, and/or on the transmission/distribution sides. On the generation side, automatic voltage regulators (AVRs) are employed for controlling voltage deviations [4], whereas on the transmission/distribution sides, various existing reactive power compensation devices can be utilized, such as flexible AC transmission systems (FACTS), transformers’ tab changers, filter devices, etc. The AVR systems have verified excellent regulation performance on the generation sides of synchronous generators. Optimized design of AVR parameters is essential for obtaining the best performance of these devices, considering the various uncertainties and parameter variations in AVRs, sensors, gains, etc. [5].

1.2. Literature Review

The integer order (IO) schemes have found wide use in AVR applications, especially the PID control scheme. It has proven better dynamic response and more stable performance, owing to its robustness and simple implementation. However, restrictions exist in the use of PID control schemes in modern control and industrial applications, due to their performance’s dependency on and sensitivity to the gains’ value-selection process [6]. Trial-and-error tuning has been presented in the literature; however, it consumes more time without ensuring an optimum control-parameter selection process. Diverse traditional tuning methods for PID controllers exist in the literature, such as the Ziegler–Nichols and Cohen–Coon methods. However, they are based on selecting a particular operating point for the tuning process when linearly evaluating the model. This, in turn, leads to the incapability to drive optimum gains and high overshoot peaks, and, accordingly, long oscillations exist in the system [7].
Additional improvements in AVR control have been proposed using fractional order (FO) control schemes. Using FO calculus with PID control in the FOPID provides enhanced performance and better degree of freedom in the control design [8]. The classical IO-based PID has three tunable control parameters (proportional (P), integral (I), and derivative (D) gains), whereas the FOPID has five tunable control parameters (P, I, and D gains in addition to the FO operators λ and μ for the I and D terms). Accordingly, FOPID-based AVR systems have better transients, more stability, and more robust performances than traditional IO-based PID AVR systems. Therefore, the FOPID has been utilized in this paper for AVR applications. The performance improvements of FOPID control can be achieved through optimizing its parameters to fulfill predefined criteria of system response (fitness/cost/objective function). Among the presented methods, metaheuristic-based optimization methods have been widely used, and are proven to have better tuning of parameters [9,10]. Another AVR control method based on PID plus the second-order derivative (PIDD2) AVR controller has been presented in [11]. A modified version of the Runge–Kutta optimization (RUN) algorithm has been proposed for optimizing the PIDD2 parameters.
Table 1 summarizes the applied optimization algorithms for tuning AVR control systems. The table shows that several algorithms have been applied in the literature for tuning PID and FOPID AVR controllers. The performance of these algorithms differs according to the operating principle of each algorithm. Recently, several new optimization algorithms have been proposed and showed improved performance in several applications. The marine predator algorithm (MPA) has been presented in [12]. It has proved better performance in the energy management of microgrids, several load-frequency control applications [13], parameter determinations, reconfiguration of PV and fuel cells [14,15], etc. Moreover, other optimizing algorithms have recently been presented for different applications. However, these algorithms have been applied in the AVR control design optimization. Therefore, their evaluation and comparisons are included in this paper.

1.3. Motivation

Based on the above-described literature review, obtaining stable and satisfactory AVR responses can be achieved with the proper designs and determinations of the best parameters. In addition, various existing optimization algorithms differ in their principal operating philosophy, which is reflected in their outputted best parameters for the optimization problem. Additionally, the FOPID-based AVR systems achieved better performance metrics than classical PID controllers. However, the evaluation of some recent optimization algorithms and the superior performance of the MPA optimizer has not been investigated. The MPA optimizer has achieved enhanced performance in several applications [12,13,14,15,47]. The main highlights of the operating principles and governing policies of MPA for achieving optimized foraging, interactions, and memories of marine predators are as follows [12]:
  • The Lévy-based strategy is used by marine predators for environments with low prey concentrations, and Brownian movement is used for areas with abundant prey. Both strategies share the same percentage of traversing various habitats within their lifetime;
  • The behavior of predators is changed by natural environment-based effects, such as eddy formations, or by human-based effects, such as fish aggregating devices (FAD), for seeking areas with prey distribution;
  • The velocity ratio (VR), which is represented by the ratio of predator velocity to prey velocity, is used for determining the best strategy. At low VR (VR = 0.1), Lévy is the best predators’ strategy, whether prey are moving using Brownian or Lévy strategies;
  • At unity VR (VR = 1), Brownian movement represents the best predator’s strategy when prey movement is through Lévy. The other scenarios depend on the system’s size.
  • At high VR (VR ≥ 10), the predator’s best strategy is not doing any movements at all, whereas prey’s movement is made through the Brownian or Lévy strategies;
  • The MPA benefits the good memory of marine predators at reminding their associates in addition to locations of successful foraging.

1.4. Paper Contributions

The paper’s contributions are summarized as follows:
  • The recent powerful MPA optimizer is presented and applied with the FOPID to improve the AVR controller. Based on authors’ knowledge at the submission date, this is the first time the MPA optimizer is presented in AVR controller design. The study is not limited to applying the MPA optimizer, it also presents a performance evaluation of MPA with several recently developed optimization algorithms are presented in this paper.
  • Additionally, the use of FOPID provides higher freedom with its additional parameters, which help improve the performance of AVR systems. The obtained optimum FOPID AVR controller is compared with the previously determined optimum FOPID using other optimizer techniques from the literature. The conducted design and analysis clearly demonstrate superior results of the proposed MPA-tuned FOPID AVR controller.
  • Better convergence performance and the determined parameters’ accuracy of the optimum FOPID AVR controller are presented in this paper using the MPA optimizer. The proposed MPA-based method is compared with recent and existing optimization methodologies. Additionally, several statistical tests are performed to make fair comparisons of optimization methods. The obtained results over 30 runs and the statistical analysis of the results confirm the superiority of MPA and its feasibility in AVR controller design.
The remaining of the paper is organized as follows: Section 2 presents state-of-the-art mathematical modeling of the AVR system, and the open loop analysis of AVR systems is presented in Section 3. A detailed description of the proposed optimization methodology is provided in Section 4. The AVR control and optimization process are detailed in Section 5. The obtained results are shown in Section 6, and paper conclusions are provided in Section 7.

2. Mathematical Description Modeling for AVR Systems

Figure 1 shows simplified general AVR modeling for this study. The control system is responsible for regulating the output voltage of the generator. It controls the excitation of the generator using a signal representing the error of subtracting the sensed output terminal voltage of the generator from the predefined reference voltage value. The operation of the AVR control can be explained at voltage drop or raise as follows:
  • At voltage drop: The output terminal voltage drops when there are increased loading conditions. In this scenario, the error between the sensed terminal voltage and the reference voltage increases with a positive value. Accordingly, the generator’s excitation increases until the sensed voltage reaches the predefined reference voltage value. This process continues until the sensed voltage equals the reference voltage. After reaching this condition, the generator’s excitation is maintained constant to preserve a stable supply voltage for all connected loads.
  • At voltage raise: At reduced load values, the output voltage increases, and hence an increase in the error signal happens, but with a negative value. Then, the generator’s excitation is reduced until it achieves equal sensed and reference terminal voltage values. Then, the generator’s excitation is maintained constant to stabilize the output terminal voltage.
The field windings of the generator have very high values for their inductance, and hence unavoidable conditions of load switching are generated in power systems. This, in turn, imposes several challenges on the AVR design process. The transfer function (TF) of the AVR system is modeled using the Laplace transform. The various AVR components (amplifier, exciter, generator, and sensor) are modeled using linearized components to facilitate TF modeling processes. The TF of AVR components ( G A ( s ) for amplifier, G E ( s ) for exciter, G G ( s ) for generator, and G S ( s ) for sensor) are represented with their reported parameters’ ranges in the literature as follows [41,48]:
G A ( s ) = K A 1 + s T A , and 10 K A 400 , 0.02 s T A 0.1 s G E ( s ) = K E 1 + s T E , and 1 K E 10 , 0.4 s T E 1 s G G ( s ) = K G 1 + s T G , and 0.7 K G 1 , 1 s T G 2 s G S ( s ) = K S 1 + s T S , and 1 K S 2 , 0.001 s T S 0.06 s
where K A , K E , K G , and K S represent gains, and T A , T E , T G , and T S represent time constants of the amplifier, exciter, generator, and sensor, respectively. The error between the measured voltage V m and its reference voltage V r e f is represented as E v in Figure 2. The AVR has to control the error E v to be minimum and to equal zero at steady state, whereas the AVR controller C ( s ) can be represented for PID C P I D ( s ) , and for FOPID C F O P I D ( s ) as follows:
G P I D ( s ) = K P + K I s + K D s G F O P I D ( s ) = K P + K I s λ + K D s μ
The complete first-order TF modeling of AVR components is shown in Figure 2. The complete system’s TF representing the output to input TF G s y s ( s ) can be expressed as follows:
G s y s ( s ) = C ( s ) G A ( s ) G E ( s ) G G ( s ) 1 + C ( s ) G A ( s ) G E ( s ) G G ( s ) G S ( s )

3. System Characteristics without Controller

The considered parameters’ values (1) are as follows: K A = 10 , K E = 1 , K G = 1 , K S = 1 , T A = 0.1 , T E = 0.4 , T G = 1 , and T S = 0.01 [42], whereas the generator’s gain K G depends on its loading level. At the no-load condition, K G equals 0.7, and it increases with loading to 1 at the nominal loading condition [42]. In this section, the analysis of the AVR system without the controller is carried out. Figure 3 shows different step responses of AVR systems at various loading conditions (varying K G at 0.7, 0.8, 0.9, and 1.0). The associated frequency response and the root locus plots of the effects of different loading conditions are shown in Figure 4 and Figure 5, respectively. It can be seen that the generator’s gain affects the system response in step time, the frequency response, and the root locus pole/zero locations.
The measured effects of the AVR system’s response are summarized in Table 2. The overshoot percentage and settling time increase at higher values of K G , whereas the rise time and steady state error decrease with the increase in K G . For the frequency response bode diagram, the gain margin and the phase margin decrease with increasing K G . In addition, the variations in K G values affect the closed-loop poles locations of the AVR system. Comparisons of the open-loop response values are included in the table.

4. MPA Optimization Methodology

The MPA represents a population-based optimizer, wherein uniform distribution over searching space exists in its initial solution. It is expressed as follows [12]:
X 0 = X m i n + r a n d ( X m a x X m i n )
where X m a x and X m i n are upper and lower boundaries of variables, and r a n d represents a uniform-random vector (between 0 and 1). The Elite, which is used for representing top predators, and prey are expressed as follows [47]:
E l i t e = X 1 , 1 I X 1 , 2 I X 1 , d I X 2 , 1 I X 2 , 2 I X 2 , d I X n , 1 I X n , 2 I X n , d I n × d
P r e y = X 1 , 1 X 1 , 2 X 1 , d X 2 , 1 X 2 , 2 X 2 , d X n , 1 X n , 2 X n , d n × d
where X I is top predators’ vector, n stands for searching agents’ number, d stands for dimensions’ number, and X i , j stands for the j t h dimension of i t h prey. The three operating scenarios of MPA can be expressed as follows:
At High VR Scenario: Within this scenario, the predators’ velocity is higher than prey’s velocity. It happens during initial optimization iterations. It is mathematically expressed as: while I t e r < 1 3 M a x I t e r [14],
S i = R B ( E l i t e i R B P r e y i ) , i = 1 , 2 , , n P r e y i = P r e y i + P . R S i
where R B is vector of random numbers for Brownian motion representation, ⊗ stands for entry-wise multiplication, P is a constant equal to 0.5, R stands for a random number vector between 0 and 1, and I t e r stands for the current iteration’s number within the maximum iterations number M a x I t e r , whereas multiplying R B with P r e y i simulates prey movements.
At Unity VR Scenario: This scenario happens while 1 3 M a x I t e r I t e r < 2 3 M a x I t e r , wherein predators use Brownian strategy and prey use Lévy strategy. Both of exploration and the exploitation happen in this scenario. Then, the population is divided between them, where prey are responsible for the exploitation part and predators are responsible for the exploration part. For the prey’s half of the population, it is expressed as [12]:
S i = R L ( E l i t e i R L P r e y i ) , i = 1 , 2 , , n / 2 P r e y i = P r e y i + P . R S i
where R L stands for Lévy distribution-based random number vector for Lévy movement representation. For the predators’ half of the population, it is expressed as [12]:
S i = R B ( R B E l i t e i P r e y i ) , i = n / 2 , , n P r e y i = E l i t e i + P . C F S i C F = ( 1 t M a x I t e r ) 2 I t e r M a x I t e r
At Low VR Scenario: This scenario happens while I t e r 2 3 M a x I t e r . It is mathematically expressed as [14]:
S i = R L ( R L E l i t e i P r e y i ) , i = 1 , 2 , , n P r e y i = E l i t e i + P . C F S i
In addition, FAD effects are mathematically expressed as [12]:
P r e y i = P r e y i + C F [ X m i n + R ( X m a x X m i n ) ] U If r F A D s P r e y i + [ F A D s ( 1 r ) + r ] ( P r e y r 1 P r e y r 2 ) If r > F A D s
where U stands for binary vector including arrays of 0 and 1, F A D s stands for FADs’ effects probability, which equals 0.2, r stands for a random number between 0 and 1, and P r e y r 1 and P r e y r 2 are random indexes within the prey matrix.

5. The AVR Control and Optimization Process

5.1. FO Theory and Representation

The FO control theory offers better flexibility and freedom compared with integer ones. The common representation of FO theory is provided by the Riemann–Liouville methodology, the Grunwald–Letnikov methodology, and the Caputo representation. For Grundwald–Letnikov, the α t h FO derivative part of function f and within a to t range limits is:
D α | a t = lim h 0 1 h α r = 0 t a h ( 1 ) r n r f ( t r h )
where h refers to step time, and the use of n is for fulfilling the ( n 1 < α < n ) condition. The coefficients in binomial are obtained as:
n r = Γ ( n + 1 ) Γ ( r + 1 ) Γ ( n r + 1 )
where Gamma refers to function and is defined as:
Γ ( n + 1 ) = 0 t x 1 e t d t
Additionally, Liouville and Riemann have provided definitions of derivatives that are capable of avoiding sums and limits. Wherein this definition uses the IO derivatives, the integrals are represented as:
D α | a t = 1 Γ ( n α ) ( d d t ) n a t f ( τ ) ( t τ ) α n + 1 d τ
The Caputo method representation defines FO derivatives as:
D α | a t = 1 Γ ( n α ) a t f ( n ) ( τ ) ( t τ ) α n + 1 d τ
whereas D α | a t is referred as:
D α | a t = α > 0 d α d t α FO derivative α < 0 t 0 t f d t α FO integral α = 0 1
Implementing of an FO-based control is commonly performed using an Oustaloup recursive approximation (ORA). It showed a suitable representation of FO derivatives for digital implementations in real-time applications. Also, it showed a better way to optimize FO control design using optimization techniques. Therefore, ORA FO representation is selected in this work to represent the FOPID AVR controller. The approximate mathematical ORA representation of the α t h derivative ( s α ) is expressed as:
s α ω h α k = N N s + ω k z s + ω k p
where ω k p and ω k z refer to pole and zero locations in ω h sequences, respectively, and they are calculated as:
ω k z = ω b ( ω h ω b ) k + N + 1 α 2 2 N + 1
ω k p = ω b ( ω h ω b ) k + N + 1 + α 2 2 N + 1
ω h α = ( ω h ω b ) α 2 k = N N ω k p ω k z
where the ORA approximated representations contains ( 2 N + 1 ) number of poles and zeros. Thence, N refers to ORA’s filter order within ( 2 N + 1 ) frequency ranges. The implemented ORA representation in this paper relies on using ( N = 5 ) in frequency range of ( ω [ ω b , ω h ] ), and accordingly selected between [ 0.001 , 1000 ] rad/s.

5.2. Proposed Optimized FOPID AVR Control

Figure 6 shows the schematic diagrams of PID- and FOPID-based AVR systems. It can be seen that the PID possesses three tunable gains, whereas the FOPID possesses five different tunable parameters. The FOPID has three gains and two FO operators, which increase its flexibility and degree of freedom compared to IO-based PID controller. The parameter optimization process is highly dependent on the objective function employed to drive the optimizer algorithm. The controller parameters are selected to minimize the predefined objective function. Usually, the measured output voltage and/or the error signals are used for constructing the objective functions.
There are four main error-based objective functions in the literature. Their general representation is as follows:
1.
Integral-squared-errors (ISE),
I S E = i = 1 m ( e i 2 ) d t
2.
Integral time-squared-errors (ITSE),
I T S E = i = 1 m ( e i 2 ) t . d t
3.
Integral-absolute-errors (IAE),
I A E = i = 1 m a b s ( e i ) d t
4.
Integral time-absolute-errors (ITAE),
I T A E = i = 1 m a b s ( e i ) t . d t
Moreover, the various error and/or performance indices can be combined in the summation and employed for constructing functions with multiple objectives. Weighting factors for each error or objective can also be added to the general form to manage the priority of each element. Table 3 summarizes the various objective function representations in the literature.
The ITAE is used as an objective function in this paper for the optimization problem due to its generality. Figure 7 shows the schematic diagram of the proposed AVR controller optimization process. The measured error voltage signal is used to calculate the objective function for each simulation run, and the MPA optimizer works to determine the optimum parameters as clarified in Algorithm 1. It can be seen that there are three tunable parameters with the PID controller and five tunable ones with the FOPID controller. The upper and lower boundaries for the proposed optimization process are as follows:
K P m i n K P K P m a x K I m i n K I K I m a x K D m i n K D K D m a x λ m i n λ λ m a x μ m i n μ μ m a x
where ( f ) m i n and ( f ) m a x are lower and upper bounds of the tunable AVR controller’s parameter, respectively, whereas K P m i n , K I m i n , K D m i n are set at 1, 0.1, and 0.1, respectively, and K P m a x , K I m a x , and K D m a x , are set at 2, 1, and 0.4, respectively. In addition, μ m i n , and λ m i n are set at 1, whereas μ m a x and λ m a x are set at 2 in the proposed optimization process. For a fair comparison, all the optimization algorithms are set at 50 iterations, and the population size is set at 8.
Algorithm 1 Pseudo-code representing proposed optimization process using MPA
Define MPA parameters (Size of population, and maximum iteration M a x I t e r
Start algorithm initialization using (4)
Construct initial E l i t e using (5)
Initialize the searching agents (Preys) population using (6)
while   I t e r < M a x I t e r   do
      Calculate objective function, and construct E l i t e matrix
      if I t e r < M a x I t e r / 3 then
         Update preys using (7)
      else if M a x I t e r / 3 I t e r < 2 M a x I t e r / 3  then
         For preys’ half of population ( i = 1 , , n / 2 ) , update preys using (8)
         For predators’ half of population ( i = n / 2 , , n ) , update preys using (9)
      else if I t e r 2 M a x I t e r / 3 then
         Update preys using (10)
      end if
      Accomplish the memory savings, and the E l i t e updates
      Applying the FADs effects using (11)
      Calculate the objective function
end while
Return the best solution

6. Results and Discussion

Performance Comparisons with Recent Algorithms

Firstly, comparisons of the performance of MPA optimizers was performed with optimizers featured in the literature recently. The considered optimizers include the Artificial Hummingbirds Algorithm (AHA) [49], arithmetics trigonometric optimizer algorithm (ATOA) [50], Differential Evolutions (DE) [51], Heap-Based Optimization (HBO) [52], Sine–Cosine Optimizer Algorithm (SCA) [53], Slime Mould Algorithm (SMA) [54], and Manta Ray-Foraging Optimization (MRFO) [55]. The proposed AVR controller and optimization algorithms are simulated using MATLAB R2021a software. The AVR system and controller are implemented in Simulink and linked with the m-file program, which contains the optimization algorithms. The optimization process is made 30 times using the same number of iterations, population size, and parameters limits. The obtained results over the 30 runs are summarized in Table 4.
The obtained best FOPID parameters for each optimizer over the 30 runs are summarized in Table 5. It also provides statistical analysis of the best value, worst value, mean value, median value, and standard deviation for each optimizer during the 30 runs. The proposed MPA gives the lowest value of the objective function compared to the studied optimizers of 0.0151, whereas the SMA optimizer provides the highest value of the objective function of 0.0233. This proves the best performance of MPA over the all runs, whereas the statistical analysis shows that the SCA provides the highest value of the worst objective of 0.0763 compared to the lowest value of worst objective of 0.0176 using the proposed MPA optimizer. Additionally, the proposed MPA provides the best results in terms of mean, median, and standard deviation over the 30 runs compared to the studied optimizers. This proves the ability and stability of proposed method over the performed runs and random processes of different optimizers.
Another test based on 100 runs of the MRFO and MPA is performed, and the obtained results are shown in Table A1 in Appendix A. The associated statistical results are summarized in Table 6. It can be seen from the table that the best objective obtained by MPA is 0.01507709, compared to 0.01534512 by the MRFO algorithm. Additionally, the mean value achieved by the MPA is 0.01564708, compared to 0.01704896 by MRFO. This, in turn, confirms the improved design performance using MPA compared to the MRFO algorithm.
An analysis of variance (ANOVA) test has been performed on the obtained results of 30 runs. The results are shown in Table 7, and their corresponding rankings are shown in Figure 8. In the ANOVA test, when the obtained value for F is higher than the obtained p-value, the null hypothesis will be true. The obtained results verify that the obtained F-value is much higher than the obtained p-value, showing a significant difference among the obtained results. From Figure 8, the proposed optimization using MPA outperforms other recently developed methods that are considered in this analysis. The proposed MPA-based design provides the smallest variance range with the lowest objective value (for our minimization problem for optimization). This confirms the resilience and accuracy of proposed method. Additionally, the Tukey Honestly Significant Difference (the Tukey HSD) post hoc analysis has been provided to provide support for the obtained ANOVA test-based results. Figure 9 provides the obtained HSD post hoc results. The proposed MPA-based optimization provides the lowest value for the mean fitness function. It is followed by the MRFO and SMA optimizers. The obtained results from the two tests confirm the superior algorithm robustness of the MPA based design for determining optimal FOPID AVR controller parameters.
The obtained optimal parameters of various considered algorithms are simulated in a closed-loop control with the AVR model. The control parameters associated with the best objective function from the 30 runs are used in the simulation study for fair comparisons. The obtained results with step change in the reference voltage at zero second and full-load condition with K G = 1 are shown in Figure 10. It can be seen that the worst response is obtained by the ATOA, SMA, and SCA optimizers, whereas the MPA and MRFO provide the best results. Although MPA appears to have a slightly higher peak, it provides faster, reaching a steady-state operating point. The measured ISE, IAE, ITSE, and ITAE for MPA are 2.7334, 2.857, 4.1826, and 4.3812, respectively, whereas the MPA provides 2.8728, 2.9025, 4.4870, and 4.6129, respectively.
In this part, a comparison of the step response of the proposed design is performed with some featured designs from the literature. The obtained step response comparisons are shown in Figure 11 at full-load condition with K G = 1. The considered FOPID designs and parameters include method 1 as in [43], method 2 as in [44], method 3 as in [37], method 4 as in [45], method 5 as in [46], and the proposed MPA design. The proposed MPA-based FOPID AVR regulator provides reduced peak values and fast tracking in this step change scenario. Although method 2 provides a lower peak overshoot value, it has a steady state error as shown in the zoomed-in results. The proposed MPA-based design achieves zero steady state error in this scenario. The method 5 design has the highest peak overshoot value, followed by the methods 3 and 4. The design of method 1 has a high rise time and oscillations. Another case study at the no-load condition with K G = 0.7 is performed and shown in Figure 12. The results confirm the superior response of the proposed method compared with studied methods from the literature at no-load condition. The proposed MPA-based design has a low peak overshoot value, a low rise time, and zero steady state error. This proves the superior performance of the proposed design method compared to existing methods in the literature.
A comparison of the convergence characteristics is provided in Figure 13. The selected algorithms for the comparison include HBO, SMA, and ATOA methods. It can be seen that the MPA optimizer has good conversion speed compared with the other algorithms. While the most important criteria is the ability of MPA to have a reduced objective function at the end of the optimization process, it can be seen that MPA has the minimum objective function.
Another check is made by adding a disturbance of 0.3 p.u. to the output voltage measurement. The obtained results are shown in Figure 14. The disturbance is injected at 6 s and removed at 10 s. It can be seen that the MPA-based FOPID is robust against disturbances, and it can recover its state very fast with lower fluctuations, whereas the other studied designs have higher transients and/or higher steady state error and response times. Another comparison is made with the PIDD2 tuned by the RUN algorithm, and the parameters are taken from [11]. The obtained results are shown in Figure 15. It can be seen that the proposed MPA-based FOPID has better performance and response characteristics compared with the existing PIDD2 method.

7. Conclusions

This paper deals with an important issue of modern power systems that is related to the AVR controller and its design optimization. The paper proposes an optimized FOPID AVR controller that is optimized with the powerful marine predator optimization algorithm (MPA). The proposed AVR controller achieves an improved transient response, while preserving system stability. Additionally, the FOPID provides additional design freedom due to the extra included fractional order operators. Compared with the integer-order PID, the FOPID has two additional design flexibilities, which are reflected it its response and optimization process. Comprehensive statistical analysis of the MPA and obtained results over 30 runs are provided and compared with featured methods from the literature. The proposed FOPID-MPA method provides the best objective function minimization with reduced variations with runs compared to the other studied methods. The ANOVA and Tukey tests are provided in addition to statistical parameters, which proved to be the best performance obtained by the proposed FOPID-MPA method. Future research includes consideration of various constraints of AVR elements and further comparisons with other control methods. Additionally, more considerations can be given to other measurements, such as noise and stator current limitations.

Author Contributions

A.M.N., M.A., M.H.A., S.Z.A. and A.S.A.; Methodology, A.M.N., M.A., S.Z.A. and E.A.M.; Software, M.A., M.H.A., A.S.A., S.Z.A. and E.A.M.; Validation, A.M.N.; Formal analysis, A.M.N., M.H.A., and S.Z.A.; Investigation, A.M.N., E.A.M., M.A., A.S.A. and S.Z.A.; Resources, A.M.N., M.H.A., A.S.A. and S.Z.A.; Data curation, M.A. and E.A.M.; Writing—original draft, M.H.A., M.A. and E.A.M.; Writing—review and editing, A.M.N., M.A., A.S.A., S.Z.A., M.H.A. and E.A.M.; Visualization, A.S.A. and M.H.A.; Supervision, A.S.A., S.Z.A., A.M.N. and E.A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded By Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia, Project Number IF2/PSAU/2022/01/23150.

Data Availability Statement

Not applicable.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number (IF2/PSAU/2022/01/23150).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Obtained objective functions for the MRFO and MPA over 100 runs.
Table A1. Obtained objective functions for the MRFO and MPA over 100 runs.
Run No.MPAMRFORun No.MPAMRFO
10.01609023091533650.0186072735390579510.01597639536272260.0164159935337242
20.01577542143860250.0157790581530151520.01546253642353440.0161561738059689
30.01638397553005650.0197874911366351530.01546739808089970.0216327342483548
40.01586143602933160.0207406428821346540.01511716657792890.0171331501582164
50.01518151415067850.0178342633107935550.01516158133448830.0177126323347514
60.01668394738080190.0199929483957633560.01562678873274160.0164519137830847
70.01708011252469330.0155652253629703570.01658679568837010.0154761204614243
80.01708618578960190.0163012176667067580.01553629433242460.0156001843471987
90.01562173779819340.0174133359370687590.01529525312652950.0191971202041615
100.01615826317461100.0164046438255427600.01540085212269800.0169098714187219
110.01704679758598540.0180406862499614610.01516324567411650.0161135094007944
120.01827223518005690.0185911187833772620.01535273561341890.0154955036509552
130.01530362844629150.0153519464778942630.01512616699517960.0160335763536030
140.01580507839639080.0153868044792599640.01520636174955590.0207565298970893
150.01527477003470150.0158201433770507650.01520062271439230.0155734274747462
160.01607245062678450.0187138746721933660.01521263994191120.0169535389239530
170.01550619488590950.0186836337099605670.01550419420547950.0153554028201223
180.01654194552597260.0185533825453101680.01550974836226220.0160391558592351
190.01532393012205920.0223799180060417690.01552296648336010.0161831344490288
200.01572485215839070.0158234562111405700.01522693942234030.0171716717715401
210.01525442994952550.0156292615025043710.01572370130823330.0182032626875426
220.01562515264036130.0160019396955112720.01577142346943870.0185312561464443
230.01591623081421320.0178834956451730730.01572369342387720.0160178288055126
240.01616224082814470.0166977069643665740.01508551547987880.0156785500585492
250.01507709127999630.0160946478966843750.01567490313366720.0183338207041531
260.01592600026292240.0202221745201756760.01526787886138150.0177098353253446
270.01550294362707590.0178024023309439770.01546485276682560.0213044716721203
280.01582564737035950.0159423598815060780.01521155317284420.0175407533787249
290.01529818244512360.0154539311942310790.01519195568093190.0154529560023485
300.01635526801012540.0168646119783646800.01848792566538290.0160860584564143
310.01530777965950180.0166950620491384810.01509635504837340.0161007964844298
320.01555075855766590.0163582696073985820.01589979723854650.0169265091642214
330.01551086493077970.0155126486258078830.01531721635775430.0158997972385465
340.01600348031045510.0183471186539829840.01540469350667510.0174072210013845
350.01536756947821600.0231118909681678850.01513925513670360.0194813295050308
360.01570116132807080.0162058702232047860.01550985812053770.0159186921213136
370.01546416008241970.0159433449512397870.01541223828258140.0166866119939868
380.01522874699115070.0158987438068760880.01551008296136500.0161010349179289
390.01538538974708700.0157847475047855890.01568970182781600.0166741742074250
400.01577033588834890.0173888037944703900.01531283348319470.0159414326243861
410.01543135963869180.0170021705087658910.01535752652658520.0168227013385350
420.01605201725011770.0157239653675728920.01545028952787780.0164106715341243
430.01538011326845940.0185347128325756930.01511503402183860.0160006097916894
440.01646149238402200.0153723037170996940.01578870277255650.0155413006166845
450.01537024217252250.0153723037170996950.01512125507398550.0155006512376765
460.01612668919646630.0153451195810031960.01555006891420180.0161683560384650
470.01538670522999260.0172731580353522970.01517890382046400.0178129214344732
480.01536352912238030.0154699708648064980.01527345205891170.0154727513622186
490.01532844436913740.0163318577343206990.01509109425306820.0174838967070398
500.01529021300499220.01563099248347811000.01563449700376470.0179487980457389

References

  1. Hassan, A.; Aly, M.; Elmelegi, A.; Nasrat, L.; Watanabe, M.; Mohamed, E.A. Optimal Frequency Control of Multi-Area Hybrid Power System Using New Cascaded TID-PIλDμN Controller Incorporating Electric Vehicles. Fractal Fract. 2022, 6, 548. [Google Scholar] [CrossRef]
  2. Amin, A.; Ebeed, M.; Nasrat, L.; Aly, M.; Ahmed, E.M.; Mohamed, E.A.; Alnuman, H.H.; Hamed, A.M.A.E. Techno-Economic Evaluation of Optimal Integration of PV Based DG with DSTATCOM Functionality with Solar Irradiance and Loading Variations. Mathematics 2022, 10, 2543. [Google Scholar] [CrossRef]
  3. Said, S.M.; Aly, M.; Hartmann, B.; Mohamed, E.A. Coordinated fuzzy logic-based virtual inertia controller and frequency relay scheme for reliable operation of low-inertia power system. IET Renew. Power Gener. 2021, 15, 1286–1300. [Google Scholar] [CrossRef]
  4. Alghamdi, S.; Sindi, H.F.; Rawa, M.; Alhussainy, A.A.; Calasan, M.; Micev, M.; Ali, Z.M.; Aleem, S.H.E.A. Optimal PID Controllers for AVR Systems Using Hybrid Simulated Annealing and Gorilla Troops Optimization. Fractal Fract. 2022, 6, 682. [Google Scholar] [CrossRef]
  5. Ghasemi, M.; Rahimnejad, A.; Gil, M.; Akbari, E.; Gadsden, S.A. A self-competitive mutation strategy for Differential Evolution algorithms with applications to Proportional–Integral–Derivative controllers and Automatic Voltage Regulator systems. Decis. Anal. J. 2023, 7, 100205. [Google Scholar] [CrossRef]
  6. Alilou, M.; Azami, H.; Oshnoei, A.; Mohammadi-Ivatloo, B.; Teodorescu, R. Fractional-Order Control Techniques for Renewable Energy and Energy-Storage-Integrated Power Systems: A Review. Fractal Fract. 2023, 7, 391. [Google Scholar] [CrossRef]
  7. Daraz, A.; Malik, S.A.; Basit, A.; Aslam, S.; Zhang, G. Modified FOPID Controller for Frequency Regulation of a Hybrid Interconnected System of Conventional and Renewable Energy Sources. Fractal Fract. 2023, 7, 89. [Google Scholar] [CrossRef]
  8. Almasoudi, F.M.; Magdy, G.; Bakeer, A.; Alatawi, K.S.S.; Rihan, M. A New Load Frequency Control Technique for Hybrid Maritime Microgrids: Sophisticated Structure of Fractional-Order PIDA Controller. Fractal Fract. 2023, 7, 435. [Google Scholar] [CrossRef]
  9. Ahmed, E.M.; Mohamed, E.A.; Elmelegi, A.; Aly, M.; Elbaksawi, O. Optimum Modified Fractional Order Controller for Future Electric Vehicles and Renewable Energy-Based Interconnected Power Systems. IEEE Access 2021, 9, 29993–30010. [Google Scholar] [CrossRef]
  10. Ahmed, E.M.; Selim, A.; Alnuman, H.; Alhosaini, W.; Aly, M.; Mohamed, E.A. Modified Frequency Regulator Based on TIλ-TDμFF Controller for Interconnected Microgrids with Incorporating Hybrid Renewable Energy Sources. Mathematics 2022, 11, 28. [Google Scholar] [CrossRef]
  11. Izci, D.; Ekinci, S.; Mirjalili, S. Optimal PID plus second-order derivative controller design for AVR system using a modified Runge–Kutta optimizer and Bode’s ideal reference model. Int. J. Dyn. Control. 2022, 11, 1247–1264. [Google Scholar] [CrossRef]
  12. Faramarzi, A.; Heidarinejad, M.; Mirjalili, S.; Gandomi, A.H. Marine Predators Algorithm: A nature-inspired metaheuristic. Expert Syst. Appl. 2020, 152, 113377. [Google Scholar] [CrossRef]
  13. Sobhy, M.A.; Abdelaziz, A.Y.; Hasanien, H.M.; Ezzat, M. Marine predators algorithm for load frequency control of modern interconnected power systems including renewable energy sources and energy storage units. Ain Shams Eng. J. 2021, 12, 3843–3857. [Google Scholar] [CrossRef]
  14. Aly, M.; Ahmed, E.M.; Rezk, H.; Mohamed, E.A. Marine Predators Algorithm Optimized Reduced Sensor Fuzzy-Logic Based Maximum Power Point Tracking of Fuel Cell-Battery Standalone Applications. IEEE Access 2021, 9, 27987–28000. [Google Scholar] [CrossRef]
  15. Yousri, D.; Babu, T.S.; Beshr, E.; Eteiba, M.B.; Allam, D. A Robust Strategy Based on Marine Predators Algorithm for Large Scale Photovoltaic Array Reconfiguration to Mitigate the Partial Shading Effect on the Performance of PV System. IEEE Access 2020, 8, 112407–112426. [Google Scholar] [CrossRef]
  16. Gaing, Z.L. A particle swarm optimization approach for optimum design of PID controller in AVR system. IEEE Trans. Energy Convers. 2004, 19, 384–391. [Google Scholar] [CrossRef]
  17. Gozde, H.; Taplamacioglu, M. Comparative performance analysis of artificial bee colony algorithm for automatic voltage regulator (AVR) system. J. Frankl. Inst. 2011, 348, 1927–1946. [Google Scholar] [CrossRef]
  18. GÜVENÇ, U.; YİĞİT, T.; IŞIK, A.H.; AKKAYA, İ. Performance analysis of biogeography-based optimization for automatic voltage regulator system. Turk. J. Electr. Eng. Comput. Sci. 2016, 24, 1150–1162. [Google Scholar] [CrossRef]
  19. Köse, E. Optimal Control of AVR System With Tree Seed Algorithm-Based PID Controller. IEEE Access 2020, 8, 89457–89467. [Google Scholar] [CrossRef]
  20. Hekimoglu, B.; Ekinci, S. Grasshopper optimization algorithm for automatic voltage regulator system. In Proceedings of the 2018 5th International Conference on Electrical and Electronic Engineering (ICEEE), IEEE, Istanbul, Turkey, 3–5 May 2018. [Google Scholar] [CrossRef]
  21. Sahu, B.K.; Panda, S.; Mohanty, P.K.; Mishra, N. Robust analysis and design of PID controlled AVR system using Pattern Search algorithm. In Proceedings of the 2012 IEEE International Conference on Power Electronics, Drives and Energy Systems (PEDES), Bengaluru, India, 16–19 December 2012; pp. 1–6. [Google Scholar] [CrossRef]
  22. Mosaad, A.M.; Attia, M.A.; Abdelaziz, A.Y. Whale optimization algorithm to tune PID and PIDA controllers on AVR system. Ain Shams Eng. J. 2019, 10, 755–767. [Google Scholar] [CrossRef]
  23. Habib, S.; Abbas, G.; Jumani, T.A.; Bhutto, A.A.; Mirsaeidi, S.; Ahmed, E.M. Improved Whale Optimization Algorithm for Transient Response, Robustness, and Stability Enhancement of an Automatic Voltage Regulator System. Energies 2022, 15, 5037. [Google Scholar] [CrossRef]
  24. Dogruer, T.; Can, M.S. Design and robustness analysis of fuzzy PID controller for automatic voltage regulator system using genetic algorithm. Trans. Inst. Meas. Control. 2022, 44, 1862–1873. [Google Scholar] [CrossRef]
  25. Bingul, Z.; Karahan, O. A novel performance criterion approach to optimum design of PID controller using cuckoo search algorithm for AVR system. J. Frankl. Inst. 2018, 355, 5534–5559. [Google Scholar] [CrossRef]
  26. Hekimoğlu, B. Sine-cosine algorithm-based optimization for automatic voltage regulator system. Trans. Inst. Meas. Control. 2018, 41, 1761–1771. [Google Scholar] [CrossRef]
  27. Ekinci, S.; Hekimoğlu, B. Improved Kidney-Inspired Algorithm Approach for Tuning of PID Controller in AVR System. IEEE Access 2019, 7, 39935–39947. [Google Scholar] [CrossRef]
  28. Bendjeghaba, O. Continuous Firefly Algorithm for Optimal Tuning of Pid Controller in AVR System. J. Electr. Eng. 2014, 65, 44–49. [Google Scholar] [CrossRef]
  29. Çelik, E.; Durgut, R. Performance enhancement of automatic voltage regulator by modified cost function and symbiotic organisms search algorithm. Eng. Sci. Technol. Int. J. 2018, 21, 1104–1111. [Google Scholar] [CrossRef]
  30. Ekinci, S.; Hekimoğlu, B.; Kaya, S. Tuning of PID Controller for AVR System Using Salp Swarm Algorithm. In Proceedings of the 2018 International Conference on Artificial Intelligence and Data Processing (IDAP), Malatya, Turkey, 28–30 September 2018; pp. 1–6. [Google Scholar] [CrossRef]
  31. Anbarasi, S.; Muralidharan, S. Enhancing the Transient Performances and Stability of AVR System with BFOA Tuned PID Controller. J. Control. Eng. Appl. Inform. 2016, 18, 20–29. [Google Scholar]
  32. DUMAN, S.; YÖRÜKEREN, N.; ALTAŞ, İ.H. Gravitational search algorithm for determining controller parameters in an automatic voltage regulator system. Turk. J. Electr. Eng. Comput. Sci. 2016, 24, 2387–2400. [Google Scholar] [CrossRef]
  33. Pradhan, R.; Majhi, S.K.; Pati, B.B. Design of PID controller for automatic voltage regulator system using Ant Lion Optimizer. World J. Eng. 2018, 15, 373–387. [Google Scholar] [CrossRef]
  34. Mohanty, P.K.; Sahu, B.K.; Panda, S. Tuning and Assessment of Proportional–Integral–Derivative Controller for an Automatic Voltage Regulator System Employing Local Unimodal Sampling Algorithm. Electr. Power Components Syst. 2014, 42, 959–969. [Google Scholar] [CrossRef]
  35. Blondin, M.; Sicard, P.; Pardalos, P. Controller Tuning Approach with robustness, stability and dynamic criteria for the original AVR System. Math. Comput. Simul. 2019, 163, 168–182. [Google Scholar] [CrossRef]
  36. Zamani, M.; Karimi-Ghartemani, M.; Sadati, N.; Parniani, M. Design of a fractional order PID controller for an AVR using particle swarm optimization. Control. Eng. Pract. 2009, 17, 1380–1387. [Google Scholar] [CrossRef]
  37. Ortiz-Quisbert, M.E.; Duarte-Mermoud, M.A.; Milla, F.; Castro-Linares, R.; Lefranc, G. Optimal fractional order adaptive controllers for AVR applications. Electr. Eng. 2016, 100, 267–283. [Google Scholar] [CrossRef]
  38. Zhang, D.L.; Tang, Y.G.; Guan, X.P. Optimum Design of Fractional Order PID Controller for an AVR System Using an Improved Artificial Bee Colony Algorithm. Acta Autom. Sin. 2014, 40, 973–979. [Google Scholar] [CrossRef]
  39. Tang, Y.; Cui, M.; Hua, C.; Li, L.; Yang, Y. Optimum design of fractional order PIλDμ controller for AVR system using chaotic ant swarm. Expert Syst. Appl. 2012, 39, 6887–6896. [Google Scholar] [CrossRef]
  40. Zeng, G.Q.; Chen, J.; Dai, Y.X.; Li, L.M.; Zheng, C.W.; Chen, M.R. Design of fractional order PID controller for automatic regulator voltage system based on multi-objective extremal optimization. Neurocomputing 2015, 160, 173–184. [Google Scholar] [CrossRef]
  41. Ayas, M.S.; Sahin, E. FOPID controller with fractional filter for an automatic voltage regulator. Comput. Electr. Eng. 2021, 90, 106895. [Google Scholar] [CrossRef]
  42. Micev, M.; Ćalasan, M.; Oliva, D. Fractional Order PID Controller Design for an AVR System Using Chaotic Yellow Saddle Goatfish Algorithm. Mathematics 2020, 8, 1182. [Google Scholar] [CrossRef]
  43. Pan, I.; Das, S. Frequency domain design of fractional order PID controller for AVR system using chaotic multi-objective optimization. Int. J. Electr. Power Energy Syst. 2013, 51, 106–118. [Google Scholar] [CrossRef]
  44. Sikander, A.; Thakur, P.; Bansal, R.; Rajasekar, S. A novel technique to design cuckoo search based FOPID controller for AVR in power systems. Comput. Electr. Eng. 2018, 70, 261–274. [Google Scholar] [CrossRef]
  45. Khan, I.A.; Alghamdi, A.S.; Jumani, T.A.; Alamgir, A.; Awan, A.B.; Khidrani, A. Salp Swarm Optimization Algorithm-Based Fractional Order PID Controller for Dynamic Response and Stability Enhancement of an Automatic Voltage Regulator System. Electronics 2019, 8, 1472. [Google Scholar] [CrossRef]
  46. Pan, I.; Das, S. Chaotic multi-objective optimization based design of fractional order PIλDμ controller in AVR system. Int. J. Electr. Power Energy Syst. 2012, 43, 393–407. [Google Scholar] [CrossRef]
  47. Ahmed, E.M.; Elmelegi, A.; Shawky, A.; Aly, M.; Alhosaini, W.; Mohamed, E.A. Frequency Regulation of Electric Vehicle-Penetrated Power System Using MPA-Tuned New Combined Fractional Order Controllers. IEEE Access 2021, 9, 107548–107565. [Google Scholar] [CrossRef]
  48. Bakir, H.; Guvenc, U.; Kahraman, H.T.; Duman, S. Improved Lévy flight distribution algorithm with FDB-based guiding mechanism for AVR system optimal design. Comput. Ind. Eng. 2022, 168, 108032. [Google Scholar] [CrossRef]
  49. Zhao, W.; Wang, L.; Mirjalili, S. Artificial hummingbird algorithm: A new bio-inspired optimizer with its engineering applications. Comput. Methods Appl. Mech. Eng. 2022, 388, 114194. [Google Scholar] [CrossRef]
  50. Devan, P.A.M.; Hussin, F.A.; Ibrahim, R.B.; Bingi, K.; Nagarajapandian, M.; Assaad, M. An Arithmetic-Trigonometric Optimization Algorithm with Application for Control of Real-Time Pressure Process Plant. Sensors 2022, 22, 617. [Google Scholar] [CrossRef]
  51. Ahmad, M.F.; Isa, N.A.M.; Lim, W.H.; Ang, K.M. Differential evolution: A recent review based on state-of-the-art works. Alex. Eng. J. 2022, 61, 3831–3872. [Google Scholar] [CrossRef]
  52. Askari, Q.; Saeed, M.; Younas, I. Heap-based optimizer inspired by corporate rank hierarchy for global optimization. Expert Syst. Appl. 2020, 161, 113702. [Google Scholar] [CrossRef]
  53. Mirjalili, S. SCA: A Sine Cosine Algorithm for solving optimization problems. Knowl.-Based Syst. 2016, 96, 120–133. [Google Scholar] [CrossRef]
  54. Li, S.; Chen, H.; Wang, M.; Heidari, A.A.; Mirjalili, S. Slime mould algorithm: A new method for stochastic optimization. Future Gener. Comput. Syst. 2020, 111, 300–323. [Google Scholar] [CrossRef]
  55. Zhao, W.; Zhang, Z.; Wang, L. Manta ray foraging optimization: An effective bio-inspired optimizer for engineering applications. Eng. Appl. Artif. Intell. 2020, 87, 103300. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of generalized AVR components.
Figure 1. Schematic representation of generalized AVR components.
Fractalfract 07 00690 g001
Figure 2. Model TFs of AVR components.
Figure 2. Model TFs of AVR components.
Fractalfract 07 00690 g002
Figure 3. Step responses of AVR systems without controller at various loading conditions.
Figure 3. Step responses of AVR systems without controller at various loading conditions.
Fractalfract 07 00690 g003
Figure 4. Frequency responses of AVR systems without controller at various loading conditions.
Figure 4. Frequency responses of AVR systems without controller at various loading conditions.
Fractalfract 07 00690 g004
Figure 5. Root locus plots of AVR systems without controller at various loading conditions.
Figure 5. Root locus plots of AVR systems without controller at various loading conditions.
Fractalfract 07 00690 g005
Figure 6. Schematic structure of AVR control schemes. (a) PID controller; (b) FOPID controller.
Figure 6. Schematic structure of AVR control schemes. (a) PID controller; (b) FOPID controller.
Fractalfract 07 00690 g006
Figure 7. Proposed MPA-based AVR controller optimization.
Figure 7. Proposed MPA-based AVR controller optimization.
Fractalfract 07 00690 g007
Figure 8. Statistical analysis using ANOVA test.
Figure 8. Statistical analysis using ANOVA test.
Fractalfract 07 00690 g008
Figure 9. Tukey test of obtained results.
Figure 9. Tukey test of obtained results.
Fractalfract 07 00690 g009
Figure 10. Step response comparisons of the various obtained optimization results at full-load condition with K G = 1.
Figure 10. Step response comparisons of the various obtained optimization results at full-load condition with K G = 1.
Fractalfract 07 00690 g010
Figure 11. Step response comparisons of the proposed AVR controller with the existing designs in the literature at full-load condition with K G = 1.
Figure 11. Step response comparisons of the proposed AVR controller with the existing designs in the literature at full-load condition with K G = 1.
Fractalfract 07 00690 g011
Figure 12. Step response comparisons of the proposed AVR controller with the existing designs in the literature at no-load condition with K G = 0.7.
Figure 12. Step response comparisons of the proposed AVR controller with the existing designs in the literature at no-load condition with K G = 0.7.
Fractalfract 07 00690 g012
Figure 13. Conversion speed evaluation of MPA method.
Figure 13. Conversion speed evaluation of MPA method.
Fractalfract 07 00690 g013
Figure 14. Step responses with load disturbance of 0.3 p.u.
Figure 14. Step responses with load disturbance of 0.3 p.u.
Fractalfract 07 00690 g014
Figure 15. Performance comparison with PIDD2 with parameters from [11].
Figure 15. Performance comparison with PIDD2 with parameters from [11].
Fractalfract 07 00690 g015
Table 1. Summary of control schemes and optimization algorithms of AVR systems in the literature.
Table 1. Summary of control schemes and optimization algorithms of AVR systems in the literature.
Ref.ControllerOptimizer Algorithm
Ref. [16]PIDParticle Swarm Optimization (PSO)
Ref. [17]PIDArtificial Bee Colony (ABC)
Ref. [18]PIDBiography-Based Optimization (BBO)
Ref. [19]PIDTree-Seed Algorithm (TSA)
Ref. [20]PIDGrasshopper Optimization Algorithm (GOA)
Ref. [21]PIDPattern Search Algorithm (PSA)
Ref. [22]PIDWhale Optimization Algorithm (WOA)
Ref. [23]PIDImproved Whale Optimization Algorithm (IWOA)
Ref. [24]PIDGenetic Algorithm (GA)
Ref. [25]PIDCuckoo Search (CS) Algorithm
Ref. [26]PIDSine–Cosine Algorithm (SCA)
Ref. [27]PIDImproved Kidney Inspired Algorithm (IKA)
Ref. [17]PIDDifferential Evolution (DE)
Ref. [28]PIDContinuous FireFly Algorithm (CFA)
Ref. [29]PIDSymbiotic Organisms Search (SOS) Algorithm
Ref. [30]PIDSalp Swarm Algorithm (SSA)
Ref. [31]PIDBacterial Foraging Optimization Algorithm (BFOA)
Ref. [32]PIDGravitational Search Algorithm (GSA)
Ref. [33]PIDAnt Lion Optimizer (ALO)
Ref. [34]PIDLocal Unimodal Sampling (LUS) Algorithm
Ref. [35]PIDAnt Colony Optimizer with Nelder–Mead (ACO-NM)
Ref. [36]FOPIDParticle-Swarm Optimization (PSO)
Ref. [37]FOPIDGenetic Algorithm (GA)
Ref. [38]FOPIDArtificial-Bee Colony Optimizer (CNC-ABC)
Ref. [39]FOPIDChaotic Ant Swarm (CAS)
Ref. [40]FOPIDMulti-Objective Extremal Optimization (MOEO)
Ref. [41]FOPIDSine–Cosine Algorithm (SCA)
Ref. [42]FOPIDChaotic Yellow Saddle Goatfish Algorithm (CYSGA)
Ref. [43]FOPID (method 1)Improved Multi-Objective NSGA-II with Henon Map
Ref. [44]FOPID (method 2)Cuckoo Search (CS)
Ref. [37]FOPID (method 3)Particle Swarm Optimization (PSO)
Ref. [45]FOPID (method 4)Salp Swarm Optimization (SSO)
Ref. [46]FOPID (method 5)Multi-Objective NSGA-II with Chaotic Map
ProposedFOPIDMarine Predator Algorithm (MPA)
Table 2. AVR system’s responses at various loading conditions.
Table 2. AVR system’s responses at various loading conditions.
TypeParametersUnits K G = 0.7 K G = 0.8 K G = 0.9 K G = 1.0
 Step Resp.OvershotPercentage %50.495655.961161.018365.7226
Rise TimeSeconds0.31720.29440.27600.2607
Settling TimeSeconds4.89805.41496.42576.9865
Steady-State Err.Per Unit (p.u.)0.12200.11580.10240.0938
Freq. Resp.Gain MargindB4.192.110.07−2
Phase MarginDegree ( )14.66.70.2−5.3
Root LocusClosed Loop
System Poles
−12.49−12.31−12.13−11.93
−99.97−99.97−99.98−99.98
−0.52 + 4.66i−0.61 + 4.47i−0.7 + 4.25i−0.8 + 4.02i
−0.52 − 4.66i−0.61 − 4.47i−0.7 − 4.25i−0.8 − 4.02i
Table 3. Existing objective functions in the literature.
Table 3. Existing objective functions in the literature.
Objective FunctionTypeRef.
O F = I A E = | e v | d t SingleRef. [38]
O F 1 = I A E , O F 2 = 1000 | E s s | , O F 3 = T s MultipleRef. [40]
O F 1 = ω c f , O F 2 = P m MultipleRef. [43]
O F = I T A E = t . | e v | d t SingleRef. [45]
O F = ( ω 1 . O S ) 2 + ω 2 . T s 2 + ω 3 d V m a x 2 SingleRef. [37]
O F = ω 1 . O S + ω 2 . T r + ω 3 . T s + ω 4 . E s s + ( ω 5 . | e v | d t + ω 6 . V f ( t ) 2 ) d t + ω 7 P m + ω 8 G m SingleRef. [36]
O F = ( 1 e β ) . ( O S + E s s ) + e β . ( T s T r ) SingleRef. [39]
O F 1 = I T S E = t . e v 2 d t , O F 2 = t . Δ u 2 ( t ) d t , O F 3 = I T S E = t . e l o a d 2 d t MultipleRef. [46]
O F = ω 1 . O S + ω 2 . T s + ω 3 . E s s + ω 4 . | e v | d t + ω 5 . u 2 d t SingleRef. [37]
OS = overshoot, T r = rise time, T s = settling time, E s s = steady-state error, e v = error voltage, u = control signal, P m = phase margin, G m = gain margin, ω c f = gain crossover-frequency, e l o a d = error signal during load disturbance, d V m a x = maximum-point of voltage signal-derivative, ω 1 ω 8 = weighting factors.
Table 4. Obtained objective functions over 30 runs of the studied optimizers.
Table 4. Obtained objective functions over 30 runs of the studied optimizers.
Run No.AHAATOADEHBOSCASMAMRFOMPA
10.02940.03250.01750.01780.05320.05070.01680.0165
20.01750.03120.01640.01700.03420.04140.01560.0154
30.01920.03580.02330.02090.02550.03670.01620.0162
40.01590.02510.02050.02020.03680.05000.01550.0151
50.02650.02880.01800.02060.03260.04840.01650.0152
60.01960.05980.01520.02050.04200.04860.01780.0154
70.01910.02530.02320.01730.05890.02960.01760.0153
80.01660.02270.01620.02330.03660.03560.01590.0162
90.02570.02370.02800.01800.05840.05930.01740.0152
100.02070.02500.02030.03140.02150.02840.01560.0152
110.03200.02830.02020.02000.06630.05120.01670.0166
120.02640.03210.02050.02650.03740.05120.01590.0151
130.02090.02640.02050.02420.03780.03580.01550.0174
140.02870.03560.02020.01620.03760.04150.01590.0170
150.02710.02040.02550.04450.03880.05460.01550.0165
160.02310.03220.01710.01660.04360.05270.01550.0152
170.02190.03210.02080.02570.04700.04780.01550.0153
180.01960.02940.01780.01610.04220.02330.01830.0176
190.02730.03240.01850.02370.03040.05500.01800.0156
200.02230.02870.01840.02280.03630.04880.01750.0151
210.01570.02000.02230.02670.06110.03750.01550.0154
220.01780.03680.01700.01870.05580.02570.01570.0156
230.02800.03380.01570.03010.03270.04930.01530.0151
240.01700.02540.02380.02670.02980.05030.01540.0165
250.02220.02160.02700.01820.03080.05580.01790.0162
260.02200.02600.02040.02000.02850.04910.01550.0151
270.01870.03260.03430.03730.03460.03630.01600.0151
280.03050.03600.01990.01850.07630.03650.01700.0172
290.02450.03870.01950.02060.02800.03280.01550.0160
300.01980.02610.01670.02510.03720.03580.01550.0159
Table 5. Best FOPID parameters and statistical analysis of the studied optimizers over the performed 30 runs.
Table 5. Best FOPID parameters and statistical analysis of the studied optimizers over the performed 30 runs.
AHAATOADEHBOSCASMAMRFOMPA
Optimum FOPID Parameters
K p 1.82952.00001.70571.52151.85151.57911.65061.7061
K i 0.81521.00000.82630.72931.00000.94170.78780.8068
K d 0.39680.40000.40000.36010.40000.38960.39320.4000
λ 1.16361.13351.11581.12691.05841.02341.12351.1286
μ 1.25121.17971.21291.20731.19111.13441.20931.2164
Statistical Parameters
Best0.01570.02000.01520.01610.02150.02330.01530.0151
Worst0.03200.05980.03430.04450.07630.05930.01830.0176
Mean0.02250.03010.02050.02280.04110.04330.01630.0158
Median0.02200.02910.02020.02060.03730.04810.01590.0155
Std. Deviation0.00470.00760.00420.00650.01310.00980.00100.0008
Table 6. Statistical analysis based on 100 runs of MPA and MRFO algorithms.
Table 6. Statistical analysis based on 100 runs of MPA and MRFO algorithms.
AlgorithmBestWorstMeanMedianStd. Deviation
MRFO0.015345120.023111890.017048960.016410320.00168091
MPA0.015077090.018487930.015647080.015466130.00059616
Table 7. Performance comparisons using ANOVA test.
Table 7. Performance comparisons using ANOVA test.
ANOVA Table
SourceSSdfMSFProb > F
Columns0.0237570.0033966.72 4.53626 × 10 52
Error0.01182320.00005
Total0.03555239
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Noman, A.M.; Almutairi, S.Z.; Aly, M.; Alqahtani, M.H.; Aljumah, A.S.; Mohamed, E.A. A Marine-Predator-Algorithm-Based Optimum FOPID Controller for Enhancing the Stability and Transient Response of Automatic Voltage Regulators. Fractal Fract. 2023, 7, 690. https://doi.org/10.3390/fractalfract7090690

AMA Style

Noman AM, Almutairi SZ, Aly M, Alqahtani MH, Aljumah AS, Mohamed EA. A Marine-Predator-Algorithm-Based Optimum FOPID Controller for Enhancing the Stability and Transient Response of Automatic Voltage Regulators. Fractal and Fractional. 2023; 7(9):690. https://doi.org/10.3390/fractalfract7090690

Chicago/Turabian Style

Noman, Abdullah M., Sulaiman Z. Almutairi, Mokhtar Aly, Mohammed H. Alqahtani, Ali S. Aljumah, and Emad A. Mohamed. 2023. "A Marine-Predator-Algorithm-Based Optimum FOPID Controller for Enhancing the Stability and Transient Response of Automatic Voltage Regulators" Fractal and Fractional 7, no. 9: 690. https://doi.org/10.3390/fractalfract7090690

APA Style

Noman, A. M., Almutairi, S. Z., Aly, M., Alqahtani, M. H., Aljumah, A. S., & Mohamed, E. A. (2023). A Marine-Predator-Algorithm-Based Optimum FOPID Controller for Enhancing the Stability and Transient Response of Automatic Voltage Regulators. Fractal and Fractional, 7(9), 690. https://doi.org/10.3390/fractalfract7090690

Article Metrics

Back to TopTop