Harmony Search and NatureInspired Algorithms for Engineering Optimization
View this Special IssueResearch Article  Open Access
Parameter Estimation of Photovoltaic Models via Cuckoo Search
Abstract
Since conventional methods are incapable of estimating the parameters of Photovoltaic (PV) models with high accuracy, bioinspired algorithms have attracted significant attention in the last decade. Cuckoo Search (CS) is invented based on the inspiration of brood parasitic behavior of some cuckoo species in combination with the Lévy flight behavior. In this paper, a CSbased parameter estimation method is proposed to extract the parameters of singlediode models for commercial PV generators. Simulation results and experimental data show that the CS algorithm is capable of obtaining all the parameters with extremely high accuracy, depicted by a low RootMeanSquaredError (RMSE) value. The proposed method outperforms other algorithms applied in this study.
1. Introduction
Photovoltaic (PV) cells, normally assembled into modules or arrays on mounting systems, are capable of producing electrons when photons strike its surface. Taking the advantages of many promising features like renewability, less pollution, and ease of installation, PV generators are envisaged to be an important energy source for the future.
Due to the high initial cost of a PVsupplied system, predictive performance tools are widely used by engineers to optimize the system performance [1, 2]. PV manufacturers normally provide limited tabular data measured under the Standard Test Conditions (STCs), which correspond to a cell temperature of 25°C and an irradiance of 1000 W/m^{2} at 1.5 air mass spectral distributions. As reported in [3], PV generators always operate under environments far from the STCs. Due to this reason, the data available in the datasheet usually fail to fulfill the engineering requirements.
PV model, with the ability to predict characteristics of PV generators under an operating environment other than the STCs, is a predictive performance tool that allows consumers to maximize the cost effectiveness of the system before installation [2]. They are generally analytical equations based on a physical description that formulate PV generated current with the most crucial technical characteristics and the environmental variables, such as the operating voltage , the ambient temperature , and the irradiance . Over the years, significant research efforts have been contributing to the development of the behavioral models [4–8]. Among numerous modeling approaches, the SingleDiode Model (SDM) is the most widely utilized PV model in the literature. A general SDM includes five parameters, namely, photocurrent (), saturation current (), diode ideality constant (), series resistance (), and shunt resistance (). In order to adapt PV model behavior to different operating conditions, de Blas et al. [9] suggested to apply the procedure described in the International Standard IEC 891 that relates current and voltage of the PV characteristics at given values of and with the corresponding values at different operating environments. The Improved Single Diode Model (ISDM) presented by De Soto et al. [5] includes the dependence of the PV parameters on operating conditions. The normal parameters at the STCs are necessary to be determined in this model. Both SDM and ISDM are adopted in this study of parameter estimation.
Analytical methods [5, 10–12] are common approaches in estimating the parameters by mathematical equations. Although having the merit of simplicity, it is hard to further reduce the errors of the estimated values. Furthermore, analytical methods utilize the curve features or semiconductor parameters that are unavailable in the datasheet. This problem often reduces its feasibility. Recently, PV parameter estimation is deemed as a multidimensional optimization problem. Several computational intelligence methods, such as Genetic Algorithms (GA) [13], Chaos Particle Swarm Optimization (CPSO) [14], Firefly [15], and Pattern Search (PS) [16], were proposed in the literature. These algorithms usually extract relevant parameters by minimizing the Root Mean Square Error (RMSE) as the objective function in the optimization process. Askarzadeh and Rezazadeh [17] reported that the optimization methods produce better results than analytical methods.
Cuckoo Search (CS) is a natureinspired optimization algorithm based on the fascinating breeding behavior such as brood parasitism of certain species of cuckoos. In [18, 19], Yang and Deb reported that the CS algorithm outperforms Particle Swarm Optimization (PSO) and GA algorithms for various standard test functions. In this paper, a CSbased parameter estimation method for the SDM and ISDM is presented. Simulation and experimental results show superior accuracy and feasibility of the proposed parameter estimation method.
The rest of the paper is organized as follows. Section 2 explains both PV models (SDM and ISDM) used in this work. The objective function formulation is given in Section 3. This is followed by results and discussions in Section 4. The results comparison is also available here, and finally the conclusions are derived in Section 5.
2. PV Modeling
2.1. Single Diode PV Model (SDM)
PV cells are made of a variety of semiconductor materials using different manufacturing processes. The working principle of PV cells is essentially on the basis of the PV effect, which refers to the generation of a potential difference at the junction in response to visible or other radiation. When a PV cell is exposed to light, the semiconductor materials absorb photons, and accordingly charged carriers are generated. Potential difference and current in the external circuit lead to the separation of carriers in the internal electric field created by the junction and collection at the electrodes. The photogenerated charge carriers can be subsequently captured in the form of an electric current, that is, electricity . Eliminated the PV effect, a PV cell behaves like a conventional diode that does not depend on any light parameters. The Shockley diode equation is generally used to describe the current flowing through the diode ():
In (1), is the normal diode current, and represents the electrical potential difference between the two ends of the diode. The ideality factor is assumed to be independent of the environment variables and . denotes the thermal voltage of the PV, and its value can be written as a function of : where and represent the Boltzmann constant (1.380650 × 10^{−23} J/K) and the electron charge (1.602176 × 10^{−19} C), respectively.
SDM assumes that the superposition principle holds; that is, the total characteristic is the sum of the dark and illuminated characteristics [3–5]. As expressed in (3) below, the terminal current is therefore equal to subtracting the current diverting through the diode and . The equivalent circuit of the SDM is shown in Figure 1.
PV module is a particular case of the PV cells connected in series. If the number of the connected cells is up to , is scaled to times. The model equation is then rewritten as
In this sense, , , , and are the corresponding parameters of a PV module.
2.2. Improved Single Diode Model (ISDM)
The traditional SDM ignores the operating conditions effect on these parameters. However, some studies have shown that the parameters, such as and , vary under different environmental conditions. These are due to changes of temperature and irradiance . Aiming to evaluate the PV behavior at environmental conditions other than the normal values and , the relations between the operating parameters and the normal parameters are studied by numerous researchers.
In [4], the value of lightgenerated is reported to be linearly dependent on the solar irradiation under the influence of temperature:
where is the lightgenerated current at the STCs. , the shortcircuit current coefficient, is one of the ISDM parameters. The difference between the standard test temperature and is denoted by .
Based on the diode theory, Messenger and Ventre [20] presented an approximate linear expression for the diode saturation current , which can be expressed as where is the material band gap. Usually, is set at a reasonable level depending on the semiconductor materials (1.12 eV for crystalline silicon and 1.75 eV for amorphous silicon) in simulation and design tools [21]. De Soto et al. [5] present an estimation method for in a wide temperature range:
where is a normal value at the STCs ( eV for silicon cells and eV for the triple junction amorphous cells).
In [3], Lo Brano et al. study how the series and shunt resistances are affected by the solar irradiance. On the basis of the experimental data, the values of and are observed varying in inverse linear modes with :
where the values of the resistances and are evaluated under the STCs.
By using the aforementioned relations, the ISDM described in [5] is able to analytically describe the characteristic of a PV generator for each generic condition of operative temperature and solar irradiance.
3. Parameter Estimation
3.1. Formulation of Parameter Estimation Problem
PV parameter estimation is a process that minimizes the difference between the measured data and the calculated current by adjusting the normal PV parameters. When the number of experimental data is up to , the objective function can be formulated by RMSE as where for SDM and for ISDM. is the homogeneous form of (4) which expresses the characteristics of the SDM:
For the case of ISDM, and satisfy the relational expressions discussed in the previous subsection, namely, (5)–(8).
3.2. Cuckoo Search
The CS algorithm [18, 19], proposed by Yang and Deb, is a natureinspired stochastic global search algorithm that follows three idealized behavior rules.(i)A cuckoo lays an egg and dumps it randomly into other bird species’ nests.(ii)The best nests with high quality eggs will be carried forward to the next generation.(iii)There are a fixed number of available host nests. If a host bird discovers that the eggs are not its own, it will either throw these alien eggs away, or it may abandon the nest and build a brand new nest at a nearby location.
Based on the three rules, the basic steps of CS can be summarized by the pseudocode shown in Pseudocode 1. In the CS algorithm, a pattern corresponds to a nest while each individual attribute of the pattern corresponds to an egg laid by the cuckoo. On the basis of randomwalk algorithms, the general system equation of the CS algorithm is given in: where and denote the generation number and the pattern number , respectively. The product means entrywise multiplications. Here is the step size scaling factor which should be related to the scales of the problem of interest [19]. The th attributes of the th pattern is initiated by using (12),

where and are the upper and lower bounds of the th attributes, respectively. In each computation step, the CS algorithm checks whether the value of an attribute exceeds the allowed search range. If this happens, the value of the related attribute will be updated with the corresponding boundary value.
Before the searching process, the CS algorithm detects the most successful pattern as pattern. Among the existing algorithms exist for generating Lévy flights in the literature, Yang and Deb [18, 19] reported that Mantegna’s algorithm [22] works well in most of the optimization problems. Accordingly the evolution phase of the pattern initialized with the detection step of , which is given by (13) [23]:
where is 1.5 in the standard software implementation of the CS algorithm [30]. denotes the gamma function.
After initialization, the evolution phase of the pattern starts by defining the donor vector , where . The required step size of the th attributes can be calculated by the following equation: where and . The function generates a uniform integer between [25]. The donor pattern is then randomly adjusted by
The CS algorithm will evaluate the fitness of the random pattern. If a better solution is caught, the pattern will be updated. The unfeasible patterns are revised by the crossover operator given in (16) as follows: where is the mutation probability value in the standard software implementation [30]). The final step of a generation is to check if the revised infeasible patterns deliver a better solution.
4. Results and Discussions
With the aim of providing a thorough evaluation of the CS algorithm in estimating the PV parameters, both SDM and ISDM are considered in this paper. Two case studies are designed to estimate the CS algorithm in model parameters estimation: (i)a commercial 57 mm diameter solar cell (R.T.C. France [26]) operating at the standard irradiance level; (ii)a PV module (KC200GT Multicrystal Photovoltaic Module) operating under varied environment conditions.
During the parameter extraction process, the objective function is minimized with respect to the parameters range. In theory, the value of is slightly larger than that of . is in a loose range from 1 eV to 2 eV. is around the value provided by the datasheet (normally less than 0.02%/°C). The is usually less than 50 μA. As stated in [27], the ideality factor ranges between 1 and 2. PV modules produced by most manufacturers have less than 0.5 Ω and between 5 and 170 Ω [8, 28]. As for PV cell, the ranges of and can be scaled by simply dividing [29].
Statistical analysis is performed to evaluate the quality of the fitted models to the experimental data. Besides RMSE, other two fundamental measures, namely, Individual Absolute Error (IAE) and the Mean Absolute Error (MAE), are applied to evaluate in this paper. Equations (17) and (18) preset the IAE and MAE, respectively:
The optimization algorithms applied in this paper are programmed in MATLAB. Similar simulation conditions, including population size, maximum generation number, and search ranges, are set to ensure a fair evaluation (population size = 25; maximum generation number = 5000).
4.1. Case Study 1: Parameter Estimation for a PV Cell at the Certain Irradiance Level
Table 1 lists the model parameters of the R.T.C France PV cell at 33°C, which are extracted from the experimental data in [26]. The parameters obtained from the CS algorithm are compared with three different parameter estimation approaches: CPSO [14], GA [13], and PS [16]. From the RMSEs of these methods, which are listed in the last row of Table 1, the CS algorithm [30] outperforms the other three optimization methods. CS obtained slightly lower RMSE, recording 0.0010 in numerical value.
During the parameter estimation process for the SDM, the values of the objective function in different optimization algorithms are shown in Figure 2. The function “ga” in MATLAB [31], whose crossover rate = 0.8 and mutation rate , is utilized for the convergence process test. As for PSO implementation [24], the algorithm parameters are set as learning factors , inertia factors , , and velocity clamping factor . In Figure 2, no further improvement by GA is observed after 500 iterations. On the contrary, the CS algorithm showed continuous improvement until the maximum generation. The CS algorithm, whose convergence speed is slightly faster than PSO, shows the best accuracy result in the minimization task after 5000 iterations.
Table 2 lists the parameters of the ISDM obtained by the CS algorithm. In order to evaluate the accuracy of the CSbased estimation, these parameters are substituted into the ISDM. Since the demonstrates nonlinear characteristics, the PV terminal current is solved by the NewtonRaphson method [32] in this paper. In Table 3, the calculated results are compared with the experimental data to observe the agreement between them. The notations IAE_{SDM} and IAE_{ISDM} denote the IAE for SDM and ISDM, respectively. Although the RMSE of the ISDM is less than that of CPSO, GA, and PS, it is similar to the RMSE of the conventional SDM under a certain environmental condition.


4.2. Case Study 2: Parameter Estimation for a PV Module under Different Environment Conditions
In this section, the validity of the CS algorithm is evaluated using KC200GT PV module operating under different environment conditions. The estimated parameters, both in the SDM and ISDM, are shown in Table 4. As illustrated in Section 1, the main application of the parameter extraction is to predict the characteristics for design purpose. It is worth pointing out that the SDM parameters can only be extracted by the experimental data measured under a certain test condition. Significant errors may occur as the experimental data are measured under varying operating conditions. In the commercial simulation tool like PSIM [21], the PV parameters of the SDM are firstly estimated at the STCs, then the equations (given in the appendix) are applied to calculate the electrical characteristics of different operating conditions. The ISDMbased parameter estimation, however, can be performed by the data measured under any conditions.
(a) SDM parameters (extracted by the CS algorithm)  
 
(b) ISDM parameters (extracted by the CS algorithm)  

Figure 3 shows the curves generated using the parameters obtained by the CS algorithm. The simulated results are compared with the experimental data, which are collected at five different irradiance levels (1000 W/m^{2}, 800 W/m^{2}, 600 W/m^{2}, 400 W/m^{2}, and 200 W/m^{2}) and three different temperature levels (25°C, 50°C, and 75°C). It can be seen that the curves of the ISDM fit the whole range of the experimental dataset. On the other hand, the errors of SDM seem larger at lower irradiance and higher temperature levels. With the experimental data, the RMSE of the current in SDM is calculated as 0.2837, while the RMSE of in ISDM is only 0.0776.
(a)
(b)
Figure 4 shows the absolute current errors of different performance predicting methods under different operating conditions. The curves denoted by the label “analytical SDM” are obtained from the analytical SDM model [4]. Ignoring the effect of incidence angle and air mass, the curves labeled by “analytical ISDM” denote the curves from De Soto’s analytical ISDM model [5]. It is evident the ISDM with the parameters extracted by the CS algorithm is more accurate than the analytical model. As for the SDM, the CS algorithm is capable of extracting a set of PV parameters with a good fit for the experimental data at the STCs. However, the SDM with the equations in the appendix does not exhibit a good prediction performance under other operating conditions.
(a)
(b)
To further validate the accuracy of the CS algorithm, the extracted parameters are compared to the ones obtained using GA in Figure 5. In general, the CS algorithm gives the better performance than GA for all cases. The Maximum Power Point (MPP), usually locating around 74% of the open circuit voltage, is an important technical data in PV modeling. However, a negative point of the GAbased ISDM is that the errors in the high voltage range are relatively high. The maximum absolute error of the GAbased ISDM is up to about 0.8 A, while the absolute error of the CS is kept below 0.2 A.
(a)
(b)
5. Conclusion
In this work, the Cuckoo Search (CS) algorithm is applied to estimate the parameters of two PV models, namely, Single Diode Model (SDM) and its improved version (ISDM). The feasibility of the proposed method has been validated by estimating the parameters of two commercial PV generators. The simulation and experimental results showed that the CS algorithm is capable of not only extracting all the parameters of the SDM under a certain condition but also successfully estimating all the parameters of ISDM under different operating conditions. In statistical analysis, CS algorithm recorded the lowest RMSE value compared to other algorithms such as GA, PSO and PS.
Appendix
PV Physical Model Adopted in PSIM
By using the parameters extracted at the STCs, the characteristics of a PV generator under nonstandard operating conditions can be calculated via the following equations:
Acknowledgments
The authors are grateful to Professor XinShe Yang for the sharing of Cuckoo Search source code online. Without his generosity, this work would not be possible. This research is supported by the National Natural Science Foundation of China under Grant 61070085.
References
 B. Amrouche, A. Guessoum, and M. Belhamel, “A simple behavioural model for solar module electric characteristics based on the first order system step response for MPPT study and comparison,” Applied Energy, vol. 91, no. 1, pp. 395–404, 2012. View at: Publisher Site  Google Scholar
 A. Orioli and A. di Gangi, “A procedure to calculate the fiveparameter model of crystalline silicon photovoltaic modules on the basis of the tabular performance data,” Applied Energy, vol. 102, pp. 1160–1177, 2013. View at: Google Scholar
 V. Lo Brano, A. Orioli, G. Ciulla, and A. di Gangi, “An improved fiveparameter model for photovoltaic modules,” Solar Energy Materials and Solar Cells, vol. 94, no. 8, pp. 1358–1370, 2010. View at: Publisher Site  Google Scholar
 M. G. Villalva, J. R. Gazoli, and E. R. Filho, “Comprehensive approach to modeling and simulation of photovoltaic arrays,” IEEE Transactions on Power Electronics, vol. 24, no. 5, pp. 1198–1208, 2009. View at: Publisher Site  Google Scholar
 W. De Soto, S. A. Klein, and W. A. Beckman, “Improvement and validation of a model for photovoltaic array performance,” Solar Energy, vol. 80, no. 1, pp. 78–88, 2006. View at: Publisher Site  Google Scholar
 K. Ishaque, Z. Salam, and H. Taheri, “Simple, fast and accurate twodiode model for photovoltaic modules,” Solar Energy Materials and Solar Cells, vol. 95, no. 2, pp. 586–594, 2011. View at: Publisher Site  Google Scholar
 V. Lo Brano, A. Orioli, and G. Ciulla, “On the experimental validation of an improved fiveparameter model for silicon photovoltaic modules,” Solar Energy Materials and Solar Cells, vol. 105, pp. 27–39, 2012. View at: Google Scholar
 A. N. Celik and N. Acikgoz, “Modelling and experimental verification of the operating current of monocrystalline photovoltaic modules using four and fiveparameter models,” Applied Energy, vol. 84, no. 1, pp. 1–15, 2007. View at: Publisher Site  Google Scholar
 M. A. de Blas, J. L. Torres, E. Prieto, and A. García, “Selecting a suitable model for characterizing photovoltaic devices,” Renewable Energy, vol. 25, no. 3, pp. 371–380, 2002. View at: Publisher Site  Google Scholar
 J. P. Charles, G. Bordure, A. Khoury, and P. Mialhe, “Consistency of the double exponential model with physical mechanisms of conduction for a solar cell under illumination,” Journal of Physics D, vol. 18, no. 11, pp. 2261–2268, 1985. View at: Publisher Site  Google Scholar
 D. S. H. Chan and J. C. H. Phang, “Analytical methods for the extraction of solarcell single and doublediode model parameters from IV characteristics,” IEEE Transactions on Electron Devices, vol. 34, no. 2, pp. 286–293, 1987. View at: Google Scholar
 J. C. H. Phang, D. S. H. Chan, and J. R. Phillips, “Accurate analytical method for the extraction of solar cell model parameters,” Electronics Letters, vol. 20, no. 10, pp. 406–408, 1984. View at: Google Scholar
 A. J. Joseph, B. Hadj, and A. L. Ali, “Solar cell parameter extraction using genetic algorithms,” Measurement Science and Technology, vol. 12, no. 11, pp. 1922–1925, 2001. View at: Publisher Site  Google Scholar
 W. Huang, C. Jiang, L. Xue, and D. Song, “Extracting solar cell model parameters based on chaos particle swarm algorithm,” in Proceedings of the International Conference on Electric Information and Control Engineering (ICEICE '11), pp. 398–402, April 2011. View at: Publisher Site  Google Scholar
 I. Fister, I. Fister Jr, X. S. Yang, and J. Brest, “A comprehensive review of firefly algorithms,” Swarm and Evolutionary Computation, 2013. View at: Publisher Site  Google Scholar
 M. F. AlHajri, K. M. ElNaggar, M. R. AlRashidi, and A. K. AlOthman, “Optimal extraction of solar cell parameters using pattern search,” Renewable Energy, vol. 44, pp. 238–245, 2012. View at: Publisher Site  Google Scholar
 A. Askarzadeh and A. Rezazadeh, “Parameter identification for solar cell models using harmony searchbased algorithms,” Solar Energy, vol. 86, pp. 3241–3249, 2012. View at: Google Scholar
 X. Yang and S. Deb, “Engineering optimisation by cuckoo search,” International Journal of Mathematical Modelling and Numerical Optimisation, vol. 1, no. 4, pp. 330–343, 2010. View at: Publisher Site  Google Scholar
 X. Yang and S. Deb, “Cuckoo search via Lévy flights,” in Proceedings of the World Congress on Nature and Biologically Inspired Computing (NABIC '09), pp. 210–214, December 2009. View at: Publisher Site  Google Scholar
 R. A. Messenger and J. Ventre, Photovoltaic Systems Engineering, CRC Press, New York, NY, USA, 2nd edition, 2004.
 PSIM User Manual, Powersim, Woburn, Mass, USA, 2001.
 R. N. Mantegna, “Fast, accurate algorithm for numerical simulation of Lévy stable stochastic processes,” Physical Review E, vol. 49, no. 5, pp. 4677–4683, 1994. View at: Publisher Site  Google Scholar
 A. H. Gandomi, X. Yang, and A. H. Alavi, “Cuckoo search algorithm: a metaheuristic approach to solve structural optimization problems,” Engineering with Computers, vol. 29, pp. 17–35, 2013. View at: Publisher Site  Google Scholar
 B. Birge, Particle Swarm Optimization Toolbox, http://www.mathworks.com/matlabcentral/fileexchange/7506particleswarmoptimizationtoolbox.
 P. Civicioglu and E. Besdok, “A conceptual comparison of the Cuckoosearch, particle swarm optimization, differential evolution and artificial bee colony algorithms,” Artificial Intelligence Review, vol. 39, pp. 315–346, 2013. View at: Publisher Site  Google Scholar
 T. Easwarakhanthan, J. Bottin, I. Bouhouch, and C. Boutrit, “Nonlinear minimization algorithm for determining the solar cell parameters with microcomputers,” International Journal of Solar Energy, vol. 4, pp. 1–12, 1986. View at: Google Scholar
 M. Bashahu and P. Nkundabakura, “Review and tests of methods for the determination of the solar cell junction ideality factors,” Solar Energy, vol. 81, no. 7, pp. 856–863, 2007. View at: Publisher Site  Google Scholar
 S. J. Jun and L. KaySoon, “Photovoltaic model identification using particle swarm optimization with inverse barrier constraint,” IEEE Transactions on Power Electronics, vol. 27, pp. 3975–3983, 2012. View at: Google Scholar
 K. Ishaque, Z. Salam, S. Mekhilef, and A. Shamsudin, “Parameter extraction of solar photovoltaic modules using penaltybased differential evolution,” Applied Energy, vol. 99, pp. 297–308, 2012. View at: Google Scholar
 X. S. Yang, Cuckoo Search Algorithm (Source Code), http://www.mathworks.com/matlabcentral/fileexchange/29809cuckoosearchcsalgorithm.
 Optimization Toolbox, The MathWorks Inc., http://www.mathworks.com/products/optimization/index.html.
 R. L. Burden and J. D. Faires, Numerical Analysis, Cengage Learning, Singapore, 2010.
Copyright
Copyright © 2013 Jieming Ma et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.