# Research on Multiaircraft Cooperative Suppression Interference Array Based on an Improved Multiobjective Particle Swarm Optimization Algorithm

Authors: Huan Zhang; Rennong Yang; Changyue Sun; Haiyan Han
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

For the problem of multiaircraft cooperative suppression interference array (MACSIA) against the enemy air defense radar network in electronic warfare mission planning, firstly, the concept of route planning security zone is proposed and the solution to get the minimum width of security zone based on mathematical morphology is put forward. Secondly, the minimum width of security zone and the sum of the distance between each jamming aircraft and the center of radar network are regarded as objective function, and the multiobjective optimization model of MACSIA is built, and then an improved multiobjective particle swarm optimization algorithm is used to solve the model. The decomposition mechanism is adopted and the proportional distribution is used to maintain diversity of the new found nondominated solutions. Finally, the Pareto optimal solutions are analyzed by simulation, and the optimal MACSIA schemes of each jamming aircraft suppression against the enemy air defense radar network are obtained and verify that the built multiobjective optimization model is corrected. It also shows that the improved multiobjective particle swarm optimization algorithm for solving the problem of MACSIA is feasible and effective.


---

## Body

## 1. Introduction

Under the condition of information warfare, the enemy air defense radar network often has a strong ability to resist “the four” [1, 2]. It increases the difficulty of combat aircraft penetration safely, so combat aircraft must have the aid of suppression interference against the enemy air defense radar network by the electronic support jamming aircraft. In this way it can provide a safe route planning space for subsequent route planning. Since the enemy air defense radar network is often deployed by a number of air defense radars, if only relying on a single electronic jamming aircraft to interfere with the enemy radar network, it is difficult to achieve the desired suppression effect due to the limited interference resources on a single jamming aircraft. Therefore, a solution of multiple electronic jamming aircraft cooperating with each other to interfere with the air defense radar network must be used. In order to make more reasonable distribution of interference resources in the process of cooperative interference, a reasonable array pattern should be used to determine the position of each jamming aircraft. This is the optimal problem of MACSIA in electronic warfare mission planning. For this problem, the domestic and foreign research are mainly focused on the operational efficiency and suppression interference effect of combat aircraft, and so on. But the study on the MACSIA is rare. In [3], Shi et al. analyzed the influence of electronic jamming on the path planning of combat aircraft. In [4], Wang et al. proposed a multiconstraint condition genetic algorithm to optimize the deployment of the enemy radar network. The calculation method of operational efficiency of combat aircraft against ground warning radar is studied under the condition of stand-off jamming [5]. Ruan et al. [6] have explored the influence of all kinds of factors on the suppression interference effect by taking the minimum interference distance as evaluation criterion. In a work by Tang et al. [7], the evaluation model of active suppression interference effect is built and the timing of interference is analyzed. In [8], Chen et al. have built a single objective optimal model of MACSIA, but the MACSIA is not considered as a multiobjective optimization problem (MOP) to solve. All the objective functions are aggregated into a single objective function by the weighted sum method, so a single objective optimization model for this problem is built. Different Pareto optimal solutions can be obtained by selecting different weight combinations. But the drawbacks of this method are obvious. Obviously the selection of weights is related to the relative importance of each objective function. If the user does not have sufficient prior knowledge of the problem, it is difficult to find the Pareto optimal solution that satisfies the decision maker. However, if the problem is regarded as a MOP to solve, it can avoid these drawbacks effectively.In this paper, for the optimal problem of solving MACSIA, firstly, the calculation method of the enemy air defense radar detection range under the terrain masking condition is given, and the active interference model of electronic warfare is built. Secondly, the concept of route planning security zone is proposed and the mathematical morphology method is used to calculate the minimum width of route planning security zone. Thirdly, fully considering the characteristics of MOP for MACSIA, the minimum width of route planning security zone and the sum of the distance between each jamming aircraft and the center of the enemy air defense radar network are taken as objective functions, so the multiobjective optimization model of MACSIA is built. Finally, multiobjective particle swarm optimization (MOPSO) [9–17] has been widely used in multiobjective optimization based on its simple, fast convergence, easy to achieve in engineering [18–20], and so on. However, due to the disadvantages of particle swarm optimization algorithm [21], the distribution of nondominated solutions generated by MOPSO along the Pareto front is not very uniform and the computation time of MOPSO is not very fast enough. In order to solve these problems, the decomposition mechanism and the proportional distribution mechanism are introduced into the MOPSO algorithm. Therefore, an improved MOPSO algorithm is proposed which is used to solve the multiobjective optimization model of MACSIA. The optimal MACSIA schemes for the limit of the minimum width of route planning security zone and ensuring the safety of the jamming aircrafts are obtained by the simulation examples.


## 2. Model Building

### 2.1. Radar Detection Model under the Terrain Masking Condition

#### 2.1.1. Digital Elevation Model (DEM)

Digital elevation model (DEM) refers to a method of storing terrestrial elevation information in the form of specific data. There are two common representations of digital elevation models, which are the grid structure and the contour map. In this paper, the grid structure is used.


#### 2.1.2. Radar Detection Model

Radar detection model is generally characterized by radar equation which is defined as(1)R=PtG2σλ24π3S/NFnkToBnL1/4.In this equation, R is the radar detection range (m), Pt is the radar emitting peak power (W), G is the radar antenna gain, σ is the radar cross section (m2), λ is the radar wavelength (m), S/N is the signal to noise ratio of radar receiver, Fn is the noise coefficient, k is the Boltzmann constant (1.38×10-23J/K), To is the radar working temperature (K), Bn is the noise bandwidth (Hz), and L is the system loss.


#### 2.1.3. Calculating the Radar Detection Range under Terrain Masking Condition

The calculation diagram of radar detection range under terrain masking condition is shown in Figure1.Figure 1
The calculation diagram of radar detection range under terrain masking condition.The height of the radar observation point ishr. The height of the target point is hf and the target point represents the combat aircraft. The polar coordinate grid P(α,3) is the furthest boundary point of radar detection range in α direction and h(α,3) is the elevation value at this point.In Figure1, not only the height angle εP(α,k) from the polar coordinate grid P(α,k) to radar observation point in α direction can be calculated, but also the height angle εo from the target point to radar observation point can be got, which are shown as(2)εPα,k=arctanhα,k-hrPα,k-Pα,0,εo=arctanhf-hrPα,n-Pα,0. Here, the values of εP(α,k)(k=0,1,2,…,k) are calculated in α direction. If there is a grid P(α,k) satisfying εP(α,k) is greater than εo; then the target point is within the shadow of the radar terrain and the radar could not detect the target.


### 2.2. Active Interference Model of Electronic Warfare

The interference power expression of active electronic warfare equipment is defined as(3)Prj=PjGjGφrjBr4π2Rj2LjBj. Here, Pj is the jammer transmission power (W), Gj is the jammer antenna gain, G(φ) is the gain of the radar antenna in the jammer interference direction, φ is the angle between the direction of the main lobe of the enemy radar antenna and jammer, rj is the jammer polarization loss, Bj is the jammer interference signal bandwidth (Hz), Br is the signal bandwidth of the enemy radar receiver (Hz), Rj is the distance between the enemy radar and the jammer (m), and Lj is the interference signal loss.The power expression of the target echo signal is defined as(4)Prt=PtG2σλ24π3R4L.The meaning of each parameter in the above formula is shown in formula (1), where the constant Prt is equivalent to the parameter S in formula (1).The expression ofG(φ) is defined as(5)Gφ=G,0≤φ≤φ0.5,K2φ0.5φ2G,φ0.5≤φ≤90°,K2φ0.590°2G,90°≤φ≤180°.In this equation, φ0.5 is the lobe width of the enemy radar antenna at the half power point, K is a constant between 0.04 and 0.1 in general, and G is the enemy radar antenna gain.The detection range expression of the enemy radar can be obtained from formulas (1), (3), and (4) in the case of interference with a single jammer, which is shown as(6)R=PtG2σλ24π3LFnkToBnKj+4πLPjGjGφrjBr/Bj2LjBj1/4. Here, Kj is the minimum suppression coefficient of the enemy radar and the remaining parameters are the same as in formulas (1), (3), and (4).When considering suppression interference with multiple jammers against the enemy air defense radar network, the total interference power received is the sum of the interference powers of each jammer, which is defined as(7)Prj_total=∑i=1nPjiGjiGφirjiBr4π2Rji2LjiBji.The detection range expression of the enemy air defense radar after cooperative suppression by multiple jammers can be calculated, which is shown as(8)R=PtG2σλ24π3LFnkToBnKj+∑i=1n4πLPjiGjiGφirjiBr/Bji2LjiBji1/4.


### 2.3. Calculation Model of Route Planning Security Zone

#### 2.3.1. Concept of Route Planning Security Zone

Electronic jamming aircraft implements active suppression interference against the enemy air defense radar network, which aims to suppress the detection range of the enemy air defense radar and expand the space scope of route planning security. So it can provide more safe and reliable planning space for subsequent route planning. The definition of route planning security zone is given here. Route planning security zone is a certain width, height range for combat aircraft flight track space which constitutes the search space for the subsequent optimal route planning of combat aircraft. For the sake of convenience, this paper mainly calculates the width of route planning security zone for combat aircraft in a certain altitude. This width refers to the minimum width within the entire security zone.


#### 2.3.2. Calculating the Width of Route Planning Security Zone Based on Mathematical Morphology

Since the detection range of the enemy air defense radar is influenced by the terrain, the detection boundary is irregular, so the model building is more difficult. If the traditional geometric method is used to calculate the width of security zone, the calculation process will be very complex and not conducive to engineering practice. Therefore, this paper will proceed from the mathematical morphology point of view so as to calculate the width of route planning security zone. Here the brief introduction to mathematical morphology will be given.The core idea of mathematical morphology is to use a probe structure element to detect a given image and gets the information about the image, so the image can be analyzed and processed [22]. Mathematical morphology involves the basic operation of a closed operation, open operation, corrosion, expansion, and so on. The effect of the corrosion and expansion is opposite. The corrosion can cause a given area to shrink inward at the same time while the expansion allows the region to expand around. Both the closed operation and open operation are complexes of corrosion and expansion. The closed operation expands the image firstly and then corrodes. As a result, a narrow fracture on a given image can be filled and the holes in the image can be removed, and so on. On the contrary, the open operation corrodes the image firstly and then expands, which usually has a smoothing effect on a given image contour and eliminates small burrs on the contour.This paper will be followed by the principle of image expansion and corrosion. Firstly, the image of radar detection range under terrain masking and electronic interference is binarized [22]. Secondly, we will do open operation on the binary image and then expand it. Thirdly, continuity check is performed while the image is being expanded and then the number of the graphic elements in the image is obtained. Finally, it is judged whether the route planning security zone satisfying the minimum width restriction is formed, and then the minimum width of security zone is calculated. The specific calculation process is shown in Figure 2.Figure 2
The calculation process of route planning security zone.


### 2.4. Multiobjective Optimization Model of MACSIA

A typical maximization MOP can be defined as follows [23]:(9)maxf1x,f2x,…,fnxs.t.lb≤x≤ubAeq∗x=beqA∗x≤b.In this equation, n is the dimension of objective space, x∈[x1,x2,…,xDim]T is decision variable, Dim represents the number of decision variables, ub and lb are the upper and lower bounds of the decision variable x, respectively, and A∗x≤b and Aeq∗x=beq are linear inequality constraints and equality constraints for decision variables, respectively.The ultimate goal of optimizing the MACSIA is to achieve a reasonable distribution of jamming aircraft positions and suppress the best against the enemy air defense radar, and the jamming aircraft themselves are not threatened. Therefore, this problem is a MOP. The expressions of single objective are defined as(10)f1=∑i=1Nradarxi-xtc2+yi-ytc2+zi-ztc21/2(11)f2=Widthsafe.Here, the objective function f1 represents the sum of the distances between our jamming aircrafts and the center of the enemy radar network, and the physical meaning is our own safety of jamming aircrafts. The objective function f2 represents the product of the quantity and width of route planning security zone, and the physical meaning is the interference suppression effect of our jamming aircrafts against the enemy air defense radar network. Nradar is the number of the enemy radars. (xi,yi,zi) represents the coordinates of the ith jamming aircraft, (xtc,ytc,ztc) is the center coordinates of the enemy radar network, and Widthsafe is the minimum width of the security zone.The following constraints and interference principles should be considered when the MACSIA is implemented. Firstly, the position of each jamming aircraft must be outside the maximum detection range of the enemy radar. Secondly, the height of each jamming aircraft should be within the given safe height range. Thirdly, the interference principle is traditional principle of many-to-one and one-to-one. The constraints are mathematically expressed as follows:(12)s.t.sqrtxi-xtj2+yi-ytj2+zi-ztj2>Rji=1,2,…,NECM,j=1,2,…,Nradarzmin≤zi≤zmax.Here, zmax and zmin are the upper and lower bounds of the given height range, respectively. The parameter (xi,yi,zi) can be seen in (10), NECM is the quantity of the jamming aircrafts, (xtj,ytj,ztj) represents the coordinates of the ith enemy radar, Nradar is the quantity of the enemy radar, and Rj is the maximum detection range of the ith enemy radar.The multiobjective optimization model of MACSIA with constraints is built by formulas (10), (11), and (12).


## 2.1. Radar Detection Model under the Terrain Masking Condition

### 2.1.1. Digital Elevation Model (DEM)

Digital elevation model (DEM) refers to a method of storing terrestrial elevation information in the form of specific data. There are two common representations of digital elevation models, which are the grid structure and the contour map. In this paper, the grid structure is used.


### 2.1.2. Radar Detection Model

Radar detection model is generally characterized by radar equation which is defined as(1)R=PtG2σλ24π3S/NFnkToBnL1/4.In this equation, R is the radar detection range (m), Pt is the radar emitting peak power (W), G is the radar antenna gain, σ is the radar cross section (m2), λ is the radar wavelength (m), S/N is the signal to noise ratio of radar receiver, Fn is the noise coefficient, k is the Boltzmann constant (1.38×10-23J/K), To is the radar working temperature (K), Bn is the noise bandwidth (Hz), and L is the system loss.


### 2.1.3. Calculating the Radar Detection Range under Terrain Masking Condition

The calculation diagram of radar detection range under terrain masking condition is shown in Figure1.Figure 1
The calculation diagram of radar detection range under terrain masking condition.The height of the radar observation point ishr. The height of the target point is hf and the target point represents the combat aircraft. The polar coordinate grid P(α,3) is the furthest boundary point of radar detection range in α direction and h(α,3) is the elevation value at this point.In Figure1, not only the height angle εP(α,k) from the polar coordinate grid P(α,k) to radar observation point in α direction can be calculated, but also the height angle εo from the target point to radar observation point can be got, which are shown as(2)εPα,k=arctanhα,k-hrPα,k-Pα,0,εo=arctanhf-hrPα,n-Pα,0. Here, the values of εP(α,k)(k=0,1,2,…,k) are calculated in α direction. If there is a grid P(α,k) satisfying εP(α,k) is greater than εo; then the target point is within the shadow of the radar terrain and the radar could not detect the target.


## 2.1.1. Digital Elevation Model (DEM)

Digital elevation model (DEM) refers to a method of storing terrestrial elevation information in the form of specific data. There are two common representations of digital elevation models, which are the grid structure and the contour map. In this paper, the grid structure is used.


## 2.1.2. Radar Detection Model

Radar detection model is generally characterized by radar equation which is defined as(1)R=PtG2σλ24π3S/NFnkToBnL1/4.In this equation, R is the radar detection range (m), Pt is the radar emitting peak power (W), G is the radar antenna gain, σ is the radar cross section (m2), λ is the radar wavelength (m), S/N is the signal to noise ratio of radar receiver, Fn is the noise coefficient, k is the Boltzmann constant (1.38×10-23J/K), To is the radar working temperature (K), Bn is the noise bandwidth (Hz), and L is the system loss.


## 2.1.3. Calculating the Radar Detection Range under Terrain Masking Condition

The calculation diagram of radar detection range under terrain masking condition is shown in Figure1.Figure 1
The calculation diagram of radar detection range under terrain masking condition.The height of the radar observation point ishr. The height of the target point is hf and the target point represents the combat aircraft. The polar coordinate grid P(α,3) is the furthest boundary point of radar detection range in α direction and h(α,3) is the elevation value at this point.In Figure1, not only the height angle εP(α,k) from the polar coordinate grid P(α,k) to radar observation point in α direction can be calculated, but also the height angle εo from the target point to radar observation point can be got, which are shown as(2)εPα,k=arctanhα,k-hrPα,k-Pα,0,εo=arctanhf-hrPα,n-Pα,0. Here, the values of εP(α,k)(k=0,1,2,…,k) are calculated in α direction. If there is a grid P(α,k) satisfying εP(α,k) is greater than εo; then the target point is within the shadow of the radar terrain and the radar could not detect the target.


## 2.2. Active Interference Model of Electronic Warfare

The interference power expression of active electronic warfare equipment is defined as(3)Prj=PjGjGφrjBr4π2Rj2LjBj. Here, Pj is the jammer transmission power (W), Gj is the jammer antenna gain, G(φ) is the gain of the radar antenna in the jammer interference direction, φ is the angle between the direction of the main lobe of the enemy radar antenna and jammer, rj is the jammer polarization loss, Bj is the jammer interference signal bandwidth (Hz), Br is the signal bandwidth of the enemy radar receiver (Hz), Rj is the distance between the enemy radar and the jammer (m), and Lj is the interference signal loss.The power expression of the target echo signal is defined as(4)Prt=PtG2σλ24π3R4L.The meaning of each parameter in the above formula is shown in formula (1), where the constant Prt is equivalent to the parameter S in formula (1).The expression ofG(φ) is defined as(5)Gφ=G,0≤φ≤φ0.5,K2φ0.5φ2G,φ0.5≤φ≤90°,K2φ0.590°2G,90°≤φ≤180°.In this equation, φ0.5 is the lobe width of the enemy radar antenna at the half power point, K is a constant between 0.04 and 0.1 in general, and G is the enemy radar antenna gain.The detection range expression of the enemy radar can be obtained from formulas (1), (3), and (4) in the case of interference with a single jammer, which is shown as(6)R=PtG2σλ24π3LFnkToBnKj+4πLPjGjGφrjBr/Bj2LjBj1/4. Here, Kj is the minimum suppression coefficient of the enemy radar and the remaining parameters are the same as in formulas (1), (3), and (4).When considering suppression interference with multiple jammers against the enemy air defense radar network, the total interference power received is the sum of the interference powers of each jammer, which is defined as(7)Prj_total=∑i=1nPjiGjiGφirjiBr4π2Rji2LjiBji.The detection range expression of the enemy air defense radar after cooperative suppression by multiple jammers can be calculated, which is shown as(8)R=PtG2σλ24π3LFnkToBnKj+∑i=1n4πLPjiGjiGφirjiBr/Bji2LjiBji1/4.


## 2.3. Calculation Model of Route Planning Security Zone

### 2.3.1. Concept of Route Planning Security Zone

Electronic jamming aircraft implements active suppression interference against the enemy air defense radar network, which aims to suppress the detection range of the enemy air defense radar and expand the space scope of route planning security. So it can provide more safe and reliable planning space for subsequent route planning. The definition of route planning security zone is given here. Route planning security zone is a certain width, height range for combat aircraft flight track space which constitutes the search space for the subsequent optimal route planning of combat aircraft. For the sake of convenience, this paper mainly calculates the width of route planning security zone for combat aircraft in a certain altitude. This width refers to the minimum width within the entire security zone.


### 2.3.2. Calculating the Width of Route Planning Security Zone Based on Mathematical Morphology

Since the detection range of the enemy air defense radar is influenced by the terrain, the detection boundary is irregular, so the model building is more difficult. If the traditional geometric method is used to calculate the width of security zone, the calculation process will be very complex and not conducive to engineering practice. Therefore, this paper will proceed from the mathematical morphology point of view so as to calculate the width of route planning security zone. Here the brief introduction to mathematical morphology will be given.The core idea of mathematical morphology is to use a probe structure element to detect a given image and gets the information about the image, so the image can be analyzed and processed [22]. Mathematical morphology involves the basic operation of a closed operation, open operation, corrosion, expansion, and so on. The effect of the corrosion and expansion is opposite. The corrosion can cause a given area to shrink inward at the same time while the expansion allows the region to expand around. Both the closed operation and open operation are complexes of corrosion and expansion. The closed operation expands the image firstly and then corrodes. As a result, a narrow fracture on a given image can be filled and the holes in the image can be removed, and so on. On the contrary, the open operation corrodes the image firstly and then expands, which usually has a smoothing effect on a given image contour and eliminates small burrs on the contour.This paper will be followed by the principle of image expansion and corrosion. Firstly, the image of radar detection range under terrain masking and electronic interference is binarized [22]. Secondly, we will do open operation on the binary image and then expand it. Thirdly, continuity check is performed while the image is being expanded and then the number of the graphic elements in the image is obtained. Finally, it is judged whether the route planning security zone satisfying the minimum width restriction is formed, and then the minimum width of security zone is calculated. The specific calculation process is shown in Figure 2.Figure 2
The calculation process of route planning security zone.


## 2.3.1. Concept of Route Planning Security Zone

Electronic jamming aircraft implements active suppression interference against the enemy air defense radar network, which aims to suppress the detection range of the enemy air defense radar and expand the space scope of route planning security. So it can provide more safe and reliable planning space for subsequent route planning. The definition of route planning security zone is given here. Route planning security zone is a certain width, height range for combat aircraft flight track space which constitutes the search space for the subsequent optimal route planning of combat aircraft. For the sake of convenience, this paper mainly calculates the width of route planning security zone for combat aircraft in a certain altitude. This width refers to the minimum width within the entire security zone.


## 2.3.2. Calculating the Width of Route Planning Security Zone Based on Mathematical Morphology

Since the detection range of the enemy air defense radar is influenced by the terrain, the detection boundary is irregular, so the model building is more difficult. If the traditional geometric method is used to calculate the width of security zone, the calculation process will be very complex and not conducive to engineering practice. Therefore, this paper will proceed from the mathematical morphology point of view so as to calculate the width of route planning security zone. Here the brief introduction to mathematical morphology will be given.The core idea of mathematical morphology is to use a probe structure element to detect a given image and gets the information about the image, so the image can be analyzed and processed [22]. Mathematical morphology involves the basic operation of a closed operation, open operation, corrosion, expansion, and so on. The effect of the corrosion and expansion is opposite. The corrosion can cause a given area to shrink inward at the same time while the expansion allows the region to expand around. Both the closed operation and open operation are complexes of corrosion and expansion. The closed operation expands the image firstly and then corrodes. As a result, a narrow fracture on a given image can be filled and the holes in the image can be removed, and so on. On the contrary, the open operation corrodes the image firstly and then expands, which usually has a smoothing effect on a given image contour and eliminates small burrs on the contour.This paper will be followed by the principle of image expansion and corrosion. Firstly, the image of radar detection range under terrain masking and electronic interference is binarized [22]. Secondly, we will do open operation on the binary image and then expand it. Thirdly, continuity check is performed while the image is being expanded and then the number of the graphic elements in the image is obtained. Finally, it is judged whether the route planning security zone satisfying the minimum width restriction is formed, and then the minimum width of security zone is calculated. The specific calculation process is shown in Figure 2.Figure 2
The calculation process of route planning security zone.


## 2.4. Multiobjective Optimization Model of MACSIA

A typical maximization MOP can be defined as follows [23]:(9)maxf1x,f2x,…,fnxs.t.lb≤x≤ubAeq∗x=beqA∗x≤b.In this equation, n is the dimension of objective space, x∈[x1,x2,…,xDim]T is decision variable, Dim represents the number of decision variables, ub and lb are the upper and lower bounds of the decision variable x, respectively, and A∗x≤b and Aeq∗x=beq are linear inequality constraints and equality constraints for decision variables, respectively.The ultimate goal of optimizing the MACSIA is to achieve a reasonable distribution of jamming aircraft positions and suppress the best against the enemy air defense radar, and the jamming aircraft themselves are not threatened. Therefore, this problem is a MOP. The expressions of single objective are defined as(10)f1=∑i=1Nradarxi-xtc2+yi-ytc2+zi-ztc21/2(11)f2=Widthsafe.Here, the objective function f1 represents the sum of the distances between our jamming aircrafts and the center of the enemy radar network, and the physical meaning is our own safety of jamming aircrafts. The objective function f2 represents the product of the quantity and width of route planning security zone, and the physical meaning is the interference suppression effect of our jamming aircrafts against the enemy air defense radar network. Nradar is the number of the enemy radars. (xi,yi,zi) represents the coordinates of the ith jamming aircraft, (xtc,ytc,ztc) is the center coordinates of the enemy radar network, and Widthsafe is the minimum width of the security zone.The following constraints and interference principles should be considered when the MACSIA is implemented. Firstly, the position of each jamming aircraft must be outside the maximum detection range of the enemy radar. Secondly, the height of each jamming aircraft should be within the given safe height range. Thirdly, the interference principle is traditional principle of many-to-one and one-to-one. The constraints are mathematically expressed as follows:(12)s.t.sqrtxi-xtj2+yi-ytj2+zi-ztj2>Rji=1,2,…,NECM,j=1,2,…,Nradarzmin≤zi≤zmax.Here, zmax and zmin are the upper and lower bounds of the given height range, respectively. The parameter (xi,yi,zi) can be seen in (10), NECM is the quantity of the jamming aircrafts, (xtj,ytj,ztj) represents the coordinates of the ith enemy radar, Nradar is the quantity of the enemy radar, and Rj is the maximum detection range of the ith enemy radar.The multiobjective optimization model of MACSIA with constraints is built by formulas (10), (11), and (12).


## 3. The Improved MOPSO Algorithm

### 3.1. Particle Swarm Optimization

Particle swarm optimization [24] is mainly used to solve single objective optimization problem. The particle position and velocity updating formulas of particle swarm optimization are defined as(13)Xi+1=Xi+Vi+1m,Vi+1=wVi+c1r1Pbesti-Xi+c2r2Gbesti-Xi,where Xi+1 is the position of the particle at the next moment and Xi represents the current position of the particle. Vi+1 is the velocity of the particle at the next moment, Vi represents the current velocity of the particle. c1 and c2 are learning factors. r1 and r2 are random numbers with values between 0 and 1. Pbesti represents the individual optimal solution of the ith particle and Gbesti represents the global optimal solution of the ith population. w is inertia weighting factor which is defined as(14)w=wmax-wmax-wminiteramaxitera.Here, wmax and wmin are the maximum and minimum values of inertia weighting factors, respectively, itera is the current population iteration number, and iteramax represents the maximum number of population iterations.


### 3.2. MOPSO Algorithm

MOPSO [25–30] algorithm uses multiple objective functions as the optimization functions and these objective functions are optimized to achieve the best state at the same time. The core idea of MOPSO can be described as follows: the particle population is initialized. Based on the idea of domination for each particle objective function value, the initial population is divided into two parts: the dominant subset A and the nondominant subset B. An external file is used to store the nondominated solution subset generated by each iteration and only the particle position and velocity in the dominant subset A are updated during the iteration. The dominance relation between the updated particles in A and B is compared. If xi belongs to A, xj belongs to B, and xi dominates xj, then xj is rejected, xi is added to B, and the external archive is updated. The end of the algorithm can be predefined by maximum number of iterations, calculation accuracy, and so on. When the algorithm is finished, the nondominated solution set is the Pareto optimal solution.The flow chart of MOPSO algorithm is shown in Figure3.Figure 3
The flow chart of MOPSO.


### 3.3. The Improved MOPSO Algorithm

#### 3.3.1. Decomposition Mechanism

The decomposition mechanism decomposes a MOP into a number of scalar optimization subproblems and optimizes them simultaneously. Each subproblem is optimized by only using information from its several neighboring subproblems, so it can reduce the computational complexity of the algorithm effectively. In this proposed algorithm, the adopted decomposition mechanism is classical Tchebycheff Approach [31], which is written in the form as follows:(15)gtexλ,z∗=max1≤i≤nλifix-zi∗subjecttox∈Ω,where λ=(λ1,…,λn)T is a weight vector and z∗=(z1∗,…,zn∗)T is the reference point; that is, zi∗=max{fi(x),x∈Ω} for each i=1,…,n.


#### 3.3.2. Proportional Distribution Mechanism

In the MOPSO algorithm, the selection of appropriate local guides from the new found nondominated solutions which aims to attain both convergence and diversity of solutions is a vital problem. If there are only a few nondominated solutions consisted in an area, more particles should be distributed to follow these nondominated solutions as a guide for searching larger variety of solutions around this area. If the situation is reversed, less particles should be distributed to avoid similar solutions gathered in this region. To maintain the diversity of new found nondominated solutions and improve the search ability of population, the proportional distribution mechanism is introduced. The proportional distribution mechanism is depicted as Figure4 [32].Figure 4
Local guides distribution for each particle.In Figure4, the circles represent the particles and the triangles are the nondominated solutions. Each particle and nondominated solution is numbered. The distance between two adjacent nondominated solutions is defined as the variable d. The process of proportional distributed mechanism for the nondominated solutions can be described as follows.Firstly, the existing nondominated solutions and all particles should be given a sorted number based on lateral axial. The number results are shown as Figure4. Secondly, calculate the density parameter Ci of each nondominated solution in the solution space. The formula of Ci is defined as(16)Ci=dxi,xi+1+dxi,xi-12,ifi≠nPF,i≠1,dxi,xi+1,ifi=1dxi,xi-1,ifi=nPF,where xi is the coordinate of nondominated solution i and nPF is the total number of nondominated solutions obtained.Thirdly, each nondominated solution is regarded as the guide distribution of particles. How the amount of particles is guided by each nondominated solution calculated is shown in (17), which is described as GbestNi.(17)GbestNi=integerCi∗Np∑i=1nCi,where Np is the total number of particles and Ci is defined as (16). The example calculation result of GbestN2 is 2∗16/(2+3+3.5+1+1+4)=2. The result indicates that three particles of numbers 1, 2, and 3 are guided by the nondominated solution 1 during the evolution process. In the proportional distribution mechanism, the sum of all GbestNi which are defined as NGbest should be equal to Np for the evolution process. Finally, the particles are guided by the nondominated solutions based on formula (17) and the guided result can be seen in Figure 4.The decomposition mechanism and the proportion distribution mechanism are joined in the MOPSO algorithm flow, which constitutes the calculation process of the improved MOPSO algorithm. Based on the calculation process of the improved MOPSO algorithm, the multiobjective optimization model of MACSIA established in Section2.4 is taken as the objective functions and then combined with the minimum width calculation method of route planning security zone in Section 2.3.2 which can be seen in Figure 2. Therefore, the MACSIA schemes can be obtained that meet the minimum width restriction of route planning security zone and the jamming aircrafts themselves are the most secure.


## 3.1. Particle Swarm Optimization

Particle swarm optimization [24] is mainly used to solve single objective optimization problem. The particle position and velocity updating formulas of particle swarm optimization are defined as(13)Xi+1=Xi+Vi+1m,Vi+1=wVi+c1r1Pbesti-Xi+c2r2Gbesti-Xi,where Xi+1 is the position of the particle at the next moment and Xi represents the current position of the particle. Vi+1 is the velocity of the particle at the next moment, Vi represents the current velocity of the particle. c1 and c2 are learning factors. r1 and r2 are random numbers with values between 0 and 1. Pbesti represents the individual optimal solution of the ith particle and Gbesti represents the global optimal solution of the ith population. w is inertia weighting factor which is defined as(14)w=wmax-wmax-wminiteramaxitera.Here, wmax and wmin are the maximum and minimum values of inertia weighting factors, respectively, itera is the current population iteration number, and iteramax represents the maximum number of population iterations.


## 3.2. MOPSO Algorithm

MOPSO [25–30] algorithm uses multiple objective functions as the optimization functions and these objective functions are optimized to achieve the best state at the same time. The core idea of MOPSO can be described as follows: the particle population is initialized. Based on the idea of domination for each particle objective function value, the initial population is divided into two parts: the dominant subset A and the nondominant subset B. An external file is used to store the nondominated solution subset generated by each iteration and only the particle position and velocity in the dominant subset A are updated during the iteration. The dominance relation between the updated particles in A and B is compared. If xi belongs to A, xj belongs to B, and xi dominates xj, then xj is rejected, xi is added to B, and the external archive is updated. The end of the algorithm can be predefined by maximum number of iterations, calculation accuracy, and so on. When the algorithm is finished, the nondominated solution set is the Pareto optimal solution.The flow chart of MOPSO algorithm is shown in Figure3.Figure 3
The flow chart of MOPSO.


## 3.3. The Improved MOPSO Algorithm

### 3.3.1. Decomposition Mechanism

The decomposition mechanism decomposes a MOP into a number of scalar optimization subproblems and optimizes them simultaneously. Each subproblem is optimized by only using information from its several neighboring subproblems, so it can reduce the computational complexity of the algorithm effectively. In this proposed algorithm, the adopted decomposition mechanism is classical Tchebycheff Approach [31], which is written in the form as follows:(15)gtexλ,z∗=max1≤i≤nλifix-zi∗subjecttox∈Ω,where λ=(λ1,…,λn)T is a weight vector and z∗=(z1∗,…,zn∗)T is the reference point; that is, zi∗=max{fi(x),x∈Ω} for each i=1,…,n.


### 3.3.2. Proportional Distribution Mechanism

In the MOPSO algorithm, the selection of appropriate local guides from the new found nondominated solutions which aims to attain both convergence and diversity of solutions is a vital problem. If there are only a few nondominated solutions consisted in an area, more particles should be distributed to follow these nondominated solutions as a guide for searching larger variety of solutions around this area. If the situation is reversed, less particles should be distributed to avoid similar solutions gathered in this region. To maintain the diversity of new found nondominated solutions and improve the search ability of population, the proportional distribution mechanism is introduced. The proportional distribution mechanism is depicted as Figure4 [32].Figure 4
Local guides distribution for each particle.In Figure4, the circles represent the particles and the triangles are the nondominated solutions. Each particle and nondominated solution is numbered. The distance between two adjacent nondominated solutions is defined as the variable d. The process of proportional distributed mechanism for the nondominated solutions can be described as follows.Firstly, the existing nondominated solutions and all particles should be given a sorted number based on lateral axial. The number results are shown as Figure4. Secondly, calculate the density parameter Ci of each nondominated solution in the solution space. The formula of Ci is defined as(16)Ci=dxi,xi+1+dxi,xi-12,ifi≠nPF,i≠1,dxi,xi+1,ifi=1dxi,xi-1,ifi=nPF,where xi is the coordinate of nondominated solution i and nPF is the total number of nondominated solutions obtained.Thirdly, each nondominated solution is regarded as the guide distribution of particles. How the amount of particles is guided by each nondominated solution calculated is shown in (17), which is described as GbestNi.(17)GbestNi=integerCi∗Np∑i=1nCi,where Np is the total number of particles and Ci is defined as (16). The example calculation result of GbestN2 is 2∗16/(2+3+3.5+1+1+4)=2. The result indicates that three particles of numbers 1, 2, and 3 are guided by the nondominated solution 1 during the evolution process. In the proportional distribution mechanism, the sum of all GbestNi which are defined as NGbest should be equal to Np for the evolution process. Finally, the particles are guided by the nondominated solutions based on formula (17) and the guided result can be seen in Figure 4.The decomposition mechanism and the proportion distribution mechanism are joined in the MOPSO algorithm flow, which constitutes the calculation process of the improved MOPSO algorithm. Based on the calculation process of the improved MOPSO algorithm, the multiobjective optimization model of MACSIA established in Section2.4 is taken as the objective functions and then combined with the minimum width calculation method of route planning security zone in Section 2.3.2 which can be seen in Figure 2. Therefore, the MACSIA schemes can be obtained that meet the minimum width restriction of route planning security zone and the jamming aircrafts themselves are the most secure.


## 3.3.1. Decomposition Mechanism

The decomposition mechanism decomposes a MOP into a number of scalar optimization subproblems and optimizes them simultaneously. Each subproblem is optimized by only using information from its several neighboring subproblems, so it can reduce the computational complexity of the algorithm effectively. In this proposed algorithm, the adopted decomposition mechanism is classical Tchebycheff Approach [31], which is written in the form as follows:(15)gtexλ,z∗=max1≤i≤nλifix-zi∗subjecttox∈Ω,where λ=(λ1,…,λn)T is a weight vector and z∗=(z1∗,…,zn∗)T is the reference point; that is, zi∗=max{fi(x),x∈Ω} for each i=1,…,n.


## 3.3.2. Proportional Distribution Mechanism

In the MOPSO algorithm, the selection of appropriate local guides from the new found nondominated solutions which aims to attain both convergence and diversity of solutions is a vital problem. If there are only a few nondominated solutions consisted in an area, more particles should be distributed to follow these nondominated solutions as a guide for searching larger variety of solutions around this area. If the situation is reversed, less particles should be distributed to avoid similar solutions gathered in this region. To maintain the diversity of new found nondominated solutions and improve the search ability of population, the proportional distribution mechanism is introduced. The proportional distribution mechanism is depicted as Figure4 [32].Figure 4
Local guides distribution for each particle.In Figure4, the circles represent the particles and the triangles are the nondominated solutions. Each particle and nondominated solution is numbered. The distance between two adjacent nondominated solutions is defined as the variable d. The process of proportional distributed mechanism for the nondominated solutions can be described as follows.Firstly, the existing nondominated solutions and all particles should be given a sorted number based on lateral axial. The number results are shown as Figure4. Secondly, calculate the density parameter Ci of each nondominated solution in the solution space. The formula of Ci is defined as(16)Ci=dxi,xi+1+dxi,xi-12,ifi≠nPF,i≠1,dxi,xi+1,ifi=1dxi,xi-1,ifi=nPF,where xi is the coordinate of nondominated solution i and nPF is the total number of nondominated solutions obtained.Thirdly, each nondominated solution is regarded as the guide distribution of particles. How the amount of particles is guided by each nondominated solution calculated is shown in (17), which is described as GbestNi.(17)GbestNi=integerCi∗Np∑i=1nCi,where Np is the total number of particles and Ci is defined as (16). The example calculation result of GbestN2 is 2∗16/(2+3+3.5+1+1+4)=2. The result indicates that three particles of numbers 1, 2, and 3 are guided by the nondominated solution 1 during the evolution process. In the proportional distribution mechanism, the sum of all GbestNi which are defined as NGbest should be equal to Np for the evolution process. Finally, the particles are guided by the nondominated solutions based on formula (17) and the guided result can be seen in Figure 4.The decomposition mechanism and the proportion distribution mechanism are joined in the MOPSO algorithm flow, which constitutes the calculation process of the improved MOPSO algorithm. Based on the calculation process of the improved MOPSO algorithm, the multiobjective optimization model of MACSIA established in Section2.4 is taken as the objective functions and then combined with the minimum width calculation method of route planning security zone in Section 2.3.2 which can be seen in Figure 2. Therefore, the MACSIA schemes can be obtained that meet the minimum width restriction of route planning security zone and the jamming aircrafts themselves are the most secure.


## 4. Experiment Analysis

432
k
m
×
432
k
m area is selected as combat scenario in our simulation experiment. The DEM of the grid structure in this area is shown in Figure 5. The DEM data resolution is 360 m. We assume that four enemy air defense radars are deployed in the area. The minimum suppression coefficient Kj of the four radars is five. Their performance parameters are shown in Table 1. The coordinates of the four radars are (240km,100km), (150km,120km), (300km,250km), and (180km,280km), respectively. In order to obtain at least 20km width of the route planning security zone, we use three electronic jamming aircrafts to implement cooperative active suppression interference against the enemy radars. The three jamming aircrafts are at a height of 2.1km and their performance parameters are shown in Table 2. The interference principle of one-to-one is used to interfere with the enemy radars.Table 1
The performance parameters list of four air defense radars.
Radar
P
t
/
k
W
G
λ
/
m
L
F
n
T
o
/
K
B
n/Hz
Number 1
4000
50
0.6
2.5
3
291
2
×
1
0
5
Number 2
3900
45
0.8
4
2.5
291
2.1
×
1
0
5
Number 3
3100
70
0.75
3
4
291
3
×
1
0
5
Number 4
2600
100
0.5
3
6
291
1.5
×
1
0
5Table 2
The parameters list of three jamming aircrafts.
Jamming aircraft
P
j
/
k
W
G
j
r
j
L
j
B
j
/
H
z
Number 1
1
×
1
0
3
4
2
20
2
×
1
0
6
Number 2
2
×
1
0
3
2
2
20
2
×
1
0
6
Number 3
2
×
1
0
3
6
1
30
2
×
1
0
6Figure 5
Digital elevation terrain map.Combined with the above method of calculating the radar detection range under the terrain masking condition, the detection range of 2100 m height for the four enemy air defense radars under the terrain masking condition is shown in Figure6. It can be seen from Figure 6 that the four radars detection range overlaps each other. So it is difficult for our combat aircraft to penetrate in the past safely. Therefore, our combat aircraft must have the aid of suppression interference against air defense radar network by the electronic support jamming aircraft. The binary image of the radar network detection range is shown in Figure 7.Figure 6
The detection range of the four enemy air defense radars under the terrain masking condition.Figure 7
The binary image of radar network detection range.In order to facilitate the comparison of the MOPSO algorithm and the improved MOPSO algorithm, the initial parameters of the two algorithms are the same which are set as follows:Dim=6, the population number xSize is 50, the maximum number of iterations Max_itera is 100, the learning factors c1 and c2 are both 0.8, and the maximum value wmax and the minimum value wmin of the inertial weighting factors are 1.2 and 0.1, respectively. The maximum number of iterations Max_itera is 100.The simulation experiment is carried out by MATLAB language and the nondominated solutions distribution obtained is shown in Figure8.Figure 8
The distribution graph of nondominated solutions obtained.It can be seen from Figure8 that the nondominated solutions obtained by the improved MOPSO algorithm obviously dominate the results of MOPSO, and the distribution of the nondominated solutions obtained by the improved MOPSO is more uniform. The number of the nondominated solutions obtained by MOPSO is 17 and the number of the nondominated solutions obtained by the improved MOPSO is 28, which constitute their Pareto front, respectively. It is obvious that the Pareto front obtained by the improved MOPSO is more smooth than that of MOPSO. So the improved MOPSO algorithm has achieved a good search result than MOPSO.For the improved MOPSO algorithm, there are 28 kinds of MACSIA schemes for the jamming aircrafts from the 28 Pareto optimal solutions. Which solution will be chosen requires the decision maker to make a choice based on the actual battlefield situation. Three typical cases are given here. Firstly, when the decision maker needs the sum of the distance between each jamming aircraft and the center of enemy radar network to be minimum and the minimum width of security zone to be maximum, that is, more focused on the suppression interference effect of the enemy radar network. So the optimal MACSIA scheme obtained by the improved MOPSO algorithm which also is called scheme one is shown in Figure10. It can be seen from Figure 10 that the three jamming aircrafts are closer to the center of enemy radar network. The binary image of suppression effect for each jamming aircraft against the enemy radar network effect in scheme one is shown in Figure 12. The minimum width of security zone obtained by the improved MOPSO algorithm is 45km. Figure 9 is the optimal MACSIA scheme obtained by MOSPO under the same condition and Figure 11 is the binary image of suppression effect. The minimum width of security zone obtained by MOPSO is 41km which is less than that obtained by the improved MOPSO algorithm. So it is proved that the improved MOPSO algorithm is better than MOPSO algorithm.Figure 9
The optimal MACSIA scheme of the jamming aircrafts under the minimum total distance and maximum security zone minimum width obtained by MOPSO.Figure 10
The optimal MACSIA scheme of the jamming aircrafts under the minimum total distance and maximum security zone minimum width obtained by improved MOPSO.Figure 11
The binary image of suppression effect for each jamming aircraft against the enemy radar network effect obtained by MOPSO.Figure 12
The binary image of suppression effect for each jamming aircraft against the enemy radar network effect in scheme one obtained by improved MOPSO.Secondly, when the decision maker needs the sum of the distance between each jamming aircraft and the center of enemy radar network to be maximum and the minimum width of security zone to be minimum, that is, more focused on the safety of each jamming aircraft itself. So the optimal MACSIA scheme obtained by the improved MOPSO algorithm which also is called scheme two is shown in Figure14. It can be seen from Figure 14 that the three jamming aircrafts are farther away from the center of enemy radar network. The binary image of suppression effect for each jamming aircraft against the enemy radar network effect in scheme two is shown in Figure 16. The minimum width of security zone obtained by the improved MOPSO algorithm is 28km. So the suppression interference effect of scheme two is worse than that of scheme one. Figure 13 is the optimal MACSIA scheme obtained by MOSPO under the same condition and Figure 15 is the binary image of suppression effect. The minimum width of security zone obtained by MOPSO is 24km which is less than that obtained by the improved MOPSO algorithm. So it is proved that the improved MOPSO algorithm is better than MOPSO algorithm.Figure 13
The optimal MACSIA scheme of the jamming aircrafts under the maximum total distance and minimum security zone minimum width obtained by MOPSO.Figure 14
The optimal MACSIA scheme of the jamming aircrafts under the maximum total distance and minimum security zone minimum width obtained by improved MOPSO.Figure 15
The binary image of suppression effect for each jamming aircraft against the enemy radar network effect obtained by MOPSO.Figure 16
The binary image of suppression effect for each jamming aircraft against the enemy radar network effect in scheme two obtained by improved MOPSO.Thirdly, when the decision maker should consider not only the safety of the jamming aircrafts themselves, but also a good suppression interference effect. Then the optimal MACSIA scheme obtained by the improved MOPSO algorithm with moderate total distance and minimum width of security zone can be chosen which also is called scheme three; it is shown in Figure18. The binary image of suppression effect for each jamming aircraft against the enemy radar network effect in scheme three is shown in Figure 20. The minimum width of security zone obtained by the improved MOPSO algorithm is 35.5km. So the suppression interference effect of scheme three is worse than that of scheme one but better than that of scheme two. Figure 17 is the optimal MACSIA scheme obtained by MOSPO under the same condition and Figure 19 is the binary image of suppression effect. The minimum width of security zone obtained by MOPSO is 32km which is less than that obtained by the improved MOPSO algorithm. So it is proved that the improved MOPSO algorithm is better than MOPSO algorithm.Figure 17
The optimal MACSIA scheme of the jamming aircrafts under the moderate total distance and moderate security zone minimum width obtained by MOPSO.Figure 18
The optimal MACSIA scheme of the jamming aircrafts under the moderate total distance and moderate security zone minimum width obtained by improved MOPSO.Figure 19
The binary image of suppression effect for each jamming aircraft against the enemy radar network effect obtained by MOPSO.Figure 20
The binary image of suppression effect for each jamming aircraft against the enemy radar network effect in scheme three obtained by improved MOPSO.It can be seen from the above that the three kinds of optimal MACSIA schemes of the jamming aircrafts obtained by the improved MOPSO algorithm have achieved the given suppression interference effect and formed the security zone that meets the given requirements.


## 5. Conclusions

(
1
) In this paper, for the problem of MACSIA in electronic warfare mission planning, the concept of route planning security zone is proposed, the solution to get the minimum width of security zone based on mathematical morphology is put forward, and the multiobjective optimization model of electronic interference array is built.(
2
) The improved MOPSO algorithm is used to solve the model.(
3
) Then the optimal MACSIA schemes can be obtained that meet the minimum width restriction of route planning security zone and the jamming aircrafts themselves are the most secure. Meanwhile, it verifies that the built model and the proposed improved MOPSO algorithm are feasible and effective.(
4
) The content of this paper is an important part of electronic warfare mission planning. It is the basis for subsequent penetration route planning which provides a safe and reliable planning space for the subsequent route planning. Therefore, this study has important practical significance.When the battlefield environment changes, research on the problem of dynamic MACSIA has more practical significance. So we will take this as an entry point and explore more efficient and feasible dynamic MOPSO algorithm to solve the problem in the future work.



---

*Source: 9843735-2017-03-05.xml*