Characteristic data from R.T.C. France (Si solar cell).
This chapter presents a comprehensive study of a new hybrid method developed for obtaining the electrical unknown parameters of solar cells. The combination of a traditional method and a recent smart swarm-based optimization method is done, with a big focus on the application of the topic of artificial intelligence algorithms into solar photovoltaic production. The combined approach was done between the traditional method, which is the noniterative Levenberg-Marquardt technic and between the recent meta-heuristic optimization technic, called Grey Wolf optimizer algorithm. For comparison purposes, some other classical solar cell parameter determination optimization-based methods are carried out, such as the numerical (iterative, noniterative) methods, the meta-heuristics (evolution, human, physic, and swarm) methods, and other hybrid methods. The final obtained results show that the used hybrid method outperforms the above-mentioned classical methods, under this study.
- solar cell
- swarm-based intelligence
Solar photovoltaic energy is becoming the most popular renewable energy used in the world, at many caring of installations. Modeling and characterization are important topics that necessitate the determination of the exact solar photovoltaic (PV) cell’s unknown parameters values and thus optimizing the PV power generated. Solar PV generator’s performances are affected by many factors, where some of them are external, related to the environmental conditions like the weather’s variations (irradiation and temperature), shading phenomenon, hotspot , dust, cell damage, wind velocity, and soiling . Others are internal, related to the electrical, physical, and mathematical modeling. After the modeling step of any PV generator, their identified parameters values are used in the established model. Therefore, it is necessary to find the accurate values of the unknown electrical PV parameters by an appropriate approach. Besides, accurate parameter values of PV cells are essential for the development of good controlling techniques for Maximum Power Point Tracking (MPPT) based power electronic converters . As shown in the Figure 1 the importance of PV parameters’ obtaining accurate values for a whole PV system.
With the complexity of the nonlinearity founded from the current-voltage relationship through the model used to represent the physical behavior of PV cells (Single Diode, Double Diode, Three Diode, and so on) . The parameters to be found become more numerous, as the complexity of models increases. The PV parameters of a Single Diode Model (SDM), which are the most determined in literature, and which are the light and the diode saturation currents, the diode ideality factor, the series, and the shunt resistances. A challenge is to obtain the values of all the PV model’s parameters value while keeping a reasonable compromise of some criteria, such as the fast speed of convergence, low implementation complexity, and so on.
Artificial intelligent (AI) algorithms have attracted attention lately, by the scientific community to be used for resolving many topic’s problems. Among them solar photovoltaics’ problems, such as modeling, identification, prediction, optimization, sizing, control, and many others . The hybrid optimization-based methods have attracted more attention to achieve more efficiency and precision. For this reason, this paper presents a combination of the traditional Levenberg-Marquardt (LM) with the recent meta-heuristic Grey Wolf optimizer (GWO) method. This hybrid LMGWO method has seemed to be the most performing, which we finally have chosen to be used in this work.
The remaining of this chapter is structured as follow. After the introduction given in this Section 1. Section 2 presents a classification of the optimization-based methods used in the literature to estimate the PV parameters values. In the Section 3, models are presented of the PV cell (single diode), and the problem of PV parameters obtaining values is presented. Section 4 gives some details about the hybrid used method to estimate the SDM parameters values. Moreover, this section provides the basic concept of the LM and the GWO. Finally, Section 5 shows some results comparing one method of each type proposed in the classification and the LMGWO. Some conclusions are given in Section 6.
2. Classification and literature review
Earliest, numerous research workings have been developed only for obtaining parasitic resistances (series and shunt) values  by the cause of their high influences on the PV’s performances. Then, it has been observed some influences of all PV’s electrical parameters on the PV’s performances , which leads the researchers, for doing a large number of studies for obtaining their accurate values. In literature, different approaches, that allow the evaluation of the PV cell’s electrical parameters values, exist [8, 9, 10, 11]. There are some analytical approaches [12, 13, 14, 15], and those based on the optimization process. By the cause of limits of the analytical method to achieve with high precision the PV parameters values, our interest is gone for the optimization-based methods. This latter can be classified as in Figure 2.
The optimization algorithms are classified into traditional, heuristic/metaheuristic, and hybrid groups. More details about each group of optimization-based are given at the following subsections.
2.1 Numeric traditional
The numeric traditional optimization-based methods are used to find the optimum of a function using gradient or hessian. These numeric-traditional methods applied for PV parameters obtaining values, are based on the reduction of the number of parameters to be evaluated, such as Kashif’s one . In this subsection’s methods, the traditional iterative Newton-Raphson (NR) approach , iterative curve-fitting , can also be used. It necessitates an iterative process with good initialization guess of PV parameters values, to converge to the best solutions. Others build a set of nonlinear transcendental equations (based on short-circuit, open-circuit, MPP, derivatives of the I-V curve) and execute an optimization problem instead of solving by numerical methods . For the noniterative method, the Levenberg-Marquardt  can be cited.
Even though with their effectiveness to get a good local search, they still have other limitations, such as the need of a convex, continuous, and differentiability of the objective function. Besides, good guessing of initial parameters values is necessary for a good converging process. Also, as the complexity of the modeling process increase, as the optimizer loses the ability for obtaining better results.
In recent times, meta-heuristic optimization-based methods, using Artificial-Intelligence (AI) inspired algorithms, have attracted the care of researchers to obtain with good precision, the unknown PV parameters values.
The metaheuristic methods use bio-inspired algorithms in the search process to identify the PV parameters values at real-time, using the errors between the real experimental data and the simulated data. These approaches are based on an experimental process and are known as identification methods . These approaches are graphically based on curve characteristics fitting.
Meta-heuristics are categorized into four main sets such as evolution-based , physic-based , immune-human-based  and swarm-based intelligence methods . Some of each category is used for obtaining PV parameters values as presented on what follow.
Evolutionary Algorithm (EA) , Differential Evolutionary (DE) , Genetic Algorithms (GAs) , Pattern Search (PS) , Simulated Annealing (SA) , Improved Shuffled Complex Evolution (ISCE) , Repaired Adaptive Differential Evolution (Rcr-IJADE) .
Electromagnetic Field Optimization (EFO), Gravitational Search Algorithm (GSA), Electromagnetism-Like Algorithm (EMA), Weighted Superposition Attraction (WSA) .
Harmony Search (HS) , Bacterial Foraging Algorithm (BFA) , Simplified Teaching-Learning-Based Optimization (STLBO) , Discrete Symbiosis Organism Search (DSOS) , Artificial Immune system (AIS) .
The swarm-based, Particle Swarm Optimization (PSO) [3, 35, 36], Bird Mating Optimization (BMO) , Artificial Bee Swarm Optimization (ABSO) . Grey Wolf Optimizer (GWO) , Chaotic Whale Optimization Algorithm (CWOA) , Cat Swarm Optimization (CSO) , and Cluster Analysis (CA) .
The metaheuristics are more attractive than the deterministic traditional methods in terms of accuracy and robustness, by the cause of their good global research achieving. Besides, they do not require a gradient or differentiable of the objective function. Besides, the initial guess of parameters values is not a necessity, but it necessitates the upper and lower limits of an interval of research.
The hybrid method combines different approaches. These methods make a mix of other methods, i.e. analytical and numeric-traditional methods ; analytical and meta-heuristics, numeric-traditional and meta-heuristics optimization; a combination of two different meta-heuristics, etc. . We can site, hybrid adaptive Nelder-Mead simplex algorithm based on eagle strategy (EHA-NMS) , Nelder-Mead simplex algorithm based on eagle strategy (EHA-NMS) , Nelder-Mead and Modified Particle Swarm Optimization (NM-MPSO) , Artificial Bee Colony-Differential Evolution (ABC-DE) , Trust-Region Reflective deterministic algorithm with the Artificial Bee Colony (ABC-TRR) , Teaching–learning–based Artificial Bee Colony (TLABC) . Our proposed Levenberg-Marquardt with Grey Wolf optimizer (LM-GWO), and so on. Those methods, which are called hybrid, have excellent performances because they restrict the universe in the search process without losing precision (without losing the optimum). They achieve outstanding results with a smaller number of iterative steps when compared with pure optimization methods.
3. Modeling and problem formulation
There are several electrical models, used by researchers, to describe the physical behaviors of PV cells. The Single Diode Model, containing the five unknown parameters, used in this paper is represented in Figure 3. By the cause of compromise between accuracy and simplicity, the SDM is selected herein.
The mathematical expressions related to the current-voltage, (I-V) relationship of the PV cell is as follow.
The overhead mathematical equation is in a nonlinear form and has a set of five unknown parameters (
Several approaches permit the formulation of the optimal nonlinear PV parameters determination problem, using the error (between real and simulated data) .
Our focus is to estimate the PV parameters values of the SDM model using RTC France data at the conditions of irradiance about 1000 W/m2 and of temperature about 300°C. We do not review the identification process as detailed on our previous work ; our focus is restricted on the third part of identification process, which is the estimation of PV parameters values. The big focus is to optimize the damping factor of LM through GWO. The characteristics of RTC France Silicon-cell data from datasheet are presented on the following Table 1.
|Characteristic data||R.T.C France|
The real experimental data used of RTC France are presented on the following Table 2.
|Measurement||V (Volts)||I (Ampere)|
4. Hybrid optimization-based method
Hybrid optimization-based algorithms have become the modern choice for resolving challenging problems [41, 42, 43]. A compromise is gotten in this work, from a combination of a traditional numeric optimization-based with a metaheuristic swarm-based method.
The estimation/identification process can be gotten in three major steps, such as the initial step of prediction through the use of least-squares mean (LSM), the getting of optimal PV parameters values through Levenberg-Marquardt (LM), and the optimization of a dominant factor through GWO as detailed below.
4.1 Least squares mean (initial phase of prediction)
Prediction of initial PV parameters values using LSM [44, 45] for the two parts of the introduced real experimental points of I-V curve characteristics as described below.
For the linear part:
The prediction in the linear part [46, 47] of the model can be obtained simply through the use of the following expressions.
For the nonlinear part:
The prediction in the nonlinear part [19, 48] of the model can be obtained with a logarithmic way through the use of the following logarithmic expression.
Once obtaining initial values of PV parameters values, we introduce them on the LM in order to optimize their values, as explained in the following subsection.
4.2 Levenberg Marquardt (get of optimal PV parameter values)
The traditional Levenberg-Marquardt approach is a gradient order from Steepest-Descent (SD) in its first step and from Gauss-Newton (GN) in its second step [48, 49, 50]. It is mainly based on an optimization of the error between real data and data from the model through the following expression.
The real and simulated data are denoted by
Evaluate the objective function
Calculus of Jacobian of
For (damping optimized) update
The dominant factor λ is considered as responsible parameters for switching from SD to GN in the LM process .
For this reason, it is important to get an optimal value of this damping factor by the use of another optimization-based method, our choice was for the recent swarm-based method called GWO, through the following idea:
In addition, it is mentioned that at each iteration of the LM process that the damping factor must be found and is considered as crucial factor for the convergence process of the algorithm. Therefore, its value must be optimized by the use of another approach such as the GWO approach.
4.3 Grey Wolf optimizer (optimization of damping factor’s value)
In this subsection, our focus is on the evolution of the function
The meta-heuristic methods are known for their simplicity, flexibility, derivation free process and the ability to find the global optimal solution. They are also appropriate for a diversity of problems without changing on their main structure. These methods can be based on a single solution or on population of solutions. The basic concepts can be obtained through exploration (exploring all of the search space and thus avoiding local optimum) and exploitation (investigating process in detail of the promising search space area).
Swarm-based intelligence (SI) methods, which derive from meta-heuristics, are based on the smart collective behavior of decentralized and self-organized swarms to ensure some biological needing such as food or security. A detailed discussion about the recent smart swarm-based algorithm, known as GWO is presented as follow.
Grey Wolf optimizer (GWO) algorithm, developed by Mirjalili in 2014, is a recent smart swarm-based meta-heuristic approach [50, 51, 52]. This algorithm mimics the leadership hierarchy and hunting process of Grey wolves in the wildlife. The following points represent the hierarchy in a wolf’s group, which is about 5 to 12 members.
The alphas wolves (α): are the leading wolves that are responsible for managing and making decisions. These are the first level of the wolves’ social hierarchical structure. This later is presented in Figure 4.
The betas wolves (β): represent the second level. Their main job is to help and support alpha’s decisions.
The deltas wolves (δ): represent the third level in the pack and are called subordinates. They use to follow alpha and beta wolves. The delta wolves can divide their tasks into five categories as follows:
Scouts: used to control the boundaries of the territory and alert the pack in case of danger.
Sentinels: protect and guarantee the safety of the pack.
Elders: among these strong and mature wolves, some of them become either alpha or beta.
Hunters: help alpha and beta in the hunting prey, providing food to the pack.
Caretakers: responsible for caring the ill, wounded and weak wolves.
The omegas wolves (ω): represent the lowest level. They have to follow alpha, beta and delta wolves.
When a pack of wolves sees a prey such as (gazelle, rabbit or a buffalo) they attack it in three steps and do not recede, Figure 5.
These three steps of the hunting process can be mentioned as follows.
Encircling, tracking, chasing, and approaching the prey (Figure 5: A, B).
Pursuing, encircling, and harassing the prey until it stops moving (Figure 5: C).
Attacking the prey (Figure 5: D, E).
The mentioned above social hierarchy and hunting process of Grey wolves have been mathematically modeled in GWO, as follows [51, 52]:
The first, second and third best solutions are considered as α, β and δ wolves, respectively.
The rest of the candidate solutions are considered as ω.
The following equations are used to model the encircling first step of Grey wolves hunting process:
To mathematically simulate the second step of the Grey wolves hunting process, we suppose that the alpha (best candidate solution), beta and delta have a better knowledge about the potential location of the prey . Therefore, the first three best solutions obtained so far are saved and oblige the other search agents (including the omegas) to update their positions according to the position of the best search agents. In this regard, the following formulas are used.
The final third step is the hunting process as attacking the prey as soon as it stops moving.
The main steps of the used hybrid LMGWO method applied for the PV parameters obtaining values are presented in Figure 6.
5. Results and discussions
The following Table 3 presents PV parameters results for the all classified optimization-based method discussed in Section 2.
|LM ||0.760782||3.166611e-07||1.479182||0.036461||53.271523||9.8680e- 4|
|PSO [35, 36]||0.76077||0.32454||1.48165||0.03636||53.8550||9.8606E−04|
From the above Table 3 it is clear that for the traditional methods, the LM is more accurate than Newton’s method, which in turn outperforms Kashif’s method. Then, for the metaheuristic methods for each of their category as follow.
It is observed that ISCE, Rcr-IJADE, and PCE outperform PS, which in turn is better than GA and SA.
All the swarm-based used outperform ABSO in terms of precision.
CSO and STBLO outperform HS, which is better than BFA.
EFO is more accurate than EMA, which is more accurate than WSA, which is more accurate then GSA.
It is mentioned that the swarm-based got the best results compared to the other metaheuristic’s category.
Finally, for the hybrid methods, it is clear that all of them have achieved the highest best optimized (minimum) values for RMSE, until now with the value of (9.8601E−04).
In addition, the hybrid methods outperform the metaheuristics, which in turn outperform the traditional methods.
The fitting obtained curves of real and simulated data, using the proposed LMGWO are illustrated in Figure 7.
The best approximation gotten from the fitted curves in Figure 7 has proved the effectiveness of our hybrid LMGWO method.
Since nonlinear and multi-parameter PV models are used to represent a PV cell’s physical behavior, classical methods are incapable of evaluating the exact parameters’ values of these models. For these reasons, the present paper presents a proposed hybrid method of obtaining the unknown electrical parameters of solar photovoltaic cells. To do so, we applied our hybrid method, the LM combined with GWO method, after having initial guess using least squares mean, and then compared it with other previous optimization-based methods. The application of LMGWO has shown high precision for the obtained solutions’ values. The LMGWO outperforms the other tested algorithms in many aspects. It is simple and accurate and converges rapidly to the optimum in every test. In addition, it has fewer parameters to set then it is easily implemented. The obtained results demonstrate the efficiency of the hybrid LMGWO approach compared to the other meta-heuristics and some of the other traditional methods.
Conflict of interest
The authors declare no conflict of interest.
Guerriero P, Daliento S. Toward a hot spot free PV module. IEEE Journal of Photovoltaics. May 2019; 9(3):796-802
Coello M, Boyle L. Simple model for predicting time series soiling of photovoltaic panels. IEEE Journal of Photovoltaics. September 2019; 9(5):1382-1387
Jieming M. Optimization approaches for parameter estimation and maximum power point tracking (MPPT) of photovoltaic systems [thesis]. Liverpool: University of Liverpool for the degree of Doctor in Philosophy; 2014
Saha C, Agbu N, Jinks R, Nazmul Huda M. Review article of the solar PV parameters estimation using evolutionary algorithms. MOJ Solar and Photoenergy Systems. 2018; 2(2):66-78
Leva S, Ogliari E. Computational intelligence in photovoltaic systems. Applied Sciences. 2019; 9:1826. DOI: 10.3390/app9091826
Bashahu M, Habyarimana A. Review and Test of Methods for Determination of the Solar Cell Series Resistance. Pergamon: Renewable Energy; 1995
Kenneth KL. Analysis of performance degradation in CdS solar cells. IEEE Transactions on Aerospace and Electronic Systems. 1969; AES-5(6):912-917
Tchoketch-Kebir S, Haddadi M, Ait-Cheikh MS. An overview of solar cells parameters extraction methods. In: IEEE Proceedings in the 3rd International Conference on Control, Engineering & Information Technology CEIT, Telmcen, Algeria, 2015
Humada Ali M, Hojabri M, Mekhilef S, Hamada Hussein M. Solar cell parameters extraction based on single and double-diode models: A review. Renewable and Sustainable Energy Reviews. 2016; 56:494-509
Pillai Dhanup S, Rajasekar N. Metaheuristic algorithms for PV parameter identification: A comprehensive review with an application to threshold setting for fault detection in PV systems. Renewable and Sustainable Energy Reviews. [Accepted: 28 October 2017]
Abbassia R, Abbassic A, Jemlic M, Chebbi S. Identification of unknown parameters of solar cell models: A comprehensive overview of available approaches. Renewable and Sustainable Energy Reviews. 2018; 90:453-474
Shongwe S, Hanif M. Comparative analysis of different single-diode PV modeling methods. IEEE Journal of Photovoltaics. 2015; 5(3)
Mehta Hitesh K, Warke H, Kukadiya K, Panchal Ashish K. Accurate expressions for single-diode-model solar cell parameterization. IEEE Journal of Photovoltaics. 2019; 9(3)
Chitti Babu B, Gurjar S. A novel simplified two-diode model of photovoltaic (PV) module. IEEE Journal of Photovoltaics. 2014; 4(4):1156-1161
Brano VL, Ciulla G. An efficient analytical approach for obtaining a five parameters model of photovoltaic modules using only reference data. Applied Energy. 2013; 111:894-903
Kashif I, Zainal S, Hamed T. Simple, fast and accurate two-diode model for photovoltaic modules. Solar Energy Materials & Solar Cells. 2011; 95:586-594
Cabestany J, Castanier L. Evaluation of solar cell parameters by nonlinear algorithms. Journal of Physics D: Applied Physics. 1983; 16:2547
Easwarakhanthan T, Bottin J, Bouhouch I, Boutrit C. Nonlinear minimization algorithm for determining the solar cell parameters with microcomputers. International Journal of Solar Energy. 1986; 4(1):1-12
Dkhichi F, Oukarfi B. Levenberg-Marquardt algorithm for parameter identification of solar cell model. In: Conference Paper June. 288 014. pp. 781-788
Tchoketch_Kebir S, Ait_Cheikh MS, Haddadi M. A detailed step-by-step electrical parameters identification method for photovoltaic generators using a combination of two approaches. ASTES Journal. 2018; 3(4):44-52
AlHajri MF, El-Naggar KM, AlRashidi MR, Al-Othman AK. Optimal extraction of solar cell parameters using pattern search. Renewable Energy. 2012; 44:238-245
Rashedi E, Nezamabadi-Pour H, Saryazdi S. GSA: A gravitational search algorithm. Information Sciences. 2009; 179:2232-2248
Chena M-H, Chang P-C, Wu J-L. A population-based incremental learning approach with artificial immune system for network intrusion detection. Engineering Applications of Artificial Intelligence. 2016
Beni G, Wang J. Swarm intelligence in cellular robotic systems. In: Robots and Biological Systems: Towards a New Bionics. Riverside, California 92521, USA: Springer; 1993. pp. 703-712
Gao X, Cui Y, Hu J, Xu G, Wang Z, Qu J. Parameter extraction of solar cell models using improved shuffled complex evolution algorithm. Energy Conversion and Management. 2018; 157:460-479
El-Naggar KM, AlRashidi MR, AlHajri MF, Al-Othman AK. Simulated annealing algorithm for photovoltaic parameters identification. Solar Energy. 2012; 86:266-274
Chen Z, Wu L, Lin P, Wu Y, Cheng S. Parameters identification of photovoltaic models using hybrid adaptive Nelder-Mead simplex algorithm based on eagle strategy. Applied Energy. 2016; 182:47-57
Gong W, Zhihua C. Parameter extraction of solar cell models using repaired adaptive differential evolution. Solar Energy. 2013; 94
Yurtkuran A, Küçükoğlu İ. Comparative study of physics-inspired meta-heuristic algorithms for the solar cell parameter identification problem. In: 16th International Conference on Clean Energy (ICCE-2018)
Askarzadeh A, Rezazadeh A. Parameter identification for solar cell models using harmony search-based algorithms. Solar Energy. 2012; 86:3241-3249
Rajasekar Neeraja N, Kumar K, Venugopala R. Bacterial foraging algorithm based solar PV parameter estimation. Solar Energy. 2013; 97:255-265
Niu Q, Zhang H, Li K. An improved TLBO with elite strategy for parameters identification of PEM fuel cell and solar cell models. International Journal of Hydrogen Energy. 2014
Bouali C, Schulte H, Mami A. A high performance optimizing method for modeling photovoltaic cells and modules array based on discrete symbiosis organism search. Energies. 2019; 12:2246. DOI: 10.3390/en12122246
Jacob B, Balasubramanian K, Sudhakar Babu T, Azharuddin SM, Rajasekar N. Solar PV modelling and parameter extraction using artificial immune system. Energy Procedia. 2015; 75:331-336
Wang X, Xu Y, Ye M. Parameter extraction of solar cells using particle swarm optimization. Journal of Applied Physics. 2009; 105:094502
Soon JJ, Low KS. Photovoltaic model identification using particle swarm optimization with inverse barrier constraint. IEEE Transactions on Power Electronics. 2012; 27:3975-3983
Askarzadeh A, Rezazadeh A. Extraction of maximum power point in solar cells using bird mating optimizer-based parameters identification approach. Solar Energy. 2013; 90:123-133
Chen X, Xu B, Mei C, Ding Y, Li K. Teaching–learning–based artificial bee colony for solar photovoltaic parameter estimation. Applied Energy. Elsevier. 2018; 212:1578-1588
Robandi DI. Photovoltaic parameter estimation using Grey Wolf optimization. In: 3rd International Conference on Control, Automation and Robotics, 2017
Oliva D, Abd El Aziz M, Ella Hassanien A. Parameter estimation of photovoltaic cells using an improved chaotic whale optimization algorithm. Applied Energy. 2017; 200:141-154
Guo L, Meng Z, Sun Y, Wang L. Parameter identification and sensitivity analysis of solar cell models with cat swarm optimization algorithm. Energy Conversion and Management. 2016; 108:520-528
Hamid NFA, Rahim NA, Selvaraj J. Solar cell parameters identification using hybrid Nelder-Mead and modified particle swarm optimization. Journal of Renewable and Sustainable Energy. 2016
Hachana O, Hemsas KE, Tina GM, Ventura C. Comparison of different metaheuristic algorithms for parameter identification of photovoltaic cell/module. Journal of Renewable and Sustainable Energy. 2013; 5
Miller SJ. The Method of Least Squares. Providence, RI: Mathematics Department Brown University. p. 02912
Lawson C, Hanson R. Solving Least Squares Problems. Society for Industrial and Applied Mathematics, Classics in applied mathematics, 15; SIAM e-books. Philadelphia, Pa, Englewood Cliffs, NJ: Prentice-Hall; 1974
Kong M, Li D, Zhang D. Research on the application of improved least square method in linear fitting. IOP Conference Series: Earth and Environmental Science. 2019; 252:052158. DOI: 10.1088/1755-1315/252/5/052158
Whittle P, Sargent T. Least-square approximation. In: Prediction and Regulation by Linear Least-Square Methods. 2nd ed. Minneapolis: University of Minnesota Press; 1983. pp. 46-55
Walling DD. Numerical methods for non-linear least squares curve fitting [thesis]. Iowa State University; 1963
Gavin HP. The Levenberg-Marquardt Method for Nonlinear Least Squares Curve-Fitting Problems. Department of Civil and Environmental Engineering, Duke University; 2016
Manolis I, Lourakis A. A Brief Description of the Levenberg-Marquardt Algorithm Implemened by levmar. Crete, Greece: Institute of Computer Science, Foundation for Research and Technology - Hellas (FORTH); 2005
Mirjalili S, Mirjalili SM, Lewis A. Grey Wolf optimizer. Advances in Engineering Software. 2014; 69:46-61
Mirjalilia S, Saremia BS, Mirjalilic BSM, Coelhod LDS. Multi-objective grey wolf optimizer: A novel algorithm for multi-criterion optimization. Expert Systems with Applications. 2016; 47:106-119
Lakshminarayanan S. Nature inspired Grey Wolf optimizer algorithm for minimizing operating cost in green smart home [thesis]. The University of Toledo; 2015