A Robust Model Predictive Control for a Photovoltaic Pumping System Subject to Actuator Saturation Nonlinearity

: In this paper, a new robust model predictive control (RMPC) for uncertain nonlinear systems subject to actuator saturation is designed to regulate the terminal voltage of a photovoltaic generator (PVG) that feeds a DC motor-pump via a buck DC–DC converter. The considered system is a combination of a PVG-converter and DC motor-pump, which possesses nonlinear behavior along with under a saturating control signal highly dependent on the operation point and climate conditions of solar radiation and temperature. As a result, the control task is complex due to the nonlinearity of the system and its dependence on climate conditions. Based on the dead-zone property, the presented paper introduces a new RMPC technique to provide an innovative and efﬁcient solution to ensure the closed-loop system’s robust stability in the presence of actuator nonlinearity. In this paper, the nonlinear system is described in polytypic form, and an appropriate linear feedback control law is designed and used to minimize an inﬁnite horizon cost function under the framework of linear matrix inequalities (LMIs). Furthermore, sufﬁcient state-feedback control law conditions are synthesized to guarantee the robust stability of the closed-loop system in the presence of polytypic uncertainties. Simulation results are provided, in which the results illustrate the effectiveness of the proposed method.


Introduction
Over the past three decades, solar energy has become increasingly popular as one of the main renewable energy sources. Several strategies have been introduced to use sunlight as a source of energy in which sunlight can be converted into heat (solar-thermal energy conversion) [1], electricity (solar-photovoltaic energy conversion) [2], solar fuel (hydrogen) generated via photocatalytic water splitting [3,4], or sunlight chemicals, such as H 2 O 2 production [5]. It can also be used for CO 2 reduction [6] and ammonia synthesis [7]. The use of photovoltaic energy has grown tremendously in a wide range of applications, where many non-electrified villages and rural regions rely on solar energy as their primary source of electricity since electrification is difficult and expensive [8,9]. In addition, photovoltaic water-pumping systems for irrigation and water supply in remote areas have been widely implemented due to their unique features of ease of installation, environment friendliness, and low maintenance costs [10,11]. For the optimization of energy, PV water-pumping systems have to operate at their maximum power point (MPP). A directly coupled PV electromechanical system operates at the intersection of current-voltage curves of the PV array and DC motor-pump set [12,13]. PV water-pumping systems are usually designed to extract as much energy as possible from the solar resource, using a simple DC-DC power converter controlled by an maximum power point tracker (MPPT) algorithm [14,15]. In the literature, many of the MPPT algorithms have been applied for obtaining the maximal power from PV systems. Perturb and observe (P&O) and incremental conductance (INC) are the most widely used methods [16][17][18]. However, the simplest and fastest algorithm is the constant voltage control approach, in which the PV array is controlled to operate at a constant voltage equal to the MPP voltage of the array at the standard test condition (STC) provided by the manufacturer, by adjusting the duty ratio of the power converter [15]. The authors in [19,20] have shown that this technique has improved stability and low dependence on solar irradiation. Elgendy et al. [13] proved in a comparative investigation that it offers significantly better energy use efficiencies (up to about 91%) compared to directly connected systems without taking the effects of insolation and temperature variations on the MPP voltage into consideration.
Generally, the dynamic behavior of power converter systems can be described as a bilinear model under a saturating control signal. The classical control techniques applied to power converters usually do not take into account input saturation, which can severely degrade the performance of the closed-loop system, thus making the closed-loop system unstable, especially if the converter is subject to large perturbation. Furthermore, when a power converter is used in solar applications, the control system becomes more complex. Many studies have lately demonstrated that PVG characteristics have a substantial effect on the dynamic behavior of power converters [21]. Since the PVG dynamic resistance is both an environmental variable and operating point dependent, major changes in the PV system might compromise its stability [15].
Despite the nonlinearity and uncertainty of these systems, the linear (state or output) feedback control of power converters based on the small signal control theory proposed by Middlebrook and Cuk [22] is the most applicable approach at present. For example, the fuzzy-logic controller proposed in [23] and the PID controller proposed in [24] are based on small signal control theory. Such linear controllers are designed based on the model linearization at a certain operation point, and the obtained model can only be useful for small variations around that specific steady point [25,26]. To cope with the deficiency of these linear controllers, several nonlinear control strategies have been developed with stability analysis, including fuzzy-logic controllers [27,28], adaptive control [29,30], neuralnetwork-based control [31], sliding mode control (SMC) [32], and feedback linearization controllers [33]. However, these methods have two major limitations: The constraints on the input actuators are not considered when using these methods. In addition, the closed-loop performance in terms of robust stability and the presence of parameter uncertainties is not optimized. In recent years, more consideration has been given to studying intelligent and robust control structures, which can maintain nonlinear systems stability over a certain operating range.
Model predictive control (MPC) is a control strategy that offers attractive solutions for controlling constrained linear or nonlinear systems [34,35]. The ability to handle hard constraints on states/outputs and inputs and time-varying behaviors make the MPC method an effective technique for a wide range of practical applications [36]. MPC algorithms are widely based on dynamic models to predict system behaviors over a prediction horizon. The reliable solution obtained is based on the minimization at each time step of an upper bound of the worst-case infinite horizon quadratic function.
In general, it is difficult to construct an accurate model of a system in which uncertainty is frequently present. It is well known that the uncertainty related to systems can be represented in the form of parametric uncertainty or bounded disturbance regions. In the presence of a significant level of uncertainty, the control law designed by MPC based on a nominal model is suboptimal and can even be infeasible. Unfortunately, traditional MPC approaches fail to explicitly handle the plant-model uncertainty. As a result, numerous researchers have recently shown great interest in the RMPC of uncertain systems, where several tools, such as the LMI optimization approach, were developed to build RMPC algorithms [37][38][39][40].
In most practical control processes, actuator saturation is a common kind of nonlinearity. The classical control techniques that ignore actuator saturation may severely reduce the performance of a closed-loop system, and the closed-loop system may lose its stability, especially if subjected to certain large perturbations [41,42]. Hence, it is important to consider the control input constraints in the RMPC design. In the literature, stability analysis and the synthesis of RMPC strategies for linear systems with actuator saturation have been considered as a popular topic in the past few decades. The common feature of most RMPC strategies is that they deal with input nonlinearity by estimating the domain of attraction in the presence of actuator saturation. Kothare et al. [37] proposed a new robust centralized constrained MPC technique using an LMI-based optimization approach to guarantee the asymptotic stability of closed-loop systems. The proposed algorithm allows the incorporation of a large class of plant uncertainty descriptions and can guarantee the stability and robustness of the controlled system. The design in [37] also includes systematic treatment of input and state constraints. Casavola et al. [38] presented a scheduling min-max MPC algorithm for linear parameter varying (LPV) systems with polytopic uncertainty subject to input saturation. Cao and Lin [43] proposed an RMPC algorithm for LPV systems with polytopic uncertainty in the presence of actuator saturation. In the approach in [43], the actuator saturation constraint was characterized in terms of the convex hull of a group of auxiliary linear feedback laws and an actual linear feedback law. Huanga et al. [44] improved the design in [43] by taking the relative weighting between the auxiliary and actual feedback laws into account. Despite the improvements obtained, this method does not consider certain issues of system nonlinearity. Another disadvantage is that the real-time implementation of this approach suffers from a heavy computational burden due to the large number of inequalities when the saturating linear feedback law is expressed on a convex hull. This may lead to the intractability of the optimization problem, especially in the case of high-dimensional systems and fast sampling applications. Moreover, numerous MPC designs for linear systems with actuator saturation have been published in [45][46][47][48][49].
During the past decade, many MPC techniques for systems subject to actuator saturation have been studied. In [50], a model predictive tracking control algorithm for the flexible air-breathing hypersonic flight vehicle was proposed to handle the flight control problem with actuator constraints and input delays. A novel, fast MPC technique with actuator saturation for large-scale structures based on the explicit expression form of the Newmark-β method and parametric variational principle was introduced in [51]. In [52], the asymptotic stability of the finite horizon MPC of nonlinear systems with incremental input constraints was proposed and investigated. The authors in [53] presented an MPC approach for systems with a dead zone and saturation, where two different control strategies based on MPC were compared: the former uses hybrid MPC, while the latter is based on dead-zone inversion and standard MPC. A novel discrete-time sliding mode predictive control technique for a tethered satellite with saturated input was proposed and investigated in [54], in which a discrete-time auxiliary controller was involved in the discrete-time MPC scheme to achieve fast and more stable performance. However, in contrast to linear systems subject to actuator saturation, there are few significant MPC designs for nonlinear systems subject to actuator saturation.
In this work, the ability of the LMI approach is exploited to accommodate an MPCbased technique for the control of nonlinear systems subject to actuator saturation, through the dead-zone property. As the main contribution in this work, a new RMPC for uncertain nonlinear systems subject to actuator saturation is proposed to control a photovoltaic pumping system where a simulation study on the PV pumping system, which comprises a PVG, a voltage-mode-controlled DC-DC buck converter, and a DC motor-pump, is illustrated to evaluate the performance of the proposed controller. The objective of the control process is to keep the PVG voltage at the MPP voltage in the presence of PVG dynamic resistance uncertainty and atmospheric condition changes. The proposed RMPC algorithm innovatively uses certain ideas from the work in [55]. In the proposed approach, a saturation model based on dead-zone nonlinearity is used to handle the saturation of the control input. The controller design is characterized as an optimization problem of the worst-case performance objective function over an infinite moving horizon. In addition, at each time step, adequate state feedback that renders the closed-loop saturated nonlinear system globally asymptotically stable to the origin is obtained via linear matrix inequalities. The main advantages of the proposed approach, compared with other well-known RMPC techniques, are its ability to consider both actuator saturation and system nonlinearity and the reduction in conservativeness by avoiding a large number of inequalities when the actuator saturation constraint is characterized in terms of the convex hull; therefore, the proposed algorithm guarantees a lower computation time compared to other methods. The proposed method can be applied to deal with complex industrial systems, such as computer numerical control (CNC) machines, active magnetic bearing (AMB), robot manipulators, overhead cranes, and DC-DC power converters. Different scenarios are used in the simulation study to illustrate the effectiveness and applicability of the proposed approach to the PV pumping system.
The rest of this paper is organized as follows: Section 2 describes the mathematical formulation of the problem to be handled. Some basic concepts concerning MPC are introduced and the saturated RMPC scheme is proposed in Sections 3 and 4, respectively. Next, several simulation results are demonstrated under the Matlab environment in Section 5 to test the validity and effectiveness of our method. Finally, the paper is concluded in Section 6.

Problem Formulation
Consider a nonlinear discrete-time dynamical system subject to actuator saturation: where k is the current time instant and x ∈ R n is the state. The control input is bounded as u(k) = sat(u(k)). sat(u(k)) ∈ R m is a saturation function of the control input u(k) ∈ R m , which is defined as: where u lim is a control input limit and u lim = u max . u max is the maximum control input limit. f ∈ C 2 is a nonlinear function of states and control inputs f (0, 0 , and the dynamic system (Equation (1)) can be reformulated as a polytopic uncertain system: where the system matrices are affine functions of a parameter vector β of r parameters (β = (β 1 , . . . , β r )). Each uncertain parameter β j is bounded between a minimum and a maximum value β j and β j . j = {1, 2, . . . , r} Moreover, we assume that: where Ω denotes the convex hull and A j B j are vertices of the convex hull. Ω can be written as: f (x(k), u(k)) is the nonlinear term obtained by differing between the nominal nonlinear model and the linear part, which is given by: Assume that f is globally Lipschitz or at least locally Lipschitz in a region D including the origin with respect to x(k) and uniformly in u(k). Therefore, we have the following Lipschitz condition: where N ∈ R n×n is a Lipschitz constant matrix. For x 0 = 0, the Lipschitz inequality (Equation (8)) can be rewritten as: where I is the identity matrix. This paper aims to find a non-saturating linear-state feedback law u(k) = Hx(k), where a dead-zone nonlinearity approach is used to handle the saturation nonlinearity. Hence, the dynamic system (Equation (3)) is written as: where H is the state feedback gain. Let us introduce the following useful lemmas for later use.

Lemma 3.
Assume that F and E are vectors or matrices with appropriate dimensions. For any positive scalar α > 0, the following inequality holds:

Robust Model-Based Predictive Control Using LMIs
Let x(k + 1|k ) ∈ X be the predicted state of the plant at time k + 1, i ≥ 0 and u(k + 1|k ) ∈ U the future control move at time k + 1, i ≥ 0, in which X and U are compact subsets of R n and R m , respectively, and both of them contain the origin as an interior point. The RMPC aims to find an efficient state feedback control law u(k + i|k) = H(k)x(k + i|k) , i ≥ 0 by minimizing the following worst-case performance function: where S ∈ R n×n and R ∈ R m×m are positive definite states and control weights, respectively. The optimization problem (Equation (17)) can be formulated as follow: Let us introduce a quadratic Lyapunov function V(x) = x T Px, P > 0 of the state x(k|k) , with V(0) = 0, at sampling time k. Suppose the following robust stability condition is satisfied: (19) Summing this inequality from i = 0 to i = ∞, we get: With x(∞|k) = 0 or V(x(∞|k)) = 0, the upper bound of the cost function is: where γ is a positive scalar and is regarded as an upper bound of the objective in Equation (17): Thus, the control action of robust MPC at time k can be obtained by solving the following optimization problem: Applying Schur complements, the condition x(k|k) T Px(k|k ) ≤ γ in Equation (21) can be expressed equivalently as the LMI: We can postulate the following theorems to construct the state feedback matrix H in order to minimize the upper bound γ. (More details can be found in Kothare et al. [37]).

Linear Systems Subject to Actuator Saturation
In this section, for the case in which the nonlinear term is not considered, a saturated controller is designed to stabilize the input-saturated linear system (Equation (25)): where u(k) ∈ R m is the saturation function of the control input defined in Equation (2). With the dead-zone nonlinearity expression of a saturating linear feedback law, as described in Lemma 1, the first result of our work is given by the following theorem.

Theorem 1.
Consider the discrete-time linear system (Equation (25)) and assume x(k|k) is the measured state of x(k) at each sample time k. The state feedback control law u(k + i|k) = H(k)x(k + i|k) , i ≥ 0 that minimizes the upper bound γ on the infinite horizon quadratic performance index J(k) is given by H(k) = YQ −1 , where Q, Q > 0, and Y are the solutions to the following LMIs: Proof of Theorem 1. See Appendix A.

Nonlinear Systems Subject to Actuator Saturation
In this section, we will extend the preceding development to the case of nonlinear systems under actuator saturation. With Theorem 1 and the Lipschitz condition (Equation (9)), the following theorem is the second result of our paper.

Theorem 2.
Consider the discrete-time nonlinear system (Equation (10)) and assume x(k|k) is the measured state of on the infinite horizon quadratic performance x(k) at each sample time k. The state feedback control law u(k + i|k) = H(k)x(k + i|k) , i ≥ 0 that minimizes the upper bound γ on the infinite horizon quadratic performance index J(k) is given by H(k) = YQ −1 , where Q, Q > 0, and Y are the solutions to the following LMIs: Subject to LMI Actually, Equation (29) is for constructing an invariant ellipsoid, and Equations (30) and (31) are for guaranteeing robust stability. The symbol " * " depicts a symmetric structure.
Proof of Theorem 2. See Appendix A.

Example 1
To demonstrate the effectiveness of the proposed RMPC approach, we consider the discrete-time double integrator [58]: which is subject to the input constraint −1 ≤ u(k) ≤ 1. The discrete-time system starts from the initial condition x(0) = 10 −5 T . We use our saturated RMPC algorithm (Equation (29)) to cope with this problem and compare the results with the saturated RMPC algorithm presented by Cao and Lin [43] and Kothare et al. [37]. Figure 1 shows the state responses and control input of the system under the three different designs. These responses demonstrate that our RMPC guarantees better stabilization performance than that of the previous designs, where the proposed controller results in a smaller overshoot, a shorter rising time, and a faster time response.
Actually, Equation (29) is for constructing an invariant ellipsoid, and Equations (30) and (31) are for guaranteeing robust stability. The symbol " * " depicts a symmetric structure.
Proof of Theorem 2. See Appendix A. □

Example 1
To demonstrate the effectiveness of the proposed RMPC approach, we consider the discrete-time double integrator [58]: which is subject to the input constraint −1 ≤ ( ) ≤ 1. The discrete-time system starts from the initial condition (0) = [ 10 −5 ] . We use our saturated RMPC algorithm (Equation (29)) to cope with this problem and compare the results with the saturated RMPC algorithm presented by Cao and Lin [43] and Kothare et al. [37]. Figure 1 shows the state responses and control input of the system under the three different designs. These responses demonstrate that our RMPC guarantees better stabilization performance than that of the previous designs, where the proposed controller results in a smaller overshoot, a shorter rising time, and a faster time response.
Moreover, Figure 2 shows that the accumulated costs associated with the proposed RMPC algorithm is lower than those of the other designs.

Example 2. DC-DC Buck-Converter-Based PV Pumping System
The photovoltaic pumping system under consideration is shown in Figure 3. It comprises a PVG, a DC-DC buck converter, and a DC motor coupled to a centrifugal pump. Moreover, Figure 2 shows that the accumulated costs associated with the proposed RMPC algorithm is lower than those of the other designs.

Example 2. DC-DC Buck-Converter-Based PV Pumping System
The photovoltaic pumping system under consideration is shown in Figure 3. It comprises a PVG, a DC-DC buck converter, and a DC motor coupled to a centrifugal pump.

Example 2. DC-DC Buck-Converter-Based PV Pumping System
The photovoltaic pumping system under consideration is shown in Figure 3. It comprises a PVG, a DC-DC buck converter, and a DC motor coupled to a centrifugal pump.

PVG Model
Generally, a PVG consists of several solar cells, which are connected in series and parallel to achieve the required voltage and current. The single-diode circuit model in Figure 4a is the most widely used one to study the behavior of a PVG, in which the PVG can be modeled as a current source to model the photo-current in parallel with a diode , an intrinsic shunt resistance , and a series resistance (more details can be found in [59,60]). A PVG has a nonlinear voltage-current ( − ) characteristic (see Figure 5) given by Equation (32), where and are the output voltage and current of the PV module, respectively. The efficiency of PVG depends on the internal characteristics of the device and and on many ambient conditions, such as the solar irradiance level, temperature, and shaded condition.
is the reverse saturation current. = is the junction thermal voltage of the array with cells connected in series.
is the ideality factor of the PV cell. = 1.60217646 10 C is the electron charge. = 1.3806503 10 J/K is the Boltzmann constant. T is the temperature in degrees Kelvin.

Since
≫ , it is possible to assume that . Thus, the photo-current for any temperature and solar irradiation can be expressed by the following equation [61]: is the short-circuit current.
is the MPP current at STC. is the short-circuit current temperature coefficient. is the solar irradiance in W/m . = 1000 W/m is the solar irradiation reference at STC. = 25 °C (298 degrees Kelvin) is the PVG temperature reference at STC.
The cell reverse saturation current, which depends on temperature, is given by: is the reverse saturation current. = 1.12 eV is the band-gap energy of the semiconductor used in the cell.

PVG Model
Generally, a PVG consists of several solar cells, which are connected in series and parallel to achieve the required voltage and current. The single-diode circuit model in Figure 4a is the most widely used one to study the behavior of a PVG, in which the PVG can be modeled as a current source to model the photo-current I ph in parallel with a diode D, an intrinsic shunt resistance R P , and a series resistance R S (more details can be found in [59,60]). A PVG has a nonlinear voltage-current (V − I) characteristic (see Figure 5) given by Equation (32), where v pv and I PV are the output voltage and current of the PV module, respectively. The efficiency of PVG depends on the internal characteristics of the device R P and R S and on many ambient conditions, such as the solar irradiance level, temperature, and shaded condition.
I 0 is the reverse saturation current. V t = N s KT q is the junction thermal voltage of the array with N s cells connected in series. n is the ideality factor of the PV cell. q = 1.60217646 × 10 −19 C is the electron charge. K = 1.3806503 × 10 −23 J/K is the Boltzmann constant. T is the temperature in degrees Kelvin. Since R p R S , it is possible to assume that I sc ≈ I mpp . Thus, the photo-current for any temperature and solar irradiation can be expressed by the following equation [61]: I sc is the short-circuit current. I mpp is the MPP current at STC. K i is the short-circuit current temperature coefficient. G 0 is the solar irradiance in W/m 2 . G re f = 1000 W/m 2 is the solar irradiation reference at STC. T re f = 25 • C (298 degrees Kelvin) is the PVG temperature reference at STC. The cell reverse saturation current, which depends on temperature, is given by: I rs is the reverse saturation current. E g = 1.12 eV is the band-gap energy of the semiconductor used in the cell.
I rs = I sc,n exp V oc,n nV t,n − 1 I sc,n is the cell short-circuit current at STC. V oc,n is the open-circuit voltage at STC. V t,n is the junction thermal voltage with N s cells connected in series at T re f . A modification on the reverse saturation current is proposed in [61] to match the open-circuit voltages of the model with the experimental data (from the manufacturer datasheet) for a wide range of temperatures K v is the open-circuit voltage temperature coefficient. The relation between R p and R S can be determined by solving Equation (33) for R S [61]: P max is the maximum power from the manufacturer datasheet. Then, R p can be derived iteratively as: To simplify the PVG model, a Thévenin equivalent circuit shown in Figure 4b can be used, with: where V TH is Thévenin equivalent voltage and R TH is Thévenin resistance. It can be observed that the Thévenin equivalent circuit parameters are not constant but depend on the environmental variables and operating point.
, is the cell short-circuit current at STC.
, is the open-circuit voltage at STC. , is the junction thermal voltage with cells connected in series at . A modification on the reverse saturation current is proposed in [61] to match the open-circuit voltages of the model with the experimental data (from the manufacturer datasheet) for a wide range of temperatures is the open-circuit voltage temperature coefficient. The relation between and can be determined by solving Equation (33) for [61]: is the maximum power from the manufacturer datasheet. Then, can be derived iteratively as: To simplify the PVG model, a Thévenin equivalent circuit shown in Figure 4b can be used, with: where is Thévenin equivalent voltage and is Thévenin resistance. It can be observed that the Thévenin equivalent circuit parameters are not constant but depend on the environmental variables and operating point.

DC Motor-Pump Model
The DC motor in Figure 3 is a DC motor with a permanent magnet (see Figure 6). This model is defined by an electrical circuit equation armature:

DC Motor-Pump Model
The DC motor in Figure 3 is a DC motor with a permanent magnet (see Figure 6). This model is defined by an electrical circuit equation armature: where e(t) is the counter electromotive force, i a is the armature current, v a is the armature voltage, and R a and L a represent the armature resistance and inductance, respectively.

DC Motor-Pump Model
The DC motor in Figure 3 is a DC motor with a permanent magnet (see Figure 6). This model is defined by an electrical circuit equation armature: where ( ) is the counter electromotive force, is the armature current, is the armature voltage, and and represent the armature resistance and inductance, respectively.  Figure 5 shows the schematic circuit diagram of a PV pumping system formed after connecting a DC motor to a DC-DC buck converter fed from a photovoltaic generator. The converter is assumed to operate in continuous conduction mode (CCM), and the inductor current is always larger than 0. Control is implemented via a pulse width modulation (PWM) approach.

Buck Converter Model
is the photovoltaic array voltage, which must be controlled through variation of the duty cycle ( ) in order to keep the array operation at the maximum power point. The diode is on inverse polarization, while represents the capacitor value, and is a power MOSFET controlled by a binary signal ( ); see Figure 7. The binary signal that triggers the switches on and off is controlled by a fixed-frequency pulse width modulation (PWM) circuit ( Figure 7). 1/ is the constant switching frequency of the PWM circuit, and is the switching period given by:  Figure 5 shows the schematic circuit diagram of a PV pumping system formed after connecting a DC motor to a DC-DC buck converter fed from a photovoltaic generator. The converter is assumed to operate in continuous conduction mode (CCM), and the inductor current is always larger than 0. Control is implemented via a pulse width modulation (PWM) approach. v pv is the photovoltaic array voltage, which must be controlled through variation of the duty cycle d(t) in order to keep the array operation at the maximum power point. The diode D is on inverse polarization, while C represents the capacitor value, and sw is a power MOSFET controlled by a binary signal S b (t); see Figure 7. The binary signal that triggers the switches on and off is controlled by a fixed-frequency pulse width modulation (PWM) circuit ( Figure 7). 1/T s is the constant switching frequency of the PWM circuit, and T s is the switching period given by: The average circuit model of the PV pumping system is shown in Figure 8. The corresponding discrete-time model is obtained by using the forward Euler approximation as follows: The total instantaneous quantities of the PV system are presented as the sum of the DC and AC components.

Buck Converter Model
where x = i a (k), v pv (k) T ∈ R 2 is the state, d(k) ∈ R is the saturated control input of the buck converter, X and D represent the equilibrium values, and x and d are the perturbed values of state and input. Using the same concept, we can obtain: Equation (46) is substituted into Equation (44), and a small-signal model is derived as follows: Equation (47) can be expressed as follows: where f is the nonlinear term that is obtained by differing between the nominal nonlinear model and the linear part, given by: The control input is subject to non-symmetric actuator saturation, and we can use the method proposed in [25] to transform it to a symmetric saturation. Thus, the saturation function of the control input is rewritten as follows: which can be expressed as If we add HX to Equation (48), we can obtain the following equality: Since the steady-state part is equal to 0, there is: In our case, the steady-state control signal is HX = 0.5, and the operating point can be calculated from: The system described in Equation (48) can be written as follows: where sat δ(k) = sat (−HX,1−HX) (H x(k)) is the new control input. Lemma 1 can be used to handle saturation nonlinearity. Thus, the new symmetric saturation model is constrained as in [62]: Based on this amplitude saturation, the decentralized dead-zone nonlinearity can be defined as: The closed-loop system is rewritten as:

Uncertainty Polytope Model
In this example, the Thévenin resistance R TH at the operating point is considered as the uncertain parameter. The discrete state-space model described in Equation (58) is expressed as: We consider r = 2 and the vector β[1/R TH ], in which: Photovoltaic generators are neither constant voltage sources nor constant current sources. In practice, the PVG is forced to operate at the boundaries of the constant current and constant voltage modes, if an MPPT is used, and the possessing dynamic resistance is equal to the resistance of its load at the maximum power transfer [63]. Based on the PVG equivalent circuit of Figure 4b, we can define a range of changes for Thévenin resistance. Under the open-circuit condition, R D is low, dominating the parallel connection with R P [15]. Therefore, we have: bounded by: With short-circuit and reference conditions, R D is high and R P dominates the parallel connection. There is: Note that R S is constant and R P is irradiation dependent [64]: where G 0 is the solar irradiance and R P,re f is the shunt resistance at STC (the solar irradiation is G re f = 1000 W/m 2 , and the PVG temperature is T re f = 25 • C). As an approximation, we have: Since the PV system matrix A(β) depends linearly on the uncertain parameter 1/R TH , we can define a polytope of r = 2 vertices, which contains all the possible values of the uncertain matrix. The uncertain matrices A j are:

Simulations Results
In this section, several simulation tests are implemented to verify the performance of the proposed control law and its robustness according to the presence of the input constraint, insolation variations, temperature variations, and dynamic resistance uncertainty. The results are compared to the RMPC method for nonlinear systems dependent on the Lipschitz bound presented in Poursafar et al. [55] and the saturated RMPC algorithm presented in Cao and Lin [43] designed by the linearized model around the equilibrium point. The parameters of the DC-DC converter and the PVG can be found in Tables 1 and 2, respectively.  Figure 9 depicts the computed input and state responses under the three different designs, where the waveforms are the PV voltage v pv , armature current i a , and duty cycle d. From Figure 9, it can be discovered that the proposed method outperforms the saturated RMPC algorithm in Cao and Lin [43] with a fairly good time response and lower fluctuations. Moreover, the PV voltage response settles to its desired value without any overshoot, whereas the approach in [55] has some unstable transient responses. In other words, our RMPC is well capable of handling the hard actuator saturation constraint. • Scenario 2: The insolation variations effect In the second case, we assume that the saturation limit is δ max = 0.5. An examination is made for various irradiances, such as G 0 = 200 W/m 2 , G 0 = 400 W/m 2 , and G 0 = 1000 W/m 2 . With the proposed control method, the power regulation response is shown in Figure 10. The increasing and decreasing nature of PV power with respect to insolation G 0 can be observed in the waveforms of PV power, which verifies the PV generator voltage at the MPP tracking. Similarly, the proposed RMPC also shows a good time response, low oscillation, and desired stability. • Scenario 3: The temperature variations effect In this scenario, we consider varying temperature with constant insolation G 0 = 1000 W/m 2 . The temperature changes between 100 • K and 298 • K at t = 100 s. Next, the temperature increases to 400 • K at t = 200 s, and it returns to 100 • K at t = 300 s, For this scenario, the saturation limit is δ max = 0.5. Figure 11 depicts the computed state responses under the proposed RMPC design. It is easily seen that the proposed RMPC guarantees the desired stability, with a good time response and low oscillation. • Scenario 4: The dynamic resistance variations effect In the last case, we consider dynamic resistance uncertainty in order to validate the robustness of the proposed controller. Figure 12 illustrates that our controller can stabilize the system on the reference PV voltage in presence of an abrupt dynamic resistance change. The overshoots and long settling time at t = 50 s and t = 120 s are the results of the aggressive move in the converter parameters, which demonstrates the effectiveness of our MPC design algorithm.  In this scenario, we consider varying temperature with constant insolation 0 = 1000 W/m 2 . The temperature changes between 100 °K and 298 °K at = 100 s . Next, the temperature increases to 400 °K at = 200 s, and it returns to 100 °K at = 300 s, For this scenario, the saturation limit is = 0.5. Figure 11 depicts the computed state responses under the proposed RMPC design. It is easily seen that the proposed RMPC guarantees the desired stability, with a good time response and low oscillation.   •

Scenario 4: The dynamic resistance variations effect
In the last case, we consider dynamic resistance uncertainty in order to validate the robustness of the proposed controller. Figure 12 illustrates that our controller can stabilize the system on the reference PV voltage in presence of an abrupt dynamic resistance change. The overshoots and long settling time at = 50 s and = 120 s are the results of the aggressive move in the converter parameters, which demonstrates the effectiveness of our MPC design algorithm.

Conclusions
This paper proposes a new RMPC design framework based on Lyapunov stability theory for a PV pumping system with actuator saturation and uncertain parameters. The input saturation effect is expressed by a model using dead-zone nonlinearity. At each time instant, the state feedback control law is obtained by minimizing the upper bound of the

Conclusions
This paper proposes a new RMPC design framework based on Lyapunov stability theory for a PV pumping system with actuator saturation and uncertain parameters. The input saturation effect is expressed by a model using dead-zone nonlinearity. At each time instant, the state feedback control law is obtained by minimizing the upper bound of the infinite horizon cost function within the framework of LMIs, which directly incorporates the input saturation. The internal stability of the closed-loop system is guaranteed in the sense of Lyapunov stability. The design of the closed-loop stabilization for a nonlinear system for the requirement of desired output responses with uncertain parameters is also investigated here. The simulation results obtained demonstrate the effectiveness and efficiency of our method. Furthermore, the novel method is examined in case the uncertainty is in the form of solar irradiance change and dynamic resistance, which proves that it is able to stabilize the system at the desired PV voltage.