Optimal design of fuzzy-PID controller for automatic generation control of multi-source interconnected power system

This paper suggests a fuzzy logic controller (FLC) structure from seven membership functions (MFs) and its input–output relationship rules to design a secondary controller to reduce load frequency control (LFC) issues. The FLC is coupled to a proportional–integral–derivative (PID) controller as the proposed FPID controller, which is tuned by an optimized water cycle algorithm (WCA). The proposed WCA: FPID scheme was implemented with two models from the literature under the integral time absolute error cost function. Initially, a two-area non-reheat unit was implemented, and the gains of PID and FPID controllers were adjusted to verify the suitability of WCA in solving LFC issues. Then, in order to identify the robustness of the closed-loop system, sensitivity analysis is carried out. Additionally, a two-area non-reheat unit was tested under the governor dead band nonlinearity. To guarantee the suitability of the proposed FPID controller, a model with a mixture of power plants, such as reheat, hydro, and gas unit in each area was carried out with and without the HVDC link, which can increase practical issues with LFC. The proposed controller's robustness was studied for all models under numerous scenarios, step load perturbations (SLP), and different objective functions. Simulation results proved that the proposed FPID controller provided superior performance compared to recently reported techniques in terms of peaks and settling time.


Introduction
Advanced power generation systems face significant challenges in order to match power generation with consumer demand. The power systems frequency is sensitive to fluctuations caused by consumer requests, making the frequency and tie-line power deviate from nominal performance. A significant deviation increases the possibility of overloading susceptible equipment, which poses a substantial threat to the stability of power systems. The primary controller is typically not adequate to return the power system to a steady state. Consequently, a secondary controller, LFC, or automatic generation control (AGC), is necessary to maintain system frequency and tie-line power within their nominal bounds during any step load perturbations (SLP) and ensure the connected generators maintain synchronization. The AGC output is the area control error (ACE), which is driven to zero to minimize the deviations in frequency and tie-line power [1,2]. Many approaches to address the concerns of LFC have been suggested, such as adaptive control, intelligent control, and optimal control [3][4][5]. Standard proportional-integral/proportional-integral-derivative (PI/PID) controllers are often used to address LFC issues because of their reliability, structural simplicity, and the satisfactory ratio between performance and reasonable cost [6]. Therefore, tuning PI/PID controllers are used in many optimization approaches, such as bacterial foraging optimization [7], differential evolution [8], firefly algorithm [9], quasi-oppositional harmony search [10], and teaching-learningbased optimization [11]. So far, due to the outcome of nonlinearity, realistic aspects of power system such as load disturbances and uncertainties such as governor dead band, generation rate control (GRC), and time delay, increase the chance of causing instability in power systems [12].
Conventional PID controllers are used at the explicit operating states. Due to the nonlinear behavior in real systems, PID controllers are not appropriate for fluctuating conditions; hence, the performance of the PID controller must be improved [13]. To overcome this weakness, the application of a fuzzy logic design is recognized as one of the most effective PI/PID controller design methods in the presence of nonlinearity [14][15][16] and including PV and wind turbine plants [17], making the FPID an expandable controller. The FLC was used to address LFC issues, to confirm that the FLC has a good dynamic response over control structures, such as optimal control, and traditional PI and PID controller [18,19]. The fuzzy plus PID (FPID) controller was designed for automatic LFC of a two-area non-reheat system to improve stability [20,21]. Previously, a FPID controller that was adjusted using the firefly algorithm (FA) for a multi-area multi-unit power system showed superior dynamic response, which was verified using traditional methods [22]. Also, for sophisticated interconnected power systems in a deregulated circumstance, including nonlinearity, a FPID controller with a derivative filter (FPIDF) was established in [23]. The design of an optimal fuzzy fractional-order PID controller for of a photovoltaic-reheat thermal is presented in [24]. More recently, for a better result, hybrid particle swarm optimization-pattern search (hPSO-PS) based on a fuzzy PI controller was implemented [25]. Additionally, an effective hybrid of the local unimodal sampling (LUS) algorithm combined with the teaching-learning-based optimization (TLBO) called LUS-TLBO [26] and a hybrid of the cuckoo optimization algorithm (COA) with a harmony search (HS) algorithm (HSCOA) [15] have been magnificently employed to adjust the gains of FPID controllers for LFC problems.
From the above discussion, the enhancement of the fuzzy plus fundamental controller performance comes from: 1. The structure of the fundamental controller such as PI, PID, PIDF, and fractional-order PID (FOPID). 2. The optimization algorithm (which plays a key role in the system enhancement), 3. The structure of fuzzy; including the membership functions (MFs) and the input-output relationship rules.
Many structures of the basic controller, optimization algorithms, and the hybrids are clearly explained in the literature. However, considering the key role in system enhancement of the MFs of input-output, its rules, and relative weights, there is little information present in the literature. There are many fuzzy structures with different numbers of MFs, that are designed to address LFC issues [2,21]. The most extensively used fuzzy structure is the five MFs, based on the prevalence in the literature [2,27]. Therefore, in this paper, design of an effective fuzzy structure in terms of MFs with enhanced rules is the hypothesized strategy in reducing LFC issues. Fuzzy structures with seven MFs and the adjusted rules plus a conventional PID (proposed FPID) controller to obtain the best dynamic performance of LFC problems are presented. As mentioned before, the optimization algorithm plays an essential role in the enhanced performance of power systems. According to the ''no free lunch'' principle, any meta-heuristic optimization algorithm might be more efficient than other algorithms in tackling a specific problem while not doing well in other problems [28]. To obtain a more desirable degree of smooth and damp fluctuations and owing to the continual development of integrated power generation, computational algorithm-based controllers give a significant motivation for research to design novel mixtures to discover the best suitable ones in LFC difficulties. One of the best optimization algorithms is the water cycle algorithm (WCA) [29], which is commonly used in different engineering challenges [30][31][32]. WCA is extensively used to combat LFC issues [1,[33][34][35]. Therefore, the optimization of the WCA was chosen to achieve the best improvements and scaling factors of the proposed FPID. Furthermore, the WCA-based FPID with five MFs (WCA: FPID) was also implemented to fairly compare its performance with the proposed FPID controller with seven MFs under the same environment. A common two-area multi-unit thermal-hydro-gas power source was also considered with and without the HVDC tie-line in test system 2 to demonstrate the effectiveness and robustness of the proposed FPID controller.
In summary, the main contributions of this study are: i. To propose a powerful scheme such as a WCA-based PID and FPID with five and seven MFs to resolve LFC issues in interconnected power systems. ii. To complete the sensitivity analysis by changing the system parameters and loading conditions. iii. To employ a two-area multi-unit system with/without an HVDC link, which is raised for the practical issues of LFC. iv. To carry out the dynamic transient responses of the interconnected power system. v. To compare the dynamic performance of the proposed FPID controller with recently reported strategies.

Power system under study
The system described herein is widely used in the literature to design and analyze the AGC of the interconnected power system. This is the first example of a system model approach that discusses and analyzes, in a straightforward manner, the limitations and issues arising from LFC. The transfer functions of test system 1 (two-area non-reheat thermal power plant) are presented in Fig. 1a [2,15]. Each non-reheat power plant consists of a speed governing system, turbine, and generator with capacity of 2000 MW , and a nominal preliminary loading of 1000 MW. The transfer functions (TF) of various blocks used in each nonreheat thermal power plant are carried out in [2,15]. The steam turbine is characterized by: Whenever the perturbation occurs, the frequency fluctuates, resulting in system instability. As a result, the controlling system detects these fluctuations and adjusts the turbine speed to achieve an optimum speed that suppresses the frequency fluctuations. The speed governor valve TF is: The general TF of the generator-load is presented by: where K PS and T PS are the power system gain and the power system time constant, respectively [11,15]. Each area has three inputs and two outputs. The inputs are the input of the controller DP ref , tie-line power deviation DP tie , and the load perturbation DP L and the outputs are the generator frequency deviation Df : The combination of frequency deviation and tie-line power error is known as the area control error (ACE). The ACEs are taken as input to the FPID controllers corresponding to each area. ACEs for the test system 1 as shown in Fig. 1a are: where DP tie is the change in tie-line power, Df 1 and Df 2 are the frequency deviations in area 1 and area 2, respectively, and B is the frequency bias factor, [36]. When the power system is subjected to a load perturbation, ACEs are used as a regulating signal to diminish DP tie and Df i to zero at the steady state. Thus, the inputs and outputs scaling factors of the FPID controller must be selected suitably to enhance the transient performance of the system. Additionally, the governor dead band (GDB) is designed in the considered system to add more challenge and to prove the robustness of the HHO/PD-PI in the presence of nonlinearity properties [15].

Controller structure
The PID is extensively used in industry because of its reliability, structural simplicity, and effective performance [37]. Therefore, employing FLC cascaded by the PID controller, referred to as a fuzzy-PID (FPID) controller, has confirmed superior reliability compared to a PID controller in overcoming the challenges in the growth of power systems [15,38]. Test system 1 being studied with the FPID was built in MATLAB/SIMULINK and incorporated the parameter values shown in Fig. 1b. The typical structure of the FPID controller is demonstrated in Fig. 2, where K 1 ; K 2 are the input scaling factors and K P ; K I; and K D are the output scaling factors. The FLC has four components.
(i) The fuzzifier which alters the value into fuzzy sets. (ii) The fuzzy inference system which performs all the logical manipulations. (iii) The rule base which contains the control rules and MFs. (iv) The fuzzy inference system's output (which is the fuzzy value that is transformed into real value by the defuzzification method) [16]. However, there are many MFs in the FLC; the triangular MFs are commonly accepted in FLC strategies because of their realtime applications, economical nature, and improved performance [39]. The fuzzy controller with scaling factors, where five linguistic variables of MFs for inputs and output of fuzzy are used, was planned in [40] and was used effectively used many times in the literature [2,15,16].  In this study, the proposed fuzzy controller structure consists of seven triangular MFs, as shown in Fig. 3. The fuzzy inputs are the ACE and the rate of change of ACE (ACE derivative). The outputs were converted into seven linguistic variables NL (Negative Large), NM (Negative Medium), NS (Negative Small), Z (Zero), PS (Positive Small), PM (Positive Medium), and PL (Positive Large). Identical MFs with a similar range are used for both the inputs and the FLC design's output to get simplicity, enhanced computational efficiency, a smaller amount of memory usage, and enhanced performance analysis [41,42]. The Mamdani fuzzy interface system with the centroid defuzzification technique was chosen in the proposed design. The proposed rules that relate the inputs and outputs of the FLC are given in Table 1. The fuzzy rules represent a significant function in the FLC performance [36]. Therefore, in this study, the input-output rules are studied extensively by examining the systems dynamic behavior and checking the rules suitability with four test systems. The surface of the fuzzy-PID controller is given in Fig. 4. A similar controller with the same gains is employed in area 1 and 2.
Designing an optimum FPID controller, the gains/factors must be appropriately tuned. The preferred dynamic response has minimum settling time, in seconds, with a trivial overshoot and undershoot, in hertz, when the system is subjected to adequate perturbation. In this study, a robust algorithm such as WCA is used to adjust the FPID controller to extract greater dynamic performance from the AGC-controlled FPID. The WCA will be discussed in the following section.

Design of objective function
The power systems frequency is sensitive to fluctuations caused by consumer requests, making the frequency and tie-line power ðDF; and DP tieÀline Þ deviate from nominal performance. AGC's desired control designations should include ability to return the frequency deviations and tieline error to predefined states, ensuring a satisfactory degree of stability for the closed-loop system, and providing a fast response with a rapidly dampened oscillation. A suitable objective function should be applied appropriately before using modern heuristic optimization-based controller schemes to achieve the best LFC performance.   [27]. Many objective functions have been applied to the AGC problems in the literature, such as the integral square error (ISE), integral of time multiplied absolute error (ITAE), integral absolute error (IAE), and integral time square error (ITSE) to restore nominal performance. Although, the ITSE-based controller decreased the contribution of large initial errors and emphasized errors later in the response [6]. This provides more time to settle compared to the ITAE. The ITAE criterion better decreased the settling time than IAE or ISE, and also decreased the peak overshoot but less-so than the ITSE criterion [6,43]. From the literature, the ITAE criterion is preferred in AGC studies. Thus, in this study, the ITAE criterion was used to adjust the scaling factors of the FPID controller. Furthermore, it was valuable to confirm the performance of the proposed controller under various performance criteria. Therefore, the corresponding value of the ITSE will be also calculated.
The ITAE and the ITSE objective functions are represented as: where DF 1 ; DF 2 are the power system frequency deviations; DP tieÀline is the incremental variation in tie-line power; t sim is the simulation time which is 10 s for test system 1 and 20 s for test system 2. The optimum performance has the lowest value of J.
Minimize J for FPID controller is subject to: where the min and max superscripts stand for the lower and the upper values of the corresponding parameter. The minimum and maximum values of the parameters are 0:0 and 2:0, respectively [2,15,36]. To examine the comparative evaluation between the proposed strategy and the published methods in the literature, the settling times and peak undershoots of the frequency deviation and tie-line power error results under different load demands have been considered. The parameters values point out the speed of the dynamic response profiles of the test systems in consideration.

Water cycle algorithm
The WCA is an advanced, nature-inspired heuristic scheme based on the water cycle procedure. The main idea of the WCA arises from the flow of rivers/streams to the sea. The procedures of the WCA are realizable, simple, and require a few parameters that are defined by the user [29]. To solve the optimization problem, the problem variables must be planned as a matrix of raindrops with dimensions N pop Â N var .
Var ¼ where and are the population size and number of design variables, respectively. The cost of the streams, i.e., objective function (ITAE), is calculated by The streams (N SR ) that provide the best cost are chosen as one sea and several rivers ðN R ) The streams flow toward the rivers or toward the sea as Allocating the streams to the rivers and/or to the sea depends on the flow intensity, as in the following equation [44]: where N sn is the number of streams. Figure 5 shows the WCA flowchart. Let us assume that d is the remoteness of the river and sea [44]. Let XðX 2 0; C Â d ð Þ C [ 1Þ is the river path flowing to the sea which is randomly connecting [45], C lies between 1 and 2; and d is the present distance among the river and the sea. Hence, the location will be updated according to the jth stream X jþ1 Stream and river X jþ1 River in the exploitation phase, which are characterized by [45] X jþ1 where a rand is a random number within 0 and 1, if the finest results achieved by the stream are better than the river to be connected, then it changes their position. Furthermore, a similar condition can be appropriate when the river's finest value is greater than the sea. Once the streams are greater than the river, the optimal solution is obtained, and they swap their positions. Evaporation is a very important factor, as it prevents fast convergence of the algorithm to local optima. From rivers and lakes, water evaporates; this process completes the cycle. The next pseudocode expresses whether a river flows to the sea. The d Max value is set constant and near zero. d Max is adaptively reduced, and its value at the ith iteration is evaluated by the following equation: Afterward evaporation, rain occurs. The new raindrops are found in different areas and create a new stream equivalent to the previous one. To specify the new streams of new areas, the following equation is applied: where LB and UB are the lower and the upper values of design variables. According to the finest/least value of the objective function (ITAE), the best stream is considered a new river, and the excess of the streams flow to the new The WCA has several benefits over other swarm systems, such as (i) avoiding entrapment at minima due to the balance of exploitation and exploration phases. (ii) Fast convergence properties, and (ii) Its adaptability, recently, the WCA is implemented many times in the interconnected power system, [45][46][47]. WCA's suitability in the presence of the nonlinearity was demonstrated in [48,49]. Also, the multi-objective WCA was applied [50,51]. Recently, WCA has shown improved performance compared to other algorithms such as the genetic algorithm, TLBO, particle swarm optimization, and differential evolution techniques when dealing with LFC [52,53].
The WCA scheme adjusts PID and the proposed FPID controllers' gains according to the ITAE criterion. The block diagram of the gain organizing controller-based ITAE is revealed in Fig. 6, where the power system's output is a combination of tie-line power error and the frequency deviation relating to the concerned area. This output is estimated using ITAE criterion. The optimizer WCA compares the obtained value of ITAE with the previous one to get the lowest ITAE value and generates a new gain for the controller parameters. It repeats this operation based on the number of iterations. The lowest obtained ITAE value corresponds to the best controller gains, which provide a trivial deviation. For WCA-based PID or FPID real-time implementation, the dominant parameter of metaheuristic algorithms is population size [1]. According to [54], ''large initial population sizes do not outperform small populations in terms of identifying the optimum solution,'' making the output results satisfied, which is suitable for real-time situations of LFC problems. Furthermore, because the number of PID or FPID controller gains is quite low, they do not require a significant number of iterations to get the best solution, allowing them to be tuned online very fast in response to the power system variations.

Simulation results and discussion
4.    Primarily, the WCA based on the PID controller utilizing the ITAE criterion was employed at the two identical non-reheat thermal power plants in test system 1, under a 10% SLP in area 1 to guarantee the suitability of WCA for the AGC issues. The values of the parameters of the WCAbased PID controller of the best run are N pop ¼ 50; d Max ¼ 10 À9 ; and N Sr ¼ 5 corresponding to the best value for ITAE ¼ 0:13065; and the gains of PID are To demonstrate the superiority of the WCA-based PID controller, the obtained results were compared to the other existing schemes like TLBO [11], DE [8], and FA [55], as shown in Table 2.
From Table 2, it is revealed that the minimum ITAE value was obtained with the WCA ( Table 2. Comparing our results and the results from the literature, the WCA had improved performance when dealing with LFC issues.

Test system 1 with 10% SLP at area 1
Initially, the WCA based on the common fuzzy five MFs in the literature (WCA: FPID) was designed [15,16,27] to fairly compare between the WCA: FPID and the WCA based on the proposed fuzzy seven MFs (referred to as proposed FPID). This was employed to accomplish the more efficient operation of the AGC issues of test system 1.
The optimal settings for all test system controller parameters are gathered in Table 3. The efficacy of the proposed FPID controller was compared to WCA: FPID, the COA into HS algorithm (HSCOA)-based FPID controller [15], foraging optimization algorithm (BFO)-based FPID controller [2], and the hybrid particle swarm optimizationpattern search (hPSO-PS) built fuzzy PI controller [25], as demonstrated in Table 4. These comparisons were made using the main indices; the objective functions like ITAE and ITSE, the settling time, and the peak undershoot. The proposed FPID controller had the lowest values of objective functions (ITAE = 0.0015, ITSE = 9.34 9 10 -7 ) compared to WCA: FPID of five MFs, and the recently published strategies such as the HSCOA: FPID controller, the BFO: FPID controller, and the hPSO-PS: FPI controller. Also, the settling time and the peak undershoot were largely reduced in both frequency and tie-line power deviations. Figure 7 shows the dynamic response at 10% SLP in area 1. As mentioned before, the ITAE and ITSE minimally conflicted in settling time and in oscillation magnitude. However, the HSCOA has a good ITAE value compared to WCA: FPID, that results in less settling time and a higher ITSE value that leads to an increased oscillation magnitude. Based on the data ( Fig. 7 and Table 4) the WCA: FPID has a similar performance to the HSCOA: FPID scheme. Therefore, there was no need to employ the WCA: FPID in further investigations. Direct comparison of the results can be done to the HSCOA: FPID and the other competitors.

Sensitivity analysis of test system 1
Evaluating the power system dynamic behavior based on variations of loading conditions and system parameters are done to identify the robustness of the closed-loop system. All the system parameters increase and decrease by 25% and 50% under DP L1 = 10% SLP at area 1 to document the stability of the system. Figure 8 shows the sensitivity analysis under variation of Tg. Table 5 displays the change in system parameters and many of the important indices like ITAE value, settling time, undershoot, and overshoot for the system transient responses, i.e., DF 1 ; DF 2 ; and DP tie . For example, the overshoots of DF 2 and DP tie remain at zero in all changes. It is clear that the effect of fluctuations on the ITAE criterion, settling time, undershoot, and overshoot acquired by the proposed FPID controller under variation of system parameters was very small. These trivial fluctuations make the proposed FPID controller robust and perform well against with up to 50% variations in the parameters.

Test system 1 with GDB nonlinearity
The presence of the nonlinearity in the interconnected power system leads to larger oscillations and lower performance of the system response. To ensure the suitability of the proposed techniques in the presence of nonlinearity, the influence of GDB is studied in test system 1. The nominal parameters of test system 1 including GDB are available in ''Appendix A.2'', [15]. The GDB nonlinearity is expressed as [15] To check the proposed scheme performance under GDB, a 1% SLP is applied to area 1, the obtained cost functions ITAE and ITSE are compared to other reported strategy HSCOA: FPID. Table 6 shows that the value of ITAE is obtained using the proposed FPID (ITAE ¼ 0:0030) compared to HSCOA: FPID (ITAE ¼ 0:0019). Although, the obtained ITAE using the proposed technique is larger than HSCOA: FPID, but the ITSE for the proposed technique (ITSE = 2.747 9 10 -6 ) is smaller than the HSCOA: FPID (ITSE = 5.107 9 10 -6 ) indicating that the calculation of many cost function is preferable for fair comparison as seen in Fig. 9. The settling time of the system response and the undershoots value are also revealed in Table 6. To further show their difference, the transient responses of the considered regulators the proposed FPID and HSCOA: FPID are checked at 2%SLP and 5% SLP in each area are shown in Table 7. Figure 10 shows the dynamic response at 5% SLP. It is noticed that the proposed FPID works efficiently to the compared scheme HSCOA: FPID in the presence of GDB nonlinearity as a practical challenge to the concerned power system.

Test system 2 with AC tie-lines
In order to add a more realistic and complex power system, the investigation was extended to a multi-source power system, as shown in Fig. 11a, [14,15,26,56]. Additionally, The SIMULINK file incorporating the parameter values is demonstrated in Fig. 11b. To investigate the multi-source power systems dynamic performance with AC tie-lines only, a 1% SLP was subjected to area 1. The gains of the proposed FPID controller utilized by WCA are depicted in Table 3, and the optimized parameters of WCA were N pop ¼ 100; d Max ¼ 10 À8 ; and N Sr ¼ 5. The transient response in terms of settling time and undershoot are represented in Table 8, and the dynamic responses are shown in Fig. 12. The most ideal values of the objective functions obtained by the proposed regulator demonstrate the suitability of the proposed approach compared to recently published approaches, such as HSCOA: FPID [15], LUS-TLBO: FPI [26], DE: FPI [26], and DE: PID [56].
To improve the investigation and differentiate between the regulators, different SLPs were applied in area 1 and area 2. The studied objective functions for this scenario are summarized in Table 9. The ITAE and ITSE values are slightly increased, which is confirmed in Figs. 13 and 14. From the above results, the proposed technique has been well done at high load perturbations, where the oscillation magnitude successfully reduced, and the transient response rapidly settled.

Test system 2 with AC-DC tie-lines
An extra study was supported by including the HVDC link in the test system 2, as shown in Fig. 11 [26,56]. The finest scaling factors obtained with WCA-based ITAE criterion for reheat thermal, hydro, and gas units are shown in Table 3. Table 10 demonstrates the time-domain analysis in terms of ITAE, ITSE, undershoot, and settling time. The deviations of frequency and tie-line power are shown in Fig. 15. The finest values of the objective functions for the proposed regulator at 1% SLP in area 1 were (ITAE = 0.0023, ITSE = 1.89 9 10 -7 ), which proves the effectiveness of the proposed FPID compared to the HSCOA: FPID [15].
Also, different SLPs were investigated, similar in AC tie-line as shown in Table 11. The deviations of frequency and tie-line power for 3% SLP in area 1 are shown in Fig. 16. From the above scenarios, it can be concluded that the main indices of the transient responses confirm the superiority of the proposed FPID scheme in all terms at different SLPs compared to recently published schemes. Also, the AC-DC link enhances the transient responses compared to AC only tie-lines.

Conclusions
In this study, an optimally tuned FPID controller with seven MFs using an ITAE criterion was calculated and its corresponding ITSE value using an optimized algorithm in a WCA power system was employed. Two models were selected, a two-area two-unit, and a two-area six-unit model to test the performance of our analysis. To have a more realistic LFC study and to demonstrate the suitability of the proposed FPID controller, the dynamic LFC response profiles of the studied test systems are compared to recently published schemes. A common system, test system 1, was initially investigated. The dynamic responses showed that the results obtained by the proposed FPID were superior to recently published schemes and the uncertainty of the system parameters was improved with the proposed FPID. Also, to authenticate the competence of the proposed FPID, the study was extended to a two-area six-unit power system with/without an HVDC link. The time-domain investigation shows the superior performance of the proposed FPID controller because it had smaller magnitudes of oscillation and less settling time in all investigated scenarios. The proposed FPID with a seven-MF controller showed superior performance for the LFC test systems and readaption was not needed in a wide range of deviations from the nominal conditions. In conclusion, the WAC based on the proposed FPID controller makes the LFC framework more robust and shows a more stable and better outcome in a wide variation of loadings and conditions. In the future, prior to establishing its robustness, the suggested approach will be tested on complex real-world applications such as LFC with integrated electric cars, wind, and PV systems with time delay nonlinearity.

Declarations
Conflict of interest The author declares no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.