Next Article in Journal
From Constant to Rough: A Survey of Continuous Volatility Modeling
Previous Article in Journal
Dynamic Behavior of a Stochastic Avian Influenza Model with Two Strains of Zoonotic Virus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Modified Gradient Search Rule Based on the Quasi-Newton Method and a New Local Search Technique to Improve the Gradient-Based Algorithm: Solar Photovoltaic Parameter Extraction

by
Bushra Shakir Mahmood
1,
Nazar K. Hussein
1,
Mansourah Aljohani
2 and
Mohammed Qaraad
3,4,*
1
Department of Mathematics, College of Computer Sciences and Mathematics, Tikrit University, Tikrit 34001, Iraq
2
College of Computer Science and Engineering, Taibah University, Yanbu 46421, Saudi Arabia
3
Computer Science Department, Faculty of Science, Abdelmalek Essaadi University, Tetouan 93000, Morocco
4
Department of Computer Science, Faculty of Science, Amran University, Amran 9677, Yemen
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(19), 4200; https://doi.org/10.3390/math11194200
Submission received: 28 August 2023 / Revised: 26 September 2023 / Accepted: 27 September 2023 / Published: 8 October 2023

Abstract

:
Harnessing solar energy efficiently via photovoltaic (PV) technology is pivotal for future sustainable energy. Accurate modeling of PV cells entails an optimization problem due to the multimodal and nonlinear characteristics of the cells. This study introduces the Multi-strategy Gradient-Based Algorithm (MAGBO) for the precise parameter estimation of solar PV systems. MAGBO incorporates a modified gradient search rule (MGSR) inspired by the quasi-Newton approach, a novel refresh operator (NRO) for improved solution quality, and a crossover mechanism balancing exploration and exploitation. Validated through CEC2021 test functions, MAGBO excelled in global optimization. To further validate and underscore the reliability of MAGBO, we utilized data from the PVM 752 GaAs thin-film cell and the STP6-40/36 module. The simulation parameters were discerned using 44 I-V pairs from the PVM 752 cell and diverse data from the STP6-40/36 module tested under different conditions. Consistency between simulated and observed I-V and P-V curves for the STM6-40/36 and PVM 752 models validated MAGBO’s accuracy. In application, MAGBO attained an RMSE of 9.8 × 10−4 for double-diode and single-diode modules. For Photowatt-PWP, STM6-40/36, and PVM 752 models, RMSEs were 2.4 × 10−3, 1.7 × 10−3, and 1.7 × 10−3, respectively. Against prevalent methods, MAGBO exhibited unparalleled precision and reliability, advocating its superior utility for intricate PV data analysis.

1. Introduction

In recent decades, interest in renewable energy sources has become very large due to environmental concerns resulting from fossil fuel pollution, global warming, solid waste, the dangers of burning coal, and the instability of fossil fuel prices. Many alternative energy sources have been found to avoid this, such as wind, water, solar energy, and others. Moreover, solar energy is considered one of the most important energy sources for generating electricity due to its high availability, safe and pollution-free manufacturing methods, and noise-free manufacturing methods. In addition, the price of solar cells has decreased in recent years.
The popularity of solar photovoltaic (PV) systems has increased nowadays because of their use of solar cells [1]. These technologies’ ability to convert solar energy into electricity is impressive. In order to regulate and optimize PV systems, it is crucial to build a precise model that incorporates observed current-voltage data, as this enables the evaluation of the PV array’s actual behavior. Several mathematical models have been developed to shed light on the performance and nonlinear features of PV systems [2], with single-diode and double-diode setups being the most common and commonly used.
The reliability of the parameters is intrinsically tied to the accuracy of the photovoltaic (PV) models. However, environmental variables such as temperature, global irradiation, material flaws, and shifting operational conditions make it difficult to consistently achieve these requirements. Therefore, accurate models that outline the complex relationship between current and voltage are required for the optimization and control of PV systems [3]. In recent years, numerous parameter identification approaches have been developed to meet this pressing requirement [4].
Recently, numerous strategies for finding the parameters of photovoltaic (PV) models have been put forward as potential solutions. They can generally be divided into the following three categories: (1) analytical methods, which are mainly based on the short circuit point, open circuit point, and maximum power point; these methods are fast, simple, and unique, but they are not accurate [5,6,7,8,9,10,11]. (2) Deterministic methods. These methods are sensitive to the initial values and require a level of calculus and convexity. However, they are still not precise enough in the same manner as Newton’s approach [12], the Newton–Raphson approach [13], or the method of the nonlinear algorithm [14]. (3) Meta-heuristic methods are a promising and powerful solution to extract the parameters of the photovoltaic model, and because they are easy to implement and most of them are inspired by the phenomena of nature, they do not need convexity conditions or the derivation of the objective function. Despite developing a set of optimization algorithms, they fail to provide satisfactory results, so meta-heuristic algorithms have been used to solve difficult and complex problems [15].
Consequently, many metaheuristic algorithms have been developed and applied in many fields and to extract the parameters of photoelectric models. The parameters of the diode model of a polycrystalline solar module cell were extracted using the moth optimization algorithm [16], the genetic algorithm (GA), which was derived from Darwin’s development theory [17], particle swarm optimization (PSO) [18], simulated annealing algorithm (SA) [19], artificial bee colony (ABC) [20], anarchic asexual reproduction (CARO) [21], multiple learning backtracking search (MLBSA) [22], cuckoo search algorithm (CS) [23], biogeography optimization (BBO) [24], MBO algorithm [25], sunflower optimization algorithm (SFO) [26], coyote optimization algorithm (COA) [27], Bacterial Foraging Optimization (BFO) [28], Harris haws optimization (HHO) [29], and modified flower algorithm (MFA) [30].
As a direct consequence of this, efforts are currently being undertaken to create more effective optimization algorithms to calculate the parameters of solar photovoltaic cells [31,32]. In the next parts, an in-depth discussion on optimizing photovoltaic systems will be held, focusing on previously developed and evaluated methods.
The algorithms discussed earlier have been applied to various topics, and their results have been found to be satisfactory. However, according to the No Free Lunch theorem [33], there is not yet a superior method of optimization in all circumstances. Furthermore, many optimization methods have benefits and drawbacks that are comparable to one another [34,35]. These approaches, for instance, have a substantial benefit in that they can handle a diverse range of systems with nonlinear fitness functions and restrictions [36]. This is just one example. The primary drawbacks include difficulty adjusting their parameters, the potential for early convergence, the inability to find global optimal solutions, the lack of diversity, and the instability in balancing exploration and exploitation in the solution space [37,38]. These drawbacks can be overcome, however, by using alternative approaches. These methods, if developed with care, have the potential to provide extraordinary performance and versatility for solving difficult optimization problems in the real world [39,40].
As per the No Free Lunch theory, existing optimization algorithms have limitations and weaknesses [41,42,43,44] and may not be suitable for all types of problems. Thus, there is still potential for improvement in algorithm efficiency. Recent studies have focused on developing new ways to improve fundamental optimization procedures by integrating a variety of initial design strategies [45,46,47], hybridizing optimization procedures [48,49], and altering search patterns [50,51,52,53,54]. The Gradient-based Algorithm (GBO), a population-based optimization technique, is an exciting new discovery made very recently. This technique has demonstrated that it is effective and precise in its optimization outcomes. During the phase of exploration and exploitation, GBO uses the information that is available to find optimal or near-optimal solutions by randomly selecting two solutions from the population to predict future relocation and direction based on the best search options. This helps GBO find the best possible options. However, it has been discovered that the new optimization method of GBO has significant drawbacks in high-dimensional and multimodal problems. These drawbacks include getting stuck in local solutions as the size of the search space increases and having poor convergence performance in difficult problems [55,56].
Additionally, the characteristics of the objective function have the potential to have a negative impact on the performance of the algorithm in specific contexts. Particularly when dealing with issues that include many modes of communication, the search method may be hampered. It is well known that the objective function for Solar Photovoltaic Parameter Estimation problems is multimodal and nonlinear with multiple local minimums. Because of this, it is challenging to develop new robust optimization techniques that can produce more accurate and faster convergence results.
This study aims to find the best values for the parameters of Photovoltaic systems. To achieve this, this study proposes a modified gradient search rule (MGSR) based on the quasi-Newton method, derives the MGSR factor that controls vector motion to improve its local and global capabilities, and enhances exploration to improve the search in the selected area. In addition, a new refresh operator (NRO) has been proposed to enhance the algorithm’s solution quality and exploration abilities, as well as a strong crossover mechanism to balance exploitation and increase population diversity. The proposed method, MAGBO, is tested for its performance on CEC2021 benchmark functions and a complex Photovoltaic system. The results are compared with various optimization methods such as Gradient-based Algorithm (GBO), Gradient-based optimization with ranking mechanisms (EGBO), Slime mould algorithm (SMA), Self-adaptive differential evolution algorithm (SADE), Equilibrium optimizer (EO), new self-organizing hierarchical PSO with jumping time-varying acceleration coefficients (HPSO_TVAC), comprehensive learning particle swarm optimizer (CL-PSO), Improved Jaya Algorithm (IJAYA), performance-guided JAYA algorithm (PGJAYA), and multiple learning backtracking search MLBSAs, which shows that MAGBO outperforms the other methods in this context. The Wilcoxon rank-sum and Friedman statistical tests are used to confirm the validity of MAGBO.
The heart of this research revolves around the advanced development and meticulous validation of the MAGBO technique, tailored for photovoltaic system parameter identification. Specifically:
  • Modified Gradient Search Rule (MGSR): utilizes the quasi-Newton method to bolster both local and global optimization.
  • Crossover Mechanism: introduced to ensure greater agent diversity and prevent premature convergence.
  • Novel Refresh Operator (NRO): enhances solution quality and strategic exploration, balancing exploration and exploitation.
  • Rigorous Validation: comprehensive evaluation of various SDMs, DDMs, and PV modules, showcasing MAGBO’s superiority over older methods.
To organize the remainder of this paper, the second section focuses on MAGBO’s related work. The original GBO and problem statement are presented in Section 3. Section 4 explains the details of the proposed MAGBO. In Section 5, experimental results and comparisons are presented. Section 6 discusses the problematic constraints and difficulties, and finally, Section 7 concludes with the conclusions and future work.

2. Related Work

In addition to analytical and deterministic methods, heuristic methods have also been proposed for determining the parameters of PV models. Heuristic methods are typically based on the principles of human reasoning and problem-solving, and they can be used to find approximate solutions to complex problems. Some of the most popular heuristic methods for determining the parameters of PV models include genetic algorithms (GA), artificial bee colonies (ABC), particle swarm optimization (PSO), ant lion optimizer (ALO), biogeography-based optimization (BBO), whale optimization algorithm (WOA), JAYA algorithm, differential evolution (DE), teaching-learning-based optimization (TLBO), harmony search (HS), flower pollination algorithm (FPA), and moth flame optimization algorithm (MFO). Each of these methods has its advantages and disadvantages, and the best method for a particular problem will depend on the specific characteristics of the problem [22,49,57,58,59]. Even though it has been shown that these methods are effective and consistent, most heuristic algorithms still struggle with locating the best possible solution for the whole problem. Estimating parameters in PV models is difficult since the objective function is multimodal, meaning it has numerous local minima. This makes parameter estimation a complex task. Because of this, locating the global minimum, which represents the optimal course of action, can be challenging. In addition, numerous heuristic methods call for particular control parameters, in addition to population size and the number of generations, which substantially influence the overall performance. Incorrect tuning can result in increased computation expenses as well as software that is insensitive to user input.
As a consequence of this, the issue of establishing a trustworthy heuristic approach for predicting PV model parameters continues to be challenging. The most recent and applicable research approaches are outlined in Table 1, which provides an overview of them. This study aims to give both a potential solution to the issue by way of an updated version of the Mas algorithm, which is a method for estimating the parameters of solar PV models, and to suggest a new version of the algorithm.

3. Preliminaries

Within this section, the algorithm for gradient-based optimization and the problem statement are broken down into their component parts.

3.1. Principles of Gradient-Based Optimizer (GBO)

The GBO algorithm [17] is one of the effective methods for solving complex optimization problems due to its ability to explore, exploit, converge, and move away from the local optimal solution. The method component has two main factors, GSR search base and LEO local escape factor.

3.1.1. Initialization

The initial population in GBO is a matrix with N rows and d columns, where N is the size of the population and d is the number of variables of the problem. The initial population is randomly generated using the following equation:
X n = L b + r a n d × U b L b
where X n represents the vector member, L b represents the lower-bounded, U b represents the upper-bounded of the problem, and r a n d denotes a random number in the range (0, 1).

3.1.2. Gradient Search Rule (GSR)

The first factor of GBO controls the vector movement to improve the search in the field and achieve better locations in the solution area. The GSR based on Newton’s method was proposed according to [67]. Therefore, the GSR is defined as:
G S R = r a n d × ρ 1 × 2 x × x n y a n m y b n m + ϵ
y a n m and y b n m are employed to search inside a population to enhance diversity and robustness. Equations (3) and (4) display the y a n m and y b n m formula.
y a n m = r a n d × (   u n + 1 x n 2 + r a n d × x )
y b n m = r a n d × (   u n + 1 x n 2 r a n d × x )
u n + 1 = x n m   r a n d × ρ 1 × 2 x × x n m   x w o r s t x b e s t + ϵ + D M )
where r a n d is a random number in the interval of [0, 1], u n + 1 signifies the new solution produced by the GBO, x w o r s t and x b e s t represent the worst and best solutions found throughout the optimization process, respectively, and ϵ denotes a small number in the [0, 0.1] range.
ρ 1 = 2 × r a n d × α α
α = β × sin 3 π 2 + sin β × 3 π 2
β = β min +   β max β min   × (   1 (   4 3   ) 3 ) 2  
where parameter ρ 1 changes based on the sine function to balance exploration and exploitation, β min and β max are 0.2 and 1.2, respectively, and x the difference between the vectors, which are define as:
x = r a n d × τ
τ = x b e s t x d 1 m + δ 2
δ = 2 × r a n d × x d 1 m + x d 2 m + x d 3 m + x d 4 m 4 x n m
where r 1 , r 2 , r 3 , and r 4 ( r 1 r 2 r 3 r 4 n ) are integers randomly chosen from { 1 , 2 , , N } , and τ denotes the step size.
Adding the ( D M ) limit will improve the performance of the algorithm by making better exploitation in the vicinity of the point x n , since ( D M ) will make the movement of the vector x n be in the direction of ( x b e s t x n ) . Therefore, this process creates a suitable local search tendency to promote the convergence speed of the GBO algorithm. The proposed D M is formulated as follows:
D M = r a n d × ρ 2 × x b e s t x n
ρ 2 = 2 × r a n d × α α

3.1.3. Generate New Solutions

The G S R and D M techniques update the population’s position. The three new species X 1 n m   , X 2 n m   , and X 3 n m are used to determine the final population location update in order to better balance the global search and local search. They are defined by the formulas in the following equations:
X 1 n m   = x n m   G S R + D M
X 1 n m   = x n m   r a n d × ρ 1 × 2 x × x n m y a n m y b n m + ϵ + r a n d × ρ 2 × x b e s t x n m
where m is the current of iteration. The problem’s parameters support the exploration process given an appropriate step size. By replacing the current position with the best position, a new value will be generated that can be generated using the following equation:
X 2 n m   = x b e s t r a n d × ρ 1 × 2 x × x n m y a n m y b n m + ϵ + r a n d × ρ 2 × x r 1 m x r 2 m
where r 1 , r 2 are two random numbers. Accordingly, based on the positions X 1 n m , X 2 n m , and the current position x n m , the new solution at the next iteration can be defined as:
X n m + 1 = r i × r j × X 1 n m + 1 r j × X 2 n m + 1 r i × X 3 n m
  X 3 n m =   x n m     ρ 1 × ( X 2 n m   X 1 n m   )
where r i and r j are two random numbers in [0, 1].

3.1.4. Local Escaping Operator (LEO)

The GBO algorithm suggests a LEO method to boost the effectiveness of handling complicated issues. With the help of this operator, the position of X n m + 1 is improved, the population’s diversity is increased, and the algorithm is assisted in leaving the local optimum. LEO starts by producing a random number, r a n d as a random number in the interval [ 0 , 1 ] . Equation (19) displays the LEO updating formula when r a n d is bigger than 0.5. Equation (20) displays the LEO updating formula when r a n d is smaller than 0.5.
X L E O m = X n m + 1 + q 1 × σ 1 × x b e s t σ 2 × x s m + q 2 × ρ 1 × ( σ 3 × X 2 n m X 1 n m + σ 2 × x r 1 m x r 2 m ) / 2
X L E O m = x b e s t + q 1 × σ 1 × x b e s t σ 2 × x s m + q 2 × ρ 1 × ( σ 3 × X 2 n m X 1 n m + σ 2 × x r 1 m x r 2 m ) / 2
where q 1 and q 2 are uniform random numbers in the interval [−1, 1] and σ 2 is a random number from a normal distribution with mean of 0 and standard deviation of 1.
σ 1 = 2 × r a n d   i f   φ 1 < 0.5 1   o t h e r w i s e
σ 2 = r a n d   i f   φ 1 < 0.5 1   o t h e r w i s e
σ 3 = r a n d   i f   φ 1 < 0.5 1   o t h e r w i s e
where r a n d is a random number in the interval of [0, 1] and φ 1 is a number in the interval [0, 1].
To determine the solution x s m in Equations (19) and (20), the following scheme is defined.
x s m = x l   i f   φ 2 < 0.5 x l m   o t h e r w i s e
x l = L b + r a n d × U b L b
where x l is a new solution and is a randomly selected solution of the population ( l { 1 , 2 , , N } ) and φ 2 and r a n d are random numbers in the interval of [0, 1].

3.2. Problem Statement

This section will describe two PV models, namely SDM and DDM, as SDM depends on the PV module (SMM) and the objective function.

3.2.1. Single-Diode Model

The SDM model is widely used because of its simple structure and high accuracy in extracting and explaining the features of solar cells [68,69]. The SDM’s equivalent circuit for the load current is depicted in Figure 1. The combination of the current source, the diode, and the shunt resistor determines the current through the load. This circuit shows both the leakage current and the electrical resistance of the battery [70]. Here is how you can write down the equation for SDM results:
I = I p h I d I s h
where I stands for output current, I p h is the photo-induced current, I d is the diode current, I s h is the shunt resistor current, and as follows:
I d = I 0 exp V + I R s α V t 1
I s h = V + I R S R s h  
where I 0 represents the diode reverse saturation current, α represents the ideal factor, R s represents the series resistance, R s h represents the shunt resistance, V represents the output voltage of the cell, and V t is the thermal conductivity potential and is given as:
V t = k × T q  
where T represents the crossroads temperature in Kelvin, k refers to the Boltzmann constant (1.3806503 ×   10 23 J/K), and q is the electron charge (1.60217646   × 10 19 C).
Substituting Equations (27) and (28) into Equation (17) gives:
I = I p h I 0 exp V + I R s α V t 1 V + I R S R s h

3.2.2. Double-Diode Model

The SDM single-diode model does not consider current recombination loss in the shading layer, which is a limitation of the model. To address this problem, researchers proposed the DDM dual diode model [71]. The two primary diodes work in parallel with the photovoltaic and shunt currents. These diodes can also be described as the first diode acting as a rectifier and the second as a designer to mimic the recombination of current charges and other factors that cannot be controlled. Figure 2 shows the DDM equivalent circuit, and the output equation can be formulated as follows [72]:
I = I p h I d 1 I d 2 I s h
I d 1 = I 01 exp V + I R s α 1 V t 1
I d 2 = I 02 exp V + I R s α 2 V t 1
where I d 1 indicates the first diode current, I d 2 is the second diode current, I 01 indicates the diffusion current, I 02 denotes the saturation current, and α 1 ,   α 2 are the first and second diode ideality factor, sequentially. Substituting Equations (32) and (33) into Equation (34), we obtain:
I = I p h I 01 exp V + I R s α 1 V t 1 I 02 exp V + I R s α 2 V t 1 V + I R S R s h

3.3. Photovoltaic Module

Figure 3 shows that we can obtain a relationship between voltage and current by using several diodes connected in series or parallel. The photoelectric model can be formulated based on SMM with the following formula [72]:
I = I p h N p I o N p exp V + I R s N s α N s V t 1 V + I R s N s / N p R s h N s / N p
where N p denoted to the number of solar cells connected in parallel and N s refers to the number of solar cells connected in series. In SMM we used the series connection so that N p equals one, and this can be expressed by the equation:
I = I p h I o exp V + I R s N s α N s V t 1 V + I R s N s R s h N s
Parameters ( I p h ,   I o , R s , R s h ,   α ) are unknown and are to be found by differentiation.

3.4. The Cost Functions

Parameter problems are often transformed into optimization problems by using optimization algorithms [73,74]. Researchers [75] used the Root Mean Square Error (RMSE) as the target function, as in the equation below:
R M S E = 1 M m = 1 M f ( V m , I m , x ) 2
The objective function for various photovoltaic models is shown in Table 2. M is the number of estimated I-V datasets and x is considered as a vector for measuring the extracted parameters.

4. Multi-Approach Gradient-Based Optimizer (MAGBO)

In this study, the GBO algorithm [17] will be improved by adding a new refresh operator and deriving GSR depending on the quasi-Newton method used to solve the optimization problem, while the NRO operator replaces the LEO factor.

4.1. Improve Gradient Search Rule ( M G S R )

To minimize the function f then the necessary condition to make x n is a minimum point of f ( x ) , where f ´ x = 0 ; in other words, x n is the solution of the equation f ´ x = 0 and Newton’s method is one of the direct methods to find the roots of equation f ´ x = 0 , which is at the same time the minimum value of f ( x ) , now, if we use the Taylor series to express the quadratic approximation of f x when x = x n , we obtain
f x = f x n + f x n ( x x n ) + f ( x n ) ( x x n ) 2 2 !
Let us make the derivative of Equation (38) equal to zero to obtain the minimum value of f x
f x = f x n + f x n x x n = 0
Now, suppose that x n is an approximation of the minimum value of f x , then from Equation (39) we can obtain the special formula for Newton’s iterative method,
x n + 1 = x n f ( x n ) f ( x n )
with the necessary condition in more real-world problems, then the function being minimized f ( x ) is not available in closed form or difficult to differentiable, then the derivatives f ( x n ) and f ( x n ) in Equation (40) can be approximated by a finite difference formula as
f x n = f x n + x + f x n x 2 f ( x n ) 2 ( x 2 )
f x n = f x n + x f x n x 2 x
where x is the small step size, so the Equations (41) and (42) are substituted in the equation to obtain
x n + 1 = x n f x n + x f x n x x 2 f x n + x + f x n x 2 f x n
The iterative process indicated by Equation (43) is known as the quasi-Newton method, where x n x represents the previous position and x n + x represents the next position.
Then, from Equation (18) we can formulated the proposed M G S R as follows:
M G S R = [ f x n + x f x n x ] x 2 [ f x n + x + f x n x 2 x n ]
In view of the ease of dealing with vectors, some modifications are necessary to deal with population-based research, so x n x with x b e s t , x n + x is replaced with x w o r s t and ( x n ), instead of its fitness f x n are inferred, and the M G S R is reformulated as follows:
M G S R = y a n m y b n m x 2 y a n m + y b n m 2 x n
where x is shown in Equation (9), y a n m and y b n m are computed using Equations (3) and (4), respectively, and x n is the current solution.

4.2. The Proposed New Refresh Operator (NRO)

Utilizing the idea of Levy flight, the new refresh operator is an approach that was developed to tackle difficult optimization problems as a local search technique. One heuristic approach to solving tough optimization issues is local search. It constantly swaps out the current answer for one of its neighbors to improve upon it. It is known that the search space neighbors of solution x are limitless. As a result, figuring out how to efficiently find a good area is crucial in local searches. It will attempt to generate a new solution at predetermined intervals by merging two random solutions with the best position found up to this point. Developing an efficient strategy for determining the optimal location to conduct the search is essential to this approach. Algorithm 1 displays the new refresh operator (NRO). During the optimization process, it will be utilized within each iteration of the proposed algorithm to enhance the search process.
Algorithm 1: The developed new refresh operator algorithm
Sort the particles in p o p m according to fitness functions
    For i = 1: Npop
      Randomly choose two solutions x d 1 m and x d 1 m from p o p m
      Calculate t h e   r a n d o m   v a l u e   z   b a s e d   o n   L e v y f l i g h t for solution n by Equation (48)
      Obtain for solution n by Equation (47)
      Calculate the new position x n m + 1 of solution n by Equation (46)
    End For
At first, in each iteration m for a population p o p m with size of N , search agents and a solution x n m = ( x n , 1 m , x n , 2 m , x n , d i m m ) will be generated by GSR as X 1 n m and X 2 n m according to the proposed algorithm and become x n m + 1 . After that, this search agent will be enhanced by NRO to become X n m + 1 by using this formula:
x n m + 1 = x n m + 1 + x b e s t x p n m + r a n d × x d 1 m x d 2 m   r 1 < 0.5 x b e s t + x b e s t x p n m + r a n d × x d 1 m x d 2 m   r 1 > 0.5
where n = 1 , 2 , , N , r a n d represent a random number in the range 0 and 1, m denotes to the current iteration, and
x p n t = x n m + 1 + z × ( r a n d × x d 1 m x d 2 m ) / 2
where x d 1 m and x d 1 m are two solutions chosen randomly from the population p o p m (expect the current particle x n m ), r a n d is a random number in the range 0 and 1, and t is the number of iterations. In addition, z is random numbers that are generated by the magenta method based on Levy distribution by the following equation:
z = 0.01 × b × ω q 1 λ
where b and q are drawn from normal distributions as b ~ N 0 , λ 2   and   q ~ N 0 , λ 2 , λ is constant with value 1.5 , and
ω = Γ 1 + λ × s i n ( π λ 2 ) Γ λ + 1 2 × λ × 2 λ 1 2 λ

4.3. Generate New Solutions

Depending on the proposed M G S R in Equation (45) and the direction of movement D M , we can generate a new position by updating the current position, and its equation can be formulated as follows:
X 1 n m = X n m M G S R + D M
X 1 n m   = x n m     y a n m y b n m x 2 y a n m + y b n m 2 x n + r a n d × ρ 2 ×   x b e s t x n m  
X 2 n m =   x b e s t M G S R + D M
X 2 n m   =   x b e s t y a n m y b n m x 2 y a n m + y b n m 2 x n + r a n d × ρ 2 ×   x r 1 x n m  
where r 1 are particles chosen randomly from the population p o p m . Accordingly, based on the positions X 1 n m   , X 2 n m   , and the current position x n m   , the new solution at the next iteration will be generated using Equation (17).
By the proposed local search, we can explore its local region and obtain the optimal neighborhood of this particle as described in Algorithm 1.
In the final stage of the proposed local search NRO, the particle x n m   will be recombined (crossover) with its neighbor x n m + 1 , which is obtained by Equation (46) to obtain the best neighbor for the particle x n m + 1 .
x n m + 1 = x n m + 1   µ < 0.1 x n m     µ > 0.1
The proposed algorithm pseudo-code is shown in Algorithm 2, and Figure 4 shows the MAGBO essential steps. Figure 5 illustrates the main idea of problem formulation.
Algorithm 2: The MAGBO
Define the set of parameters, N, Tmax, dim, lb and ub bounds of feasible region.
Initialize a population of N random solutions.
Calculate Fitness, BestX, Best_Cost, WorstX, Worst_Cost
While (Termination condition not satisfied) Do
   Compute  α by Equation (7) and β by Equation (8)
  For i = 1: N
    Compute  X 1 n m by Equation (51)
    Compute  X 2 n m by Equation (53)
      For j = 1: dim
        Calculate  X n m + 1 by Equation (17)
       If  r a n d < 0.5
        Apply local search technique as Algorithm (1)
        End if
    Apply Crossover operator as equation Equation (54)
     Update Fitness, BestX, Best_Cost, WorstX, Worst_Cost
   End for
   End while

5. Computer Results and Simulations

In this section, the primary objective centers around an in-depth examination of the newly proposed MAGBO algorithm. To ensure a comprehensive analysis, we juxtapose MAGBO against several state-of-the-art algorithms, utilizing the challenging functions from the IEEE CEC 2021 test suite. This rigorous benchmarking offers a solid foundation for assessing MAGBO’s optimization capabilities. To supplement the numerical findings, the Wilcoxon signed-rank test, a robust statistical method, is incorporated, lending further weight to the comparative performance claims. Additionally, a qualitative review delves into the intricacies and unique advantages inherent in the MAGBO algorithm. Further fortifying our analysis, MAGBO is applied to the practical task of solar photovoltaic model parameter identification, the results of which are scrutinized through statistical outputs and convergence curve presentations. In the subsequent Section 6, we transition to an introspective reflection on potential challenges and limitations faced during our research, ensuring that the reader is provided a balanced and holistic view of our work, while also hinting at avenues for future improvements.
All of the evaluations for this investigation were carried out on a personal computer operating under the Windows operating system and including a 2.50 GHz Intel Core i5-7300U processor along with 8 gigabytes of random-access memory (RAM). Python was utilized in the implementation of both MAGBO and its rivals.
MAGBO’s efficiency was evaluated concerning shifted, rotated, and biased functions using the IEEE CEC 2021 test suite in 20 dimensions. This was done so that MAGBO’s performance could be thoroughly investigated by comparing it experimentally to that of rival algorithms. Additional details on the IEEE CEC 2021 test suite may be found in [76].
The experimental data used for the statistical analysis were obtained by carrying out 30 sets of replicated runs, each consisting of 2500 iterations. The MAGBO algorithm will terminate once the maximum number of iterations has been completed. We evaluate the algorithm’s performance by comparing the average solution (Avg), the median solution (Med), and the standard deviation. All of this is done without compromising generality (std).
In order to investigate the results of both the MAGBO and the competitors, a pair of non-parametric statistical hypothesis tests called the Friedman test and the Wilcoxon signed-rank test are utilized. The Wilcoxon signed-rank test is used to compare the performance of MAGBO and its competitors across several different metrics, such as how effectively they handle specific tasks, to draw conclusions about which of the two is the superior option. In addition, the Friedman test’s final rankings for algorithms on all functions can be used to examine the substantial differences in overall performance between different algorithms.
The effectiveness of the MAGBO is analyzed and compared to that of eight other algorithms throughout the course of this research project. These techniques include GBO [17], EGBO [77], SMA [78], PSO [79], EO [80], SADE [81], CL-PSO [82], and HPSO_TVAC [83]. In this section, we have decided to adhere to the suggestions made by the authors of those other publications concerning the important aspects of competition, and we have outlined them in Table 3.

5.1. Comparison of MAGBO with State-of-the-Art Competitors Using the IEEE CEC 2021 Test Suite

MAGBO was thoroughly evaluated using the IEEE CEC 2021 test suite, and its results were compared to those of eight other algorithms. The results, shown in Table 4, indicate that MAGBO performed well in terms of average values for most of the test functions, except for f5, f7, and f10. In terms of standard deviation, MAGBO performed significantly better than the other algorithms in almost all cases. These results suggest that MAGBO is a stable, robust, and scalable optimization algorithm that performs well compared to the other algorithms included in the comparison.
The results of the Wilcoxon signed-rank test, which was conducted with a significance level of 0.05, are summarized in Table 5. These results show that MAGBO outperformed several other algorithms, including HPSO_TVAC, CL-PSO, and PSO, in all CEC 2021 test functions. In addition, MAGBO outperformed EGBO, GBO, SMA, and EO on nine functions and SADE on seven functions. These results suggest that MAGBO is a strong performer among the algorithms included in the comparison.
According to Table 4, the convergence rate of the algorithm is improved by using a GBO with a crossover mechanism, the NRO local search technique, and a modified gradient search rule based on the quasi-Newton. Together, the presented methods enable the MAGBO algorithm to generate novel solutions by fusing the properties of two or more existing solutions. This increases the likelihood of finding superior solutions by introducing variety into the solution pool. Keeping the algorithm from settling on a locally optimal solution also facilitates exploration and exploitation. Crossover also makes it easier to probe a wider region of the search space. This is of utmost significance in optimization problems involving complicated and rough landscapes, in which the global optimum may be hidden in a previously undiscovered region. Additionally, the MAGBO algorithm seeks to find a middle ground between exploitation (enhancing the best-known solutions) and exploration (looking for new and perhaps superior answers). While some operators, such as mutation, aid in exploitation, the suggested NRO aids exploration by producing more candidate solutions. In addition, the MAGBO algorithm’s convergence can be hastened with the help of a modified gradient search rule based on the quasi-Newton, which may ultimately lead to the generation of superior solutions.
The Wilcoxon signed-rank test showed that MAGBO is a highly reliable optimization technique compared to eight other algorithms. Figure 6 illustrates the convergence behaviors of MAGBO and its competitors on the CEC 2021 test functions. It can be observed that MAGBO significantly outperforms the other algorithms in both initial and final iterations, indicating that it has strong exploration (using swarm intelligence) and exploitation (using quadratic functions) capabilities. In addition, as shown in Figure 6, MAGBO can find the global optimal solution and avoid getting stuck in local optima thanks to its new local search feature. This feature allows the algorithm to explore the solution space more thoroughly and find the best possible solution.

5.2. Qualitative Analysis of MAGBO

The researchers in this article combined a multi-strategy—the modified gradient search rule (MGSR) derived from the quasi-Newton method to improve its local and global capabilities, a new refresh operator (NRO) to enhance solution quality and algorithm exploration abilities, and a crossover mechanism to balance exploitation and boost the population’s variety—to create a novel approach (MAGBO) that overcomes the shortcomings of the standard GBO.
MAGBO is a new method for dealing with GBO-related problems such as uneven exploration and exploitation, premature convergence, and population heterogeneity. Functions 1, 2, 3, 8, and 9 are chosen as exemplary examples from the CEC 2021 set for this study. These procedures were picked because they effectively demonstrate both unimodal and multimodal features.
As part of a qualitative study of the system’s performance, Figure 7 gives an in-depth look at how MAGBO acts when looking for solutions in unimodal and multimodal functions. The graph makes it easy to see where MAGBO stands and how fit it is at any given time during the search. This shows the average global fitness level of MAGBO over the exploration and exploitation phases. MAGBO’s location in the first dimension is also shown on the graph, along with its evolution through time. The variations in MAGBO’s global best fitness (mean) throughout the iterations show how its normal fitness level varies. This visual depiction accurately depicts the research and development stages of MAGBO throughout the entire refinement process.
MAGBO’s primary purpose of exploration is best illustrated by the fact that its first-dimensional trajectory can represent multiple regions of the search space. The MAGBO particle can quickly and correctly find the optimal solution thanks to its rapid oscillation during the prophase and moderate oscillation during the anaphase. Figure 7b shows that MAGBO’s location curve has a fairly large amplitude early on, which could cover as much as half of the exploration zone. The position amplitude of a MAGBO particle diminishes with time if the function is smooth. However, the position amplitude also shows substantial changes when the function amplitudes are in flux. This exemplifies the flexibility and dependability of MAGBO in many contexts. Depending on the point of view, the changes in amplitude can seem either spectacular or subtle. MAGBO’s substantial early fluctuations attest to the strength of its search skills, while the subtle alterations later on attest to the persistence with which it seeks to locate the ideal solution.
Figure 7c is a graphical representation of MAGBO’s discovery and development processes. The graph has two curves, one blue for the exploration phase and one orange for the exploitation phase of the method. Due to a smaller percentage of exploration relative to exploitation, the initial MAGBO displays a larger exploration ratio, as illustrated in Figure 7c. However, as the iteration count rises, the algorithm begins to shift its attention from exploration to exploitation for the vast majority of the chosen functions. This continuous pattern guarantees that MAGBO always keeps the right balance between exploration and exploitation.
Figure 7d, which depicts an average global fitness curve, shows that MAGBO’s fitness varies when using the iterative method. The curve has a lot of wiggle room, as can be seen by looking at it closely. The average fitness value drops monotonically with increasing iterations, as does the frequency of oscillation. This indicates that MAGBO searches exhaustively during the anaphase and quickly converges on a solid solution.
The outcomes of diversity analysis on the CEC2021 dataset using the MAGBO algorithm are presented in Figure 7e. The depicted graph illustrates the relationship between the number of iterations and the diversity measure, with the horizontal axis representing the former and the vertical axis representing the latter. The process commences by initializing the population with a random generation, leading to increased levels of diversity during the initial phases. As the iterations go, diversity gradually decreases. It is important to acknowledge that the CEC2021 dataset has a wide range of functions, including unimodal, multimodal, hybrid, and composition functions.
The visual representation presented in Figure 7e clearly demonstrates that MAGBO regularly displays a more pronounced decrease in average diversity across different functions. The aforementioned finding indicates the overall accelerated convergence rate of the MAGBO algorithm. Moreover, the findings highlight the general dominance of MAGBO in terms of its performance, making it a viable option for a diverse range of purposes. The primary cause of this phenomenon can be traced to the integration of a “Crossover mechanism” as outlined in Equation (54) within the MAGBO algorithm.
The inclusion of this mechanism in the MAGBO algorithm grants it the capacity to produce innovative solutions through the amalgamation of attributes derived from two or more parent solutions. As a result, the introduction of variation within the population of solutions serves to prevent early convergence and increases the likelihood of finding superior solutions. Furthermore, the integration of crossover facilitates a thorough investigation of the solution space. This characteristic is of great importance in optimization issues that are characterized by complex and uneven terrains, where the global optimum solution may be located in a distant and unexplored area.
However, it is crucial to emphasize that the effectiveness of crossover depends on the particular problem being addressed and the design of the crossover operator. In specific cases, the use of unsuitable crossover operators or parameter settings can hinder the effectiveness of bio-inspired algorithms. Therefore, it is crucial to carefully design and refine every aspect of the algorithm, such as the crossover mechanism, to correspond with the unique characteristics of the optimization problem being addressed.

5.3. The Suggested Approach for Identification of Solar Photovoltaic Model Parameters

In this study, MAGBO is applied to three different models to show its efficiency: the single, double, and PV models. The benchmark data used in this study from [84] were collected from 36 polycrystalline PV cells and a monocrystalline STM6-40/36 module, both under conditions of 1000 W/m2 irradiation and at different temperatures. These data have been widely used to evaluate different methods for estimating PV model parameters. In previous research, the parameter ranges for PV cells and modules have been kept constant across all studies to ensure a consistent search space. As shown in Table 6, the parameter ranges for PV cells and modules are provided. Eight algorithms were selected for comparison; the population size N is 30, the maximum frequency is 600 (18,000 maximum evaluation), and the number of independent runs is 30.

5.3.1. The Single-Diode Model

The single-diode module’s current-voltage (I-V) and voltage-current (P-V) curves are displayed. Figure 8 illustrates the MAGBO mistake that was present on the diode. As shown in Table 7, MAGBO achieved the best results compared to other algorithms; therefore, it is reasonable to conclude that MAGBO possesses the potential to be an effective instrument for the identification of SDM. In addition, Figure 9 illustrates the absolute and relative differences in current value that are shown by making a comparison between the data that were simulated and the data that were observed.
The values of five parameters and the RME are presented in Table 8, demonstrating that the suggested method generates superior results compared to other algorithms already in use. The findings above suggest that the proposed method could be utilized as an efficient approach to the localization of single-diode models (SDM).

5.3.2. Double-Diode Module Results (DDM)

Figure 10 compares the actual observed data and the estimated model properties that MAGBO provided. Figure 11 presents a graphical representation of a double-diode module’s RE and IAE values. The simulated current-voltage (I-V) and power-voltage (P-V) characteristic curves are an excellent match for the experimental data collected. The simulation results are presented in Table 9, which includes the current, power, IAE, and RE numbers. The results presented in Table 9 demonstrate that MAGBO’s double-diode module successfully reproduces the primary properties of solar cells. Table 10 presents the results of comparisons between the performance of MAGBO and that of alternative algorithms. Table 10 presents the RMSE comparison findings and the values of the seven extracted parameters taken out of the model. Compared to the other seven methods, the Root Mean Square Error (RMSE) for the suggested MAGBO is as low as possible. Table 10 presents the results of comparisons between the performance of MAGBO and that of alternative algorithms. Table 10 presents the RMSE comparison findings and the values of the seven extracted parameters taken out of the model. When contrasted with the other seven approaches, the RMSE for the proposed MAGBO is the lowest it can be.

5.3.3. PV Module (PV)

Figure 12 illustrates the high degree of consistency between the simulated and collected data, with I-V and P-V curves that are realistic approximations of the properties of the Photowatt-PWP 201 module model. Figure 13 displays the experimentally gathered and computer-generated representations of observed currents for the IAE and RE of the PV module model. Table 11 presents the findings of this experiment, including the current, power, IAR, and RE values, which support the idea that the values proposed by MAGBO for the model parameters of the PV modules are reliable. Table 12 lists the best RMSE score and the five extracted parameter values from each of the seven techniques’ total of 30 tests. In conclusion, the proposed MAGBO performs better than other methods in forecasting the parameters of PV module models compared to other methods.

5.3.4. Statistical Results and Convergence Curves

Table 13 provides a summary of the statistical information on competing algorithms. (after carrying out 30 separate runs, each of which consisted of 600 iterations) The statistical data are presented in Table 13 as best, worst, Avg, std, and Rank values. By examining these numbers, we can deduce that MAGBO achieves the best results within 30 independent trials with a Rank (Friedman test) of 1. On a single-diode module, a double-diode module, and a PV module, MAGBO performs at a generally competitive level with the most recent robust algorithm. Figure 14 displays the convergence curves of the MAGBO method and those of its competitors.

5.4. Tests MAGBO on a Wide Selection of Solar Cells and Modules

Data from the PVM 752 GaAs thin-film cell at 25 °C and total irradiation (1000 W/m2) and the STP6-40/36 module are used to prove the validity and reliability of MAGBO, a method for detecting the properties of both the SDM and DDM of photovoltaic modules. The simulation parameters were determined using 44 pairs of I-V points from the PVM 752 cell [85] and data from the STP6-40/36 module, which consisted of 36 series-connected polycrystalline cells and were tested at temperatures of 51 °C and 55 °C, along with measurements at varying temperatures and irradiance levels. Table 6 details how we used actual data from STM6-40/36 and PVM 752 GaAs thin-film modules to determine the simulation experiment’s settings. The parameters of the compared algorithms are listed in Table 3. Each issue runs using the same maximum number of iterations (600) across all methods.

5.4.1. Results Obtained Using the STM6-40/36 Simulation Model

In order to evaluate how well the proposed method compares to preexisting algorithms, we publish the results from the STM6-40/36 model [86]. Figure 15 shows that the I-V and P-V curves generated by the simulation and the observations are highly consistent with one another. This supports the idea that MAGBO’s projections about the STM6-40/36 model’s characteristics are accurate. The STM6-40/36 model’s IAE and RE are also depicted in Figure 16. The complete experimental results, including current, power, and corresponding IAR and RE values, are shown in Table 14. It was found that IAE is less than 6.08803600 × 10 3 , which verifies that MAGBO’s determination of the STM6-40/36 model’s parameters is accurate. Table 15 displays the best RMSE values for each of the seven approaches and the five most significant extracted parameter values based on a total of 30 tests. Table 15 provides the RMSE values that are the most accurate, in addition to the five extracted parameter values that are the most significant, for each of the seven methods evaluated using a total of 30 tests. Table 15 shows that MAGBO produces the best results with an RMSE value of 1.72981457 × 10 3 , the lowest possible value. Compared to other methods, the proposed MAGBO method is superior in terms of its performance when estimating the parameters of the STM6-40/36 model.

5.4.2. Results Obtained Using the PVM 752 GaAs Thin-Film Simulation Model

Using the SDM, we analyzed the suggested MAGBO technique to learn more about the characteristics of a PVM 752 GaAs thin-film cell. Table 16 and Table 17 detail the experiment outcomes, such as the current, power, IAE, and RE measured values. Table 17 displays the simulated current data alongside the IAE and RE error measurements. The optimization method was found to be more challenging when applied to calculating the parameters of the PVM 752 thin-film cell than the RTC France solar cell. Table 16 shows that despite this, the MAGBO approach is effective. Table 17 displays the outcomes of different algorithms’ attempts to estimate the SDM parameters, and it can be observed that MAGBO fared best, with the lowest RMSE value. Figure 17 shows the IAE and RE currents, both experimentally and as measured, for the PVM 752 thin-film cell model. Figure 18 demonstrates that the I-V and P-V curves generated by the simulation and the measurements are highly consistent with one another, proving that MAGBO’s predictions for the PVM 752 thin-film model are correct.

6. Problem Limitation

Even if it performs better than certain other algorithms already on the market for the chosen tests, MAGBO has a few drawbacks. The following are some of our ideas regarding the primary restrictions imposed by the algorithm: to begin, the research does not address a large-scale optimization problem that involves more than 500 choice variables related to scalability. Second, the suggested method can successfully obtain optimal global solutions and is utilized to optimize the process of extracting parameters from static models of solar photovoltaic systems. In addition, the report completely ignores the difficulties associated with the dynamic model. Finally, while selecting multi-objective real-world issues, it is a good idea to consider a mix of computing time, complexity, stability, and scalability. This is because these factors all play a role in the solution.

7. Conclusions and Future Work

Solar energy is one of the most important forms of renewable energy that has the potential to replace fossil fuels. Unlike fossil fuels, solar energy does not release any harmful gaseous emissions that are a major factor in the degradation of the natural environment and a consequent risk to the well-being of all forms of life. The vital activity of parameter extraction is a component of the optimization process for solar power systems. In this research, optimal solutions can be discovered by utilizing a more recent and improved algorithm known as MAGBO. The modified gradient search rule (MGSR), which was derived from the quasi-Newton method to improve its local and global capabilities, a new refresh operator (NRO), which was designed to improve solution quality and algorithm exploration abilities, and a crossover mechanism, which was designed to balance exploitation and boost the population’s variety, were all combined by MAGBO to create a novel approach that overcomes the shortcomings of the standard GBO. A precise evaluation of MAGBO’s usefulness was made possible thanks to the successful extraction of parameters from various PV models, such as single and double-diode (DDM), PV modules, as well as the CEC2021 test suite. The results revealed that the recently developed MAGBO algorithm performed better than its contemporaries, including the first version of the GBO algorithm. According to the findings, MAGBO can compete with numerous other algorithms already in use, establishing it as a viable alternative approach for determining the parameters of solar models. Further study will present additional strategies to complement the suggested algorithm to address other energy-related optimization difficulties, particularly the dynamic models. This will be done to improve the overall quality of the proposed solution. In addition, the technique that has been described can be utilized to address large-scale global optimization issues without compromising diversity. In addition, MAGBO has not been considered from a medical point of view at any step in this process. As a direct consequence, the MAGBO that has been proposed will conduct an experiment on this subject.

Author Contributions

Conceptualization, M.Q. and N.K.H.; methodology, M.Q. and B.S.M.; software, M.Q. and M.A.; validation, M.A., M.Q. and N.K.H.; formal analysis, M.A.; investigation, B.S.M. and M.A.; resources, M.A.; data curation, B.S.M., M.A. and M.Q.; writing—original draft preparation, N.K.H.; writing—review and editing, M.Q., B.S.M., M.A. and N.K.H.; visualization, M.Q., B.S.M., M.A. and N.K.H.; supervision, M.Q., B.S.M., M.A. and N.K.H.; project administration, M.Q. and N.K.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia, project number 445-9-394.

Data Availability Statement

We have made the code for this research project available on GitHub. at https://github.com/MohammedQaraad/MAGBO-algorithm, (accessed on 22 September 2023). Readers interested in accessing the code can download it from this repository.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia, for funding this research work through project number 445-9-394.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mcelroy, M.B.; Chen, X. Wind and Solar Power in the United States: Status and Prospects. CSEE J. Power Energy Syst. 2017, 3, 1–6. [Google Scholar] [CrossRef]
  2. Alam, D.F.; Yousri, D.A.; Eteiba, M.B. Flower Pollination Algorithm Based Solar PV Parameter Estimation. Energy Convers. Manag. 2015, 101, 410–422. [Google Scholar] [CrossRef]
  3. Fathy, A.; Rezk, H. Parameter Estimation of Photovoltaic System Using Imperialist Competitive Algorithm. Renew. Energy 2017, 111, 307–320. [Google Scholar] [CrossRef]
  4. Chin, V.J.; Salam, Z.; Ishaque, K. Cell Modelling and Model Parameters Estimation Techniques for Photovoltaic Simulator Application: A Review. Appl. Energy 2015, 154, 500–519. [Google Scholar] [CrossRef]
  5. Batzelis, E.I.; Papathanassiou, S.A. A Method for the Analytical Extraction of the Single-Diode PV Model Parameters. IEEE Trans. Sustain. Energy 2015, 7, 504–512. [Google Scholar]
  6. Louzazni, M.; Aroudam, E.H. An Analytical Mathematical Modeling to Extract the Parameters of Solar Cell from Implicit Equation to Explicit Form. Appl. Sol. Energy 2015, 51, 165–171. [Google Scholar]
  7. Brano, V.L.; Ciulla, G. An Efficient Analytical Approach for Obtaining a Five Parameters Model of Photovoltaic Modules Using Only Reference Data. Appl. Energy 2013, 111, 894–903. [Google Scholar]
  8. Stornelli, V.; Muttillo, M.; De Rubeis, T.; Nardi, I. A New Simplified Five-Parameter Estimation Method for Single-Diode Model of Photovoltaic Panels. Energies 2019, 12, 4271. [Google Scholar]
  9. Atay, B.K.; Eminoğlu, U. A New Approach for Parameter Estimation of the Single-Diode Model Forphotovoltaic Cells/Modules. Turkish J. Electr. Eng. Comput. Sci. 2019, 27, 3026–3039. [Google Scholar] [CrossRef]
  10. Qaraad, M.; Amjad, S.; Hussein, N.K.; Elhosseini, M.A. Addressing Constrained Engineering Problems and Feature Selection with a Time-Based Leadership Salp-Based Algorithm with Competitive Learning. J. Comput. Des. Eng. 2022, 9, 2235–2270. [Google Scholar]
  11. Mehta, H.K.; Warke, H.; Kukadiya, K.; Panchal, A.K. Accurate Expressions for Single-Diode-Model Solar Cell Parameterization. IEEE J. Photovolt. 2019, 9, 803–810. [Google Scholar]
  12. Lun, S.; Du, C.; Guo, T.; Wang, S.; Sang, J.; Li, J. A New Explicit I–V Model of a Solar Cell Based on Taylor’s Series Expansion. Sol. Energy 2013, 94, 221–232. [Google Scholar]
  13. 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]
  14. Cabestany, J.; Castaner, L. Evaluation of Solar Cell Parameters by Nonlinear Algorithms. J. Phys. D Appl. Phys. 1983, 16, 2547. [Google Scholar] [CrossRef]
  15. Li, S.; Gong, W.; Gu, Q. A Comprehensive Survey on Meta-Heuristic Algorithms for Parameter Extraction of Photovoltaic Models. Renew. Sustain. Energy Rev. 2021, 141, 110828. [Google Scholar] [CrossRef]
  16. Ismail, M.S.; Moghavvemi, M.; Mahlia, T.M.I. Characterization of PV Panel and Global Optimization of Its Model Parameters Using Genetic Algorithm. Energy Convers. Manag. 2013, 73, 10–25. [Google Scholar]
  17. Ahmadianfar, I.; Bozorg-Haddad, O.; Chu, X. Gradient-Based Optimizer: A New Metaheuristic Optimization Algorithm. Inf. Sci. 2020, 540, 131–159. [Google Scholar] [CrossRef]
  18. Hamid, N.F.A.; Rahim, N.A.; Selvaraj, J. Solar Cell Parameters Extraction Using Particle Swarm Optimization Algorithm. In Proceedings of the 2013 IEEE Conference on Clean Energy and Technology (CEAT), Langkawi, Malaysia, 18–20 November 2013; pp. 461–465. [Google Scholar]
  19. El-Naggar, K.M.; AlRashidi, M.R.; AlHajri, M.F.; Al-Othman, A.K. Simulated Annealing Algorithm for Photovoltaic Parameters Identification. Sol. Energy 2012, 86, 266–274. [Google Scholar]
  20. Oliva, D.; Cuevas, E.; Pajares, G. Parameter Identification of Solar Cells Using Artificial Bee Colony Optimization. Energy 2014, 72, 93–102. [Google Scholar] [CrossRef]
  21. Yuan, X.; He, Y.; Liu, L. Parameter Extraction of Solar Cell Models Using Chaotic Asexual Reproduction Optimization. Neural Comput. Appl. 2015, 26, 1227–1239. [Google Scholar]
  22. Yu, K.; Liang, J.J.; Qu, B.Y.; Cheng, Z.; Wang, H. Multiple Learning Backtracking Search Algorithm for Estimating Parameters of Photovoltaic Models. Appl. Energy 2018, 226, 408–422. [Google Scholar] [CrossRef]
  23. Kang, T.; Yao, J.; Jin, M.; Yang, S.; Duong, T. A Novel Improved Cuckoo Search Algorithm for Parameter Estimation of Photovoltaic (PV) Models. Energies 2018, 11, 1060. [Google Scholar]
  24. Niu, Q.; Zhang, L.; Li, K. A Biogeography-Based Optimization Algorithm with Mutation Strategies for Model Parameter Estimation of Solar and Fuel Cells. Energy Convers. Manag. 2014, 86, 1173–1185. [Google Scholar] [CrossRef]
  25. Wang, G.G.; Deb, S.; Cui, Z. Monarch Butterfly Optimization. Neural Comput. Appl. 2015, 31, 1995–2014. [Google Scholar] [CrossRef]
  26. Qais, M.H.; Hasanien, H.M.; Alghuwainem, S. Identification of Electrical Parameters for Three-Diode Photovoltaic Model Using Analytical and Sunflower Optimization Algorithm. Appl. Energy 2019, 250, 109–117. [Google Scholar]
  27. Qais, M.H.; Hasanien, H.M.; Alghuwainem, S.; Nouh, A.S. Coyote Optimization Algorithm for Parameters Extraction of Three-Diode Photovoltaic Models of Photovoltaic Modules. Energy 2019, 187, 116001. [Google Scholar]
  28. Subudhi, B.; Pradhan, R. Bacterial Foraging Optimization Approach to Parameter Extraction of a Photovoltaic Module. IEEE Trans. Sustain. Energy 2017, 9, 381–389. [Google Scholar] [CrossRef]
  29. Chen, H.; Jiao, S.; Wang, M.; Heidari, A.A.; Zhao, X. Parameters Identification of Photovoltaic Cells and Modules Using Diversification-Enriched Harris Hawks Optimization with Chaotic Drifts. J. Clean. Prod. 2020, 244, 118778. [Google Scholar]
  30. Benkercha, R.; Moulahoum, S.; Taghezouit, B. Extraction of the PV Modules Parameters with MPP Estimation Using the Modified Flower Algorithm. Renew. Energy 2019, 143, 1698–1709. [Google Scholar]
  31. Houssein, E.H.; Zaki, G.N.; Abualigah, L.; Younis, E.M.G. Metaheuristics for Parameter Estimation of Solar Photovoltaic Cells: A Comprehensive Review. Stud. Comput. Intell. 2022, 1038, 149–179. [Google Scholar]
  32. Ridha, H.M.; Hizam, H.; Mirjalili, S.; Othman, M.L.; Ya’acob, M.E.; Ahmadipour, M. Parameter Extraction of Single, Double, and Three Diodes Photovoltaic Model Based on Guaranteed Convergence Arithmetic Optimization Algorithm and Modified Third Order Newton Raphson Methods. Renew. Sustain. Energy Rev. 2022, 162, 112436. [Google Scholar] [CrossRef]
  33. Wolpert, D.H.; Macready, W.G. No Free Lunch Theorems for Optimization. IEEE Trans. Evol. Comput. 1997, 1, 67–82. [Google Scholar] [CrossRef]
  34. Aliramezani, M.; Koch, C.R.; Shahbakhti, M. Modeling, Diagnostics, Optimization, and Control of Internal Combustion Engines via Modern Machine Learning Techniques: A Review and Future Directions; Elsevier: Amsterdam, The Netherlands, 2022. [Google Scholar] [CrossRef]
  35. Baburaj, E. Comparative Analysis of Bio-Inspired Optimization Algorithms in Neural Network-Based Data Mining Classification. Int. J. Swarm Intell. Res. (IJSIR) 2022, 13, 1–25. [Google Scholar]
  36. Brookes, D.H.; Aghazadeh, A.; Listgarten, J. On the Sparsity of Fitness Functions and Implications for Learning. Proc. Natl. Acad. Sci. USA 2022, 119, e2109649118. [Google Scholar] [CrossRef]
  37. Deng, W.; Shang, S.; Cai, X.; Zhao, H.; Zhou, Y.; Chen, H.; Deng, W. Quantum Differential Evolution with Cooperative Coevolution Framework and Hybrid Mutation Strategy for Large Scale Optimization; Elsevier: Amsterdam, The Netherlands, 2021. [Google Scholar]
  38. LaTorre, A.; Molina, D.; Osaba, E.; Poyatos, J.; Del Ser, J.; Herrera, F. A Prescription of Methodological Guidelines for Comparing Bio-Inspired Optimization Algorithms; Elsevier: Amsterdam, The Netherlands, 2021. [Google Scholar]
  39. Deng, W.; Zhang, X.; Zhou, Y.; Liu, Y.; Zhou, X.; Chen, H.; Zhao, H. An Enhanced Fast Non-Dominated Solution Sorting Genetic Algorithm for Multi-Objective Problems; Elsevier: Amsterdam, The Netherlands, 2022. [Google Scholar]
  40. Wu, E.Q.; Zhou, M.; Hu, D.; Zhu, L.; Tang, Z.; Qiu, X.Y.; Deng, P.Y.; Zhu, L.M.; Ren, H. Self-Paced Dynamic Infinite Mixture Model for Fatigue Evaluation of Pilots’ Brains. IEEE Trans. Cybern. 2020, 52, 5623–5638. [Google Scholar]
  41. Mirjalili, S. Moth-Flame Optimization Algorithm: A Novel Nature-Inspired Heuristic Paradigm. Knowl.-Based Syst. 2015, 89, 228–249. [Google Scholar] [CrossRef]
  42. Li, M.D.; Zhao, H.; Weng, X.W.; Han, T. A Novel Nature-Inspired Algorithm for Optimization: Virus Colony Search. Adv. Eng. Softw. 2016, 165, 1339–1351. [Google Scholar] [CrossRef]
  43. Hashim, F.A.; Houssein, E.H.; Mabrouk, M.S.; Al-Atabany, W.; Mirjalili, S. Henry Gas Solubility Optimization: A Novel Physics-Based Algorithm. Futur. Gener. Comput. Syst. 2019, 101, 646–667. [Google Scholar] [CrossRef]
  44. Qaraad, M.; Aljadania, A.; Elhosseini, M. Large-Scale Competitive Learning-Based Salp Swarm for Global Optimization and Solving Constrained Mechanical and Engineering Design Problems. Mathematics 2023, 11, 1362. [Google Scholar] [CrossRef]
  45. 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]
  46. 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]
  47. 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]
  48. Nadimi-Shahraki, M.H.; Taghian, S.; Zamani, H.; Mirjalili, S.; Elaziz, M.A. MMKE: Multi-Trial Vector-Based Monkey King Evolution Algorithm and Its Applications for Engineering Optimization Problems. PLoS ONE 2023, 18, e0280006. [Google Scholar] [CrossRef]
  49. Nadimi-Shahraki, M.H.; Taghian, S.; Mirjalili, S.; Abualigah, L. Binary Aquila Optimizer for Selecting Effective Features from Medical Data: A COVID-19 Case Study. Mathematics 2022, 10, 1929. [Google Scholar] [CrossRef]
  50. Nadimi-Shahraki, M.H.; Moeini, E.; Taghian, S.; Mirjalili, S. DMFO-CD: A Discrete Moth-Flame Optimization Algorithm for Community Detection. Algorithms 2021, 14, 314. [Google Scholar] [CrossRef]
  51. Daoud, M.S.; Shehab, M.; Al-Mimi, H.M.; Abualigah, L.; Zitar, R.A.; Shambour, M.K.Y. Gradient-Based Optimizer (GBO): A Review, Theory, Variants, and Applications. Arch. Comput. Methods Eng. 2022, 1, 1–19. [Google Scholar]
  52. Alorf, A. A Survey of Recently Developed Metaheuristics and Their Comparative Analysis. Eng. Appl. Artif. Intell. 2023, 117, 105622. [Google Scholar] [CrossRef]
  53. 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]
  54. Zhang, J.; Sanderson, A.C. JADE: Adaptive Differential Evolution with Optional External Archive. IEEE Trans. Evol. Comput. 2009, 13, 945–958. [Google Scholar] [CrossRef]
  55. Liang, J.; Qiao, K.; Yu, K.; Ge, S.; Qu, B.; Xu, R.; Li, K. Parameters Estimation of Solar Photovoltaic Models via a Self-Adaptive Ensemble-Based Differential Evolution. Sol. Energy 2020, 207, 336–346. [Google Scholar] [CrossRef]
  56. Hassan, M.H.; Kamel, S.; El-Dabah, M.A.; Rezk, H. A Novel Solution Methodology Based on a Modified Gradient-Based Optimizer for Parameter Estimation of Photovoltaic Models. Electronics 2021, 10, 472. [Google Scholar] [CrossRef]
  57. Jiang, Y.; Luo, Q.; Zhou, Y. Improved Gradient-Based Optimizer for Parameters Extraction of Photovoltaic Models. IET Renew. Power Gener. 2022, 16, 1602–1622. [Google Scholar] [CrossRef]
  58. 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. 2023, 14, 7109–7131. [Google Scholar]
  59. Zhou, W.; Wang, P.; Heidari, A.A.; Zhao, X.; Turabieh, H.; Chen, H. Random Learning Gradient Based Optimization for Efficient Design of Photovoltaic Models. Energy Convers. Manag. 2021, 230, 113751. [Google Scholar] [CrossRef]
  60. Yu, S.; Chen, Z.; Heidari, A.A.; Zhou, W.; Chen, H.; Xiao, L. Parameter Identification of Photovoltaic Models Using a Sine Cosine Differential Gradient Based Optimizer. IET Renew. Power Gener. 2022, 16, 1535–1561. [Google Scholar] [CrossRef]
  61. Ramadan, A.; Kamel, S.; Hassan, M.H.; Tostado-Véliz, M.; Eltamaly, A.M. Parameter Estimation of Static/Dynamic Photovoltaic Models Using a Developed Version of Eagle Strategy Gradient-Based Optimizer. Sustainability 2021, 13, 13053. [Google Scholar] [CrossRef]
  62. Qaraad, M.; Amjad, S.; Hussein, N.K.; Badawy, M.; Mirjalili, S.; Elhosseini, M.A. Photovoltaic Parameter Estimation Using Improved Moth Flame Algorithms with Local Escape Operators. Comput. Electr. Eng. 2023, 106, 108603. [Google Scholar] [CrossRef]
  63. Özban, A.Y. Some New Variants of Newton’s Method. Appl. Math. Lett. 2004, 17, 677–682. [Google Scholar] [CrossRef]
  64. Liu, Y.; Chong, G.; Heidari, A.A.; Chen, H.; Liang, G.; Ye, X.; Cai, Z.; Wang, M. Horizontal and Vertical Crossover of Harris Hawk Optimizer with Nelder-Mead Simplex for Parameter Estimation of Photovoltaic Models. Energy Convers. Manag. 2020, 223, 113211. [Google Scholar]
  65. Fan, Y.; Wang, P.; Heidari, A.A.; Zhao, X.; Turabieh, H.; Chen, H. Delayed Dynamic Step Shuffling Frog-Leaping Algorithm for Optimal Design of Photovoltaic Models. Energy Rep. 2021, 7, 228–246. [Google Scholar]
  66. Yu, K.; Chen, X.; Wang, X.; Wang, Z. Parameters Identification of Photovoltaic Models Using Self-Adaptive Teaching-Learning-Based Optimization. Energy Convers. Manag. 2017, 145, 233–246. [Google Scholar]
  67. 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]
  68. Li, S.; Gong, W.; Yan, X.; Hu, C.; Bai, D.; Wang, L.; Gao, L. Parameter Extraction of Photovoltaic Models Using an Improved Teaching-Learning-Based Optimization. Energy Convers. Manag. 2019, 186, 293–305. [Google Scholar] [CrossRef]
  69. 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]
  70. Ridha, H.M.; Gomes, C.; Hizam, H.; Ahmadipour, M.; Heidari, A.A.; Chen, H. Multi-Objective Optimization and Multi-Criteria Decision-Making Methods for Optimal Design of Standalone Photovoltaic System: A Comprehensive Review. Renew. Sustain. Energy Rev. 2021, 135, 110202. [Google Scholar]
  71. Zhang, H.; Heidari, A.A.; Wang, M.; Zhang, L.; Chen, H.; Li, C. Orthogonal Nelder-Mead Moth Flame Method for Parameters Identification of Photovoltaic Modules. Energy Convers. Manag. 2020, 211, 112764. [Google Scholar]
  72. Abbassi, R.; Abbassi, A.; Heidari, A.A.; Mirjalili, S. An Efficient Salp Swarm-Inspired Algorithm for Parameters Identification of Photovoltaic Cell Models. Energy Convers. Manag. 2019, 179, 362–372. [Google Scholar] [CrossRef]
  73. Special Session on Real. Available online: https://www3.ntu.edu.sg/home/epnsugan/index_files/CEC2021/CEC2021-2.htm (accessed on 23 December 2022).
  74. Ahmadianfar, I.; Gong, W.; Heidari, A.A.; Golilarz, N.A.; Samadi-Koucheksaraee, A.; Chen, H. Gradient-Based Optimization with Ranking Mechanisms for Parameter Identification of Photovoltaic Systems. Energy Rep. 2021, 7, 3979–3997. [Google Scholar] [CrossRef]
  75. Li, S.; Chen, H.; Wang, M.; Heidari, A.A.; Mirjalili, S. Slime Mould Algorithm: A New Method for Stochastic Optimization. Future Gener. Comput. Syst. 2020, 111, 300–323. [Google Scholar] [CrossRef]
  76. 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]
  77. Faramarzi, A.; Heidarinejad, M.; Stephens, B.; Mirjalili, S. Equilibrium Optimizer: A Novel Optimization Algorithm. Knowl.-Based Syst. 2020, 191, 105190. [Google Scholar] [CrossRef]
  78. Qin, A.K.; Huang, V.L.; Suganthan, P.N. Differential Evolution Algorithm with Strategy Adaptation for Global Numerical Optimization. IEEE Trans. Evol. Comput. 2009, 13, 398–417. [Google Scholar] [CrossRef]
  79. 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]
  80. Ghasemi, M.; Aghaei, J.; Hadipour, M. New Self-Organising Hierarchical PSO with Jumping Time-Varying Acceleration Coefficients. Electron. Lett. 2017, 53, 1360–1362. [Google Scholar] [CrossRef]
  81. Bakır, H. Comparative Performance Analysis of Metaheuristic Search Algorithms in Parameter Extraction for Various Solar Cell Models. Environ. Chall. 2023, 11, 100720. [Google Scholar] [CrossRef]
  82. Rezaee Jordehi, A. Enhanced Leader Particle Swarm Optimisation (ELPSO): An Efficient Algorithm for Parameter Estimation of Photovoltaic (PV) Cells and Modules. Sol. Energy 2018, 159, 78–87. [Google Scholar] [CrossRef]
  83. Tong, N.T.; Pora, W. A Parameter Extraction Technique Exploiting Intrinsic Properties of Solar Cells. Appl. Energy 2016, 176, 104–115. [Google Scholar] [CrossRef]
  84. Yu, S.; Heidari, A.A.; He, C.; Cai, Z.; Althobaiti, M.M.; Mansour, R.F.; Liang, G.; Chen, H. Parameter Estimation of Static Solar Photovoltaic Models Using Laplacian Nelder-Mead Hunger Games Search. Sol. Energy 2022, 242, 79–104. [Google Scholar] [CrossRef]
  85. Long, W.; Jiao, J.; Liang, X.; Xu, M.; Tang, M.; Cai, S. Parameters Estimation of Photovoltaic Models Using a Novel Hybrid Seagull Optimization Algorithm. Energy 2022, 249, 123760. [Google Scholar] [CrossRef]
  86. Liang, J.; Qiao, K.; Yu, K.; Qu, B.; Yue, C.; Guo, W.; Wang, L. Utilizing the Relationship Between Unconstrained and Constrained Pareto Fronts for Constrained Multiobjective Optimization. IEEE Trans. Cybern. 2023, 53, 3873–3886. [Google Scholar] [CrossRef]
  87. Yu, K.; Zhang, D.; Liang, J.; Chen, K.; Yue, C.; Qiao, K.; Wang, L. A Correlation-Guided Layered Prediction Approach for Evolutionary Dynamic Multiobjective Optimization. IEEE Trans. Evol. Comput. 2022. [Google Scholar] [CrossRef]
Figure 1. Equivalent circuit diagram of the SDM model.
Figure 1. Equivalent circuit diagram of the SDM model.
Mathematics 11 04200 g001
Figure 2. Equivalent circuit diagram of the DDM model.
Figure 2. Equivalent circuit diagram of the DDM model.
Mathematics 11 04200 g002
Figure 3. Equivalent circuit diagram of the SMM model.
Figure 3. Equivalent circuit diagram of the SMM model.
Mathematics 11 04200 g003
Figure 4. Problem formulation.
Figure 4. Problem formulation.
Mathematics 11 04200 g004
Figure 5. MAGBO Flow chart.
Figure 5. MAGBO Flow chart.
Mathematics 11 04200 g005
Figure 6. MAGBO convergence curves vs. some basics and advanced algorithms during 2500 iterations.
Figure 6. MAGBO convergence curves vs. some basics and advanced algorithms during 2500 iterations.
Mathematics 11 04200 g006aMathematics 11 04200 g006b
Figure 7. (a) MAGBO qualitative analysis (Cec2021_functions). (b) MAGBO qualitative analysis (first-dimensional trajectory). (c) MAGBO qualitative analysis (phase of exploration and exploitation). (d) MAGBO qualitative analysis (global best fitness—average). (e) MAGBO qualitative analysis (diversity analysis—average).
Figure 7. (a) MAGBO qualitative analysis (Cec2021_functions). (b) MAGBO qualitative analysis (first-dimensional trajectory). (c) MAGBO qualitative analysis (phase of exploration and exploitation). (d) MAGBO qualitative analysis (global best fitness—average). (e) MAGBO qualitative analysis (diversity analysis—average).
Mathematics 11 04200 g007aMathematics 11 04200 g007bMathematics 11 04200 g007c
Figure 8. Experimental and simulated current-voltage I–V (a) and power-voltage (P–V) characteristics (b) for a single diode by MAGBO.
Figure 8. Experimental and simulated current-voltage I–V (a) and power-voltage (P–V) characteristics (b) for a single diode by MAGBO.
Mathematics 11 04200 g008
Figure 9. Experimental and simulated current data error-index values of SDM: (a) IAE; (b) RE.
Figure 9. Experimental and simulated current data error-index values of SDM: (a) IAE; (b) RE.
Mathematics 11 04200 g009
Figure 10. Experimental and simulated current-voltage I–V (a) and power-voltage (P–V) characteristics (b) for a double diode by MAGBO.
Figure 10. Experimental and simulated current-voltage I–V (a) and power-voltage (P–V) characteristics (b) for a double diode by MAGBO.
Mathematics 11 04200 g010
Figure 11. Experimental and simulated current data error-index values of SDM: (a) IAE; (b) RE.
Figure 11. Experimental and simulated current data error-index values of SDM: (a) IAE; (b) RE.
Mathematics 11 04200 g011
Figure 12. Experimental and simulated current-voltage I-V (a) and power-voltage (P-V) characteristics (b) for a PV module model by MAGBO.
Figure 12. Experimental and simulated current-voltage I-V (a) and power-voltage (P-V) characteristics (b) for a PV module model by MAGBO.
Mathematics 11 04200 g012
Figure 13. Experimental and simulated current data error-index values of a PV module model: (a) IAE; (b) RE.
Figure 13. Experimental and simulated current data error-index values of a PV module model: (a) IAE; (b) RE.
Mathematics 11 04200 g013
Figure 14. MAGBO and other competitors’ convergence curves for three models: (a) Single-diode SD; (b) Double-diode DD; (c) PV module.
Figure 14. MAGBO and other competitors’ convergence curves for three models: (a) Single-diode SD; (b) Double-diode DD; (c) PV module.
Mathematics 11 04200 g014
Figure 15. Experimental and simulated current-voltage I–V (a) and power-voltage (P–V) characteristics (b) for an STM6-40/36 model by MAGBO.
Figure 15. Experimental and simulated current-voltage I–V (a) and power-voltage (P–V) characteristics (b) for an STM6-40/36 model by MAGBO.
Mathematics 11 04200 g015
Figure 16. Experimental and simulated current data error-index values of an STM6-40/36 model: (a) IAE; (b) RE.
Figure 16. Experimental and simulated current data error-index values of an STM6-40/36 model: (a) IAE; (b) RE.
Mathematics 11 04200 g016
Figure 17. Experimental and simulated current data error-index values of a PVM 752 GaAs thin-film model: (a) IAE; (b) RE.
Figure 17. Experimental and simulated current data error-index values of a PVM 752 GaAs thin-film model: (a) IAE; (b) RE.
Mathematics 11 04200 g017
Figure 18. Experimental and simulated current-voltage I–V (a) and power–voltage (P–V) characteristics (b) for PVM 752 GaAs thin-film model by MAGBO.
Figure 18. Experimental and simulated current-voltage I–V (a) and power–voltage (P–V) characteristics (b) for PVM 752 GaAs thin-film model by MAGBO.
Mathematics 11 04200 g018aMathematics 11 04200 g018b
Table 1. The methodologies of the recent modification/hybridization of GBO.
Table 1. The methodologies of the recent modification/hybridization of GBO.
Algorithm and Ref. MethodologyType of Improvement Problem FormulationStatistical Test
Modified Gradient-Based Optimizer (MGBO) [60]The strategy is to adjust the value of the direction of movement (DM) by introducing a new number D that ranges from 1 to 2 and goes up in increments of 1.ModificationThree Photovoltaics module (single, double-diode), and PV module.none
Improved gradient-based optimizer (IBGO) [61]Two different approaches have been suggested as ways to enhance the functionality of meta-heuristic algorithms. The first tactic employs adaptive weights to move closer to the optimal solution and stay out of local optima during the various phases of the method. The second tactic involves substituting chaotic behavior for the algorithm’s inherent unpredictability to boost the method’s convergence speed and accuracy.ModificationThree Photovoltaics module (single, double-diode), and PV module.none
Opposition decided gradient-based optimizer (OBGBO) [62]OBGBO is an enhancement over GBO, which employs an opposition-based process to produce more precise answers. OBGBO uses the same technique to generate solutions, but it does so in a different way.ModificationThree Photovoltaics module (single, double-diode), and PV module.none
Random learning gradient based optimization (RLGBO) [63]To solve the issue of the GBO being stuck in local optima and to make the convergence process faster and more accurate, a new learning mechanism based on randomization was added to the original version of the GBO. The random learning mechanism allows individuals to maintain constant communication with one another, ultimately resulting in improved performance.Modificationthree various PV models (KC200GT, ST40, and SM55) Wilcoxon ranked test
Sine cosine differential gradient based optimizer (SDGBO) [64]The suggested method utilizes the mutation crossover of the differential evolution algorithm in conjunction with the sine cosine method. With either the crossover operator or the sine cosine approach, the new search agent is encouraged to determine the absolute worst place to settle. By doing so, the algorithm is freed from the confines of whatever local optimums it might have been experiencing.HybridizationFour various PV models, and (KC200GT) Wilcoxon ranked test
Developed Version of Eagle Strategy Gradient-Based Optimizer (ESCGBO) [65]Estimating the parameters of static and dynamic photovoltaic models led to developing an updated version of the GBO algorithm known as ESCGBO. The ESCGBO employs a strategy known as the chaotic eagle. The authors carried out assessments to establish that ESCGBO is successful on both a macro and micro level.Hybridization23 benchmark functions and three various PV modelsNone
Improved moth flame algorithms (IMFOL) [66]IMFOL stands for Improved Moth Flame with Local escape operators, which is the name of the algorithm that has been proposed. The MFO algorithm is enhanced by adding the Local Escape Operators method, which ultimately leads to improved performance in terms of accuracy and efficiency. Additionally, the LEO strategy broadens the demographic range of the population and improves the MFO algorithm’s ability to explore uncharted territories efficiently.Hybridizationfive various PV modelsNone
Table 2. Objective function and variables of different PV models.
Table 2. Objective function and variables of different PV models.
PV ModelObjective FunctionParameters
f ( V ,   I ,   x ) x
SDM I p h I 0 exp V + I R s α V t 1 V + I R S R s h I ( I p h ,   I o , R s , R s h ,   α )
DDM I = I p h I 01 exp V + I R s α 1 V t 1 I 02 exp V + I R s α 2 V t 1 V + I R S R s h I ( I p h ,   I o 1 , I o 2 , R s , R s h ,   α 1 , α 2 )
SMM I = I p h I o exp V + I R s N s α N s V t 1 V + I R s N s R s h N s I ( I p h ,   I o , R s , R s h ,   α )
Table 3. Parameter settings for the competitors.
Table 3. Parameter settings for the competitors.
AlgorithmParameter Setting
GBOeps = 0.005 × 10−3 × random();
EGBOLC = 0.7, eps = 5 × 10−20 × random()
SMAz = 0.03
PSOVmax = 6, wMax = 0.9, wMin = 0.2, c1 = c2 = 2
EOV = 1, a1 = 2, a2 = 1, GP = 0.5
SADEprobability = 50, cr = 5m, crm = 0.5, p1 = 0.5
JADE miu _ f = 0.5 ,   miu _ cr = 0.5 ,   p = 0.1 ,   c = 0.1
CL-PSOc_local = 1.2, w_min = 0.4, w_max = 0.9, max_flag = 7
HPSO_TVACCi = 0.5, cf = 0.0
MFO T = [ - 1 , 1 ] ,   b = 1
LNMHGSVC2 = 0.03, NP = 50
SSA c 1 = [ 2 / e ,   2 ]
MAGBOeps = 0.005 × 10−3 × random();
Table 4. Comparison of MAGBO with basics and advanced algorithms.
Table 4. Comparison of MAGBO with basics and advanced algorithms.
A l g o r i t h m Criteria F1F2F3F4F5F6F7F8F9F10Friedman Avg/Rank
GBOAvg9.217 × 10 2 1.122 × 10 3 1.183 × 10 5 1.907 × 10 3 1.946 × 10 4 2.115 × 10 3 7.997 × 10 3 2.332 × 10 3 2.642 × 10 3 3.216 × 10 3 4.444.000
Std 1.415 × 10 3 6.650 × 10 6.442 × 10 5 0.366 × 10 1.119 × 10 4 5.474 × 10 2 5.292 × 10 3 0.809 × 10 1.117 × 10 2 6.647 × 10
Med 4.134 × 10 2 1.101 × 10 3 7.000 × 10 2 1.906 × 10 3 1.944 × 10 4 2.013 × 10 3 7.122 × 10 3 2.334 × 10 3 2.600 × 10 3 3.213 × 10 3
EGBOAvg 1.499 × 10 2 1.398 × 10 9 1.183 × 10 5 1.908 × 10 3 1.280 × 10 4 3.107 × 10 3 4.270 × 10 4 2.337 × 10 3 2.622 × 10 3 3.208 × 10 3 4.023.000
Std 2.704 × 10 2 7.657 × 10 9 6.442 × 10 5 5.142 × 100 1.259 × 10 4 3.569 × 10 3 9.986 × 10 4 0.4241 × 10 8.419 × 10 8.238 × 10
Med 1.000 × 10 2 1.100 × 10 3 7.000 × 10 2 1.907 × 10 3 7.881 × 10 3 2.341 × 10 3 8.697 × 10 3 2.336 × 10 3 2.600 × 10 3 3.183 × 10 3
SMAAvg 7.679 × 10 3 9.045 × 10 5 9.370 × 10 3 1.904 × 10 3 1.073 × 10 5 2.309 × 10 3 5.924 × 10 4 2.333 × 10 3 2.614 × 10 3 3.160 × 10 3 5.716.000
Std 3.384 × 10 3 7.804 × 10 5 1.662 × 10 4 0.1525 × 10 4.806 × 10 3 8.531 × 10 2 2.737 × 10 4 0.6499 × 10 4.880 × 10 4.418 × 10
Med 8.933 × 10 3 6.705 × 10 5 1.305 × 10 3 1.903 × 10 3 1.053 × 10 5 2.006 × 10 3 5.601 × 10 4 2.300 × 10 3 2.601 × 10 3 3.149 × 10 3
PSOAvg 2.137 × 10 3 2.590 × 10 5 2.991 × 10 4 1.905 × 10 3 5.173 × 10 4 1.788 × 10 3 1.203 × 10 4 2.362 × 10 3 2.715 × 10 3 3.278 × 10 3 4.845.000
Std 2.187 × 10 3 3.880 × 10 5 1.418 × 10 5 0.1161 × 10 2.961 × 10 4 2.248 × 10 2 6.035 × 10 3 5.505 × 10 3.634 × 10 2 1.597 × 10 2
Med 1.453 × 10 3 1.351 × 10 5 7.000 × 10 2 1.903 × 10 3 4.306 × 10 4 1.677 × 10 3 1.063 × 10 4 2.342 × 10 3 2.600 × 10 3 3.232 × 10 3
EOAvg 4.280 × 10 7 6.700 × 10 9 1.897 × 10 9 1.907 × 10 3 8.561 × 10 4 6.070 × 10 3 2.166 × 10 5 2.320 × 10 3 3.014 × 10 3 3.035 × 10 3 6.528.00
Std 1.346 × 10 8 1.369 × 10 10 2.781 × 10 9 1.179 × 10 4.972 × 10 4 2.584 × 10 3 6.129 × 10 5 1.135 × 10 5.436 × 10 2 7.140 × 10
Med 4.233 × 10 3 2.509 × 10 6 5.883 × 10 7 1.903 × 10 3 6.680 × 10 4 5.359 × 10 3 6.654 × 10 4 2.324 × 10 3 2.738 × 10 3 3.017 × 10 3
SADEAvg 1.000 × 10 2 1.100 × 10 3 7.000 × 10 2 1.905 × 10 3 5.156 × 10 4 1.672 × 10 3 2.068 × 10 4 2.318 × 10 3 2.615 × 10 3 3.164 × 10 3 3.452.000
Std 0.000 0.000 0.000 7.25 × 10 1 2.661 × 10 4 1.763 × 10 2 1.239 × 10 4 1.442 × 10 3.009 × 10 4.642 × 10
Med 1.000 × 10 2 1.100 × 10 3 7.000 × 10 2 1.905 × 10 3 4.570 × 10 4 1.604 × 10 3 1.719 × 10 4 2.328 × 10 3 2.600 × 10 3 3.158 × 10 3
CL-PSOAvg 6.292 × 10 6 1.753 × 10 8 1.913 × 10 3 5.718 × 10 4 2.920 × 10 3 2.040 × 10 8 2.328 × 10 3 2.739 × 10 3 3.423 × 10 3 6.709.000
Std 2.358 × 10 7 5.861 × 10 8 3.641 × 10 8 3.100 × 10 2.007 × 10 4 1.306 × 10 3 7.703 × 10 5 0.2165 × 10 3.289 × 10 2 9.664 × 101
Med 7.389 × 10 4 9.251 × 10 5 7.717 × 10 6 1.906 × 10 3 5.562 × 10 4 2.649 × 10 3 1.648 × 10 4 2.328 × 10 3 2.602 × 10 3 3.443 × 10 3
HPSO_TVACAvg 1.754 × 10 4 1.575 × 10 8 1.554 × 10 3 1.926 × 10 3 5.280 × 10 4 2.163 × 10 3 1.678 × 10 4 2.339 × 10 3 2.647 × 10 3 3.298 × 10 3 6.037.000
Std 6.309 × 10 4 8.617 × 10 8 3.196 × 10 3 1.612 × 10 2.411 × 10 4 7.765 × 10 2 8.566 × 10 3 0.8417 × 10 1.382 × 10 2 1.214 × 10 2
Med 2.327 × 10 3 1.646 × 10 4 7.000 × 10 2 1.922 × 10 3 4.808 × 10 4 2.000 × 10 3 1.326 × 10 4 2.338 × 10 3 2.600 × 10 3 3.281 × 10 3
MAGBOAvg 1.000 × 10 2 1.100 × 10 3 7.000 × 10 2 1.905 × 10 3 2.000 × 10 4 1.659 × 10 3 8.714 × 10 3 2.300 × 10 3 2.614 × 10 3 3.237 × 10 3 2.961.000
Std 0.000 0.000 0.000 0.2648 × 10 1.103 × 10 4 1.329 × 10 2 5.893 × 10 3 1.059 × 10 4.880 × 10 7.765 × 10
Med 1.000 × 10 2 1.100 × 10 3 7.000 × 10 2 1.904 × 10 3 2.080 × 10 4 1.606 × 10 3 7.324 × 10 3 2.329 × 10 3 2.601 × 10 3 3.222 × 10 3
Table 5. Wilcoxon rank-sum of the MAGBO vs. basic and advanced recent algorithms on CEC2021.
Table 5. Wilcoxon rank-sum of the MAGBO vs. basic and advanced recent algorithms on CEC2021.
AlgorithmsF1F2F3F4F5F6F7F8F9F10
MAGBO vs. GBO<0.050<0.050<0.050<0.0500.222171<0.0500.680259<0.050<0.0500.266174
MAGBO vs. EGBO<0.050<0.0500.993796<0.050<0.050<0.0500.379592<0.050<0.0500.070029
MAGBO vs. SMA<0.050<0.050<0.0500.137504<0.050<0.050<0.050<0.0501<0.05
MAGBO vs. PSO<0.050<0.0500.80952<0.050<0.050<0.050<0.050<0.050<0.0500.498735
MAGBO vs. EO<0.050<0.050<0.0500.479204<0.050<0.050<0.050<0.050<0.050<0.050
MAGBO vs. SADE1110.331072<0.050<0.050<0.0500.063115<0.050<0.050
MAGBO vs. CL-PSO<0.050<0.050<0.050<0.050<0.050<0.050<0.0500.3795920.114461<0.050
MAGBO vs. HPSO_TVAC<0.050<0.050<0.050<0.050<0.050<0.050<0.050<0.050<0.050<0.050
Table 6. The parameter search ranges.
Table 6. The parameter search ranges.
P a r a m e t e r S i n g l e / D o u b l e   D i o d e P V   M o d u l e S T M 6 - 40 / 36 M o d u l e
L o w e r U p p e r L o w e r U p p e r L o w e r U p p e r
I p h ( A ) 010202
I s d , I s d 1 , I s d 2 ( µ A ) 01050050
R s ( Ω ) 00.50200.36
R s h ( Ω ) 01000200001000
n , n 1 , n 2 12150160
Table 7. Absolute MAGBO error (IAE) on the single-diode SDM.
Table 7. Absolute MAGBO error (IAE) on the single-diode SDM.
Item M e a s u r e d   D a t a S i m u l a t e d   C u r r e n t   D a t a S i m u l a t e d   P o w e r   D a t a
V ( V ) I ( A ) I s i m ( A ) I A E I ( A ) P s i m ( W ) I A E P ( W )
1−0.20570.76400.76408770400.0000877040−0.15717284070.0000180407
2−0.12910.76200.76266308600.0006630860−0.09845980440.0000856044
3−0.05880.76050.76135530700.0008553070−0.04476769210.0000502921
40.00570.76050.76015399100.00034600900.00433287770.0000019723
50.06460.76000.75905520800.00094479200.04903496640.0000610336
60.11850.75900.75804234500.00095765500.08982801790.0001134821
70.16780.75700.75709165300.00009165300.12703997940.0000153794
80.21320.75700.75614136400.00085863600.16120933880.0001830612
90.25450.75550.75508687200.00041312800.19216960890.0001051411
100.29240.75400.75366387800.00033612200.22037131790.0000982821
110.32690.75050.75139096600.00089096600.24562970680.0002912568
120.35850.74650.74735385100.00085385100.26792635560.0003061056
130.38730.73850.74011722200.00161722200.28664740010.0006263501
140.41370.72800.72738222500.00061777500.30091802650.0002555735
150.43730.70650.70697265200.00047265200.30915914070.0002066907
160.45900.67550.67528015300.00021984700.30995359020.0001009098
170.47840.63200.63075827500.00124172500.30175475880.0005940412
180.49600.57300.57192836200.00107163800.28367646760.0005315324
190.51190.49900.49960702400.00060702400.25574883560.0003107356
200.52650.41300.41364879900.00064879900.21778609270.0003415927
210.53980.31650.31751011900.00101011900.17139196220.0005452622
220.55210.21200.21215495100.00015495100.11713074840.0000855484
230.56330.10350.10225132600.00124867400.05759817190.0007033781
240.5736−0.0100−0.00871752400.0012824760−0.00500037180.0007356282
250.5833−0.1230−0.12550739100.0025073910−0.07320846120.0014625612
260.5900−0.2100−0.20847230200.0015276980−0.12299865820.0009013418
Table 8. A comparison of various algorithms applied to the model of a single-diode SDM.
Table 8. A comparison of various algorithms applied to the model of a single-diode SDM.
Algorithm I p h ( A ) I s d ( µ A ) R s ( Ω ) R s h ( Ω ) n RMSE
MAGBO 7.60775530 × 10 7 3.23020770 × 10 7 3.63770933 × 10 2 5.37185214 × 10 0.148118358 × 10 9.860218779287832 × 10 4
PSO 7.60830725 × 10 1 2.97833118 × 10 7 3.67049082 × 10 2 5.15209651 × 10 0.147305453 × 10 9.981852770215333 × 10 4
PGJAYA 7.60777374 × 10 1 3.22628317 × 10 7 3.63818340 × 10 2 5.36737322 × 10 0.148106137 × 10 9.860248538038482 × 10 4
IJAYA 7.60656504 × 10 1 3.29562058 × 10 7 3.62293998 × 10 2 5.37054978 × 10 0.148323279 × 10 9.963586798743776 × 10 4
MLBSA 7.60774976 × 10 1 3.23425701 × 10 7 3.63719322 × 10 2 5.37466051 × 10 0.148130985 × 10 9.860248926943814 × 10 4
CL-PSO 7.60360878 × 10 1 5.44600932 × 10 7 3.48181218 × 10 2 9.60050380 × 10 0.153544041 × 10 1.6631313822752311 × 10 3
HPSO_TVA 7.60711254 × 10 1 3.48500520 × 10 7 3.60688873 × 10 2 5.60201199 × 10 0.148886781 × 10 9.967934554405249 × 10 4
PSO_W 7.60768649 × 10 1 3.24287287 × 10 7 3.63625558 × 10 2 5.38923057 × 10 0.148157655 × 10 9.860565785062050 × 10 4
GBO 7.60776465 × 10 1 3.22245709 × 10 7 3.63871340 × 10 2 5.36607682 × 10 0.148094154 × 10 9.860328319688564 × 10 4
Table 9. Absolute MAGBO error (IAE) on the double-diode DDM.
Table 9. Absolute MAGBO error (IAE) on the double-diode DDM.
Item M e a s u r e d   D a t a S i m u l a t e d   C u r r e n t   D a t a S i m u l a t e d   P o w e r   D a t a
V ( V ) I ( A ) I s i m ( A ) I A E I ( A ) P s i m ( W ) I A E P ( W )
1−0.20570.76400.76398343600.0000165640−0.15715139280.0000034072
2−0.12910.76200.76260411600.0006041160−0.09845219140.0000779914
3−0.05880.76050.76133771400.0008377140−0.04476665760.0000492576
40.00570.76050.76017380000.00032620000.00433299070.0000018593
50.06460.76000.75910769000.00089231000.04903835680.0000576432
60.11850.75900.75812142600.00087857400.08983738900.0001041110
70.16780.75700.75718861700.00018861700.12705624990.0000316499
80.21320.75700.75624360600.00075639400.16123113680.0001612632
90.25450.75550.75517729700.00032270300.19219262210.0000821279
100.29240.75400.75372234500.00027765500.22038841370.0000811863
110.32690.75050.75139912100.00089912100.24563237270.0002939227
120.35850.74650.74730142500.00080142500.26790756090.0002873109
130.38730.73850.74001063900.00151063900.28660612050.0005850705
140.41370.72800.72724693100.00075306900.30086205540.0003115446
150.43730.70650.70685027900.00035027900.30910562700.0001531770
160.45900.67550.67521053000.00028947000.30992163330.0001328667
170.47840.63200.63076075400.00123924600.30175594470.0005928553
180.49600.57300.57199473800.00100526200.28370939000.0004986100
190.51190.49900.49970614700.00070614700.25579957660.0003614766
200.52650.41300.41373368900.00073368900.21783078730.0003862873
210.53980.31650.31754622400.00104622400.17141145170.0005647517
220.55210.21200.21212301600.00012301600.11711311710.0000679171
230.56330.10350.10216329900.00133670100.05754858630.0007529637
240.5736−0.0100−0.00879172100.0012082790−0.00504293120.0006930688
250.5833−0.1230−0.12554339500.0025433950−0.07322946230.0014835623
260.5900−0.2100−0.20837153500.0016284650−0.12293920570.0009607943
Table 10. A comparison of various algorithms applied to the model of a double-diode DDM.
Table 10. A comparison of various algorithms applied to the model of a double-diode DDM.
Algorithm I p h ( A ) I s d 1 ( µ A ) I s d 2 ( µ A ) R s ( Ω ) R s h ( Ω ) n 1 n 2 RMSE
MAGBO7.60781094   ×   10 1 2.25969455   ×   10 7 3.67404535 × 10 2 5.54852834 × 10 0.145101498 × 10 7.49386194 × 10 7 0.200000000 × 10 9.82484851996943 × 10 4
PSO 7.60818668 × 10 1 1.92236886 × 10 7 3.70214586 × 10 2 5.49831144 × 10 0.143704875 × 10 9.69884458 × 10 7 0.199958878 × 10 9.860071193600016 × 10 4
PGJAYA 7.60751928 × 10 1 2.91525289 × 10 7 3.65535114 × 10 2 5.51135412 × 10 0.176996546 × 10 2.27698207 × 10 7 0.145498866 × 10 9.852611722895091 × 10 4
IJAYA 7.60787349 × 10 1 1.66868484 × 10 7 3.65910355 × 10 2 5.36663167 × 10 0.170854351 × 10 2.41944586 × 10 7 0.145998738 × 10 9.860914205884810 × 10 4
MLBSA 7.60780345 × 10 1 2.86499791 × 10 7 3.64993221 × 10 2 5.42564433 × 10 0.147096561 × 10 2.63125416 × 10 7 0.199999784 × 10 9.840561693493248 × 10 4
CL-PSO 7.62937590 × 10 1 1.80424338 × 10 7 3.55917202 × 10 2 4.77950471 × 10 0.144959311 × 10 4.20720089 × 10 7 0.167893978 × 10 1.6430906963831881 × 10 4
HPSO_TVA 7.60769273 × 10 1 3.83641383 × 10 7 3.56754579 × 10 2 5.80624538 × 10 0.149871283 × 10 0.00000000 0.178885865 × 10 1.0402286091084594 × 10 4
PSO_W 7.60770205 × 10 1 3.23613029 × 10 7 3.63706958 × 10 2 5.38189735 × 10 0.148136716 × 10 3.29665828 × 10 12 0.186219871 × 10 9.860319283422937 × 10 4
GBO 7.60777595 × 10 1 3.23396195 × 10 7 3.63727942 × 10 2 5.37457057 × 10 0.148129998 × 10 1.33799831 × 10 12 0.190873826 × 10 9.860264229026566 × 10 4
Table 11. Absolute MAGBO error (IAE) on the PV module model.
Table 11. Absolute MAGBO error (IAE) on the PV module model.
Item M e a s u r e d   D a t a S i m u l a t e d   C u r r e n t   D a t a S i m u l a t e d   P o w e r   D a t a
V ( V ) I ( A ) I s i m ( A ) I A E I ( A ) P s i m ( W ) I A E P ( W )
10.124801.031501.02911916200.00238083800.12843407140.0023808380
21.809301.030001.02738107400.00261892601.85884057720.0026189260
33.351101.026001.02574179800.00025820203.43736333930.0002582020
44.762201.022001.02410715600.00210715604.87700309830.0021071560
56.053801.018001.02229180600.00429180606.18875013520.0042918060
67.236401.015501.01993068200.00443068207.38062638720.0044306820
78.318901.014001.01636310700.00236310708.45502305080.0023631070
89.309701.010001.01049615200.00049615209.40741602630.0004961520
910.216301.003501.00062897100.002871029010.22272575640.0028710290
1011.044900.988000.98454837900.003451621010.87423839120.0034516210
1111.801800.963000.95952167700.003478323011.32408292760.0034783230
1212.492900.925500.92283881800.002661182011.52893306940.0026611820
1313.123100.872500.87259966300.000099663011.45121263750.0000996630
1413.698300.807500.80727426300.000225737011.05828503690.0002257370
1514.222100.726500.72833647700.001836477010.35847420950.0018364770
1614.699500.634500.63713799800.00263799809.36561000160.0026379980
1715.134600.534500.53621306000.00171306008.11537017790.0017130600
1815.531100.427500.42951132100.00201132106.67078327760.0020113210
1915.892900.318500.31877447800.00027447805.06625090140.0002744780
2016.222900.208500.20738950000.00111050003.36445911960.0011105000
2116.524100.101000.09616716400.00483283601.58907583470.0048328360
2216.79870−0.00800−0.00832539500.0003253950−0.13985581300.0003253950
2317.04990−0.11100−0.11093649200.0000635080−1.89145609500.0000635080
2417.27930−0.20900−0.20924727600.0002472760−3.61564645620.0002472760
2517.48850−0.30300−0.30086359800.0021364020−5.26165303360.0021364020
Table 12. A comparison of various algorithms applied to a single-diode DDM model.
Table 12. A comparison of various algorithms applied to a single-diode DDM model.
Algorithm I p h ( A ) I s d ( µ A ) R s ( Ω ) R s h ( Ω ) n RMSE
MAGBO1.03051430   ×   10 3.48226289   ×   10 6 0.120127101 × 10 9.81982284 × 10 2 4.86428348 × 10 2.425074868100019 × 10 3
PSO 0.143559362 × 10 0.00000000 0.00000000 1.46242721 × 10 0.100000000 × 10 2.742507777740578100 × 10 1
PGJAYA 0.103079770 × 10 3.30962024 × 10 6 0.120648712 × 10 9.30403403 × 10 2 4.84489889 × 10 2.4291524524006703 × 10 3
IJAYA 0.103041596 × 10 3.50647703 × 10 6 0.120112652 × 10 1.00002657 × 10 3 4.86687509 × 10 2.4256537082842510 × 10 3
MLBSA 0.103051756 × 10 3.48211561 × 10 6 0.120128062 × 10 9.81736099 × 10 2 4.86426751 × 10 2.4250751724230210 × 10 3
CL-PSO 0.102600504 × 10 4.74880673 × 10 6 0.124732241 × 10 7.81089319 × 10 2 4.99392076 × 10 1.49624138566545101 × 10 2
HPSO_TVA 0.105173347 × 10 5.05107355 × 10 8 0.153633461 × 10 2.05631823 × 10 2 3.64672902 × 10 1.11282483130807390 × 10 2
PSO_W 0.103032474 × 10 3.54993442 × 10 6 0.119961187 × 10 1.01326705 × 10 2 4.87162092 × 10 2.4258587696614703 × 10 3
GBO 0.103066014 × 10 3.42682201 × 10 6 0.120289483 × 10 9.60567724 × 10 2 4.85816052 × 10 2.4255399509107550 × 10 3
Table 13. MAGBO and other competitors for various PV models.
Table 13. MAGBO and other competitors for various PV models.
PV ModuleAlgorithmRMSESignificant Wilcoxon Test Friedman Test/Rank
M i n M a x A v g . S t d
Single-diode model of RTC FranceMAGBO 9.8602 × 10 4 9.8602 × 10 4 9.86022 × 10 4 8.41862 × 10 17 1.2068971
PSO 9.94813 × 10 4 3.8151316 × 10 2 2.34396 × 10 2 5.40938 × 10 2 +8.1034488
PGJAYA 9.86025 × 10 4 1.033627 × 10 3 9.94721 × 10 4 1.26863 × 10 5 +3.4827593
IJAYA 9.96359 × 10 4 1.482295 × 10 3 1.13907 × 10 3 1.48842 × 10 4 +5.5862076
MLBSA 9.86022 × 10 4 1.00735 × 10 3 9.88786 × 10 4 4.85071 × 10 6 +2.2413792
CL-PSO 1.663131 × 10 3 3.606517 × 10 3 2.46597 × 10 3 5.03603 × 10 4 +8.4482769
HPSO_TVA 9.96793 × 10 4 2.170638 × 10 3 1.49135 × 10 3 3.40663 × 10 4 +6.4482767
PSO_W 9.86057 × 10 4 3.8151316 × 10 2 2.47509 × 10 3 6.74863 × 10 3 +5.5172415
GBO 9.86033 × 10 4 3.015336 × 10 3 1.11602 × 10 3 4.22836 × 10 4 +3.9655174
Double-diode model of RTC FranceMAGBO 9.8248 × 10 4 9.9667 × 10 4 9.84729 × 10 4 2.94950 × 10 6 1.6551721
PSO 9.86007 × 10 4 3.8151316 × 10 2 1.23483 × 10 2 1.56943 × 10 2 +8.2413799
PGJAYA 9.85261 × 10 4 1.368534 × 10 3 1.03260 × 10 3 8.03099 × 10 5 +3.5517244
IJAYA 9.86451 × 10 4 1.492106 × 10 3 1.18338 × 10 3 1.73319 × 10 4 +55
MLBSA 9.84056 × 10 4 1.130151 × 10 3 1.00554 × 10 3 4.14601 × 10 5 +2.8275862
CL-PSO 1.643091 × 10 3 4.408675 × 10 3 2.61184 × 10 3 7.56568 × 10 4 +8.1379318
HPSO_TVA 1.040229 × 10 3 2.686124 × 10 3 1.76918 × 10 3 4.45671 × 10 4 +7.1034487
PSO_W 9.86032 × 10 4 0.03339187 2.38006 × 10 3 5.86805 × 10 3 +5.4827596
GBO 9.82733 × 10 4 2.073028 × 10 3 1.06864 × 10 3 2.28955 × 10 4 +33
PV module model of Photo Watt-PWP 201MAGBO 2.42507 × 10 3 2.42507 × 10 3 2.42507 × 10 3 6.15834 × 10 16 11
PSO 2.74250778 × 10 1 7.83911954 × 10 1 3.19336 × 10 1 9.37310 × 10 2 +8.0172419
PGJAYA 2.42507 × 10 3 3.17454 × 10 3 2.46388 × 10 3 1.37087 × 10 4 +2.9655173
IJAYA 2.425654 × 10 3 2.662246 × 10 3 2.46927 × 10 3 5.20336 × 10 5 +3.7586204
MLBSA 2.42507 × 10 3 2.74250778 × 10 1 2.06983 × 10 2 6.89281 × 10 2 +2.6551722
CL-PSO 1.4962413 × 10 2 311203745 × 10 1 1.18115 × 10 1 8.93019 × 10 2 +5.4827585
HPSO_TVA 1.1128248 × 10 2 2.74250778 × 10 1 2.32335 × 10 1 9.53834 × 10 2 +6.7931036
PSO_W 2.425859 × 10 3 4.42058515 × 10 1 2.86844 × 10 1 9.92052 × 10 2 +7.4827588
GBO 2.42554 × 10 3 2.74250778 × 10 1 2.47187 × 10 1 8.2581072 × 10 2 +6.8448277
Table 14. Absolute MAGBO error (IAE) on the STM6-40/36 model.
Table 14. Absolute MAGBO error (IAE) on the STM6-40/36 model.
Item M e a s u r e d   D a t a S i m u l a t e d   C u r r e n t   D a t a S i m u l a t e d   P o w e r   D a t a
V ( V ) I ( A ) I s i m ( A ) I A E I ( A ) P s i m ( W ) I A E P ( W )
101.6631.6634582580.0004582580.00000000000.0000000000
20.1181.6631.663252310.000252310.19626377260.0000297726
32.2371.6611.6595508090.0014491913.71241515970.0032418403
45.4341.6531.6539146990.0009146998.98737247440.0049704744
57.261.651.6505659140.00056591411.98310853560.0041085356
69.681.6451.6454306050.00043060515.92776825640.0041682564
711.591.641.6392335370.00076646318.99871669380.0088833062
812.61.6361.6337126960.00228730420.58477996960.0288200304
913.371.6291.6272858080.00171419221.75681125300.0229187470
1014.091.6191.6183135760.00068642422.80203828580.0096717142
1114.881.5971.6030900450.00609004523.85397986960.0906198696
1215.591.5811.5815883760.00058837624.65696278180.0091727818
1316.41.5421.5423305910.00033059125.29422169240.0054216924
1416.711.5241.5211926340.00280736625.41912891410.0469110859
1516.981.51.4991947440.00080525625.45632675310.0136732469
1617.131.4851.485275270.0002752725.44276537510.0047153751
1717.321.4651.4656542420.00065424225.38513147140.0113314714
1817.911.3881.3875893680.00041063224.85172558090.0073544191
1919.081.1181.1183913760.00039137621.33890745410.0074674541
2021.020−0.0000248140.000024814−0.00052159030.0005215903
Table 15. A comparison of a variety of algorithms applied to the model of an STM6-40/36.
Table 15. A comparison of a variety of algorithms applied to the model of an STM6-40/36.
Algorithm I p h ( A ) I s d ( µ A ) R s ( Ω ) R s h ( Ω ) n RMSE
MAGBO 1.6639078 1.73865682 × 10 6 1.53855769 × 10 1 5.73418580 × 10 2 5.57630767 × 10 1.72981457 × 10 3
SEDE 1.66390478 1.73865696 × 10 6 1.53855761 × 10 1 5.73418588 × 10 2 5.57630770 × 10 1.72981371 × 10 3
MLBSA 1.66390478 1.73865694 × 10 6 1.53855762 × 10 1 5.73418594 × 10 2 5.57630770 × 10 1.72981371 × 10 3
IJAYA 1.66820359 2.71000974 × 10 7 3.53241396 × 10 1 3.85497843 × 10 2 4.91485938 × 10 4.69534212 × 10 3
JADE 1.66390478 1.73870000 × 10 6 4.27377120 × 10 3 1.59282939 × 10 1.52030292 1.72981371 × 10 3
PGJAYA 1.66390506 1.73854219 × 10 6 1.53863151 × 10 1 5.73403658 × 10 2 5.57628110 × 10 1.72981372 × 10 3
SSA 1.86482388 3.13491034 × 10 7 1.35532499 × 10 1 3.31630801 × 10 5.11425353 × 10 1.78578316 × 10 1
MFO 1.86239751 0.00000000 0.00000000 3.28634387 × 10 2.70424691 × 10 3.10757409 × 10 1
LNMHGS [87] 1.66345338 2.08416563 × 10 6 1.32775299 × 10 1 6.08102259 × 10 2 5.65036041 × 10 1.78155207 × 10 3
HSOA [85] 1.66298248 2.51123031 × 10 6 1.10621018 × 10 1 6.49583949 × 10 2 5.72859543 × 10 1.93518826 × 10 3
Table 16. Absolute MAGBO error (IAE) on the PVM 752 GaAs thin-film model.
Table 16. Absolute MAGBO error (IAE) on the PVM 752 GaAs thin-film model.
Item M e a s u r e d   D a t a S i m u l a t e d   C u r r e n t   D a t a S i m u l a t e d   P o w e r   D a t a
V ( V ) I ( A ) I s i m ( A ) I A E I ( A ) P s i m ( W ) I A E P ( W )
1−0.16590.10010.1033130.003213−0.016610.00053298
2−0.12810.10000.1029350.002935−0.012810.00037599
3−0.08880.09990.1025430.002643−0.008870.00023467
4−0.04900.09990.1021450.002245−0.00490.00010999
5−0.01020.09990.1017570.001857−0.001020.00001894
60.02750.09980.101380.001580.0027450.00004345
70.06950.09990.100960.001060.0069430.00007364
80.10610.09980.1005940.0007940.0105890.00008426
90.14600.09980.1001950.0003950.0145710.00005769
100.18280.09970.0998280.0001280.0182250.00002332
110.22300.09970.0994260.0002740.0222330.00006120
120.26000.09960.0990560.0005440.0258960.00014145
130.30010.09970.0986540.0010460.029920.00031383
140.34060.09960.0982490.0013510.0339240.00046004
150.37890.09950.0978660.0016340.0377010.00061913
160.41680.09940.0974860.0019140.041430.00079786
170.45830.09940.0970670.0023330.0455550.00106936
180.49490.09930.0966940.0026060.0491440.00128985
190.53700.09930.0962550.0030450.0533240.00163542
200.57530.09920.0958370.0033630.057070.00193487
210.61230.09900.0953990.0036010.0606180.00220516
220.65460.09880.0948080.0039920.0646740.00261315
230.69180.09830.0941280.0041720.0680040.00288618
240.73180.09770.0930390.0046610.0714970.00341085
250.77020.09630.0912970.0050030.074170.00385325
260.80530.09370.0885590.0051410.0754570.00413981
270.83290.09000.085110.004890.0749610.00407271
280.85500.08550.0811370.0043630.0731030.00373058
290.87380.07990.0767120.0031880.0698170.00278586
300.88870.07430.0723880.0019120.066030.00169898
310.90160.06830.0680570.0002430.0615790.00021924
320.91410.06180.0632240.0014240.0564910.00130153
330.92480.05550.0586130.0031130.0513260.00287915
340.93440.04930.0540890.0047890.0460660.00447473
350.94450.04220.0489380.0067380.0398580.00636414
360.95330.03570.0440510.0083510.0340330.00796108
370.96180.02910.0389870.0098870.0279880.00950938
380.97020.02220.0336850.0114850.0215380.01114278
390.97780.01570.0286010.0129010.0153510.01261445
400.98520.00920.0233470.0141470.0090640.01393780
410.99260.00260.0177380.0151380.0025810.01502608
420.99990.00400.0054650.0014650.0040.00146486
431.00460.00850.006810.0153070.0085390.01537696
441.00890.01240.01930.0317030.012510.03198466
Table 17. A comparison of various algorithms applied on the model of a PVM 752 GaAs thin-film.
Table 17. A comparison of various algorithms applied on the model of a PVM 752 GaAs thin-film.
Algorithm I p h ( A ) I s d ( µ A ) R s ( Ω ) R s h ( Ω ) n RMSE
MAGBO 1.021541340 × 10 1 4.912951570 × 10 10 5.000000000 × 10 1 1.000000000 × 10 2 2.000000000 7.498895770 × 10 3
SEDE 1.021541340 × 10 1 4.912951570 × 10 10 5.000000000 × 10 1 1.000000000 × 10 2 2.000000000 7.498895770 × 10 3
MLBSA 1.021541340 × 10 1 4.912951570 × 10 10 5.000000000 × 10 1 1.000000000 × 10 2 2.000000000 7.498928370 × 10 3
IJAYA 1.021541340 × 10 1 4.912951570 × 10 10 5.000000000 × 10 1 1.000000000 × 10 2 2.000000000 7.498895770 × 10 3
JADE 1.021541340 × 10 1 4.912951570 × 10 10 5.000000000 × 10 1 1.000000000 × 10 2 2.000000000 7.498895770 × 10 3
PGJAYA 1.021541340 × 10 1 4.912951560 × 10 10 5.000000000 × 10 1 1.000000000 × 10 2 2.000000000 7.498895770 × 10 3
SSA 1.025574810 × 10 1 3.993709460 × 10 10 1.706274330 × 10 2 6.503079280 × 10 1.972598610 9.017017310 × 10 3
MFO 1.021541340 × 10 1 4.912951570 × 10 10 5.000000000 × 10 1 1.000000000 × 10 2 2.000000000 7.498895770 × 10 3
LNMHGS [87] 1.021541340 × 10 1 4.912951560 × 10 10 5.000000000 × 10 1 1.000000000 × 10 2 2.000000000 7.498895770 × 10 3
HSOA [85] 1.021541340 × 10 1 4.912951570 × 10 10 5.000000000 × 10 1 1.000000000 × 10 2 2.000000000 7.498895770 × 10 3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mahmood, B.S.; Hussein, N.K.; Aljohani, M.; Qaraad, M. A Modified Gradient Search Rule Based on the Quasi-Newton Method and a New Local Search Technique to Improve the Gradient-Based Algorithm: Solar Photovoltaic Parameter Extraction. Mathematics 2023, 11, 4200. https://doi.org/10.3390/math11194200

AMA Style

Mahmood BS, Hussein NK, Aljohani M, Qaraad M. A Modified Gradient Search Rule Based on the Quasi-Newton Method and a New Local Search Technique to Improve the Gradient-Based Algorithm: Solar Photovoltaic Parameter Extraction. Mathematics. 2023; 11(19):4200. https://doi.org/10.3390/math11194200

Chicago/Turabian Style

Mahmood, Bushra Shakir, Nazar K. Hussein, Mansourah Aljohani, and Mohammed Qaraad. 2023. "A Modified Gradient Search Rule Based on the Quasi-Newton Method and a New Local Search Technique to Improve the Gradient-Based Algorithm: Solar Photovoltaic Parameter Extraction" Mathematics 11, no. 19: 4200. https://doi.org/10.3390/math11194200

APA Style

Mahmood, B. S., Hussein, N. K., Aljohani, M., & Qaraad, M. (2023). A Modified Gradient Search Rule Based on the Quasi-Newton Method and a New Local Search Technique to Improve the Gradient-Based Algorithm: Solar Photovoltaic Parameter Extraction. Mathematics, 11(19), 4200. https://doi.org/10.3390/math11194200

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