Research Article  Open Access
A Multiobjective Particle Swarm Optimization Algorithm Based on Multipopulation Coevolution for WeaponTarget Assignment
Abstract
Weapontarget assignment (WTA) is critical to command and decision making in modern battlefields and is a typical nondeterministic polynomial complete problem. To solve WTA problems with multiple optimization objectives, a multipopulation coevolutionbased multiobjective particle swarm optimization (MOPSO) algorithm is proposed to realize the rapid search for the globally optimal solution. The algorithm constructs a masterslave population coevolution model. Each slave population corresponds to an objective function and is used to search for noninferior solutions. The master population receives all the noninferior solutions from the slave populations, repairs the gaps between the noninferior solutions, and generates a relatively optimal Pareto optimal solution set. In addition, to accelerate the slave populations searching for noninferior solutions and master population repairing the gaps between noninferior solutions, the particle velocity update method is improved. The simulation results show that the proposed algorithm has higher computational efficiency and achieves better solutions than existing algorithms capable of providing a good solution. The method is suitable for rapidly solving multiobjective WTA (MOWTA) problems.
1. Introduction
Weapontarget assignment (WTA) [1, 2] is critical to operational command and directly affects the progression and outcome of operations. WTA is an important military problem studied by numerous military powers. The core concept of WTA is to rapidly and accurately assign weapons to targets and to obtain better solutions to satisfy operational goals (i.e., optimization objectives). These optimization objectives mainly include maximum combat effectiveness, minimum weapon consumption, minimum potential threat of target, minimum target surplus, and minimum incidental damage. According to the number of optimization objectives in the problem, WTA problems can be divided into singleobjective WTA (SOWTA) problems, with only one optimization objective, and MOWTA problems, with at least two optimization objectives. Currently, research on the SOWTA problem is more indepth, and respective optimization algorithms [3–7] can quickly obtain a better allocation scheme. For the SOWTA problem with the background of air defense intercept, Liu et al. [4] proposed a firepower unit correlation matrix and designed a hybrid optimized algorithm based on the PSO and tabu search algorithm. The simulation results show that the method is more efficient than the existing methods and can output optimization results at any time. For the SOWTA problem with the background of againstground targets, Wang et al. [5] improved the particle initialization and inertia weight selection methods of the PSO algorithm and effectively improved the optimization efficiency and allocation results of the largescale SOWTA problem. For the dynamic SOWTA problem, Cho et al. [6] constructed a static SOWTA model with several constraints and designed an improved greedy algorithm with phased optimization, which effectively improved the optimization speed of the problem, but it was still solving a static SOWTA problem; Mei et al. [7] constructed a dynamic SOWTA model based on the killing region of weapon platform and proposed a combinatorial algorithm derived from heuristic algorithm and receding horizon control. The experimental results show that the proposed model can describe the combat scene accurately and the algorithm can improve the speed of solving the dynamic SOWTA problem. But there is still a gap from the actual application. Compared with the SOWTA problem, the MOWTA problem considers multiple optimization objectives simultaneously, and the allocation scheme satisfies the multiple operational goals of the decision makers, therein being more in line with the combat and decision makers’ needs. Research on the MOWTA problem is lacking, and the optimization speed is slow. This study researches the MOWTA problem and builds a fast optimization algorithm for such problems.
The multiobjective evolutionary algorithm based on Pareto theory [8–12] is widely used for optimization in the MOWTA problem. For the MOWTA problem with the background of againstground targets, Li et al. [10] proposed a modified Pareto AC optimization (MPACO) algorithm based on newly designed operators, including a dynamic heuristic information calculation approach, improved movement probability rule, dynamic evaporation rate strategy, and global updating rule for pheromones. The improved MPACO algorithm has high computational efficiency and accuracy; however, the Pareto optimal solution set generated by the improved MPACO algorithm has a poor distribution and cannot satisfactorily cover the relatively optimal solutions of the WTA problem. For the MOWTA problem with the background of air defense intercept, Xia et al. [11] proposed an improved MOPSO algorithm called the MOPSOII algorithm. By introducing a random decomposition strategy and a cooperative framework for the variables, that algorithm decomposes a largescalevariable problem into several variable groups of moderate dimension, called subpopulations. The entire problem is then solved by cooperation between the subpopulations. Compared with the improved MPACO algorithm, the MOPSOII algorithm has a higher optimization speed; however, the solution accuracy is slightly worse. For the dynamic MOWTA problem, Chang et al. [12] divided the dynamic process into several stages and proposed an improved ABC algorithm based on ranking selection and elite guidance. The rulebased inspiration was designed for the poor initial population of the algorithm. The simulation results show that the proposed method can improve the optimization speed and the quality of the solution, but it was still solving the static MOWTA problem. Note that the application of machine learning methods to the MOWTA problem has yet to be studied. A few scholars have studied machine learning methods used to solve the SOWTA problem [13, 14]; however, the optimization performance is not satisfactory.
The PSO algorithm is easily implemented and has a high optimization speed, also obtaining satisfactory results in the optimization of the MOWTA problem. However, the PSO algorithm is prone to premature convergence and becoming trapped in locally optimal solutions, mainly caused by a rapid loss of population diversity. To maintain population diversity, researchers have introduced operations such as random mutation [15] and stochastic perturbation [16] into the PSO algorithm or adjusted its parameters (e.g., acceleration factors [17] and inertia weights [18]). Although they maintain satisfactory population diversity and improve the global search capability of the algorithm, these improvement strategies lack an effective theoretical basis for parameter adjustment and must rely on simulation experiments, which significantly weakens their adaptability. Multipopulation strategies are another important technique to improve population diversity. Zhao et al. [19] proposed a stratified multipopulation coevolution strategy, which evenly divides the population into several subpopulations, and each subpopulation evolves in a relatively independent manner. This strategy maintains the population diversity and can effectively balance the global and local search capabilities of the algorithm. Sorkhabi et al. [20] divided the population into infeasible and feasible solution populations. These two populations evolved in parallel and independently, and the particles that evolved into feasible particles in the infeasible solution population were moved into the feasible solution population. This strategy slows the decrease in population diversity and strengthens the global search capability of the algorithm. These two multipopulation cooperation strategies provide ideas for the maintenance of population diversity and have achieved good results. However, the rate of convergence is degraded, which is unfavourable to improving the operational efficiency.
To quickly optimize the MOWTA problem, considering the advantages and disadvantages of the PSO algorithm and multipopulation cooperation strategy, a multipopulation coevolution MOPSO (MPCMOPSO) algorithm is proposed. The MPCMOPSO algorithm constructs a masterslave population coevolution model. A slave population corresponds to an objective function for searching noninferior solutions and speeding up the search for noninferior solutions by randomly selecting several noninferior solutions from other populations to replace particles with lower fitness values in the population. The master population receives all the noninferior solutions from the slave populations, repairs the gaps between noninferior solutions, and generates a relatively optimal Pareto optimal solution set. Through the masterslave population coevolution model, the master and slave populations can simultaneously develop different search areas of the problem. This not only maintains the good diversity of the population but also achieves the efficient coordination and global search of multiple populations. Additionally, to accelerate the process by which the populations search for noninferior solutions, repair the gaps between the noninferior solutions, and further improve the global search speed of the MPCMOPSO algorithm, a multidirectional competition factor is introduced into the particle velocity update to guide the particles to competitively search in multiple directions. Simulation results show that, compared with the improved MPACO algorithm proposed by Li et al. [10] and the MOPSOII algorithm proposed by Xia et al. [11], the MPCMOPSO algorithm has higher computational efficiency and produces better solutions.
In modern and future battlefields, the benefit of assigning smart weapons to targets heavily relies on the assignment of sensors; thus, the sensorweapontarget assignment (SWTA) problem is derived from the WTA problem. To effectively solve the SWAT problem, Bogdanowicz et al. [21] constructed a model that seeks to rationally assign sensors and weapons to targets for maximum performance. Based on Bogdanowicz et al., Wang et al. [22] modelled the damage probability as the product of the probability of target recognition by a sensor and the damage probability of a weapon paired with the sensor. Xin et al. [23] improved the damage probability model of Wang et al. and modelled the damage probability as the product of the weapon’s probability of kill and the sensor’s probability of detection. They then proposed a marginalreturnbased constructive heuristic (MRBCH) algorithm, which achieved good optimization results. However, the MRBCH algorithm only considers one optimization objective.
The remainder of this study is organized as follows. In Section 2, the theoretical basis of the MOWTA mode, multiobjective optimization and Pareto set, and the MOPSO algorithm are introduced and analysed. In Section 3, an MPCMOPSO algorithm for solving the MOWTA problem is proposed. The results of employing different algorithms to solve the MOWTA problem are presented and analysed in Section 4. Section 5 concludes this study.
2. Theoretical Basis
2.1. Construction of the MOWTA Model
To correctly construct the MOWTA model, the following assumptions are made:
(1) The probability of weapon damage to the target is the comprehensive damage probability, considering the weapon’s penetration probability, target hit probability, target damage probability, etc.
(2) The sequence of a target strike and the maximum projection ability of a wave of strikes are not considered.
(3) When multiple weapons strike a target, the damage probability of each weapon for the target is unchanged, and each weapon does not affect each other; that is, each strike is independent of each other.
Based on the above assumptions, the MOWTA problem can be described as follows: N types of weapons are used to attack M targets. The quantities of the N types of weapons are N_{1}, N_{2}, …, N_{N}; the combat value coefficients of the M targets are V_{1}, V_{2}, …, V_{M}, respectively; and the comprehensive damage probability of the i^{th} type of weapon against the j^{th} target is p_{ij}, i=1, 2, …, N, j=1, 2, …, M. There are two constraints on the number of weapons used. Assuming that the number of the i^{th} type of weapon against the j^{th} target is m_{ij}, the constraints can be described as follows.
N_{i} is the maximum number of weapons of the type used, i.e.,The attacker can use a maximum of C_{j} weapons on the target, i.e.,The MOWTA problem has multiple optimization objectives, resulting in numerous combinations of optimization objectives. The differences in the optimization objectives and quantities of optimization objectives reflect the different operational goals of the decision makers; however, these factors do not affect the model construction technique or the optimization performance of algorithms. Therefore, this study selects two commonly used and representative optimization objectives (i.e., maximum combat effectiveness and minimum weapon consumption) to construct the model. The former seeks to maximally reduce the combat value of the target, whereas the latter seeks to destroy the target with lower cost as much as possible.
(1) Combat Effectiveness. When m_{ij} weapons of the type attack the target, the kill probability (S_{ij}) isThus, when all weapons of the N types attack the target, the kill probability (S_{j}) isLet V_{j} be the combat value coefficient of the target. Thus, the combat effectiveness (F_{1j}) of all weapons of the N types in attacking the target isThe combat effectiveness (F_{1}) of all weapons of the N types in attacking M targets isFor the targets, F_{1} is the value damage; that is, the value damage (F_{2}) of M targets isTo rapidly and clearly compare the combat effectiveness optimization results, the coefficients V_{1}, V_{2}, …, V_{M} are normalized, and their conversion satisfies the following condition:Thus, the combat effectiveness (F_{1}) of the weapons is converted toTo facilitate the use of the MOPSO algorithm in the optimization of multiple objective functions, the minimum of the objective function (f_{1}=1F_{1}) is used to represent the maximum F_{1}, i.e.,
(2) Number of Weapons Used. If m_{ij} weapons of the type are used to attack the target, then the number (m_{i}) of weapons of the type used to attack M targets isThus, the number (f_{2}) of weapons of the N types used to attack M targets isBased on the functions f_{1} and f_{2} as well as the constraints on the number of weapons used, the MOWTA model is constructed as follows:
2.2. Multiobjective Optimization and Pareto Set
In contrast to a singleobjective optimization problem, the multiobjective optimization (MOO) problem does not have a unique optimal solution; its optimal solution becomes a set that may contain an infinite number of solutions that do not dominate one another (i.e., a Pareto optimal solution set). An MOO problem can be described as follows:In (14), is a decision vector (); is an objective function vector consisting of n objective functions (); and is a constraint vector consisting of k equations () or inequalities (). Given an MOO problem with a feasible region , its set of Pareto optimal solutions is as follows:
2.3. The MOPSO Algorithm
The MOPSO algorithm is the name of the PSO algorithm used to solve MOO problems. The MOPSO and PSO algorithms update the particle velocity and location in the same way. The PSO algorithm updates the velocity and location of a particle through two “guides”: the optimal solution that the particle has found (i.e., the individual guide ) and the optimal solution that the entire population has found thus far (i.e., the global guide ). In the basic variant of the PSO algorithm, the velocity and location updates satisfy the following conditions:In (16) and (17), in the ddimensional direction of the search space, is the current velocity of the particle, is the current location of the particle, is the updated velocity of the particle, is the updated location of the particle, is the optimal solution that the particle has found, is the optimal solution that the entire population has found thus far, w is the inertia weight, c_{1} and c_{2} are acceleration coefficients, and r_{1} and r_{2} are random numbers within the interval .
The MOPSO and PSO algorithms also differ in several respects: first, the update and selection methods of the individual and global guides are different. In addition, the solution obtained by the MOPSO algorithm is a Pareto optimal solution set. Moreover, the MOPSO algorithm requires a storage set for storing the noninferior solutions found by the population to be established. Numerous researchers have conducted indepth studies on the updating and selection of individual and global guides, the establishment of a storage set, and the parameter settings in the MOPSO algorithm and have obtained some results [24–26]. On largescale problems, the MOPSO algorithm is still prone to premature convergence and faces difficulties in efficiently searching for the global optimal solution.
3. The MPCMOPSO Algorithm
To maintain high population diversity and perform a rapid and global search for the Pareto optimal solutions of the MOWTA problem, the MPCMOPSO algorithm constructs a masterslave population coevolution model, which allows the algorithm to simultaneously develop various search regions for the problem and perform an efficient multipopulation cooperative global search. Additionally, to accelerate the process in which the populations search for noninferior solutions, repair the gaps between the noninferior solutions, and further improve the global search speed of the MPCMOPSO algorithm, a multidirectional competition factor is introduced into the particle velocity update to guide the particles to competitively search in multiple directions.
3.1. Integer Coding Based on the Attacking Target
In solving the MOWTA problem, the PSO algorithm first needs to encode its particles. The coding method should not only represent the solution of the problem and attempt to satisfy the constraints set in the model but also reduce the length of the code as much as possible. For this purpose, an integer coding method based on the attacking target was designed. The coding structure is shown in Figure 1. In Figure 1, the length of the code is N×M; that is, the dimension is D= N×M, where, P_{1}, P_{2}, …, and P_{N×1} represent the allocation of all types of weapons for target 1; P_{N+1}, P_{N+2}, …, and P_{N×2} represent the allocation of all types of weapons for target 2; and P_{N×(M1)+1}, P_{N×(M2)+1}, …, and P_{N×M} represent the allocation of all types of weapons for target M. This coding method does not distinguish the order of attack on the same target for the same type of weapon, and its weapon usage is represented by only one code, which greatly reduces the length of the code.
3.2. Particles’ Multidirectional Competitive Search for Optimal Solutions
A principal factor that leads the PSO and MOPSO algorithm towards premature convergence and their low global search efficiency is described as follows. The search space contains multiple local optimal regions. The populations are guided by a rule that they tend to evolve in the direction of the optimal particle in the populations. Under this rule, the particles in the populations rapidly exchange information and tend to move in a single direction. Consequently, the populations rapidly converge to a local optimal region, and the population diversity rapidly decreases, resulting in premature convergence. If the region does not contain the globally optimal solution, the algorithm becomes trapped in local optima. The algorithm requires a substantial amount of time to jump out of local optimal solutions, resulting in low global search efficiency.
If the particles can also evolve in the directions of the optimal particles in their locally optimal regions while tending to evolve in the direction of the optimal particle, it is possible to prevent the populations from rapidly accumulating in a certain locally optimal region, thereby maintaining satisfactory population diversity, effectively preventing the algorithm from becoming trapped in locally optimal solutions and increasing the algorithm’s global search speed. Hence, the particle update method is improved. A multidirectional competition factor is introduced into the particle velocity update to guide the particles to competitively search in multiple directions. The particles may still evolve in the directions of the optimal particles in their locally optimal regions even while evolving in the direction of the optimal solution. The improvement of the particle update method includes two parts: determination of the local guides and introduction of a multidirectional competition factor.
3.2.1. Determination of Local Guides
It is difficult to obtain the locally optimal regions in the search space by mathematical operations. Instead, clustering techniques can rapidly obtain more accurate locally optimal regions. Zhang et al. [27] proposed a clustering technique based on Pareto domination that can rapidly obtain accurate locally optimal regions. In this section, the Paretodominationbased clustering technique is used to obtain locally optimal regions. “Local guides” are particles that tend to evolve in locally optimal regions and are the optimal particles in the locally optimal regions. In an MOO problem, the optimal particles are often not unique, and there may be multiple optimal particles that do not dominate one another. As a result, multiple particles may become local guides. Particle diversity can easily suffer when particles evolve in a single direction. To prevent the rapid loss of particle diversity in a locally optimal region caused by the selection of a fixed particle as the guide for all particles within a region, one particle is randomly selected from multiple optimal particles as the local guide for each particle.
3.2.2. Introduction of a Multidirectional Competition Factor into the Velocity Update
To prevent premature convergence and low global search efficiency caused by the tendency of particles to evolve in a single direction (the direction of the optimal particle in a population), a multidirectional competition factor is introduced into the particle velocity update to allow particles to evolve in the directions of the optimal particles within their locally optimal regions. This technique can realize a multidirectional competitive search by the particles, effectively preventing the algorithm from experiencing premature convergence or becoming trapped in local optima and thus improving the global search efficiency. If is the local guide for particle i, then after introducing a multidirectional competition factor, the velocity update shown in (16) is improved as follows:In (18), c_{3} is an acceleration coefficient (), r_{3} is a random number in the interval , and l_{id} is the component of in the dimension of the search space.
3.3. MasterSlave Population Coevolution Model
The MPCMOPSO algorithm constructs a masterslave population coevolution model, which includes a master population and multiple slave populations. The number of slave populations is determined by the number of objective functions. One slave population corresponds to one objective function. By implementing the PSO algorithm containing a multidirectional competition factor, the noninferior solutions of the objective function can be searched. The master population receives all noninferior solutions from the slave populations, repairs the gaps between the noninferior solutions by implementing the MOPSO algorithm containing a multidirectional competition factor, and generates a relatively optimal Pareto optimal solution set. Figure 2 shows the masterslave population coevolution model.
3.3.1. Slave Population Evolution Model
For an MOO problem with as the objective functions, the slave population evolution model constructs n populations and encodes the particles in the populations using the encoding method proposed in Section 3.1. Each slave population corresponds to an objective function (i=1, 2, …, n) and searches for an optimal solution to the corresponding objective function by using the PSO algorithm containing a multidirectional competition factor. The particles in a slave population are the decision vectors of not only the singleobjective function but also the multiobjective function and will generate the noninferior solutions needed by the master population during the iterative optimization process. Therefore, after each iterative search of a slave population, it is necessary to determine the noninferior solutions in the population and transfer them to the master population.
The slave populations operate in parallel and independently. Since the rates of evolution of the slave populations are different, some slave populations will inevitably have converged while others will have not yet converged. For slave populations that have converged, new noninferior solutions will not be generated due to the loss of population diversity. To avoid this deficiency and increase the rate of convergence of the slave populations that have not converged, a cooperative search strategy for the slave populations is proposed as follows. After each iteration of a slave population, several noninferior solutions are randomly selected from other slave populations to replace the particles with lower fitness values. The cooperative search strategy for the slave populations not only avoids the situation whereby some slave populations cannot generate new noninferior solutions due to rapid convergence but also enables the slave populations to assist one another and jointly search for noninferior solutions of singleobjective functions.
3.3.2. Master Population Evolution Model
After each iteration, the master population receives all the noninferior solutions from the slave populations. Then, the population searches for more and better noninferior solutions by the MOPSO algorithm proposed, thus repairing the gap between the noninferior solutions and generating a good Pareto optimal solution set.
Since the master population continuously receives noninferior solutions from the slave populations, its scale and particle density increase, resulting in a decrease in the optimization speed. Densely distributed particles not only are unfavourable to the maintenance of population diversity but also significantly reduce the search efficiency of the populations. Therefore, it is necessary to eliminate some densely distributed particles.
The more densely the particles are distributed, the closer the particles are to each other, and the “crowding distance” is used to measure the density of the particles. The crowding distance of a particle is the sum of its distances from its two nearest particles. If particles j and k are the particles closest to particle i, the crowding distance (D_{i}) of particle i isBased on (19), the larger the crowding distance is, the more scattered the particles and the lower the particle density are. Additionally, the crowding distances of all particles in the master population can be obtained by using (19). If the population scale is limited to S, then the master population only retains the top S particles with the largest crowding distances.
3.4. The MPCMOPSO Algorithm Flow
Figure 3 shows the flowchart of the MPCMOPSO algorithm. The steps of the MPCMOPSO algorithm are as follows:
(1) Initialize the slave populations. The number of slave populations is determined by the number of objective functions. First, encode the particles of the population following the method proposed in Section 3.1. Then, initialize the particles.
(2) Search for noninferior solutions. The slave population evolution model is used to search for noninferior solutions to the singleobjective functions and move them to the master population.
(3) Generate a Pareto optimal solution set. The master population receives all the noninferior solutions from the slave populations and uses the master population evolution model to repair the gaps between the noninferior solutions. Then, a good Pareto optimal solution set is generated.
(4) Determine whether the termination condition is satisfied. If the termination condition is met, the algorithm stops searching; otherwise, the algorithm returns to Step (2).
4. Simulation and Analysis of the MOWTA Problem
4.1. MOWTA Problem Cases
The case background is as follows. There are ten different types of weapons available, and twelve targets will be attacked. The number of each type of weapon, the combat value coefficient of each target, and the maximum number of weapons that can be used on each target are shown in Table 1. The comprehensive damage probability of each type of weapon for the target is shown in Table 2. Objectives with the maximum combat effectiveness and the minimum weapon consumption must be optimized.


Based on the above background, five different cases were established. The weapon types and targets in each case are shown in Table 3.

4.2. Simulation and Analysis of the MOWTA Problem
The improved MPACO algorithm proposed by Li et al. [10], the MOPSOII algorithm proposed by Xia et al. [11], and the MPCMOPSO algorithm proposed in this study were implemented to solve the MOWTA problem described above using the same software and hardware (64bit Windows 7, Intel® Xeon® E31240 (Version 3) CPU @3.4 GHz, and 16 GB of RAM). The time (unit: s) required by each of the three algorithms to perform the optimization and the optimal assignment determined by each of the three algorithms were recorded.
4.2.1. Optimization Speeds of the Algorithms
Table 4 summarizes the time consumption statistics of the three algorithms after 20 simulations in 5 cases. and are the shortest and longest time of 20 simulations, respectively. and are the average time and the standard deviation of the time spent on 20 simulations, respectively.

A comparison of , , and for the optimizations by the three algorithms shows that the , , and of the MPCMOPSO algorithm are much lower than those of the other two algorithms, suggesting that the MPCMOPSO algorithm can significantly increase the optimization speed on MOWTA problems. Additionally, as the scale of the problem increases, of the MPCMOPSO algorithm does not increase significantly compared to that of the other two algorithms, indicating that the MPCMOPSO algorithm is suitable for solving largescale MOWTA problems. A comparison of the of the three algorithms shows that the of the MPCMOPSO algorithm is much lower than that of the other two algorithms, indicating that the MPCMOPSO algorithm has good optimization stability.
4.2.2. Results Obtained Using the Algorithms
Figures 4–8 visualize the results of the Pareto optimal solution set (Pareto front, i.e., optimal numerical relationships between the maximum combat effectiveness and the minimum weapons consumption) obtained by the three algorithms in Cases 1 to 5, respectively.
Figures 4–8 show that, compared with the improved MPACO algorithm and the MOPSOII algorithm, the Pareto front optimized by the MPCMOPSO algorithm is closer to the true Pareto front. The comparison demonstrates that the MPCMOPSO algorithm is accurate. Figures 4–8 also show that, with increasing problem scale, the MPCMOPSO algorithm can still obtain good results, which shows that the MPCMOPSO algorithm can be used to solve largescale MOWTA problems.
For the targets, if additional weapons are put into the strike sequence, the targets will be assigned more suitable weapons. For the weapons, if additional targets are added to the strike mission, the weapons will choose more suitable targets to attack. Both of the above cases will increase the value damage of targets. To further evaluate the effectiveness, whether the MPCMOPSO algorithm can assign more suitable weapons to the targets for the above two cases is determined. Tables 5 and 6 record the value damage of the targets under the allocation scheme obtained by the MPCMOPSO algorithm in Case 1 to Case 3 and Case 3 to Case 5, respectively. The “N.A.” in Table 5 indicates that there is no calculated value due to an insufficient number of weapons. The “N.A.” in Table 6 indicates that there is no calculated value due to the weapon usage exceeding the maximum weapon usage to targets.


Table 5 shows that, under the same weapon consumption, the value damage of targets in Case 3 is higher than that in Case 1 and Case 2, and the value damage of targets in Case 2 is higher than that in Case 1. This shows that if additional weapons are available, the MPCMOPSO algorithm can assign weapons with a higher damage probability to the target and obtain a better allocation scheme. Table 6 shows that, under the same weapon consumption, the value damage of targets in Case 5 is higher than that in Case 3 and Case 4, and the value damage of targets in Case 4 is higher than that in Case 3. This shows that if additional targets are added to the strike mission, the MPCMOPSO algorithm can assign targets with greater operational value to the weapons.
The comparative analysis in Section 4.2 demonstrates that the MPCMOPSO algorithm can significantly improve the optimization speed on the MOWTA problem and generate relatively optimal assignments. Furthermore, the MPCMOPSO algorithm can be used to solve largescale MOWTA problems.
5. Conclusion
To quickly search for the optimal solution of the MOWTA problem, the MPCMOPSO algorithm is proposed. The MPCMOPSO algorithm constructs a masterslave population coevolution model and introduces a multidirectional competition factor into the particle velocity update. The simulation results show that the MPCMOPSO algorithm has higher computational efficiency and obtains better solutions than the improved MPACO algorithm or the MOPSOII algorithm. The MPCMOPSO algorithm can significantly increase the optimization speed of the MOWTA problem and generate relatively optimal assignments. Furthermore, the MPCMOSPO algorithm can be used for rapidly solving largescale MOWTA problems.
The MPCMOPSO algorithm can also solve the SWTA problem. Compared with the MRBCH algorithm, which can solve the SWTA problem well, the MPCMOPSO algorithm performs better in solving the SWTA problem with multiple optimization objectives.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was jointly supported by the National Natural Science Foundation for Young Scientists of China (Grants Nos. 61403397 and 61202332) and the Natural Science Foundation of Shanxi Province, China (Grant No. 2015JM6313).
References
 M. F. Hocaoğlu, “Weapon target assignment optimization for land based multiair defense systems: a goal programming approach,” Computers & Industrial Engineering, 2019. View at: Google Scholar
 X. Li, D. Zhou, Q. Pan et al., “Weapontarget assignment problem by multiobjective evolutionary algorithm based on decomposition,” Complexity, vol. 2018, Article ID 8623051, 19 pages, 2018. View at: Publisher Site  Google Scholar
 X. Hu, P. Luo, X. Zhang et al., “Improved ant colony optimization for weapontarget assignment,” Mathematical Problems in Engineering, vol. 2018, Article ID 6481635, 14 pages, 2018. View at: Publisher Site  Google Scholar
 Z. Liu, Z. Shi, L. Wu et al., “Solving cooperative antimissile weapontarget assignment problems using hybrid algorithms based on particle swarm and tabu search,” International Conference on Computer Science and Application Engineering, pp. 898–906, 2017. View at: Google Scholar
 W. Shunhong, Q. Yang, R. Wang et al., “Particle Swarm Optimization based weapontarget assignment for attacking ground targets,” Electronics Optics & Control, vol. 24, no. 3, pp. 36–40, 2017. View at: Google Scholar
 D.H. Cho and H.L. Choi, “Greedy maximization for assetbased weapontarget assignment with timedependent rewards,” Cooperative Control of MultiAgent Systems: Theory and Applications, pp. 115–139, 2017. View at: Google Scholar  MathSciNet
 Z. Mei, Z. Peng, and X. Zhang, “Optimal dynamic weapontarget assignment based on receding horizon control heuristic,” in Proceedings of the 13th IEEE International Conference on Control and Automation, ICCA 2017, pp. 876–881, Macedonia, Balkans, July 2017. View at: Google Scholar
 C. Feng, P. Yao, Y. Jing et al., “Weapontarget assignment based on improved quantuminspired immune clonal multiobjective optimization algorithm,” Computer Engineering and Science, vol. 39, no. 12, pp. 2314–2319, 2017. View at: Google Scholar
 Y. Zhang, Z. Qiao, and J. Jng, “Research on weapontarget allocation based on genetic algorithm,” Fuzzy Systems and Data Mining, pp. 260–266, 2016. View at: Google Scholar
 Y. Li, Y. Kou, Z. Li et al., “A modified pareto ant colony optimization approach to solve biobjective weapontarget assignment problem,” International Journal of Aerospace Engineering, vol. 2017, Article ID 1746124, 14 pages, 2017. View at: Google Scholar
 X. Xinxue, L. Yangtao, W. Fan et al., “Weapontarget assignment with an improved multiobjective particle swarm optimization algorithm,” Acta Armamentarii, vol. 37, no. 11, pp. 2085–2093, 2016. View at: Google Scholar
 T. Chang, D. Kong, N. Hao et al., “Solving the dynamic weapon target assignment problem by an improved artificial bee colony algorithm with heuristic factor initialization,” Applied Soft Computing, pp. 845–863, 2018. View at: Google Scholar
 H. Mouton, J. Roodt, and H. Le Roux, “Applying reinforcement learning to the weapon assignment problem in air defence,” Scientia Militaria: South African Journal of Military Studies, vol. 39, no. 2, pp. 123–140, 2011. View at: Google Scholar
 Y. Wang, J. Sun, M. Xiao et al., “Research of dynamic weapontarget assignment problem based on type2 interval fuzzy knearest neighbors classifier,” Systems Engineering and Electronics, vol. 38, no. 6, pp. 1314–1319, 2016. View at: Google Scholar
 X. Zhang, X. Wang, Q. Kang et al., “Differential mutation and novel social learning particle swarm optimization algorithm,” Information Sciences, vol. 480, pp. 109–129, 2019. View at: Google Scholar
 Z. Ouyang, Y. Liu, S. Ruan et al., “An improved particle swarm optimization algorithm for reliabilityredundancy allocation problem with mixed redundancy strategy and heterogeneous components,” Reliability Engineering & System Safety, vol. 181, pp. 62–74, 2019. View at: Google Scholar
 F. Zhou and S. Chen, “DVHop node localization algorithm based on improved particle swarm optimization,” Communications, Signal Processing, and Systems, vol. 423, pp. 541–550, 2018. View at: Google Scholar
 A. Nickabadi, R. Safabakhsh, and M. Ebadzadeh, “Inertia weight control strategies for PSO algorithms,” in Swarm Intelligence, vol. 1 of Principles, Current Algorithms and Methods, pp. 169–198, 2018. View at: Google Scholar
 N. Zhao and J. Zhao, “Selfadaptive particle swarm algorithm based on stratified multipopulation,” Journal of Lanzhou University of Technology, vol. 43, no. 4, pp. 108–111, 2017. View at: Google Scholar
 S. A. Ebrahim, A. M. Deljavan, and K. A. Reza, “Duality evolution: an efficient approach to constraint handling in multiobjective particle swarm optimization,” Soft Computing, vol. 21, no. 24, pp. 7251–7267, 2017. View at: Google Scholar
 Z. Bogdanowicz and N. Coleman, “Sensortarget and weapontarget pairings based on auction algorithm,” in Proceedings of the 11th WSEAS International Conference on Applied Mathematics, pp. 92–96, 2007. View at: Google Scholar
 J. Wang and C. Chen, “Sensorweapon joint management based on improved genetic algorithm,” in Proceedings of the 34th Chinese Control Conference, pp. 2738–2742, 2015. View at: Google Scholar
 B. Xin, Y. Wang, and J. Chen, “An efficient marginalreturnbased constructive heuristic to solve the sensorweapontarget assignment problem,” IEEE Transactions on Systems, Man, and Cybernetics: Systems, pp. 1–12, 2018. View at: Google Scholar
 A. Monir, A. AlMuta'a Ebtsam, and A. Sanabani Maher, “Integrated MOPSO algorithms for task scheduling in cloud computing,” Journal of Intelligent & Fuzzy Systems, pp. 1–14, 2018. View at: Google Scholar
 W. Zhao, Z. Luan, and C. Wang, “Parameter optimization design of vehicle EHHPS system based on an improved MOPSO algorithm,” Advances in Engineering Software, vol. 123, pp. 51–61, 2018. View at: Publisher Site  Google Scholar
 M. Najimi and S. H. R. Pasandideh, “Modelling and solving a biobjective single period problem with incremental and all unit discount within stochastic constraints: NSGAII and MOPSO,” International Journal of Services and Operations Management, vol. 30, no. 4, pp. 520–541, 2018. View at: Publisher Site  Google Scholar
 Y. Zhang and G. Gong, “Advanced multiobjective particle swarm optimization theory and its application,” Beijing City: Science Press, pp. 7778, 2016. View at: Google Scholar
Copyright
Copyright © 2019 Guangyuan Fu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.