Next Article in Journal
A Review of Raman-Based Technologies for Bacterial Identification and Antimicrobial Susceptibility Testing
Next Article in Special Issue
Enhanced Absorption in InP Nanodisk Arrays on Ultra-Thin-Film Silicon for Solar Cell Applications
Previous Article in Journal
All-Optical Phase-Change Memory with Improved Performance by Plasmonic Effect
Previous Article in Special Issue
Passivation Effect of CsPbI3 Quantum Dots on the Performance and Stability of Perovskite Solar Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Chaotic Second Order Oscillation JAYA Algorithm for Parameter Extraction of Photovoltaic Models

School of Optical Electrical and Computer Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China
*
Author to whom correspondence should be addressed.
Photonics 2022, 9(3), 131; https://doi.org/10.3390/photonics9030131
Submission received: 20 January 2022 / Revised: 21 February 2022 / Accepted: 22 February 2022 / Published: 25 February 2022
(This article belongs to the Special Issue Advances in Photovoltaic Technologies from Atomic to Device Scale)

Abstract

:
In order to identify the parameters of photovoltaic (PV) cells and modules more accurately, reliably and efficiently, a chaotic second order oscillation JAYA algorithm (CSOOJAYA) is proposed. Firstly, both logical chaotic map and the mutation strategy are brought in enhancing the population diversity and improving exploitation. Secondly, in order to better balance exploration and exploitation, the second order oscillation factor is added, which not only improves the diversity of the population, but also has strong exploration at the beginning of the iteration and strong exploitation at the end of the iteration. In order to balance the two abilities, the self-adaptive weight is introduced into the CSOOJAYA algorithm to regulate individuals the tendency of moving toward the optimal solution and escaping from the worst solution, so as to enhance the search efficiency and exploitation. In order to validate the behavior of CSOOJAYA, it is employed to the parameter identification problem of PV models. Finally, the experimental results show that CSOOJAYA delivers excellent behavior in the aspects of convergence, reliability and accuracy.

1. Introduction

In order to alleviate the greenhouse effect and solve problems such as global warming, the application of renewable energy has attracted much attention [1]. Among the renewable energy sources, solar energy is considered as one of the most promising renewable energy source because of its green, clean, widely distributed and free energy [2]. Solar PV systems can directly convert solar energy into electrical energy, so it has become one of the most widely used power generation technologies [3]. However, PV systems are susceptible to external factors, especially irradiance and temperature [4]. PV arrays tend to degrade over time under harsh environmental conditions, significantly affecting the behavior and utilization efficiency of PV systems [5]. Therefore, in order to accurately control PV systems, it is significant to adopt the experimentally measured current-voltage (I-V) data to estimate the actual behavior of PV systems at work [6]. The single diode model (SDM) and the double diode model (DDM) are widely employed and can accurately reflect the nonlinear behavior of PV panels [7]. The parameter accuracy of these models is essential for modeling, control and optimization of PV panels [8]. However, these parameters are susceptible to various environmental factors and become unavailable [9]. PV models are implicitly nonlinear transcendental equations [10]. It is very significant and challenging to employ an efficient and reliable methods to identify the PV parameters.
Many methods have been proposed to solve PV parameter identification. In principle, these methods can be divided into analytical methods and optimization methods [11]. The analytical methods mainly employ several key data (short-circuit current, open-circuit voltage, and maximum power) provided by the manufacturer, and use mathematical equations to deduce the parameters [12]. The analytical methods have the advantages of simple, fast and direct calculation. However, these data are obtained under standard test conditions [13]. Therefore, the extracted parameters above cannot precisely predict the I-V curves at different temperatures and solar irradiance [14]. Furthermore, this approach depends heavily on the key points of the selected I-V characteristics. If the crucial points are wrongly selected, the computational accuracy will be significantly reduced [11].
Instead of choosing a few key points, the optimization method takes into account all the actual measured current and voltage data. From the algorithm point of view, optimization methods can be divided into meta-heuristic methods and deterministic methods [15]. The optimization methods transform the PV parameter identification problem into an optimization problem, and then the parameters are identified according to all the reference points on the I-V characteristic curve. Deterministic methods, such as Newton–Raphson method [16,17], iterative curve fitting [18,19] and Lambert W-function [20,21], require the objective function to be continuous, convex, and differentiable [22]. The I-V curve of the PV cell model is non-linear, multi-peak and multi-mode. And the deterministic method is easily affected by the initial conditions and gradient information, so it is free to become trapped in a local optimum when resolving such complicated problems [23].
The heuristic methods have no strict requirements on the form of the optimization problem and are not affected by initial conditions and gradient information [24]. Therefore, many heuristic methods have been used to identify the parameters of PV models in the past ten years. These include the classified perturbation mutation based particle swarm optimization algorithm (CPMPSO) [1], efficient teaching-learning-based optimization algorithm (MTLBO) [16], multiswarm spiral leader particle swarm optimization algorithm (MSLPSO) [18], enhanced adaptive differential evolution algorithm (EJADE) [22], memetic adaptive differential evolution (MADE) [23], modified Rao-1 optimization algorithm (MRAO-1) [25], improved gaining-sharing knowledge algorithm (IGSK) [26] and chaos induced coyote algorithm (CICA) [27]. The characteristics of the three existing methods are summarized in Table 1.
Compared with deterministic methods, heuristic methods can acquire more precise and robust PV parameter identification results. However, PV models are nonlinear and multi-mode, so some heuristic methods may become trapped in local optima during the optimization process [8]. In addition to the basic parameters of population size and termination condition, most heuristic algorithms have specific parameters related to their own mechanisms. The behavior of the heuristic algorithm depends to some extent on these specific parameters. Choosing inappropriate parameters will not only increase the computational burden of the algorithm, but it may also converge prematurely or fall into a local optimum [28].
Rao proposed a new and efficient heuristic algorithm in 2016, namely the JAYA algorithm [29]. The structure of the JAYA algorithm is relatively simple, with only two parameters: population size and termination conditions [28]. Therefore, the advantage of the JAYA algorithm is that it can reduce the time of the optimization process and avoid the difficulties caused by parameter adjustment. However, the individual position is only affected by the current optimal individual and the worst individual in the JAYA algorithm, which cannot effectively maintain the diversity of the population [8]. Therefore, JAYA algorithm is free to trap in local optimum when solving complicated problems such as multi-peak and multi-mode.
Although the performance of some improved JAYA algorithms such as the performance-guided JAYA algorithm or PGJAYA [8], comprehensive learning JAYA algorithm (CLJAYA) [10], improved JAYA algorithm (IJAYA) [28], logistic chaotic JAYA algorithm (LCJAYA) [29], and elite opposition-based JAYA algorithm (EO-JAYA) [30] have been improved to some extent, the imbalance between exploration and exploitation is not fully considered. The corresponding summary of the extraction parameters of existing heuristic methods is shown in Table 2.
In order to improve the performance of JAYA algorithm, this paper proposes a chaotic second order oscillation JAYA algorithm (CSOOJAYA), which can identify parameters more accurately and stably. Firstly, the introduction of logical chaotic mapping mechanism improves the population diversity and exploration. Secondly, the introduction of the second order oscillation mechanism not only improves the diversity of the population, but also has strong exploration ability due to the oscillation convergence of the solution vector in the early iteration. The solution vector converges asymptotically in the later iteration and has strong exploitation ability. The self-adaptive weight mechanism is introduced to adjust the tendency of individuals to approach the optimal solution and escape the worst solution, improve the search efficiency of the population and the exploitation. The algorithm improves and balances the global search ability and local optimization ability as a whole. The introduction of mutation mechanisms ensures that individuals avoid getting stuck in local optima. The structure of CSOOJAYA is similar to the JAYA algorithm, and it only has two parameters: population size and termination condition. In order to verify the effectiveness of the proposed CSOOJAYA algorithm, CSOOJAYA and other novel algorithms are used to extract the parameters of three PV models. The experimental results show that CSOOJAYA has the best performance in terms of identification accuracy, stability and convergence speed.

2. PV Models and Objective Function

2.1. Formula of SDM

Figure 1 presents the circuit diagram of the SDM. Obviously, the framework of the SDM is relatively simple. There are photogenerated current I p h , output current I L , diode current I d , and shunt current I s h in the circuit. According to Kirchhoff current law, I L can be get by Equation (1). According to the diode Shockley equation and Ohm law, I d can be get by Equation (2), and I s h can be obtained by Equation (3):
I L   =   I p h I d I s h
I d = I s d · exp q · V L + R S · I L n · k · T 1  
I s h = V L + R S · I L R s h
I L = I p h I d · exp q · V L + R S · I L n · k · T 1 V L + R S · I L R s h
R S refers to the series resistance, R s h stands for the parallel resistance, V L stands for the output voltage, n is the ideality factor of the diode. k refers to the Boltzmann constant (1.3806503 × 10 23 J/K), q stands for the elementary charge (1.60217646 × 10 19 C), and T refers to the absolute temperature. Therefore, we can substitute Equations (2) and (3) into Equation (1) to get Equation (4).
Obviously, SDM needs to identify five unknown different parameters ( I p h , I d , R S , R s h , n).

2.2. Formula of DDM

Taking into account the inherent drawbacks of SDM, DDM can more accurately illustrate the relationship between voltage and current. Figure 2 presents the circuit diagram of the DDM. It can be seen from the Figure 2, DDM differs from SDM in that it has one more diode in parallel with the current source than SDM. I L can be obtained by Equation (5):
I L = I p h I d 1 · exp q · V L + R S · I L n 1 · k · T 1 I d 2 · exp q · V L + R S · I L n 2 · k · T 1 V L + R S · I L R s h
I d 1 refers to the diode diffusion current, and I d 2 stands for the saturation current. n 1 and n 2 represent the ideal saturation factor of two diodes, respectively. Obviously, the DDM needs to extract seven unknown different parameters ( I p h , I d 1 , I d 2 , R S , R s h , n 1 , n 2 ).

2.3. Formula of PV Module

Figure 3 presents the circuit diagram of a PV module, which consists of several PV cells in parallel or in series. I L can be obtained by Equation (6):
I L = N P · I p h N P · I d · exp q · V L / N S + R S · I L / N P n · k · T 1 N P · V L / N S + R S · I L R s h
N S and N P refer to the number of serial and parallel connections of PV cells, respectively. This paper adopts the SDM of PV modules. Therefore, the PV module needs to identify five unknown different parameters ( I p h , I d , R S , R s h , n).

2.4. Objective Function

An objective function is determined to quantitatively evaluate the difference between the actual measured value and the calculated value. Equations (7)–(9) represent the error functions for the experimental and calculated data points of the SDM, DDM, and PV module, respectively:
f S D V L , I L , X = I L I p h + I d · exp q · V L + R S · I L n · k · T 1 + V L + R S · I L R s h   X = I p h , I d , R S , R s h , n
f D D V L , I L , X = I L I p h + I d 1 · exp q · V L + R S · I L n 1 · k · T 1 + I d 2 · exp q · V L + R S · I L n 2 · k · T 1 + V L + R S · I L R s h X = I p h , I d 1 , I d 2 , R S , R s h , n 1 , n 2
f M D V L , I L , X = I L N P · I p h + N P · I d · exp q · V L / N S + R S · I L / N P n · k · T 1 + N P · V L / N S + R S · I L R s h X = I p h , I d , R S , R s h , n
RMSE X = 1 N k = 1 N f k ( V L , I L , X ) 2
However, the error function cannot reflect all the differences between the calculated and experimental data. The root mean square error (RMSE) defined by Equation (10) is used as the objective function to quantitatively evaluate the overall difference between the experimental and calculated data. References like [10,11,14] generally use the objective function. X represents the solution composed of different unknown parameters, and N represents the number of actual measured data.

3. JAYA Algorithm

JAYA is a new population-based intelligent optimization algorithm [29]. In each generation of the JAYA algorithm, the solution vector is optimized by approaching the optimal solution while staying away from the worst solution. JAYA does not need to tune algorithm-specific parameters. The only two parameters are population size and maximum number of evaluation functions ( m a x F E S ).
For the objective function with d unknown variables (j = 1,2,3...d), x i , j represent the value of the j-th variable of the i-th unknown solution, then X i = ( x i , 1 , x i , 2 ...... x i , d ) represents the i-th unknown solution. If an individual can get the worst/best value of f(X) among all individuals, then this is the worst/best individual, expressed as x b e s t   = ( x b e s t , 1 , x b e s t , 2 ,… x b e s t , d )/ x w o r s t   = ( x w o r s t , 1 , x w o r s t , 2 ,... x w o r s t , d ). Each solution x i , j of the JAYA algorithm can be updated by Equation (11):
x n e w i , j = x i , j + r a n d 1 · x b e s t , j x i , j r a n d 2 · x w o r s t , j x i , j
x i , j represents the absolute value of x i , j , and x n e w i , j represents the updated value of x i , j . r a n d 1 and r a n d 2 represent random numbers in [0,1]. The second and third terms in the above Equation (11) represent the tendency of the solution vector x to move towards the optimal solution and escape from the worst solution, respectively. In order to retain a better solution vector, the updated solution x n e w i = ( x n e w i , 1 , x n e w i , 2 ,…, x n e w i , d ) is received if it can provide a better function value.

4. Chaotic Second Order Oscillation JAYA Algorithm

4.1. Motivation

At present, the identification of PV parameters is considered to be a multi-peak and multi-modal problem, that is, there are multiple local optimal solutions. When the JAYA algorithm deals with such problems, it easily falls into a local optimum because it cannot effectively explore different regions of the space. How to effectively enhance the behavior of the algorithm while balancing the exploitation and exploration? Therefore, a chaotic second order oscillation JAYA (CSOOJAYA) algorithm is proposed in this research.

4.2. Logistic Chaotic Map Strategy

In the optimization process of the JAYA algorithm, it is hard to use two uniformly distributed random numbers to maintain the population diversity. However, the number of chaotic sequences generated by the logical chaotic mapping strategy is aperiodic and does not converge to a specific value. It enhances the population diversity of the heuristic algorithm, improves the global search ability, and avoids falling into the local optimum. Therefore, two random numbers are replaced by chaotic sequence numbers, and Equation (12) represents the new update solution:
x n e w i , j = x i , j + C 1 i , j · x b e s t , j x i , j C 2 i , j · x w o r s t , j x i , j
where C 1 i , j and C 2 i , j represent two chaotic sequence numbers, which are generated by Equation (13):
C n + 1 = 4 C n · 1 C n
C n represents the chaotic value of the n-th iteration. Its value range is [0,1], and C 1 is usually set to 0.8.

4.3. Second Order Oscillation Strategy

According to the update equation of the JAYA algorithm, When the i-th solution is the best/worst solution, one term on the right side of Equation (11) is 0, which cannot maintain the diversity of the population. There are only two factors that affect the population position of the JAYA algorithm, namely the current optimal individual and the worst individual. Although the algorithm can speed up the convergence speed and improve the exploitation, it cannot effectively enhance the diversity of the population and weaken the exploration ability in the optimization process.
Considering the above two situations, the population individuals may fall into a local optimum in the JAYA algorithm. It can be seen from Equations (11) and (12) that x best and x worst have a significant impact on the performance of JAYA. Therefore, taking advantages of population information and enhancing exploration capabilities are the available methods to improve the behavior of the JAYA. How to balance exploitation and exploration while improving population diversity. Therefore, the second order oscillation factor is introduced into Equation (12), that is, two unequal second order oscillation factors are added respectively to the two rightmost terms of Equation (12), as shown in Equation (14). The difference between the current optimal individual and the optimal individual of the previous generation and the difference between the current worst individual and the worst individual of the previous generation can guide the individual to update a better range, maintain population diversity and improve exploration ability. After adopting the second order oscillation strategy, because the solution vector has oscillation convergence in the early iteration, the algorithm has strong global search ability; The solution vector converges asymptotically in the later stage of iteration, and the algorithm has strong local optimization ability. This strategy improves the convergence accuracy of the algorithm as a whole:
x n e w i , j = x i , j + C 1 i , j · 1 + k 1 · x b e s t , j k 1 · x b e s t p , j x i , j C 2 i , j · 1 + k 2 ·   x w o r s t , j k 2 · x w o r s t p , j x i , j
where x b e s t p , j and x w o r s t p , j represent the optimal solution and the worst solution of the jth variable of the previous generation of population individuals, respectively. k 1 and k 2 are two randomly generated numbers in [0,1], and k 1 is not equal to k 2 .

4.4. Self-Adaptive Weight

Equation (12) enhances the population diversity and exploration, and Equation (14) not only improves exploitation and exploration, but also effectively balances the two abilities of the algorithm. However, how to balance the two abilities as a whole is the core problem of the algorithm. Considering that the strategy is dedicated to improving exploitation ability, the self-adaptive weight strategy is introduced into Equation (14) to obtain Equation (15). w is obtained from Equation (16). It can be seen from Equation (16) that the value of w increases gradually with the increase of the number of function evaluations. Because the difference between the optimal solution and the worst solution becomes smaller and smaller as the search process progresses. Therefore, in the optimization process of the JAYA algorithm, the population can approach the promising solution at the beginning of the iteration and conduct local optimization in the area of the promising solution at the end of the iteration. Thus, it enhances the search efficiency of the population and the exploitation of the algorithm:
x n e w i , j = x i , j + C 1 i , j · 1 + k 1 · x b e s t , j k 1 · x b e s t p , j x i , j w · C 2 i , j · 1 + k 2 ·     x w o r s t , j k 2 · x w o r s t p , j x i , j  
w = ( f x b e s t f x w o r s t ) 2 ,   i f   f x w o r s t 0 1   ,   o t h e r w i s e
f x b e s t and f x w o r s t represent the values of the optimal and worst individual of the objective function, respectively. Furthermore, the self-adaptive weight is automatically determined, so there is no need to adjust parameters.

4.5. Mutation Mechanism

The three position update equations based on the JAYA algorithm not only improve the overall exploration and exploitation, but also effectively balance the two abilities. However, the optimal individual may still have a very small probability of being in a local optimal state to solve the identification of PV parameters. In order to make the individual jump out of the local optimum, the mutation mechanism is proposed. If from the first function evaluation to a quarter of the maximum function evaluation time, and so on, if the RMSE values of the two are equal, the population individual is reinitialized. And if the best RMSE of the previous generation is better than the current optimal RMSE, the RMSE of the previous generation is accepted.

4.6. Framework of CSOOJAYA

The overall update equation of the improved algorithm based on the above strategy is shown in Equation (17). p is the probability, which represents random numbers in [0,1]:
x n e w i , j = { x i , j + C 1 i , j · ( ( 1 + k 1 ) · x b e s t , j k 1 · x b e s t p , j | x i , j | ) C 2 i , j · ( ( 1 + k 1 ) · x w o r s t , j k 2 · x w o r s t p , j | x i , j | ) , i f   0 p < 1 / 2 x i , j + C 1 i , j · ( ( 1 + k 1 ) · x b e s t , j k 1 · x b e s t p , j | x i , j | ) C 2 i , j · w · ( ( 1 + k 2 ) · x w o r s t , j k 2 · x w o r s t p , j | x i , j | ) , i f 1 2 p < 3 / 4 x i , j + C 1 i , j · ( x b e s t , j | x i , j | ) C 2 i , j · ( x w o r s t , j | x i , j | ) , i f   3 4 p 1
The pseudocode of the CSOOJAYA algorithm can be summarized by the Algorithm 1. In addition, Figure 4 shows the flow chart of the specific implementation of CSOOJAYA. Similar to the JAYA algorithm, the CSOOJAYA does not need to adjust the specific algorithm parameters.
Algorithm 1: CSOOJAYA algorithm.
  • Set the population size (NP) and m a x F E S ;
  • Randomly initialize the population and RMSE values of all individuals are calculated;
  • FES = NP;
  • While FES < m a x F E S do
  •   Identify the best xi and worst xi of all individuals;
  •     For i = 1 to NP do
  •       if (0 ≤ p < 1/2) then
  •         Update the xi using Equation (14);
  •       else if (1/2 ≤ p < 3/4) then
  •         Calculate the w using Equation (15);
  •         Update the xi using Equation (16)
  •       else
  •         Update the xi using Equation (12)
  •       End if
  •        Obtain the new population and Calculate the RMSE for the updated individual x n e w i ;
  •       Receive the updated solution if it is more excellent than the pre-update one
  •     End For
  •     If the conditions of the mutation strategy are met
  •      Randomly initialize the population and Save the current optimal solution and RMSE
  •      End if
  • End While

5. Experimental Results and Discussion

In order to estimate the effective behavior of CSOOJAYA, it was applied to the identification of PV parameters. The benchmark data are extracted from reference [11]. where the data for SDM and DDM were obtained from a commercial RTC silicon PV cell (1000 w/m2 at 33 °C, Photowatt, Bourgoin-Jallieu, France) with a diameter of 57 mm, and the data for the PV module was obtained from a Photowatt-PWP201 (1000 w/m2 at 45 °C) solar module consisting of 36 polycrystalline silicon cells connected in series (Photowatt, Bourgoin-Jallieu, France). In order to ensure fairness, the search space of the solution vector is the same, and the value range of the identification parameter is the same as that adopted in the previous literature. Table 3 shows the value ranges of the parameters corresponding to different PV models.
CSOOJAYA is compared with other novel algorithms to verify its superior performance. These algorithms are DE/WOA, MLBSA, EJADE, MTLBO, MPPCEDE, MADE, ITLBO, CPMPSO, DERAO, JAYA, IJAYA, LCJAYA and PGJAYA. Because these optimization algorithms have better performance in PV model parameter identification, they are used for comparison with CSOOJAYA. Table 4 presents the corresponding parameter settings for the selected algorithms, which are extracted from their respective references. For fair comparison, set m a x F E S to 50,000 for all algorithms and run 30 times independently. The results of DE/WOA [11], MLBSA [5], EJADE [22], MTLBO [16], MPPCEDE [15], MADE [23], ITLBO [2], CPMPSO [1], DERAO [25], JAYA [31], IJAYA [28], LCJAYA [29] and PGJAYA [8] were obtained directly from the corresponding references. The selected algorithms are all executed in MATLAB R2016b on a PC configured with an Intel(I) Core (TM) i5-7200u CPU operating at 3.1 GHz and equipped with 4 Gb of RAM.

5.1. Results on the SDM

For SDM, Table 5 illustrates the best parameters and RMSE obtained for all selected algorithms. The overall best RMSE values of all compared algorithms are indicated in bold. It is obvious from Table 5 that CSOOJAYA, DE/WOA, MLBSA, EJADE, MTLBO, MPPCEDE, MADE, ITLBO, CPMPSO, DERAO, LCJAYA and PGJAYA obtained the smallest RMSE (9.8602 × 10−4), followed by IJAYA. Although accurate parameter value information cannot be obtained, it can be seen from the literature [1,2] that the RMSE can be employed to represent the accuracy. Furthermore, the smaller the RMSE, the more precise the identified parameters.
To further confirm the performance of CSOOJAYA, the IAE of power and current between the measured data and the calculated data was used. The I A E c u r r e n t can be calculated by Equation (18), and the I A E p o w e r can be obtained by Equation (19):
I A E c u r r e n t = I m e a s u r e d I e s t i m a t e d
I A E p o w e r = V I m e a s u r e d V I e s t i m a t e d
where I e s t i m a t e d represents the output current estimated by different models through Equations (4), (5) or (6), and I m e a s u r e d represents the actual measured current. V represents the actual measured voltage. I A E c u r r e n t represents the individual absolute error of current, and I A E p o w e r represents the individual error of power.
The IAE of power and current are shown in Table 6. All values of I A E c u r r e n t are less than 2.50 × 10−3 and all values of I A E p o w e r are less than 1.46 × 10−3, proving the accuracy of extracted parameters. To further verify the accuracy of the results, Figure 5 illustrates the best results obtained using CSOOJAYA to plot the I-V and P-V curves. Obviously, the data calculated by CSOOJAYA are in good agreement with the measured data across the entire voltage range.

5.2. Results on the DDM

For DDM, the extraction of seven parameters increases the difficulty of parameter identification. Table 7 illustrates the best parameters and RMSE obtained for all selected algorithms. The overall best RMSE value of all compared algorithms are indicated in bold. It is obvious from Table 7 that CSOOJAYA, DE/WOA, EJADE, MTLBO, ITLBO, CPMPSO, DERAO have the smallest RMSE (9.8248 × 10−4), followed by MLBSA and MPPCEDE. To further verify the accuracy of the results, Figure 6 illustrates the best results obtained using CSOOJAYA to plot the I-V and P-V curves. Obviously, the data calculated by CSOOJAYA are in good agreement with the measured data across the entire voltage range.
The IAE of power and current are shown in Table 8. All values of I A E c u r r e n t are less than 2.54 × 10−3 and all values of I A E p o w e r are less than 1.48 × 10−3, which demonstrates the accuracy of CSOOJAYA parameter identification.

5.3. Results on the PV Module

The effectiveness of CSOOJAYA is further verified using the PhotoWatt-PWP201 PV module. Table 9 illustrates the best parameters and RMSE obtained for all selected algorithms. The overall best RMSE value of all the compared algorithms are marked in bold. Among all the comparison algorithms, CSOOJAYA, DE/WOA, MLBSA, EJADE, MTLBO, MPPCEDE, MADE, ITLBO, CPMPSO, PGJAYA, IJAYA, LCJAYA and DERAO obtained the best RMSE value, followed by JAYA. To further verify the accuracy of the results, Figure 7 illustrates the best results obtained using CSOOJAYA to plot the I-V and P-V curves. Obviously, the data calculated by CSOOJAYA are in good agreement with the measured data across the entire voltage range. The IAE of power and current are shown in Table 10. All values of I A E c u r r e n t are less than 4.83 × 10−3 and all values of I A E p o w e r are less than 7.99 × 10−2, which verifies the accuracy of CSOOJAYA parameter identification again.

5.4. Statistical Results and Convergence Curve

The above three sections demonstrate the superior accuracy of CSOOJAYA in solving parameter identification on different PV models. However, convergence and robustness should also be considered. Therefore, this section compares the convergence curves and statistical results of the three PV models.
The maximum RMSE (Max), minimum RMSE(Min), average RMSE (Mean) and standard deviation (SD) of RMSE obtained by running all the comparison algorithms 30 times independently on three different PV models are shown in Table 11. The overall best RMSE values for all compared algorithms are marked in bold. where Max represents the worst accuracy, Min is the best accuracy, and Mean is the average accuracy, and SD represents the stability of algorithm. It can be seen from the Table 11:
(1)
For the best RMES value, only CSOOJAYA, DE/WOA, EJADE, MTLBO, ITLBO, CPMPSO and DERAO get the best RMSE values for all PV models.
(2)
For the worst RMSE value, only CSOOJAYA and DE/WOA can obtain the optimal worst value of all PV models
(3)
For the average RMSE value, only CSOOJAYA can obtain the optimal average value for all PV models.
(4)
Taking into account the SD, it can reflect the reliability of the algorithm. Although the SD of MTLBO, MPPCEDE, ITLBO, CPMPSO and DERAO are slightly better than those of the CSOOJAYA algorithm in the SDM and the PV module, the SD of CSOOJAYA is far superior to other algorithms in the DDM. On the whole, the CSOOJAYA algorithm has stronger reliability than other algorithms. Based on the superior stability and accuracy of three PV models, CSOOJAYA can obtain the best parameter identification results.
In order to validate the convergence speed of CSOOJAYA, Figure 8, Figure 9 and Figure 10 plot the iterative curves of several algorithms on the PV models. Obviously, CSOOJAYA not only get the best results among these three models, but also converges faster than other algorithms.
In summary, the above comparison shows that CSOOJAYA has faster convergence, better accuracy and robustness in extracting the parameters of PV models. Furthermore, the behavior of CSOOJAYA is competitive when compared to other algorithms.

6. Conclusions

This paper proposes a chaotic second order oscillation JAYA algorithm to provide an accurate and stable way to identify the PV parameters. The algorithm introduces the logical chaotic map strategy, and the generated random sequence numbers are non-periodic and non-convergent, which can enhance the diversity of the population and improve the exploration. In order to improve and balance the above two abilities, the second order oscillation strategy is introduced into the algorithm. The strategy can not only maintain the diversity of the population, but also has strong exploration ability due to the oscillation convergence of the solution vector in the early iteration; The solution vector converges asymptotically in the later iteration and has strong exploitation ability. In order to balance the exploitation and exploration as a whole, the self-adaptive weight mechanism is introduced to adjust the tendency of individuals to approach the optimal solution and escape the worst solution and improve the search efficiency of the population and the exploitation. The introduction of mutation mechanisms ensures that individuals avoid getting stuck in local optima. In order to validate the behavior of CSOOJAYA, which is employed to the parameter identification problem of PV models. Finally, the experimental results present that CSOOJAYA delivers excellent behavior in the aspects of convergence, reliability and accuracy.
The data of the PV cells used in this article are measured under certain temperature and light irradiance conditions. In order to further verify the applicability of CSOOJAYA algorithm to parameter identification under different PV cell or module operating scenarios. The next step of the research group will be to study the parameter identification of PV cells or components in multiple scenarios.

Author Contributions

Conceptualization, X.J.; methodology, X.J.; software, Y.C.; validation, Y.C.; writing—original draft preparation, X.J.; writing—review and editing, Y.C.; supervision, X.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (NSFC), grant number 11774017.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liang, J.; Ge, S.; Qu, B. Classified perturbation mutation based particle swarm optimization algorithm for parameters extraction of photovoltaic models. Energy Convers. Manag. 2020, 203, 112138.1–112138.20. [Google Scholar] [CrossRef]
  2. Li, S.; Gong, W.; Yan, X. Parameter extraction of photovoltaic models using an improved teaching-learning-based optimization. Energy Convers. Manag. 2019, 186, 293–305. [Google Scholar] [CrossRef]
  3. 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] [CrossRef]
  4. Abd, M.; Oliva, D. Parameter estimation of solar cells diode models by an improved opposition-based whale optimization algorithm. Energy Convers. Manag. 2018, 171, 1843–1859. [Google Scholar]
  5. Yu, K.; Liang, J.; Qu, Y. Multiple learning backtracking search algorithm for estimating parameters of photovoltaic models. Appl. Energy 2018, 226, 408–422. [Google Scholar] [CrossRef]
  6. Abdel, M.; Mohamed, R.; Mirjalili, S. Solar photovoltaic parameter estimation using an improved equilibrium optimizer. Sol. Energy 2020, 209, 694–708. [Google Scholar] [CrossRef]
  7. Xiong, G.; Zhang, J.; Shi, D. Parameter extraction of solar photovoltaic models using an improved whale optimization algorithm. Energy Convers. Manag. 2018, 174, 388–405. [Google Scholar] [CrossRef]
  8. Yu, K.; Qu, B.; Yue, C. A performance-guided JAYA algorithm for parameters identification of photovoltaic cell and module. Appl. Energy 2019, 237, 241–257. [Google Scholar] [CrossRef]
  9. Premkumar, M.; Jangir, P.; Sowmya, R.; Elavarasan, R.M.; Kumar, B.S. Enhanced chaotic JAYA algorithm for parameter estimation of photovoltaic cell/module. ISA Trans. 2021, 116, 139–166. [Google Scholar] [CrossRef]
  10. Zhang, Y.; Ma, M.; Jin, Z. Comprehensive learning Jaya algorithm for parameter extraction of photovoltaic models. Energy 2020, 211, 118644. [Google Scholar] [CrossRef]
  11. Xiong, G.; Zhang, J.; Yuan, X. Parameter extraction of solar photovoltaic models by means of a hybrid differential evolution with whale optimization algorithm. Sol. Energy 2018, 176, 742–761. [Google Scholar] [CrossRef]
  12. Wei, D.; Wei, M.; Cai, H. Parameters extraction method of PV model based on key points of I-V curve. Energy Convers. Manag. 2020, 180, 112656.1–112656.8. [Google Scholar] [CrossRef]
  13. Ibrahim, H.; Anani, N. Evaluation of Analytical Methods for Parameter Extraction of PV modules. Energy Procedia 2017, 134, 69–78. [Google Scholar] [CrossRef]
  14. Abdel, M.; El-Shahat, D.; Chakrabortty, K. Parameter estimation of photovoltaic models using an improved marine predators algorithm. Energy Convers. Manag. 2021, 227, 113491. [Google Scholar] [CrossRef]
  15. Song, Y.; Wu, D.; Deng, W. MPPCEDE: Multi-population parallel co-evolutionary differential evolution for parameter optimization. Energy Convers. Manag. 2021, 228, 113661. [Google Scholar] [CrossRef]
  16. Abdel-Basset, M.; Mohamed, R.; Chakrabortty, R.K.; Sallam, K.; Ryan, M.J. An efficient teaching-learning-based optimization algorithm for parameters identification of photovoltaic models: Analysis and validations. Energy Convers. Manag. 2021, 227, 113614. [Google Scholar] [CrossRef]
  17. Pardhu, B.G.; Kota, V.R. Radial movement optimization based parameter extraction of double diode model of solar photovoltaic cell. Sol. Energy 2021, 213, 312–327. [Google Scholar] [CrossRef]
  18. Nunes, H.G.; Silva, P.N.; Pombo, J.A.; Mariano, S.J.; Calado, M.R. Multiswarm spiral leader particle swarm optimization algorithm for PV parameter identification. Energy Convers. Manag. 2020, 225, 113388. [Google Scholar] [CrossRef]
  19. Yu, K.; Chen, X.; Wang, X. Parameters identification of photovoltaic models using self-adaptive teaching-learning-based optimization. Energy Convers. Manag. 2017, 145, 233–246. [Google Scholar] [CrossRef]
  20. Tao, Y.; Bai, J.; Pachauri, K. Parameter extraction of photovoltaic modules using a heuristic iterative algorithm. Energy Convers. Manag. 2020, 224, 113386. [Google Scholar] [CrossRef]
  21. Gao, X.; Cui, Y.; Hu, J.; Xu, G.; Yu, Y. Lambert W-function based exact representation for double diode model of solar cells: Comparison on fitness and parameter extraction. Energy Convers. Manag. 2016, 127, 443–460. [Google Scholar] [CrossRef]
  22. Li, S.; Gu, Q.; Gong, W. An enhanced adaptive differential evolution algorithm for parameter extraction of photovoltaic models. Energy Convers. Manag. 2020, 205, 112443. [Google Scholar] [CrossRef]
  23. Li, S.; Gong, W.; Yan, X. Parameter estimation of photovoltaic models with memetic adaptive differential evolution. Sol. Energy 2019, 190, 465–474. [Google Scholar] [CrossRef]
  24. Lekouaghet, B.; Boukabou, A.; Chabane, B. Estimation of the photovoltaic cells/Modules parameters using an improved Rao-based chaotic optimization technique. Energy Convers. Manag. 2021, 229, 113722. [Google Scholar] [CrossRef]
  25. Jian, X.; Zhu, Y. Parameters identification of photovoltaic models using modified Rao-1 optimization algorithm. Optik 2021, 231, 166439. [Google Scholar] [CrossRef]
  26. Sallam, M.; Hossain, A.; Chakrabortty, K. An Improved Gaining-Sharing Knowledge Algorithm for Parameter Extraction of Photovoltaic Models. Energy Convers. Manag. 2021, 237, 114030. [Google Scholar] [CrossRef]
  27. Khan, A.; Ahmad, S.; Sarwar, A. Chaos Induced Coyote Algorithm (CICA) for Extracting the Parameters in a Single, Double, and Three Diode Model of a Mono-Crystalline, Polycrystalline, and a Thin-Film Solar PV Cell. Electronics 2021, 10, 2094. [Google Scholar] [CrossRef]
  28. Yu, K.; Liang, J.; Qu, Y. Parameters identification of photovoltaic models using an improved JAYA optimization algorithm. Energy Convers. Manag. 2017, 150, 742–753. [Google Scholar] [CrossRef]
  29. Jian, X.; Weng, Z. A logistic chaotic JAYA algorithm for parameters identification of photovoltaic cell and module models. Optik 2019, 203, 164041. [Google Scholar] [CrossRef]
  30. Wang, L.; Huang, C. A novel Elite Opposition-based Jaya algorithm for parameter estimation of photovoltaic cell models. Optik 2020, 203, 164034. [Google Scholar] [CrossRef]
  31. Venkata, R. Jaya: A simple and new optimization algorithm for solving constrained and unconstrained optimization problems. Int. J. Ind. Eng. Comput. 2016, 7, 19–34. [Google Scholar] [CrossRef]
Figure 1. The equivalent circuit of SDM.
Figure 1. The equivalent circuit of SDM.
Photonics 09 00131 g001
Figure 2. The equivalent circuit of DDM.
Figure 2. The equivalent circuit of DDM.
Photonics 09 00131 g002
Figure 3. The equivalent circuit of a PV module.
Figure 3. The equivalent circuit of a PV module.
Photonics 09 00131 g003
Figure 4. The flowchart of CSOOJAYA.
Figure 4. The flowchart of CSOOJAYA.
Photonics 09 00131 g004
Figure 5. Experimental and calculated data of CSOOJAYA for SDM.
Figure 5. Experimental and calculated data of CSOOJAYA for SDM.
Photonics 09 00131 g005
Figure 6. Experimental and calculated data of CSOOJAYA for DDM.
Figure 6. Experimental and calculated data of CSOOJAYA for DDM.
Photonics 09 00131 g006
Figure 7. Experimental and calculated data of CSOOJAYA for PV module.
Figure 7. Experimental and calculated data of CSOOJAYA for PV module.
Photonics 09 00131 g007
Figure 8. Convergence curves of the seven algorithms on SDM.
Figure 8. Convergence curves of the seven algorithms on SDM.
Photonics 09 00131 g008
Figure 9. Convergence curves of the seven algorithms on DDM.
Figure 9. Convergence curves of the seven algorithms on DDM.
Photonics 09 00131 g009
Figure 10. Convergence curves of the seven algorithms on PV.
Figure 10. Convergence curves of the seven algorithms on PV.
Photonics 09 00131 g010
Table 1. The characteristics of the three methods.
Table 1. The characteristics of the three methods.
MethodCharacteristicInstance
Analytical method
-
Advantages such as simple, fast and direct calculation.
-
If the initial value is not selected well, it is necessary to repeatedly assign a new initial value for a new round of solution, and the error in the solution process will also increase with the increase of the number of identification parameters, which requires a large amount of calculation and a long time.
Deterministic method
-
Convert the PV parameter identification problem into an optimization problem
-
The objective function has characteristics such as continuity, convexity and differentiability. Easily affected by initial conditions and gradient information.
-
It is easy to get stuck in local optima when dealing with complex multimodal problems.
Newton–Raphson method [16,17], iterative curve fitting [18,19] and Lambert W-function [20,21]
Heuristic method
-
Overcome the shortcomings of the above two methods, the solution speed is significantly improved, and the accuracy is improved
-
contains some control parameters
CPMPSO [1], MTLBO [16], MSLPSO [18], EJADE [22], MADE [23], MRAO-1 [25]
Table 2. Summary of some existing heuristic methods.
Table 2. Summary of some existing heuristic methods.
AlgorithmSummary
CPMPSO [1]
-
Introducing the classification perturbation mutation strategy to improve the local search ability and population diversity of the algorithm.
PGJAYA [8]
-
The performance of each individual is quantified by probability.
An adaptive chaotic perturbation mechanism is introduced to improve the quality of the population.
CLJAYA [10]
-
Introducing a comprehensive learning mechanism to improve the global search ability of the algorithm.
MTLBO [16]According to the score level of each learner, each teaching stage is divided into three levels to update the population.
MSLPSO [18]
-
Combine multiple populations with different search mechanisms to maintain a balance between exploration and exploitation.
-
Guided by several different spiral trajectories, the algorithm maintains the diversity of the population while exploring different regions of the multidimensional space.
EJADE [22]
-
Introduce crossover rate sorting mechanism to retain excellent individuals.
-
Introduce a dynamic population reduction strategy to improve convergence speed and balance exploration and exploitation.
MADE [23]
-
Adaptive differential evolution strategy based on success-history improves global search ability.
-
A ranking-based elimination strategy is proposed to retain promising individuals in the population.
MRAO-1 [25]
-
A bidirectional update strategy is proposed to improve population diversity.
IGSK [26]
-
Introduce an adaptive mechanism to automatically adjust the knowledge rate parameter.
-
A constraint handling method is proposed to improve the convergence speed and maintain the balance between exploration and exploitation.
-
Automatically select a binding processing strategy from three algorithm pools to improve the performance of the algorithm.
CICA [27]
-
Chaos map can solve problems such as easy to fall into local optimum and slow algorithm convergence speed.
IJAYA [28]
-
Experience-based learning strategies maintain the diversity of the population and improve the exploration ability of the population.
-
A chaotic elite learning method is proposed to improve the quality of the best solution in each generation.
LCJAYA [29]
-
Introduces a chaotic mutation strategy to improve and balance exploration and exploitation.
EO-JAYA [30]
-
Introducing an elite opposition learning strategy to increase population diversity.
Table 3. Parameters ranges of PV models.
Table 3. Parameters ranges of PV models.
ParameterSDM/DDMPV Module
LowerUpperLowerUpper
0102
01050
00.502
010002000
12150
Table 4. Parameter settings of the compared algorithms.
Table 4. Parameter settings of the compared algorithms.
AlgorithmParameters
DE/WOANp = 40, F = rand (0.1,1), Cr = rand (0,1)
MLBSANp = 50
EJADE Np max = 50, Np min = 4
MTLBONp = 50
MPPCEDENp = 40; F = rand (0.1,1); Cr = rand (0,1)
MADENP = 20; 0.05 for the SDM and PV module; 0.01 for the DDM
IJAYANP = 20
ITLBONp = 50
CPMPSONP = 50, w = 0.729, c1 = 1.49455, c2 = 1.49455
DERAONP = 20
JAYANP = 20
PGJAYANP = 20
LCJAYANP = 20
CSOOJAYANP = 20
Table 5. Comparison of the extracted optimal parameters of the SDM.
Table 5. Comparison of the extracted optimal parameters of the SDM.
Algorithm I p h (A)   I d ( μ A) R S ( Ω ) R s h ( Ω )nRMSE
DE/WOA0.7607760.3230210.03637753.7185241.4811849.860219 × 104
MLBSA0.76080.323020.036453.71851.48129.8602 × 10−4
EJADE0.76080.32300.036453.71851.48129 8602 × 10−4
MTLBO0.760775530.323000.0363770953.71852511.481183599 860219 × 10−4
MPPCEDE0.76080.323020.036453.71851.48129.8602 × 10−4
MADE0.76080.32300.036453.71851.48129.8602 × 10−4
IJAYA0.76080.32280.036453.75951.48119.8603 × 10−4
ITLBO0.76080.32300.036453.71851.48129.8602 × 10−4
CPMPSO0.7607760.3230210.03637753.7185221.4811849.860219 × 10−4
DERAO0.7607760.3230210.03637753.7185221.4811359.860219 × 10−4
JAYA0.76080.32810.036454.92981.48289.8946 × 10−4
PGJAYA0.76080.32300.036453.71851.48129.8602 × 10−4
LCJAYA0.76080.32300.036453.71851.48199.8602 × 10−4
CSOOJAYA0.7607760.3230210.03637753.7185251.4811849.860219 × 10−4
Table 6. IAE of CSOOJAYA for the SDM.
Table 6. IAE of CSOOJAYA for the SDM.
Item V m e a s u r e d / V I m e a s u r e d / A I e s t i m a t e d / A I A E c u r r e n t   P e s t i m a t e d / W I A E p o w e r
1−0.20570.76400.76408817470.0000881747−0.15717293750.0000181375
2−0.12910.76200.76266355700.0006635570−0.09845986520.0000856652
3−0.05880.76050.76135577800.0008557780−0.04476771970.0000503197
40.00570.76050.76015446180.00034553810.00433288040.0000019695
50.06460.76000.75905567970.00094432020.04903499690.0000610030
60.11850.75900.75804281610.00095718380.08982807370.0001134262
70.16780.75700.75709212480.00009212480.12704005850.0000154585
80.21320.75700.75614183580.00085816410.16120943930.0001829606
90.25450.75550.75508734390.00041265600.19216972900.0001050209
100.29240.75400.75366434990.00033565000.22037145590.0000981440
110.32690.75050.75139143920.00089143920.24562986150.0002914115
120.35850.74650.74735432600.00085432600.26792652580.0003062758
130.38730.73850.74011769970.00161769970.28664758510.0006265351
140.41370.72800.72738270750.00061729240.30091822610.0002553738
150.43730.70650.70697313900.00047313900.30915935370.0002069037
160.45900.67550.67528064250.00021935740.30995381490.0001006850
170.47840.63200.63075875910.00124124080.30175499030.0005938096
180.49600.57300.57192882310.00107117680.28367669630.0005313036
190.51190.49900.49960743170.00060743170.25574904420.0003109442
200.52650.41300.41364911060.00064911060.21778625670.0003417567
210.53980.31650.31751027880.00101027880.17139204850.0005453485
220.55210.21200.21215489450.00015489450.11713071720.0000855172
230.56330.10350.10225098790.00124901200.05759798150.0007035684
240.5736−0.0100−0.00871821290.0012817870−0.00500076690.0007352330
250.5833−0.1230−0.12550850170.0025085017−0.07320910900.0014632090
260.5900−0.2100−0.20847376600.0015262339−0.12299952190.0009004780
Table 7. Comparison of the extracted optimal parameters of the DDM.
Table 7. Comparison of the extracted optimal parameters of the DDM.
Algorithm I p h (A) I d 1 ( μ A) I d 2 ( μ A) R S ( Ω ) R s h ( Ω )n1n2RMSE
DE/WOA0.7607810.2259740.7493460.03674055.4854371.45101729.824849 × 10−4
MLBSA0.76080.227280.738350.036755.46121.451529.8249 × 10−4
EJADE0.76080.22600.74930.036755.48541.451029.8248 × 10−4
MTLBO0.76078100.225970.74930.0367404355.4854471.45101629.824849 × 10−4
MPPCEDE0.76080.227280.738350.036755.46121.451529.82485 × 10−4
MADE0.76080.22460.73940.036855.43291.45051.99639.8261 × 10−4
IJAYA0.76010.00504450.750940.037677.85191.21861.62479.8293 × 10−4
ITLBO0.76080.22600.74930.036755.48541.451029.8248 × 10−4
CPMPSO0.760780.225970.749350.0367455.485441.4510229.824849 × 10−4
DERAO0.7607810.2259590.7496670.03674055.4860451.45096329.824849 × 10−4
JAYA0.76070.00607630.315070.036452.65751.47881.84369.8934 × 10−4
PGJAYA0.76080.210310.885340.036855.81351.445029.8263 × 10−4
LCJAYA0.76080.225960.746400.036755.48151.45182.00009.8250 × 10−4
CSOOJAYA0.7607810.2259740.7493450.03674055.4854371.4510172.0000009.824849 × 10−4
Table 8. IAE of CSOOJAYA for the DDM.
Table 8. IAE of CSOOJAYA for the DDM.
Item V m e a s u r e d / V I m e a s u r e d / A I e s t i m a t e d / A I A E c u r r e n t   P e s t i m a t e d / W I A E p o w e r
1−0.20570.76400.76398341180.0000165881−0.15715138780.0000034121
2−0.12910.76200.76260409570.0006040957−0.09845218870.0000779887
3−0.05880.76050.76133769800.0008376980−0.04476665660.0000492566
40.00570.76050.76017378770.00032621220.00433299050.0000018594
50.06460.76000.75910768010.00089231980.04903835610.0000576438
60.11850.75900.75812141980.00087858010.08983738820.0001041117
70.16780.75700.75718861360.00018861360.12705624930.0000316493
80.21320.75700.75624360680.00075639310.16123113690.0001612630
90.25450.75550.75517730150.00032269840.19219262320.0000821267
100.29240.75400.75372235350.00027764640.22038841610.0000811838
110.32690.75050.75139913420.00089913420.24563237690.0002939269
120.35850.74650.74730144320.00080144320.26790756740.0002873174
130.38730.73850.74001066070.00151066070.28660612880.0005850788
140.41370.72800.72724695250.00075304740.30086206420.0003115357
150.43730.70650.70685029770.00035029770.30910563520.0001531852
160.45900.67550.67521054250.00028945740.30992163900.0001328609
170.47840.63200.63076075760.00123924230.30175594640.0005928535
180.49600.57300.57199473290.00100526700.28370938750.0004986124
190.51190.49900.49970613520.00070613520.25579957060.0003614706
200.52650.41300.41373367280.00073367280.21783077870.0003862787
210.53980.31650.31754620570.00104620570.17141144180.0005647418
220.55210.21200.21212299600.00012299600.11711310610.0000679061
230.56330.10350.10216327650.00133672340.05754857360.0007529763
240.5736−0.0100−0.00879175090.0012082490−0.00504294830.0006930516
250.5833−0.1230−0.12554343500.0025434350−0.07322948560.0014835856
260.5900−0.2100−0.20837158910.0016284108−0.12293923760.0009607623
Table 9. Comparison of the extracted optimal parameters of the PV module.
Table 9. Comparison of the extracted optimal parameters of the PV module.
Algorithm I p h (A) I d ( μ A) R S ( Ω ) R s h ( Ω )nRMSE
DE/WOA1.0305143.4822631.201271981.98214348.6428352.425075 × 10−3
MLBSA1.03053.48231.2013981.982348.64282.4251 × 10−3
EJADE1.03053.48231.2013981.982448.64282.4251 × 10−3
MTLBO1.03051433.48231.2012710981.982373248.64283492425075 × 10−3
MPPCEDE1.0305143.4822631.201271981.98214348.6428352.42507 × 10−3
MADE1.03053.48231.2013981.982348.64282.4250 × 10−3
IJAYA1.03053.47031.2016977.375248.62982.4251 × 10−3
ITLBO1.03053.48231.2013981.982348.64282.4251 × 10−3
CPMPSO1.0305143.48231.201271981.982348.642842.425075 × 10−3
DERAO1.0305143.48231.201271981.982148.641312.425075 × 10−3
JAYA1.03023.49311.20141022.548.65312.427785 × 10−3
PGJAYA1.03053.48181.2013981.854548.64242.425075 × 10−3
LCJAYA1.03053.48231.2024981.982848.66842.425075 × 10−3
CSOOJAYA1.0305143.4822631.201271981.98224348.6428342.425075 × 10−3
Table 10. IAE of CSOOJAYA for PV module.
Table 10. IAE of CSOOJAYA for PV module.
Item V m e a s u r e d / V I m e a s u r e d / A I e s t i m a t e d / A I A E c u r r e n t   P e s t i m a t e d / W I A E p o w e r
10.12481.03151.02911886220.00238113770.12843403400.0002971659
21.80931.03001.02738077400.00261922591.85884003450.0047389654
33.35111.02601.02574149780.00025850213.43736233330.0008662666
44.76221.02201.02410685560.00210685564.87700166800.0100332680
56.05381.01801.02229150540.00429150546.18874831540.0259799154
67.23641.01551.01993038160.00443038167.38062421380.0320600138
78.31891.01401.01636280630.00236280638.45502054960.0196559496
89.30971.01001.01049585170.00049585179.40741323120.0046162312
910.21631.00351.00062866980.002871330110.22272267940.0293343705
1011.04490.98800.98454807790.003451922010.87423506640.0381261335
1111.80180.96300.95952137460.003478625311.32407935920.0410540407
1212.49290.92550.92283851520.002661484711.52892928660.0332496633
1313.12310.87250.87259935810.000099358111.45120863650.0013038865
1413.69830.80750.80727395660.000226043311.05828083960.0030964103
1514.22210.72650.72833616800.001836168010.35846981610.0261141661
1614.69950.63450.63713768680.00263768689.36560542790.0387726779
1715.13460.53450.53621274640.00171274648.11536543200.0259217320
1815.53110.42750.42951100440.00201100446.67077836100.0312331110
1915.89290.31850.31877415840.00027415845.06624582270.0043571727
2016.22290.20850.20738917840.00111082153.36445390350.0180207464
2116.52410.10100.09616683970.00483316021.58907047680.0798636231
2216.7987−0.0080−0.00832572170.0003257217−0.13986130110.0054717011
2317.0499−0.1110−0.11093682170.0000631782−1.89146171730.0010771826
2417.2793−0.2090−0.20924760820.0002476082−3.61565219700.0042784970
2517.4885−0.3030−0.30086393220.0021360677−5.26165887990.0373566200
Table 11. Comparison of statistical results for PV models.
Table 11. Comparison of statistical results for PV models.
ModelAlgorithmRMSE
MinMaxMeanSD
SDMDE/WOA9.860219 × 10−49.860219 × 10−49.860219 × 10−43.545178 × 10−17
MLBSA9.8602 × 10−49.8602 × 10−49.8602 × 10−49.1461 × 10−12
EJADE9.8602 × 10−49.8602 × 10−49.8602 × 10−45.13 × 10−17
MTLBO9.860219 × 10−49.860219 × 10−49.860219 × 10−41.9092749 × 10−17
MPPCEDE9.86022 × 10−49.86022 × 10−49.86022 × 10−42.44932 × 10−17
MADE9.8602 × 10−49.8602 × 10−49.8602 × 10−42.74 × 10−15
IJAYA9.8603 × 10−41.0622 × 10−39.9204 × 10−41.4033 × 10−5
ITLBO9.8602 × 10−49.8602 × 10−49.8602 × 10−42.19 × 10−17
CPMPSO9.86022 × 10−49.86022 × 10−49.86022 × 10−42.17556 × 10−17
DERAO9.860219 × 10−49.860219 × 10−49.860219 × 10−43.642031 × 10−17
JAYA9.8946 × 10−41.4783 × 10−31.1617 × 10−31.8796 × 10−4
PGJAYA9.8602 × 10−49.8603 × 10−49.8602 × 10−41.4485 × 10−9
LCJAYA9.8602 × 10−49.8602 × 10−49.8602 × 10−45.6997 × 10−16
CSOOJAYA9.860219 × 10−49.860219 × 10−49.860219 × 10−44.717305 × 10−17
DDMDE/WOA9.824849 × 10−49.860377 × 10−49.829703 × 10−49.152178 × 10−7
MLBSA9.8249 × 10−49.8518 × 10−49.8798 × 10−41.3482 × 10−6
EJADE9.8248 × 10−49.8602 × 10−49.8363 × 10−41.36 × 10−6
MTLBO9.824849 × 10−49.825026 × 10−49.824855 × 10−43.30000 × 10−9
MPPCEDE9.82485 × 10−49.82908 × 10−49.82504 × 10−48.02951 × 10−8
MADE9.8261 × 10−49.8786 × 10−49.8608 × 10−48.02 × 10−5
IJAYA9.8293 × 10−41.4055 × 10−31.0269 × 10−39.8325 × 10−5
ITLBO9.8248 × 10−49.8812 × 10−49.8497 × 10−41.54 × 10−6
CPMPSO9.82485 × 10−49.83137 × 10−49.86022 × 10−41.3398 × 10−6
DERAO9.824841 × 10−49.824897 × 10−49.824845 × 10−41.01560 × 10−9
JAYA9.8934 × 10−41.4793 × 10−31.1767 × 10−31.9356 × 10−4
PGJAYA9.8263 × 10−49.9499 × 10−49.8582 × 10−42.5375 × 10−6
LCJAYA9.8250 × 10−49.8602 × 10−49.8308 × 10−41.3118 × 10−6
CSOOJAYA9.824849 × 10−49.824849 × 10−49.824849 × 10−45.576332 × 10−17
PV module DE/WOA2.425075 × 10−32.425442 × 10−32.425092 × 10−36.270718 × 10−8
MLBSA2.425075 × 10−32.425084 × 10−32.425312 × 10−34.336794 × 10−8
EJADE2.4251 × 10−32.4251 × 10−32.4251 × 10−33.27 × 10−17
MTLBO2.425075 × 10−32.425075 × 10−32.425075 × 10−31.3070107 × 10−17
MPPCEDE2.42507 × 10−32.42507 × 10−32.42507 × 10−32.23156 × 10−17
MADE2.4250 × 10−32.4251 × 10−32.4251 × 10−33.07 × 10−17
IJAYA2.4251 × 10−32.4393 × 10−32.4289 × 10−33.7755 × 10−6
ITLBO2.4251 × 10−32.4251 × 10−32.4251 × 10−31.27 × 10−17
CPMPSO2.425075 × 10−32.425075 × 10−32.425075 × 10−31.515959 × 10−17
DERAO2.425075 × 10−32.425075 × 10−32.425075 × 10−31.716113 × 10−17
JAYA2.427785 × 10−32.595873 × 10−32.453710 × 10−33.456290 × 10−5
PGJAYA2.425075 × 10−32.426764 × 10−32.425144 × 10−33.071420 × 10−7
LCJAYA2.425075 × 10−32.425075 × 10−32.425075 × 10−32.415229 × 10−16
CSOOJAYA2.425075 × 10−32.425075 × 10−32.425075 × 10−32.699858 × 10−17
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jian, X.; Cao, Y. A Chaotic Second Order Oscillation JAYA Algorithm for Parameter Extraction of Photovoltaic Models. Photonics 2022, 9, 131. https://doi.org/10.3390/photonics9030131

AMA Style

Jian X, Cao Y. A Chaotic Second Order Oscillation JAYA Algorithm for Parameter Extraction of Photovoltaic Models. Photonics. 2022; 9(3):131. https://doi.org/10.3390/photonics9030131

Chicago/Turabian Style

Jian, Xianzhong, and Yougang Cao. 2022. "A Chaotic Second Order Oscillation JAYA Algorithm for Parameter Extraction of Photovoltaic Models" Photonics 9, no. 3: 131. https://doi.org/10.3390/photonics9030131

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