Next Article in Journal
Project-Based STEM Learning Using Educational Robotics as the Development of Student Problem-Solving Competence
Next Article in Special Issue
Computational Modeling of Latent Heat Thermal Energy Storage in a Shell-Tube Unit: Using Neural Networks and Anisotropic Metal Foam
Previous Article in Journal
MACA: Multi-Agent with Credit Assignment for Computation Offloading in Smart Parks Monitoring
Previous Article in Special Issue
Numerical Modelling and Experimental Validation of Novel Para Winglet Tape for Heat Transfer Enhancement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Photovoltaic Models’ Parameter Extraction Using New Artificial Parameterless Optimization Algorithm

1
Department of Electrical Engineering, College of Engineering, Jouf University, Sakaka 72388, Saudi Arabia
2
Department of Electrical Engineering, College of Engineering, Northern Border University, Ar’Ar 73222, Saudi Arabia
3
Department of Computer Engineering and Networks, College of Computer and Information Sciences, Jouf University, Sakaka 72388, Saudi Arabia
4
Centre for Smart Systems, AI and Cybersecurity, Staffordshire University, Stoke-on-Trent ST4 2DE, UK
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(23), 4617; https://doi.org/10.3390/math10234617
Submission received: 11 November 2022 / Revised: 28 November 2022 / Accepted: 1 December 2022 / Published: 6 December 2022

Abstract

:
Identifying parameters in photovoltaic (PV) cell and module models is one of the primary challenges of the simulation and design of photovoltaic systems. Metaheuristic algorithms can find near-optimal solutions within a reasonable time for such challenging real-world optimization problems. Control parameters must be adjusted with many existing algorithms, making them difficult to use. In real-world problems, many of these algorithms must be combined or hybridized, which results in more complex and time-consuming algorithms. This paper presents a new artificial parameter-less optimization algorithm (APLO) for parameter estimation of PV models. New mutation operators are designed in the proposed algorithm. APLO’s exploitation phase is enhanced by each individual searching for the best solution in this updating operator. Moreover, the current best, the old best, and the individual’s current position are utilized in the differential term of the mutation operator to assist the exploration phase and control the convergence speed. The algorithm uses a random step length based on a normal distribution to ensure population diversity. We present the results of a comparative study using APLO and well-known existing parameter-less meta-heuristic algorithms such as grey wolf optimization, the salp swarm algorithm, JAYA, teaching-learning based optimization, colliding body optimization, as well as three major parameter-based algorithms such as differential evolution, genetic algorithm, and particle swarm optimization to estimate the parameters of PV the modules. The results revealed that the proposed algorithm could provide excellent exploration–exploitation balance and consistency during the iterations. Furthermore, the APLO algorithm shows high reliability and accuracy in identifying the parameters of PV cell models.

1. Introduction

1.1. Background

Because of solar energy’s outstanding environmental, technical, and economic properties, increasing the integration of solar photovoltaic systems with electric utilities is inevitable [1]. Solar radiation is abundant in most areas of the world. These systems can be used for energy generation, enabling customers to invest quickly in their electrical systems. Photovoltaic (PV) systems convert solar energy into electricity. Solar energy’s potential for generating electricity depends on various factors, including temperature and solar radiation [2,3]. Therefore, it is vital to assess how PV systems perform in operation to be modeled, managed, and optimized for future operations [4].
The inherent characteristics of PV cells give rise to a nonlinear power–voltage (P–V) curve, which is highly influenced by the environment. Ideally, a PV panel should operate at maximum efficiency at the peak of the curve. A proper maximum power point tracking (MPPT) technique enables PV systems to be used more efficiently in different environmental conditions, requiring accurate and reliable data regarding PV parameters [5,6]. Unknown parameters are the most significant contributors to PV models. It is, therefore, crucial to identify the model’s unknown parameters early using a feasible optimization algorithm, regardless of the model used. Accordingly, a robust optimization algorithm to accurately estimate the parameters of PV models and track the maximum power point under different conditions is urgent.

1.2. Related Works

Generally, solar PV systems can be optimized using deterministic and metaheuristic algorithms. Due to their reliance on gradient information and sensitivity to initial points, deterministic algorithms are unreliable. Moreover, because of their nonlinearity, these classical algorithms also have trouble capturing local optima in the nonconvex space the equivalent PV circuits created. The result may be an inaccurate estimation of parameters and, consequently, a failure to track the maximum power point [7,8,9].
Metaheuristic optimization techniques are considered modern and straightforward alternatives to deterministic algorithms. In general, metaheuristic algorithms fall into two main categories: single-solution-based and multiple-solution-based algorithms [10]. It is understood that the former algorithm uses an iterative process to achieve a superior solution by starting with a randomly selected candidate solution and moving and improving it in a promising search space iteratively. A common single-solution metaheuristic algorithm is simulation annealing (SA) [11]. Multiple-solution algorithms employ several random solutions to enhance their performance.
In classifications of population-based metaheuristic algorithms, evolutionary algorithms, physics algorithms, chemistry algorithms, and swarm-based algorithms are all included [10,12]. Based on evolution in nature, evolutionary-based metaheuristic algorithms move populations based on improvements and movements. The physics-based metaheuristic algorithm enhances the initial population through search space by using principles established based on the physics’ lows, such as mechanics, relativity, gravity, electrodynamics, electromagnetism, and optics. A chemical reaction and molecule characteristics are utilized to develop chemistry-based algorithms. Living organisms, including birds, ants, swarms, schools, and so on, are modeled using swarm-based algorithms. Some known population-based metaheuristic algorithms are introduced in Table 1.
Unlike deterministic algorithms, metaheuristic algorithms find near-optimal solutions within reasonable time for challenging real-world optimization problems. As a result, many optimization problems in science and engineering have been solved thanks to easy implementation and efficiency (see sample [49,50,51,52,53,54]). Hence, researchers have been motivated to develop successful algorithms inspired by natural and artificial processes, to solve complex optimization problems.
In recent years, metaheuristic algorithms have been used to optimize PV systems more accurately and flexibly. The development of parameterless metaheuristic algorithms for optimal parameter identification of solar PV cells has been the subject of considerable research. For example, for the extraction of parameters from PV cell-based single and double diode models [55], the salp swarm algorithms (SSA) were used. The JAYA algorithm was developed by Rao [56] as a powerful heuristic for solving optimization problems. It has been demonstrated by [57] that JAYA can be used to estimate PV cell and module parameters based on performance-guided criteria. Through a chaotic learning process, ref. [58] proposed an improved JAYA algorithm (IJAYA) to find PV model parameters reliably and accurately. In Ref. [59], TLBO was numerically simulated interactively and applied to various solar cells. TLBO was improved and simplified by [60] using an elite strategy and a local search to identify the parameters of the solar PV cells. GOTLBO (generalized oppositional TLBO) was derived from generalized opposition-based learning and was employed to identify solar cell models’ parameters [61]. The modified salp swarm algorithm (MSSA) has been used as an efficient metaheuristic for identifying PV model parameters [62]. However, according to the results of this modification, sufficient robust solutions haven’t been achieved for all PV models.
Examples of metaheuristic algorithms with control parameters for optimal parameter estimation of PV cells include DE, GA, and PSO. An estimation method for solar PV module parameters based on penalized differential evolution (P-DE) was proposed by [63]. Adding new scaling factors and crossover rates to adaptive DE improved the parameter estimation [64]. The genetic algorithm was utilized by [65] to determine solar cells’ I-V characteristics. PSOs have been successfully applied in some modified and hybrid forms to identify the parameters of solar cells. For instance, the chaotic heterogeneous comprehensive learning particle swarm optimizer was proposed by [66] for dynamic and static PV models’ parameter identification. In Ref. [67], the flexible particle swarm optimizer was presented for parameter extraction of different PV models. A classified perturbation mutation-based particle swarm optimization was introduced by [68]. In Ref. [69], a niche particle swarm optimization in parallel computing was proposed and applied to identify the unknown parameters of PV cells. A fractional chaotic ensemble particle swarm optimizer was utilized for estimating the parameters of three models of PV cells [70]. Exploitation and exploration phases were balanced effectively to mitigate the premature convergence associated with PSO [71]. A variation of opposition-based GOA (OBGOA) has been proposed by [72] to identify the electrical parameters of various PV models. A modified spotted hyena optimization algorithm MSHOA was proposed by [73] to boost the optimal solution’s performance through an accelerated function. In Ref. [74], a hybrid optimization algorithm called hARS-PS was presented that uses adaptive rat swarm optimization (ARSO) and pattern search (PS) to extract PV parameters.

1.3. Motivation

Estimating the parameters of PV systems has been widely done using metaheuristic algorithms. Nonetheless, researchers must develop algorithms that efficiently take advantage of these two factors. The first issue with most metaheuristic algorithms in PV system optimization is determining their particular control parameters. For instance, DE needs the crossover probability and scaling factor, and GA needs crossover and mutation rates. In PSO, the inertia weight and cognitive parameters should be defined. Choosing an incorrect parameter set can distract the user from the main problem and disrupt the algorithms. A parameter-free algorithm, such as TLBO, JAYA, SSA, colliding bodies optimization (CBO), and GWO, makes optimization more effortless and efficient by avoiding adjusting parameters. Note that a parameter-less algorithm in this paper must only determine standard parameters, including the number of iterations (or function evaluations) and initial population size.
Secondly, no single algorithm can solve all optimization problems, according to Wolpert and Macready’s “No Free Lunch” theorem [75]. Since the theorem was published 30 years ago, researchers have significantly improved metaheuristic algorithms and developed new ones. Nevertheless, many articles related to improving the performance of these original algorithms show their weakness in effectively and reliably detecting the parameters of different photovoltaic models. Consequently, developing new ideas that result in simple and efficient metaheuristic algorithms without requiring additional parameter settings and modifications for practical optimization problems is necessary.

1.4. Contribution

For the reasons outlined above, this paper proposes a metaheuristic algorithm called APLO for the parameter estimation of PV models. The algorithm is efficient and straightforward. In each iteration of the algorithm, each individual is updated using the current best, last best, and individual’s current solutions. The individual moves to a new position if a better solution is obtained or remains unchanged. The process is iterated until all individuals converge on the best solution or a specified criterion is reached. The performance of the proposed APLO for the problems of optimal parameters estimation of PV models is evaluated and compared with five well-known parameterless metaheuristic algorithms, i.e., TLBO, JAYA, SSA, CBO, GWO, and three conventional algorithms such as DE, GA, and PSO. Moreover, in some cases, the related results are also reported from the literature and are compared.
The main contributions of this paper are:
  • A novel, simple, and efficient parameterless optimization algorithm is proposed;
  • For parameter estimation of PV models, APLO is tested in a series of experiments;
  • High accuracy and reliability in finding the PV models’ unknown parameters;
  • Reasonable performance of the proposed algorithm compared with other original, improved, and hybrid metaheuristic algorithms.

1.5. Paper Structure

The remainder of the paper is organized as follows. The mathematical formulation of the proposed APLO algorithm is presented in Section 2. Detailed descriptions of the single-diode, double-diode, and PV module models are presented in Section 3. Section 4 simulates and evaluates the results of the experiment. Lastly, Section 5 summarizes some concluding remarks.

2. APLO Algorithm

2.1. Mathematical Model

The main steps of the proposed APLO algorithm are mathematically described in this section.

2.1.1. Initialization

In most metaheuristic algorithms, an initial population is randomly generated. Then, this population is improved in the problem’s solution space during the iterative process using a proper evolutionary mechanism. APLO is an artificial population-based metaheuristic algorithm that follows the general procedure of evolutionary algorithms such as the DE algorithm. The population consists of npop individuals, which i th individual is represented by P i ( t ) = [ P i , 1 ( t ) ,   ,   P i , D ( t ) ] | i = 1 ,   2 ,   ,   n p o p at iteration t. The initial population can be created as follows:
P i , j ( t ) = p min , j + r a n d . ( p max , j p min , j ) ;   t = 1 ,   i = 1 ,   2 ,   ,   n p o p ; j = 1 ,   2 ,   ,   D
where p min , j and p max , j define the lower and upper bounds of the problem’s decision variable j . r a n d is a uniform distribution value, and D is the dimension of the problem.

2.1.2. Search Operator

After generating the initial population, they are committed to searching around the current best, P B ( t ) , using the mean knowledge obtained from the last best, P L ( t ) , to find a better solution. Therefore, for each individual i , the new position of its arbitrary elements, i.e., Υ i , P i , j ( t + 1 ) , is updated using the following equation:
P i , j ( t + 1 ) = P B , j ( t ) + r 1 , i r 2 , i ( P B , j ( t ) P L , j ( t ) + P i , j ( t ) 2 ) ;   i = 1 ,   2 ,   ,   n p o p ; j Υ i
where r 1 , and r 2 are chosen randomly from the normal distribution function. The random value r 1 / r 2 is the main factor for balancing the exploration and exploitation phases of the proposed APLO algorithm. An example of this random value over the 1000 iterations is shown in Figure 1. Large values cause P i ( t + 1 ) to leave a feasible space for the problem. Hence, by applying restrictions related to each variable’s upper and lower bounds, the violent variables are replaced by new feasible solutions. This preserves diversity in the population during the algorithm process. Thus, it helps the algorithm to search the solutions globally (exploitation phase). In contrast, small random values assist the algorithm in searching locally around the current best solution, helping the exploration phase. It is worth mentioning that an important factor in creating stable behavior in the proposed algorithm is that the updating operation is applied only to a part of the variables of each individual.
A representative example of updating the equation in 2D space is shown in Figure 2.
Some elements of the test vector may exceed their allowed range. Various algorithms may be used to face this challenge. In this article, a simple algorithm is considered for it, such that every element of the new vector that is out of the lower or upper range is replaced with a random value as follows:
P i , j ( t + 1 ) = p min , j + r a n d . ( p max , j p min , j ) ;   if   P i , j ( t + 1 ) < p min , j P i , j ( t + 1 ) = p min , j + r a n d . ( p max , j p min , j ) ;   if   P i , j ( t + 1 ) < p min , j
This algorithm can effectively create diverse solutions during the algorithm process. In addition, random long steps created by r 1 / r 2 can contribute to this feature of the algorithm and thus prevent its premature convergence.

2.1.3. Selection

The selection mechanism choices the better solution from the old position P i ( t ) and the newly generated position P i ( t + 1 ) based on their fitness values, i.e., f ( . ) . Hence, for instance, the following selection operator is utilized for the minimization problem:
P i ( t + 1 ) = { P i ( t + 1 ) ; if   f ( P i ( t + 1 ) ) < f ( P i ( t ) ) P i ( t ) ; else .
P B ( t ) and P L ( t ) continually are updated after each function evaluation. The procedure above continues until the termination conditions are met.
The flowchart of the proposed APLO algorithm is outlined in Figure 3, and the pseudocode of the APLO algorithm (Algorithm 1) is defined as follows:
Algorithm 1: The pseudo-code of APLO
Input:  M a x I t e r ,   n p o p
Output: the best solution
1:Initialize the n p o p population randomly using Equation (1)
2:Calculate the fitness values of all individuals
3:Determine the current best and last best solutions
4: t 1
5:While  t < M a x ¯ I t e r  do
6:for i = 1 : n p o p  do
7:Update some arbitrary elements of individual i using Equations (2) and (3)
8:Calculate the fitness of individual i
9:Accept the updated solution if it is better than the old one using Equation (4)
10:Update the last best and current best solutions
11: t t + 1
12:end
13:end

2.2. Algorithm Complexity

An algorithm’s complexity plays a vital role in assessing its performance. As with all metaheuristic algorithms, APLO requires O ( n × n p o p ) times to initialize each population time, where n indicates the number of objective functions and n p o p represents the number of populations. Every algorithm has an O ( M a x ¯ I t e r × f c ) complexity, where M a x ¯ I t e r is the predefined maximum number of iterations. For a given problem, f c represents the complexity of the evaluation function. Simulating the entire process would require O ( N ) . Accordingly, the algorithm has a computational complexity of O ( N × M a x ¯ I t e r × f c × n × n p o p ) .

3. The Problem of PV Models’ Parameter Extraction

This section introduces three standard models of photovoltaic cells, i.e., SDM, DDM, and PV module models. Then, the mathematical model of the problem of finding the optimal parameters of these three PV models is expressed.

3.1. Single-Diode Model (SDM)

PV systems must be mathematically modeled, considering their practicalities to show their real-time characteristics. PV arrays can be modeled from their basic unit, which is the cell. Due to its simplicity and ease of implementation, the SDM is popular. A parallel resistor, a series resistor, a diode, and a current source make up the equivalent circuit for the single diode model, as shown in Figure 4. To determine the output current, I P V , the following formula can be used [65]:
I P V = I p h ( I s h + I d )
where the photogenerated, the diode, and the shunt resistor currents, respectively, are represented by I p h , I d , and I s h . Shockley’s equation and Kirchhoff’s voltage law (KVL) can be used to calculate I d and I s h as follows:
I d = I s d [ exp ( V P V + R s I P V m v ) 1 ]
I s h = V P V + R s I P V R s h
where I s d indicates the diode reverse saturation current; I P V denotes the cell output voltage; R s h and R s represent the shunt and series resistances, respectively; moreover, the non-physical diode ideality factor is defined by m . v in Equation (6) is the junction thermal voltage that can be expressed as follows:
v = k T q
Boltzmann’s constant, k , is 1.8865033 × 10 23 J/K, T is junction temperature, and q is electron charge ( 1.60217646 × 10 19 C). The output current I P V can be expressed in the following manner by combining Equations (5)–(8):
I P V = I p h V P V + R s I P V R s h I s d [ exp ( V P V + R s I P V m v ) 1 ]
The single-diode model requires the identification of five unknown parameters, that includes I p h , I s d , R s h , R s , and m .

3.2. Double Diode Model (DDM)

In SDM, recombination current is ignored in the depletion region despite being widely used to simulate PV cells [37]. DDM solves this problem by having three components: a photo-generated current source, a shunt resistance, two rectifying diodes, and a series resistance, as shown in Figure 5.
As a result of KCL, the output current I P V in DDM can be calculated as follows:
I P V = I p h ( I s h + I d 2 + I d 1 )
As shown in Figure 4, I d 1 is the current flowing through the first diode (D1), and I d 2 is the current flowing through the second diode (D2). It is also possible to express the magnitude of these currents in terms of the Shockley diode equation as follows:
I d 1 = I s d , 1 [ exp ( V P V + R s I P V m 1 v ) 1 ]
I d 2 = I s d , 2 [ exp ( V P V + R s I P V m 2 v ) 1 ]
There are two ideality factors for diodes, m 1 and m 2 , as well as diffusion and saturation currents, I s d , 1 and I s d , 2 , respectively. Hence, Equation (10) can be rewritten as follows by substituting Equations (7), (11) and (12):
I P V = I p h V P V + R s I P V R s h I s d , 2 [ exp ( V P V + R s I P V m 2 v ) 1 ] I s d , 1 [ exp ( V P V + R s I P V m 1 v ) 1 ]
Compared to SDM, DDM requires more parameters to be identified, i.e., I p h , I s d , 2 , I s d , 1 , R s h , R s , m 1 , and m 2 .

3.3. PV Module Model

Multiple PV cells are arranged parallel or in series in PV modules to increase current and voltage. The equivalent circuit model of the PV module is shown in Figure 6. The PV module model output current can be calculated as follows:
I P V = N p I p h V P V + [ N s / N p ] R s I P V [ N s / N p ] R s h N p I s d [ exp ( V P V + [ N s / N p ] R s I P V m v ) 1 ]
where N p represents the number of solar cells in parallel and N s represents the number of solar cells in series. The PV module has five unknown parameters, similar to the SDM ( I p h , I s d , R s h , R s , and m ).

3.4. Problem Formulation of PV Models’ Parameters Extraction

Based on measurements of I-V from actual PV cells and PV modules, mathematical models for PV models aim to estimate unknown parameters with remarkable accuracy. Experimental and built I-V data differences are commonly minimized with an optimization technique. It is, therefore, common to consider the minimization of root mean square error (RMSE) between the estimated current ( I P V , n ) and experiment current ( I ^ P V , n ) as an objective function:
Minimize   RMSE = ( 1 N n = 1 N ( I ^ P V , n I P V , n ) 2 ) 1 / 2
S.t.
I ¯ p h I p h I ¯ p h
R ¯ s R s R ¯ s
R ¯ s h R s h R ¯ s h
I ¯ s d I s d I ¯ s d
m ¯ m m ¯
I ¯ s d I s d , i I ¯ s d
m ¯ m i m ¯
where N is the number of experimental data. I ^ P V , n and I P V , n are the nth measured samples and the calculated value of PV output current based on each model. Constraints (16)–(20) indicate the upper and lower limits on the decision variables (unknown parameters) for the SDM and PV module model. Constraints (16)–(18), (21) and (22) indicate the upper and lower limits on the decision variables for the DDM.
In Equation (15), the estimated PV output current at each measured sample n can be calculated using Equations (23)–(25) for SDM, DDM, and PV module models, respectively.
I P V , n = I p h V ^ P V , n + R s I ^ P V , n R s h I s d [ exp ( V ^ P V , n + R s I ^ P V , n m v ) 1 ]
I P V , n = I p h V ^ P V , n + R s I ^ P V , n R s h I s d , 2 [ exp ( V ^ P V , n + R s I ^ P V , n m 1 v ) 1 ] I s d , 1 [ exp ( V ^ P V , n + R s I ^ P V , n m 2 v ) 1 ]
I P V , n = N p I p h V ^ P V , n + [ N s / N p ] R s I ^ P V , n [ N s / N p ] R s h N p I s d [ exp ( V ^ P V , n + [ N s / N p ] R s I ^ P V , n m v ) 1 ]

4. Experimental Results

In this section, we evaluate APLO’s effectiveness for parameter estimation with three types of PV models: SDM ( N s = N p = 1 ), DDM ( N s = 1 ,   N p = 2 ), and the PV module ( N s = 36 ,   N p = 1 ). As a standard experiment for SDM and DDM, current–voltage data were collected on silicon solar cells with a diameter of 57 mm (R.T.C. France) [76]. The PV cell characteristics are as follows: V o c = 0.5728   V ,   I s c = 0.7603   A ,   V m = 0.4507   V ,   and   I m = 0.6894   A . In addition, a PV module (Photo Watt-PWP 201) with 36 polycrystalline PV cells is used under 1000 W/m2 irradiance [76]. This PV module’s characteristics are as follows: V o c = 16.778   V ,   I s c = 1.030   A ,   V m = 12.649   V ,   and   I m = 0.912   A . A wide range of algorithms has been developed to estimate the parameters of PV models based on experimental data. Table 2 provides the minimum and maximum limits for PV model parameters [58].
Moreover, eight well-known algorithms, including DE [17], GA [13], PSO [33], the original parameterless algorithms such as GWO [44], TLBO [77], JAYA [56], SSA [78], CBO [79], are selected to validate and verify the effectiveness of APLO. It is assumed that the population size and the maximum number of iterations are set to 50 and 1000 (i.e., 50,000 evaluations of each function), respectively. The other parameters of the algorithm are maintained as they were in the original literature. To perform the statistical analysis, each algorithm is run 30 times independently in MATLAB 2021b.

4.1. Exploration and Exploitation Analysis

One of the effective factors in creating a balance between exploration and exploitation is ensuring sufficient diversity among individuals. This can prevent an algorithm from getting trapped in locally optimal solutions and result in a better solution. In this section, some experiments are performed to evaluate the exploration–exploitation and diversity in the nine applied metaheuristics on the SDM problem. The percentage of exploration and exploitation, visualizing the two abilities and population diversity in the individuals of the competitive algorithms through the iterations, is shown in Figure 7 and Figure 8, respectively. These numerical measures are calculated based on the procedure reported in [80].
It can be observed from Figure 7 that, in GWO and SSA algorithms, unlike other algorithms, the percentage of exploration during iterations is much higher than the exploitation. This is evidenced by Figure 9, where the average values of the exploration-exploitation in these two algorithms are 72%:28% and 83%:17%, respectively. As a result, these two algorithms are explorative. Moreover, regarding population diversity, SSA exhibits the highest values during the iterations (see Figure 8). GWO gives the second-highest diversity among the applied algorithms. In contrast, the most exploitation and most minor exploration capabilities are provided by CBO and GA, with mean values of 97% and 3%, respectively. This is further evidenced by the diversity measures illustrated in Figure 8. These two algorithms couldn’t provide good diversity throughout the iterations. Hence, premature convergence is the main weakness of these algorithms.
Moreover, all algorithms except SSA and GWO, as shown in Figure 7, are explorative at the beginning. After a few iterations, they were threatened as exploitative algorithms. Additionally, the diversity of the population in these algorithms is high initially. After a few iterations, it drops and remains approximately consistent.
The exploration–exploitation and diversity measurements alone cannot show the superiority of an algorithm in solving an optimization problem compared to other algorithms. A better understanding of this goal can be gained by analyzing the characteristic of convergence of these algorithms as shown in Figure 10. The exploration–exploitation balance of 18%:82% in APLO provides the best convergence performance. The final best solution obtained by the proposed algorithm is 9.86 × 10−4, while the best solution obtained by the TLBO algorithm is 1.20 × 10−3. TLBO can provide an exploration–exploitation ratio of 17%:83%, which is very close to the proposed algorithm. The explorative algorithms perform the worst performance, i.e., SSA and GWO, which fall into the local optimal after a few iterations.

4.2. Population Size Analysis

The performance of the metaheuristics in solving specific optimization problems can be affected by the population size. Six different populations of 10, 20, 30, 40, 50, and 60 are evaluated for our algorithm while solving the SDM, DDM, and PV module models’ parameter identification. The APLO algorithm was run 30 times for each population size independently, and the stop criteria were set at 50,000 function evaluations. The statistical results of this experiment in terms of minimum (Min), average (Mean), maximum (Max), and standard deviation (SD) are summarized in Table 3. Moreover, a Freidman rank test is also applied to compare the algorithm’s performance in different cases of the population size. As can be seen form Table 3, the best-suited population size on SDM and DDM problems is npop = 40, with mean rank in Freidman of 2.433 and 3.1, respectively. While the population size of 50 exhibits the best when solving the PV module problem with a mean rank of 3.1. Overall, from the last row of Table 3, the sum of Freidman tests over three PV models indicates that the population size of 40 results in the best performance to solve these problems. The population size of 30 and 50 are the second-best and third-best options. However, to perform a fair comparison with other state-of-the-art algorithms, the population size of 50 is adapted in the following sections.

4.3. Results of Parameter Extraction Based on SDM

In order to investigate the silicon solar cell model developed by RTC France, we analyzed it using SDM. We solved it competitively using nine individual algorithms. Statistics over 30 runs are shown in Figure 11, along with the min, mean, max and SD in Table 4. The bolded values indicate the best results among the applied algorithms.
A glance at the results of Table 4 and the boxplots provided in Figure 11 reveals that the proposed algorithm can achieve better statistical results than other algorithms. For example, the value of the best solution found by the proposed algorithm is 9.860218 × 10−4, while the best solution of the following algorithm, PSO, is 9.86901 × 10−4. TLBO algorithm achieves the solutions close to PSO. Also, from the perspective of the robustness of the solutions, as shown in Figure 6 and the standard deviation of Table 4, the superiority of the proposed algorithm over other competitive algorithms is visible. The SSA algorithm is shown to have the worst performance among the applied algorithms. As shown in the last column of Table 4, the proposed algorithm demonstrates superiority over other algorithms when compared using Wilcoxon’s rank sum test at a 5% confidence level.
Furthermore, Table 5 contains the best-estimated parameters obtained by APLO and other applied algorithms. As a result of these optimal parameters obtained by APLO, the estimated and measured values of current and power, as well as their individual absolute errors (IAEs), are shown in Figure 7. APLO’s simulation provides I-V and P-V characteristics that are highly similar to those of standard data. Figure 12 shows, for instance, that the IAEA of the current range from 2.51 × 10 3 to 8.77 × 10 5 , while the IAEA of the power range from 1.46 × 10 3 to 1.97 × 10 6 , demonstrating that the APLO estimates are highly accurate.

4.4. Results of Parameter Extraction Based on DDM

Based on DDM, nine algorithms are used to solve the silicon solar cell model of RTC France. The summary of the statistical results over 30 runs can be found in Figure 13 and Table 6. Among the applied algorithms, the most desirable values are highlighted in bold. As seen in Table 6 and the boxplot provided in Figure 8, the proposed approach outperforms its competitors in terms of statistical results. For example, the value of the best solution found by the proposed algorithm is 9.83065 × 10−4, while the value of the best solution found by the second-best algorithm, i.e., PSO, is 9.84872 × 10−4. As shown in Figure 13 and the fifth column of Table 6, the proposed algorithm has an advantage over other competitive algorithms from a robustness standpoint. Despite this, the DE and SSA show reasonable robustness but fail to find optimal solutions since they fall into the local optimum at every run. A comparison of Table 6 reveals that the SSA has the least effective performance. Moreover, Table 6 demonstrates the superiority of the proposed algorithm over the other algorithms based on a pairwise comparison using Wilcoxon’s rank sum test at a 5% confidence level.
Furthermore, Table 7 shows the best-estimated parameters for DDM from APLO and other algorithms. Figure 9 shows the calculated and measured values of current and power based on these optimal parameters. According to APLO’s simulation, the I-V and P–V characteristics are highly similar to those in the measured data. Figure 14 shows the high accuracy of the estimated parameters by APLO for DDM. It shows the maximum IAE index of current equaling 0.0025 and the maximum IAE of power equaling 0.0015.

4.5. PV Module Model-Based Photo Watt-PWP 201

This section uses the PV module model to model the Photo Watt-PWP 201 solar cell. In contrast, nine algorithms are used to solve the model. Based on 30 runs of this model, Figure 15 and Table 8 illustrate the statistical results of these algorithms. Figure 15 provides box plots showing how the proposed algorithm outperforms other competitor algorithms in terms of statistical performance. In this study, the APLO algorithm provided the best RMSE at 2.42507486809 × 10−3, followed by PSO at 2.47057268706 × 10−3. Based on the SD results in the fifth column of Table 8, the proposed algorithm appears more reliable than other competitive algorithms. DE obtained the second-best SD value, 7.13076666106 × 10−4, while APLO obtained the best SD value, 5.96208271747 × 10−17.
In contrast to other applied algorithms, this reasonable difference in SD indicates that the APLO algorithm is significantly more robust than different algorithms. The GA has the worst SD performance of the applied algorithms, and the SSA has the worst Min value performance. Moreover, SSA exhibits the worst performance among all algorithms. Table 8’s last column shows that the proposed algorithm is superior to other algorithms when compared pairwise using Wilcoxon’s rank sum test at a 5% confidence level.
Further analysis of the best-estimated parameters of the PV module model generated by APLO and other applied algorithms is presented in Table 9. The estimated and measured current and power, along with their IAE, are displayed in Figure 16 using these optimal parameters. Based on the simulation results obtained by APLO, the I-V and P-V characteristics are very similar to those measured. In this case, the IAE of current and power are less than 0.006 and 0.0799, respectively. This demonstrates that the estimated parameters by APLO for the PV module model are highly accurate.

4.6. Comprehensive Comparison

4.6.1. Convergence Characteristics

A comparison of the convergence curves of different algorithms is shown in Figure 17, Figure 18 and Figure 19 for SDM, DDM, and PV module models under the best run condition (e.g., achieving the “Min” solution). When solving the SDM and PV module models, APLO, CBO, and GA show the fastest convergence rates in the first 50 iterations. Unlike the proposed algorithm (APLO), CBO, and GA come close to obtaining locally optimal solutions after iteration 200. The convergence performance of CBO and GA is insufficient regarding DDM. PSO and APLO, on the other hand, provide better convergence performance than other algorithms used for DDM. However, APLO can converge to a better final solution than PSO in this model after iteration 500. Within the first 500 iterations, the convergence speeds of GA and PSO for the PV module model are faster than APLO. However, APLO can achieve a better final solution than GA and PSO.
It appears that APLO offers an appropriate balance between exploration and exploitation based on the convergence curves for the three models. It should be noted that this comparison is based on the convergence curves presented by the best runs of the applied algorithms. On the other runs, other algorithms are more likely to get stuck in local optimal points and exhibit significantly weaker convergence performance than the proposed algorithm.

4.6.2. Computational Time

There is little difference between APLO, CBO, GWO, JAYA, and SSA regarding their computational time when solving three models, as shown in Figure 20, Figure 21 and Figure 22. They are also faster than DE, GA, and PSO algorithms. The CBO offers minor computational costs, whereas the TLBO is the most complex parameterless algorithm. It is necessary to teach and learn at different times for each generation of TLBO, which results in two different function evaluations (FEs) for each participant. TLBO is thus more expensive to compute in a single generation than an algorithm with one FE per generation. Compared to parameterless algorithms, GA, DE, and PSO require more computational time. PSO is ranked second in this regard, and GA consumes the most time.

4.6.3. Wilcoxon and Friedman Tests

Different algorithms are commonly tested for significant differences in results using the Wilcoxon signed-rank test [13,50,52]. Based on 30 runs of each algorithm, Wilcoxon signed-rank tests were performed. Table 10, Table 11 and Table 12 summarize the results.
In these tables, “mean_PRs”/”sum_PRs”/”N_PRs” and “mean_NRs”/” sum_NRs”/”N_NRs”, respectively, indicate the mean/sum/number of ranks for the problem in which Algorithm I outperformed Algorithm II. The mean/sum/number of ranks for the problem in which Algorithm II outperformed Algorithm I. p-values represent the degree of significance, where a smaller value indicates a more significant difference. Based on the comparison of all models considered in Table 13, it is evident that the proposed APLO differs significantly from the compared algorithms.
Based on the Freidman test, APLO takes the top ranking for three models, demonstrating that the proposed algorithm is superior to those compared. PSO and TLBO rank second and third in average rankings, respectively. This algorithm offers a promising approach to extracting accurate parameters for various solar cells. Aside from this, SSA ranks lowest.

4.6.4. Advantages and Disadvantages of Applied Algorithms

Based on the studies conducted so far, the advantages and disadvantages of the used algorithms are compared from the global search ability, convergence speed, local entrapment probability, exploration/exploitation capability, diversity of individuals, balancing in exploration–exploitation, and computational time points of view. The summary of these comparisons is represented in Table 14. As seen from this table, the proposed APLO shows good performance. However, this algorithm is at the beginning of the way. It can still be examined from different aspects and on various problems so that its challenges are well-known and solved.

4.7. Comparison with the State-of-the-Art Algorithms

Several state-of-the-art improved and hybrid algorithms, including PGJAYA [57], IJAYA [58], STLBO [60], GOTLBO [61], MSSA [62], hARS-PS [74], BLPSO [81], CLPSO [82], TLABC [83], DE/BBO [84], and CMM-DE/BBO [85], are used to validate the performance of the proposed APLO algorithm in identifying PV models’ parameters. Table 15 shows the max, min, mean, and SD of the parameters identified by the algorithms on each model corresponding to RMSE. Bold highlights indicate the best results. APLO yields the best results for the single diode and PV module models in terms of max, min, mean, and SD. It should be noted that PGJAYA and TLABC can only obtain the best results in terms of Min for SDM. Moreover, PGJAYA, TLABC, STLBO, MSSA, and CMM-DE/BBO can achieve the best solution in terms of the min value in calculating the parameters of the PV module model. It should be noted that, the results of hARS-PS on SDM model seems incorrectly reported in [74]. For DDM, hARS-PS is the best algorithm for all statistical values. Regarding the Min value for this model, PGJAYA is the second-best, STLBO is the third best, and IJAYA is the fourth-best algorithm. However, in terms of the mean, max and SD, proposed APLO is the third-best algorithm. The results indicate that the performance of the proposed basic algorithm is acceptable compared to other combined and improved algorithms.

5. Conclusions and Future Directions

This paper proposed a novel parameterless algorithm to estimate parameters without specifying any control parameters, called artificial parameterless optimization (APLO). As part of the proposed APLO, a novel mutation operator was designed. To advance the exploration phase of the APLO, this operator required all participants to move around the best available solution. The current best, the old best, and the individual’s current position were incorporated into the differential term of the mutation operator to assist the exploitation phase and maintain convergence speed. Furthermore, a random multiplication term using a normal distribution was proposed to ensure population diversity through the iteration process of the algorithm. Comprehensive studies were established to evaluate the performance of the proposed algorithm to estimate the parameters of the PV cells. The results revealed that the proposed algorithm could provide an excellent exploration–exploitation balance and consistency during the iterations. Two main factors were responsible for this: (a) the endorsement of positive feedback from individuals who had already achieved improvements, and (b) the presence of well-representative individuals ensured that the entire population was consistent in its diversity.
In addition, some comparisons were made in terms of statistical analysis. Based on the experimental parameter estimation results for the SDM, DDM, and PV module models, the APLO algorithm was more accurate and reliable than other original and well-known algorithms. A comparison of the calculated and standard data of the V-I and P-V curves was also conducted to determine the accuracy of the estimated parameters using the proposed algorithm. Moreover, the results indicated the proper performance of the proposed basic algorithm compared to other state-of-the-art combined and improved algorithms. Regarding the four statistical metrics, the proposed algorithm outperforms others for SDM and PV module models. Additionally, in the case of DDM, it showed good performance compared to the comparative algorithms. However, the results can be further improved by improving or combining the proposed algorithm with other algorithms.
Since APLO is simple and efficient, it can also be used to solve more complex engineering optimization problems. A modification of the APLO can also speed up convergence and reduce computational costs. Moreover, we intend to develop binary and multiobjective versions of APLO algorithms. Furthermore, APLO can be used to optimize support vector machines or kernel extreme learning machines.

Author Contributions

All the authors contributed to formulating the research idea, algorithm design, result analysis, writing and reviewing the research. Conceptualization, M.A. and A.A. (Abdulaziz Alanazi); methodology, M.A.; software, A.A. (Ahmad Almadhor); validation, H.T.R. and A.A. (Abdulaziz Alanazi); formal analysis, M.A.; investigation, A.A. (Ahmad Almadhor); resources, A.A. (Abdulaziz Alanazi); data curation, H.T.R.; writing—original draft preparation, M.A.; writing—review and editing, H.T.R.; visualization, A.A. (Ahmad Almadhor); supervision, M.A.; project administration, A.A. (Abdulaziz Alanazi); funding acquisition, M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Deanship of Scientific Research at Jouf University under Grant Number DSR2022-RG-0113.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not Applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sheng, R.; Du, J.; Liu, S.; Wang, C.; Wang, Z.; Liu, X. Solar Photovoltaic Investment Changes across China Regions Using a Spatial Shift-Share Analysis. Energies 2021, 14, 6418. [Google Scholar] [CrossRef]
  2. Leitão, D.; Torres, J.P.N.; Fernandes, J.F.P. Spectral Irradiance Influence on Solar Cells Efficiency. Energies 2020, 13, 5017. [Google Scholar] [CrossRef]
  3. Wang, P.; Yu, P.; Huang, L.; Zhang, Y. An Integrated Technical, Economic, and Environmental Framework for Evaluating the Rooftop Photovoltaic Potential of Old Residential Buildings. J. Environ. Manag. 2022, 317, 115296. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, H.; Jiao, S.; Heidari, A.A.; Wang, M.; Chen, X.; Zhao, X. An Opposition-Based Sine Cosine Approach with Local Search for Parameter Estimation of Photovoltaic Models. Energy Convers. Manag. 2019, 195, 927–942. [Google Scholar] [CrossRef]
  5. Rasheduzzaman, M.; Fajri, P.; Kimball, J.; Deken, B. Modeling, Analysis, and Control Design of a Single-Stage Boost Inverter. Energies 2021, 14, 4098. [Google Scholar] [CrossRef]
  6. Zhang, W.; Zheng, Z.; Liu, H. A Novel Droop Control Method to Achieve Maximum Power Output of Photovoltaic for Parallel Inverter System. CSEE J. Power Energy Syst. 2021, 8, 1636–1645. [Google Scholar] [CrossRef]
  7. Mehta, H.K.; Warke, H.; Kukadiya, K.; Panchal, A.K. Accurate Expressions for Single-Diode-Model Solar Cell Parameterization. IEEE J. Photovoltaics 2019, 9, 803–810. [Google Scholar] [CrossRef]
  8. Hejri, M.; Mokhtari, H.; Azizian, M.R.; Ghandhari, M.; Söder, L. On the Parameter Extraction of a Five-Parameter Double-Diode Model of Photovoltaic Cells and Modules. IEEE J. Photovoltaics 2014, 4, 915–923. [Google Scholar] [CrossRef]
  9. Chin, V.J.; Salam, Z. A New Three-Point-Based Approach for the Parameter Extraction of Photovoltaic Cells. Appl. Energy 2019, 237, 519–533. [Google Scholar] [CrossRef]
  10. Akbari, M.A.; Zare, M.; Azizipanah-Abarghooee, R.; Mirjalili, S.; Deriche, M. The Cheetah Optimizer: A Nature-Inspired Metaheuristic Algorithm for Large-Scale Optimization Problems. Sci. Rep. 2022, 12, 10953. [Google Scholar] [CrossRef]
  11. Kirkpatrick, S.; Gelatt, C.D.; Vecchi, M.P. Optimization by Simulated Annealing. Science 1983, 220, 671–680. [Google Scholar] [CrossRef] [PubMed]
  12. Ghasemi, M.; Akbari, M.-A.; Jun, C.; Bateni, S.M.; Zare, M.; Zahedi, A.; Pai, H.-T.; Band, S.S.; Moslehpour, M.; Chau, K.-W. Circulatory System Based Optimization (CSBO): An Expert Multilevel Biologically Inspired Meta-Heuristic Algorithm. Eng. Appl. Comput. Fluid Mech. 2022, 16, 1483–1525. [Google Scholar] [CrossRef]
  13. Holland, J.H. Genetic Algorithms. Sci. Am. 1992, 267, 66–73. [Google Scholar] [CrossRef]
  14. Yao, X.; Liu, Y.; Lin, G. Evolutionary Programming Made Faster. IEEE Trans. Evol. Comput. 1999, 3, 82–102. [Google Scholar]
  15. Koza, J.R. Genetic Programming as a Means for Programming Computers by Natural Selection. Stat. Comput. 1994, 4, 87–112. [Google Scholar] [CrossRef]
  16. Simon, D. Biogeography-Based Optimization. IEEE Trans. Evol. Comput. 2008, 12, 702–713. [Google Scholar] [CrossRef] [Green Version]
  17. Storn, R.; Price, K. Differential Evolution—A Simple and Efficient Heuristic for Global Optimization over Continuous Spaces. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  18. Rashedi, E.; Nezamabadi-Pour, H.; Saryazdi, S. GSA: A Gravitational Search Algorithm. Inf. Sci. 2009, 179, 2232–2248. [Google Scholar] [CrossRef]
  19. Kaveh, A.; Talatahari, S. A Novel Heuristic Optimization Method: Charged System Search. Acta Mech. 2010, 213, 267–289. [Google Scholar] [CrossRef]
  20. Rabanal, P.; Rodríguez, I.; Rubio, F. Using River Formation Dynamics to Design Heuristic Algorithms. In International Conference on Unconventional Computation; Springer: Berlin/Heidelberg, Germany, 2007; pp. 163–177. [Google Scholar]
  21. Erol, O.K.; Eksin, I. A New Optimization Method: Big Bang–Big Crunch. Adv. Eng. Softw. 2006, 37, 106–111. [Google Scholar] [CrossRef]
  22. Boettcher, S.; Percus, A.G. Optimization with Extremal Dynamics. Complexity 2002, 8, 57–62. [Google Scholar] [CrossRef]
  23. Shah-Hosseini, H. Principal Components Analysis by the Galaxy-Based Search Algorithm: A Novel Metaheuristic for Continuous Optimisation. Int. J. Comput. Sci. Eng. 2011, 6, 132–140. [Google Scholar]
  24. Formato, R.A. Central Force Optimization. Prog. Electromagn. Res. 2007, 77, 425–491. [Google Scholar] [CrossRef] [Green Version]
  25. Zhao, W.; Zhang, Z.; Wang, L. Manta Ray Foraging Optimization: An Effective Bio-Inspired Optimizer for Engineering Applications. Eng. Appl. Artif. Intell. 2020, 87, 103300. [Google Scholar] [CrossRef]
  26. Eskandar, H.; Sadollah, A.; Bahreininejad, A.; Hamdi, M. Water Cycle Algorithm—A Novel Metaheuristic Optimization Method for Solving Constrained Engineering Optimization Problems. Comput. Struct. 2012, 110, 151–166. [Google Scholar] [CrossRef]
  27. Hosseini, H.S. Problem Solving by Intelligent Water Drops. In Proceedings of the 2007 IEEE Congress on Evolutionary Computation, Singapore, 25–28 September 2007; pp. 3226–3231. [Google Scholar]
  28. Li, B.; Jiang, W. Chaos Optimization Method and Its Application. Control Theory Appl. 1997, 14, 613–615. [Google Scholar]
  29. Birbil, Ş.İ.; Fang, S.-C. An Electromagnetism-like Mechanism for Global Optimization. J. Glob. Optim. 2003, 25, 263–282. [Google Scholar] [CrossRef]
  30. Alatas, B. ACROA: Artificial Chemical Reaction Optimization Algorithm for Global Optimization. Expert Syst. Appl. 2011, 38, 13170–13180. [Google Scholar] [CrossRef]
  31. Irizarry, R. LARES: An Artificial Chemical Process Approach for Optimization. Evol. Comput. 2004, 12, 435–459. [Google Scholar] [CrossRef]
  32. Abdechiri, M.; Meybodi, M.R.; Bahrami, H. Gases Brownian Motion Optimization: An Algorithm for Optimization (GBMO). Appl. Soft Comput. 2013, 13, 2932–2946. [Google Scholar] [CrossRef]
  33. Kennedy, J.; Eberhart, R. Particle Swarm Optimization. In Proceedings of the ICNN’95—International Conference on Neural Networks, Perth, WA, Australia, 27 November–1 December 1995; Volume 4, pp. 1942–1948. [Google Scholar]
  34. Yang, X.-S.; Deb, S. Cuckoo Search via Lévy Flights. In Proceedings of the 2009 World Congress on Nature & Biologically Inspired Computing (NaBIC), Coimbatore, India, 9–11 December 2009; pp. 210–214. [Google Scholar]
  35. Mirjalili, S. The Ant Lion Optimizer. Adv. Eng. Softw. 2015, 83, 80–98. [Google Scholar] [CrossRef]
  36. Pham, D.T.; Ghanbarzadeh, A.; Koç, E.; Otri, S.; Rahim, S.; Zaidi, M. The Bees Algorithm—A Novel Tool for Complex Optimisation Problems. In Intelligent Production Machines and Systems; Elsevier: Amsterdam, The Netherlands, 2006; pp. 454–459. [Google Scholar]
  37. Eusuff, M.M.; Lansey, K.E. Optimization of Water Distribution Network Design Using the Shuffled Frog Leaping Algorithm. J. Water Resour. Plan. Manag. 2003, 129, 210–225. [Google Scholar] [CrossRef]
  38. Yang, X.; Gandomi, A.H. Bat Algorithm: A Novel Approach for Global Engineering Optimization. Eng. Comput. 2012, 29, 464–483. [Google Scholar] [CrossRef]
  39. Mirjalili, S. Moth-Flame Optimization Algorithm: A Novel Nature-Inspired Heuristic Paradigm. Knowl.-Based Syst. 2015, 89, 228–249. [Google Scholar] [CrossRef]
  40. Panigrahi, B.K.; Ravikumar Pandi, V. Bacterial Foraging Optimisation: Nelder–Mead Hybrid Algorithm for Economic Load Dispatch. IET Gener. Transm. Distrib. 2008, 2, 556. [Google Scholar] [CrossRef]
  41. Gandomi, A.H.; Alavi, A.H. Krill Herd: A New Bio-Inspired Optimization Algorithm. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 4831–4845. [Google Scholar] [CrossRef]
  42. Mirjalili, S.; Lewis, A. The Whale Optimization Algorithm. Adv. Eng. Softw. 2016, 95, 51–67. [Google Scholar] [CrossRef]
  43. Drigo, M. The Ant System: Optimization by a Colony of Cooperating Agents. IEEE Trans. Syst. Man Cybern. B 1996, 26, 1–13. [Google Scholar] [CrossRef] [Green Version]
  44. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey Wolf Optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
  45. Yang, X.-S. Firefly Algorithms for Multimodal Optimization. In International Symposium on Stochastic Algorithms; Springer: Berlin/Heidelberg, Germany, 2009; pp. 169–178. [Google Scholar]
  46. Akay, B.; Karaboga, D. Artificial Bee Colony Algorithm for Large-Scale Problems and Engineering Design Optimization. J. Intell. Manuf. 2012, 23, 1001–1014. [Google Scholar] [CrossRef]
  47. Pan, W.-T. A New Fruit Fly Optimization Algorithm: Taking the Financial Distress Model as an Example. Knowl.-Based Syst. 2012, 26, 69–74. [Google Scholar] [CrossRef]
  48. Krishnanand, K.N.; Ghose, D. Glowworm Swarm Optimization for Simultaneous Capture of Multiple Local Optima of Multimodal Functions. Swarm Intell. 2009, 3, 87–124. [Google Scholar] [CrossRef]
  49. Nadimi-Shahraki, M.H.; Taghian, S.; Mirjalili, S.; Abualigah, L.; Abd Elaziz, M.; Oliva, D. EWOA-OPF: Effective Whale Optimization Algorithm to Solve Optimal Power Flow Problem. Electronics 2021, 10, 2975. [Google Scholar] [CrossRef]
  50. Nadimi-Shahraki, M.H.; Taghian, S.; Mirjalili, S.; Faris, H. MTDE: An Effective Multi-Trial Vector-Based Differential Evolution Algorithm and Its Applications for Engineering Design Problems. Appl. Soft Comput. 2020, 97, 106761. [Google Scholar] [CrossRef]
  51. Nadimi-Shahraki, M.H.; Taghian, S.; Mirjalili, S. An Improved Grey Wolf Optimizer for Solving Engineering Problems. Expert Syst. Appl. 2021, 166, 113917. [Google Scholar] [CrossRef]
  52. Eslami, M.; Shareef, H.; Mohamed, A.; Khajehzadeh, M. Optimal Location of PSS Using Improved PSO with Chaotic Sequence. In Proceedings of the International Conference on Electrical, Control and Computer Engineering 2011 (InECCE), Kuantan, Malaysia, 21–22 June 2011; pp. 253–258. [Google Scholar]
  53. Eslami, M.; Shareef, H.; Mohamed, A.; Khajehzadeh, M. Coordinated Design of PSS and SVC Damping Controller Using CPSO. In Proceedings of the 2011 5th International Power Engineering and Optimization Conference, Shah Alam, Malaysia, 6–7 June 2011; pp. 11–16. [Google Scholar]
  54. Khajehzadeh, M.; Taha, M.R.; Eslami, M. Multi-Objective Optimisation of Retaining Walls Using Hybrid Adaptive Gravitational Search Algorithm. Civ. Eng. Environ. Syst. 2014, 31, 229–242. [Google Scholar] [CrossRef]
  55. Messaoud, R. Ben Extraction of Uncertain Parameters of Single and Double Diode Model of a Photovoltaic Panel Using Salp Swarm Algorithm. Measurement 2020, 154, 107446. [Google Scholar] [CrossRef]
  56. Rao, R. Jaya: A Simple and New Optimization Algorithm for Solving Constrained and Unconstrained Optimization Problems. Int. J. Ind. Eng. Comput. 2016, 7, 19–34. [Google Scholar] [CrossRef]
  57. Yu, K.; Qu, B.; Yue, C.; Ge, S.; Chen, X.; Liang, J. A Performance-Guided JAYA Algorithm for Parameters Identification of Photovoltaic Cell and Module. Appl. Energy 2019, 237, 241–257. [Google Scholar] [CrossRef]
  58. Yu, K.; Liang, J.J.; Qu, B.Y.; Chen, X.; Wang, H. Parameters Identification of Photovoltaic Models Using an Improved JAYA Optimization Algorithm. Energy Convers. Manag. 2017, 150, 742–753. [Google Scholar] [CrossRef]
  59. Patel, S.J.; Panchal, A.K.; Kheraj, V. Extraction of Solar Cell Parameters from a Single Current–Voltage Characteristic Using Teaching Learning Based Optimization Algorithm. Appl. Energy 2014, 119, 384–393. [Google Scholar] [CrossRef]
  60. Niu, Q.; Zhang, H.; Li, K. An Improved TLBO with Elite Strategy for Parameters Identification of PEM Fuel Cell and Solar Cell Models. Int. J. Hydrogen Energy 2014, 39, 3837–3854. [Google Scholar] [CrossRef]
  61. Chen, X.; Yu, K.; Du, W.; Zhao, W.; Liu, G. Parameters Identification of Solar Cell Models Using Generalized Oppositional Teaching Learning Based Optimization. Energy 2016, 99, 170–180. [Google Scholar] [CrossRef]
  62. Yaghoubi, M.; Eslami, M.; Noroozi, M.; Mohammadi, H.; Kamari, O.; Palani, S. Modified Salp Swarm Optimization for Parameter Estimation of Solar PV Models. IEEE Access 2022, 10, 110181–110194. [Google Scholar] [CrossRef]
  63. Ishaque, K.; Salam, Z.; Mekhilef, S.; Shamsudin, A. Parameter Extraction of Solar Photovoltaic Modules Using Penalty-Based Differential Evolution. Appl. Energy 2012, 99, 297–308. [Google Scholar] [CrossRef]
  64. Jiang, L.L.; Maskell, D.L.; Patra, J.C. Parameter Estimation of Solar Cells and Modules Using an Improved Adaptive Differential Evolution Algorithm. Appl. Energy 2013, 112, 185–193. [Google Scholar] [CrossRef]
  65. AlRashidi, M.R.; AlHajri, M.F.; El-Naggar, K.M.; Al-Othman, A.K. A New Estimation Approach for Determining the I-V Characteristics of Solar Cells. Sol. Energy 2011, 85, 1543–1550. [Google Scholar] [CrossRef]
  66. Yousri, D.; Allam, D.; Eteiba, M.B.; Suganthan, P.N. Static and Dynamic Photovoltaic Models’ Parameters Identification Using Chaotic Heterogeneous Comprehensive Learning Particle Swarm Optimizer Variants. Energy Convers. Manag. 2019, 182, 546–563. [Google Scholar] [CrossRef]
  67. Ebrahimi, S.M.; Salahshour, E.; Malekzadeh, M.; Gordillo, F. Parameters Identification of PV Solar Cells and Modules Using Flexible Particle Swarm Optimization Algorithm. Energy 2019, 179, 358–372. [Google Scholar] [CrossRef]
  68. Liang, J.; Ge, S.; Qu, B.; Yu, K.; Liu, F.; Yang, H.; Wei, P.; Li, Z. Classified Perturbation Mutation Based Particle Swarm Optimization Algorithm for Parameters Extraction of Photovoltaic Models. Energy Convers. Manag. 2020, 203, 112138. [Google Scholar] [CrossRef]
  69. Lin, X.; Wu, Y. Parameters Identification of Photovoltaic Models Using Niche-Based Particle Swarm Optimization in Parallel Computing Architecture. Energy 2020, 196, 117054. [Google Scholar] [CrossRef]
  70. Yousri, D.; Thanikanti, S.B.; Allam, D.; Ramachandaramurthy, V.K.; Eteiba, M.B. Fractional Chaotic Ensemble Particle Swarm Optimizer for Identifying the Single, Double, and Three Diode Photovoltaic Models’ Parameters. Energy 2020, 195, 116979. [Google Scholar] [CrossRef]
  71. Nunes, H.G.G.; Silva, P.N.C.; Pombo, J.A.N.; Mariano, S.; Calado, M.R.A. Multiswarm Spiral Leader Particle Swarm Optimisation Algorithm for PV Parameter Identification. Energy Convers. Manag. 2020, 225, 113388. [Google Scholar] [CrossRef]
  72. Premkumar, M.; Jangir, P.; Elavarasan, R.M.; Sowmya, R. Opposition Decided Gradient-Based Optimizer with Balance Analysis and Diversity Maintenance for Parameter Identification of Solar Photovoltaic Models. J. Ambient Intell. Humaniz. Comput. 2021, 1–23. [Google Scholar] [CrossRef]
  73. Gafar, M.; El-Sehiemy, R.A.; Hasanien, H.M.; Abaza, A. Optimal Parameter Estimation of Three Solar Cell Models Using Modified Spotted Hyena Optimization. J. Ambient Intell. Humaniz. Comput. 2022, 1–12. [Google Scholar] [CrossRef]
  74. Eslami, M.; Akbari, E.; Seyed Sadr, S.T.; Ibrahim, B.F. A Novel Hybrid Algorithm Based on Rat Swarm Optimization and Pattern Search for Parameter Extraction of Solar Photovoltaic Models. Energy Sci. Eng. 2022, 10, 2689–2713. [Google Scholar] [CrossRef]
  75. Wolpert, D.H.; Macready, W.G. No Free Lunch Theorems for Optimization. IEEE Trans. Evol. Comput. 1997, 1, 67–82. [Google Scholar] [CrossRef] [Green Version]
  76. Easwarakhanthan, T.; Bottin, J.; Bouhouch, I.; Boutrit, C. Nonlinear Minimization Algorithm for Determining the Solar Cell Parameters with Microcomputers. Int. J. Sol. Energy 1986, 4, 1–12. [Google Scholar] [CrossRef]
  77. Rao, R.V.; Savsani, V.J.; Vakharia, D.P. Teaching–Learning-Based Optimization: A Novel Method for Constrained Mechanical Design Optimization Problems. Comput. Des. 2011, 43, 303–315. [Google Scholar] [CrossRef]
  78. Mirjalili, S.; Gandomi, A.H.; Mirjalili, S.Z.; Saremi, S.; Faris, H.; Mirjalili, S.M. Salp Swarm Algorithm: A Bio-Inspired Optimizer for Engineering Design Problems. Adv. Eng. Softw. 2017, 114, 163–191. [Google Scholar] [CrossRef]
  79. Kaveh, A.; Mahdavi, V.R. Colliding Bodies Optimization: A Novel Meta-Heuristic Method. Comput. Struct. 2014, 139, 18–27. [Google Scholar] [CrossRef]
  80. Hussain, K.; Salleh, M.N.M.; Cheng, S.; Shi, Y. On the Exploration and Exploitation in Popular Swarm-Based Metaheuristic Algorithms. Neural Comput. Appl. 2019, 31, 7665–7683. [Google Scholar] [CrossRef]
  81. Chen, X.; Tianfield, H.; Mei, C.; Du, W.; Liu, G. Biogeography-Based Learning Particle Swarm Optimization. Soft Comput. 2017, 21, 7519–7541. [Google Scholar] [CrossRef]
  82. Liang, J.J.; Qin, A.K.; Suganthan, P.N.; Baskar, S. Comprehensive Learning Particle Swarm Optimizer for Global Optimization of Multimodal Functions. IEEE Trans. Evol. Comput. 2006, 10, 281–295. [Google Scholar] [CrossRef]
  83. Chen, X.; Xu, B.; Mei, C.; Ding, Y.; Li, K. Teaching–Learning–Based Artificial Bee Colony for Solar Photovoltaic Parameter Estimation. Appl. Energy 2018, 212, 1578–1588. [Google Scholar] [CrossRef]
  84. Gong, W.; Cai, Z.; Ling, C.X. DE/BBO: A Hybrid Differential Evolution with Biogeography-Based Optimization for Global Numerical Optimization. Soft Comput. 2010, 15, 645–665. [Google Scholar] [CrossRef] [Green Version]
  85. Chen, X.; Tianfield, H.; Du, W.; Liu, G. Biogeography-Based Optimization with Covariance Matrix Based Migration. Appl. Soft Comput. 2016, 45, 71–85. [Google Scholar] [CrossRef]
Figure 1. A sample of random value over 1000 iterations.
Figure 1. A sample of random value over 1000 iterations.
Mathematics 10 04617 g001
Figure 2. A schematic of position updating in the APLO algorithm.
Figure 2. A schematic of position updating in the APLO algorithm.
Mathematics 10 04617 g002
Figure 3. Flowchart of APLO algorithm.
Figure 3. Flowchart of APLO algorithm.
Mathematics 10 04617 g003
Figure 4. Equivalent circuit of SDM.
Figure 4. Equivalent circuit of SDM.
Mathematics 10 04617 g004
Figure 5. Equivalent circuit of DDM.
Figure 5. Equivalent circuit of DDM.
Mathematics 10 04617 g005
Figure 6. Equivalent circuit of PV module model.
Figure 6. Equivalent circuit of PV module model.
Mathematics 10 04617 g006
Figure 7. Exploration and exploitation of the competitive algorithms through the iterations on the SDM problem.
Figure 7. Exploration and exploitation of the competitive algorithms through the iterations on the SDM problem.
Mathematics 10 04617 g007aMathematics 10 04617 g007b
Figure 8. Diversity in individuals of the competitive algorithms on SDM problem.
Figure 8. Diversity in individuals of the competitive algorithms on SDM problem.
Mathematics 10 04617 g008
Figure 9. Mean values of exploration and exploitation of competitive algorithms on SDM problem.
Figure 9. Mean values of exploration and exploitation of competitive algorithms on SDM problem.
Mathematics 10 04617 g009
Figure 10. Convergence characteristics of the applied algorithms on SDM problem.
Figure 10. Convergence characteristics of the applied algorithms on SDM problem.
Mathematics 10 04617 g010
Figure 11. Box chart of different algorithms on SDM parameter extraction over 30 runs.
Figure 11. Box chart of different algorithms on SDM parameter extraction over 30 runs.
Mathematics 10 04617 g011
Figure 12. Measured and calculated data of the RTC France silicon solar cell based on SDM by APLO.
Figure 12. Measured and calculated data of the RTC France silicon solar cell based on SDM by APLO.
Mathematics 10 04617 g012
Figure 13. Boxplot comparison of different algorithms for DDM.
Figure 13. Boxplot comparison of different algorithms for DDM.
Mathematics 10 04617 g013
Figure 14. Measured and calculated data of the RTC France silicon solar cell based on DDM by APLO.
Figure 14. Measured and calculated data of the RTC France silicon solar cell based on DDM by APLO.
Mathematics 10 04617 g014
Figure 15. Boxplot comparison of different algorithms for PV module.
Figure 15. Boxplot comparison of different algorithms for PV module.
Mathematics 10 04617 g015
Figure 16. Comparison between the measured and calculated data yielded by APLO for the PV module model based on Photo Watt-PWP 201.
Figure 16. Comparison between the measured and calculated data yielded by APLO for the PV module model based on Photo Watt-PWP 201.
Mathematics 10 04617 g016
Figure 17. Convergence curves of the competitor algorithms for SDM.
Figure 17. Convergence curves of the competitor algorithms for SDM.
Mathematics 10 04617 g017
Figure 18. Convergence curves of the competitor algorithms for DDM.
Figure 18. Convergence curves of the competitor algorithms for DDM.
Mathematics 10 04617 g018
Figure 19. Convergence curves of the competitor algorithms for the PV module model.
Figure 19. Convergence curves of the competitor algorithms for the PV module model.
Mathematics 10 04617 g019
Figure 20. CPU time comparison in solving parameter extraction of SDM by different algorithms over 30 runs.
Figure 20. CPU time comparison in solving parameter extraction of SDM by different algorithms over 30 runs.
Mathematics 10 04617 g020
Figure 21. CPU time comparison in solving parameter extraction of DDM by different algorithms over 30 runs.
Figure 21. CPU time comparison in solving parameter extraction of DDM by different algorithms over 30 runs.
Mathematics 10 04617 g021
Figure 22. CPU time comparison in solving parameter extraction of PV module model by different algorithms over 30 runs.
Figure 22. CPU time comparison in solving parameter extraction of PV module model by different algorithms over 30 runs.
Mathematics 10 04617 g022
Table 1. The classification of population-based metaheuristic algorithms.
Table 1. The classification of population-based metaheuristic algorithms.
CategoryMost Popular Algorithms and Abbreviations
Evolutionary-basedGenetic algorithm (GA) [13]
Evolutionary programming (EP) [14]
Genetic programming (GP) [15]
Biogeography-based optimizer (BBO) [16]
Differential evolution (DE) [17]
Evolutionary strategy (ES)
Physics-basedGravitational Ssearch Algorithm (GSA) [18]
Charged system search (CSS) [19]
River formation dynamics algorithm (RFDA) [20]
Big bang–big crunch (BB-BC) [21]
Extremal optimization (EO) [22]
Galaxy-based search algorithm (GBSA) [23]
Central force optimization (CFO) [24]
Ray optimization (RO) [25]
Water cycle algorithm (WCA) [26]
Intelligent water drops (IWD) [27]
Chaos optimization algorithm (COA) [28]
Electromagnetism-like mechanism (EM) [29]
Chemistry-basedArtificial chemical reaction optimization algorithm (ACROA) [30]
Artificial chemical process (ACP) [31]
Gases Brownian motion optimization (GBMO) [32]
Swarm-basedParticle swarm optimization (PSO) [33]
Cuckoo search (CS) [34]
Ant lion optimizer (ALO) [35]
Bees algorithm (BA) [36]
Shuffled frog-leaping algorithm (SFLA) [37]
Bat algorithm (BA) [38]
Moth–flame optimization (MFO) [39]
Bacterial foraging algorithm (BFA) [40]
Krill herd (KH) [41]
Whale optimization algorithm (WOA) [42]
Ant colony algorithms (ACO) [43]
Grey wolf optimizer (GWO) [44]
Firefly algorithm (FA) [45]
Artificial bee colony (ABC) [46]
Fruit fly optimization algorithm (FOA) [47]
Glowworm swarm optimization (GSO) [48]
Table 2. Parameters’ upper and lower ranges.
Table 2. Parameters’ upper and lower ranges.
ModelParameters’ Limits
I ¯ p h   ( A ) I ¯ p h   ( A ) I ¯ s d   ( μ A ) I ¯ s d   ( μ A ) m ¯ m ¯ R ¯ s   ( Ω ) R ¯ s   ( Ω ) R ¯ s h   ( Ω ) R ¯ s h   ( Ω )
SDM01011200.50100
DDM01011200.50100
PV module020501500202000
Table 3. Statistical analysis of population size on the performance of APLO for different PV models (the significant values are bolded).
Table 3. Statistical analysis of population size on the performance of APLO for different PV models (the significant values are bolded).
ProblemMeasurenpop = 10npop = 20npop = 30npop = 40npop = 50npop = 60
SDMMin9.860219 × 10−49.860219 × 10−49.860219 × 10−49.860219 × 10−49.860219 × 10−49.860219 × 10−4
Mean9.860219 × 10−49.860219 × 10−49.860219 × 10−49.860219 × 10−49.860219 × 10−49.860219 × 10−4
Max9.860219 × 10−49.860219 × 10−49.860219 × 10−49.860219 × 10−49.860219 × 10−49.860219 × 10−4
SD4.905627 × 10−124.149942 × 10−168.666882 × 10−171.126621 × 10−169.397154 × 10−178.574883 × 10−17
Mean rank in Freidman5.966673.616672.733332.433333.13.15
Sum rank Freidman179108.582739394.5
DDMMin9.829117 × 10−49.826894 × 10−49.824958 × 10−49.824849 × 10−49.827377 × 10−49.857074 × 10−4
Mean1.008130 × 10−31.006467 × 10−31.003314 × 10−31.011153 × 10−31.013609 × 10−31.022566 × 10−3
Max1.191914 × 10−31.131532 × 10−31.191268 × 10−31.402166 × 10−31.338409 × 10−31.189575 × 10−3
SD4.099853 × 10−53.637408 × 10−54.636394 × 10−57.838701 × 10−57.170668 × 10−55.765975 × 10−5
Mean rank in Freidman3.933333.433333.166673.13.233334.13333
Sum rank Freidman118103959397124
PV moduleMin9.825837 × 10−49.825044 × 10−49.825115 × 10−49.827521 × 10−49.832545 × 10−49.825992 × 10−4
Mean4.845660 × 10−39.993596 × 10−49.966059 × 10−41.001461 × 10−31.011433 × 10−31.029418 × 10−3
Max1.149226 × 10−11.060000 × 10−31.051849 × 10−31.132161 × 10−31.225400 × 10−31.445858 × 10−3
SD2.079038 × 10−22.217082 × 10−51.878370 × 10−53.269111 × 10−55.876697 × 10−59.490092 × 10−5
Mean rank in Freidman4.666673.233333.133333.33.13.56667
Sum rank Freidman14097949993107
Sum rankMean rank in Freidman14.5666710.283339.033338.833339.4333310.85
Sum rank Freidman437308.5271265283325.5
Table 4. The statistical results of RMSE for SDM (the significant values are bolded).
Table 4. The statistical results of RMSE for SDM (the significant values are bolded).
AlgorithmMinMeanMaxSDSignificance
APLO9.860218778914× 10−49.860218778916× 10−49.860218778922× 10−41.599419351161× 10−16
CBO1.148573645374 × 10−37.208263932384 × 10−38.576997450439 × 10−31.645251727229 × 10−3
DE5.548414731899 × 10−36.915313985034 × 10−37.850952678115 × 10−36.541614750785 × 10−4
GA3.817138390879 × 10−36.837349221993 × 10−28.169353790970 × 10−21.780787047086 × 10−2
GWO1.158182213682 × 10−22.158244241773 × 10−12.228964974584 × 10−13.857526966905 × 10−2
JAYA1.782507253178 × 10−34.846317911802 × 10−39.666422102947 × 10−31.496853195820 × 10−3
PSO9.869017992242 × 10−42.588954874116 × 10−34.826880699904 × 10−31.103903013051 × 10−3
SSA2.228658847342 × 10−12.230928389927 × 10−12.239170090380 × 10−12.800981095230 × 10−4
TLBO9.887536713543 × 10−43.863042069493 × 10−38.677457995748 × 10−32.079581847494 × 10−3
† Indicates APLO has a significant advantage over its competitor when Wilcoxon’s rank sum test is performed at 5% confidence.
Table 5. The best-identified parameters for SDM using the competitor algorithms (the significant values are bolded).
Table 5. The best-identified parameters for SDM using the competitor algorithms (the significant values are bolded).
Algorithm I p h   ( A ) I s d     ( μ A ) R s h   ( Ω ) R s   ( Ω ) m RMSE
APLO0.76077550.323020853.718524000.03637711.48118559.8602188 × 10−4
CBO0.76063650.439183763.846879860.03512731.51276691.1485736 × 10−3
DE0.76338403.5631230100.00000000.02390901.77080805.5484140 × 10−3
GA0.76179081.823216199.999927360.02805091.67920313.8171384 × 10−3
GWO0.76668610.813449414.695651250.02737661.58712881.1581822 × 10−2
JAYA0.75930100.5978732100.00000000.03410351.54568571.7825073 × 10−3
PSO0.76076640.330157954.310319010.03628961.48338919.8690180 × 10−4
SSA0.83619820.00000001.1550933330.00000002.00000002.2286588 × 10−1
TLBO0.76070490.331494554.987919230.03627151.48377859.8875367 × 10−4
Table 6. Statistical results of algorithms for parameter extraction of DDM (the significant values are bolded).
Table 6. Statistical results of algorithms for parameter extraction of DDM (the significant values are bolded).
AlgorithmMinMeanMaxSDSignificance
APLO9.830657938188 × 10−41.019874828873 × 10−31.342343716280 × 10−37.797062995961 × 10−5
CBO6.445682677224 × 10−37.474969722257 × 10−38.275185158026 × 10−35.353607495130 × 10−4
DE6.761178531429 × 10−38.130913835822 × 10−39.217501072530 × 10−36.044793459888 × 10−4
GA1.048914576651 × 10−34.277390976996 × 10−29.079203590158 × 10−23.753851916526 × 10−2
GWO7.723436258756 × 10−31.804748010930 × 10−12.228793281927 × 10−18.624502969932 × 10−2
JAYA2.798207817458 × 10−35.479908398673 × 10−39.234026428131 × 10−31.396786107898 × 10−3
PSO9.848728345415 × 10−42.231732852337 × 10−33.903661545904 × 10−39.893579067548 × 10−4
SSA2.228624210303 × 10−12.234264679038 × 10−12.243772831227 × 10−15.008454232877 × 10−4
TLBO1.006361630316 × 10−36.071268460472 × 10−31.954518527698 × 10−24.166155543866 × 10−3
† Indicates APLO has a significant advantage over its competitor when Wilcoxon’s rank sum test is performed at 5% confidence.
Table 7. The best-identified parameters for DDM using the competitor algorithms (the significant values are bolded).
Table 7. The best-identified parameters for DDM using the competitor algorithms (the significant values are bolded).
Algorithm I p h   ( A ) I s d , 1     ( μ A ) I s d , 2   ( μ A ) R s   ( Ω ) R s h   ( Ω ) m 1 m 2 RMSE
APLO0.76077570.46979110.2609657790.036580754.904627801.99996411.46314859.8306579382 × 10−4
CBO0.76388700.24351628.0406156020.023444899.982253741.53646642.00000006.4456826772 × 10−3
DE0.76490501.70980704.766027700.023261099.999990001.72903701.98675106.7611785310 × 10−3
GA0.76057980.38845310.000000000.035644860.423842291.49997941.47541311.0489145767 × 10−3
GWO0.77063570.15719630.000000000.041898730.382969581.41074971.67473187.7234362588 × 10−3
JAYA0.75731800.05568820.3152246340.0356605100.00000001.47772751.49795132.7982078175 × 10−3
PSO0.76077950.29575040.1452544620.036461753.997354151.47381471.92989069.8487283454 × 10−4
SSA0.83772290.00000000.0000000000.00000001.1451146811.00000002.00000002.2286242103 × 10−1
TLBO0.76068200.35433790.3328655550.036086156.6489828861.49068291.00000001.0063616303 × 10−3
Table 8. Statistical results of algorithms for parameter extraction of PV module (the significant values are bolded).
Table 8. Statistical results of algorithms for parameter extraction of PV module (the significant values are bolded).
AlgorithmMinMeanMaxSDSignificance
APLO2.42507486809 × 10−32.42507486810 × 10−32.42507486810 × 10−35.96208271747 × 10−17
CBO2.59323307598 × 10−39.58546210390 × 10−31.44726212656 × 10−21.62949321579 × 10−3
DE6.67237575292 × 10−37.74335013025 × 10−39.59388984237 × 10−37.13076666106 × 10−4
GA5.22983393188 × 10−31.72642387255 × 10−13.51706994835 × 10−11.17731216363 × 10−1
GWO9.30534798242 × 10−33.03524536723 × 10−21.09058464220 × 10−12.40958661843 × 10−2
JAYA3.42058453137 × 10−31.63022106280 × 10−27.64341081931 × 10−22.27085786595 × 10−2
PSO2.47057268706 × 10−35.06269373279 × 10−36.76745897635 × 10−31.11580692510 × 10−3
SSA5.36253332910 × 10−21.45064970492 × 10−12.75832201568 × 10−18.99178958837 × 10−2
TLBO2.81460265874 × 10−34.17330355891 × 10−38.12934998713 × 10−31.16062006801 × 10−3
† Indicates APLO has a significant advantage over its competitor when Wilcoxon’s rank sum test is performed at 5% confidence.
Table 9. The best-identified parameters for PV module using the competitor algorithms (the significant values are bolded).
Table 9. The best-identified parameters for PV module using the competitor algorithms (the significant values are bolded).
Algorithm I p h I s d R s h R s m RMSE
APLO1.03051433.4822628927.27728450.03336861.35119162.4250749 × 10−3
CBO1.02870094.8445513942.71957030.03238661.38721842.5932331 × 10−3
DE1.030615622.16902381999.999590.02649981.58264586.6723758 × 10−3
GA1.028403615.19658371427.988760.02817021.52906815.2298339 × 10−3
GWO1.032952714.2377500798.6394840.02711841.51911829.3053480 × 10−3
JAYA1.02478857.900936722000.000000.03094921.44417783.4205845 × 10−3
PSO1.02944894.1160674833.98298670.03286841.36918682.4705727 × 10−3
SSA1.137514650.00000001.641100860.00252541.74592205.3625333 × 10−2
TLBO1.02887045.5023333847.48944270.03185491.40163052.8146027 × 10−3
Table 10. Wilcoxon ranks test results for SDM.
Table 10. Wilcoxon ranks test results for SDM.
Algorithm IAlgorithm IImena_NRsmean_PRssum_NRssum_PRsN_NRsN_PRsp-Value
APLOCBO15.5NaN46503001.8627 × 10−9
APLODE15.5NaN46503001.8627 × 10−9
APLOGA15.5NaN46503001.8627 × 10−9
APLOGWO15.5NaN46503001.8627 × 10−9
APLOJAYA15.5NaN46503001.8627 × 10−9
APLOPSO15.5NaN46503001.8627 × 10−9
APLOSSA15.5NaN46503001.8627 × 10−9
APLOTLBO15.5NaN46503001.8627 × 10−9
CBODE17.1151203457231.9660 × 10−2
CBOGA16.51.546232829.3132 × 10−9
CBOGWO15.5NaN46503001.8627 × 10−9
CBOJAYA7.017.2354305257.9945 × 10−6
CBOPSO2.016.544612281.3039 × 10−8
CBOSSA15.5NaN46503001.8627 × 10−9
CBOTLBO9.016.0184472284.7125 × 10−7
DEGA16.51.546232829.3132 × 10−9
DEGWO15.5NaN46503001.8627 × 10−9
DEJAYA13.515.6274382282.3488 × 10−6
DEPSONaN15.504650301.8627 × 10−9
DESSA15.5NaN46503001.8627 × 10−9
DETLBO3.817.3154504262.5518 × 10−7
GAGWO16.0146412913.7253 × 10−9
GAJAYA1.016.014641293.7253 × 10−9
GAPSONaN15.504650301.8627 × 10−9
GASSA15.5NaN46503001.8627 × 10−9
GATLBONaN15.504650301.8627 × 10−9
GWOJAYANaN15.504650301.8627 × 10−9
GWOPSONaN15.504650301.8627 × 10−9
GWOSSA15.5NaN46503001.8627 × 10−9
GWOTLBONaN15.504650301.8627 × 10−9
JAYAPSO2.516.454602281.8627 × 10−8
JAYASSA15.5NaN46503001.8627 × 10−9
JAYATLBO17.414.91223437232.2100 × 10−2
PSOSSA15.5NaN46503001.8627 × 10−9
PSOTLBO17.711.235311220101.2050 × 10−2
SSATLBONaN15.504650301.8627 × 10−9
Table 11. Wilcoxon ranks test results for DDM.
Table 11. Wilcoxon ranks test results for DDM.
Algorithm IAlgorithm IImena_NRsmean_PRssum_NRssum_PRsN_NRsN_PRsp-Value
APLOCBO15.5NaN46503001.8627 × 10−9
APLODE15.5NaN46503001.8627 × 10−9
APLOGA15.5NaN46503001.8627 × 10−9
APLOGWO15.5NaN46503001.8627 × 10−9
APLOJAYA15.5NaN46503001.8627 × 10−9
APLOPSO17.33.75450152642.5518 × 10−7
APLOSSA15.5NaN46503001.8627 × 10−9
APLOTLBO16.0146412913.7253 × 10−9
CBODE17.09.5408572461.2334 × 10−4
CBOGA18.893759020102.5600 × 10−3
CBOGWO16.0146412913.7253 × 10−9
CBOJAYA10.016.1304353273.7905 × 10−6
CBOPSONaN15.504650301.8627 × 10−9
CBOSSA15.5NaN46503001.8627 × 10−9
CBOTLBO13.216.51193469211.8530 × 10−2
DEGA20.88.535411117131.1300 × 10−2
DEGWO16.0146412913.7253 × 10−9
DEJAYA4.516.394562286.1467 × 10−8
DEPSONaN15.504650301.8627 × 10−9
DESSA15.5NaN46503001.8627 × 10−9
DETLBO10.917.5983679214.6600 × 10−3
GAGWO16.25.5454112821.0245 × 10−7
GAJAYA11.316.0344313276.9179 × 10−6
GAPSO2.016.024631295.5879 × 10−9
GASSA15.5NaN46503001.8627 × 10−9
GATLBO8.717.2524136247.0568 × 10−5
GWOJAYANaN15.504650301.8627 × 10−9
GWOPSONaN15.504650301.8627 × 10−9
GWOSSA16.0146412913.7253 × 10−9
GWOTLBO2.016.024631295.5879 × 10−9
JAYAPSONaN15.504650301.8627 × 10−9
JAYASSA15.5NaN46503001.8627 × 10−9
JAYATLBO16.514.524721815157.7657 × 10−1
PSOSSA15.5NaN46503001.8627 × 10−9
PSOTLBO16.87437282642.7623 × 10−6
SSATLBONaN15.504650301.8627 × 10−9
Table 12. Wilcoxon ranks test results for the PV module model.
Table 12. Wilcoxon ranks test results for the PV module model.
Algorithm IAlgorithm IImena_NRsmean_PRssum_NRssum_PRsN_NRsN_PRsp-Value
APLOCBO15.5NaN46503001.86265 × 10−9
APLODE15.5NaN46503001.86265 × 10−9
APLOGA15.5NaN46503001.86265 × 10−9
APLOGWO15.5NaN46503001.86265 × 10−9
APLOJAYA15.5NaN46503001.86265 × 10−9
APLOPSO15.5NaN46503001.86265 × 10−9
APLOSSA15.5NaN46503001.86265 × 10−9
APLOTLBO15.5NaN46503001.86265 × 10−9
CBODE15.515.5314342284.42192 × 10−6
CBOGA16.214295.5454112821.02446 × 10−7
CBOGWO15.5NaN46503001.86265 × 10−9
CBOJAYA23.71429131662997230.17719
CBOPSO815.7586284571294.65661 × 10−8
CBOSSA15.5NaN46503001.86265 × 10−9
CBOTLBO215.9655224631295.58794 × 10−9
DEGA16.51.546232829.31323 × 10−9
DEGWO15.5NaN46503001.86265 × 10−9
DEJAYA20.213.1520226310200.54253
DEPSONaN15.504650301.86265 × 10−9
DESSA15.5NaN46503001.86265 × 10−9
DETLBO11614641293.72529 × 10−9
GAGWO4.8571418.73913344317236.91786 × 10−6
GAJAYA7.516.07143154502282.55182 × 10−7
GAPSO11614641293.72529 × 10−9
GASSA1219.519227316140.41613
GATLBONaN15.504650301.86265 × 10−9
GWOJAYA20.3333314.291671223436240.0221
GWOPSONaN15.504650301.86265 × 10−9
GWOSSA15.93103346232919.31323 × 10−9
GWOTLBONaN15.504650301.86265 × 10−9
JAYAPSO6.8333317.66667414246241.82446 × 10−5
JAYASSA16.51.546232829.31323 × 10−9
JAYATLBO816.33333244413271.41934 × 10−6
PSOSSA15.5NaN46503001.86265 × 10−9
PSOTLBO12.2516.68182983678220.00466
SSATLBONaN15.504650301.86265 × 10−9
Table 13. Friedman ranks test results for three models (the significant values are bolded).
Table 13. Friedman ranks test results for three models (the significant values are bolded).
AlgorithmSDMDDMPV Module Model
Mean RankSum RankMean RankSum RankMean RankSum Rank
APLO1.0000301.1667351.000030
CBO5.53331665.16671555.7000171
DE5.10001535.90001774.7667143
GA6.86672065.96671798.0000240
GWO7.96672397.86672367.0667212
JAYA3.96671193.76671134.6000138
PSO2.4667742.0333613.000090
SSA9.00002708.96672698.4333253
TLBO3.1000934.16671252.433373
Table 14. Advantages and disadvantages of applied algorithms on optimal parameter extraction of PV cell and module.
Table 14. Advantages and disadvantages of applied algorithms on optimal parameter extraction of PV cell and module.
AlgorithmGlobal Search AbilityConvergence SpeedLocal Entrapment ProbabilityExplorative/ExploitativeDiversityExploration-Exploitation BalanceComputational Complexity
APLOHighMediumLowExploitativeAdequateGoodLow
CBOLowHighMediumExploitativeInadequateMediumLow
DELowHighHighExploitativeInadequateMediumHigh
GALowHighHighExploitativeInadequateWeakHigh
GWOLowHighHighExplorativeHighWeakLow
JAYAMediumLowMediumExploitativeAdequateMediumLow
PSOMediumMediumMediumExploitativeAdequateMediumHigh
SSALowHighHighExplorativeHighWeakMedium
TLBOHighLowLowExploitativeAdequateGoodMedium
Table 15. Statistic results of the APLO and state-of-the-art algorithms for three PV models (the significant values are bolded).
Table 15. Statistic results of the APLO and state-of-the-art algorithms for three PV models (the significant values are bolded).
SDM
AlgorithmPGJAYAIJAYASTLBOGOTLBOTLABCMSSA
Min9.8602 × 10−49.8603 × 10−49.8602 × 10−49.8856 × 10−49.8602 × 10−49.86 × 10−4
Mean9.8602 × 10−49.9204 × 10−49.8607 × 10−41.0450 × 10−39.9417 × 10−49.86 × 10−4
Max9.8603 × 10−41.0622 × 10−39.8655 × 10−41.2067 × 10−31.0308 × 10−39.87 × 10−4
SD1.4485 × 10−91.4033 × 10−51.8602 × 10−55.0218 × 10−51.1896 × 10−53.01 × 10−7
AlgorithmCLPSOBLPSODE/BBOCMM-DE/BBOAPLOhARS-PS
Min9.9633 × 10−41.0272 × 10−39.9922 × 10−49.8605 × 10−49.8602 × 10−49.84 × 10−4
Mean1.0581 × 10−31.3139 × 10−31.2948 × 10−31.0486 × 10−39.8602 × 10−49.85 × 10−4
Max1.3196 × 10−31.7928 × 10−32.2258 × 10−31.3475 × 10−39.8602 × 10−49.87 × 10−4
SD7.4854 × 10−52.1166 × 10−42.5074 × 10−48.1679 × 10−51.5994 × 10−163.01 × 10−7
DDM
AlgorithmPGJAYAIJAYASTLBOGOTLBOTLABCMSSA
Min9.8263 × 10−49.8293 × 10−49.8252 × 10−49.8742 × 10−41.0012 × 10−39.83 × 10−4
Mean9.8582 × 10−41.0269 × 10−31.0585 × 10−31.1475 × 10−31.2116 × 10−39.94 × 10−4
Max9.9499 × 10−41.4055 × 10−32.4480 × 10−31.3947 × 10−31.9826 × 10−39.99 × 10−4
SD2.5375 × 10−69.8325 × 10−52.8978 × 10−41.1330 × 10−42.1100 × 10−41.49 × 10−6
AlgorithmCLPSOBLPSODE/BBOCMM-DE/BBOAPLOhARS-PS
Min9.9894 × 10−41.0628 × 10−31.0255 × 10−31.0088 × 10−39.8307 × 10−49.82 × 10−4
Mean1.1458 × 10−31.4821 × 10−31.5571 × 10−31.5487 × 10−31.0199 × 10−39.84 × 10−4
Max1.5494 × 10−31.7411 × 10−32.4042 × 10−32.0589 × 10−31.3423 × 10−39.87 × 10−4
SD1.4367 × 10−41.7789 × 10−43.6297 × 10−42.9413 × 10−47.7971 × 10−51.45 × 10−7
PV Module Model
AlgorithmPGJAYAIJAYASTLBOGOTLBOTLABCMSSA
Min2.425075 × 10−32.425129 × 10−32.425075 × 10−32.426583 × 10−32.425075 × 10−32.42 × 10−3
Mean2.425144 × 10−32.428855 × 10−32.055293 × 10−22.475386 × 10−32.425464 × 10−32.54 × 10−3
Max2.426764 × 10−32.439269 × 10−32.742508 × 10−12.563849 × 10−32.428731 × 10−32.78 × 10−3
SD3.071420 × 10−73.775523 × 10−66.896273 × 10−22.938836 × 10−58.746462 × 10−71.75 × 10−5
AlgorithmCLPSOBLPSODE/BBOCMM-DE/BBOAPLOhARS-PS
Min2.428064 × 10−32.425236 × 10−32.428255 × 10−32.425075 × 10−32.425075 × 10−32.42 × 10−3
Mean2.454903 × 10−32.437873 × 10−32.461623 × 10−32.425175 × 10−32.425075 × 10−32.43 × 10−3
Max2.543269 × 10−32.488348 × 10−32.525560 × 10−32.426796 × 10−32.425075 × 10−32.50 × 10−3
SD2.580951 × 10−51.372409 × 10−52.925123 × 10−53.554783 × 10−75.962083 × 10−171.38 × 10−5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alanazi, M.; Alanazi, A.; Almadhor, A.; Rauf, H.T. Photovoltaic Models’ Parameter Extraction Using New Artificial Parameterless Optimization Algorithm. Mathematics 2022, 10, 4617. https://doi.org/10.3390/math10234617

AMA Style

Alanazi M, Alanazi A, Almadhor A, Rauf HT. Photovoltaic Models’ Parameter Extraction Using New Artificial Parameterless Optimization Algorithm. Mathematics. 2022; 10(23):4617. https://doi.org/10.3390/math10234617

Chicago/Turabian Style

Alanazi, Mohana, Abdulaziz Alanazi, Ahmad Almadhor, and Hafiz Tayyab Rauf. 2022. "Photovoltaic Models’ Parameter Extraction Using New Artificial Parameterless Optimization Algorithm" Mathematics 10, no. 23: 4617. https://doi.org/10.3390/math10234617

APA Style

Alanazi, M., Alanazi, A., Almadhor, A., & Rauf, H. T. (2022). Photovoltaic Models’ Parameter Extraction Using New Artificial Parameterless Optimization Algorithm. Mathematics, 10(23), 4617. https://doi.org/10.3390/math10234617

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop