Simulation-Based Coyote Optimization Algorithm to Determine Gains of PI Controller for Enhancing the Performance of Solar PV Water-Pumping System

: In this study, a simulation-based coyote optimization algorithm (COA) to identify the gains of PI to ameliorate the water-pumping system performance fed from the photovoltaic system is presented. The aim is to develop a stand-alone water-pumping system powered by solar energy, i.e., without the need of electric power from the utility grid. The voltage of the DC bus was adopted as a good candidate to guarantee the extraction of the maximum power under partial shading conditions. In such a system, two proportional-integral (PI) controllers, at least, are necessary. The adjustment of (Proportional-Integral) controllers are always carried out by classical and tiresome trials and errors techniques which becomes a hard task and time-consuming. In order to overcome this problem, an optimization problem was reformulated and modeled under functional time-domain constraints, aiming at tuning these decision variables. For achieving the desired operational characteristics of the PV water-pumping system for both rotor speed and DC-link voltage, simultaneously, the proposed COA algorithm is adopted. It is carried out through resolving a multiobjective optimization problem employing the weighted-sum technique. Inspired on the Canis latrans species, the COA algorithm is successfully investigated to resolve such a problem by taking into account some constraints in terms of time-domain performance as well as producing the maximum power from the photovoltaic generation system. To assess the e ﬃ ciency of the suggested COA method, the classical Ziegler–Nichols and trial–error tuning methods for the DC-link voltage and rotor speed dynamics, were compared. The main outcomes ensured the e ﬀ ectiveness and superiority of the COA algorithm. Compared to the other reported techniques, it is superior in terms of convergence rapidity and solution qualities.


Introduction
Thanks to the huge availability of the sun supply in comparison with other energy sources, the solar PV systems seem one of the most encouraging applications of renewable energy today and in the forthcoming years. Moreover, the cost of the PV modules is reducing increasing a result of the progress in the manufacturing technology of the solar modules [1]. Standalone and grid-connected systems are the two main types of solar systems. Indeed, Sun availability is taking as one of the most important motivation for depending on the solar system as a promising solution. As being a type of renewable energy resources, the solar PV has gained a growing importance for uses in the electric power domain which has several advantages like (1) being harmless energy, (2) being suitable for isolated sites and (3) being cheap and easy maintenance.
Indeed, the efficiency and cost of the PV pumping system are not attractive yet to spread widely, but, many literature works are devoted to improving the total effectiveness and reducing the total cost of the solar pumping generation system [1,2].
A major drawback of using this sort of vector control lies in the fact that the effectiveness of the water-pumping system is in the heavy dependence on the appropriate adjustment of the PI parameters. Commonly, the set of these parameters is based on the classical and tiresome trials and error methods [3][4][5][6].
Adjustment method turned out to be complex and time-consuming. Moreover, the conventional PI tuning approach such as the symmetrical optimum [7], Ziegler-Nichols [8], Tyreus-Luyben [9] and Cohen-Coon [10] techniques requires that the designer has a good insight about the characteristics and dynamics of the controlled system. Hence, for the sake of overcoming these drawbacks, a systematic process to identify best gains of PI regulators has been provided in [3] and the modern optimization is considered as an auspicious solution.
In [11] the authors considered two optimization metaheuristics: Ant colony optimization (ACO) and differential evolution (DE). These strategies are reported to provide an optimized adjustment (PI) regulators in the direct torque control-space vector modulation (DTC-SVM) control loops, like rotor speed, electromagnetic torque, stator flux linkage and estimation of the linkage stator flux. Compared to other metaheuristics methods, those considered approaches require few tuning parameters to provide a fast convergence rate. Simulation and practical experimental on three phase induction motor (TIM) controlled by DTC-SVM are carried out. On the other hand, classical procedures such as frequency and root locus have been applied.
In [12], the study reported a novel optimization algorithm that combines all these three techniques: genetic algorithm (GA), artificial neural networks (ANN) and fuzzy logic, which was named as GNFPID. This technique aims to outperform the performance of the used PID controller. In [13], a genetic algorithm is applied as the adequate adjustment method of PI controllers applied to DTC drive for doubly fed asynchronous machine. Following in the same strategy, particle-swarm optimization (PSO), bat algorithm (BA) and genetic algorithm (GA) are considered for the adequate tuning of the gains of PI for brushless DC motor in [14].
Metaheuristics are an iterative processes for solving different optimization problems by taking into account simple trial operations [15,16]. Particle-swarm optimization [17], genetic algorithm [18], Imperialist competitive algorithm [19], ant colony optimization [20,21], differential evolution [22,23] are commonly treated. In [24], water cycle algorithm was adopted to adjust the parameters of PI controller for speed control and reactive power control loops in RSC (rotor side converter) and DC-link regulation voltage loop in GSC (grid side converter).
In the same way, [25] has referred to a thermal exchange optimization (TEO), this metaheuristic method was applied for the tuning of PI controller gains for the external loops in the conventional vector control strategy of a doubly fed induction generator based on wind turbine system. For assessing the efficiency and superiority of the proposed control strategy, other well-known metaheuristic procedures such as grasshopper optimization algorithm (GOA), particle-swarm optimization (PSO), Energies 2020, 13, 4473 3 of 17 water cycle algorithm (WCA) and harmony search algorithm (HSA), are taken into account for comparison purposes.
To best of our knowledge, for first time, the coyote optimization (COA) algorithm is used to identify the parameters of PI controller to improve the performance of solar PV water-pumping system under partial shading condition. Therefore, in this study, COA algorithm is investigated and further, to prove the superiority of the suggested COA algorithm, the obtained results are compared with the classical Ziegler-Nichols method, trial-error tuning method and PSO algorithm.
The remainder of this study is arranged as follows: the mathematical model of the solar photovoltaic water-pumping system (SPWPS) is described in Section 2. Section 3 is devoted to the formulation of the outer-loops PI controllers' design problem given as a constrained optimization problem. The proposed COA method is given in detail, aimed to resolve such a problem. In Section 4, the demonstrative results of the COA-tuned PI controllers tuning are carried out. Concluding remarks are summarized in Section 5.

Modeling of the Solar Photovoltaic Water-Pumping System
Referring to the solar pumping system, several control methods were developed to provide an effective functioning. These include speed control, distinct maximum power point tracking (MPPT) techniques, voltage control and frequency regulation under load variation [26].
This part studies and simulates a new structure of a solar pumping system. The overall structure of the proposed PV pumping system is described, as well as its various subsystems.
Vector control also called field-oriented control, is a method that permits independent control of the flux and torque components. It is equipped with high-performance drives. In this section, the interest is focused on the type of control.
The proposed system includes the following parts: PV panel, three-phase voltage inverter, DC-DC boost converter and induction engine-driven pump unit.
A schematic diagram of the new structure including indirect field oriented control (IFOC) and MPPT strategies, is illustrated in Figure 1. Concerning the control of this sort of system, different techniques were considered for effective functioning in terms of speed control, voltage control, distinct MPPT techniques and frequency regulation under load variations. The most important parts of this system are indicated in Figure 1 and are recapitulated as follows: (1) Shaded PV modules represent the power supply source of the induction engine via a three-phase voltage source inverter (VSI); (2) A hysteresis pulse-width modulation (PWM) based current control; (3) A boost converter for forcing the PV panel to operate at the MPP under different partial shading conditions [27]; (4) Flux weakening component is required to produce the reference current (I * ds ) and the speed controller output gives (I * qs ); (5) Vector control of asynchronous machine for driving centrifugal pump [28]. The reference speed is a measure of the generated PV power. The PI controllers' design and the suitable values of the gains regulators (Kps, Kis, Kpdc and Kidc) are typically tuned employing trials and errors technique, as mentioned before. Regarding the control design stage of VSI, this hard and challenging task turns into more difficult and time-consuming. A recent global meta-heuristic algorithm is an effective solution to tackle this matter.

Model of the Induction Machine
In this context, the Park transformation is adopted when modeling the asynchronous machine to omit the inductance variation with time [29].  The reference speed is a measure of the generated PV power. The PI controllers' design and the suitable values of the gains regulators (K ps , K is , K pdc and K idc ) are typically tuned employing trials and errors technique, as mentioned before. Regarding the control design stage of VSI, this hard and challenging task turns into more difficult and time-consuming. A recent global meta-heuristic algorithm is an effective solution to tackle this matter.

Model of the Induction Machine
In this context, the Park transformation is adopted when modeling the asynchronous machine to omit the inductance variation with time [29].
where σ = 1 − L m 2 L s L r V ds , V qs , V dr and V qr indicate the stator and the rotor voltages, respectively.
where τ r = L r R r In steady states conditions I ds = I mdr (8) The instantaneous value of the angle θ s is determined by this equation: The sliding speed ω sl can be calculated using this equation [29]: In steady states conditions (I mdr = I ds ) and thus ω sl = I qs τ r I ds (12) Equations (13) and (14) are needed to determine the reference flux and reference stator current, respectively [29]: where ϕ rn denotes nominal value of flux. ω mn denotes nominal value of mechanical speed L m : is the mutual inductance The correlation process between the torque and the speed of the pump model is described as follows [30]: The power of the pump has a cubic relationship with the motor speed w m , thus driving the centrifugal pumps [31,32].
To evaluate one of the reference speed components, an affinity law shall be applied for this purpose. Under different solar irradiation, the induction machine operates and seeks to determine the centrifugal pump flow rates. Moreover, to optimize the produced PV power, a MPPT control was investigated. It enables the extraction of the global MPP. The generated power serves to determine the first component of the reference speed as follows [33]: The PV array power is converted in terms of speed via the constant k.
The second component of the reference speed is evaluated by means of the DC-link voltage controller. In this system, a comparison between the measured bus voltage (V dc ) and the desired bus voltage (V * dc ), is carried out, which conducts to a voltage error (∆V dc ), described as [34]: The output speed can be obtained by the following equation [34]: where K pdc and K idc are, respectively the proportional and integral parameters of the DC-linkvoltage controller. The reference speed of the induction motor is obtained according to the following equation [33]: where ω * m is the reference speed and w m1 is the first part of reference speed, determined by means of Equation (17). w m2 is the second part of reference speed, determined by means of Equation (19).

Pump Model
The mechanical model of the induction motor actuating the centrifugal pump is provided by the following formula [35]: T r and T e represent the load torque and the electromagnetic torque, respectively, J is the moment of inertia and f is the damping coefficient which links the friction torque C f to the motor speed by The following equation used to represent the pump torque [33]: where A p denotes torque constant and estimates as following [33]: where P n and w mn indicate rated power and nominal value of mechanical speed, respectively.

Three-Phase Voltage Inverter
To achieve an optimum control of the power supply between the PV source and the asynchronous engine, a voltage source inverter is required as a crucial part. The employed inverter consists of three independent arms, as illustrated in Figure 2. DC-AC converters are dealing with in this section, which already applied in the solar energy sector as well as in the three-phase voltage inverters. The state of the switches (k n ) has a considerable impact on the inverter modeling. The produced voltages are provided by the following matrix form [36]: Energies 2020, 13, x FOR PEER REVIEW 8 of 18 The current provided by the DC source obeys to; A PWM hysteresis strategy is employed to control these power converters.

DC-DC Boost Converter
An ideal scheme of the DC-DC boost converter is illustrated in Figure 3. The relation of proportionality between the input voltage (E) and the output voltage, is described as follows: where u denotes the duty cycle. The DC converter can be integrated with MPPT strategy to extract the maximum power from the PV array. The current provided by the DC source obeys to; A PWM hysteresis strategy is employed to control these power converters.

DC-DC Boost Converter
An ideal scheme of the DC-DC boost converter is illustrated in Figure 3. The relation of proportionality between the input voltage (E) and the output voltage, is described as follows: where u denotes the duty cycle.

Hysteresis Based PWM Current Controller
In most studies, the PWM technique is the most common strategy since it perfectly adapts with the use of the voltage inverter. It is recommended also in case the vector control is adopted.
The principle of this command as presented in Figure 4 is to generate a PWM signal based directly on the variable to be controlled, through decisions based on the on-off control. In addition, it consists of maintaining the variation of the measuring currents within a band at a self-adjusting width, centered on the reference currents. The hysteresis comparator of any PWM structure receives as an input signal: the difference between the reference and measurement currents in order to produce at the output the control signal for the power switches [36]. In fact, this control is simple to operate as it does not require precise knowledge of the machine but suffers from the lack of control of the switching frequency.

COA Design for Optimal PI Parameters Tuning
The coyote optimization algorithm (COA) was recently proposed by Pierezan and Coelho [37]. The core idea of COA optimizer is based on Canis latrans species that reside mainly North America [37]. The algorithm is adapted to consider the social organization of the agents called coyotes and it has been serving with a different algorithmic construction. The greatest advantage of this method is to preserve a balance between the exploitation and exploration mechanisms during that optimizing process. By contrast with the gray wolf optimization (GWO), the COA does not center on how these dominant norms have been followed by these animals and the social hierarchy. Additionally, the COA depends not only on the hunting preys which take place in the GWO, but also on the social structure and regular experience interchange which are carried out by the coyotes. They are distinguished by cooperative functionalities as they head toward the prey in the close chain while they have a strong sense of smell which makes it possible to identify the location of the prey. Regarding the hunting process, coyotes attack in groups, this action forces the agents to update their positions to improve them. When coyotes' striking their opponents, they are well-prepared with a The DC converter can be integrated with MPPT strategy to extract the maximum power from the PV array.

Hysteresis Based PWM Current Controller
In most studies, the PWM technique is the most common strategy since it perfectly adapts with the use of the voltage inverter. It is recommended also in case the vector control is adopted.
The principle of this command as presented in Figure 4 is to generate a PWM signal based directly on the variable to be controlled, through decisions based on the on-off control. In addition, it consists of maintaining the variation of the measuring currents within a band at a self-adjusting width, centered on the reference currents. The hysteresis comparator of any PWM structure receives as an input signal: the difference between the reference and measurement currents in order to produce at the output the control signal for the power switches [36]. In fact, this control is simple to operate as it does not require precise knowledge of the machine but suffers from the lack of control of the switching frequency.

Hysteresis Based PWM Current Controller
In most studies, the PWM technique is the most common strategy since it perfectly adapts with the use of the voltage inverter. It is recommended also in case the vector control is adopted.
The principle of this command as presented in Figure 4 is to generate a PWM signal based directly on the variable to be controlled, through decisions based on the on-off control. In addition, it consists of maintaining the variation of the measuring currents within a band at a self-adjusting width, centered on the reference currents. The hysteresis comparator of any PWM structure receives as an input signal: the difference between the reference and measurement currents in order to produce at the output the control signal for the power switches [36]. In fact, this control is simple to operate as it does not require precise knowledge of the machine but suffers from the lack of control of the switching frequency.

COA Design for Optimal PI Parameters Tuning
The coyote optimization algorithm (COA) was recently proposed by Pierezan and Coelho [37]. The core idea of COA optimizer is based on Canis latrans species that reside mainly North America [37]. The algorithm is adapted to consider the social organization of the agents called coyotes and it has been serving with a different algorithmic construction. The greatest advantage of this method is to preserve a balance between the exploitation and exploration mechanisms during that optimizing process. By contrast with the gray wolf optimization (GWO), the COA does not center on how these dominant norms have been followed by these animals and the social hierarchy. Additionally, the COA depends not only on the hunting preys which take place in the GWO, but also on the social structure and regular experience interchange which are carried out by the coyotes. They are distinguished by cooperative functionalities as they head toward the prey in the close chain while they have a strong sense of smell which makes it possible to identify the location of the prey. Regarding the hunting process, coyotes attack in groups, this action forces the agents to update their positions to improve them. When coyotes' striking their opponents, they are well-prepared with a probability of threat and move away from their current position by an excessive random distance. Further, the COA method is initialized by an initializing population of coyotes whose size can be

COA Design for Optimal PI Parameters Tuning
The coyote optimization algorithm (COA) was recently proposed by Pierezan and Coelho [37]. The core idea of COA optimizer is based on Canis latrans species that reside mainly North America [37]. The algorithm is adapted to consider the social organization of the agents called coyotes and it has been serving with a different algorithmic construction. The greatest advantage of this method is to preserve a balance between the exploitation and exploration mechanisms during that optimizing process. By contrast with the gray wolf optimization (GWO), the COA does not center on how these dominant norms have been followed by these animals and the social hierarchy. Additionally, the COA depends not only on the hunting preys which take place in the GWO, but also on the social structure and regular experience interchange which are carried out by the coyotes. They are distinguished by cooperative functionalities as they head toward the prey in the close chain while they have a strong Energies 2020, 13, 4473 9 of 17 sense of smell which makes it possible to identify the location of the prey. Regarding the hunting process, coyotes attack in groups, this action forces the agents to update their positions to improve them. When coyotes' striking their opponents, they are well-prepared with a probability of threat and move away from their current position by an excessive random distance. Further, the COA method is initialized by an initializing population of coyotes whose size can be defined as the multiplication of packs and coyotes. The number of agents called coyotes in each pack is supposed to be equal and constant. In that regard, to facilitate the users' understanding, each coyote can be considered as a candidate solution for the optimization problem and its social condition presents the fitness function. The COA optimization algorithm started by a population of coyotes that are created randomly within the search space. More details about physical meaning and mathematical representation can be found in Ref [37,38]. Figure 5 presents the optimization process of COA.
Energies 2020, 13, x FOR PEER REVIEW 10 of 18 defined as the multiplication of packs and coyotes. The number of agents called coyotes in each pack is supposed to be equal and constant. In that regard, to facilitate the users' understanding, each coyote can be considered as a candidate solution for the optimization problem and its social condition presents the fitness function. The COA optimization algorithm started by a population of coyotes that are created randomly within the search space. More details about physical meaning and mathematical representation can be found in Ref [37,38]. Figure 5 presents the optimization process of COA.

Define input parameters and initialization
Estimate the corresponding fitness function  Regarding the PI design control, the suitable parameters of this type of controllers are usually tuned by trials and error techniques [3][4][5][6]. This non-systematic process turns into more complicated and into more time-consuming, particularly for the complex systems' design. Accordingly, the idea of transforming the identification process of these parameters (Kps, Kis, Kpdc and Kidc) to an optimization problem, is a promising alternative. By means of some advanced meta-heuristic algorithms, such optimization problems can be effectively managed [3]. Within the framework provided, two interactive PI regulators for speed and DC-link voltage loops are being introduced. These PI regulators are demanded to achieve an optimized value of this control using such COA algorithm. Figure 6 provides the suggested tuning process of PI controllers for the solar photovoltaic waterpumping system (SPVWPS) based on optimization approach. Regarding the PI design control, the suitable parameters of this type of controllers are usually tuned by trials and error techniques [3][4][5][6]. This non-systematic process turns into more complicated and into more time-consuming, particularly for the complex systems' design. Accordingly, the idea of transforming the identification process of these parameters (K ps , K is , K pdc and K idc ) to an optimization problem, is a promising alternative. By means of some advanced meta-heuristic algorithms, such optimization problems can be effectively managed [3]. Within the framework provided, two interactive PI regulators for speed and DC-link voltage loops are being introduced. These PI regulators are demanded to achieve an optimized value of this control using such COA algorithm. Figure 6 provides the suggested tuning process of PI controllers for the solar photovoltaic water-pumping system (SPVWPS) based on optimization approach.
During the optimization process, the gains of PI controllers for the DC-link voltage and speed loops are considered as the decision variables of the problem, these parameters are denoted as follows: , , T ps is pdc idc (27) where Kpdc and Kidc indicate the proportional and integral gains of DC-link regulator and Kps and Kis denote the proportional and integral gains of speed regulator. Under time-domain, such constraints like rise time, maximum overshoot, settling time and steady-state error of the system step-response, the performance criterion is minimized according to the suitable parameters of the problem [3].
The tuning process based on optimization is defined as follows: In formulating the optimization problem, cost function are outlined as follows: Which are adopted the Integral of Absolute error (IAE) and maximum overshoot (MO) criteria, is the initial search space. During the optimization process, the gains of PI controllers for the DC-link voltage and speed loops are considered as the decision variables of the problem, these parameters are denoted as follows: where K pdc and K idc indicate the proportional and integral gains of DC-link regulator and K ps and K is denote the proportional and integral gains of speed regulator.
Under time-domain, such constraints like rise time, maximum overshoot, settling time and steady-state error of the system step-response, the performance criterion is minimized according to the suitable parameters of the problem [3].
The tuning process based on optimization is defined as follows: In formulating the optimization problem, cost function are outlined as follows: Which are adopted the Integral of Absolute error (IAE) and maximum overshoot (MO) criteria, Under such circumstances and to show superior performance, some constraints had to be considered in the optimization framework, h j : R 4 + → R, j ∈ {1, 2} presented the inequality constraints which are considered for this optimization process.
The terms δ dc and δ speed are the overshoot of the DC-link voltage and rotor speed, respectively. The optimization process described in Equation (28) is considered as a multi-objective optimization type; an aggregate of all costs functions to form one single objective function is determined by means of an aggregation function [3].
where w i > 0, it denotes the weighting coefficients of the aggregation cost function. The objective function befitting for COA-tuned PI regulators for SPVWPS under shading condition, reformulated as follows:

Results and Discussion
The considered system in this research work contains a three-phase induction motor drive of 2200 W, 230 V, used to feed the water-pumping system. The main source is PV array with rating of 2400 W. The PV array composed of 234 PV module. The number of series PV modules interconnected is series is 18 modules whereas number of parallel string is 13 strings. The open circuit voltage and short circuit current of PV modules are 21.6 V and 0.64 A, respectively whereas the voltage and current at maximum power point are 17.6 V and 0.58 A. The case study considered the partial shading condition, three different level of the solar radiation; 1.0 kW/m 2 , 0.8 kW/m 2 and 0.5 kW/m 2 are subjected to the PV array.
In order to tackle the problem of time-consuming, the PI control tune by COA is performed, which further highlights the contribution of the proposed adjustment approach based on metaheuristic versus the given classical Ziegler-Nichols and Trial and Error method. Table 1 provides a summary of the obtained gains of the PI controllers for the classical based methods and the proposed optimization meta-heuristic called the coyote optimization algorithm (COA). Furthermore, a comparison is carried out between the obtained results for the proposed COA-tuned PI controllers and those obtained by the well-known Ziegler-Nichols and Trial and Error method for the PI controllers' design of DC-link voltage (Table 2 and Figure 7).   According to these results, COA-tuned PI controllers with the IAE criteria, outperforms these considered classical based methods in terms of settling time, overshoot and steady-state error indices. Indeed, the DC-link voltage performance using COA for designing the PI controller has a negligible steady-state error and also a negligible overshoot and it provides a shorter rise time than that in the case of a conventional adjustment method which presents a higher overshoot at the transient operation. For such reformulated PI design issues and in terms of offering the quality of solutions and non-premature convergence as well as fastness, the COA algorithm outperforms all other methods. To reinforce the effectiveness of the proposed COA method, the homologous PSO algorithm is considered for comparison purposes. Figure 8 represents the dynamic response of the DC-link voltage based on the PI parameters design by COA and PSO. It can be observed from this figure, that COA can maintain the DC-link voltage to constant value faster than that of the other meta-heuristic (PSO). This confirms the higher global convergence of COA due to the impact of exploration and exploitation mechanisms that are developed by that algorithm. From Figure 8, the tracking dynamic indicates a considerable oscillating using the PSO algorithm for adjustment PI controller design and the COA tuned PI controller enables us to reach rapidly the steady-state. To confirm the reliability of COA compared to PSO, both algorithms executed 30 times. The population size and maximum of iterations are 20 and 50, respectively. The results of statistical comparison are presented in Table 3. According to these results, COA-tuned PI controllers with the IAE criteria, outperforms these considered classical based methods in terms of settling time, overshoot and steady-state error indices. Indeed, the DC-link voltage performance using COA for designing the PI controller has a negligible steady-state error and also a negligible overshoot and it provides a shorter rise time than that in the case of a conventional adjustment method which presents a higher overshoot at the transient operation. For such reformulated PI design issues and in terms of offering the quality of solutions and non-premature convergence as well as fastness, the COA algorithm outperforms all other methods. To reinforce the effectiveness of the proposed COA method, the homologous PSO algorithm is considered for comparison purposes. Figure 8 represents the dynamic response of the DC-link voltage based on the PI parameters design by COA and PSO. It can be observed from this figure, that COA can maintain the DC-link voltage to constant value faster than that of the other meta-heuristic (PSO). This confirms the higher global convergence of COA due to the impact of exploration and exploitation mechanisms that are developed by that algorithm. From Figure 8, the tracking dynamic indicates a considerable oscillating using the PSO algorithm for adjustment PI controller design and the COA tuned PI controller enables us to reach rapidly the steady-state. To confirm the reliability of COA compared to PSO, both algorithms executed 30 times. The population size and maximum of iterations are 20 and 50, respectively. The results of statistical comparison are presented in Table 3.  The detailed performance (power, speed, DC-link voltage, load torque, electromagnetic torque, etc. ….) of the PV generation system, obtained with COA algorithm used for tuning PI controllers design of the DC-link voltage and rotor speed, are illustrated in Figure 9. Considering Figure 9, it can be observed that the desired and the detected values of the rotor speed converge simultaneously to the nominal value that demonstrate the ameliorated efficiency of the speed control. Considering the heavy dependence between the rotor speed components and the load torque, the water spilling out of the pump was reduced due to sudden torque variation in the  The detailed performance (power, speed, DC-link voltage, load torque, electromagnetic torque, etc. . . . .) of the PV generation system, obtained with COA algorithm used for tuning PI controllers design of the DC-link voltage and rotor speed, are illustrated in Figure 9.  The detailed performance (power, speed, DC-link voltage, load torque, electromagnetic torque, etc. ….) of the PV generation system, obtained with COA algorithm used for tuning PI controllers design of the DC-link voltage and rotor speed, are illustrated in Figure 9. Considering Figure 9, it can be observed that the desired and the detected values of the rotor speed converge simultaneously to the nominal value that demonstrate the ameliorated efficiency of the speed control. Considering the heavy dependence between the rotor speed components and the load torque, the water spilling out of the pump was reduced due to sudden torque variation in the Considering Figure 9, it can be observed that the desired and the detected values of the rotor speed converge simultaneously to the nominal value that demonstrate the ameliorated efficiency of the speed control. Considering the heavy dependence between the rotor speed components and the load torque, the water spilling out of the pump was reduced due to sudden torque variation in the suggested regulation method. Figure 9c illustrates the rotor speed behavior of the PV system due to the increased proportionally of the load torque, however, it achieves steady-state shortly at 0.2 s.
However, the proposed system uses an MPPT algorithm that provides the reference speed as an output signal, to be fed into indirect-field oriented control (IFOC) strategy. It is worth noting that the reference and measured values of the DC bus voltage converge towards the desired value, which has demonstrated an enhanced efficiency of the voltage regulation. The reference value (set-point) of the DC bus voltage must be constant and its value was selected to 388 V. The outer PI controller loop is adopted to keep the DC bus voltage constant. The actual value of the DC bus voltage is the value that is measured by the voltage sensor.
In other words, this Set point (reference speed), which must be applied to the machine, comes from a new maximum power tracking algorithm based on the PSO technique (PSO-MPPT) and a pump affinity law (Equation (16)). This process results in an optimization of photovoltaic power under partial shading conditions. Further, it can be seen that the load torque and the rotor speed follow the PV power. This is found in Figure 9. It was found that the inaccuracy of the speed has an important effect on the efficiency and reliability of the pump drive (The relationship of proportionality was indicated in Equation (16)).
In sum, the COA algorithm can identify the gains of PI controllers of DC link voltage and rotor speed in a reduced computational time and under functional constraints compared to the reported classical techniques which usually provide local solutions for the control problem with constraints.

Conclusions
The optimal gains of PI controllers were determined using coyote optimization algorithm (COA) in order to enhance the effectiveness of the solar photovoltaic water-pumping system. In fact, the controllers' adjustment is reformulated as a constraint optimization problem under the functional time domain. A comparative study of the dynamic response of the DC-link voltage and rotor speed for the COA-tuning PI controllers is carried out with the well-known classical Ziegler-Nichols and Trial and Error tuning procedures. Similarly, the proposed meta-heuristic COA succeeds in the tune of PI controllers for DC-link voltage and rotor speed more effectively than the considered meta-heuristic PSO algorithm. The COA provides good responses of the DC-link voltage and rotor speed, during the transient and steady state cycle. The main outcomes confirmed the performance superiority of the developed COA-tuned PI controllers, in terms of the rapidity of convergence and solutions quality. To reduce the cost of the controller, a senseless-MPPT will be considered in the future work.