An Effective Method for Parameter Estimation of a Solar Cell

Parameter extraction of the photovoltaic cell is a highly nonlinear complex optimization problem. This article proposes a new hybrid version of whale optimization and particle swarm optimization algorithm to optimize the photovoltaic cell parameters. The exploitation ability of particle swarm optimization with adaptive weight function is implemented in the pipeline mode with a whale optimization algorithm to improve its exploitation capability and convergence speed. The performance of the proposed hybrid algorithm is compared with six different optimization algorithms in terms of root mean square error and rate of convergence. The simulation result shows that the proposed hybrid algorithm produces not only optimized parameters at different irradiation levels (i.e., 1000 W/m2, 870 W/m2, 720 W/m2, and 630 W/m2) but also estimates minimum root mean square error even at a low level of irradiations. Furthermore, the statistical analysis validates that the average accuracy and robustness of the proposed algorithm are better than other algorithms. The best values of root mean square error generated by the proposed algorithm are 7.1700×10−4 and 9.8412×10−4 for single-diode and double-diode models. It is observed that the estimated parameters based on the optimization process are highly consistent with the experimental data.


Introduction
The depletion of fossil fuel resources and resulting environmental impact due to their usages embarks the need for alternate energy resources [1]. Solar energy is one of the most promising alternative sources for fossil fuel. The free access to the energy of sunlight can be extracted employing the photovoltaic (PV) panels. The rapid adoption of solar energy by the domestic and industrial sector makes it a vital source to be explored [2]. Despite the very low operational and maintenance cost, there are various limitations for efficient energy generation. An enormous amount of research has been performed and carried out to better the power output from the PV panels [3,4]. The major limitation in the execution and implementation of the solar PV power plants is the very high capital cost for installation [5]. PV cells are having nonlinear current-voltage (I-V) and power-voltage (P-V) characteristics curves with some operational limitations [6]. This non-linearity makes it difficult for any probability and approximation to increase efficiency. Every PV panel can operate at maximum efficiency, as defined by the manufacturer, only if the practical parameters (voltage-current) are somewhat close to or coinciding with the maximum power point (MPP). The real behaviour of PV panels rather different from the optimal conditions, due to the non-linearity of I-V characteristics of solar cells makes it essential to determine the MPP in each moment. It could be done through simulation techniques for better operational efficiency [7]. This technology is ensured by the model of the equivalent extraction of solar PV cell. In this approach, the hybridization of two algorithms are implemented, the firefly and pattern search. The exploration phase is completed by the firefly algorithm during the first half iteration and then the pattern search algorithm takes control of the population for the exploitation phase. A new opposition-based learning approach is incorporated with whale optimization and shuffled complex evolutionary algorithm for optimization of solar cell parameters [35,36]. This approach is tested on unimodal as well as on multimodal benchmark functions and simulation results clearly show the robustness of the algorithms.
The whale optimization algorithm (WOA) [37] and particle swarm optimization (PSO) [38] are the two most prominent used metaheuristics techniques as available in the literature. However, they differ from each other in the search mechanism for the best solution in a defined search space. WOA mimics the social behaviour of humpback whales while PSO mimics the searching behaviour of the birds in a group. It is shown by many previous research studies that WOA is good at exploring [39] the search space but suffers from a slow convergence rate due to low exploitation ability while PSO don't have good capability in exploring [40] the search space but have good local search capability. the convergence speed of the algorithm. In [41], the author proposed a chaotic WOA (CWOA) to improve maps utilized their dynamic behavior to prevent an optimization algorithm to trap in local optima and improves its global search capability. In [42], the author proposed Levy flight trajectory based WOA (LWOA) to improve the accuracy and convergence speed of the algorithm. Levy flight allowed for the algorithm to get rid of local optima and prevents premature convergence.
There are certain complex and non-convex optimization problems that are not solved by continuous metaheuristic therefore, in [43], the author proposed binary WOA (BWOA). In [44], the author proposed a modified WOA that includes whale memory and new random search agent to enhance the exploitation capability of the algorithm. In [45], the author improved the exploration capability of WOA and proposed three modified WOA which are based on opposition-based learning, exponentially decreasing parameters, and re-initialization of the worst particles. The hybridization of metaheuristic algorithms is another approach to improve the exploration and exploitation capability of population based stochastic algorithm. Furthermore, researchers have proposed hybrid approach grey wolf optimization (HAGWO) [46], WOA-CBO (colliding bodies optimization) [47], memetic-WOA (MWOA) [48], WOA-SA (simulated annealing) [39], WOA-MFO (moth flame optimization) [49], Sine-Cosine (SC-WOA) [50], WOA-PS (Pattern Search) [51], and Brain Storm (BS-WOA) [52][53][54] to improve the global and local search capability of WOA.
According to the literature survey, WOAPSO has not yet been implemented for the parameter extraction of the solar cell (and it cannot be used to establish a PV parameter estimation technique that can overcome all existing techniques). Therefore, this research paper aims to anticipate a new parameter estimation algorithm for solar cell/module. The novelty of the proposed study is that the exploitation capability of WOA is significantly improved by incorporating the exploitation capability of PSO with adaptive weight in sequential mode. As a result, equivalent circuit parameters converge equally good to the true values with minimum error. The proposed WOAPSO algorithm's performance is measured based on convergence analysis, robustness, reliability, and statistical analysis for three PV models at diverse operating conditions. The manuscript is organized as follows: the problem formulation and mathematical model for solar PV cell/module are presented in Section 2. Section 3 gives a brief introduction of the WOA, PSO, and proposed WOAPSO algorithm and discussed its implementation to estimate the optimized value of unknown parameters of a single-diode, double-diode, and PV module model. In Section 4, the simulation results of the WOAPSO algorithm are discussed and compared with pre-existing metaheuristic algorithms. Finally, Section 5 provides a conclusive remark to summarize the paper.

Methodology
In this section, the equivalent circuits of a photovoltaic solar cell are formulated using a single-diode and double-diode models. These equivalent circuit models are used to describe the current-voltage characteristics of a solar cell.

PV Panel Model
The equivalent circuit of PV panel module is shown in Figure 1. The relation between current and voltage at output terminal for the PV panel module is expressed as: where N s and N p represents the number of solar cells connected in series and parallel respectively. It is clearly depicted from Figure 1 that only five parameters (I p , I SD , a 1 , R s and R sh ) are needed to be estimated for minimum value of the RMSE.

Objective Function
The key purpose of this work is to optimize the unknown parameters for both the models (SDM and DDM) and to reduce the error between experimental and estimated data. The objective function for error used here is same as the authors have used previously in as: where V l and I l are the measured voltage and current of PV module. The parameter 'k' stands for the number of experimental data set. The best solution found by WOAPSO is represented by a vector X. For the single-diode model: For the double-diode model: − I l X = I p , I SD1 , I SD2 , a 1 , a 2 , R s , R sh For the PV panel module model:

Hybrid Algorithm
The hybridization of the metaheuristic algorithm plays a vital role in improving their performance. The fundamental principle of hybridization is to blend the best features of two or more metaheuristic algorithms to improve search capability, accuracy, and convergence speed of an individual algorithm. A hybrid algorithm is also known as a memetic algorithm. In the last few years, researchers have proposed different strategies for hybridizing metaheuristic algorithms. The three most explored methodologies of hybridization are multi-stage, sequential and parallel.
In the multi-stage methodology, one of the algorithms globally explores the search space and the second algorithm locally discovers the optimal solution. In sequential search, both the algorithms run sequentially and find the optimal solution in the search space. In the parallel mode, both the algorithms run parallel on the same population of the defined problem.

Particle Swarm Optimization (PSO)
Particle swarm optimization is a nature inspired stochastic optimization technique proposed by J. Kennedy and R. C. Eberhard in 1995. It is a population-based computationally inexpensive technique that is inspired by the social behaviour of fish schooling and bird flocking. The methodology of the algorithm is that the swarm of particles fly in the search space and finds the optimal solution by updating their own best solution and the best solution obtained by the swarms. The swarm is randomly initialized as particles in N-dimensional search space with position x i and velocity v i . The position of the particles represents the probable solution, and the velocity represents the rate of change of position of the particle concerning the current position. The particles change their positions with respect to the positions of the best particle. The velocity update equations are given by: where v d i (t) and x d i (t) represents the velocity and position of ith particle in dth dimension at tth iteration, v d i (t + 1) and x d i (t + 1) is the velocity and position of the ith particle in dth dimension at (t + 1)th iteration. pbest d i represents the current best position of the particles and gbest d represents the best position among all the particles in dth dimension, c 1 and c 2 are the acceleration parameter, r 1 and r 2 are the random number in the range [0, 1] and w is the inertial weight vector which maintains balance between exploration and exploitation.

Whale Optimization Algorithm (WOA)
The whale optimization algorithm is a population-based optimization algorithm that mimics the social behaviour of humpback whales and was proposed by Mirjalili and Lewis in 2016. Humpback whales are long in size and have an interesting food searching capability: they attack their prey (krill and small fishes) by a bubble-net hunting strategy. WOA is inspired by this hunting behaviour, and works in three phases. First, it searches for prey then encircles the prey and lastly, attacks the prey. Humpback whales swim around the prey either following a shrinking path or through a spiral movement. A probability factor p assumed to be 50% simultaneously choose either of the two movements.

Shrinking Movement
Initially in the exploration phase, humpback whales search around a prey chosen randomly in the search space with the following mathematical model: where t is the current iteration and (t + 1)th is the next iteration, → X rand is the random position of the prey, → A and → C are the coefficient vectors and is defined as: where → a is decreased from 2 to 0 over the course of iterations and → r is the random number in the range [0, 1]. In the exploitation phase the position of whales are updated based on the position of the best search prey → X * . Mathematically it is defined as:

Spiral Movement
In the spiral movement of the humpback whale, first the distance is evaluated between the whale located at (X, Y) and best search prey located at (X*, Y*). Once the distance is evaluated then the helix-shaped movement of whale around the prey is defined with following mathematical equation: where is the distance between the whale and best searched prey, b is the constant which maintains the shape of the logarithmic spiral and l is the random number defined in the range [-1, 1].
In WOA, coefficient vector 'A' maintains the balance in exploration and exploitation, when the value of p < 0.5 and A > 1 then the positions are updated by Equations (9) and (13) while when p < 0.5 and A < 1 the positions are updated by Equations (13) and (14) and when the p ≥ 0.5 then the positions are updated using Equation (14).

Hybrid WOAPSO Algorithm
In this section, the principle of the proposed hybrid WOAPSO algorithm is briefly addressed. In general, the performance of any optimization technique while solving any NLP problem is affected by premature convergence and slow rate of convergence. Some algorithms better explore the search space and have a slow convergence rate while some algorithms less diversely explore the search space and did not find the optimal solution. Maintaining the balance between exploration and exploitation is a critical issue in any optimization algorithm. WOA has good exploration capability but exploitation depends on evaluating the distance between the whale and the best position of the prey, and if the distance is large then it takes more time to converge. While PSO has fast rate of convergence but it is prone to premature convergence due to weakness in global search capability. Since in PSO, if the global best solution gets trapped in local optima, then the rest of the particles do not explore the search space and follow the global best solution, and become trapped in local optima. Therefore, it can be concluded that WOA is good at exploring the search space, but suffers from a slow convergence rate while PSO doesn't have good capability in exploring the search space but have good local search capability. The aim of the proposed hybrid algorithm is to enhance the exploitation capability of WOA by embedding the PSO algorithm to find an optimal solution around the region explored by WOA. The proposed approach is mixed, co-evolutionary in which PSO is used as a component of WOA and thus the hybrid approach utilizes the strength of both the algorithms to avoid the premature convergence and local optima. Figure 2 depicts the process flow chart of the proposed algorithm. The mathematical model of the proposed algorithm is illustrated in the following steps: Step 1: Initialize the random population of search agents with position and velocity defined as: Step 2: Calculate the fitness of each search agent. If the problem is the minimization problem, then → X * is the position corresponding to the minimum fitness and for maximization problem → X * is the position corresponding to the maximum fitness. → X * is the best search agent.
Step 3: Update the constant parameters A, C, using Equations (10) and (11) and l lying between [-1, 1] and p is the probability between 0 and 1.
Step 4: If p < 0.5 and |A|≥1, then select the random position of search agent (X*) in search space and update the position of search agent using Equations (9) and (13).
Else if p < 0.5 and |A|<1, then update the position of search agent using Equations (13) and (14).
Else p > 0.5, then update the position of search agent using Equation (14).
Step 5: Update the velocity of search agent based on the best position of search agent (X*) in the search space using the following equation: Step 6: Update the position of the particles using Equation (17).
Step 7: Go to step 3 until the termination criteria is met. The algorithm terminates when either maximum number of iterations or minimum error criteria is attained.
Step 8: In the last iteration the returned value of → X * represents the global minimum and the position corresponding to it represents the solution of the problem.

Implementation of WOAPSO for Parameter Extraction Single-Diode Model
Initialize the population of search agents of fifth order dimension in the search space. The fifth order dimension represents the photovoltaic current (I p ), series resistance (R s ), shunt resistance (R sh ), diode saturation current (I SD ) and diode ideality factor (a 1 ). The range of these parameters are [0-1, 0.001-0.5, 0-100, 0.01-0.5, 1-2].
Regulate the fitness of all agents in the search space using Equation (3). Update the position of the agents at every iteration using WOAPSO. The algorithm is designed to work in the minimization mode thus the location of particles that acquire minimum cost represents the optimized parameters of SDM with minimum RMSE.
Regulate the fitness of all agents in the search space using Equation (3). Update the position of the agents at every iteration using WOAPSO. The algorithm is designed to work in the minimization mode thus the location of particles that acquire minimum cost represents the optimized parameters of SDM with minimum RMSE.
Regulate the fitness of all agents in the search space using Equation (4).
Update the position of all agents at every iteration using WOAPSO. The algorithm is designed to work in the minimization mode. Thus, the location of particles having minimum cost represents the parameters of the double-diode model with minimum RMSE.

Results
In this section, the feasibility of the proposed new hybrid WOAPSO was tested and evaluated using mainly two types of PV devices: one PV cell (R.T.C France solar cell) and one PV module (SS2018P) at different solar irradiation. As a result, the retrieved PV cell and module parameters were monitored and used to create simulated I-V data for each device type. The accuracy and reliability of the WOAPSO were assessed by comparing the techniques published in the literature with the existing art. The efficiency of the proposed method is evaluated based on distinct empirical tools such as the individual absolute error Regulate the fitness of all agents in the search space using Equation (4). Update the position of all agents at every iteration using WOAPSO. The algorithm is designed to work in the minimization mode. Thus, the location of particles having minimum cost represents the parameters of the double-diode model with minimum RMSE.

Results
In this section, the feasibility of the proposed new hybrid WOAPSO was tested and evaluated using mainly two types of PV devices: one PV cell (R.T.C France solar cell) and one PV module (SS2018P) at different solar irradiation. As a result, the retrieved PV cell and module parameters were monitored and used to create simulated I-V data for each device type.  Table 1, which are the same as those being used by investigators in [27]. The proposed WOAPSO algorithm is implemented on MATLAB 2018a platform with Intel ® core ™ i7-HQ CPU, 2.4 GHz, 16 GB RAM Laptop. In order to conduct the experiment, the sample size, and the estimated number of objective function evaluations are set at 30 and 50,000, respectively. Furthermore, a minimum of 30 separate runs are carried out to prevent the contingency.

Parameter Estimation of Single-Diode Model Using WOAPSO
Only five parameters (I p , I sd , a, R s , R sh ) are required to be estimated for a single-diode model. Table 2 signifies the values of parameters optimized by WOAPSO and RMSE for the comparison. The WOAPSO algorithm provides the lowest RMSE of 7.1700 × 10 −4 than others (Table 2 and Table S2). Here RMSE values are acquired as the index for the evaluation of results with previously existing algorithms implemented by the researchers. The characteristics curve of current-voltage and power-voltage for a single-diode model is redrawn based on the best optimized parameters obtained by implementing the WOAPSO algorithm and depicted in Figure 3. It is observed that the calculated data obtained by the WOAPSO is very effectively in keeping with the experimental data set, under S.T.C (i.e., 1000 w/m 2 and 33 • C), all over the voltage range. The error relating the measurement results for each of 20 pair points is determined by IAE and RE, which is calculated by using Equations (18) and (19), respectively.

WOAPSO for Parameter Estimation of Double-Diode Model
In the case of DDM, the seven parameters (Ip, Isd, Isd1, a1, a2, Rs, Rsh) are required to be optimized. The values of optimized parameters and minimum of RMSE are presented in Table 3. The characteristics curve in terms of current-voltage and power-voltage for the double-diode model is redrawn based on the best optimized parameters (Figure 4). It can be observed that the estimated data based on optimized parameters are in keeping with the experimental data set.

WOAPSO for Parameter Estimation of Double-Diode Model
In the case of DDM, the seven parameters (I p , I sd , I sd1 , a 1 , a 2 , R s , R sh ) are required to be optimized. The values of optimized parameters and minimum of RMSE are presented in Table 3. The characteristics curve in terms of current-voltage and power-voltage for the double-diode model is redrawn based on the best optimized parameters (Figure 4). It can be observed that the estimated data based on optimized parameters are in keeping with the experimental data set.

WOAPSO for Parameter Estimation of Double-Diode Model
In the case of DDM, the seven parameters (Ip, Isd, Isd1, a1, a2, Rs, Rsh) are required to be optimized. The values of optimized parameters and minimum of RMSE are presented in Table 3. The characteristics curve in terms of current-voltage and power-voltage for the double-diode model is redrawn based on the best optimized parameters (Figure 4). It can be observed that the estimated data based on optimized parameters are in keeping with the experimental data set.

WOAPSO for Parameter Estimation of SS2018P PV Module
In order to further evaluate the efficiency of the proposed WOAPSO algorithm, parameters for SS2018P PV module were also estimated at different level of irradiance by utilizing the SDM model. The optimal value of five parameters (I p , I sd , a, R s , R sh ) for SDM of solar PV module at distinct levels of irradiance and constant temperature of 25 • C is presented in Table 4 and Tables S7-S9. The characteristics curve of current-voltage and power-voltage for solar PV module is redrawn based on best optimized parameters obtained by implementing the WOAPSO algorithm at a different level of irradiance, i.e., 1000 W/m 2 , 870 W/m 2 , 720 W/m 2 , and 630 W/m 2 and is depicted in Figure 5. It is found that the calculated data obtained by the WOAPSO is very effectively in keeping with the experimental data set. The curve of IAE between experimental and estimated values at 1000 W/m 2 , 870 W/m 2 , 720 W/m 2 , and 630 W/m 2 , is shown in Figure 6.

Convergence Analysis
To analyze the computational competence of WOAPSO, the convergence curves of the single-diode model, double-diode model, and PV module is presented in Figure 7. It is depicted that the proposed WOAPSO algorithm outperforms the GSA, SCA, GWO, PSO, WOA, PSOGSA algorithms in terms of convergence speed and generates a precise solution for the identical number of function evaluations (i.e., 50,000).

Convergence Analysis
To analyze the computational competence of WOAPSO, the convergence curves of the single-diode model, double-diode model, and PV module is presented in Figure 7. It is depicted that the proposed WOAPSO algorithm outperforms the GSA, SCA, GWO, PSO, WOA, PSOGSA algorithms in terms of convergence speed and generates a precise solution for the identical number of function evaluations (i.e., 50,000).

Robustness and Statistical Analysis
This section presents statistical evaluation based on mean, minimum, maximum, and standard deviation of RMSE for all previously implemented methods, and a comparison with respect to precision and consistency of the distinct algorithms in a total of thirty runs and depicted in Table 5. The mean of RMSE is calculated to evaluate the precision of algorithms, and the standard deviation is calculated to evaluate the consistency of the parameter estimation methods. The statistical results presented in Table 5 indicate that WOAPSO is the most accurate and reliable parameter optimization technique. As shown in Table 6, based on the Friedman ranking test result, the best ranking is obtained by the WOAPSO, followed by WOA, GWO, GSA, PSOGSA, SCA, and PSO. Also, Figure 8 shows the distribution of results (i.e., RMSE) obtained from the distinct algorithms in 30 runs in the form of a boxplot graph for the SDM, DDM, and PV module. It can be anticipated from Figure 8 that the proposed WOAPSO algorithm delivers the best results in terms of accuracy and reliability compared to the other six algorithms.

Discussion
To evaluate the reliability of the WOAPSO, the proposed hybrid algorithm is compared with six well established metaheuristics algorithms, i.e., GSA [56], SCA [57], GWO [58], PSO [59], WOA [37], PSOGSA [60] as well as other algorithms existing in the literature. It is observed that the estimated parameters based on the optimization process are highly consistent with the experimental data for SDM, DDM, and SS2018P PV module.
For SDM, the hybrid WOAPSO algorithm generates the lowest RMSE values (7.1700 × 10 −4 ) compared to the GSA, SCA, GWO, PSO, and WOA, PSOGSA algorithms ( Table 2). The RMSE of the proposed WOAPSO algorithm is also compared with previously studied algorithms (Table S2). It is noted that the hybrid WOAPSO algorithm provides the lowest RMSE values than that of others. Table S3 represents the absolute IAE for SDM analysis. The magnitude of IAE for different observations is less than 0.0018 (Table  S3), which indicates that the parameters optimized by the WOAPSO are very precise.

Discussion
To evaluate the reliability of the WOAPSO, the proposed hybrid algorithm is compared with six well established metaheuristics algorithms, i.e., GSA [56], SCA [57], GWO [58], PSO [59], WOA [37], PSOGSA [60] as well as other algorithms existing in the literature. It is observed that the estimated parameters based on the optimization process are highly consistent with the experimental data for SDM, DDM, and SS2018P PV module.
For SDM, the hybrid WOAPSO algorithm generates the lowest RMSE values (7.1700 × 10 −4 ) compared to the GSA, SCA, GWO, PSO, and WOA, PSOGSA algorithms ( Table 2). The RMSE of the proposed WOAPSO algorithm is also compared with previously studied algorithms (Table S2). It is noted that the hybrid WOAPSO algorithm provides the lowest RMSE values than that of others. Table S3 represents the absolute IAE for SDM analysis. The magnitude of IAE for different observations is less than 0.0018 (Table S3), which indicates that the parameters optimized by the WOAPSO are very precise.
For the SS2018P PV module, the hybrid WOAPSO algorithm produces the lowest RMSE values compared to the GSA, SCA, GWO, PSO, WOA, and PSOGSA algorithms. The IAE magnitudes for different observations (at 1000 W/m 2 ) are less than 0.0018 (Table  S6). More importantly, the computational time for WOAPSO is less than other algorithms ( Table 4). The average execution time of each algorithm on the three PV models is calculated and illustrated in Figure 9. The WOAPSO algorithm requires less time (about 26.1 s) than GWO, PSO, SCA, WOA, and PSOGSA, while GSA has the worst execution time of approximately 52 s. Table S5. It is noticed that the IAE values are less than 0.0097, which demonstrates the accuracy of optimized parameters produced by WOAPSO.
For the SS2018P PV module, the hybrid WOAPSO algorithm produces the lowest RMSE values compared to the GSA, SCA, GWO, PSO, WOA, and PSOGSA algorithms. The IAE magnitudes for different observations (at 1000 W/m 2 ) are less than 0.0018 (Table  S6). More importantly, the computational time for WOAPSO is less than other algorithms ( Table 4). The average execution time of each algorithm on the three PV models is calculated and illustrated in Figure 9. The WOAPSO algorithm requires less time (about 26.1 s) than GWO, PSO, SCA, WOA, and PSOGSA, while GSA has the worst execution time of approximately 52 s. Furthermore, the Friedman ranking test is also performed for all algorithms and depicted in Table 6. Table 6 shows that the proposed WOAPSO algorithm significantly outperforms the GSA, SCA, GWO, PSO, WOA, PSOGSA algorithms for all three models, i.e., single-diode, double-diode, and PV module models.

Conclusions
In this study, the hybridization of whale optimization and particle swarm optimization algorithm (WOAPSO) is anticipated. The exploitation ability of PSO is only implemented in pipeline mode when WOA stops to improve the best-found solution. The collaboration of both metaheuristic algorithms can establish an effective balance between exploitation and exploration ability. The proposed technique is further used to estimate the parameter of three PV cell models, i.e., single-diode, double-diode, and SS108P PV panel module model at different operating conditions. It should be noted that this suggested technique is, for the first time, intended to track the estimation of parameters for photovoltaic models reliably. The major conclusions are classified as follows: • The proposed WOAPSO is relatively accurate and reliable at delivering the solution in terms of RMSE as compared with other algorithms such as GSA, SCA, GWO, PSO, WOA, PSOGSA, and existing algorithms in the literature.

•
The I-V and P-V characteristic curves and IAE results indicate that WOAPSO can generate the optimized value of estimated parameters for all the models of solar PV cell as compared with other algorithms.

•
The statistical analysis clearly depicts the robustness of the proposed WOAPSO technique on parameter estimation problem at different operating conditions. Furthermore, the Friedman ranking test is also performed for all algorithms and depicted in Table 6. Table 6 shows that the proposed WOAPSO algorithm significantly outperforms the GSA, SCA, GWO, PSO, WOA, PSOGSA algorithms for all three models, i.e., single-diode, double-diode, and PV module models.

Conclusions
In this study, the hybridization of whale optimization and particle swarm optimization algorithm (WOAPSO) is anticipated. The exploitation ability of PSO is only implemented in pipeline mode when WOA stops to improve the best-found solution. The collaboration of both metaheuristic algorithms can establish an effective balance between exploitation and exploration ability. The proposed technique is further used to estimate the parameter of three PV cell models, i.e., single-diode, double-diode, and SS108P PV panel module model at different operating conditions. It should be noted that this suggested technique is, for the first time, intended to track the estimation of parameters for photovoltaic models reliably. The major conclusions are classified as follows:

•
The proposed WOAPSO is relatively accurate and reliable at delivering the solution in terms of RMSE as compared with other algorithms such as GSA, SCA, GWO, PSO, WOA, PSOGSA, and existing algorithms in the literature.

•
The I-V and P-V characteristic curves and IAE results indicate that WOAPSO can generate the optimized value of estimated parameters for all the models of solar PV cell as compared with other algorithms.

•
The statistical analysis clearly depicts the robustness of the proposed WOAPSO technique on parameter estimation problem at different operating conditions. • The convergence curves demonstrate that the best values of estimated parameters are obtained by WOAPSO, and RMSE is 7.1700 × 10 −4 and 9.8412 × 10 −4 in the case of single-and double-diode respectively. • At different irradiation levels (i.e., 1000 W/m 2 , 870 W/m 2 , 720 W/m 2 , and 630 W/m 2 ), the proposed WOAPSO algorithm is best in producing optimized parameters (I p , I sd , a, R s , R sh ) and minimum value of RMSE for PV module even at a low level of irradiation (630 W/m 2 ).
The WOAPSO is an efficient and robust technique to estimate the unknown optimized parameters of the solar PV model at different operating conditions. For future study, the implementation of proposed WOAPSO to solve the other problems related to energy optimization such as economic load dispatch, energy scheduling and optimization of PV array configuration may also be interesting for scientists and research scholars.