Abstract

In this paper, we propose IJADE-TSO, a novel hybrid algorithm in which an improved adaptive differential evolution with optional external archive (IJADE) has been combined with the tuna swarm optimization (TSO). The proposed algorithm incorporates the spiral foraging search and parabolic foraging search of TSO into the mutation strategy in IJADE to improve the exploration ability and population diversity. Additionally, to enhance the convergence efficiency, crossover factor (CR) ranking, CR repairing, top α r1 selection, and population linear reduction strategies have been included in the algorithm. To evaluate the superiority of the proposed algorithm, IJADE-TSO has been benchmarked with its state-of-the-art counterparts using the CEC 2014 test set. Finally, to check the validity of IJADE-TSO, we apply it to photovoltaic (PV) parameter identification and compare its performance with those of other recently developed well-known algorithms. The statistical results reveal that IJADE-TSO outperforms the other compared algorithms.

1. Introduction

There are an increasing number of optimization problems regarding engineering design and industrial production calling for more effective optimization methods in the real world [1, 2]. One efficient approach is the stochastic heuristic algorithm, which can be divided into three categories [3]: evolution-based algorithms [4, 5], physical-phenomena-based algorithms [6, 7] and swarm-based algorithms [8, 9]. Differential evolution (DE) [10], one of the most excellent heuristic algorithms, has been widely used to solve complex optimization problems such as mechanical engineering [11], optimal power flow [12], parameter optimization [13], neural network training [14, 15], and complex control system [16]. The simplicity and efficiency of the DE algorithm has been attracting extensive exploration of it since its proposition by Storn and Price in 1995.

The performance of DE algorithm is sensitive to the mutation strategy and parameters and prone to fall into a local optimum. To fix these defects, a series of DE variants were proposed over the past two decades. These DE variants can be easily traced through the result of the competition organized by IEEE Congress on Evolutionary Computation (CEC), where DE-based variant algorithms have ranked the top 3 in successive years except for 2013. According to the ways of improvement, current improvement methods can be classified into three main categories: optimizing mutation scheme, adaptive setting of control parameters, and hybridization.

1.1. Optimizing Mutation Scheme

Zhang and Sanderson [17] proposed JADE which implements the current-to-best mutation strategy for the first time with adaptively set CR and mutant factors (F) achieved through evolving the mutation factors and crossover probabilities based on their historical record of success. Based on current-to-pbest strategy, V. Stanovov et al. [17] developed current-to-pbest/r, a new mutation strategy in LSHADE-RSP, which performs better in comparison with the alternative algorithms.

1.2. Adaptive Setting of Control Parameters

As an improvement upon the robustness of JADE, SHADE [18] enhances the adaptive performance of control parameters by adopting a diverse set of parameters to guide control parameter adaptation based on the historical memory of successful parameter settings of JADE. As an extension of the SHADE algorithm, LSHADE [19] incorporates the linear population size reduction strategy to balance the exploration and exploitation. iL-SHADE, an improved version of LSHADE, was proposed by J. Brest et al. [20], in which the memory update mechanism is modified from LSHADE and the initial value of CR is set to 0.8 to enhance population diversity.

1.3. Hybridization

Hybridization is another important strategy to modify DE as plentiful hybridized algorithms concerning DE have been developed. C. Zhang et al. [21] proposed DE-PSO algorithm, where DE is hybridized with PSO for enhancing the exploration ability with most new individuals generated by the modified DE operator. However, the algorithm was only tested on the basic benchmark functions. S. Y. Du and Z. G. Liu [22] proposed PSOJADE by introducing a multicrossover operation and JADE to enhance the global exploration and the local exploitation in PSO. W. Gong et al. [23] proposed DE/BBO, a hybrid DE with BBO, which combines the exploration of DE with the exploitation of BBO effectively, enabling it to generate the promising candidate solutions. A. W. Mohamed et al. [24] proposed a hybridization framework in which CMA-ES is hybridized with LSHADE through adaptive approach, and semiparameter adaptation is adopted to effectively adapt the values of F.

Nevertheless, conventional current-to-pbest mutation strategy in the above hybrid algorithms can only capture the interindividual differential information to form the trial vector, which is single and prone to stagnation. Therefore, it is judicious to integrate it with other algorithms with strong global exploitation capabilities to improve its performance. However, roughly hybridizing swarm-based algorithms with DE rarely performs well. Hence, it is crucial to find an effective way to hybridize the DE algorithm with other algorithms.

Tuna swarm optimization [25] is one of the excellent swarm-based algorithms proposed by L. Xie et al. It shows a strong global search ability due to spiral foraging and parabolic foraging search behavior. This paper introduces TSO into the mutation strategy instead of simply hybridizing TSO with IJADE. In the iteration of each generation, individuals in the population except the optimal individual are selected with a certain probability to generate offspring individuals by IJADE or TSO, respectively. Similarly, the mutation vectors generated by TSO are subject to crossover operations. In addition, a series of improvements are made to improve the convergence efficiency.

The main contributions of this paper are as follows: (1) A novel hybrid mutation strategy is proposed. In the mutation operator, spiral foraging and parabolic foraging are incorporated into classical “current-to-pbest” mutation strategy. (2) To enhance the convergence efficiency of JADE, CR sorting mechanism and CR repairing are introduced. (3) A top α r1 selection strategy is proposed to accelerate convergence efficiency. (4) The linear population size reduction strategy is introduced to approximate the optimal individual in the final phase of the search.

2. Brief Review of DE and JADE

Optimization problems are common in the real world. The common form is shown as follows [10]:

We have that is the optimization vector, D is the dimension of the optimization problem, and and are the boundaries.

2.1. Differential Evolution

At the beginning of the optimization problem, DE stochastically generates populations in the search space. The individuals in the population create the next generation in an evolutionary manner. When the individual explores a new location with a better fitness value, the individual moves to that location. There are four main operators in DE, that is, initialization, mutation, crossover, and selection operations. These operations will be discussed in detail as follows.

2.1.1. Initialization

Like other swarm intelligence optimization algorithms, the population is initialized first.where is the ith individual, j represents the jth dimension, and NP is the initial population number.

2.1.2. Mutation

For each individual in the population, the mutant vector of DE is generated in the following way:

In the above formula, , , and are randomly selected in the population and . Parameter F is used to control the amplification of the difference vector and .

In addition, the mutant vector can be generated in other ways [26].

2.1.3. Crossover

The trial vector is obtained through replacing some components of target vector with corresponding mutant vector .

randi(1, D) generate a random integer between 0 and D. is the crossover factor which decides the proportion of replaced components in .

2.1.4. Selection

In the selection operation, according to the greedy strategy, the individual of next generation is selected by comparing the trail vector and the target vector in DE. The selection method is as follows:

2.2. JADE

The JADE algorithm is a variant of DE that was proposed by Zhang J. et al. [27] in 2009. In recent years, many excellent algorithms have been established based on JADE, and their advantages are as follows.

2.2.1. DE/Current-to-pbest Strategy

The most important improvement is that the new mutation strategy DE/current-to-pbest is implemented with the optional external archive, which uses historical information to provide evolution direction information.

Mutation strategy without external archive is

Mutation strategy with external archive is

The difference between the above two mutation strategies is the element . In the strategy with external archive, there is a certain possibility for to be selected from either the current group or the external archive, which enhances the diversity of the population.

2.2.2. Adaptive Control Parameters of F and CR

The control parameters F and CR are updated adaptively to improve the optimization performance.

CR is generated according to the normal distribution with mean and standard deviation equal to 0.1, and F is generated according to the Cauchy distribution with mean and standard deviation equal to 0.1. and are initially set to 0.5.

SCR is a collection for storing the CR value of each generation of successful individuals. Similarly, SF is a collection for storing the F value of each generation of successful individuals.

After each iteration, the values of and are updated as follows:

We have that c is the learning rate, which is set to 0.1, meanA(.) refers to the arithmetic mean, and meanL(.) represents the Lehmer mean, which is defined as follows:

In general, the control parameters of the DE algorithm are not adaptive [28] and the convergence performance is undesired. Although the convergence speed of JADE is high [29], the single mutation strategy results in poor population diversity and liability of falling into local optimum.

3. IJADE-TSO

3.1. Tuna Swarm Optimization

The JADE algorithm uses a single mutation strategy, giving rise to a greater possibility of the algorithm falling into a local optimum. To compensate for this, two foraging search strategies in tuna swarm optimization (TSO) [25] have been introduced into the mutation operation of IJADE. The two mutation strategies account for a respective percentage of the population to improve population diversity and avoid local optimum.

Tuna swarm optimization is a novel swarm-based metaheuristic algorithm for global optimization. The inspiration for TSO comes from the cooperative foraging behavior of tuna swarm. It mainly consists of two foraging behaviors of tuna swarm: spiral foraging and parabolic foraging. Its global exploration capacity excels its exploitation capacity.

3.1.1. Spiral Foraging

When foraging, the tuna group will form a spiral formation to round up the target in the face of the rapid position change of the target fish group. If a small part of tuna group moves firmly towards a certain direction, the surrounding tuna group will gradually adjust their directions and follow to round up. Each tuna will follow the previous tuna while chasing the target, forming a chain of information transmission. The position of the most important lead tuna is updated by the mutation strategy of IJADE, which augments the foraging efficiency. The mathematical model of the spiral foraging strategy is as follows:

We have that is the ith individual of the +1 generation, is the current optimal individual (lead tuna), and are weight coefficients that control the tendency of individuals to move towards either the optimal individual or the previous individual, a is a constant used to determine the extent to which the tuna follows the optimal individual or the previous individual in the initial phase, denotes the number of current iterations, is the maximum number of iterations, and b is a random number uniformly distributed between 0 and 1.

In the initial stage of foraging, the prey tracked by the lead fish is not necessarily the fish group that provides the most abundant food. Therefore, the whole tuna group needs to scatter to find the location of the fish group with the most abundant food. It is more efficient to stochastically select a point in the space and spirally search around it. The specific mathematical model is described as follows:

We have that is a randomly generated reference point in the search space.

Metaheuristic algorithms usually perform extensive global exploration in the early stage and then gradually transition to precise local exploitation. Therefore, TSO changes the reference points of spiral foraging from random individuals to optimal individuals with the increase in iteration. In summary, the final mathematical model of the spiral foraging strategy is as follows:

3.1.2. Parabolic Foraging

To prevent preys from escaping, in addition to spiral foraging, tunas also adopt parabolic foraging. Tuna forms a parabolic formation with prey as a reference point. Meanwhile, tuna forage preys by searching areas around themselves. The two approaches are performed simultaneously, with the assumption that the selection probabilities are 50% for both. The specific mathematical model is described as follows:where TF is a random number with a value of 1 or −1.

Tuna work cooperatively with the two foraging strategies to hunt for the prey. In each iteration, each individual randomly executes one of the two foraging strategies.

3.2. Crossover Rate Sorting Mechanism and CR Repairing

In order to establish the relationship between CR and the individual fitness values, the CR sorting mechanism [30] is introduced. Firstly, the CR values are generated by Gaussian distribution and then are sorted in an ascending order. This is shown as follows:

A variant of JADE is introduced in [31], which modifies the crossover factor according to the real crossover rate of each generation. In this way, the adaption of the crossover factor is improved. Its crossover operation becomes

By sorting the CR values, the individuals with better fitness are given smaller CR, so their next generation can retain more characteristics of the parent individuals. Meanwhile the poor individuals will be given larger CR and a larger proportion of components will be replaced by the mutated individuals. This helps to improve the exploration efficiency.

With the CR sorting mechanism, the issue of neglecting the relationship between CR values and the individual fitness values in JADE can be effectively addressed. Through CR correction, the actual crossover rate is introduced into the calculation of the crossover factor in the next generation.

3.3. Top α r1 Selection

In LSHADE-RSP [17], a ranking-based approach was proposed for the selection of r1 and r2. In the JADE algorithm, the selection of the r1 individual is random. To improve the convergence efficiency of the algorithm, we use the top α r1 selection strategy. The selection of r1 is shown as follows:

We have that is the number of candidates for the selection of r1, and rand is a random value selected in [0, 1]. The individuals with better fitness values will be more likely to be selected. In this way, it is easier to form a difference vector that evolves towards the current optimal individual and expedite convergence.

3.4. Linear Population Size Reduction

The linear population size reduction strategy is introduced to accomplish depth exploration around the optimal individual and facilitate algorithm convergence in the final stage of the search. The mathematical formula is as follows:

NPmin and NPinit represent the minimum population and the number of initial populations, respectively. However, the increasing number of populations will result in a decrease in population diversity, making it liable to be trapped in local optimum. Therefore, rather than setting NPmin to 4 in [19], NPmin is set to 50 in this paper.

The algorithm pseudocode is shown as Table 1 and the corresponding flowchart of IJADE-TSO algorithm is shown in Figure 1.

4. Numerical Experiment and Discussion

4.1. CEC 2014

In this section, we use the 2014 IEEE CEC test suite to verify the performance of the IJADE-TSO by comparing it with state-of-the-art algorithms, including LSHADE [19], iLSHADE [20], SPS-LSHADE-EIG [32], CPI-JADE [33], GEDGWO [34], and DOLTLBO [35]. Among them, LSHADE won a prize at the 2014 IEEE CEC. The first four are all improved algorithms based on JADE. GEDGWO is a variant of GWO [8] applying the Gauss probability model to estimate the distribution. DOLTLBO is an enhanced teaching-learning-based optimization algorithm proposed in 2019. The parameters are the same as the recommended settings in the original works, as reported in Table 2.

The definitions and optimal values of CEC 2014 test suite are provided in [36]. The number of maximum function evaluations (MaxFEs) is set to D × 10000, and D is set to 30 in this paper. Each function is evaluated 51 times independently for statistics. The error value between the best obtained solution and the optimum solution is recorded as the result.

All experiments are performed on a computer with AMD R7 4800U (1.80 GHz) processor and 16 GB of RAM. The programs are implemented by MATLAB 2016B platform.

From Table 3, It is apparent that, for unimodal functions, IJADE-TSO outperforms all other comparison algorithms, and the optimal value can be stably obtained, which confirms the excellent performance of IJADE-TSO in solving unimodal functions. As for the multimodal and hybrid functions, IJADE-TSO ranked in the top four except F5 and F12. IJADE-TSO ranks in the top two except F26 in solving composition functions. It is worth noting that the latter two groups of test functions are more complex, so the superiority of IJADE-TSO in solving complex optimization problems is better demonstrated. SPS-LSHADE-EIG performs well in multimodal functions except for F14, unimodal and multimodal functions, and hybrid functions but performs poorly in composition functions except F26. LSHADE and iLSHADE are similar to SPS-LSHADE-EIG, but LSHADE shows poor performance in F19 and iLSHADE performs poorer than SPS-LSHADE-EIG for composition functions. The performances of CPI-JADE and GEDGWO are poor. DOLTLBO performs well in F23, F25, F26, and F28 of composition functions but poor in other functions.

Figure 2 shows the ranking of the algorithm in each test function. The closer to the center of the circle, the higher the ranking, and vice versa. The area enclosed by the curves shows the overall ranking of the algorithms.

4.1.1. Friedman Test

The Friedman test is used [37] to analyze the performance difference of algorithms on the 30 test functions. Table 4 shows the average ranking results of each algorithm.

The ranking value of IJDAE-TSO is 2.92, ranking the top among all algorithms. The rankings of the other six algorithms are as follows: SPS-LSHADE-EIG, iLSHADE, LSHADE, GEDGWO, CPI-JADE, and DOLTLBO. Note that, compared with the improved variant of JADE, CPI-JADE, which ranks sixth, IJDAE-TSO ranks the first, proving the effectiveness of the improvement strategies. The p value calculated by Iman-Davenport test is 4.43e-11, showing that there is a significant difference between algorithms.

4.1.2. Wilcoxon Signed-Rank Test

Reference [38] shows that it is insufficient to analyze and compare the algorithms only according to the average value. Because, in some tests, the results produced by the algorithms are slightly different, and the occurrence of these differences may be accidental, to verify whether the statistical difference between IJADE-TSO and the comparison algorithms is accidental or not, the Wilcoxon signed-rank test is used to analyze the performances of the IJADE-TSO and the other six algorithms. Table 5 lists the Wilcoxon signed-rank test [39] results of each competitor and IJADE-TSO with the significance level set to 0.05. In Table 5, “+” indicates that the optimization result of IJADE-TSO is better than the competitor, “−” denotes the opposite, and “=” represents that the optimization results are similar. “R+” denotes the magnitude with which IJADE-TSO surpasses the comparison algorithm, and “R–” indicates the opposite result. The statistical results of the test are also given in the last row of the table. According to the statistical results in Table 5, among the 30 test functions of CEC 2014 test set, IJADE-TSO shows a performance that is similar to those of SPS-LSHADE-EIG, LSHADE, and iLSHADE and precedes CPI-JADE in all test functions except F5, F7, F12, and F14.

4.1.3. Convergence Curve of Algorithms

To further illustrate the convergence performance of IJADE-TSO, the average error convergence curves in solving the test function are given in Figure 3.

For most unimodal functions, simple multimodal functions, and hybrid functions, IJADE-TSO shows better convergence accuracy and higher convergence speed compared to GEDGWO, CPI-JADE, and DLOTLBO and slower convergence speed compared to LSHADE, iLSHADE, and SPS-LSHADE-EIG. For the composite function, the convergence curve of IJADE-TSO drops slightly at the initial stage of iteration and then falls into local optimization, but it breaks through local optimization at the later stage of iteration and shows a significant decline. For composition functions, LSHADE, SPS-LSHADE-EIG, and iLSHADE converge quickly, but their accuracy is significantly worse than that of IJADE-TSO, and they are apt to fall into the local optimum. As mentioned above, IJADE-TSO possesses both strong convergence ability and the capacity to break through local optimization.

4.1.4. Box Graph of the Results

The following box graphs are drawn to analyze the distribution characteristics of the results for CEC 2014 test functions.

As shown in Figure 4, the center mark of each box represents the median value of the results of 51 runs, and the bottom and upper edges of each box represent the third and fourth quartiles, respectively. Outliers are drawn separately with “+.” IJADE-TSO has no outliers in 13 test functions (F2-F3, F5–F7, F9, F13-F14, F17, F19, F24, F26, and F29). The distribution of the IJADE-TSO solutions is more centralized, with rarely produced outliers, except for F30. Thus, it is safe to conclude that IJADE-TSO is more robust than the competing algorithms.

4.2. Parameters Identification of Photovoltaic Models

Global warming has become an increasingly serious issue threatening the survival of mankind nowadays. In response to the environmental impact of global warming, solar energy has become widely exploited as a clean energy source [40]. Solar energy has many advantages, such as being nonpolluting, renewable, cheap, and easy to obtain. Photovoltaic systems are commonly used in daily life as an important device for collecting solar energy. PV models are crucial for PV systems, as they are used to imitate the response of real PV cells, fitting measured current-voltage (I-V) data under all operating conditions. The most widely used PV systems are single-diode model (SDM), double-diode model (DDM), and PV module model [41]. It is necessary to identify the parameters in the established photovoltaic models to make the model more consistent with the measured data. In this section, the proposed IJADE-TSO algorithm is used to identify parameters of photovoltaic models.

4.2.1. Problem Statement

In [42], there are three PV models, single-diode model, double-diode model, and PV module model, which are not repeated in this paper. The PV parameter identification models referenced from [41, 42] are as follows.

In general, parameters extraction problem is usually transformed into a class of optimization problem to leverage the optimization algorithms. Similar to [42], the root mean square error (RMSE) is also used as the objective function in this study, which is defined as follows:where N is the number of measured I-V data samples and x is a vector that concludes the unknown parameters to be extracted. For different PV models, their objective functions are represented as follows.

For single-diode model,

For double-diode model,

For PV module based on the SDM,

4.2.2. Experimental Results and Analysis

To evaluate the performance of IJADE-TSO, we employ it to extract parameters of SDM, DDM, and PV module model. The I-V data of SDM and DDM are obtained from [43], which are measured on a 57 mm diameter commercial (R.T.C France) silicon solar cell under 1000 W/m2 at 33°C.

(1) Results for the Three Models. In order to further show the quality of the results, we have introduced individual absolute error (IAE).

Im is the measured current and Ic is the simulated current. IAE represents the absolute error between them.

The I-V characteristics obtained by IJADE-TSO and the individual absolute error (IAE) between the experimental data and simulated data are shown in Tables 68 and the I-V characteristics are plotted in Figure 5. In addition, Tables 68 present the individual absolute error (IAE) between the measured data and calculated data. From Table 6, all the IAE values of single-diode model are smaller than 2.5434E − 03. From Table 7, all the IAE values of double-diode model are smaller than 1.6284E − 03. From Table 8, all the IAE values of double-diode model are smaller than 4.8328E − 03. The results show that the calculated data obtained by IJADE-TSO are remarkably in accordance with the measured data over the whole voltage range, which indicates that IJADE-TSO identifies highly accurate parameters.

(2) Statistical Results and Convergence Speed. To validate the superior performance of IJADE-TSO, the comparisons are carried out with other algorithms including MLHADE [44], improved JAYA [45], JADE [27], crow search algorithm (CSA) [46], and multiple learning backtracking search algorithm (MLBSA) [47]. Among the above algorithms, MLSHAPE, MLBSA, IJAYA, and CSA are excellent algorithms that have been applied in the field of photovoltaic parameter identification in recent years. JADE is the prototype of the improved algorithm in this paper. For fair comparison, the value of function evaluations of all algorithms is set to 50000 for all problems. Meanwhile, each algorithm is tested 30 times independently. The population sizes of all compared algorithms are set to 50. Other algorithm parameter settings are the same as the references.

The statistical results are reported in Table 9. For single-diode model, IJADE-TSO has the same accuracy as MLSHADE and MLBSA. The Min RMSE searched by all algorithms have no difference except IJAYA. As for standard deviation, IJADE-TSO is better than its competitors. For double-diode model, only IJADE-TSO and MLSHADE algorithms obtain the same Min RMSE, but the Mean and SD RMSE of IJADE-TSO ranks the first. For PV module model, the minimum RMSE of PV model obtained by all algorithms are no different. Only MLBSA, IJADE-TSO, and JADE ranked the first in the Mean RMSE. The SD value obtained by IJADE-TSO ranks the first. Overall, IJADE-TSO has better search ability and robustness than the comparison algorithms and can steadily obtain the minimum RMSE value in all three PV parameter identification models.

The convergence curves among all competitors for three PV models are shown in Figure 6. Obviously, for the single-diode model, the descent curve of the IJADE-TSO algorithm is significantly higher compared to the other algorithms in the initial parts, and then it declines fastest and obtains the best results after 10,000 function evaluations. Although MLSHADE and MLBSA also obtain the best results, they are significantly slower than IJADE-TSO. For the double-diode model, IJADE-TSO achieves the best results with the least number of function evaluations. For the PV model, IJADE-TSO, MLBSA, and JADE all achieve the best results, but IJADE-TSO uses the least number of function evaluations and its iteration curve declines the fastest. Among the three PV models, the descent curve of the IJADE-TSO declines the fastest compared with other algorithms and achieves the best results.

In order to comprehensively compare the performances of all algorithms in the three PV models, we conducted a Friedman test with a significance of 0.05, and the results are shown in Table 10.

In terms of Min, Mean, and SD of RMSE, IJADE-TSO excels all other algorithms, ranking the first, which shows that IJADE-TSO has an excellent performance in terms of accuracy, robustness, and computational efficiency.

5. Conclusions

This paper proposed a novel hybrid algorithm, named IJADE-TSO. The spiral foraging search and parabolic foraging search of TSO are introduced into the mutation strategy in IJADE to improve the exploration ability and population diversity. Meanwhile, this paper integrates JADE with three strategies: CR sorting mechanism and CR repairing, top α r1 selection, and population linear reduction.

The proposed IJADE-TSO is evaluated using CEC 2014 benchmarks with dimensionality of 30 and photovoltaic parameter identification problem. Although the statistical results of the CEC 2014 benchmarks demonstrate that IJADE-TSO is slightly superior to SPS-LSHADE-EIG and iLSAHDE, IJADE-TSO is significantly better than other comparison algorithms such as MLHADE in solving photovoltaic parameter identification problem. This shows that the two mutant strategies combined contributed significantly to expanding the search scope and avoiding local optimization, and three other strategies are conducive to improve the convergence efficiency. IJADE-TSO achieves a good balance in exploration and exploitation and demonstrates great value in solving practical engineering problems.

The current spiral foraging and parabolic foraging strategies are selected with a stationary probability, which means the selected strategy may not be optimal and an adaptive selection strategy needs to be developed in the future.

Data Availability

The codes used to support the findings of this study can be downloaded at https://github.com/lailaidaren/Matlab-code-of-IJADE-TSO.git.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors acknowledge funding received from the following science foundations: the National Natural Science Foundation of China (no. 62101590), the Science Foundation of the Shanxi Province, China (2021JM-223, 2022JQ-584, and 2021JM-224), and Military Scientific Research Project (KJ2020C009002).