# Dynamic Multiaircraft Cooperative Suppression Interference Array Optimization by Dynamic MOPSO Algorithm

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

---

## Abstract

Dynamic multiaircraft cooperative suppression interference array (MACSIA) optimization problem is a typical dynamic multiobjective optimization problem. In this paper, the sum of the distance between each jamming aircraft and the enemy air defense radar network center and the minimum width of the safety area for route planning are taken as the objective functions. The dynamic changes in the battlefield environment are reduced to two cases. One is that the location of the enemy air defense radar is mobile, but the number remains the same. The other is that the number of the enemy air defense radars is variable, but the original location remains unchanged. Thus, two dynamic multiobjective optimization models of dynamic MACSIA are constructed. The dynamic multiobjective particle swarm optimization algorithm is used to solve the two models, respectively. The optimal dynamic MACSIA schemes which satisfy the limitation of the given suppression interference effect and ensure the safety of the jamming aircraft themselves are obtained by simulation experiments. And then verify the correctness of the constructed dynamic multiobjective optimization model, as well as the feasibility and effectiveness of the dynamic multiobjective particle swarm optimization algorithm in solving dynamic MACSIA problem.


---

## Body

## 1. Introduction

Under the conditions of high-tech warfare, the air defense radar network on the enemy air defense position usually has a strong ability to resist ‘the four’, that is, antilow altitude intrusion, antijamming, antistealth, and antiradiation missile [1, 2]. It is usually difficult for combat aircraft to try to break through the enemy air defense radar network position safely. At this time the combat aircraft must rely on multiple jamming aircraft to implement cooperative suppression interference against the enemy air defense radar network. So that it can provide a secure and planned space for the follow-up route planning of the combat aircraft. In order to make more rational use of limited interference resources, a reasonable array mode should be considered for the location of each jamming aircraft in the process of cooperative interference. Since the actual battlefield environment is constantly changing, it must be considered how to adjust the location of each jamming aircraft based on the actual situation to achieve the best array mode when the battlefield changes. And this problem is the dynamic MACSIA optimization problem in electronic warfare mission planning.For this problem, in [3], a multiconstraint genetic algorithm is used to solve the deploy optimization problem of the enemy radar network. In [4], Ruan et al. analyzed the various factors that affect the effect of suppression interference and selected the minimum interference distance as a measure for the effect of suppression interference. A weighted sum method was adopted for constructing a single-objective MACSIA model to study the optimal electronic warfare array against radar network [5]. However, the shortcoming of the single-objective optimization model is obvious. Firstly, the determination of the weight depends on the relative importance of each subobjective. Secondly, the experimenter needs to set different weights by multiple experiments to obtain different solutions and the computational overhead is larger. In a work by Zhang et al. [6], a multiobjective optimization model of MACSIA was constructed under static conditions. A modified multiobjective particle swarm optimization algorithm is used to solve the static MACSIA optimization problem [7]. Literatures 6 and 7 are our previous research results, but only the MACSIA optimization problem under static conditions is studied. The MACSIA optimization problem is not solved as a dynamic multiobjective optimization problem [8, 9] and the impact of changes in the battlefield environment is ignored.In view of the above problems, in this paper, from the perspective of the real battlefield environment, the dynamic changes in the battlefield environment are reduced to two cases. One is that the location of the enemy air defense radar is mobile, but the number remains the same. The other is that the number of the enemy air defense radars is variable, but the original location remains unchanged. As a basis, on the basis of fully considering the dynamic multiobjective optimization characteristics of the dynamic MACSIA, the sum of the distance between each jamming aircraft and the enemy air defense radar network center and the minimum width of the safety area for route planning are taken as the objective functions. So two dynamic multiobjective optimization models of dynamic MACSIA are constructed. Dynamic multiobjective particle swarm optimization (DMOPSO) [10–18] algorithm has a better effect on the optimization of dynamic environment. It can obtain the optimal solution quickly and effectively. As the algorithm is simple, so the implementation of the project is not difficult and it is widely applied in dynamic multiobjective optimization field [19–22]. Therefore, this paper uses DMOPSO algorithm to solve the two dynamic multiobjective optimization models of dynamic MACSIA. The optimal dynamic MACSIA schemes of each jamming aircraft against the enemy air defense network are calculated, which can provide some reference for the decision-maker to make the next decision.The novelty of the paper is mainly reflected in the following aspects.Firstly, to the best of our knowledge, this is the first work that deals with the MACSIA optimization problem as a dynamic multiobjective optimization problem. However, the previous work has considered the problem as a single-objective optimization or static multiobjective optimization problem to solve. The impact of changes in the battlefield environment was ignored.Secondly, this paper starts with the characteristics of the real battlefield environment and considers two different environmental changes. One is that the location of the enemy air defense radar is mobile, but the number remains the same. The other is that the number of the enemy air defense radars is variable, but the original location remains unchanged. Thus, two dynamic multiobjective optimization models of dynamic MACSIA are constructed.Thirdly, the DMOPSO algorithm is used to solve these two models, and the results obtained meet the requirements of the mission. At the same time, the correctness of the model is verified and the feasibility of the DMOPSO algorithm to solve the problem is verified.Finally, when the battlefield environment changes, it can provide a basis for decision-maker to make decisions. Therefore, the research on the dynamic MACSIA optimization problem has more practical significance.


## 2. Model Constructing

### 2.1. Dynamic Multiobjective Optimization Model of Dynamic MACSIA

A typical maximized dynamic multiobjective optimization problem can be described as the following mathematical expression [23, 24].(1)maxfx,t=f1x,t,f2x,t,⋯,fmx,ts.t.gx,t>0,hx,t=0In this equation, f(x,t) is the set of objective functions with the time variable t, x∈[x1,x2,⋯,xDim]T represents the decision variable, Dim is the number of decision variables, m is the number of objective functions, and g(x,t) and h(x,t) represent the constraint functions with the time variable t.The ultimate goal of dynamic MACSIA optimization is to achieve a reasonable array of jamming aircraft when the battlefield environment changes, which achieve the best effect of suppression interference against the enemy radar network and ensure the safety of the jamming aircraft themselves. So we can get two dynamic multiobjective optimization models of dynamic MACSIA under two kinds of battlefield environment are as follows.The first case: the location of the enemy air defense radar is mobile, but the number remains the same.(2)f1=∑i=1NECMxi-xtct2+yi-ytct2+zi-ztct21/2f2=Widthsafes.t.sqrtxi-xtjt2+yi-ytjt2+zi-ztjt2>Rji=1,2,⋯,NECM;j=1,2,⋯,Nradarzmin≤zi≤zmaxHere, the objective function f2 represents the minimum width of the safety area for route planning, and the physical meaning is the suppression interference effect of jamming aircraft against the enemy radar network. The objective function f1 represents the sum of the distance between each jamming aircraft and the enemy air defense radar network center, and the physical meaning characterizes the safety of jamming aircraft themselves. (xtc(t),ytc(t),ztc(t)) indicates the coordinate of the enemy air defense radar center; (xi,yi,zi) indicates the coordinate of the i-th jamming aircraft. In the first case, the location of the enemy air defense radar is mobile, so the coordinate of the enemy air defense radar network is variable and the function of time variable t. Widthsafe is the minimum width of route planning safety area. In the constraints, the location of each jamming aircraft should be outside the maximum detection range of each enemy air defense radar, the height of each jamming aircraft should be within the required safety altitude range, and zmax and zmin represent the maximum and minimum values for a given height, respectively. NECM indicates the number of jamming aircraft, (xtj(t),xtj(t),xtj(t)) is the coordinate of the i-th enemy air defense radar, and the function of time variable t. Rj indicates the maximum detection range of the i-th enemy air defense radar; Nradar represents the number of radars in the enemy air defense radar network.The second case: the number of the enemy air defense radar is variable, but the original location remains unchanged.In the second case, as the number of the enemy air defense radars is variable, so the number of our jamming aircraft will change accordingly and the location of the enemy air defense radar network center is also changing.(3)f1=∑i=1NECMtxi-xtct2+yi-ytct2+zi-ztct21/2f2=Widthsafes.t.sqrtxi-xtjt2+yi-ytjt2+zi-ztjt2>Rji=1,2,⋯,NECMt;j=1,2,⋯,Nradartzmin≤zi≤zmaxIn this equation, NECM(t) indicates the number of jamming aircraft and is the function of time variable t. Nradar(t) represents the number of radars in the enemy air defense radar network and is the function of time variable t. The meaning of the remaining parameters is the same as in equation (1).


### 2.2. Calculation of Route Planning Safety Area Width

The definition of route planning safety area refers to a range of width and height of the combat aircraft within the flight path space and constitutes the follow-up optimal route planning search space for the combat aircraft. In order to facilitate the study, this paper still mainly considers the calculation of route planning safe area width in a certain height of the scope for the combat aircraft. And the width refers to the minimum width of the entire route planning safety area. The method of solving the width of the safe area is still based on the mathematical morphology [25] method. Due to the detection boundary shape of the detection range of the enemy air defense radar network is irregular under terrain shading and electronic jamming conditions, so the traditional geometric method is difficult to calculate the safety area width. The core of mathematical morphology is that you can use a probe structure element to detect the image and obtain the relevant information, so as to achieve the purpose of image analysis and processing. Firstly, the images of the enemy air defense radar network are binarized according to the terrain shelter and the electronic jamming conditions. Secondly, using the principle of graphic compression and expansion in mathematical morphology, the binary image of the radar detection range is opened and then the expansion operation is carried out, and the image connectivity is continuously checked. Finally, determine whether the minimum width of the route planning safety zone satisfying the given width limit is formed and the minimum width is calculated. On this part of the content which can also be seen in our previous research results [6, 7], the specific calculating steps are shown in Figure 1.Figure 1
The minimum width calculation process of route planning safety area.The calculation model of the range of enemy air defense radar under terrain shading and electronic jamming in the process is as follows [6, 7].The detection range calculation of the enemy air defense radar under terrain shading is as follows.Figure2 shows the detection range calculation diagram of the enemy air defense radar under terrain shading.Figure 2
The detection range calculation diagram of the enemy air defense radar under terrain shading.In Figure2, the combat aircraft is the target point, the height of which is expressed by hf. In α direction, the farthest boundary point of the enemy radar is represented by the polar coordinate grid P(α,3), where the elevation value is denoted by h(α,3). hr represents the height of radar observation point. The height angle ε0 from the enemy radar observation point to the target and the height angle εP(α,k) from the enemy radar observation to P(α,k) can be obtained in Figure 2, respectively. As shown in (4), calculate εP(α,k)(k=0,1,2⋯k) in turn and compare with εo, respectively. If there is a P(α,k) value satisfying εP(α,k)>εo, then the enemy radar observation line of sight is blocked by the terrain and the target can not be detected.(4)εPα,k=arctan⁡hα,k-hrPα,k-Pα,0εo=arctan⁡hf-hrPα,n-Pα,0The detection range calculation of the enemy air defense radar under electronic jamming is as follows.The detection range expression of multiple jamming aircraft cooperative suppression interference against the enemy air defense radar network is shown in(5)R=PtG2σλ24π3LFnkToBnKj+∑i=1n4πLPjiGjiGφirjiBr/Rji2LjiBjiKj1/4Here, Rji is the distance between the enemy radar and the i-th jamming aircraft and its unit is m, Lji represents the interference loss of the i-th jamming aircraft, Bji indicates the interference signal bandwidth of the i-th jamming aircraft and its unit is Hz, and Pji is the transmit power of the i-th jamming aircraft and its unit is W. Gji represents the antenna gain of the i-th jamming aircraft, rji indicates the polarization loss of the i-th jamming aircraft, Br represents the signal bandwidth received by the enemy radar receiver and its unit is Hz, L is the system loss of the enemy radar, Bn indicates the noise bandwidth of the enemy radar, T0 is the current working environment temperature of the enemy radar and its unit is K, k represents the Boltzmann constant which is 1.38×10-23J/K, Fn indicates the noise coefficient of the enemy radar, λ represents the wavelength of the enemy radar and its unit is m, σ indicates the target cross-sectional area of the enemy radar and its unit is m2, G is the antenna gain of the enemy radar, Pt represents the emission peak power of the enemy radar and its unit is W, and n indicates the number of the jamming aircraft. When n is 1, formula (5) shows the detection distance of single jamming aircraft suppression interference against the enemy radar. Kj is the minimum suppression factor of the enemy radar, G(φi) represents the gain of the enemy radar antenna in the interference direction of the i-th jamming aircraft, and the expression is shown in (6)Gφ=G,0≤φ≤φ0.5K2φ0.5φ2G,φ0.5≤φ≤90°K2φ0.590°2G,90°≤φ≤180°In the above equation, the value range of K is generally 0.04 to 0.1 and φ0.5 indicates the lobe width of the enemy radar antenna at the half power point.


## 2.1. Dynamic Multiobjective Optimization Model of Dynamic MACSIA

A typical maximized dynamic multiobjective optimization problem can be described as the following mathematical expression [23, 24].(1)maxfx,t=f1x,t,f2x,t,⋯,fmx,ts.t.gx,t>0,hx,t=0In this equation, f(x,t) is the set of objective functions with the time variable t, x∈[x1,x2,⋯,xDim]T represents the decision variable, Dim is the number of decision variables, m is the number of objective functions, and g(x,t) and h(x,t) represent the constraint functions with the time variable t.The ultimate goal of dynamic MACSIA optimization is to achieve a reasonable array of jamming aircraft when the battlefield environment changes, which achieve the best effect of suppression interference against the enemy radar network and ensure the safety of the jamming aircraft themselves. So we can get two dynamic multiobjective optimization models of dynamic MACSIA under two kinds of battlefield environment are as follows.The first case: the location of the enemy air defense radar is mobile, but the number remains the same.(2)f1=∑i=1NECMxi-xtct2+yi-ytct2+zi-ztct21/2f2=Widthsafes.t.sqrtxi-xtjt2+yi-ytjt2+zi-ztjt2>Rji=1,2,⋯,NECM;j=1,2,⋯,Nradarzmin≤zi≤zmaxHere, the objective function f2 represents the minimum width of the safety area for route planning, and the physical meaning is the suppression interference effect of jamming aircraft against the enemy radar network. The objective function f1 represents the sum of the distance between each jamming aircraft and the enemy air defense radar network center, and the physical meaning characterizes the safety of jamming aircraft themselves. (xtc(t),ytc(t),ztc(t)) indicates the coordinate of the enemy air defense radar center; (xi,yi,zi) indicates the coordinate of the i-th jamming aircraft. In the first case, the location of the enemy air defense radar is mobile, so the coordinate of the enemy air defense radar network is variable and the function of time variable t. Widthsafe is the minimum width of route planning safety area. In the constraints, the location of each jamming aircraft should be outside the maximum detection range of each enemy air defense radar, the height of each jamming aircraft should be within the required safety altitude range, and zmax and zmin represent the maximum and minimum values for a given height, respectively. NECM indicates the number of jamming aircraft, (xtj(t),xtj(t),xtj(t)) is the coordinate of the i-th enemy air defense radar, and the function of time variable t. Rj indicates the maximum detection range of the i-th enemy air defense radar; Nradar represents the number of radars in the enemy air defense radar network.The second case: the number of the enemy air defense radar is variable, but the original location remains unchanged.In the second case, as the number of the enemy air defense radars is variable, so the number of our jamming aircraft will change accordingly and the location of the enemy air defense radar network center is also changing.(3)f1=∑i=1NECMtxi-xtct2+yi-ytct2+zi-ztct21/2f2=Widthsafes.t.sqrtxi-xtjt2+yi-ytjt2+zi-ztjt2>Rji=1,2,⋯,NECMt;j=1,2,⋯,Nradartzmin≤zi≤zmaxIn this equation, NECM(t) indicates the number of jamming aircraft and is the function of time variable t. Nradar(t) represents the number of radars in the enemy air defense radar network and is the function of time variable t. The meaning of the remaining parameters is the same as in equation (1).


## 2.2. Calculation of Route Planning Safety Area Width

The definition of route planning safety area refers to a range of width and height of the combat aircraft within the flight path space and constitutes the follow-up optimal route planning search space for the combat aircraft. In order to facilitate the study, this paper still mainly considers the calculation of route planning safe area width in a certain height of the scope for the combat aircraft. And the width refers to the minimum width of the entire route planning safety area. The method of solving the width of the safe area is still based on the mathematical morphology [25] method. Due to the detection boundary shape of the detection range of the enemy air defense radar network is irregular under terrain shading and electronic jamming conditions, so the traditional geometric method is difficult to calculate the safety area width. The core of mathematical morphology is that you can use a probe structure element to detect the image and obtain the relevant information, so as to achieve the purpose of image analysis and processing. Firstly, the images of the enemy air defense radar network are binarized according to the terrain shelter and the electronic jamming conditions. Secondly, using the principle of graphic compression and expansion in mathematical morphology, the binary image of the radar detection range is opened and then the expansion operation is carried out, and the image connectivity is continuously checked. Finally, determine whether the minimum width of the route planning safety zone satisfying the given width limit is formed and the minimum width is calculated. On this part of the content which can also be seen in our previous research results [6, 7], the specific calculating steps are shown in Figure 1.Figure 1
The minimum width calculation process of route planning safety area.The calculation model of the range of enemy air defense radar under terrain shading and electronic jamming in the process is as follows [6, 7].The detection range calculation of the enemy air defense radar under terrain shading is as follows.Figure2 shows the detection range calculation diagram of the enemy air defense radar under terrain shading.Figure 2
The detection range calculation diagram of the enemy air defense radar under terrain shading.In Figure2, the combat aircraft is the target point, the height of which is expressed by hf. In α direction, the farthest boundary point of the enemy radar is represented by the polar coordinate grid P(α,3), where the elevation value is denoted by h(α,3). hr represents the height of radar observation point. The height angle ε0 from the enemy radar observation point to the target and the height angle εP(α,k) from the enemy radar observation to P(α,k) can be obtained in Figure 2, respectively. As shown in (4), calculate εP(α,k)(k=0,1,2⋯k) in turn and compare with εo, respectively. If there is a P(α,k) value satisfying εP(α,k)>εo, then the enemy radar observation line of sight is blocked by the terrain and the target can not be detected.(4)εPα,k=arctan⁡hα,k-hrPα,k-Pα,0εo=arctan⁡hf-hrPα,n-Pα,0The detection range calculation of the enemy air defense radar under electronic jamming is as follows.The detection range expression of multiple jamming aircraft cooperative suppression interference against the enemy air defense radar network is shown in(5)R=PtG2σλ24π3LFnkToBnKj+∑i=1n4πLPjiGjiGφirjiBr/Rji2LjiBjiKj1/4Here, Rji is the distance between the enemy radar and the i-th jamming aircraft and its unit is m, Lji represents the interference loss of the i-th jamming aircraft, Bji indicates the interference signal bandwidth of the i-th jamming aircraft and its unit is Hz, and Pji is the transmit power of the i-th jamming aircraft and its unit is W. Gji represents the antenna gain of the i-th jamming aircraft, rji indicates the polarization loss of the i-th jamming aircraft, Br represents the signal bandwidth received by the enemy radar receiver and its unit is Hz, L is the system loss of the enemy radar, Bn indicates the noise bandwidth of the enemy radar, T0 is the current working environment temperature of the enemy radar and its unit is K, k represents the Boltzmann constant which is 1.38×10-23J/K, Fn indicates the noise coefficient of the enemy radar, λ represents the wavelength of the enemy radar and its unit is m, σ indicates the target cross-sectional area of the enemy radar and its unit is m2, G is the antenna gain of the enemy radar, Pt represents the emission peak power of the enemy radar and its unit is W, and n indicates the number of the jamming aircraft. When n is 1, formula (5) shows the detection distance of single jamming aircraft suppression interference against the enemy radar. Kj is the minimum suppression factor of the enemy radar, G(φi) represents the gain of the enemy radar antenna in the interference direction of the i-th jamming aircraft, and the expression is shown in (6)Gφ=G,0≤φ≤φ0.5K2φ0.5φ2G,φ0.5≤φ≤90°K2φ0.590°2G,90°≤φ≤180°In the above equation, the value range of K is generally 0.04 to 0.1 and φ0.5 indicates the lobe width of the enemy radar antenna at the half power point.


## 3. Dynamic Multiobjective Particle Swarm Optimization Algorithm

### 3.1. Particle Swarm Optimization

Particle swarm optimization (PSO) [26] algorithm is widely used in solving single target optimization problem. The updating location and speed expressions of particles in PSO are shown in (7)Xi+1=Xi+Vi+1Vi+1=wVi+c1r1Pbesti-Xi+c2r2Gbesti-XiIn this equation, Vi indicates the current speed of the particle, Vi+1 represents the speed at the next moment, Xi indicates the current location of the particle, the position of the particle at the next moment is denoted by Xi+1, Gbesti is the global optimal solution for the i-th generation of population, Pbesti represents the individual optimal solution for the i-th generation of particle, and r1 and r2 are random numbers with values between 0 and 1. Both c1 and c2 represent learning factors, w indicates the inertia weight factor, which is a variable that reflects the effect of the current velocity of the particle on the velocity of the particle at the next moment, and the expression is (8)w=wmax-wmax-wminiteramax∗iteraHere, itera is the current number of population iterations, iteramax represents the maximum number of population iterations, and wmin and wmax represent the minimum and maximum values of the inertia weighting factor, respectively.


### 3.2. DMOPSO Algorithm

DMOPSO algorithm is to optimize the function with multiple objective functions in the dynamic environment, while optimizing the objective functions to achieve the best condition. In the dynamic environment, the optimal solution of the individual and the optimal solution of the population may change over time, and the particle is easy to fall into the optimization state of the former environment and stagnate. Therefore, the core of the dynamic multiobjective optimization problem is to adapt to environmental changes and can detect changes in the environment quickly and accurately and how to make adjustments to adapt to environmental changes. In this paper, the dynamic environment monitoring method based on sentinel particles [20] is used. The basic idea of this method is to generate a part of the sentinel particles randomly as the population is initialized at the same time. The fitness values of the sentinel particles are calculated in each iteration. When the fitness values of the sentinel particles change, it is determined that the current environment has changed. The sentinel particles do not participate in each iteration, so the fitness values of the sentinel particles should remain constant in the static environment. When monitoring the environment changes, the DMOPSO algorithm needs to adjust the current optimization results. Firstly, the information in the nondominated solution memory is recalculated under the current environmental conditions, the dominated solution is removed, and then a part of the particle location in the search space is reinitialized. Generally, thirty percent of the particles are selected. The specific flow chart of the dynamic MOPSO algorithm based on sentinel particle monitoring is depicted in Figure 3.Figure 3
The flow chart of DMOPSO algorithm based on sentinel particle monitoring.Based on the process of DMOPSO algorithm, formulae(1) and (3) are used as the objective functions, and then the calculation method of the route planning safety area minimum width in Section 3 is combined. So the dynamic MACSIA array schemes can be obtained in the case of two kinds of battlefield environment dynamic changes, which satisfy the limitation of the given suppression interference effect and ensure the safety of the jamming aircraft themselves.


## 3.1. Particle Swarm Optimization

Particle swarm optimization (PSO) [26] algorithm is widely used in solving single target optimization problem. The updating location and speed expressions of particles in PSO are shown in (7)Xi+1=Xi+Vi+1Vi+1=wVi+c1r1Pbesti-Xi+c2r2Gbesti-XiIn this equation, Vi indicates the current speed of the particle, Vi+1 represents the speed at the next moment, Xi indicates the current location of the particle, the position of the particle at the next moment is denoted by Xi+1, Gbesti is the global optimal solution for the i-th generation of population, Pbesti represents the individual optimal solution for the i-th generation of particle, and r1 and r2 are random numbers with values between 0 and 1. Both c1 and c2 represent learning factors, w indicates the inertia weight factor, which is a variable that reflects the effect of the current velocity of the particle on the velocity of the particle at the next moment, and the expression is (8)w=wmax-wmax-wminiteramax∗iteraHere, itera is the current number of population iterations, iteramax represents the maximum number of population iterations, and wmin and wmax represent the minimum and maximum values of the inertia weighting factor, respectively.


## 3.2. DMOPSO Algorithm

DMOPSO algorithm is to optimize the function with multiple objective functions in the dynamic environment, while optimizing the objective functions to achieve the best condition. In the dynamic environment, the optimal solution of the individual and the optimal solution of the population may change over time, and the particle is easy to fall into the optimization state of the former environment and stagnate. Therefore, the core of the dynamic multiobjective optimization problem is to adapt to environmental changes and can detect changes in the environment quickly and accurately and how to make adjustments to adapt to environmental changes. In this paper, the dynamic environment monitoring method based on sentinel particles [20] is used. The basic idea of this method is to generate a part of the sentinel particles randomly as the population is initialized at the same time. The fitness values of the sentinel particles are calculated in each iteration. When the fitness values of the sentinel particles change, it is determined that the current environment has changed. The sentinel particles do not participate in each iteration, so the fitness values of the sentinel particles should remain constant in the static environment. When monitoring the environment changes, the DMOPSO algorithm needs to adjust the current optimization results. Firstly, the information in the nondominated solution memory is recalculated under the current environmental conditions, the dominated solution is removed, and then a part of the particle location in the search space is reinitialized. Generally, thirty percent of the particles are selected. The specific flow chart of the dynamic MOPSO algorithm based on sentinel particle monitoring is depicted in Figure 3.Figure 3
The flow chart of DMOPSO algorithm based on sentinel particle monitoring.Based on the process of DMOPSO algorithm, formulae(1) and (3) are used as the objective functions, and then the calculation method of the route planning safety area minimum width in Section 3 is combined. So the dynamic MACSIA array schemes can be obtained in the case of two kinds of battlefield environment dynamic changes, which satisfy the limitation of the given suppression interference effect and ensure the safety of the jamming aircraft themselves.


## 4. Simulation Experiment

Simulation environment is in Windows 7 32-bit system, and the processor is the Intel (R) Core (TM) i5-4590 CPU @ 3.3GHz. The programming language is programmed with MATLAB 2010a. A432km∗432km size area is selected as the simulation experiment of the combat scene. The digital elevation of this region is shown in Figure 4, it uses a grid structure, and the elevation data resolution is 360m.Figure 4
Digital elevation terrain map.The first case: the location of the enemy air defense radar is mobile, but the number remains the same.Assuming that four enemy air defense radars are deployed in the combat area, the minimum suppression factorsKj of the four enemy radars are 5, as shown in Table 1 for the list of performance parameters of the four enemy radars. The initial coordinate locations of the four enemy radars are (240km,100km), (150km,120km), (300km,250km), and (180km,280km).Table 1
The list of four enemy air defense radar performance parameters.
Radar P t / k W G λ / m L F n T o / K B n / H z 1# 4000 50 0.6 2.5 3 291 2 × 10 5 2# 3900 45 0.8 4 2.5 291 2.1 × 10 3# 3100 70 0.75 3 4 291 3 × 10 5 4# 2600 100 0.5 3 6 291 1.5 × 10Initially, the detection range of the four enemy air defense radars under terrain shading at a height of2.1km is shown in Figure 5. In the figure, the detection range of the four enemy air defense radars overlaps, so it is very difficult for our combat aircraft to penetrate the defense safely. At this time the combat aircraft must rely on the electronic support jamming aircraft to implement cooperative suppression interference against the enemy air defense radar network, so as to obtain a given limit to meet the conditions of the route planning safety area width. The minimum width of the route planning safety area is limited to 20km. We intend to send three electronic support jamming aircraft to implement cooperative suppression interference against the enemy air defense radar network. The three jamming aircraft are at a height of 2.1km and the performance parameters are shown in Table 2.Table 2
The list of three jamming aircrafts parameters.
Jamming aircraft P j / k W G j r j L j B j / H z 1# 1 × 10 3 4 2 20 2 × 10 6 2# 2 × 10 3 2 2 20 2 × 10 6 3# 2 × 10 3 6 1 30 2 × 10 6Figure 5
The detection range of the four enemy air defense radars under terrain shading at the beginning.The initial parameters of DMOPSO algorithm are as follows: the population sizexSize is 100, the dimension of decision variables Dim is 6, the learning factor c1=c2=0.75, the maximum and minimum values of the inertia weighting factor are wmax=1.1, wmin=0.2, respectively, and the maximum number of iterations Max_itera is 500. The number of environment changes is set to T=3, the number of iterations of DMOPSO algorithm is 150 for each time the environment changes, and the threshold of environment change is δ=0.5; set the initial environment t=0. As the moving law of the enemy air defense radars is difficult to obtain, therefore, in this paper, for the sake of simplified study, the coordinates (in km) of the four enemy air defense radars after three environment changes are shown in Table 3. Here is the detection range of the four air defense radars under terrain shading after the third environment change, as shown in Figure 6.Table 3
The coordinates of the four enemy air defense radars after three environment changes.
Radar After the first environment change After the Second environment change After the third environment change 1# (240, 100) (230, 100) (240, 105) 2# (150, 120) (160, 130) (160, 120) 3# (300, 250) (300, 250) (300, 220) 4# (170, 270) (170, 270) (180, 250)Figure 6
The detection range of the four air defense radars under terrain shading after the third environment change.The experimental results were run independently by DMOPSO algorithm 30 times, and the optimal result was selected from them. The calculated nondominated solutions distribution is shown in Figure7.Figure 7
The obtained nondominated solutions distribution in the first case.It can be seen from Figure7 that DMOPSO algorithm has a uniform distribution of nondominated solution sets, which constitute four Pareto front surfaces, and it is concluded that the search effect of DMOPSO algorithm is still good when the environment changes. In this paper, the Pareto front surface generated after the third environmental change is taken as an example. In order to facilitate the decision-maker to make decisions, three points are selected from the Pareto front surface: the intermediate point S2 and the two endpoints S1 and S3, respectively, where S1 is biased to optimize the f2 objective and S3 is biased towards the optimized f1 objective, while S2 takes into account both objective f2 and f1. The optimal dynamic MACSIA schemes of the three points are shown in Figures 8, 9, and 10, respectively.Figure 8
The optimal dynamic MACSIA scheme of pointS1 in the first case.Figure 9
The optimal dynamic MACSIA scheme of pointS3 in the first case.Figure 10
The optimal dynamic MACSIA scheme of pointS2 in the first case.When the decision-maker who develops the combat plan needs to minimize the sum of the distance between each jamming aircraft and the enemy air defense radar network center and maximize the minimum width of the route planning safety area, the emphasis is the suppression interference effect on the enemy air defense radar network. And the optimal dynamic MACSIA scheme of pointS1 can be chosen; as shown in Figure 8, the three jamming aircraft are closer to the enemy air defense radar network center in this scheme. When the decision-maker who develops the combat plan needs to maximize the sum of the distance between each jamming aircraft and the enemy air defense radar network center and minimize the minimum width of the route planning safety area, that is to say, the focus is on the safety of each jamming aircraft itself, the optimal dynamic MACSIA scheme of point S3 can be chosen; as shown in Figure 9, the three jamming aircraft are more far away from the enemy air defense radar network center in this scheme. When the decision-maker who develops the combat plan not only needs to consider the safety of each jamming aircraft itself, but also needs to achieve a given suppression interference effect, the optimal dynamic MACSIA scheme of point S2 may be the best choice, as it can be seen in Figure 10. At this time the total distance between each jamming aircraft and the enemy air defense radar network center is moderate in this scheme.Figures11, 12, and 13 are the binary images of suppression interference effect for the jamming aircraft against the enemy air defense radar network in three schemes. The obtained minimum safety widths are 45km, 28km, and 37.5km, respectively.Figure 11
The binary image of suppression interference effect of pointS1 for the jamming aircraft against the enemy air defense radar network in the first case.Figure 12
The binary image of suppression interference effect of pointS3 for the jamming aircraft against the enemy air defense radar network in the first case.Figure 13
The binary image of suppression interference effect of pointS2 for the jamming aircraft against the enemy air defense radar network in the first case.The second case: the number of the enemy air defense radar is variable, but the original location remains unchanged.At first, assuming that three enemy air defense radars are deployed in the combat area, the minimum suppression factorsKj of the three enemy air defense radars are 5. The performance parameters are the same as the number 1, 2, and 3 radar in Table 1. And the initial coordinate locations of the three enemy air defense radars are (220km,105km), (290km,230km), and (190km,270km).Initially, the detection range of the three enemy air defense radars under terrain shading at a height of2.1km is shown in Figure 14. It can be seen from Figure 13 that the detection range of the three enemy air defense radars overlaps, so it is very difficult for our combat aircraft to penetrate the defense safely. At this time the combat aircraft must rely on the electronic support jamming aircraft to implement cooperative suppression interference against the enemy air defense radar network, so as to obtain a given limit to meet the conditions of the route planning safety area width. The minimum width of the route planning safety area is limited to 20km. At the beginning, we intend to send two electronic support jamming aircraft to implement cooperative suppression interference against the enemy air defense radar network. The two jamming aircraft are at a height of 2.1km and their performance parameters are the same as the number 1 and 2 jamming aircraft in Table 2.Figure 14
The detection range of the three enemy air defense radars under terrain shading at the beginning.The initial parameters of DMOPSO algorithm are as follows: the maximum number of iterationsMax_itera is 800, the number of environment changes is set to T=2, the number of iterations of DMOPSO algorithm is 250 for each time the environment changes, the threshold of environment change is δ=0.6, and the remaining parameters are consistent with the simulation parameters in the first case. In order to simplify the study, at the same time more in line with the real battlefield environment, the number of enemy air defense radar is limited to 3 to 5, and the number of enemy air defense radars is increased in the process of environment change. The coordinates of the additional enemy air defense radar after the two environment changes are (150km,210km) and (300km,130km), respectively. Here is given the detection range of the enemy air defense radar network by terrain shading after the second environment change, which is shown in Figure 15.Figure 15
The detection range of the five air defense radars under terrain shading after the second environment change.The experimental results were run independently by DMOPSO algorithm 30 times, and the optimal result was selected from them. The calculated nondominated solutions distribution is shown in Figure16.Figure 16
The obtained nondominated solutions distribution in the second case.As can be seen from Figure16, when the environment changes, the search effect of DMOPSO algorithm is still relatively good, and three more uniform Pareto front surfaces are generated. Then the Pareto front surface generated after the second environment change is taken as an example. Similarly, in order to facilitate decision-maker to make decisions, three points called S1, S2, and S3 are selected on the Pareto frontier surface. The optimal dynamic MACSIA schemes of the three points are shown in Figures 17, 18, and 19, respectively.Figure 17
The optimal dynamic MACSIA scheme of pointS1 in the second case.Figure 18
The optimal dynamic MACSIA scheme of pointS3 in the second case.Figure 19
The optimal dynamic MACSIA scheme of pointS2 in the second case.As shown in Figure17, each jamming aircraft is relatively close to the enemy air defense radar network center in this scheme, the total distance is the smallest and the minimum width of the route planning safe area is the largest. The choice of point S1 corresponding to the scheme means that more emphasis is on the suppression interference effect on the enemy air defense radar network. Figure 18 shows that the four jamming aircraft are far from the enemy air defense radar network center in this scheme, correspondingly, the sum of the distance between each jamming aircraft and the enemy air defense radar network center is the largest, and the minimum width of the route planning safe area is the smallest. If the focus is on the safety of the jamming aircraft themselves, the point S3 corresponding to the scheme can be selected. As shown in Figure 19, at this time, the total distance between the four jamming aircraft and the enemy air defense radar network center is moderate. When the decision-maker who develops the combat plan needs to consider the safety of the jamming aircraft themselves and achieve the given suppression interference effect, the point S2 corresponding to the scheme can be selected.Figures20, 21, and 22 are the binary images of suppression interference effect for the jamming aircraft against the enemy air defense radar network in three schemes. The obtained minimum safety widths are 32km, 22km, and 27.5km, respectively.Figure 20
The binary image of suppression interference effect of pointS1 for the jamming aircraft against the enemy air defense radar network in the second case.Figure 21
The binary image of suppression interference effect of pointS3 for the jamming aircraft against the enemy air defense radar network in the second case.Figure 22
The binary image of suppression interference effect of pointS2 for the jamming aircraft against the enemy air defense radar network in the second case.It can be concluded from the above that the six kinds of dynamic MACSIA schemes have achieved good suppression interference effect in the two cases and formed the route planning safety area which satisfies the restriction condition and ensures the safety of the jamming aircraft themselves. The decision-maker who develops the combat plan can select the corresponding scheme according to the actual operational needs.


## 5. Conclusions

In this paper, aiming at the problem of dynamic MACSIA optimization in electronic warfare mission planning, based on two different environment changes, two kinds of dynamic multiobjective optimization models of the dynamic MACSIA are constructed by regarding the sum of the distance between each jamming aircraft and the enemy air defense radar network center and the minimum width of the route planning safety area as the objective functions. The DMOPSO algorithm is used to solve the two models, and the optimal dynamic MACSIA schemes are calculated in the case of two environment changes. Simultaneously, the feasibility and validity of the proposed model and the proposed method are verified. The study content of this paper is a key part of the dynamic planning in electronic warfare mission planning. When the battlefield environment changes, it can provide a basis for decision-maker to make decisions and also provide a safe and reliable planning space for the follow-up route planning of the combat aircraft. Therefore, this paper has a relatively strong significance of actual combat.



---

*Source: 4521701-2018-07-16.xml*