Output-Feedback Nonlinear Adaptive Control Strategy of the Single-Phase Grid-Connected Photovoltaic System

This paper addresses the problem of controlling the single-phase grid connected to the photovoltaic system through a full bridge inverter with LCL-filter. The control aims are threefold: (i) imposing the voltage in the output of PV panel to track a reference provided by theMPPT block; (ii) regulating the DC-link voltage to guarantee the power exchange between the source and AC grid; (iii) ensuring a satisfactory power factor correction (PFC).The problem is dealt with using a cascade nonlinear adaptive controller that is developed making use of sliding-mode technique and observers in order to estimate the state variables and grid parameters, by measuring only the grid current, PV voltage, and the DC bus voltage. The control problem addressed by this work involves several difficulties, including the uncertainty of some parameters of the system and the numerous state variables are inaccessible to measurements. The results are confirmed by simulation under MATLAB\Simulink\SimPowerSystems, which show that the proposed regulator is robust with respect to climate changes.


Introduction
The global concern about climate change and the growing energy demand of industrialized countries have necessarily led to exploring other new sources like the renewable energy.The main advantages of this type of renewable energies reside in the reduction of pollution caused by the production of greenhouse gases.Among different types of these energies, the photovoltaic energy has obtained a great attention.
The photovoltaic energy systems are classified according to their use.The two principal classifications are gridconnected systems and stand-alone systems.The first one is connected to the grid through a three-phase or single-phase inverter; this category is used to deliver the power directly to utility grid and must be properly controlled according to power electrical legislations.The second one is used with a battery bank for electrifying remote rural areas.
The power factor correction, the DC output voltage regulation, and the maximization of the power provided by the PV modules are the main control objectives for allowing high power quality to the grid.To meet these requirements, various control methods have been proposed in [1][2][3].Indeed, [4] proposed a passivity based control.However, this technique requires being in a passive state.In [5], the proportional resonant (PR) controller is designed.The latter provides acceptable dynamic performance and eliminates the steady-state error.In [6,7] a sliding-mode controller (SMC) is used to have an excellent robustness and a very good steady-state performance as well as a fast-dynamic response.In [8] the authors propose a control law based on backstepping [9] and Lyapunov function so as to stabilize the global system.
In the literature, a few papers dealt with state-feedback control [10] and state observer design at the same time.In [11] a state-feedback control law is combined with observer to enhance disturbance rejection capability of a grid-connected photovoltaic inverter.In order to achieve the third goal, it is important that the PV system operates at its optimal power point, and for this task a maximum power point tracking is required.
In this work, we seek a control strategy that meets the following three control objectives simultaneously: (i) Perfect power factor correction (PFC): the grid phase currents and its corresponding voltages must be in phase.(ii) DC output voltage regulation: this voltage must be tightly regulated to a constant reference value to ensure the power exchange between AC grid and the DC bus.
(iii) Maximization of the power provided by the PV models.
To achieve the above objectives, a cascaded nonlinear adaptive controller is designed.The latter is constituted by a PV voltage loop and grid current loop.The first one is designed to extract the maximum power from the PV array by regulating the voltage provided by the PV generator.The second one includes the inner loop and aims to regulate the grid current to meet the PFC, and the outer loop is intended to enhance the power exchange, between the source and the grid, by regulating the DC-link voltage.These loops are designed based on sliding-mode technique combined with a Luenberger and extended Kalman filter type.Compared to previous works, the contribution of the present study enjoys several interesting features including the following: (i) Several control objectives are simultaneously considered (MPPT, DC Regulation, and PFC) while only some of these objectives have been tackled in previous works [12,13].
(ii) The nonlinearity of the controlled system was preserved [14] in order to keep all the properties of the studied system.
(iii) The grid voltage is not accessible to measurement and the internal impedance is assumed to be unknown, unlike previous works which assumed that voltage is available and the grid impedance is null or known [15][16][17].(iv) The present nonlinear adaptive control system does not necessitate many sensors for the measurement of some needed variables unlike previous works [18].
The paper is structured as follows: in Section 2, a mathematical model and description of all system stage are described.In Section 3, the design of the cascade nonlinear adaptive controller is presented.Section 4 covers the simulation results and discussion about the results.The conclusion is in the end.

System Description
This section describes the modelling of photovoltaic system connected to the grid.The power plant under study is shown in Figure 1.It consists of a PV panel, a DC-DC boost converter that drains the energy from the photovoltaic module and feed the DC bus capacitor, and a full bridge single-phase inverter with LCL-filter used at the output of the converter to achieve a satisfactory total harmonic distortion of the injected current.
By analyzing the circuit and applying the well-known Kirchhoff laws, the equations describing the dynamics of the system of Figure 1 are given below: or V pv and  pv are the voltage and the current generated by PV array, V dc is DC-link voltage,   and V  are the current and voltage of the grid,  0 designates the input current chopper, and  0 and  1 are the switching functions given by (2) The above instantaneous model (1a)-(1f) cannot be used directly for controller design as it involves the binary inputs, namely,  0 and  1 .To overcome this problem, let us use the averaging model (3a)-(3f).The state variables   , V  ,  1 , V pv ,  0 , and V dc are replaced by their average values  1 ,  2 ,  3 ,  4 ,  5 , and  6 over a cutting period.The control inputs  0 and  1 denote the average values of  0 and  1 , respectively.
The supply net voltage V  is considered inaccessible to measurement.Equation (3a) is completed with the internal model of the grid voltage signal V  =   sin(  ).In particular, where   and   denote, respectively, the amplitude and the angular frequency of V  .The inaccessible states and unknown parameters are presented in Table 1.

Control Strategies of the System
In this section, an output-feedback nonlinear controller will be synthetized.As represented by the averaging model (3a)-(3f) the system has two control inputs ( 0 ,  1 ).The controllers (Figure 2) will be designed to achieve the three main objectives mentioned previously.In addition, the observers are designed to estimate the values of unmeasurable states.The first task is dedicated to the design of the observers and the second task is devoted to development of an output-feedback nonlinear controllers.

State Observer Design.
The purpose of the present subsection is to design the observers, which provide accurate estimates of states variables and use them later to develop an output-feedback controller such that the estimation errors converge to zero.For that, an adaptive observer and Luenberger observer [19] are designed to estimate the state variables ( 1 ,  2 ,  3 ,  4 , and  5 ) based only on the measurement of the states ( 1 ,  4 ).
The model described by (3a)-(3f) can be given in the following two subsystems, denoted by Σ 1 and Σ 2 with where Obviously  1 =  1 and  2 =  4 are the measured output of the PV system connected to the grid.
Parameter estimator Equation (35b) The design strategy consists in synthesizing separately an observer for each one of subsystems (6a)-(6b).In the first step, a linear observer is designed for subsystem Σ 1 .In the second step, an adaptive observer is designed for subsystem Σ 2 .
3.1.1.Linear Observer.The form of system Σ 1 suggests the following Luenberger observer for the estimation of the unknown state variables ( 1 ,  2 ,  3 ): where ] are the observer gain.
From (6a) and ( 8) one obtains the estimation error dynamics: with The gain vector  is selected to make  0 =  1 −   1 a Hurwitz matrix, which will guarantee the asymptotic error convergence.
Introduce the following Lyapunov function candidate: where  is asymmetric positive definite matrix.The derivative of   along the trajectory of  1 is given by The  matrix is chosen as   0  +  0 = −, where  denotes the identity matrix, and this choice leads to 3.1.2.Adaptive Observer.The system Σ 2 is state-affine in the sense that all unknown states come in linearly.According to the [18] methodology the following adaptive observer is developed to estimate the state variables for subsystem Σ 2 : where   2 is the solution of the fallowing equation: The matrix   2 is ensured bounded positive definite provided the following persistent excitation condition holds: for some constants ( 1 ,  1 , ) > 0, where   0 () denotes the transition matrix for the subsystem Ẋ2 =  2 ( 0 ) 2 ,  2 =  2  2 , and Σ is a positive definite bounded matrix.The system can be seen as a linear time-varying system parameterized by initial conditions as soon as the function  0 is fixed.
To study the convergence of the proposed observer (13), defining the estimation error as its dynamics are given by Proof.To analyze the error system (17), the following Lyapunov function candidate is considered: Its time derivative is given by where  = min() and  are any positive constants.
Proposition 1.Under condition (15), the estimation error ( 17) is exponentially vanishing; that is, the estimate X2 converges exponentially to its true value  2 with a rate driven by  [19].

Output-Feedback Controller.
The PV output voltage  4 is regulated by controlling the switching device of the boost converter, while the DC-link voltage  6 and the grid current  1 are adjusted through the switching devices of the inverter.

Control of Boost Chopper (PV Voltage Controller).
The control objective is to enforce the voltage provided by the PV panel to track the desired signal in order to achieve maximum power point.This regulator consists of two loops: a loop for seeking of the nominal power point, in which we used the IncCond algorithm, and a loop for regulating the voltage x4 .Then, to reach these aims, we seek a control law of slidingmode type.This law takes the following form: where  0eq is nominal control, it keeps the system on the sliding surface, and  0 is discontinuous control, and it allows to reach the sliding surface.
To design a controller for subsystem (3d)-(3e), the error is defined as follows: Its dynamic are In order to stabilize subsystem (3d)-(3e), the slidingmode control technique is used.The sliding surface is defined as follows: where  is the parameter of the sliding surface.
The dynamics of the surface are given by The equivalent command  0eq is calculated by means of the method of "Utkin"; setting the sliding surface Ṡ 1 () = 0, the nominal control law can be defined as and, to elaborate the discontinuous control, consider the following Lyapunov function candidate: and its dynamics are given by Using ( 20) and ( 25), the dynamics of the Lyapunov function are written as: One seeks  0 of the type This choice guaranties the negativity of the dynamics of Lyapunov function: where  is positive parameter.The final control law is given by Proposition 2. Consider the closed-loop control system, consisting of system (6b) in closed loop with the control law (31) and the state adaptive observer error (17).Its dynamic behaviour, expressed in ( 1 ,  1 , X2 )-coordinates, is governed by the following equations: System ( 32) is globally asymptotically stable with respect to the Lyapunov functions ( 18) and (26).

Control of Single-Phase Inverter.
To guarantee high performance transmission of the power and good functioning of the system, the current and voltage grid should be in phase.
Hence, there is a necessity for regulator that enforces the estimate current x1 to track a given reference current x * 1 =  x1 / θ2 , where  is computed from the output of the outer loop voltage.The proposed controller uses a cascaded loop: an outer voltage loop and an inner current loop.The former compares the sensing DC bus voltage in the link capacitor with the given reference, whereas the latter uses slidingmode controller to regulate the grid current allowing the synchronization of the currents   with utility grid voltages V  .
(i) Network Voltage and Impedance Observer.It is readily checked that equations (3a) and (5a)-(5b) could be given the following compact form: with Then, (33) suggests the following state observer [19] for the estimation of the inaccessible state vector   and the unknown parameter vector : The notation x ∈ R 3 and θ ∈ R 2 , respectively, denote the estimate of the state variables   and the estimate of the unknown parameters .  and   are symmetric positive definite matrices with   (0) > 0 and   (0) > 0. Furthermore,   and   are positive constants that ensure the rapid convergence of the observers.The Λ dynamic is an auxiliary system that can be seen as filter.
The convergence properties of the adaptive observer are analyzed based on the following error system dynamics: with the errors   = x − Λ θ, θ =  − θ, and x =   − x , and, following the same idea as in [19], we indeed get One can choose the Lyapunov function (38), to analyze the convergence properties of the observers (35a)-(35e): Its time derivation leads to the following inequality: knowing that and this implies that with  min = min(  ,   ).
The stability results are summarized in the proof which can be found in [11].Therefore, the error system is exponentially stable with a rate driven by  min .
(ii) Power Factor Correction (Inner Current Loop).We seek a law control of sliding-mode type.For the PV system connected to the single-phase grid represented by the average model of (3a)-(3c) and (5a), the sliding-mode control which ensures that the error  2 ( 1 ) tends asymptotically to zero in finite time can be written as follows: To design this controller, one defines the sliding surface: where  = 3 is the relative degree of the system, and  2 ( x1 ) = x1 −  * 1 is the error between the signal and its reference  * 1 .The surface is given by Consider the first derivative of The command  1 appears in the first derivative of the sliding surface.Then the equivalent command To elaborate the discontinuous control, consider the following Lyapunov function candidate: Its dynamics are given by Differentiating  1 with respect to time, using (42) and ( 45), the dynamics of the Lyapunov function are written as Moreover, (48) shows that, to ensure the stability of the closed-loop system, a choice for  1 is of the form  The choice of (49) guarantees the negativity of the dynamics of Lyapunov function, so where  is positive parameter and sgn(⋅) is the sign function defined as follows: sgn The overall law is given by which is closer to the function sgn(⋅) but it is continuous especially at the point zero.The real positive constant  is selected sufficiently small to better approach the function sgn(⋅).Figure 3 shows that, for smaller values of , the function tanh(⋅) behaves close to a function sgn().Figure 4 depicts the errors between sgn() and tanh(/), which shows that the error decreases for small values of .

(iii) DC Bus Proportional-Integral Control (Outer Voltage Loop
).The objective of this loop is to generate the signal  so that the square of the DC bus voltage  = ( 6 ) 2 is regulated to a reference value  * = ( * 6 ) 2 .Considering the fact that  and its time derivatives should be available, a filtered control law of type PI is retained; namely, with where "s" denotes the Laplace variable and ( 3 ,   ,   ) are any positive real constants.

Simulation
In this section, the controller that has been designed in the above section using output nonlinear feedback technique will be tested.The simulation results have been obtained under normal conditions ( = 1000 w/m 2 ;  = 25 ∘ C) and other different conditions.The whole system (power system and controllers) is simulated with the numerical values of Tables 2 and 3.The parameters of the controllers are determined based on the desired specifications of the closed-loop system.In this paper, the controller parameters are designed so that the current control loop has much faster response than that of the DC-link voltage control loop.In this work the value  = 0.1 is used.

Tracking Performance in the Presence of Meteorological
Constant .show the results of simulation, which are selected to demonstrate the most significant aspects of the PV system grid-connected behaviour.Figure 5 shows that the observers do well.Also, Figures 5(a)-5(e) show that the estimated state provided by the observers converges rapidly due to its true values.Figure 6 shows the grid current spectrogram, where the THD value of this current is very low equal to 1.37%, which is still below the international standards.The performances of the grid adaptive observer are illustrated by Figures 7-10.Indeed in Figures 7 and 8 it is observed that the estimated parameters θ1 and θ2 , provided by the grid adaptive observer, converge rapidly to their true values, respectively, ( 1 +   )/( 1 +   ) and 1/( 1 +   ) after a short transient phase.Figure 11 illustrates the response of the PV voltage in standard climatic conditions.It can be clearly seen that, in steady state, the PV generator provides the maximal power, which is equal to 250 V. Figure 12 presents a view of the injected current in the grid with its reference; as it is seen in the figure that the current follows its reference with error nearly 0. Figure 13 shows the regulation of the DC-link voltage.As shown in the figure, the voltage is maintained at a constant level (650 V).Consequently, the real power extracted from the PV generator can be totally transferred to the grid.In the last figure, Figure 14 shows the current injected to the grid.As can be seen in this figure, the current and the grid voltage are in phase and sinusoidal.As a result, a unit power factor is achieved.
Different powers of the system are given by Figures 15-17.Figure 15 shows the power provided by the PV generators.According to Figures 16 and 17, the system injects into the grid active power  grid , where the reactive power  grid is kept around zero.

Tracking Performance in the Presence of Meteorological
Variation.The robustness of the nonlinear adaptive controller is checked under a time variation climate conditions.The purpose of this simulation is to test the closed-loop system under a change in the PV output that occurred by a sudden change in climate change.In the following simulations, the DC-link capacitor voltage is kept constant equal to 650 V.
The levels of the irradiance are illustrated in Figure 18, which shows that the irradiation has increased from 900 W/m 2 to 1000 W/m 2 at 0.3 s.
In the following simulations, the output-feedback controller performances are illustrated by Figure 19.The curves (a-c) show that the tracking quality of the proposed observers is quite satisfactory for all supposed unknown states.irradiance, the inner loop ensures a perfect asymptotic tracking of the current reference signal.Figure 19(c) shows that the DC bus voltage regulation is recovered after a short transient period following each change of the irradiation.Finally, Figure 19(d) shows that the correction of the power factor is preserved even during variations of irradiation.

The Controller's Ability to Compensate for Variation in
Network Impedance Parameters.The aim of this test is to check the grid current controller performance in the presence of variations in the network impedance.The elements of the impedances (  ) and (  ) are modified, according to the protocol presented in Figure 20, to generate variations at  1 and  2 .Note that all system parameters are those used for previous tests.The following figures present the simulation results in the case of variations undergone by the network impedance.Figures 21 and 22 show that unknown parameters and network current estimates converge rapidly to their actual trajectories.Finally, Figure 23 illustrates that the adaptive controller provides satisfactory power factor correction during the uncertainty interval.

Tracking Performance in the Presence of Grid Faults.
This test aims at evaluating the performance of the system against the change of the amplitude and the frequency of the grid.The amplitude and the reference of the grid are modified according to the protocol presented in Figure 24.
Figures 25 and 26 show that, in spite of the variation of the amplitude and the frequency of the grid, the regulator has the capacity to track the reference provided by the outer loop and to inject a current in phase with the grid voltage.
Note that, for all simulations, it is clear that the proposed controller reacts in a quick manner to reach the reference and to remove the steady-state error quickly to keep the stability of the system.

Conclusion
An output-feedback nonlinear control strategy for a singlephase grid-connected PV system is proposed in this paper.The system is described by 6th order nonlinear averaged  designed.Specifically, it is shown that all control objectives are achieved, including PFC requirement, extracting a maximum power from the PV array and DC-link voltage regulation without requiring a lot of current and voltage sensors.

Figure 1 :
Figure 1: Photovoltaic system tied to the single-phase grid.
Photovoltaic generator with BOOST converter

Figure 5 :
Figure 5: (a) Grid current and its estimate.(b) Capacitor filter voltage and its estimate.(c) Inverter output current and its estimate.(d) PV voltage and its estimate.(e) Current of inductance boost and its estimate.

Figures 9
and 10 show that the estimated signals of the grid

Table 1 :
State variables and unknown parameters.

Table 2 :
PV system and single-phase grid characteristics.