Abstract

In the past few decades, a lot of optimization methods have been applied in estimating the parameter of photovoltaic (PV) models and obtained better results, but these methods still have some deficiencies, such as higher time complexity and poor stability. To tackle these problems, an enhanced success history adaptive DE with greedy mutation strategy (EBLSHADE) is employed to optimize parameters of PV models to propose a parameter optimization method in this paper. In the EBLSHADE, the linear population size reduction strategy is used to gradually reduce population to improve the search capabilities and balance the exploitation and exploration capabilities. The less and more greedy mutation strategy is used to enhance the exploitation capability and the exploration capability. Finally, a parameter optimization method based on EBLSHADE is proposed to optimize parameters of PV models. The different PV models are selected to prove the effectiveness of the proposed method. Comparison results demonstrate that the EBLSHADE is an effective and efficient method and the parameter optimization method is beneficial to design, control, and optimize the PV systems.

1. Introduction

In recent years, in order to solve the problems of environmental pollution and burn out, the utilization of solar, wind, hydropower, nuclear, and so on has been increasing attention [1]. Among them, solar energy is considered as one of the most promising alternatives to inexhaustible and clean sources. At the present time, the PV systems play a very important role in power system, because solar energy can directly be converted into electric energy to supply power with humans. Therefore, solar PV systems have been widely applied in the whole world and have continued growing. However, the PV systems are exposed to the external environment and their PV arrays are prone to aging, which seriously affects production efficiency of PV panels and are harmful to the work efficiency of solar energy [2, 3]. Hence, in order to effectively design, simulate, estimate, control, and optimize PV systems, it is paramount to estimate the performance using exact model in operation. The most widely used mathematical models for describing the nonlinear behaviour and performance are the single and double diode models and the PV model. However, the accurateness of these models is dependent on the parameters’ values of models. The parameters are the reflection of the intrinsic characteristics of the PV model, and the I-V equation can be determined by identifying the PV parameters in order to predict the output power of PV array. The parameter optimization problem of the PV model is to fastly and accurately identify the parameters of the PV model in order to obtain better output power prediction and maximum power point tracking. Due to the aging, failure, breakdown, and unstable working conditions, it is very difficult to determine these parameters. Therefore, it is an extremely essential work to deeply study an effective method to optimize the parameters of PV and improve the solar energy utilizing efficiency.

Usually, the estimated problem of unknown parameters for the PV models is considered as an optimization objective function [4, 5]. Due to the measured current data and voltage data involving noise, the constructed objective function is a nonlinear and multimodal function with multiple local optimums. Some researchers have studied and presented many methods to optimize the parameters of PV models in recent years, such as analytical method, deterministic method, and heuristic method. The heuristic method is a promising alternative to analytical method and deterministic method [6]. Due to no strict restrictions for the objective function, a lot of heuristic methods have obtained more and more attention for optimizing the parameters of PV models, such as simulated annealing (SA), genetic algorithm (GA), differential evolution (DE), particle swarm optimization (PSO), artificial bee colony (ABC), bacterial foraging algorithm (BFA), teaching-learning-based optimization (TLBO), whale optimization algorithm (COA), bird mating optimizer (BMO), month flame optimizer (MFO), and backtracking search algorithm (BSA) [720]. Recently, gaining-sharing knowledge-based algorithm has been proposed by Mohamed et al. [21]. Zagrouba et al. [22] used genetic algorithms to identify the parameters of PV solar cells and modules. Ishaque et al. [23] presented a penalty-based DE to extract the parameters of solar PV models. Merchaoui et al. [24] presented an improved PSO with adaptive mutation strategy to extract the parameters of PV solar cell/module. Yu et al. [25] presented a performance-guided JAYA algorithm to identify the parameters of PV cell and modules. Chen et al. [26] presented a teaching-learning-based ABC algorithm to estimate the parameters of the solar PV model. Li et al. [27] presented an improved teaching-learning-based optimization algorithm to extract the parameters of PV models. Oliva et al. [28] used ABC algorithm to identify the parameters of solar cells. Wu et al. [29] presented an improved ant lion optimizer to identify the parameters of the PV cell model. Ram et al. [30] presented a new hybrid bee pollinator flower pollination algorithm to estimate the parameters of solar PV. Oliva et al. [31] presented an improved chaotic WOA to estimate the parameters of PV cells. Xu and Wang [32] presented a hybrid flower pollination algorithm to estimate the parameters of PV modules. Chen et al. [33] presented a diversification-enhanced Harris Hawks optimizer to identify the parameters of PV modules. Long et al. [34] presented a new hybrid algorithm based on a grey wolf optimizer and cuckoo search to extract the parameters of solar PV models. Chen et al. [35] presented a novel opposition-based sine cosine approach with local search to estimate the parameters of PV models. Cai et al. [36] employed an orthogonal experiment method to optimize the parameters of the dust absorbing structure for PV panels. The results reported by these heuristic methods have achieved satisfied results, which indicate that these heuristic methods are promising alternative to identify the parameters of PV models. But because the parameter identification of PV models is a nonlinear and multimodal problem, their accuracy and reliability need to be further improved. In addition, many heuristic methods have their own parameters to be experimentally tuned to improve their efficacy, accuracy, reliability, and scalability.

Differential evolution (DE), a random evolution algorithm based on population evolution, was proposed by Storn and Price in 1997 [37], which has been regarded a simple and efficient optimization algorithm. Due to the fast convergence, robustness, and search ability, various advanced DE variants have been presented and widely applied in various fields of engineering design, operations research, biomedicine, and so on. Most of the various advanced DE variants are used to optimize the parameters of different PV models. Chang [38] presented an ant direction hybrid differential evolution algorithm. Gong and Cai [39] presented an improved adaptive differential evolution with crossover rate repairing technique and ranking-based mutation, as R-cr-IJADE. Muhsen et al. [40] presented an improved differential evolution with adaptive mutation per iteration algorithm (DEAM). Chellaswamy and Ramesh [41] presented adaptive differential evolution algorithm. Mohamed and Abdulaziz [42] proposed a differential evolution with novel mutation and adaptive crossover strategies to solve large-scale global optimization problem. Zhang et al. [43] presented a gradient decent-based multiobjective cultural differential evolution (GD-MOCDE) to improve the optimal efficiency. Mohamed [44] proposed an enhanced adaptive differential evolution algorithm to solve large-scale global optimization problems. Rashidi and Khorshidi [45] presented a multiobjective differential evolution algorithm. Mohamed [44] presented a multiobjective self-adaptive differential evolution (MOSaDE) algorithm. Mohamed [46] proposed a new approach to differential evolution algorithm for solving stochastic programming problem. Ramli et al. [47] presented a multiobjective self-adaptive differential evolution algorithm. Mohamed and Suganthan [48] proposed enhanced fitness-adaptive differential evolution algorithm with novel mutation to solve real-parameter unconstrained optimization problem. Mohamed and Mohamed [49] proposed an adaptive guided directed differential evolution algorithm to solve unconstrained optimization problems. Tey et al. [50] presented an improved global search space differential evolution algorithm. Xiong et al. [51] presented an effective hybrid method based on the exploration of DE with the exploitation of WOA as DE/WOA. Hadi et al. [52] proposed a LSHADE-SPA memetic framework to solve large-scale optimization problem. Mohamed and Mohamed [53] proposed an enhanced AGDE algorithm for real-parameter unconstrained optimization problem. Li et al. [54] presented a memetic adaptive differential evolution as MADE algorithm. Mohamed et al. [55] proposed an enhanced DE algorithm (EDDE) that utilizes the information given by good individuals and bad individuals in the population. The new mutation scheme maintains effectively the exploration/exploitation balance. Essiet et al. [56] presented an improved enhanced differential evolution algorithm for implementing demand response between the aggregator and consumer. In the team of professor Mohamed, a lot of work has been done to improve the DE algorithm. Mohamed et al. [53] proposed an enhanced directed differential evolution algorithm for solving constrained nonlinear integer and mixed-integer global optimization problems. Mohamed et al. [57] proposed an EBLSHADE algorithm based on novel mutation strategy. In this work, two mutation operators are introduced, ord_best and ord_pbest, which are versions of the classical DE/current-to-best/1 scheme. The proposed mutations were incorporated into SHADE and LSHADE algorithms in order to enhance their performances.

Due to the better robustness, stability, and quality of the solution of LSHADE, an enhanced SHADE (EBLSHADE) algorithm is applied to propose a parameter optimization method to optimize the parameters of PV models quickly, accurately, and reliably. In the EBLSHADE, a less and more greedy mutation strategy is used to enhance the exploitation capability and the exploration capability. A linear population size reduction strategy is used to gradually reduce population for improving the convergence speed, balancing the exploration and exploitation capabilities in the process of evolution. The EBLSHADE is employed to propose a parameter optimization method to optimize the unknown parameter estimation of the single diode model, double diode model, and PV model. The experimental results demonstrate that the parameter optimization method can exactly and reliably optimize the unknown parameters of different PV models and provide highly competitive results compared with other algorithms.

The main contributions of this paper are described as follows:(i)An EBLSHADE is applied to effectively optimize the parameters of PV models. Based on quantified performance, the proper update strategy can adaptively be selected for individuals to immeasurably elevate the related searching performance.(ii)The less and more greedy mutation strategy in the EBLSHADE is used to enhance the exploitation capability and the exploration capability, respectively.(iii)The linear population size reduction strategy is employed to gradually reduce population to further improve the convergence speed, balance the exploration and exploitation capabilities, and avoid to falling into local optimum in the process of evolution.(iv)The performance of EBLSHADE has been extensively investigated by the parameter estimation of different PV models.

The rest of this paper is arranged as follows. Section 2 describes different PV models and their objective functions. The differential evolution algorithm is briefly described in Section 3. Section 4 introduces an enhanced SHADE with multistrategies in detail. The experimental results on different PV models are shown and analyzed in Section 5. Finally, the conclusion is given in Section 6.

2. Differential Evolution

The main operations of DE contain initialization population, mutation operation, crossover operation, selection operation, and so on. Its main thoughts are to differentiate and scale between two different individual vectors in the same population and add a third individual vector in this population to obtain a mutation individual vector, which is crossed with the parent individual vector with a certain probability to generate an attempted individual vector. Finally, the attempted individual vector and the parent individual vector are executed greedy selection, and the better individual vector is saved to the next generation population. The basic evolution processes of e DE are described.

2.1. Abbreviations and Acronyms

The DE uses NP D-dimension vectors as the initial solution. Setting population number N, each individual can be expressed as . The initial population is generated in .where G represents the Gth generation, represents the maximum search space value, represents the minimum search space value, and represents a random number.

2.2. Mutation Operation

It generates a mutation vector for , namely, target vector. For each generated target vector, the mutation strategy is used to obtain a corresponding mutation vector. The mutation operation is the most important operation of DE. According to the different generation methods of the mutation individuals, several different mutation strategies of DE are formed. The most common mutation strategies are described as follows:(1)DE/best/2/bin(2)DE/rand/2/bin(3)DE/current-to-best/1/bin(4)DE/current-to-rand/1/bin

2.3. Crossover Operation

Each pair of target vectors and the corresponding mutation vectors are crossed to obtain a test vector . In the basic DE, it uses a binomial crossover to define.where the CR is a constant between 0 and 1, which is used to control the duplicated proportion from the mutation vector. is a selected integer within [1, D], .

2.4. Selection Operation

Comparing the objective function value of each test vector, the objective function value of the test vector is less than the corresponding target vector, and then the target vector is replaced by the experimental vector. The selection operation is given as follows:

2.5. Implementation of DE

The DE algorithm evolves generation by generation until the result or ending condition has been met. The flow chart of DE algorithm is shown in Figure 1.

3. EBLSHADE

3.1. SHADE

The SHADE is one of the most successful variants of DE [57]. In the SHADE, a historical memory with H entries is used. It is made up of MCR and MF, which can adaptively control parameters CR and F. In each iteration, each individual (xi) has its own Fi and CRi to generate a new test vector ui. The two parameters are expressed as follows:where ri is a random integer on [1, H], randn is the Gaussian normal distribution, and the scale parameter is 0.1, and randc is the Cauchy distribution and the variance is 0.1.

The CR and F of the generated test vectors are recorded as SCR and SF. Their average values are stored in MCR and MF. The SHADE keeps H parameters to guide the control parameters in order to achieve self-adaptive search. Even if the SCR and SF of some offsprings contain a set of poor values, the stored parameters of the previous generation will not be affected. The control parameters are shown in Table 1.

Then, the MCR and MF are updated according to the following expressions:where t is an index to determine the saved position and k is the set 1. When a pair new SCR and SF are added in the history, k is incremented by 1. MeanWA (SCR) and MeanWL (SF) are weight mean values of SCR and weight Lehmer mean value of SF.

3.2. LSHADE

An improved SHADE based on the reduction strategy of the linear population size, namely, LSHADE is developed. In the LSHADE, the population size is gradually lessened. Therefore, the linear function is described as follows:where round returns the nearest integer number. NFE is the fitness optimization, MaxNFE is the optimization with the maximum iterations, NP0 is the initial population size, and NPmin is the possible minimum population size, which is the minimum number of individuals (NPmin = 4).

3.3. EBLSHADE
3.3.1. A Mutation Strategy

The DE/current-to-best/1 strategy can find the best solution in the evolution. But it may deteriorate or lose the population diversity and exploration capability. To overcome these shortcomings, a variant of DE/current-to-best/1, namely, DE/current-to-or_best/1 is used to balance the local exploitation and the global exploration abilities and enhance the convergence speed. In the strategy, all individuals are sorted to divide into three vectors. The best vector is referred as xor_best, G, the median vector is referred as xor_median, G, and the worst vector is referred as xor_worst, G. Therefore, the trial vector is described as follows:

As can be seen from the new mutation strategy equation (12), the added objective function value has two advantages. The difference vector of the best vector and the target vector is the first perturbation part of the new mutation strategy, which can substantially avoid prematurity and accelerate convergence. The second perturbation part of the new mutation strategy is the difference vector of the median vector and the worst vector. Therefore, the DE/current-to-or_best/1 can get the global optimal solution.

When the population size reduction strategy can improve optimization performance of the algorithm, the initial population size will be reduced to 18 dimensions in the LSHADE. The increased population size will affect DE/current-to-or_best/1; the probability of xor_best, G, will be decreased to be the global best solution. Therefore, the behaviour of the DE/current-to-or_best/1 will approximate to DE/rand/1. To solve this problem, an enhanced version of DE/current-to-or_pbest/1 is used. In this mutation strategy, one vector from top p best vectors is included. The other two vectors are chosen randomly. Then, three vectors are sorted. The best vector is referred as xor_pbest, G, the median vector is referred as xor_pmedian, G, and the worst vector is referred as xor_pworst, G. Therefore, the trial vector is obtained as follows:

3.3.2. Model of EBLSHADE

The flow of the EBLSHADE is shown in Figure 2.

4. Modeling of PV Models

4.1. SDM

This model can describe the behaviour of solar cell effectively [17]. Equivalent circuit is given in Figure 3.

In this model, it is made up of a leakage current shunt resistance, a parallel current source with a diode, and several resistors for the related load current loss. Therefore, the output current is described as follows:

From the abovementioned equations, the parameters of Iph, Id, Rs, Rsh, and α need to be optimized. Therefore, the objective functions can be formulated:

4.2. DDM

It is used to take the recombination current loss effect. Equivalent circuit is given in Figure 4.

The output current of this model is described as follows:

From the abovementioned equations, the parameters of Iph, Isd1, Isd2, Rs, Rsh, α1, and α2 need to be estimated. Therefore, the objective functions can be formulated:

4.3. PVM

It is usually based on series or parallel solar cells. Its equivalent circuit is described in Figure 5. This circuit can be extended to the N-diode PV model.

The output current of this model is formulated.

Due to the used PV model with series in the experiment, there is NP = 1. The abovementioned equation can be reformulated:

From the abovementioned equations, the parameters of Iph, Id, Rs, Rsh, and α will be optimized. The objective functions of the PV model can be formulated:

5. Experiment Results and Analysis

5.1. Data and Environment

The benchmark experimental current-voltage data of SDM (Ns = Np = 1) and DDM (Ns = 1, Np = 2) are used from reference [58]. The current-voltage data came from references [59, 60]. The experimental platform is a PC with an Intel i7-7700HQ, 16 GB, under Windows 10.

5.2. Parameter Setting

In order to guarantee the comparative fairness, the lower bound (LB) and upper bound (UP) are described in Table 2.

To prove the superior ability of the EBLSHADE, some existing algorithms are selected in Table 3. The parameters of the EBLSHADE are set as follows: NP = 40, H = 100,  = 0.5ones (1, H),  = 0.5ones (1, H),  = 0.2,  = 0.6, pmin = 0.05, and pmax = 0.2. The settings of parameters are the same as the corresponding literatures for all compared algorithms. These algorithms are executed in Matlab2018b and 30 runs.

5.3. Root Mean Square Error

Root mean square error (RMSE) is the ratio square root of the square of the deviation between the predicted value and the actual value and the number of observations n. In the actual measurement, the observation number n is always limited, and the true value can only be replaced by the most reliable value. It is used to measure the deviation between the observed value and the true value and illustrates the dispersion degree of the data. The root mean square error (RMSE) is defined as follows:where N is the experiment number and x is a vector.

5.4. Experiment Results
5.4.1. Experimental Results of SDM

Here, the RMSE is the obtained computation value according to the equation of root mean square error (equation 20) by each independent run of the algorithm. MRMSE is the mean value of RMSE of different set of 30 runs; BRMSE is the best value of RMSE among different set of 30 runs. Therefore, the compared results of five parameters, statistical results of the best RMSE (BRMSE), the mean RMSE (MRMSE), the standard deviation (SD), and the least computing resources (NFE) are given in Table 3, where the obtained best results are bold. The boxplot of RMSE is presented in Figure 6.

From Table 3 and Figure 6, the TLABC, MLBSA, JADE, SHADE, MADE, and EBLSHADE obtain the BRMSE value (9.8602E − 04). Especially for MLBSA, SHADE, MADE, and EBLSHADE, the best, worst, and mean RMSE values are the same value (9.8602E − 04). Due to the unavailable information, the RMSE is usually used to express the more accurate parameters. Although the second BRMSE value (9.8603E − 04) of IJAYA is infinitely close to the BRMSE value (9.8602E − 04), it is significant to reduce the order of objective function and improve the true value of parameters. For standard deviation, the EBLSHADE obtains the third best standard deviation value, which is close to the best standard deviation value. For max NFE, the EBLSHADE required max NFE = 4,000, which is much less than the other compared algorithms.

To further prove the effectiveness of the EBLSHADE, the detailed results of the individual absolute error of current (IAEI) and power (IAEP) between the experimental data and the measured data are shown in Table 4. See references for measurement data. The I-V and P-V characteristics curve obtained by EBLSHADE are shown in Figures 710. Note that Ic is calculated by the optimized parameters using EBLSHADE.

From Table 4 and Figures 710, all IAEI values are lesser than 2.5075E − 03 and all IAEP values are lesser than 1.4626E − 03. It can be seen that the EBLSHADE can accurately optimize these parameters. It is also evident that the obtained experimental data by EBLSHADE are highly consistent with the measured data, which effectively reflects the optimized parameters to be accurate enough. Therefore, the EBLSHADE can be considered as a significant method for the parameter optimization of SDM.

5.4.2. Experimental Results of DDM

Seven optimized parameters will increase the optimization difficulty. The EBLSHADE is compared with the selected methods in Table 5. The results are shown in Table 6. The boxplot of RMSE is shown in Figure 11.

From Table 5 and Figure 11, the SHADE only achieved the BRMSE value (9.8248E04), trailed by SATLBO, IJAYA, MLBSA, MADE, and EBLSHADE. The obtained BRMSE value (9.8295E04) by EBLSHADE is close to the obtained BRMSE value(9.8248E04). By comparison, the RMSE values of the other algorithm are lower, especially the obtained result by BLSPO is worst. While the SHADE consumed 20000 NFE, although the EBLSHADE did no’t obtain the BRMSE value, it only consumed 10000 NFE to obtain the better competitive results. Nevertheless, other compared algorithms need much more computing resources except for MADE. It is worth knowing that the MLBSA, JADE, TLABC, and MADE did not reach the BRMSE values. For standard deviation, the EBLSHADE obtained the best SD value (1.2825E06), which has much less computing resources than the other selected algorithms, except for MADE.

To further prove the accuracy of the parameters optimized by EBLSHADE, detailed results of the IAEI and IAEP are shown in Table 6. The obtained I-V and P-V curves are shown in Figures 1215. See references for measurement data.

From Table 6 and Figures 1215, all IAEI values are lesser than 2.5070E03 and all IAEP values are lesser than 1.4623E03, which prove that the EBLSHADE can obtain the optimized parameters with higher accuracy. The experimental results show that the experimental data by EBLSHADE are highly in agreement with the measured data. So, the EBLSHADE is regarded as a significant algorithm to optimize the parameters of DDM.

5.4.3. Experimental Results of PVM

For the PV models, to further prove the performance of the EBLSHADE, photowatt-PWP201, STM6-40/36, and STP6-120/36 are used in here. The experimental results of the PVMs are shown in Tables 79. The boxplot of RMSE are presented in Figures 1618.

From Table 7 and Figure 16 on the photowatt-PWP201 model, all compared algorithms can obtain the same BRMSE (2.4251E03), while the EBLSHADE only consumed the least 5,000 NFE. Especially for the SHADE, MADE, PGJAYA, and EBLSHADE, the best, worst, and mean RMSE values are the same value (2.4251E03), which indicate that the SHADE, MADE, PGJAYA, and EBLSHADE take on better stability. For standard deviation, the EBLSHADE obtained the second best standard deviation value (2.8821E17), which is close to the best standard deviation value (2.0700E17). Although the EBLSHADE did no’t obtain the best standard deviation value, it only consumed 5,000 NFE to obtain the better optimization results.

From Table 7 and Figure 17 on the STM6-40/36 model, except for IJAYA and JADE, the other compared algorithms can get the same BRMSE value (1.7298E03). Especially for MADE and EBLSHADE, the best, worst, and mean RMSE values are the same value (1.7298E03), which indicate that the MADE and EBLSHADE take on better stability. For standard deviation, the EBLSHADE obtained the best standard deviation value (6.40591E14). The max NFE of EBLSHADE is 10,000 NFE, which utilizes less computational resources compared to other compared algorithms except for MADE. Although the EBLSHADE is not less NFE for the STM6-40/36 model, it consumed 10,000 NFE to obtain the better optimization results.

From Table 8 and Figure 18 on the STP6-120/36 model, except for IJAYA and JADE, the other compared algorithms can get the same BRMSE value (1.6601E02). Especially for SATLBO, MADE, and EBLSHADE, the best, worst, and mean RMSE values are the same value (1.6601E02), which indicate that the SATLBO, MADE, and EBLSHADE take on better stability. For standard deviation, the EBLSHADE obtained the best standard deviation value (8.0544E16). The max NFE of EBLSHADE is 15,000 NFE, which are much less computational resources than the other compared algorithms except for MADE. Although the EBLSHADE is not less NFE for the STP6-120/36 model, it consumed 15,000 NFE to obtain the better optimization results.

Therefore, it is clear that the EBLSHADE can effectively and consistently provide better results of integrating the RMSE and NFE considerations by comparing with other methods.

To show the accuracy of the parameters optimized by EBLSHADE, detailed results of the IAEI and IAEP for three PVM are shown in Tables 1012. The obtained I-V and P-V curves are given in Figures 1930. See references for measurement data.

From Tables 1012 and Figures 1930, all IAEI values are lesser than 4.8329E03 and all IAEP values are lesser than 3.7363E02 for the photowatt-PWP201 PV model. For the STM6-40/36 PV model, all IAEI values are lesser than 6.0901E03 and all IAEP values are lesser than 9.0620E01. For the STP6-120/36 PV model, all IAEI values are lesser than 4.4318E02 and all IAEP values are lesser than 7.1264E01. In general, the EBLSHADE can effectively optimize parameters of the PV models with higher accuracy. It is also clear that the experimental data by EBLSHADE are highly in agreement with the measured data. Therefore, the EBLSHADE is considered as a significant method to optimize the parameters of the PV models.

5.4.4. Discussion of the Results

As the experiment results are given in Section 5.3, the EBLSHADE is applied in the parameter optimization of PV models, respectively. The superiority of the EBLSHADE has been proved by comparing with some other algorithms. It can be seen from the I-V curves and P-V curves that the experimental results of EBLSHADE are highly in agreement with the measured data for all datasets, which effectively reflects the optimized parameters to be accurate enough by using EBLSHADE. All the compared results containing the optimal parameters, the BRMSE, WRMSE, MRMSE, SD, and NFE indicate that the EBLSHADE has better capacities of the exploration and exploitation. The reason is that the novel mutation strategies of DE/current-to-or_best/1 and DE/current-to-or_pbest/1 for EBLSHADE can enhance and balance the exploitation capability and exploration capability and enhance the convergence. Rregarding the whole result, the EBLSHADE is considered as a significant method for the parameter optimization of SDM, DDM, and PVM. It is also expected to be applied for optimizing parameters of other PV models.

6. Conclusion

In this paper, an enhanced success history adaptive DE(EBLSHADE) has been applied to propose a parameter optimization method for deferent PV models. The effectiveness of the EBLSHADE and parameter optimization method has been verified on estimating parameters of different PV models. The experimental and statistical results with the reliability, accuracy, and computing efficiency show that the EBLSHADE is superior to the other compared algorithms, and the parameter optimization method is an effective method and can design, control, and optimize the PV systems. Therefore, the EBLSHADE is considered as a significant method for optimizing parameters of other PV models.

Data Availability

The data are available in this paper.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (61771087), the China Education Ministry Humanities and Social Science Research Youth Fund Project (18YJCZH192), the CERNET Innovation Project (NGII20190605), Yantai Key Research and Development Program (2020YT06000970), Sichuan Science and Technology Program (2019YFG0216 and 2019ZYZF0169), and the A Ba Achievements Transformation Program (19CGZH0006).