# Parameter Estimation of the Photovoltaic System Using Bald Eagle Search (BES) Algorithm

Authors: Ndongmo Fotsa Nicaire; Perabi Ngoffe Steve; Ndjakomo Essiane Salome; Abessolo Ondoua Grégroire
Journal: International Journal of Photoenergy
 (2021)

---

## Abstract

The global demand for renewable energy is growing, and one of the proposed solutions to this energy crisis is the use of photovoltaic systems. So far, they are a reliable solution, as they are nonpolluting and can be used almost anywhere on the planet. However, the design and development of more efficient photovoltaic cells and modules require an accurate extraction of their intrinsic parameters. Up to date, metaheuristic algorithms have proven to be the best methods to obtain accurate values of these intrinsic parameters. Hence, to extract these parameters reliably and accurately, this paper presents an optimization method based on the principle of bald eagle search (BES) during fish hunting. This search is divided into three steps: in the first stage (space selection), the eagle selects the space with the largest number of prey; in the second stage (space search), the eagle moves into the selected space to search for prey; in the third stage (dive), the eagle swings from the best position identified in the second stage and determines the best point to hunt. Thus, we used the proposed BES algorithm to determine the parameters of the single-diode model (SDM), the double-diode model (DDM), and the PV modules. This algorithm converges very quickly and gives a root mean square error (RMSE) of9.8602e−04for the single-diode model and 9.8248e−4 for the dual-diode model. The results obtained show that the proposed algorithm is more efficient than the other methods available in the literature, in terms of the better accuracy of the results obtained. The good harmony of the I-V and P-V characteristic curve of the calculated parameters with that of the measured data from a PV module/cell data sheet proves that the proposed BES should be used among the methods provided in the literature for the identification of PV solar cell parameters.


---

## Body

## 1. Introduction

The energy demand of almost every country in the world is increasing due to its large-scale industrial expansion, population growth, and the continuous growth of per capita energy consumption. It is worth noting that most of the energy needs are in the form of electricity. In contrast, the use of fossil fuel-based electricity generation has reached saturation levels due to increased environmental concerns and limited resources. Thus, the gaps between demand and production in the future must be filled by renewable energy sources [1]. However, as the solar energy obtained from a solar PV module is not constant, a major challenge is therefore to maximise the use of solar energy due to the unpredictability of the power output of PV modules caused by the resulting variations in irradiance levels and cell temperature [2]. Thus, the competition to optimize and increase the efficiency of photovoltaic cells has led researchers to find methods to determine the intrinsic parameters of these cells. In the literature review, several methods have been proposed for the extraction of the parameters; each of these methods has drawbacks, either in terms of complexity of use and accuracy or in terms of convergence and speed. These methods are classified into three categories: analytical, numerical, and metaheuristic methods [3].In the analytical method, a set of transcendental equations is solved to estimate the parameters of the solar cell. The main advantage of the analytical method is the speed of calculation and relatively accurate results. Analytical methods are simple, with a short computation time. Sometimes, a single iteration is sufficient to obtain the result [4]. Although this approach is very popular, it is not always easy to apply. In addition, they need many data points of the I-V curve, which in turn complicates the computation [5]. Analytical methods include Lambert’s W function [6], Taylor series expansion [7], and Chebyshev polynomials [8, 9]. Some of these methods estimate 5 parameters, and others extract only the series and shunt resistance. The main weakness of analytical methods is that they are only suitable for standard conditions; consequently, they have poor results with variable ambient conditions [10]. Numerical methods with curve-fitting techniques are better than analytical methods. The algorithms of these methods provide accurate results by evaluating all points of the PV-IV curves using the algorithm [10].Metaheuristic algorithms are global optimization techniques which do not impose any restrictions on the problem formulation and have the ability to solve various complex problems [11]. In the literature, many metaheuristic algorithms have been suggested for extracting the parameters of PV solar cell models, such as the Genetic Algorithm (GA) [12, 13], the Cuckoo Search (CS) Algorithm [11, 14], Particle Swarm Optimization (PSO) [15–17], Differential Evolution (DE) algorithm [18], Artificial Bee Colony (ABC) algorithm [19–21], Artificial Algorithm of Bee Swarm Optimization (ABSO) [22], Bacterial Foraging Optimization (BFO) algorithm [23, 24], Biogeography-Based Optimization (BBO) algorithm [25], Floral Pollination Algorithm (FPA) [26, 27], Jaya Optimization Algorithm (JAYA) [28, 29], Salp Swarm Algorithm (SSA) [30], Bird Mating Optimization (BMO) algorithm [31], Teaching-Learning-Based Algorithm (TLBO) [20, 32–34], Whale Optimization Algorithm (WOA) [35–37], Backtracking Search Algorithm (BSA) [38], Sine-Cosine Algorithm (SCA) [39], Imperialist Competitive Algorithm (ICA) [40, 41], Multiverse Optimizer (MVO) algorithm [42], Ant-Lion Optimizer (ALO) algorithm [43, 44], Eagle Strategy (ES) [45], Cat Swarm Optimization (CSO) [46], Harmony Search (HS) [47], Firefly Algorithm (FA) [48], Simplified Swarm Optimization (SSO) [49], Moth-Flame Optimization (MFO) algorithm [50], Water Cycle Algorithm (WCA) [51], Enhanced Vibration of Particles System (EVPS) [52], Harris Hawks Optimization (HHO) [53], Shuffled Frog Leaping (SFL) algorithm [54], Metaphor-Free Dynamic Spherical Evolution (DSCE) [55], enhanced metaphor-free gradient-based optimizer (EGBO) [56], delayed dynamic step shuffling frog-leaping algorithm (DDSFLA) [57], evolutionary shuffled frog leaping with memory pool (SFLBS) [58], enhanced spherical evolution (LCNMSE) [59], random reselection particle swarm optimization (PSOCS) [60], boosting slime mould algorithm (CNMSMA) [61], and Teaching-Learning Artificial-Based Bee Colony (TLABC) [20]. Comparative with analytical and numerical approaches, these metaheuristic algorithms were able to provide satisfactory results for the extraction of PV model parameters. However, these metaheuristic algorithms still have inherent drawbacks. For example, HS is very sensitive to the initial population, PSO is easily subjected to premature convergence, ABC is poor when the system is in operation, and CS suffers from slow convergence [11].This paper presents a comprehensive study on the estimation of design parameters for SDM, DDM, and PV modules using the bald eagle search (BES) algorithm [62]. The paper is organised as follows: Section 1, introduction; Section 2, mathematical models and analysis of different models for cells, modules, and PV; Section 3, illustration of the formulation and optimization problem of models; in Section 4, we present the methodology for applying the optimization algorithm to estimate the model parameters; Section 4 presents the results and discussions and finally we conclude the paper.The main contribution of our work is that we are proposing a new algorithm (BES) to extract the parameters of different cell models (single diode, double diode) and photovoltaic panels. Comparisons were made between the BES and the different algorithms on the extraction of the parameters of the three models. Comparing the results obtained using BES with those obtained through other methods, the accuracy and the reliability of the results can be clearly observed. Thus, BES can be an effective alternative for the parameter extraction from PV models.


## 2. Photovoltaic Model Description

In practical applications, a single-diode model (SDM) and double-diode model (DDM) are commonly used to describe the nonlinear voltage–current characteristics of photovoltaic systems. This section describes the properties of each of these two models [3, 52, 57].


### 2.1. Single-Diode Model

The single-diode model is shown in Figure1. This model is the most used in many researches and better than the two-diode models (described after) due to its behaviour which is closer to a PV cell than the series resistance model (simplified) and its simplicity for the mathematical calculation [29].
(1)I=Iph−Id−Ish,(2)Ish=V+Rs×IRsh,(3)I=Iph−IsdexpqV+Rs×In×K×T−1−V+Rs×IRsh,Figure 1
Single-diode equivalent model.withId as the saturation current of the diode, q the charge of one electron q=1.60217646×10−19C, n the ideality factor of the diode, K Boltzmann’s constant K=1.3806503×10−23J/K, and T the temperature in Kelvin.The equation model (3) is also called the implicit model with five unknown parameters: Iph, Isd, Rs, Rsh, and n.


### 2.2. Two-Diode Model

Due to its simplicity and accuracy, the above-mentioned single-diode model has been widely used to describe the static characteristics of the photovoltaic cell. However, the single-diode model has inherent drawbacks as it assumes that the ideality factor of the diode remains constant throughout the range of output voltage variation [29]. Currently, the closest electrical model to a photovoltaic cell is the two-diode (double exponential) model, where the cell is of course presented as an electrical current generator whose behaviour is equivalent to a current source with two diodes in parallel. The two-diode model is shown in Figure 2 [18, 52].
(4)I=Iph−Id1−Id2−Ish,Ish=V+Rs×IRsh,I=Iph−Isd1expqV+Rs×In1×K×T−1−Isd2expqV+Rs×In2×K×T−1−V+Rs×IRsh.Figure 2
Equivalent model with two diodes.The parameters of this double-diode model to be estimated areIph, Isd1, Isd2, Rs, Rsh, n1, and n2.


### 2.3. The PV Model

The PV module model is shown in Figure3 [52, 63, 64].
(5)I=Iph×Np−Isd×NpexpqV/Ns+Rs/Np×In×K×T−1−VNp/Ns+Rs×IRsh,Figure 3
Equivalent model of a PV panel.whereNp and Ns represent the number of solar cells in parallel and in series, respectively. Thus, for this PV model, five unknown parameters (Iph, Isd, Rs, Rsh, and n) have to be identified.


### 2.4. Objective Function

For a single-diode model, the objective function is expressed as(6)fV,I,X=Iph−IsdexpqV+Rs×In×K×T−1−V+Rs×IRsh,X=Iph,Isd,Rs,Rsh,n.For the double-diode model, the objective function is(7)fV,I,X=Iph−Isd1expqV+Rs×In1×K×T−1−Isd2expqV+Rs×In2×K×T−1−V+Rs×IRsh,X=Iph,Isd1,Isd2,Rs,Rsh,n1,n2.The parameters can be estimated by minimizing the objective functionRMSEX, i.e., by searching for the solution vector x [11, 16, 29, 65, 66].
(8)RMSEX=1N∑k=1NfV,I,X2,whereX represents the parameters composed of the solution vector. V and I are the measured voltage and current, respectively. N represents the number of experiments. Hence, to estimate the parameters is equivalent to search the X in the range which minimizes the objective function.


## 2.1. Single-Diode Model

The single-diode model is shown in Figure1. This model is the most used in many researches and better than the two-diode models (described after) due to its behaviour which is closer to a PV cell than the series resistance model (simplified) and its simplicity for the mathematical calculation [29].
(1)I=Iph−Id−Ish,(2)Ish=V+Rs×IRsh,(3)I=Iph−IsdexpqV+Rs×In×K×T−1−V+Rs×IRsh,Figure 1
Single-diode equivalent model.withId as the saturation current of the diode, q the charge of one electron q=1.60217646×10−19C, n the ideality factor of the diode, K Boltzmann’s constant K=1.3806503×10−23J/K, and T the temperature in Kelvin.The equation model (3) is also called the implicit model with five unknown parameters: Iph, Isd, Rs, Rsh, and n.


## 2.2. Two-Diode Model

Due to its simplicity and accuracy, the above-mentioned single-diode model has been widely used to describe the static characteristics of the photovoltaic cell. However, the single-diode model has inherent drawbacks as it assumes that the ideality factor of the diode remains constant throughout the range of output voltage variation [29]. Currently, the closest electrical model to a photovoltaic cell is the two-diode (double exponential) model, where the cell is of course presented as an electrical current generator whose behaviour is equivalent to a current source with two diodes in parallel. The two-diode model is shown in Figure 2 [18, 52].
(4)I=Iph−Id1−Id2−Ish,Ish=V+Rs×IRsh,I=Iph−Isd1expqV+Rs×In1×K×T−1−Isd2expqV+Rs×In2×K×T−1−V+Rs×IRsh.Figure 2
Equivalent model with two diodes.The parameters of this double-diode model to be estimated areIph, Isd1, Isd2, Rs, Rsh, n1, and n2.


## 2.3. The PV Model

The PV module model is shown in Figure3 [52, 63, 64].
(5)I=Iph×Np−Isd×NpexpqV/Ns+Rs/Np×In×K×T−1−VNp/Ns+Rs×IRsh,Figure 3
Equivalent model of a PV panel.whereNp and Ns represent the number of solar cells in parallel and in series, respectively. Thus, for this PV model, five unknown parameters (Iph, Isd, Rs, Rsh, and n) have to be identified.


## 2.4. Objective Function

For a single-diode model, the objective function is expressed as(6)fV,I,X=Iph−IsdexpqV+Rs×In×K×T−1−V+Rs×IRsh,X=Iph,Isd,Rs,Rsh,n.For the double-diode model, the objective function is(7)fV,I,X=Iph−Isd1expqV+Rs×In1×K×T−1−Isd2expqV+Rs×In2×K×T−1−V+Rs×IRsh,X=Iph,Isd1,Isd2,Rs,Rsh,n1,n2.The parameters can be estimated by minimizing the objective functionRMSEX, i.e., by searching for the solution vector x [11, 16, 29, 65, 66].
(8)RMSEX=1N∑k=1NfV,I,X2,whereX represents the parameters composed of the solution vector. V and I are the measured voltage and current, respectively. N represents the number of experiments. Hence, to estimate the parameters is equivalent to search the X in the range which minimizes the objective function.


## 3. Problem Formulation

The problem can be set as an optimization problem with the objective to minimize the difference between the measured and estimated current. The objective function (OF) is defined as the root mean square error (RMSE), where the error function is defined as the difference between the estimated and experimental currents. It is expressed as follows [11, 16, 29, 65, 66]:
(9)MinRMSEX=1N∑k=1NIk,mes−Ik,extX2,whereRMSEX is the objective function to minimize, N is the number of points measured, Ik,mes is the measured current, and Ik,extX is the estimated current.For a single-diode model, the fitness function is expressed as(10)MinRMSEX=1N∑k=1NIk,mes−Iph+Isd1expqVk,mes+Rs×Ik,mesn1×K×T−1+Vk,mes+Rs×Ik,mesRsh2,withX=Iph,Isd,Rs,Rsh,n the parameters to be estimated.For the double-diode model, the fitness function is(11)MinRMSEX=1N∑k=1NIk,mes−Iph+Isd1expqVk,mes+Rs×Ik,mesn1×K×T−1+Isd2expqVk,mes+Rs×Ik,mesn2×K×T−1+Vk,mes+Rs×Ik,mesRsh2,withX=Iph,Isd1,Isd2,Rs,Rsh,n1,n2 the parameters to be estimated.In this paper, the approach of bald eagle search (BES) is used for the optimization of results of equations (10) and (11).


## 4. Bald Eagle Search (BES) Algorithm [62]

Bald eagles are occasional predators and are at the top of the food chain only because of their size. Furthermore, bald eagles are considered scavengers that feast on any available, easy, and protein-rich food. Bald eagles are an opportunistic forager that mainly select fish (alive or dead), especially salmon, as the primary food. Bald eagles frequently hunt from perch but may also hunt while in flight. They are capable of spotting fish at enormous distances because obtaining fish from water is difficult. When they start to search for food over a water spot, these eagles set off in a specific direction and select a certain area to begin the search. Accordingly, finding the search space is achieved by self-searching and tracking other birds with the concentration of fish (dead or alive).The proposed BES algorithm mimics the behaviour of bald eagles during hunting to justify the consequences of each hunting step. This algorithm is divided into three parts, namely, search space selection, search in the selected search space, and swooping.


### 4.1. Selection Stage

In the selection stage, bald eagles identify and select the best area (in terms of amount of food) within the selected search space where they can hunt for prey. Equation (12) presents this behaviour mathematically.
(12)Pi,new=Pbest+α∗Pmean−Pi∗r,whereα is the position change control parameter that takes a value between 1.5 and 2 and r is a random number that takes a value between 0 and 1; Pi,new and Pi are updated position and old position, respectively, at time i. In the selection step, the bald eagles select an area based on the information available from the previous step. The eagles randomly select another search area that differs from the previous search area but is located nearby. Pbest denotes the search area that is currently selected by the bald eagles based on the best position identified in their previous search. The eagles randomly search all points near the previously selected search space. Meanwhile, Pmean indicates that these eagles have used all the information from the previous points.


### 4.2. Search Stage

In the search stage, bald eagles search for prey within the selected search space and move in different directions within a spiral space to accelerate their search. The best position for the swoop is mathematically expressed in(13)Pi,new=Pi+yi∗Pi−Pi+1+xi∗Pi−Pmean,(14)xi=xrimaxxr,(15)yi=yrimaxyr,(16)xri=ri∗sinθi,(17)yri=ri∗cosθi,(18)θi=a∗π∗rand,(19)ri=θi+R∗rand,wherea is a parameter that takes a value between 5 and 10 for determining the corner point search in the central point and R takes a value between 0.5 and 2 for determining the number of search cycles.This algorithm uses the polar graph property to mathematically represent this movement. This property also allows the BES algorithm to discover new spaces and increase diversification by multiplying the difference between the current point and the next point with the polar point in they-axis and adding the difference between the current point and the center point with the polar point in the x-axis. We use the average solution in the search point because all search points move towards the centre point. All points in the polar plot take a value between -1 and 1, and we use a special equation for the shape of the spiral (20–22).


### 4.3. Swooping Stage

In the swooping stage, bald eagles swing from the best position in the search space to their target prey. All points also move towards the best point. Equation (14) mathematically illustrates this behaviour.
(20)Pi,new=rand∗Pbest+x1i∗Pi−c1∗Pmean+y1i∗Pi−c2∗Pbest,x1i=xrimaxxr,y1i=yrimaxyr,xri=ri∗sinhθi,yri=ri∗coshθi,θi=a∗π∗rand,ri=θi,wherec1,c2∈1,2.


### 4.4. Complete BES Algorithm

The previous steps have presented the main components of BES, which include the selection, search, and swooping steps. To describe the remaining operations and facilitate the implementation of BES, the flowchart algorithm is described in Figure4:Figure 4
Flowchart of BES algorithm.


## 4.1. Selection Stage

In the selection stage, bald eagles identify and select the best area (in terms of amount of food) within the selected search space where they can hunt for prey. Equation (12) presents this behaviour mathematically.
(12)Pi,new=Pbest+α∗Pmean−Pi∗r,whereα is the position change control parameter that takes a value between 1.5 and 2 and r is a random number that takes a value between 0 and 1; Pi,new and Pi are updated position and old position, respectively, at time i. In the selection step, the bald eagles select an area based on the information available from the previous step. The eagles randomly select another search area that differs from the previous search area but is located nearby. Pbest denotes the search area that is currently selected by the bald eagles based on the best position identified in their previous search. The eagles randomly search all points near the previously selected search space. Meanwhile, Pmean indicates that these eagles have used all the information from the previous points.


## 4.2. Search Stage

In the search stage, bald eagles search for prey within the selected search space and move in different directions within a spiral space to accelerate their search. The best position for the swoop is mathematically expressed in(13)Pi,new=Pi+yi∗Pi−Pi+1+xi∗Pi−Pmean,(14)xi=xrimaxxr,(15)yi=yrimaxyr,(16)xri=ri∗sinθi,(17)yri=ri∗cosθi,(18)θi=a∗π∗rand,(19)ri=θi+R∗rand,wherea is a parameter that takes a value between 5 and 10 for determining the corner point search in the central point and R takes a value between 0.5 and 2 for determining the number of search cycles.This algorithm uses the polar graph property to mathematically represent this movement. This property also allows the BES algorithm to discover new spaces and increase diversification by multiplying the difference between the current point and the next point with the polar point in they-axis and adding the difference between the current point and the center point with the polar point in the x-axis. We use the average solution in the search point because all search points move towards the centre point. All points in the polar plot take a value between -1 and 1, and we use a special equation for the shape of the spiral (20–22).


## 4.3. Swooping Stage

In the swooping stage, bald eagles swing from the best position in the search space to their target prey. All points also move towards the best point. Equation (14) mathematically illustrates this behaviour.
(20)Pi,new=rand∗Pbest+x1i∗Pi−c1∗Pmean+y1i∗Pi−c2∗Pbest,x1i=xrimaxxr,y1i=yrimaxyr,xri=ri∗sinhθi,yri=ri∗coshθi,θi=a∗π∗rand,ri=θi,wherec1,c2∈1,2.


## 4.4. Complete BES Algorithm

The previous steps have presented the main components of BES, which include the selection, search, and swooping steps. To describe the remaining operations and facilitate the implementation of BES, the flowchart algorithm is described in Figure4:Figure 4
Flowchart of BES algorithm.


## 5. Simulation Results and Analysis

The BES algorithm is applied to extract the SDM, DDM, and PV module parameters. To examine in more detail the accuracy of the data obtained by the BES method for the optimized parameters, the current was calculated from the values estimated on the basis of the different models and compared with that obtained from the experimental measurements. The error in the measured values for each of the models was evaluated by IAE (individual absolute error) and RE (relative error), calculated as shown in equations (21) and (22), respectively.
(21)IAE=Imeasured−Iestimated,(22)RE=Imeasured−IestimatedImeasured.The lower and upper bounds are expressed in Table1 [18, 63, 67].Table 1
Parameter range of different PV models.
ParametersSingle and double diodePhotowatt-PWP201STM6-40/36STP6-120/36Lower boundUpper boundLower boundUpper boundLower boundUpper boundLower boundUpper boundIph (A)01020208Id,Id1,Id2 (μA)01050050050Rs (Ω)00.50200.3600.36Rsh (Ω)0100020000100001500n,n1,n212150160150


### 5.1. Case Study 1: Single-Diode Model

In this case study, the BES optimization algorithm is used to extract the parameters of the two proposed models of the R.T.C. France solar cell. The measured data of the characteristic curves (I-V) of the R.T.C. France solar cell are reported in [63, 67]. Table 2 includes the results of the parameters estimated based on BES and those estimated based on other optimization techniques such as ILSA [68], GAMS [3], ITLBO [69], IMFO [67], IJAYA [63], MADE [18], EVPS [52], GOTLBO [64], DDSFLA [57], EGBO [56], and CLPSO [70]. From Table 2, it can be seen that for the SDM, the application of the proposed BES algorithm results in the minimum RMSE value which is equal to 9.8602e−04.Table 2
Comparison of the results obtained from single diode model R.T.C. France solar cell with other methods in the literature.
AlgorithmIph (A)Isd (μ A)Rs (Ω)Rsh (Ω)nRMSEBES0.76070.32300.036453.71851.48129.8602e−04ILSA [68]0.76070.32290.036453.72041.48119.8602e−04GAMS [3]0.76070.32300.036453.71851.48129.8602e−04ITLBO [69]0.76070.32300.036453.71841.48129.8602e−04IMFO [67]0.76070.32340.036353.76081.48139.8602e−04IJAYA [63]0.76080.32280.036453.75951.48119.8603e−04MADE [18]0.76070.32300.036353.71851.48119.8602e−04EVPS [52]0.76070.32500.036353.89601.48219.8609e−04GOTLBO [64]0.76080.32970.036353.36641.48339.8856e−04TLABC [20]0.76080.32300.036453.71641.48129.8602e−04CLPSO [70]0.76080.34300.036154.19651.48739.9633e−04TLBO [71]0.76070.32940.036354.30151.48319.8733e−03TPTLBO [72]0.76080.32300.036453.71851.48129.8602e−04DSCE [55]0.76070.32300.036353.71851.48119.8602e−04EGBO [56]0.76080.32300.036453.71851.48119.8602e−04DDSFLA [57]0.76080.31910.03653.37701.48009.8630e−04SFLBS [58]0.76070.32300.036353.71851.48119.8602e−04LCNMSE [59]0.76070.32300.036353.711.48119.8602e−04PSOCS [60]0.76070.32300.036353.7191.48129.8602e−04CNMSMA [61]0.76070.32300.036353.71821.48119.8602e−04From Figure5, it can be clearly seen that the I-V and P-V curves of the simulated data found by BES are very compatible with the experimental data.Figure 5
Characteristics of the measured and estimated curve with one diode: (a)I-V and (b) P-V.
(a)
I-V characteristic curve(b)
P-V characteristic curveFrom Table3, it is evident that all IAE values are lower than 3.399e−03 and the RE values are between −6.91e−03 and 2.65897e−01, demonstrating the high efficiency identified by BES for the single-diode model.Table 3
The individual absolute error and relative error of a single-diode model obtained by BES.
ItemTest dataSimulated currentPowerV (V)I (A)Isim (A)IAE (A)Psim (W)RE1−0.20570.76400.7641120.000112-0.157178-0.0001462−0.12910.76200.7626870.000687-0.098463-0.0009023−0.05880.76050.7613790.000879-0.044769-0.00115640.00570.76050.7601780.0003220.0043330.00042350.06460.76000.7590790.0009210.0490370.00121160.11850.75900.7580660.0009340.0898310.00123070.16780.75700.7571160.0001160.127044-0.00015380.21320.75700.7561660.0008340.1612140.00110290.25450.75550.7551110.0003890.1921760.000515100.29240.75400.7536890.0003110.2203790.000413110.32690.75050.7514170.0009170.245638-0.001222120.35850.74650.7473830.0008830.267937-0.001183130.38730.73850.7401520.0016520.286661-0.002238140.41370.72800.7274300.0005700.3009380.000783150.43730.70650.7070430.0005430.309190-0.000769160.45900.67550.6753890.0001110.3100030.000165170.47840.63200.6309250.0010750.3018350.001700180.49600.57300.5721800.0008200.2838010.001432190.51190.49900.4999710.0009710.255935-0.001945200.52650.41300.4141570.0011570.218054-0.002802210.53980.31650.3181940.0016940.171761-0.005352220.55210.21200.2130460.0010460.117623-0.004934230.56330.10350.1033750.0001250.0582310.001209240.5736−0.0100-0.0073410.002659-0.0042110.265897250.5833−0.1230-0.1238500.000850-0.072242-0.006910260.5900−0.2100-0.2066010.003399-0.1218950.016186Total IAE0.023977Table4 and Figure 6 show the statistical results and the convergence curve, respectively. From Table 4, it can be seen that BES obtains the best minimum value of RMSE.Table 4
Statistical results for R.T.C. France PV cell single diode.
AlgorithmRMSEMinMeanMaxStd.BES9.8602e−049.8602e−049.8602e−042.6314e−13ITLBO9.8602e−049.8602e−049.8602e−42.19e−17IMFO9.8602e−049.8767e−049.9641e−042.1810e−06IJAYA9.8606e−041.0261e−031.1223e−034.160e−05GOTLBO9.8602e−041.4388e−031.0289e−031.01e−04CLPSO9.9455e−041.0507e−031.1865e−044.6730e−05TPTLBO9.8602e−049.8602e−049.8602e−042.28e−17MADE9.8602e−049.8602e−049.8602e−042.47e−15TLABC9.8602e−049.9852e−041.2358e−031.86e−05TLBO9.8722e−041.0476e−041.0397e−036.59e−05DSCE9.8602e−049.8602e−049.8602e−041.1320e−09EGBO9.8602e−049.9500e−041.1161e−042.62e−05DDSFLA9.8630e−041.0819e−031.3056e−038.6464e−05SFLBS9.8602e−049.8602e−049.8602e−041.4301e−14PSOCS9.8602e−049.8602e−049.8603e−041.7459e−09Figure 6
Convergence curve during the parameter extraction for one diode.


### 5.2. Case Study 2: Two-Diode Model

Table5 lists the results obtained from the application of the BES technique to extract the DDM parameters of the R.T.C. France solar cell. To validate the applied technique, the table also presents the results from the application of other techniques of GAMS [3], ITLBO [69], IMFO [67], IJAYA [63], MADE [18], EVPS [52], GOTLBO [64], and CLPSO [70]. The table shows that the BES optimization technique applied gives the best results with the minimum objective function of RMSE being 0.000982484851801148.Table 5
Comparison of the results obtained from the double-diode model R.T.C. France solar cell with other methods in the literature.
AlgorithmIph (A)Isd1 (μA)Isd2 (μA)Rs (Ω)Rsh (Ω)n1n2Best RMSEBES0.76080.22590.74930.036755.48541.45102.0009.8248e−4GAMS [3]0.76080.22600.75940.036755.48541.21861.62479.8293e−4ITLBO [69]0.76080.22600.74930.036755.48541.45102.00009.8248e−4IMFO [67]0.76070.23500.68370.036755.29971.45372.00009.8252e−4IJAYA [63]0.76010.00500.75090.037677.85191.21861.62479.8293e−4MADE [18]0.76080.22460.73940.036855.43291.45051.99639.8261e−4EVPS [52]0.76070.29750.25040.036355.88271.47491.97269.8510e−4GOTLBO [64]0.76080.27170.25950.036653.61871.46681.91619.9544e−4CLPSO [70]0.76060.28750.26860.036655.28951.95861.46529.9224e−4TLABC [20]0.76080.42390.24010.036754.66801.90751.45679.8415e−4TLBO [71]0.76100.29470.13730.036653.12101.47301.99381.0069e−03TPTLBO [72]0.76080.74340.22660.036755.48312.00001.45139.8248e−04DSCE [55]0.76080.69800.23180.036755.37501.99991.45539.8250e−04EGBO [56]0.76080.2250.7490.036755.48551.45102.00009.8248e−04DDSFLA [57]0.76080.29310.22710.036554.37101.47302.00009.8434e−04SFLBS [58]0.76070.77590.22850.036755.54962.00001.44989.8249e−04LCNMSE [59]0.76070.74930.22590.036755.48542.00001.45109.8248e−04PSOCS [60]0.76071.00000.19810.036856.1722.00001.44019.8297e−04CNMSMA [61]0.76070.22590.75060.036755.48541.45101.99999.8249e−04From Figure7, it can be clearly seen that the I-V and P-V curves of the simulated data found by BES are very compatible with the experimental data.Figure 7
Characteristics of the measured and estimated curve with two diodes: (a)P-V and (b) I-V.
(a)
P-V characteristic(b)
I-V characteristicFrom Table6, it is evident that all IAE values are lower than 2.836e−03 and the RE values are between −1.1111e−02 and 2.3090e−01, demonstrating the high efficiency identified by BES for the two-diode model.Table 6
The individual absolute error and relative error of double-diode model obtained by BES.
ItemTest dataSimulated currentPowerV (V)I (A)Isim (A)IAE (A)Psim (W)RE1−0.20570.76400.7640030.000003-0.157155-0.0000042−0.12910.76200.7626240.000624-0.098455-0.0008183−0.05880.76050.7613570.000857-0.044768-0.00112740.00570.76050.7601930.0003070.0043330.00040350.06460.76000.7591270.0008730.0490400.00114860.11850.75900.7581410.0008590.0898400.00113270.16780.75700.7572080.0002080.127060-0.00027580.21320.75700.7562630.0007370.1612350.00097390.25450.75550.7551970.0003030.1921980.000400100.29240.75400.7537440.0002560.2203950.000340110.32690.75050.7514230.0009230.245640-0.001230120.35850.74650.7473320.0008320.267919-0.001114130.38730.73850.7400540.0015540.286623-0.002104140.41370.72800.7273130.0006870.3008890.000943150.43730.70650.7069550.0004550.309151-0.000644160.45900.67550.6753760.0001240.3099970.000184170.47840.63200.6310100.0009900.3018750.001567180.49600.57300.5723510.0006490.2838860.001132190.51190.49900.5001880.0011880.256046-0.002381200.52650.41300.4143530.0013530.218157-0.003277210.53980.31650.3183050.0018050.171821-0.005702220.55210.21200.2130150.0010150.117605-0.004787230.56330.10350.1031720.0003280.0581170.003166240.5736−0.0100-0.0076910.002309-0.0044120.230909250.5833−0.1230-0.1243670.001367-0.072543-0.011111260.5900−0.2100-0.2071640.002836-0.1222260.013507Total IAE0.023441Table7 and Figure 8 show the statistical results and the convergence curve, respectively. From Table 7, it can be seen that BES obtains the best minimum value of RMSE.Table 7
Statistical results for R.T.C. France PV cell two diodes.
AlgorithmRMSEMinMeanMaxStd.BES9.8248e−049.9518e−041.1881e−034.6013e−05ITLBO9.8248e−049.8812e−049.8497e−041.54e−06IMFO9.8252e−049.9737e−041.1409e−033.2939e−05IJAYA9.8380e−041.0240e−031.3507e−038.5647e−05GOTLBO9.8407e−041.4380e−031.0453e−031.01e−04CLPSO9.9224e−041.0522e−031.1462e−034.3141e−05TPTLBO9.8248e−049.8602e−049.8363e−049.31e−07MADE9.8261e−049.8608e−049.8786e−048.02e−05TLABC9.8415e−041.0555e−031.5048e−031.54e−06TLBO1.0069e−031.1598e−031.5206e−031.56e−04EGBO9.8248e−049.8484e−049.8681e−041.66e−04DDSFLA9.8434e−041.1071e−031.4225e−031.3014e−04SFLBS9.8249e−049.8541e−049.8787e−041.7882e−06PSCOS9.8297e−041.0286e−031.4133e−039.9217e−05Figure 8
Convergence curve during the parameter extraction for the two-diode models.


### 5.3. Case Study 3: Photowatt-PWP201 PV Module

To further validate the BES technique and show its effectiveness in estimating the optimal parameters of different models, we used this algorithm on the Photowatt-PWP201 PV module, which consists of 36 silicon cells connected in series under operating conditions of 1000 W/m2 of solar irradiation and a cell temperature of 45°C. The results obtained were compared with those reported in the literature based on other techniques.The results have been listed in Table8. This table also presents a comparison with the results of other techniques from the literature ITLBO [69], IMFO [67], IJAYA [63], and MADE [18]. The comparison validated the effectiveness of BES compared to other techniques. The RMSE based on BES application to extract the parameters of PV model is equal to 0.00242507 which is better. From Figure 9, it can be clearly seen that the I-V and P-V curves of the simulated data found by BES are very compatible with the experimental data.Table 8
Comparison of the results obtained from Photowatt-PWP201 model with other methods in the literature.
AlgorithmIph (A)Isd (μ A)Rs (Ω)Rsh (Ω)nRMSEBES1.03053.48231.2013981.982448.64282.42507e−03GAMS [3]1.03203.26811.2062828.29281.34452.4426e−03ITLBO [69]1.03053.48231.2013981.982348.64282.4251e−03IMFO [67]1.03053.47831.2013980.467248.63852.4251e−03IJAYA [63]1.03023.47031.2011984.876048.64822.4251e−03MADE [18]1.03053.48231.2013981.982348.64282.4251e−03EVPS [52]1.03183.26791.2066845.7591.34452.4267e−03GOTLBO [64]1.03053.49911.2008989.688948.66112.4251e−03TLABC [20]1.03063.47151.2017972.935748.63132.4251e−03TLBO [71]1.03053.48721.2011984.876048.64822.4251e−03CLPSO [70]1.03043.61311.19781017.048.78472.4280e−03TPTLBO [72]1.03053.48231.2013981.982248.64282.4251e−03EGBO [56]1.03053.481.2013981.982248.64282.4151e−03DDSFLA [57]1.03063.44731.2023971.250048.60402.4252e−03SFLBS [58]1.03053.48221.2012981.980448.64282.4251e−03LCNMSE [59]1.03153.48221.2013981.974148.64282.4251e−03PSOCS [60]1.03053.48231.2013981.9848.6432.4251e−03Figure 9
Characteristics of the measured and estimated curve with one diode: (a)I-V and (b) P-V.
(a)
I-V characteristic curve(b)
P-V characteristic curveFrom Table9, it is evident that all IAE values are lower than 4.418e−03 and the RE values are between −6.44e−03 and 1.3167e−01, demonstrating the high efficiency identified by BES for the PV model.Table 9
The individual absolute error and relative error of Photowatt-PWP201 obtained by BES.
ItemTest dataSimulated currentPowerV (V)I (A)Isim (A)IAE (A)Psim (W)RE10.12481.03151.0291050.0023950.1284320.00232221.80931.03001.0273670.0026331.8588150.00255733.35111.02601.0257280.0002723.4373150.00026644.76221.02201.0240930.0020934.876936-0.00204856.05381.01801.0222780.0042786.188666-0.00420267.23641.01551.0199180.0044187.380532-0.00435078.31891.01401.0163520.0023528.454931-0.00232089.30971.01001.0104900.0004909.407354-0.000485910.21631.00351.0006310.00286910.2227480.0028591011.04490.98800.9845670.00343310.8744390.0034751111.80180.96300.9595670.00343311.3246160.0035651212.49290.92550.9229260.00257411.5300230.0027811313.12310.87250.8727480.00024811.453153-0.0002841413.69830.80750.8075040.00000411.061434-0.0000051514.22210.72650.7286690.00216910.363202-0.0029851614.69950.63450.6375920.0030929.372288-0.0048741715.13460.53450.5368060.0023068.124337-0.0043131815.53110.42750.4302530.0027536.682303-0.0064401915.89290.31850.3196740.0011745.080544-0.0036862016.22290.20850.2084500.0000503.3816630.0002402116.52410.10100.0973910.0036091.6092940.0357362216.7987−0.0080-0.0069470.001053-0.1166940.1316702317.0499−0.1110-0.1094040.001596-1.8653230.0143802417.2793−0.2090-0.2075660.001434-3.5865890.0068632517.4885−0.3030-0.2990420.003958-5.2297960.013063Total IAE0.054687Table10 and Figure 10 show the statistical results and the convergence curve, respectively. From Table 10, it can be seen BES obtains the best minimum value of RMSE.Table 10
The statistical results for Photowatt-PWP201.
AlgorithmRMSEMinMeanMaxStd.BES2.42507e−32.42507e−32.42507e−32.4532e−17ITLBO2.4251e−032.4251e−032.4251e−031.27e−17IMFO2.4251e−032.4294e−032.5005e−031.3831e−05IJAYA2.4251e−032.4393e−032.4289e−033.78e−06GOTLBO2.4251e−032.4852e−032.4419e−031.38e−05CLPSO2.4280e−032.4549e−032.5432e−032.5809e−05TPTLBO2.4251e−032.4251e−032.4251e−031.20e−17MADE2.4250e−032.4251e−032.4251e−031.96e−16TLABC2.4251e−032.4265e−032.4458e−034.00e−06TLBO2.4251e−032.4383e−032.5475e−032.43e−05EGBO2.4251e−032.4251e−032.4251e−032.38e−17DDSFLA2.4252e−032.4974e−032.9281e−031.1000e−04SFLBS2.4251e−032.4251e−032.4251e−033.9417e−17PSOCS2.4251e−032.4252e−032.4282e−035.9113e−07Figure 10
Convergence curve during the parameter extraction for the two-diode models.


### 5.4. Case Study 4: Schutten Solar STM6-40/36 Monocrystalline PV Module

Here, we use the BES algorithm to extract the parameters of the Schutten Solar STM6-40/36 PV module. It contains 36 polycrystalline cells (size156mm×156mm) connected in series. The data set contains 20 data points measured at T=51°C [66]. For the STM6-40/36 PV module model, Table 11 shows the results of the parameters obtained from ITLBO, IJAYA, GOTLBO, TPTLBO …. From the results, it can be seen that the BES provides a better RMSE: 0.00172981370994066.Table 11
Comparison of the results obtained from Solar STM6-40/36 PV model with other methods in the literature.
AlgorithmIph (A)Isd (μA)Rs (Ω)Rsh (Ω)nRMSEBES1.66391.73870.004315.92831.52031.7298e−03ITLBO [69]1.66391.73870.004315.92831.52031.7298e−03IJAYA [63]1.66371.83530.004015.94491.52631.7548e−03MADE [18]1.66391.13870.004315.92831.52031.7298e−03GOTLBO [64]1.66391.73870.004315.92831.52031.7298e−03TPTLBO [72]1.66391.73870.004315.92831.52031.7298e−03TLABC [20]1.66391.73870.004315.92831.52031.7298e−03TLBO [71]1.66381.73070.004315.99551.51981.7305e−03EGBO [56]1.66391.730.004315.92831.52031.7298e−03In addition, to confirm the accuracy of the extracted parameters, Figure11 shows the I-V and P-V curves. It is evident that the simulated data from the BES match well with the measured data in the voltage range for both I-V and P-V curves.Figure 11
Characteristics of the measured and estimated curve Schutten Solar STM6-40/36 PV model: (a)I-V and (b) P-V.
(a)
I-V characteristic curve(b)
P-V characteristic curveIn addition, the IAE and RE (relative error) are given in Table12. The IAE describes the error between the extracted parameter and the measured data. In other words, the extracted parameters are better when the IAE is small. According to Table 12, the sum of the IAE is less than 2.50E−02, which indicates that the measured and extracted data coincide well.Table 12
The individual absolute error and relative error of Schutten Solar STM6-40/36 PV model obtained by BES.
ItemTest dataSimulated currentPowerV (V)I (A)Isim (A)IAE (A)Psim (W)RE10.0001.6631.6634530.0004530.000000-0.00027320.1181.6631.6632470.0002470.196263-0.00014932.2371.6611.6595460.0014543.7124040.00087545.4341.6531.6539100.0009108.987346-0.00055057.261.651.6505610.00056111.983074-0.00034069.681.6451.6454260.00042615.927728-0.000259711.591.641.6392320.00076818.9986960.000468812.61.6361.6337140.00228620.5848020.001397913.371.6291.6272930.00170721.7569020.0010481014.091.6191.6183280.00067222.8022470.0004151114.881.5971.6031200.00612023.854427-0.0038321215.591.5811.5816420.00064224.657796-0.0004061316.41.5421.5424300.00043025.295849-0.0002791416.711.5241.5213180.00268225.4212170.0017601516.981.51.4993470.00065325.4589140.0004351617.131.4851.4854450.00044525.445677-0.0003001717.321.4651.4658490.00084925.388510-0.0005801817.911.3881.3878870.00011324.8570560.0000811919.081.1181.1190650.00106521.351757-0.0009522022.0200.0023860.0023860.050143—Total IAE0.02486996To further prove the reliability of the BES, the statistical results containing the minimum (Min), maximum (Max), mean value (Mean), and standard deviation (Std.) are analyzed. Table13 shows the statistical results, and Figure 12 shows the convergence curve of Photowatt STM6-40/36.Table 13
The statistical results for Photowatt STM6-40/36.
AlgorithmRMSEMinMeanMaxStd.BES1.7298e−031.7298e−031.7298e−035.6525e−18ITLBO1.7298e−031.7298e−031.7298e−037.13e−03IJAYA17548e−032.5223e−031.9305e−031.91e−04GOTLBO1.7298e−031.1244e−034.2347e−036.41e−02TPTLBO1.7298e−031.7298e−031.7298e−034.96e−18MADE1.7298e−031.7298e−031.7298e−038.49e−14TLABC1.7298e−032.1827e−036.5053e−039.22e−04TLBO1.7305e−034.3487e−032.0593e−023.45e−03EGBO1.7298e−031.7298e−031.7298e−038.22e−18Figure 12
Convergence curve during the parameter extraction for Solar STM6-40/36 PV model.


### 5.5. Case Study 5: STP6-120/36 Module

The polycrystalline STP6-120/36 has 36 cells connected in series and is measured under 1000 W/m2 at 55°C. The current-voltage data was obtained from [72]. The results have been listed in Table 14. This table also presents a comparison with the results of other techniques from the literature ITLBO [69], IMFO [67], IJAYA [63], TLABC [20], GOTLBO [64], and MADE [18]. The comparison validated the effectiveness of BES compared to other techniques. The RMSE based on BES application to extract the parameters of PV model equals 0.0167828544285882. From Figure 13, it can be clearly seen that the I-V and P-V curves of the simulated data found by BES are very compatible with the experimental data.Table 14
Comparison of the results obtained from Solar STP6-120/36 PV model with other methods in the literature.
AlgorithmIph (A)Isd (μ A)Rs (Ω)Rsh (Ω)nRMSEBES7.47162.32180.004623.02651.25961.6782e−03ITLBO [69]7.47252.33500.004622.21991.26011.6601e−02IJAYA [63]7.46722.25360.004627.59251.25711.6731e−02MADE [18]7.47252.33500.004622.21991.26011.6601e−02GOTLBO [64]7.47252.33500.004622.21991.26011.6601e−02TPTLBO [72]7.47573.01000.004622.21991.26013.0343e−01TLABC [20]7.47252.33490.004622.21171.26011.6601e−02TLBO [71]7.47821.91940.004713.26881.24401.6892e−02EGBO [56]7.47252.330.004622.21991.26011.6601e−02Figure 13
Characteristics of the measured and estimated curve SPT6-120/36: (a)I-V and (b) P-V.
(a)
I-V characteristic curve(b)
P-V characteristic curveIn addition, the IAE and RE (relative error) are given in Table15. The IAE describes the error between the extracted parameter and the measured data. In other words, the extracted parameters are better when the IAE is small. According to Table 15, the sum of the IAE is less than 4.64E−02, which indicates that the measured and extracted data coincide well.Table 15
The individual absolute error and relative error of Solar STP6-120/36 PV model obtained by BES.
ItemTest dataSimulated currentPowerV (V)I (A)Isim (A)IAE (A)Psim (W)RE107.487.4701030.00989791.0073840.00132329.067.457.4520810.00208193.400729-0.00027939.477.427.4489540.02895495.858207-0.003902410.327.447.4387160.00128498.0869360.000173511.177.417.4200200.010020100.437238-0.001352611.817.387.3957830.015783101.453230-0.002139712.367.377.3630570.006943101.7408930.000942812.747.347.3312970.008703101.0794300.001186913.167.297.2840580.00594299.7392770.0008151013.597.237.2175820.01241897.2267420.0017181114.177.17.0880200.01198094.5287000.0016871214.586.976.9583830.01161788.4607380.0016671314.936.836.8145270.01547386.0062300.0022651415.396.586.5678640.01213681.8815400.0018441515.716.366.3487760.01122478.4099650.0017651616.0866.0464390.04643974.377396-0.0077401716.345.755.7851100.03511067.754636-0.0061061816.765.275.2780870.0080870.274440-0.0015351916.95.075.0891260.01912691.007384-0.0037722017.14.794.7883940.00160693.4007290.0003352117.254.564.5455050.01449595.8582070.0031792217.414.294.2721080.01789298.0869360.0041712317.653.833.8387900.008790100.437238-0.0022952419.2100.0142860.014286101.453230—Total IAE0.330286Table16 and Figure 14 show the statistical results and the convergence curve, respectively.Table 16
The statistical results for Photowatt STP6-120/36 obtained by BES.
AlgorithmRMSEMinMeanMaxStd.BES1.6782e−021.69034e−021.7223e−021.1219e−04ITLBO1.6601e−021.6601e−021.6601e−027.22e−17IJAYA1.6731e−021.6891e−021.7304e−021.12e−04GOTLBO1.6601e−022.9588e−021.8099e−013.05e−02MADE1.6601e−021.6601e−021.6601e−021.69e−15TLABC1.6601e−021.6963e−022.1497e−029.47e−04TLBO1.6892e−023.6690e−022.1604e−023.51e−02EGBO1.6601e−021.6601e−041.6601e−021.47e−16Figure 14
Convergence curve during the parameter extraction for Solar STP6-120/36 PV model.


## 5.1. Case Study 1: Single-Diode Model

In this case study, the BES optimization algorithm is used to extract the parameters of the two proposed models of the R.T.C. France solar cell. The measured data of the characteristic curves (I-V) of the R.T.C. France solar cell are reported in [63, 67]. Table 2 includes the results of the parameters estimated based on BES and those estimated based on other optimization techniques such as ILSA [68], GAMS [3], ITLBO [69], IMFO [67], IJAYA [63], MADE [18], EVPS [52], GOTLBO [64], DDSFLA [57], EGBO [56], and CLPSO [70]. From Table 2, it can be seen that for the SDM, the application of the proposed BES algorithm results in the minimum RMSE value which is equal to 9.8602e−04.Table 2
Comparison of the results obtained from single diode model R.T.C. France solar cell with other methods in the literature.
AlgorithmIph (A)Isd (μ A)Rs (Ω)Rsh (Ω)nRMSEBES0.76070.32300.036453.71851.48129.8602e−04ILSA [68]0.76070.32290.036453.72041.48119.8602e−04GAMS [3]0.76070.32300.036453.71851.48129.8602e−04ITLBO [69]0.76070.32300.036453.71841.48129.8602e−04IMFO [67]0.76070.32340.036353.76081.48139.8602e−04IJAYA [63]0.76080.32280.036453.75951.48119.8603e−04MADE [18]0.76070.32300.036353.71851.48119.8602e−04EVPS [52]0.76070.32500.036353.89601.48219.8609e−04GOTLBO [64]0.76080.32970.036353.36641.48339.8856e−04TLABC [20]0.76080.32300.036453.71641.48129.8602e−04CLPSO [70]0.76080.34300.036154.19651.48739.9633e−04TLBO [71]0.76070.32940.036354.30151.48319.8733e−03TPTLBO [72]0.76080.32300.036453.71851.48129.8602e−04DSCE [55]0.76070.32300.036353.71851.48119.8602e−04EGBO [56]0.76080.32300.036453.71851.48119.8602e−04DDSFLA [57]0.76080.31910.03653.37701.48009.8630e−04SFLBS [58]0.76070.32300.036353.71851.48119.8602e−04LCNMSE [59]0.76070.32300.036353.711.48119.8602e−04PSOCS [60]0.76070.32300.036353.7191.48129.8602e−04CNMSMA [61]0.76070.32300.036353.71821.48119.8602e−04From Figure5, it can be clearly seen that the I-V and P-V curves of the simulated data found by BES are very compatible with the experimental data.Figure 5
Characteristics of the measured and estimated curve with one diode: (a)I-V and (b) P-V.
(a)
I-V characteristic curve(b)
P-V characteristic curveFrom Table3, it is evident that all IAE values are lower than 3.399e−03 and the RE values are between −6.91e−03 and 2.65897e−01, demonstrating the high efficiency identified by BES for the single-diode model.Table 3
The individual absolute error and relative error of a single-diode model obtained by BES.
ItemTest dataSimulated currentPowerV (V)I (A)Isim (A)IAE (A)Psim (W)RE1−0.20570.76400.7641120.000112-0.157178-0.0001462−0.12910.76200.7626870.000687-0.098463-0.0009023−0.05880.76050.7613790.000879-0.044769-0.00115640.00570.76050.7601780.0003220.0043330.00042350.06460.76000.7590790.0009210.0490370.00121160.11850.75900.7580660.0009340.0898310.00123070.16780.75700.7571160.0001160.127044-0.00015380.21320.75700.7561660.0008340.1612140.00110290.25450.75550.7551110.0003890.1921760.000515100.29240.75400.7536890.0003110.2203790.000413110.32690.75050.7514170.0009170.245638-0.001222120.35850.74650.7473830.0008830.267937-0.001183130.38730.73850.7401520.0016520.286661-0.002238140.41370.72800.7274300.0005700.3009380.000783150.43730.70650.7070430.0005430.309190-0.000769160.45900.67550.6753890.0001110.3100030.000165170.47840.63200.6309250.0010750.3018350.001700180.49600.57300.5721800.0008200.2838010.001432190.51190.49900.4999710.0009710.255935-0.001945200.52650.41300.4141570.0011570.218054-0.002802210.53980.31650.3181940.0016940.171761-0.005352220.55210.21200.2130460.0010460.117623-0.004934230.56330.10350.1033750.0001250.0582310.001209240.5736−0.0100-0.0073410.002659-0.0042110.265897250.5833−0.1230-0.1238500.000850-0.072242-0.006910260.5900−0.2100-0.2066010.003399-0.1218950.016186Total IAE0.023977Table4 and Figure 6 show the statistical results and the convergence curve, respectively. From Table 4, it can be seen that BES obtains the best minimum value of RMSE.Table 4
Statistical results for R.T.C. France PV cell single diode.
AlgorithmRMSEMinMeanMaxStd.BES9.8602e−049.8602e−049.8602e−042.6314e−13ITLBO9.8602e−049.8602e−049.8602e−42.19e−17IMFO9.8602e−049.8767e−049.9641e−042.1810e−06IJAYA9.8606e−041.0261e−031.1223e−034.160e−05GOTLBO9.8602e−041.4388e−031.0289e−031.01e−04CLPSO9.9455e−041.0507e−031.1865e−044.6730e−05TPTLBO9.8602e−049.8602e−049.8602e−042.28e−17MADE9.8602e−049.8602e−049.8602e−042.47e−15TLABC9.8602e−049.9852e−041.2358e−031.86e−05TLBO9.8722e−041.0476e−041.0397e−036.59e−05DSCE9.8602e−049.8602e−049.8602e−041.1320e−09EGBO9.8602e−049.9500e−041.1161e−042.62e−05DDSFLA9.8630e−041.0819e−031.3056e−038.6464e−05SFLBS9.8602e−049.8602e−049.8602e−041.4301e−14PSOCS9.8602e−049.8602e−049.8603e−041.7459e−09Figure 6
Convergence curve during the parameter extraction for one diode.


## 5.2. Case Study 2: Two-Diode Model

Table5 lists the results obtained from the application of the BES technique to extract the DDM parameters of the R.T.C. France solar cell. To validate the applied technique, the table also presents the results from the application of other techniques of GAMS [3], ITLBO [69], IMFO [67], IJAYA [63], MADE [18], EVPS [52], GOTLBO [64], and CLPSO [70]. The table shows that the BES optimization technique applied gives the best results with the minimum objective function of RMSE being 0.000982484851801148.Table 5
Comparison of the results obtained from the double-diode model R.T.C. France solar cell with other methods in the literature.
AlgorithmIph (A)Isd1 (μA)Isd2 (μA)Rs (Ω)Rsh (Ω)n1n2Best RMSEBES0.76080.22590.74930.036755.48541.45102.0009.8248e−4GAMS [3]0.76080.22600.75940.036755.48541.21861.62479.8293e−4ITLBO [69]0.76080.22600.74930.036755.48541.45102.00009.8248e−4IMFO [67]0.76070.23500.68370.036755.29971.45372.00009.8252e−4IJAYA [63]0.76010.00500.75090.037677.85191.21861.62479.8293e−4MADE [18]0.76080.22460.73940.036855.43291.45051.99639.8261e−4EVPS [52]0.76070.29750.25040.036355.88271.47491.97269.8510e−4GOTLBO [64]0.76080.27170.25950.036653.61871.46681.91619.9544e−4CLPSO [70]0.76060.28750.26860.036655.28951.95861.46529.9224e−4TLABC [20]0.76080.42390.24010.036754.66801.90751.45679.8415e−4TLBO [71]0.76100.29470.13730.036653.12101.47301.99381.0069e−03TPTLBO [72]0.76080.74340.22660.036755.48312.00001.45139.8248e−04DSCE [55]0.76080.69800.23180.036755.37501.99991.45539.8250e−04EGBO [56]0.76080.2250.7490.036755.48551.45102.00009.8248e−04DDSFLA [57]0.76080.29310.22710.036554.37101.47302.00009.8434e−04SFLBS [58]0.76070.77590.22850.036755.54962.00001.44989.8249e−04LCNMSE [59]0.76070.74930.22590.036755.48542.00001.45109.8248e−04PSOCS [60]0.76071.00000.19810.036856.1722.00001.44019.8297e−04CNMSMA [61]0.76070.22590.75060.036755.48541.45101.99999.8249e−04From Figure7, it can be clearly seen that the I-V and P-V curves of the simulated data found by BES are very compatible with the experimental data.Figure 7
Characteristics of the measured and estimated curve with two diodes: (a)P-V and (b) I-V.
(a)
P-V characteristic(b)
I-V characteristicFrom Table6, it is evident that all IAE values are lower than 2.836e−03 and the RE values are between −1.1111e−02 and 2.3090e−01, demonstrating the high efficiency identified by BES for the two-diode model.Table 6
The individual absolute error and relative error of double-diode model obtained by BES.
ItemTest dataSimulated currentPowerV (V)I (A)Isim (A)IAE (A)Psim (W)RE1−0.20570.76400.7640030.000003-0.157155-0.0000042−0.12910.76200.7626240.000624-0.098455-0.0008183−0.05880.76050.7613570.000857-0.044768-0.00112740.00570.76050.7601930.0003070.0043330.00040350.06460.76000.7591270.0008730.0490400.00114860.11850.75900.7581410.0008590.0898400.00113270.16780.75700.7572080.0002080.127060-0.00027580.21320.75700.7562630.0007370.1612350.00097390.25450.75550.7551970.0003030.1921980.000400100.29240.75400.7537440.0002560.2203950.000340110.32690.75050.7514230.0009230.245640-0.001230120.35850.74650.7473320.0008320.267919-0.001114130.38730.73850.7400540.0015540.286623-0.002104140.41370.72800.7273130.0006870.3008890.000943150.43730.70650.7069550.0004550.309151-0.000644160.45900.67550.6753760.0001240.3099970.000184170.47840.63200.6310100.0009900.3018750.001567180.49600.57300.5723510.0006490.2838860.001132190.51190.49900.5001880.0011880.256046-0.002381200.52650.41300.4143530.0013530.218157-0.003277210.53980.31650.3183050.0018050.171821-0.005702220.55210.21200.2130150.0010150.117605-0.004787230.56330.10350.1031720.0003280.0581170.003166240.5736−0.0100-0.0076910.002309-0.0044120.230909250.5833−0.1230-0.1243670.001367-0.072543-0.011111260.5900−0.2100-0.2071640.002836-0.1222260.013507Total IAE0.023441Table7 and Figure 8 show the statistical results and the convergence curve, respectively. From Table 7, it can be seen that BES obtains the best minimum value of RMSE.Table 7
Statistical results for R.T.C. France PV cell two diodes.
AlgorithmRMSEMinMeanMaxStd.BES9.8248e−049.9518e−041.1881e−034.6013e−05ITLBO9.8248e−049.8812e−049.8497e−041.54e−06IMFO9.8252e−049.9737e−041.1409e−033.2939e−05IJAYA9.8380e−041.0240e−031.3507e−038.5647e−05GOTLBO9.8407e−041.4380e−031.0453e−031.01e−04CLPSO9.9224e−041.0522e−031.1462e−034.3141e−05TPTLBO9.8248e−049.8602e−049.8363e−049.31e−07MADE9.8261e−049.8608e−049.8786e−048.02e−05TLABC9.8415e−041.0555e−031.5048e−031.54e−06TLBO1.0069e−031.1598e−031.5206e−031.56e−04EGBO9.8248e−049.8484e−049.8681e−041.66e−04DDSFLA9.8434e−041.1071e−031.4225e−031.3014e−04SFLBS9.8249e−049.8541e−049.8787e−041.7882e−06PSCOS9.8297e−041.0286e−031.4133e−039.9217e−05Figure 8
Convergence curve during the parameter extraction for the two-diode models.


## 5.3. Case Study 3: Photowatt-PWP201 PV Module

To further validate the BES technique and show its effectiveness in estimating the optimal parameters of different models, we used this algorithm on the Photowatt-PWP201 PV module, which consists of 36 silicon cells connected in series under operating conditions of 1000 W/m2 of solar irradiation and a cell temperature of 45°C. The results obtained were compared with those reported in the literature based on other techniques.The results have been listed in Table8. This table also presents a comparison with the results of other techniques from the literature ITLBO [69], IMFO [67], IJAYA [63], and MADE [18]. The comparison validated the effectiveness of BES compared to other techniques. The RMSE based on BES application to extract the parameters of PV model is equal to 0.00242507 which is better. From Figure 9, it can be clearly seen that the I-V and P-V curves of the simulated data found by BES are very compatible with the experimental data.Table 8
Comparison of the results obtained from Photowatt-PWP201 model with other methods in the literature.
AlgorithmIph (A)Isd (μ A)Rs (Ω)Rsh (Ω)nRMSEBES1.03053.48231.2013981.982448.64282.42507e−03GAMS [3]1.03203.26811.2062828.29281.34452.4426e−03ITLBO [69]1.03053.48231.2013981.982348.64282.4251e−03IMFO [67]1.03053.47831.2013980.467248.63852.4251e−03IJAYA [63]1.03023.47031.2011984.876048.64822.4251e−03MADE [18]1.03053.48231.2013981.982348.64282.4251e−03EVPS [52]1.03183.26791.2066845.7591.34452.4267e−03GOTLBO [64]1.03053.49911.2008989.688948.66112.4251e−03TLABC [20]1.03063.47151.2017972.935748.63132.4251e−03TLBO [71]1.03053.48721.2011984.876048.64822.4251e−03CLPSO [70]1.03043.61311.19781017.048.78472.4280e−03TPTLBO [72]1.03053.48231.2013981.982248.64282.4251e−03EGBO [56]1.03053.481.2013981.982248.64282.4151e−03DDSFLA [57]1.03063.44731.2023971.250048.60402.4252e−03SFLBS [58]1.03053.48221.2012981.980448.64282.4251e−03LCNMSE [59]1.03153.48221.2013981.974148.64282.4251e−03PSOCS [60]1.03053.48231.2013981.9848.6432.4251e−03Figure 9
Characteristics of the measured and estimated curve with one diode: (a)I-V and (b) P-V.
(a)
I-V characteristic curve(b)
P-V characteristic curveFrom Table9, it is evident that all IAE values are lower than 4.418e−03 and the RE values are between −6.44e−03 and 1.3167e−01, demonstrating the high efficiency identified by BES for the PV model.Table 9
The individual absolute error and relative error of Photowatt-PWP201 obtained by BES.
ItemTest dataSimulated currentPowerV (V)I (A)Isim (A)IAE (A)Psim (W)RE10.12481.03151.0291050.0023950.1284320.00232221.80931.03001.0273670.0026331.8588150.00255733.35111.02601.0257280.0002723.4373150.00026644.76221.02201.0240930.0020934.876936-0.00204856.05381.01801.0222780.0042786.188666-0.00420267.23641.01551.0199180.0044187.380532-0.00435078.31891.01401.0163520.0023528.454931-0.00232089.30971.01001.0104900.0004909.407354-0.000485910.21631.00351.0006310.00286910.2227480.0028591011.04490.98800.9845670.00343310.8744390.0034751111.80180.96300.9595670.00343311.3246160.0035651212.49290.92550.9229260.00257411.5300230.0027811313.12310.87250.8727480.00024811.453153-0.0002841413.69830.80750.8075040.00000411.061434-0.0000051514.22210.72650.7286690.00216910.363202-0.0029851614.69950.63450.6375920.0030929.372288-0.0048741715.13460.53450.5368060.0023068.124337-0.0043131815.53110.42750.4302530.0027536.682303-0.0064401915.89290.31850.3196740.0011745.080544-0.0036862016.22290.20850.2084500.0000503.3816630.0002402116.52410.10100.0973910.0036091.6092940.0357362216.7987−0.0080-0.0069470.001053-0.1166940.1316702317.0499−0.1110-0.1094040.001596-1.8653230.0143802417.2793−0.2090-0.2075660.001434-3.5865890.0068632517.4885−0.3030-0.2990420.003958-5.2297960.013063Total IAE0.054687Table10 and Figure 10 show the statistical results and the convergence curve, respectively. From Table 10, it can be seen BES obtains the best minimum value of RMSE.Table 10
The statistical results for Photowatt-PWP201.
AlgorithmRMSEMinMeanMaxStd.BES2.42507e−32.42507e−32.42507e−32.4532e−17ITLBO2.4251e−032.4251e−032.4251e−031.27e−17IMFO2.4251e−032.4294e−032.5005e−031.3831e−05IJAYA2.4251e−032.4393e−032.4289e−033.78e−06GOTLBO2.4251e−032.4852e−032.4419e−031.38e−05CLPSO2.4280e−032.4549e−032.5432e−032.5809e−05TPTLBO2.4251e−032.4251e−032.4251e−031.20e−17MADE2.4250e−032.4251e−032.4251e−031.96e−16TLABC2.4251e−032.4265e−032.4458e−034.00e−06TLBO2.4251e−032.4383e−032.5475e−032.43e−05EGBO2.4251e−032.4251e−032.4251e−032.38e−17DDSFLA2.4252e−032.4974e−032.9281e−031.1000e−04SFLBS2.4251e−032.4251e−032.4251e−033.9417e−17PSOCS2.4251e−032.4252e−032.4282e−035.9113e−07Figure 10
Convergence curve during the parameter extraction for the two-diode models.


## 5.4. Case Study 4: Schutten Solar STM6-40/36 Monocrystalline PV Module

Here, we use the BES algorithm to extract the parameters of the Schutten Solar STM6-40/36 PV module. It contains 36 polycrystalline cells (size156mm×156mm) connected in series. The data set contains 20 data points measured at T=51°C [66]. For the STM6-40/36 PV module model, Table 11 shows the results of the parameters obtained from ITLBO, IJAYA, GOTLBO, TPTLBO …. From the results, it can be seen that the BES provides a better RMSE: 0.00172981370994066.Table 11
Comparison of the results obtained from Solar STM6-40/36 PV model with other methods in the literature.
AlgorithmIph (A)Isd (μA)Rs (Ω)Rsh (Ω)nRMSEBES1.66391.73870.004315.92831.52031.7298e−03ITLBO [69]1.66391.73870.004315.92831.52031.7298e−03IJAYA [63]1.66371.83530.004015.94491.52631.7548e−03MADE [18]1.66391.13870.004315.92831.52031.7298e−03GOTLBO [64]1.66391.73870.004315.92831.52031.7298e−03TPTLBO [72]1.66391.73870.004315.92831.52031.7298e−03TLABC [20]1.66391.73870.004315.92831.52031.7298e−03TLBO [71]1.66381.73070.004315.99551.51981.7305e−03EGBO [56]1.66391.730.004315.92831.52031.7298e−03In addition, to confirm the accuracy of the extracted parameters, Figure11 shows the I-V and P-V curves. It is evident that the simulated data from the BES match well with the measured data in the voltage range for both I-V and P-V curves.Figure 11
Characteristics of the measured and estimated curve Schutten Solar STM6-40/36 PV model: (a)I-V and (b) P-V.
(a)
I-V characteristic curve(b)
P-V characteristic curveIn addition, the IAE and RE (relative error) are given in Table12. The IAE describes the error between the extracted parameter and the measured data. In other words, the extracted parameters are better when the IAE is small. According to Table 12, the sum of the IAE is less than 2.50E−02, which indicates that the measured and extracted data coincide well.Table 12
The individual absolute error and relative error of Schutten Solar STM6-40/36 PV model obtained by BES.
ItemTest dataSimulated currentPowerV (V)I (A)Isim (A)IAE (A)Psim (W)RE10.0001.6631.6634530.0004530.000000-0.00027320.1181.6631.6632470.0002470.196263-0.00014932.2371.6611.6595460.0014543.7124040.00087545.4341.6531.6539100.0009108.987346-0.00055057.261.651.6505610.00056111.983074-0.00034069.681.6451.6454260.00042615.927728-0.000259711.591.641.6392320.00076818.9986960.000468812.61.6361.6337140.00228620.5848020.001397913.371.6291.6272930.00170721.7569020.0010481014.091.6191.6183280.00067222.8022470.0004151114.881.5971.6031200.00612023.854427-0.0038321215.591.5811.5816420.00064224.657796-0.0004061316.41.5421.5424300.00043025.295849-0.0002791416.711.5241.5213180.00268225.4212170.0017601516.981.51.4993470.00065325.4589140.0004351617.131.4851.4854450.00044525.445677-0.0003001717.321.4651.4658490.00084925.388510-0.0005801817.911.3881.3878870.00011324.8570560.0000811919.081.1181.1190650.00106521.351757-0.0009522022.0200.0023860.0023860.050143—Total IAE0.02486996To further prove the reliability of the BES, the statistical results containing the minimum (Min), maximum (Max), mean value (Mean), and standard deviation (Std.) are analyzed. Table13 shows the statistical results, and Figure 12 shows the convergence curve of Photowatt STM6-40/36.Table 13
The statistical results for Photowatt STM6-40/36.
AlgorithmRMSEMinMeanMaxStd.BES1.7298e−031.7298e−031.7298e−035.6525e−18ITLBO1.7298e−031.7298e−031.7298e−037.13e−03IJAYA17548e−032.5223e−031.9305e−031.91e−04GOTLBO1.7298e−031.1244e−034.2347e−036.41e−02TPTLBO1.7298e−031.7298e−031.7298e−034.96e−18MADE1.7298e−031.7298e−031.7298e−038.49e−14TLABC1.7298e−032.1827e−036.5053e−039.22e−04TLBO1.7305e−034.3487e−032.0593e−023.45e−03EGBO1.7298e−031.7298e−031.7298e−038.22e−18Figure 12
Convergence curve during the parameter extraction for Solar STM6-40/36 PV model.


## 5.5. Case Study 5: STP6-120/36 Module

The polycrystalline STP6-120/36 has 36 cells connected in series and is measured under 1000 W/m2 at 55°C. The current-voltage data was obtained from [72]. The results have been listed in Table 14. This table also presents a comparison with the results of other techniques from the literature ITLBO [69], IMFO [67], IJAYA [63], TLABC [20], GOTLBO [64], and MADE [18]. The comparison validated the effectiveness of BES compared to other techniques. The RMSE based on BES application to extract the parameters of PV model equals 0.0167828544285882. From Figure 13, it can be clearly seen that the I-V and P-V curves of the simulated data found by BES are very compatible with the experimental data.Table 14
Comparison of the results obtained from Solar STP6-120/36 PV model with other methods in the literature.
AlgorithmIph (A)Isd (μ A)Rs (Ω)Rsh (Ω)nRMSEBES7.47162.32180.004623.02651.25961.6782e−03ITLBO [69]7.47252.33500.004622.21991.26011.6601e−02IJAYA [63]7.46722.25360.004627.59251.25711.6731e−02MADE [18]7.47252.33500.004622.21991.26011.6601e−02GOTLBO [64]7.47252.33500.004622.21991.26011.6601e−02TPTLBO [72]7.47573.01000.004622.21991.26013.0343e−01TLABC [20]7.47252.33490.004622.21171.26011.6601e−02TLBO [71]7.47821.91940.004713.26881.24401.6892e−02EGBO [56]7.47252.330.004622.21991.26011.6601e−02Figure 13
Characteristics of the measured and estimated curve SPT6-120/36: (a)I-V and (b) P-V.
(a)
I-V characteristic curve(b)
P-V characteristic curveIn addition, the IAE and RE (relative error) are given in Table15. The IAE describes the error between the extracted parameter and the measured data. In other words, the extracted parameters are better when the IAE is small. According to Table 15, the sum of the IAE is less than 4.64E−02, which indicates that the measured and extracted data coincide well.Table 15
The individual absolute error and relative error of Solar STP6-120/36 PV model obtained by BES.
ItemTest dataSimulated currentPowerV (V)I (A)Isim (A)IAE (A)Psim (W)RE107.487.4701030.00989791.0073840.00132329.067.457.4520810.00208193.400729-0.00027939.477.427.4489540.02895495.858207-0.003902410.327.447.4387160.00128498.0869360.000173511.177.417.4200200.010020100.437238-0.001352611.817.387.3957830.015783101.453230-0.002139712.367.377.3630570.006943101.7408930.000942812.747.347.3312970.008703101.0794300.001186913.167.297.2840580.00594299.7392770.0008151013.597.237.2175820.01241897.2267420.0017181114.177.17.0880200.01198094.5287000.0016871214.586.976.9583830.01161788.4607380.0016671314.936.836.8145270.01547386.0062300.0022651415.396.586.5678640.01213681.8815400.0018441515.716.366.3487760.01122478.4099650.0017651616.0866.0464390.04643974.377396-0.0077401716.345.755.7851100.03511067.754636-0.0061061816.765.275.2780870.0080870.274440-0.0015351916.95.075.0891260.01912691.007384-0.0037722017.14.794.7883940.00160693.4007290.0003352117.254.564.5455050.01449595.8582070.0031792217.414.294.2721080.01789298.0869360.0041712317.653.833.8387900.008790100.437238-0.0022952419.2100.0142860.014286101.453230—Total IAE0.330286Table16 and Figure 14 show the statistical results and the convergence curve, respectively.Table 16
The statistical results for Photowatt STP6-120/36 obtained by BES.
AlgorithmRMSEMinMeanMaxStd.BES1.6782e−021.69034e−021.7223e−021.1219e−04ITLBO1.6601e−021.6601e−021.6601e−027.22e−17IJAYA1.6731e−021.6891e−021.7304e−021.12e−04GOTLBO1.6601e−022.9588e−021.8099e−013.05e−02MADE1.6601e−021.6601e−021.6601e−021.69e−15TLABC1.6601e−021.6963e−022.1497e−029.47e−04TLBO1.6892e−023.6690e−022.1604e−023.51e−02EGBO1.6601e−021.6601e−041.6601e−021.47e−16Figure 14
Convergence curve during the parameter extraction for Solar STP6-120/36 PV model.


## 6. Conclusion

In this paper, we have presented a new and very recent algorithm based on the metaheuristic technique, called the bald eagle search (BES) algorithm to extract the best values of cell and panel parameters. To demonstrate the performance of the algorithm, many cases were implemented using the single-diode, double-diode, and PV panel models. The current-voltage and power-voltage characteristics of the measured and estimated data show the good accuracy of the proposed method. Simulation result after 20 tests and comparisons with other methods show the accuracy and validity of the method for extracting the parameters of a PV cell and module. It has the advantage of producing stable results of each test result and converging rapidly (in less than 50 iterations). The method is verified using practical data from various manufacturers. Its accuracy is confirmed by comparing its RMSE with many metaheuristic methods. In all considered scenarios, a high level of accuracy is obtained. Therefore, the excellent correspondence of the simulatedI-V and P-V curves with the measured characteristics confirms the accuracy of the BES and its applicability to parameter estimation and for solving the optimization problems of other power systems.



---

*Source: 4343203-2021-10-31.xml*