Very Fast and Accurate Procedure for the Characterization of Photovoltaic Panels from Datasheet Information

In recent years several numerical methods have been proposed to identify the five-parameter model of photovoltaic panels from manufacturer datasheets also by introducing simplification or approximation techniques. In this paper we present a fast and accurate procedure for obtaining the parameters of the five-parameter model by starting from its reduced form. The procedure allows characterizing, in few seconds, thousands of photovoltaic panels present on the standard databases. It introduces and takes advantage of further important mathematical considerations without any model simplifications or data approximations. In particular the five parameters are divided in two groups, independent and dependent parameters, in order to reduce the dimensions of the search space. The partitioning of the parameters provides a strong advantage in terms of convergence, computational costs, and execution time of the present approach. Validations on thousands of photovoltaic panels are presented that show how it is possible to make easy and efficient the extraction process of the five parameters, without taking care of choosing a specific solver algorithm but simply by using any deterministic optimization/minimization technique.


Introduction
The one-diode model for the photovoltaic (PV) panel characterization has been widely used within both specific software toolboxes for the estimation and the prediction of the electrical power produced by PV plants [1][2][3][4] and algorithms for the Maximum Power Point Tracking [5][6][7] or irradiance measurements [8,9].Indeed, it guarantees a good tradeoff between accuracy and complexity for its setup [10,11].On the other hand, the extraction of the five-parameter model at standard reference conditions (SRC) (i.e., an inverse problem) has been widely faced in the literature.Although two approaches are generally the most adopted (the one that uses the datasheet information and the other one that exploits the experimental data on I-V curves), the use of only data provided by manufacturer on datasheet appears more interesting because it does not require a specific experimental study on PV module.Nevertheless, the approaches proposed the in literature differ between them and often it is difficult to understand what is the best one to be used.Indeed, on one hand several works proposed different equations/approaches for the extraction of the five parameters; on the other hand, almost any kinds of optimization techniques have been presented to solve the inverse problem of the extraction of the five parameters.This is essentially due to the nature of the involved equations which are transcendental and hard to manage.Just to give some references within the wide literature regarding this issue, hereafter some of the more recent works are briefly reported.Regarding the techniques for finding the inverse problem solutions, in [12] an improved differential evolution algorithm is presented for the extraction of five parameters from both synthetic data and experimental  data, in [13] penalty differential evolution is used in a similar way, in [14,15] pattern search and Bacterial Foraging Algorithm are used, respectively, and so on (see the reference within these works for further journal articles).Regarding the alternative analytical approaches, in [16] an explicit - model of a solar cell which uses Padé approximation is presented, that is, the exponential function is approximated by means of a rational function; in [17] the Taylor series is instead used; in [18][19][20] the - relations are more explicitly written by means of the  Lambert function [21], and then the extraction of the five parameters is performed by numerical techniques.Although the aim of these works is to effectively solve the problem, they suffer from unsuitable mathematical approximations (which lead to errors in the results) or complicated implementations and high computational costs.As a consequence these approaches are not so easy to be applied.
In this paper we present a fast and accurate procedure for obtaining the parameters of the five-parameter model by starting from its reduced form [22] which allows the identification of thousands of PV panels available on the standard databases (such as Californian Energy Commission database [23]).Indeed, by using suitable initial guesses it is possible to fully characterize thousands of PV panel in few seconds without any model simplifications or data approximations, simply by introducing further important mathematical considerations about the five-parameter model.The paper is structured as follows: in Section 2 the traditional onediode model and the problem related to its characterization are presented; in Section 3 the reduced form of the fiveparameter model is described; the validation results obtained on thousands of PV panels are shown in Section 4; finally, Section 5 is for the conclusions.

The One-Diode Model and Its Reduced Form
The equivalent circuit of the one-diode model is shown in Figure 1.The relation between current  and voltage  for a PV array/panel of arbitrary dimension (  parallel connected strings of   series-connected PV cells) at the equivalent port is [10] where  irr is the irradiance current (photocurrent),  0 is the cell reverse saturation current (diode saturation current),  is the electron charge ( = 1.602 × 10  [11]; other slightly different relations are presented in several works, such as [24,25], but their use does not affect the effectiveness and validity of the presented procedure): In ( 5)   = 1.17 − 4.73 × 10 −4 × ( 2 /( + 636)) is the bandgap energy for silicon in eV.Thus, there are five unknown parameters at SRC,  ref ,  ,ref ,  irr,ref ,  0,ref , and  SH,ref to be found within (2)- (6).Then, by starting from their values and by using the above relations, it is possible to write the - curves for every temperature and irradiance values.In order to determine these five unknowns we need five independent equations based on datasheet information.Usually the PV panel manufacturers provide several information on datasheet at standard reference conditions (SRC), that is, for the irradiance  ref and the temperature  ref : the values of the short-circuit current ( SC,ref ) and the open-circuit voltage ( OC,ref ), the current and voltage values at the maximum power point ( mp,ref and  mp,ref ).In addition, the datasheets report the temperature coefficients (or percentage) of both the short-circuit current (  or  %  ) and the open-circuit voltage (  or  %  ).On the basis of the three characteristic points, open-circuit, short-circuit, and maximum power points at SRC, it is possible to write the first four equations of the five-parameter model [10,11,15,22,26]: indeed the first equation arises by writing (1) for the open-circuit (OC) condition; the second equation arises by using (1) for the short-circuit (SC) condition; the third equation arises by exploiting the current and voltage values at the maximum power point (MPP) condition.The fourth equation is written by imposing the slope of the - curve (power versus voltage) over the MPP equal to zero, ( ⋅ )/ = 0, that can be also expressed in terms of  mp,ref / mp,ref ratio.The last fifth equation used to complete the five-parameter model is written by exploiting the temperature dependence of ( 1) at the open-circuit condition and irradiance  =  ref by using the previously stated temperature coefficients (  and   ) [10][11][12].Before showing the five equations, it is useful to briefly recall the constants specified in (7) adopted in order to simplify the writing of the mathematical expressions.Furthermore, the temperature-dependent factor is used in ( 5) and the shunt conductance  SH,ref =  −1 SH,ref is adopted as unknown in (6) instead of  SH,ref : Thus, the five equations are the following: In ( 12), a value of  =  ref ±10 K is used, even if variations of temperature Δ belonging to the range [−10, +10] with respect to  ref return very similar solutions [10,11].
The five-parameter model is thus defined by a system of five equations, ( 8)- (12), with the five unknowns (parameters),  ref ,  ,ref ,  irr,ref ,  0,ref , and  SH,ref .Due to the presence of transcendental equations this problem is not so simple to manage and it can be only solved by means of numerical methods.Since it is practically an inverse problem, many minimization algorithms can be used and almost any kinds of computing techniques have been tested in the literature: for example, in [27] the comparison between several techniques to extract the five parameters is presented and compared by using the criteria of applicability, convergence, stability, calculation speed, and error on various types of  data.In addition, due to its nonlinear nature, the system returns solutions that are very sensitive to the choice of the initial guesses [22,26,27].As the following section shows, this problem can be easily overwhelmed by using a reduced form of the model employing only a set of two equations in two unknowns.For the reader's convenience, the list of the technical parameters used in this work is reported at the end of the paper.

Reduction to a Two-Parameter Model
In [22] it has been proven that the five-parameter model can be reduced to a two-parameter model improving the efficiency of the algorithm finding the solution.By using this reduced form of the system instead of the original one, it is also possible to demonstrate: (i) the uniqueness of the solution for the problem; (ii) the existence of a unique solution without physical meaning for some PV panels; (iii) the matter of the optimal choice of the initial guesses that make easy and effective the solution of the inverse problem.
As first thing, let us show the way to reduce the five-parameter model to a two-parameter model.This is obtained by simple algebraic manipulations of three of the five equations ( 8)- (11).Indeed from the first equation (8) By substituting ( 13) in (10), it is possible to write that can be written also as On the other hand, from (11) we can also obtain  0,ref Furthermore, by posing the expression (15) equal to the (16), we can write from which it is possible to obtain  SH,ref as a function of  ref and  ,ref , Now, substituting (18) in (15) or ( 16) also allows expressing  0,ref as a function of  ref and  ,ref ,  International Journal of Photoenergy Finally, (19) and ( 18) are utilized together with (13) so that  irr,ref can be written as a function of  ref and  ,ref , In ( 18), (19), and (20)

Solution of the Reduced form of the Five-Parameter Model.
Although the two equations ( 21) of the reduced form of the five-parameter model are transcendental equations, they are quite affordable that simple and fast numerical methods can be utilized to find the solutions instead of more complex and expensive algorithms in terms of computational costs [28][29][30][31].On the other hand, the effectiveness of the reduced form with respect to the original system based on five equations is evident since it returns the same unique solution also by using different numerical methods.Moreover, ( 18)-( 21) allows making several important considerations about the solutions of the five-parameter model.
(i) Since the values of  irr,ref ,  0,ref , and  SH,ref must be positive in order to obtain a physical meaning for these three parameters, it is possible to find the conditions for the range of the independent unknowns  ref and  ,ref from ( 18), (19), and (20).In particular, the maximum admissible value for  ,ref is a function of  ref , according to the following relations [22]: with  max ,ref as a function of  ref whereas the Lambert  function [8] in (23) has been used: this special mathematical function allows obtaining a closed form representation for the - curves and it is often successfully used for the analysis of PV modules [14,15].
(ii) Thus it is possible to individuate the feasible domain for the two remaining independent unknowns by assuming (iii) It is also possible to graph the 2D functional used for solving the system (21): where The presence of only one minimum (i.e., one global minimum) makes the problem of finding the solution of the system (21) a convex problem allowing the use of simple initial guesses without choosing specific optimization algorithms.Indeed, as also discussed in [26] where an empirical approach is adopted, the choice of the initial guesses is one of the more critical aspects regarding the identification of the fiveparameter model [22].
(iv) As a consequence, it is also possible to state that the system (21), that is, the reduced form of the original five-equation system, has a unique solution that corresponds to the one with physical meaning.
(v) By observing the graph of the functional ( ref ,  S,ref ) it is also possible to verify that some PV panels have the minimum (i.e., the solution of the problem) lying outside the feasible domain.This means that the solution still exists but it is not physical (i.e., at least one among the five parameters is negative).5 clearly prove the uniqueness of the solutions of the four panels.Furthermore, the quasi-monotonic behaviours of the functionals allow for an easy search of the solution (convex optimization).The initial guesses chosen for the search procedure of the solutions were the following:  with  guess = 1.0 for multi-Si and mono-Si PV panels and  guess = 2.0 for thin film PV panels.It can be noted by observing Figure 5 that the PV panel BP Q Series 230 W does not provide for physical solutions of the five parameters model (i.e., the solution exists but it lies outside the feasible domain and then at least one among the dependent parameters is lower than zero).It is worth noting that all the PV panels of BP Q series cannot be modelled by using the five-parameter model.The issue about the existence or not of the solution International Journal of Photoenergy  is really complex and nothing can be said a priori by simply observing the datasheets of the PV panels.

Tests on California Energy Commission Database
In order to prove the effectiveness of the proposed procedure, aimed at the identification of the five-parameter model simply by starting from PV datasheets, in this section a statistical validation is presented.In particular, the tests have involved around 11000 PV panels belonging to the California Energy Commission (CEC) database [23] (updated monthly).Table 1 shows the number of tested PV panels (# PV panels) from CEC database grouped by type of technology.With the aim to demonstrate both the robustness and the fastness of the proposed approach, several initial guesses have been chosen and three different numerical algorithms/functions have been used in Matlab: fsolve (suitable for systems with nonlinear equations), fminsearch (generic unconstrained minimization function), and lsqnonlin (function aimed to solve least squares problems).The simulations showed very accurate results (i.e., with functional values less than 1 − 20) which are independent of the adopted algorithm and the unique solutions have been found at first launch for almost all panels.The execution time for the extraction of the five parameters of all the 11764 PV shown in Table 1 performed on an Intel i5 core 2.5 GHz based notebook with 4 GB of RAM was around 90 seconds for the most efficient algorithm (fsolve with trust-region-dogleg) and 400 seconds for the slowest algorithm (fsolve with trust-region-reflective).This means that, also for the worst cases, the herein proposed extraction procedure of the five parameters spent less than 30 msec for each panel.The comparisons of the performance achievable by using various algorithms and initial guesses are reported in Tables 2, 3, and 4: in particular the results are expressed in terms of average number (mean # steps) and standard deviation (std # steps) of iterations steps, average number (mean # FEs) and standard deviation (std # FEs) of function evaluations (FEs) (in Table 2 the initial guesses are the ones proposed in [22]).
It is worth noting that the proposed initial guess (26) allows obtaining effective results also by using one of the most generic solvers for minimization problems, fminsearch, which employs the Nelder-Mead simplex method discussed in Lagarias et al. [32].By using instead more effective Matlab functions (such as fsolve or lsqnonlin) and algorithms (such as Levenberg-Marquardt [33] or trust-region-dogleg [34] algorithms) the number of iterations and FEs becomes extremely low (around 7 for the average number of iterations and 45 for the one of FEs).Consequently the computational costs of the proposed procedure is quite negligible, as the various soft computing based approaches, like the ones in [12][13][14][15], typically require thousands of FEs.In addition, it is worth noting that the obtained results have physical meaning for more than 97% of the total number of panels (11488 on 11764 PV panels).The unphysical solutions could be due to the impossibility of identifying the five-parameter model as was for the previous case of BP Q series PV panels.Nevertheless, since we do not have any information about how the datasheets are loaded into the CEC database, we assumed they were correct and no check was made about the exactness of data.Thus, some unphysical solutions could be also due to this last matter and caused by the presence of some errors within the CEC database.Finally, with the aim to show the importance of adopting an accurate 5-parameter model rather than approximated ones (such as for example the 3parameter and 4-parameter models [35]) the - and - curves at SRC for BP 3 235T module have been considered as last test.The 3-parameter and 4-parameter models seem to be very similar to the 5-parameter one.Nevertheless, with the aim to simplify the characterization problem,  ,ref = 0 and  SH,ref = inf conditions are used for the 3-parameter model and  SH,ref = inf condition is used for the 4-parameter model.The - and - curves for the three models are shown in Figures 6 and 8, whereas the close-up around the maximum power point (MPP) is shown in Figures 7 and 9.It is evident, by observing the - curves, that (1) the 5parameter and 4-parameter models are both accurate in the evaluation of MPP, whereas the 3-parameter model is not; (2) the 4-parameter model overestimates the power around the MPP, causing possible errors in the prediction of electric power produced by a PV plant.

Conclusions
In this paper a fast and accurate procedure has been presented for the characterization of thousands of photovoltaic modules in few seconds, by starting from the manufacturer datasheets.The proposed procedure utilizes the five-parameter model and takes advantage from its reduced form [22] in order to decrease the dimensions of the search space.Indeed, the reduced system provides a strong advantage in terms of convergence, computational costs, and execution time of the present approach (less than 30 msec for each panel was spent on a simple Intel i5 core 2.5 GHz based notebook).In particular, it allows (1) choosing suitable initial guesses within a well-defined feasible domain; (2) using very simple and standard numerical algorithms for finding parameters; (3) proving the existence, or not, of the unique physical solution of the five-parameter model for each PV panel; (4) proving the existence of only unphysical solutions for the cases in which the five-parameter model cannot be identified.The results of the tests performed on around 11.000 photovoltaic modules belonging to the CEC database demonstrated both the fastness and the effectiveness of the proposed method.

Figure 2 :
Figure 2: Graph of the functional (25) and the feasible domain for the BP 3 235 T PV panel (multi-Si technology).A unique physical solution exists.

Figure 3 :
Figure 3: Graph of the functional (25) and feasible domain for the Sharp NT175 panel (mono-Si technology).A unique physical solution exists.

Figure 4 :
Figure 4: Graph of the functional (25) and feasible domain for the PV panel Xunlight XR12-88 panel (thin film technology).A unique physical solution exists.

Figure 5 :
Figure 5: Graph of the functional (25) and feasible domain for the PV panel BP Q 230 panel (multi-Si technology).A unique unphysical solution exists (it lies outside the feasible domain); that is, this panel cannot be modelled by using the five-parameter model.

𝑓 1 (
ref ,  ,ref ) = 0 and  2 ( ref ,  ,ref ) = 0 represent the first and the second equation of the system (21), respectively.The graph of the functional ( ref ,  ,ref ) is smooth and free from local minima (some examples are shown in the next section).
Graphs.In order to prove the above considerations, in Figures 2, 3, 4, and 5 the results of four different panel modules are shown: a mono-Si PV panel (Sharp NT-175UC1), a multi-Si PV panel (BP 3235 T), a thin film PV panel (Xunlight XR12-88), and another multi-Si PV panel (BP Q Series 230 W).Each figure shows the graph of the functional ( ref ,  S,ref ) together with its feasible domain (i.e., the set of points of the two independent parameters  ref and  ,ref for which the dependent parameters  0,ref ,  irr,ref and  SH,ref , have physical meaning) and the solution (minimum of the functional) of the reduced system (21).Figures 2-

Figure 7 :
Figure 7: Close-up around maximum power point of Current-Voltage curves at SRC of BP 3 235T, traced by using 3/4/5 parameter models.

Figure 9 :
Figure 9: Close-up around maximum power point of Power-Voltage curves at SRC of BP 3 235T, traced by using 3/4/5 parameter models.
SH represent the cell series and shunt resistance, respectively.In (1), the governing variables ,   ,  irr ,  0 , and  SH can be assumed as dependent or not on the irradiance () and the temperature () but in any case they are in function of certain reference (ref) parameters at SRC ( ref = 1000 W/m 2 ,  ref = 25 ∘ C) (hereafter we present and utilize the relations proposed byDe Soto et al. in −19C),  is the cell ideality factor,  is the Boltzmann constant ( = 1.3806503 × 10 −23 J/K),  is the cell temperature, and   and we have written  SH,ref ,  0,ref , and  irr,ref as functions of  ref and  ,ref , respectively.This means that now there are only two independent unknowns, called  ref and  ,ref , to be found by using the two equations coming out from the other two conditions not yet utilized in the previous steps.They are (9) and (12) obtained from the evaluation of (1) at short-circuit ( = 0,  =  SC,ref ) condition and for   ≈ ( OC −  OC,ref )/( −  ref ) condition, respectively.Thus, the reduced form of the original five-equation system is the following:

Table 1 :
Number of tested PV panels from CEC database grouped by type of technology.

Table 2 :
Results obtained with Matlab fsolve for

Table 3 :
Results obtained with Matlab fsolve for

Table 4 :
Results obtained with different Matlab functions with ∘ C at SRC  ref : 1000 W/m 2 at SRC  ref : I d e a l i t y f a c t o r a t S R C  ,ref : Series resistance at SRC  irr,ref : P h o t o c u r r e n ta tS R C  0,ref : Reverse saturation current at SRC  SH,ref =  −1 ,ref : Shunt resistance at SRC  OC : O p e nc i r c u i tv o l t a g e  SC : Short circuit current  mp : M a x i m u mp o w e rv o l t a g e  mp : Maximum power current  OC,ref : O p e nc i r c u i tv o l t a g ea tS R C  SC,ref : Short circuit current at SRC  mp,ref : Maximum power voltage at SRC  mp,ref : Maximum power current at SRC   : Temperature coeff.for  SC  %  : Percentage temperature coeff.for  SC   : Temperature coeff.for  OC  %  : Percentage temperature coeff.for  OC  1 :  ref /  2 :  OC,ref /   1  3 :  mp,ref /   1  4 :  mp,ref /   1  5 :  SC,ref /   1 .