MPPT Algorithm for Photovoltaic Panel Based on Augmented Takagi-Sugeno Fuzzy Model

This paper deals with the Maximum Power Point Tracking (MPPT) for photovoltaic energy system. It includes photovoltaic array panel, DC/DC converter, and load. The operating point for photovoltaic energy system depends on climatic parameters and load. For each temperature and irradiation pair, there exists only one optimal operating point which corresponds to themaximum power transmitted to the load. The photovoltaic energy system is described by nonlinear equations. It is transformed into an augmented system which is described with a Takagi-Sugeno (T-S) fuzzy model. The proposed MPPT algorithm which permits transfering the maximum power from the panel to the load is based on Parallel Distributed Compensationmethod (PDC).The control parameters have been computed based on Linear Matrix Inequalities tools (LMI). The Lyapunov approach has been used to prove the stability of the system. Some reliable simulation results are provided to check the efficiency of the proposed algorithm.


Introduction
In the most recent years, photovoltaic (PV) energy has been the subject of several research projects. It is well known that the PV array power panel depends on climatic variables such as temperature and irradiation as shown in Figures 2 and 3. Actually, the operating point of the PV array panel depends on three parameters such as temperature, irradiation, and the load. In fact, the operating point results from the intersection of the I-V characteristic and the load characteristic as shown in Figure 4. In most cases, the value of load is constant and the climatic parameters vary in the day, so the load characteristic remains fixed and the characteristic of the panel varies according to climatic variables. Consequently, the operating point is variable and the load cannot extract maximum power from the panel. To overcome this disadvantage, a DC/DC converter is inserted between the panel and the load. In this way, the load value seen by the PV panel can be changed by varying the duty cycle. In this context, several studies have been developed.
Most of papers dealing with the MPPT control algorithms are based on perturb and observe (P&O) [1][2][3][4], Incremental Conductance [3,5], Mamdani type fuzzy logic controller (FLC) [4,6], and some different approaches as neural network controller (NNC) [7]. In this paper, the PV array panel has been modelled by fuzzy system approach. At every time, the desired state variables have been computed based on the measurement of temperature and irradiation. Also, the MPP tracker algorithm has been developed based on Parallel Distributed Compensation (PDC) method which was designed for fuzzy system. Most of papers which have used fuzzy system applied Mamdani method. Whereas in this work, the main contribution which deals with Maximum Power Point Tracking for photovoltaic panel, consists on developing a new algorithm based on augmented T-S fuzzy system. It is well known that the Mamdani fuzzy system includes three blocks which are fuzzification, fuzzy inference rules and defuzzification, whereas the T-S fuzzy system needs only two blocks such as fuzzification and fuzzy inference rules. The block number has been reduced. However, we decided to choose the T-S type fuzzy system. In this paper, the T-S fuzzy system has been used in modelling stage of PV system and in the control stage. The optimal duty cycle, which permits extracting the maximum power from the photovoltaic array panel, is computed based on PDC techniques. This paper is organised as follows. In Section 2, we describe the problem statement and we show the influence of climatic parameters such as temperature and irradiation on the electrical characteristic of the PV array panel. Then, we recall the model of photovoltaic panel. At the end of the second section, we describe the photovoltaic energy system by a state model. The third section is reserved to present the control strategy. In the first part, we recall the T-S fuzzy system; then we describe the photovoltaic energy system by an augmented T-S fuzzy model. In the second part of the third section, we describe the T-S fuzzy reference model and the MPPT algorithm which is based on T-S fuzzy system. Section four is devoted to the stability analysis of the closed loop system. The feedback gains have been computed by solving LMIs expressions. The simulation results of photovoltaic energy system showing performances of the proposed MPPT algorithm tracker are discussed in Section 5. Conclusions are drawn in the final section.

Problem Statement
The photovoltaic power depends on climatic parameters such as temperature and irradiation as shown in Figures 2 and 3. In fact, the photovoltaic power, which is transmitted to the load, is function of the impedance of the load and the climatic parameters as shown in Figure 4. However, to change the impedance seen by the panel, it is necessary to insert a DC/DC converter. The photovoltaic system consists of a photovoltaic array panel connected to a DC-DC converter which provides energy to the load, as shown in Figure 1. Figures 2 and 3 show the evolution of the generated power curves as a function of voltage, respectively, for a given constant irradiation and different values of temperature and then for a given constant temperature and different values of irradiation.
In conclusion, we can say that the PV array panel is nonlinear and time-variant system. From Figures 2 and 3, it is clear that the temperature affects essentially the voltage and the irradiation affects fundamentally the intensity of the PV array panel. Also, we can conclude that the output power generated by the PV array panel depends on the climatic parameters " and . " In fact, the power increases with an increase in solar radiation and decreases with an increase in temperature. For each given pair of parameters ( , ), there exists only one Maximum Power Point (MPP). The operating point is determined by the intersection of the panel current-voltage characteristic and the load current-voltage characteristic.   However, a specific algorithm tracker should be used to search the optimal operating point which permits to extract the maximum power from the PV array panel.
In the next part of this section, we recall the most popular model which is developed by Singer et al. [8]; then we describe the modeling of the overall photovoltaic system energy.
The electrical equivalent circuit of the PV cell is given by Figure 5.
It consists of a current generator which depends on irradiation ( ) and temperature ( ), in parallel with a diode, and connected to an internal parallel and series resistor, namely, respectively, sh and .
The PV cell model is described by the following equations: (1) The generated current by the photovoltaic panel varies with temperature and irradiation; its expression is given by the following equation: ph, is the rated current generated by the PV panel under standard condition of temperature and irradiation ( = 25 ∘ C and = 1000 w/m 2 ): where is a reverse saturation current: where oc is the open circuit voltage and sc is the short circuit current.
The overall photovoltaic system energy can be represented by the scheme illustrated in Figure 6.
The average dynamic model of the photovoltaic system given by Figure 6 can be expressed in continuous conduction by the following equations: In the continuous conduction, the average value of PV current is equal to the average value of current. It is very clear that the system can be described as the form oḟ( where ( ) = [ pv 2 ] is the state vector, ( , ) is the state matrix, is the input vector, and is the duty ratio. However, the state matrix is nonlinear.

Control Strategy
The control strategy that we propose is given by Figure 7.
The control strategy consists of three blocks: reference model, controller, and plant, which is defined by PV system (see Figure 1).

T-S Fuzzy Model.
Several studies have been proving that the Takagi-Sugeno fuzzy system can describe the behavior of continuous nonlinear system. However, we use in this work the T-S fuzzy system to describe the nonlinear energy conversion system. The fuzzy model is described by fuzzy rules where each rule represents input-output relations of linear local model. The th rule of the fuzzy model has the following form: IF 1 is 1 and 2 is 2 and, . . . , and is where { } are the fuzzy sets, ( ) is the state vector, ( ) is the input vector, is the state matrix, is the input matrix, 1 ( ), . . . , ( ) are the premise variables, ( ) ∈ , ∈ × , ∈ , ( ) ∈ , and is the number of fuzzy rules. The global fuzzy model of the system has the following form:̇( .  Figure 6: Photovoltaic system.

Photovoltaic energy conversion system
T-S fuzzy reference model T-S fuzzy controller (PDC) For each rule attributed a weight ( ( )) which depends on grade of membership function of premise variables ( ) in fuzzy sets : ( ( )) is the grade of membership of ( ) to the fuzzy set : The polytopic form of state equation iṡ

Augmented T-S Fuzzy Model of Photovoltaic System.
In the first stage of this subsection, we transform the average dynamic model of the photovoltaic system described by (5) into an augmented model. However, a new state variable must be added to the state vector. In other words an integrator is included previous to the real input and leṫ= . Therefore, becomes new state variable and is the new control input of the augmented system. Then, the system can be described as a nonlinear system in the form oḟ where ] is the state vector, ( , ) is the state matrix, and is the input matrix: ] .
It is clear that (12) is nonlinear. To obtain the T-S fuzzy model, we choose the following three premises variables: PV , PV , and 2 . However, eight local models have been obtained to describe the T-S fuzzy model. The state matrix of each local model has the following structure: where each of the variables = pv / pv , = 2 , and = must be replaced, respectively, by the appropriate value according to fuzzy rule base ( min or max ), ( min or max ), and ( min or max ).

Reference Model.
It is well known that the maximum power produced by panel, also the corresponding optimal Photovoltaic energy conversion system T-S fuzzy reference model MPPT algorithm based on T-S fuzzy controller tension MP and optimal current MP , depends on temperature and irradiation . However, these latest are used as premise variables, 1 = and 2 = , to compute the reference model based on T-S fuzzy system. The th rule of the fuzzy reference model has the following form: where { } are the fuzzy sets, ( ) = [ pv 2 ] is the state reference variable vector, ∈ × is the local reference state matrix, and { 1 ( ), . . . , ( )} are the premise variables. ( ) ∈ is the output vector, is the number of fuzzy rules. Then, the T-S fuzzy reference model is given by the following equation: with ∑ =1 ( ( )) = 1 and is the local reference state matrix:

MPPT Algorithm Based on T-S Fuzzy Model.
The crucial function of the MPPT algorithm is to search the coordinates of the optimal operating point. In this work, the MPPT algorithm provided at each time the appropriate duty ratio based on the desired and measured state variables. However, the MPPT algorithm which represents a controller is based on Parallel Distributed Compensation (PDC) method. The T-S fuzzy controller is designed from the T-S fuzzy model of photovoltaic energy conversion system. The T-S fuzzy controller is based on Parallel Distributed Compensation (PDC) technique, which is proposed by Wang et al. [9]. The particularity of the PDC controller is that it shares the same fuzzy sets in the premise parts with the fuzzy model. Based on the T-S fuzzy models, the PDC fuzzy controller is designed as follows. th controller rule: The global fuzzy controller is represented by

Stability Analysis
In the previous section, we discussed the main blocks of the control strategy and we have developed an MPPT algorithm ( Figure 8). However, the system of energy conversion transmits the maximum PV generated power to the load and it is crucial to ensure that the closed-loop system is stable.
Proof. The state tracking error is given by The following quadratic Lyapunov candidate function which is positive definite, has been used to verify the system stability and compute the feedback gains : The system is asymptotically stable if we prove thaṫ( ) < 0: We note in this analysis that the feedback gains and stability conditions will be transformed to an LMI problem. However, the inequalitẏ( ) < 0 is guaranteed when the following inequalities are satisfied: Using Schur complement [10], the inequalities given in (24) can be written as Since coupled elements, such as , have been enclosed in these inequalities, then we have BiLMIs inequalities. However, we must transform them to the LMIs using a congruence transformation by diag [ −1 ] to (29) and considering = −1 , = −1 , we obtain the following matrices in the LMI form:

Simulation Results
This section is reserved for presenting the main results. We use MATLAB to simulate the behavior of the energy conversion system. The main characteristics of the PV array panel are given by Table 1.
The resolution of the LMIs gives the following matrices and feedback gains, respectively, , , 1 , 2 , 3 , 4 , 5 , 6 , 7 , and 8 .  To demonstrate the performance of the proposed MPPT control approach, we apply a sudden variation of temperature or solar irradiation as shown in Figures 9 and 10. In Figure 9, we have applied a sudden change of temperature, although it is impossible to have a really dramatic change.
In this test, we have chosen four pairs of irradiation and temperature. We know that for each pair there exists only one optimal operating point which can be determined from the power-voltage characteristics of the PV array panel which is not always available for each pair ( , ). It is important to mention that it is not possible to know the appropriate coordinates of the ideal optimal operating point ( MPP , MPP ) for all pairs ( , ) as there are infinitely of pairs ( , ).
In Table 2   and irradiation and the computed values ( MPP , MPP ) by our algorithm. Figures 11, 12, 13, 14, 15, 16, 17, 18, and 19 show, respectively, the evolution of MPP voltage,error of MPP voltage, output   voltage of converter, error of output voltage of converter, panel current, error of panel current, delivered power, error of delivered power, and the duty cycle. In Figures 12, 13, 14, 15, 16, 17, and 18, we observe momentary peaks; they are due to sudden and significant change in temperature and irradiation. The changes in temperature and irradiation are not made as that way given in Figures 9 and 10, but we have used it to show the performance of the proposed algorithm. It is clear that at the steady state, the errors tend toward zero and the state variables reach the reference one. Also, it is visible that the computed coordinates, of optimal operating point, based on the proposed algorithm, are almost the same as the ideal optimal operating point. This analysis allows demonstrating the performance of the proposed algorithm.

Conclusion
In this paper, a new algorithm strategy based on the augmented Takagi-Sugeno type fuzzy system has been proposed for the MPPT of a PV energy system. All the PV system has been modeled by T-S fuzzy system. Based on the measurement of temperature and irradiation, we deduce the coordinates of the desired optimal operating point which corresponds to the maximum power. The MPPT algorithm is based on an augmented T-S fuzzy model and PDC method. The controller parameters have been computed based on the LMI tools. The stability of system has been proved based on Lyapunov approach. The simulation results show that the proposed algorithm tracks quickly the optimal operating point despite sudden variations of temperature and irradiation.