# Rapid Penetration Path Planning Method for Stealth UAV in Complex Environment with BB Threats

Authors: Zhe Zhang; Jian Wu; Jiyang Dai; Cheng He
Journal: International Journal of Aerospace Engineering
 (2020)

---

## Abstract

This paper presents the flight penetration path planning algorithm in a complex environment with Bogie or Bandit (BB) threats for stealth unmanned aerial vehicle (UAV). The emergence of rigorous air defense radar net necessitates efficient flight path planning and replanning for stealth UAV concerning survivability and penetration ability. We propose the improved A-Star algorithm based on the multiple step search approach to deal with this uprising problem. The objective is to achieve rapid penetration path planning for stealth UAV in a complex environment. Firstly, the combination of single-base radar, dual-base radar, and BB threats is adopted to different threat scenarios which are closer to the real combat environment. Besides, the multistep search strategy, the prediction technique, and path planning algorithm are developed for stealth UAV to deal with BB threats and achieve the penetration path replanning in complex scenarios. Moreover, the attitude angle information is integrated into the flight path which can meet real flight requirements for stealth UAV. The theoretical analysis and numerical results prove the validity of our method.


---

## Body

## 1. Introduction

### 1.1. Background and Motivation

In recent years, the emergence of new military equipment and systems, the stealth unmanned aerial vehicles (UAV) are gaining much attention due to the air-defense radar net is increasingly in the field of close and low altitude combat. Meanwhile, the survivability and combat effectiveness of stealth UAV will have a great influence on modern air warfare [1, 2]. The nose, fuselage, and tail of stealth UAV have relatively small radar cross-section (RCS), which aims to reduce the detection probability of radar system and improve the penetration ability and survivability of the stealth UAV [3]. However, previous works have mainly concentrated on the path planning in the static environment, there are many theoretical challenges and practical problems for stealth UAV on the problem of penetration path planning, such as the constraint of attitude angle and control of track point, the efficiency of path planning algorithm, and real-time performance. Besides, with regard to the Bogie or Bandit (BB) threats in the real combat environment, it is difficult for stealth UAV to achieve penetration path planning which is rapidly and safely. Therefore, the kinematics model and path planning algorithm for stealth UAV should be studied in a complex environment.


### 1.2. Literature Review

More recently, a lot of research has been done in the area of penetration path planning for aircraft, and so many algorithms had been proposed. In Reference [4], a new concept on the aircraft path planning for the variation of radar cross-section (RCS) was firstly presented, which indicates that the stealth UAV can significantly reduce the RCS by adjusting the attitude angle. In Reference [5], a research framework of low detectability tactical trajectory planning was proposed in the presence of multiple radar-detected threat environments; the solution of the optimal trajectory is transformed into the optimal control problem with the minimum radar detection probability. A penetration trajectory optimization method considering the influence of aircraft radar RCS was proposed in Reference [6]. In Reference [7] and Reference [8], a hybrid heuristic adaptive pseudospectrum method was proposed to address the low-detectability trajectory planning problem. In contrast, these methods usually adopt a simplified kinematic model, which aims to obtain a rough reference track, and the constraints and control conditions of radar attitude are rarely considered in the reference trajectories. In Reference [9], a path planning method for stealth UAV based on the RCS ellipsoid model was proposed to find the optimal trajectory during flight. Inversely, this model requires a large amount of computation in addressing the problem of path planning. In Reference [10], a trajectory planning method based on dynamic game theory was presented to analyze the relationship between radar and UAV. In Reference [11], a real-time air battle trajectory optimization and game model for aircraft based on a rolling time-domain control strategy were proposed in a complex environment. However, two kinds of shortages exist in the research. One was that only can be applied to the static environment with known radar position, and it rarely achieves the path planning for stealth UAV in the dynamic environment. The other was that the calculation process is relatively complicated, and the real-time performance of the path planning method is poor. Erlandsson proposed a path planning model, and the expected path cost of air mission was discussed in a hostile environment; the particle swarm optimization algorithm (PSO) was adopted to address the optimization path planning problem [12, 13]. Meanwhile, many scholars are focus on analyzing a series of optimization algorithms to address the penetration path planning problem for unmanned aerial vehicle, such as A-Star algorithm [14], genetic algorithms (GA) [15], sparse A-Star algorithm (SAS) [16], rapidly exploring random-tree algorithm (RRT) [17, 18], particle swarm optimization algorithm (PSO) [19, 20], black hole algorithm [21], and other combinatorial optimization algorithms [22]. However, these optimization algorithms are only applied to address the stealth path planning problem in the static combat environment, the single-base radar is usually adopted in the combat scenario that the position of each radar is known, and the number is relatively small. Furthermore, the efficiency of the algorithm and the safety of the path are rarely optimal.


### 1.3. Contribution

In this paper, we focus on analyzing the penetration path planning for stealth UAV based on an improved A-Star algorithm, which aims to achieve rapid penetration path planning in the dynamic combat environment with BB threats. The novelty of this method is summarized as follows: firstly, the main idea of the model-based predictive control (MPC) and learning real-time A-Star algorithm (LRTA) is integrated into the path to devise the improved A-Star algorithm. Additionally, for the BB threats in a complex environment, the improved A-Star algorithm is applied for the penetration path replanning with fixed altitude to improve the survivability and combat effectiveness of stealth UAV. Moreover, numerical simulation results are applied to demonstrate the validity of the improved A-Star algorithm in the presence of the radar net with BB threats.


### 1.4. Paper Organization

The paper is organized as follows. Section2 describes the mathematical model of stealth UAV. Section 3 discusses the detection probability calculation of multiple radar net, which is close to the real combat environment. Section 4 presents the improvement of the path planning algorithm, includes improved A-Star algorithm, LRTA-Star algorithm, and D-Star algorithm. Our numerical results and evaluation are presented in Section 5. Section 6 states our conclusions and future work.


## 1.1. Background and Motivation

In recent years, the emergence of new military equipment and systems, the stealth unmanned aerial vehicles (UAV) are gaining much attention due to the air-defense radar net is increasingly in the field of close and low altitude combat. Meanwhile, the survivability and combat effectiveness of stealth UAV will have a great influence on modern air warfare [1, 2]. The nose, fuselage, and tail of stealth UAV have relatively small radar cross-section (RCS), which aims to reduce the detection probability of radar system and improve the penetration ability and survivability of the stealth UAV [3]. However, previous works have mainly concentrated on the path planning in the static environment, there are many theoretical challenges and practical problems for stealth UAV on the problem of penetration path planning, such as the constraint of attitude angle and control of track point, the efficiency of path planning algorithm, and real-time performance. Besides, with regard to the Bogie or Bandit (BB) threats in the real combat environment, it is difficult for stealth UAV to achieve penetration path planning which is rapidly and safely. Therefore, the kinematics model and path planning algorithm for stealth UAV should be studied in a complex environment.


## 1.2. Literature Review

More recently, a lot of research has been done in the area of penetration path planning for aircraft, and so many algorithms had been proposed. In Reference [4], a new concept on the aircraft path planning for the variation of radar cross-section (RCS) was firstly presented, which indicates that the stealth UAV can significantly reduce the RCS by adjusting the attitude angle. In Reference [5], a research framework of low detectability tactical trajectory planning was proposed in the presence of multiple radar-detected threat environments; the solution of the optimal trajectory is transformed into the optimal control problem with the minimum radar detection probability. A penetration trajectory optimization method considering the influence of aircraft radar RCS was proposed in Reference [6]. In Reference [7] and Reference [8], a hybrid heuristic adaptive pseudospectrum method was proposed to address the low-detectability trajectory planning problem. In contrast, these methods usually adopt a simplified kinematic model, which aims to obtain a rough reference track, and the constraints and control conditions of radar attitude are rarely considered in the reference trajectories. In Reference [9], a path planning method for stealth UAV based on the RCS ellipsoid model was proposed to find the optimal trajectory during flight. Inversely, this model requires a large amount of computation in addressing the problem of path planning. In Reference [10], a trajectory planning method based on dynamic game theory was presented to analyze the relationship between radar and UAV. In Reference [11], a real-time air battle trajectory optimization and game model for aircraft based on a rolling time-domain control strategy were proposed in a complex environment. However, two kinds of shortages exist in the research. One was that only can be applied to the static environment with known radar position, and it rarely achieves the path planning for stealth UAV in the dynamic environment. The other was that the calculation process is relatively complicated, and the real-time performance of the path planning method is poor. Erlandsson proposed a path planning model, and the expected path cost of air mission was discussed in a hostile environment; the particle swarm optimization algorithm (PSO) was adopted to address the optimization path planning problem [12, 13]. Meanwhile, many scholars are focus on analyzing a series of optimization algorithms to address the penetration path planning problem for unmanned aerial vehicle, such as A-Star algorithm [14], genetic algorithms (GA) [15], sparse A-Star algorithm (SAS) [16], rapidly exploring random-tree algorithm (RRT) [17, 18], particle swarm optimization algorithm (PSO) [19, 20], black hole algorithm [21], and other combinatorial optimization algorithms [22]. However, these optimization algorithms are only applied to address the stealth path planning problem in the static combat environment, the single-base radar is usually adopted in the combat scenario that the position of each radar is known, and the number is relatively small. Furthermore, the efficiency of the algorithm and the safety of the path are rarely optimal.


## 1.3. Contribution

In this paper, we focus on analyzing the penetration path planning for stealth UAV based on an improved A-Star algorithm, which aims to achieve rapid penetration path planning in the dynamic combat environment with BB threats. The novelty of this method is summarized as follows: firstly, the main idea of the model-based predictive control (MPC) and learning real-time A-Star algorithm (LRTA) is integrated into the path to devise the improved A-Star algorithm. Additionally, for the BB threats in a complex environment, the improved A-Star algorithm is applied for the penetration path replanning with fixed altitude to improve the survivability and combat effectiveness of stealth UAV. Moreover, numerical simulation results are applied to demonstrate the validity of the improved A-Star algorithm in the presence of the radar net with BB threats.


## 1.4. Paper Organization

The paper is organized as follows. Section2 describes the mathematical model of stealth UAV. Section 3 discusses the detection probability calculation of multiple radar net, which is close to the real combat environment. Section 4 presents the improvement of the path planning algorithm, includes improved A-Star algorithm, LRTA-Star algorithm, and D-Star algorithm. Our numerical results and evaluation are presented in Section 5. Section 6 states our conclusions and future work.


## 2. The System Modeling

### 2.1. Kinematics

It is assumed that the stealth UAV can be defined as a particle with attitude information and ignores the influence of wind and other external conditions [1]. We are focus on studying the control of the flight path. The stealth UAV moves in a horizontal plane at a constant altitude h. Therefore, the kinematics model can be given by
(1)x•=vcosφy•=vsinφφ•=uv,where x and y are the Cartesian coordinates of the aircraft, φ is the heading angle, v is the constant speed, u is the input signal, and is the acceleration normal to the flight path vector.Specify radar locations in the planeh=0 (on the ground) by the coordinates xi,yi,0. The range from the ith radar to the stealth UAV is given by
(2)Ri=x−xi2+y−yi2+h2.Additionally, let(3)αi=arctany−yix−xiλi=αi−φ+πθ=arctanug,be the azimuth, aspect, and roll angles, respectively, measured to the ith radar, g is the acceleration of gravity.


### 2.2. Dynamic RCS Features

Radar cross-section (RCS) is a measure of ability for UAV to reflect electromagnetic radiation emitted by radar [23]; the RCS of the nose, fuselage, and tail of the stealth UAV is smaller than conventional UAV, and the RCS value will be changed due to different radar bands and radar polarization modes. It plays an important role in the penetration path of stealth UAV in combat. For the path planning problem, RCS is usually defined as a constant, but it is not reasonable in practical application. By convention [24], it is modeled as a function of the aspect angle and roll angle as viewed from the radar. Therefore, the dynamic RCS model is adopted to an ellipsoid [1], which is given by
(4)σλ,θ=πa2b2c2a2sin2λcos2θ+b2sin2λsin2θ+c2cos2λ2,where σ represents the RCS value of UAV, λ represents aspect angle, θ represents roll angle, a is relatively small frontal RCS, b is a larger beam aspect RCS, and c is relatively large RCS when viewed from above or below.


## 2.1. Kinematics

It is assumed that the stealth UAV can be defined as a particle with attitude information and ignores the influence of wind and other external conditions [1]. We are focus on studying the control of the flight path. The stealth UAV moves in a horizontal plane at a constant altitude h. Therefore, the kinematics model can be given by
(1)x•=vcosφy•=vsinφφ•=uv,where x and y are the Cartesian coordinates of the aircraft, φ is the heading angle, v is the constant speed, u is the input signal, and is the acceleration normal to the flight path vector.Specify radar locations in the planeh=0 (on the ground) by the coordinates xi,yi,0. The range from the ith radar to the stealth UAV is given by
(2)Ri=x−xi2+y−yi2+h2.Additionally, let(3)αi=arctany−yix−xiλi=αi−φ+πθ=arctanug,be the azimuth, aspect, and roll angles, respectively, measured to the ith radar, g is the acceleration of gravity.


## 2.2. Dynamic RCS Features

Radar cross-section (RCS) is a measure of ability for UAV to reflect electromagnetic radiation emitted by radar [23]; the RCS of the nose, fuselage, and tail of the stealth UAV is smaller than conventional UAV, and the RCS value will be changed due to different radar bands and radar polarization modes. It plays an important role in the penetration path of stealth UAV in combat. For the path planning problem, RCS is usually defined as a constant, but it is not reasonable in practical application. By convention [24], it is modeled as a function of the aspect angle and roll angle as viewed from the radar. Therefore, the dynamic RCS model is adopted to an ellipsoid [1], which is given by
(4)σλ,θ=πa2b2c2a2sin2λcos2θ+b2sin2λsin2θ+c2cos2λ2,where σ represents the RCS value of UAV, λ represents aspect angle, θ represents roll angle, a is relatively small frontal RCS, b is a larger beam aspect RCS, and c is relatively large RCS when viewed from above or below.


## 3. Radar Detecting Probability

When the acquired target information meets the common track criteria, the radar system confirms the detection of the target and transmits the acquired target information to the information fusion center. The detection probability of target in a single-base radar, dual-base radar, and a radar net system is depicted as follows.


### 3.1. Single-Base Radar

For the single-base radar, the detection probability of the UAV is only related to the distance from the UAV to the radar center. In a period, the radar instantaneous detection probability can be expressed as follows:(5)Pt=11+c2Ri4/σc1,where Pt is the instantaneous detection probability of the single base radar, Ri is given by Eq.(2), and c1 and c2 are the performance parameters of specific radar, respectively.


### 3.2. Dual-Base Radar

The dual-base radar system is considered to be effective systems to counter the four major threats due to their technical characteristics of transmitter-receiver separation and passive reception. For the performance of the radar system, it has the advantages that single-base radar can not match, such as a larger detection range and higher accuracy. Dual-base radars are usually deployed in modern air defense systems; therefore, we focus on analyzing the dual-base radar system.The transmitter is also employed as a receiver for dual-base radar systems, which are T/R-R systems. A transmitter can work synchronously with other splitter receivers in the dual-base radar system. The receiver synchronizes with the transmitter in time, phase, and space. The configuration of the dual-base radar system is depicted in Figure1.Figure 1
The configuration of the dual-base radar system.wherePx,y represents the position of the target, T represents the transmitter, R represents the receiver, Db represents the distance between the transmitter and receiver, Dr represents the distance between the target and transmitter, Dt represents the distance between the target and receiver, and βd represents the dual-base angle.Additionally,D is defined as the detection range of the bistatic radar, that is D=DtDr. where Dt=x2+y2 and Dr=x−Db2+y2. D is a constant when the radar parameters are given. Hence, the maximum detection range of the bistatic radar Dmax can be obtained by
(6)Dmax=Pτσλc2G4π3KT0LS/N1/4,where Pτ represents the peak power of the transmitter, G represents the antenna gain, L represents the system loss, K represents Boltzmann constant, T0 is the noise temperature of the receiver, S/N is the signal to noise ratio of the receiver, and λc represents the wavelength. Similarly, the detection probability of the bistatic radar under exponential distribution is given by
(7)PD=e−D2/D2max0≤D≤Dmax0D>Dmax.


### 3.3. Detection Probability of Networked Radar System

For a complete netted radar system, the radar net can improve the detection probability of UAV, significantly. The detection probability of a networked radar system mainly refers to the target detection probability calculated by the information fusion center. The rank K fusion rule is widely adopted in modern networked radar systems. Therefore, we adopt the rank K fusion rule to analyze the detection probability of networked radar [25, 26].It is assumed that the number of radars in the netted radar system isM, which includes single-base radar and bistatic radar. When the number of radar detected in the system exceeds the detection threshold according to the rank K fusion rule, that means, the target has been detected by the radar system, and the approximate value of the optimal detection threshold K0 can be depicted as follows.
(8)K0=1.5M.The radar makes a local judgment based on the detection of the stealth UAV, and the judgment results are either 0 or 1, which depends on whether the local threshold detection target exists or not.H0 and H1 are binary assumptions, where H0 represents target does not exist and H1 represents target exists. Therefore, the decision value of the jth radar j=1,2,⋯,M) can be expressed as follows:
(9)dj=0ifdecisionresultisH01ifdecisionresultisH1.Additionally, the local decision result is passed to the information fusion center of the radar system to form a global decision matrixDc, that is, Dc=d1,d2,⋯,dM. And the radar information fusion rule for the network is denoted as R; the decision rule RDc for rank K fusion is given by
(10)RDc=1ifthe∑j=1Mdj≥K00ifthe∑j=1Mdj<K0.Meanwhile, the total detection probability of the network radar system on the UAV is given by(11)PNet=∑DcRDc∏dj∈S1Pdj∏dj∈S01−Pdj,where S1 is a set of local decision vectors that make the fusion center judge “1”, S0 is a set of local decision vectors that make the fusion center judge “0”, and Pdj is the discovery probability of the first radar in the radar net.


## 3.1. Single-Base Radar

For the single-base radar, the detection probability of the UAV is only related to the distance from the UAV to the radar center. In a period, the radar instantaneous detection probability can be expressed as follows:(5)Pt=11+c2Ri4/σc1,where Pt is the instantaneous detection probability of the single base radar, Ri is given by Eq.(2), and c1 and c2 are the performance parameters of specific radar, respectively.


## 3.2. Dual-Base Radar

The dual-base radar system is considered to be effective systems to counter the four major threats due to their technical characteristics of transmitter-receiver separation and passive reception. For the performance of the radar system, it has the advantages that single-base radar can not match, such as a larger detection range and higher accuracy. Dual-base radars are usually deployed in modern air defense systems; therefore, we focus on analyzing the dual-base radar system.The transmitter is also employed as a receiver for dual-base radar systems, which are T/R-R systems. A transmitter can work synchronously with other splitter receivers in the dual-base radar system. The receiver synchronizes with the transmitter in time, phase, and space. The configuration of the dual-base radar system is depicted in Figure1.Figure 1
The configuration of the dual-base radar system.wherePx,y represents the position of the target, T represents the transmitter, R represents the receiver, Db represents the distance between the transmitter and receiver, Dr represents the distance between the target and transmitter, Dt represents the distance between the target and receiver, and βd represents the dual-base angle.Additionally,D is defined as the detection range of the bistatic radar, that is D=DtDr. where Dt=x2+y2 and Dr=x−Db2+y2. D is a constant when the radar parameters are given. Hence, the maximum detection range of the bistatic radar Dmax can be obtained by
(6)Dmax=Pτσλc2G4π3KT0LS/N1/4,where Pτ represents the peak power of the transmitter, G represents the antenna gain, L represents the system loss, K represents Boltzmann constant, T0 is the noise temperature of the receiver, S/N is the signal to noise ratio of the receiver, and λc represents the wavelength. Similarly, the detection probability of the bistatic radar under exponential distribution is given by
(7)PD=e−D2/D2max0≤D≤Dmax0D>Dmax.


## 3.3. Detection Probability of Networked Radar System

For a complete netted radar system, the radar net can improve the detection probability of UAV, significantly. The detection probability of a networked radar system mainly refers to the target detection probability calculated by the information fusion center. The rank K fusion rule is widely adopted in modern networked radar systems. Therefore, we adopt the rank K fusion rule to analyze the detection probability of networked radar [25, 26].It is assumed that the number of radars in the netted radar system isM, which includes single-base radar and bistatic radar. When the number of radar detected in the system exceeds the detection threshold according to the rank K fusion rule, that means, the target has been detected by the radar system, and the approximate value of the optimal detection threshold K0 can be depicted as follows.
(8)K0=1.5M.The radar makes a local judgment based on the detection of the stealth UAV, and the judgment results are either 0 or 1, which depends on whether the local threshold detection target exists or not.H0 and H1 are binary assumptions, where H0 represents target does not exist and H1 represents target exists. Therefore, the decision value of the jth radar j=1,2,⋯,M) can be expressed as follows:
(9)dj=0ifdecisionresultisH01ifdecisionresultisH1.Additionally, the local decision result is passed to the information fusion center of the radar system to form a global decision matrixDc, that is, Dc=d1,d2,⋯,dM. And the radar information fusion rule for the network is denoted as R; the decision rule RDc for rank K fusion is given by
(10)RDc=1ifthe∑j=1Mdj≥K00ifthe∑j=1Mdj<K0.Meanwhile, the total detection probability of the network radar system on the UAV is given by(11)PNet=∑DcRDc∏dj∈S1Pdj∏dj∈S01−Pdj,where S1 is a set of local decision vectors that make the fusion center judge “1”, S0 is a set of local decision vectors that make the fusion center judge “0”, and Pdj is the discovery probability of the first radar in the radar net.


## 4. Path Planning Algorithm

A-star algorithm is an effective search method to solve the shortest path in a static road network, which is widely used to settle the path planning problems of many agents. The closer the estimated distance is to the actual value, the faster the final search speed is [27–29]. This paper discusses the improvement method of the A-Star algorithm to address the penetration path planning problem for stealth UAV.


### 4.1. The A-Star Algorithm

The main idea of the standard A-Star algorithm is as follows. Firstly, select the appropriate heuristic function, estimate the generation value of the extensible search points in the search area, comprehensively. Moreover, compare the different cost values of each point. Moreover, consider the operation time and distance cost of the track point search. Finally, find an optimal path. In the A-Star algorithm, the operation of the OPEN list and CLOSE list are usually performed to achieve the storage and update of track points.However, the standard A-Star algorithm that is applied to find the optimal path for stealth UAV has many disadvantages: (a) the route obtained by the A-star algorithm only has the position of stealth UAV, and it can not reflect the characteristics of the dynamic RCS and attitude angle information. (b) Unknown path cost estimation is rarely calculated in the process of path planning accurately, and the computation time is too long to find a globally optimal path. (c) the performance of real-time is poor, and the algorithm can not deal with BB threats. Therefore, we are focus on improving the performance of the A-Star algorithm.


### 4.2. The D-Star Algorithm

The D-Star algorithm is developed by the A-Star algorithm and the Dijkstra algorithm [30, 31], which is suitable for addressing the path planning problem in unknown environments. The main idea of the algorithm is to search the reverse path from the destination to the start, and the heuristic function expression of the D-Star algorithm is given by
(12)fX,E=hX+gX,E,where hX represents the actual cost of the path from the destination to state X, and gX,E is the estimated cost of the path from the state X to the starting point.The main steps of the D-Star algorithm are given as follows. Step 1: the value of all states are set to NEW,f∗ and h∗ are set to infinity for all states, hG is set to zero, and the destination G is added to the OPEN list. Step 2: keep performing f∗ the search until the location of E which is removed from the OPEN list, if the CLOSE list contains the tag value of the state X, the complete path sequence will be obtained. In contrast, path planning has failed in the threat scenario. Step 3: if the path exists, the state X can be employed to point to the destination G by the back pointer. Besides, if the f∗ changes, the current path cost C∗ will be immediately recalculated. Step 4: insert the affected state which near the new threat into the OPEN list, and then return to Step 2. The pseudocodes of the D-Star algorithm are depicted as follows Algorithm 1.Algorithm 1
Input: the starting point E, destination GOutput: paths1X←Min_State2While~isemptyOPEN&&tX~=G3ifX=null,thenretuntoStep14tag←tagmin,deleteX5iftag<hX6ifhY≤tag&&hX>hY+CY,X7Calculateh∗andf∗8end9end10iftag=hX11iftY=NEWor12bY=X&hY≠hX+CX,Yor13bY≠X&hY>hX+CX,Ythen14bY←X,InseertY,hX+CX,Y15end16else17ifbY~=X&&hY>hX+CX,Y18InsertX,hX19end20else21InsertY,hY,CX,Y←cval22end23iftX=E&&tX~=CLOSED24InsertX,hX;Calculatef∗25else26returnstep227end28end


### 4.3. Improved A-Star Algorithm

We are focus on analyzing the learning real-time A-Star algorithm (LRTA-Star) which satisfies the requirements of real-time planning in a dynamic environment [32, 33]. In contrast, the path obtained by adopting the LRTA-Star algorithm is composed of a series of the tortuous track; it is difficult to achieve accurate track tracking control due to the limited maneuverability. Besides, in the flight process of stealth UAV, LRTA-Star algorithm is easy to fall into a local dead loop, which leads to the failure of path searching. Therefore, further combined with the idea of model-based predictive control (MPC) [34, 35], an improved A-Star algorithm with the multistep optimal search is proposed in this section.MPC is an optimization control method, which is mainly composed of model prediction, feedback correction, and rolling optimization. The principle of the MPC system is depicted in Figure2.Figure 2
The principle of the MPC system.The closed-loop output predictionypk+1 is given by
(13)ypk+1=ymk+1+h1yk−ymk,where uk is the actual control quantity on the system at moment k, yrk is the reference track softened by the input filter, ymk+1 is the predicted output value of the model, ypk+1 is the closed-loop output prediction, and h1 is the error correction coefficient.According to Eq. (13), the optimal control law for multistep prediction is given by
(14)J=∑i=1Nqiypk+i−yrk+i2+∑j=1Wδjuk+j−12,where N is the length of the prediction domain; W is the length of the control domain; qi is the output prediction error, and δj is the weighting coefficient of the control variable.In the process of a multistep optimization search, the track cost of each predicted track is considered as the cumulative value of the costs ofN predicted track nodes. Therefore, the objective function of UAV flight path planning can be obtained according to the optimal control law of the MPC system. The cost of the predicted flight path in N steps from the kth node is given by
(15)Jk=∑j=kk+Nb1Tj∣kB1b1j∣k+εjPNetj∣k+∑j=kk+WδjuTj∣kuj∣k,where PNetj∣k is the radar threat cost at the jth prediction point on the current predicted track segment, which can be obtained by Eq. (11). εj is the weight of threat cost; B1 is the distance cost weighted matrix, uj∣k is the control sequence, and b1j∣k is the distance cost between the coordinate of destination xe,ye,he and the coordinate of jth track point xj∣k,yj∣k,hj∣kT and b1j∣k=xj∣k−xe,yj∣k−ye.The stealth UAV is limited by its maneuverability in the process of a multistep search, whereu=u0,u1,u2,⋯uN−1 is sought by the current track point and heading angle, and the cost function Jk for the node has a minimum value. Therefore, add the parameter ω0≤ω≤1 to the heuristic function of the original algorithm, which aims to simplify the operation, reduce the search time and ensure the optimality of path planning, and the new heuristic function expression in the improved A-Star algorithm is given by
(16)fns=hns+1+ηiμikn,ns.Whereηi/μi is the ratio of the path cost of the current state to the path cost of an unknown region, n is the current node, ns is the adjacent extended node, and kn,ns is the cost from the current node to the adjacent node, and kn,ns=gns−gn. The pseudocode of the improved A-Star algorithm is depicted as follows Algorithm 2, and the flow chart of the improved A-Star algorithm is depicted in Figure 3.Algorithm 2
Input: the starting point E, destination GOutput: paths1n←x0,y0,ns←ns0n→OPEN2While~isemptyOPEN3fns←hns+1+ηi/μikn,ns4ifBBthreats=05ifPNet<Pc6ifns=G&&ns~=OPEN7n←ns,n→OPEN,outputpaths8else9n←nsi,reutrntoStep310else11ns→CLOSE,updatenandreturntoStep312else13Multi−StepSearch14hn←maxhn,minns∈Jkn,ns+hns15Updatethensthen16fns←hns+1+ηi/μikn,ns17CalculateσandPNet18n←ns,updatefnandreturntoStep514end15endFigure 3
The flow chart of the improved A-Star algorithm.


## 4.1. The A-Star Algorithm

The main idea of the standard A-Star algorithm is as follows. Firstly, select the appropriate heuristic function, estimate the generation value of the extensible search points in the search area, comprehensively. Moreover, compare the different cost values of each point. Moreover, consider the operation time and distance cost of the track point search. Finally, find an optimal path. In the A-Star algorithm, the operation of the OPEN list and CLOSE list are usually performed to achieve the storage and update of track points.However, the standard A-Star algorithm that is applied to find the optimal path for stealth UAV has many disadvantages: (a) the route obtained by the A-star algorithm only has the position of stealth UAV, and it can not reflect the characteristics of the dynamic RCS and attitude angle information. (b) Unknown path cost estimation is rarely calculated in the process of path planning accurately, and the computation time is too long to find a globally optimal path. (c) the performance of real-time is poor, and the algorithm can not deal with BB threats. Therefore, we are focus on improving the performance of the A-Star algorithm.


## 4.2. The D-Star Algorithm

The D-Star algorithm is developed by the A-Star algorithm and the Dijkstra algorithm [30, 31], which is suitable for addressing the path planning problem in unknown environments. The main idea of the algorithm is to search the reverse path from the destination to the start, and the heuristic function expression of the D-Star algorithm is given by
(12)fX,E=hX+gX,E,where hX represents the actual cost of the path from the destination to state X, and gX,E is the estimated cost of the path from the state X to the starting point.The main steps of the D-Star algorithm are given as follows. Step 1: the value of all states are set to NEW,f∗ and h∗ are set to infinity for all states, hG is set to zero, and the destination G is added to the OPEN list. Step 2: keep performing f∗ the search until the location of E which is removed from the OPEN list, if the CLOSE list contains the tag value of the state X, the complete path sequence will be obtained. In contrast, path planning has failed in the threat scenario. Step 3: if the path exists, the state X can be employed to point to the destination G by the back pointer. Besides, if the f∗ changes, the current path cost C∗ will be immediately recalculated. Step 4: insert the affected state which near the new threat into the OPEN list, and then return to Step 2. The pseudocodes of the D-Star algorithm are depicted as follows Algorithm 1.Algorithm 1
Input: the starting point E, destination GOutput: paths1X←Min_State2While~isemptyOPEN&&tX~=G3ifX=null,thenretuntoStep14tag←tagmin,deleteX5iftag<hX6ifhY≤tag&&hX>hY+CY,X7Calculateh∗andf∗8end9end10iftag=hX11iftY=NEWor12bY=X&hY≠hX+CX,Yor13bY≠X&hY>hX+CX,Ythen14bY←X,InseertY,hX+CX,Y15end16else17ifbY~=X&&hY>hX+CX,Y18InsertX,hX19end20else21InsertY,hY,CX,Y←cval22end23iftX=E&&tX~=CLOSED24InsertX,hX;Calculatef∗25else26returnstep227end28end


## 4.3. Improved A-Star Algorithm

We are focus on analyzing the learning real-time A-Star algorithm (LRTA-Star) which satisfies the requirements of real-time planning in a dynamic environment [32, 33]. In contrast, the path obtained by adopting the LRTA-Star algorithm is composed of a series of the tortuous track; it is difficult to achieve accurate track tracking control due to the limited maneuverability. Besides, in the flight process of stealth UAV, LRTA-Star algorithm is easy to fall into a local dead loop, which leads to the failure of path searching. Therefore, further combined with the idea of model-based predictive control (MPC) [34, 35], an improved A-Star algorithm with the multistep optimal search is proposed in this section.MPC is an optimization control method, which is mainly composed of model prediction, feedback correction, and rolling optimization. The principle of the MPC system is depicted in Figure2.Figure 2
The principle of the MPC system.The closed-loop output predictionypk+1 is given by
(13)ypk+1=ymk+1+h1yk−ymk,where uk is the actual control quantity on the system at moment k, yrk is the reference track softened by the input filter, ymk+1 is the predicted output value of the model, ypk+1 is the closed-loop output prediction, and h1 is the error correction coefficient.According to Eq. (13), the optimal control law for multistep prediction is given by
(14)J=∑i=1Nqiypk+i−yrk+i2+∑j=1Wδjuk+j−12,where N is the length of the prediction domain; W is the length of the control domain; qi is the output prediction error, and δj is the weighting coefficient of the control variable.In the process of a multistep optimization search, the track cost of each predicted track is considered as the cumulative value of the costs ofN predicted track nodes. Therefore, the objective function of UAV flight path planning can be obtained according to the optimal control law of the MPC system. The cost of the predicted flight path in N steps from the kth node is given by
(15)Jk=∑j=kk+Nb1Tj∣kB1b1j∣k+εjPNetj∣k+∑j=kk+WδjuTj∣kuj∣k,where PNetj∣k is the radar threat cost at the jth prediction point on the current predicted track segment, which can be obtained by Eq. (11). εj is the weight of threat cost; B1 is the distance cost weighted matrix, uj∣k is the control sequence, and b1j∣k is the distance cost between the coordinate of destination xe,ye,he and the coordinate of jth track point xj∣k,yj∣k,hj∣kT and b1j∣k=xj∣k−xe,yj∣k−ye.The stealth UAV is limited by its maneuverability in the process of a multistep search, whereu=u0,u1,u2,⋯uN−1 is sought by the current track point and heading angle, and the cost function Jk for the node has a minimum value. Therefore, add the parameter ω0≤ω≤1 to the heuristic function of the original algorithm, which aims to simplify the operation, reduce the search time and ensure the optimality of path planning, and the new heuristic function expression in the improved A-Star algorithm is given by
(16)fns=hns+1+ηiμikn,ns.Whereηi/μi is the ratio of the path cost of the current state to the path cost of an unknown region, n is the current node, ns is the adjacent extended node, and kn,ns is the cost from the current node to the adjacent node, and kn,ns=gns−gn. The pseudocode of the improved A-Star algorithm is depicted as follows Algorithm 2, and the flow chart of the improved A-Star algorithm is depicted in Figure 3.Algorithm 2
Input: the starting point E, destination GOutput: paths1n←x0,y0,ns←ns0n→OPEN2While~isemptyOPEN3fns←hns+1+ηi/μikn,ns4ifBBthreats=05ifPNet<Pc6ifns=G&&ns~=OPEN7n←ns,n→OPEN,outputpaths8else9n←nsi,reutrntoStep310else11ns→CLOSE,updatenandreturntoStep312else13Multi−StepSearch14hn←maxhn,minns∈Jkn,ns+hns15Updatethensthen16fns←hns+1+ηi/μikn,ns17CalculateσandPNet18n←ns,updatefnandreturntoStep514end15endFigure 3
The flow chart of the improved A-Star algorithm.


## 5. Numerical Results

Numerical simulations of penetration path planning are performed by employing the improved A-Star algorithm, LRTA-Star algorithm, and D-Star algorithm in different threat scenarios, which aims to verify the effectiveness of the improved A-Star algorithm. The simulation experiments are conducted with MATLAB2020Ra software and a Windows 10 system. The parameters of flight are depicted in Table1, the value K0 is determined by the number of radar M, and it will change when there is a BB threat in the threat scenario. Additionally, it is considered as the probability state of radar high detection when PNet exceeds 0.4, and Ph is performed to evaluate the safety degree of the track segment. The RCS data of a certain type of stealth UAV is depicted in Figure 4.Table 1
Parameters of path planning.
ParametersValuev500 (km/h)h1 (km)N4φmaxπ/3r6 (km)Pc0.5a,b,c(0.317,0.178,1.000)U9.8 (m/s2)θmaxπ/4Figure 4
The RCS data of a certain type of stealth UAV.


### 5.1. The First Threat Scenario

The threat region has a range of50km×50km; the geographical location and types of threat sources are presented in Table 2. The coordinate of the starting point is (5,5) km, and the coordinate of the target point is (90,90) km. The numerical results of penetration path planning which is performed by different algorithms in the first scenario are depicted in Figure 5. The variation of heading angle and roll angle which is applied by the improved A-Star algorithm in the first scenario is depicted in Figure 6, the detection probability of netted radar in the first threat scenario is depicted in Figure 7, and the statistical result of flight in the first scenario is presented in Table 3.Table 2
The geographical location and types of threat sources in the first threat scenario.
ThreatGeographic locations (km)TypeThreat module 1(15,10)Known threatThreat module 2(25,25)Known threatThreat module 3(40,35)Known threatThreat module 4(30,35)BB threatFigure 5
Numerical results of penetration path planning in the first scenario.
(a)
Improved A-Star algorithm(b)
LRTA-Star algorithm(c)
D-Star algorithmFigure 6
The variation of attitude angle that performed by improved A-Star algorithm in the first scenario.
(a)
The variation of heading angle(b)
The variation of roll angleFigure 7
The detection probability of netted radar in the first threat scenario.
(a)
Original planning(b)
ReplanningTable 3
The statistical result of flight in the first scenario.
SituationAlgorithmDistance (km)PhRun timeOriginal planningImproved A-Star62.3720.4085.372LRTA-Star67.7430.4188.276D-Star61.8120.4318.749ReplanningImproved A-Star67.4150.4177.508LRTA-Star70.8330.42914.624D-Star65.1520.44213.391Figure5 demonstrates the penetration paths for stealth UAV which is performed by employing three different algorithms in a scenario. From Figure 7 and Table 3, we can infer that the stealth UAV can achieve penetration path planning by employing in the first scenario in which a static or dynamic environment. However, compared with the other two algorithms, the paths obtained by employing the improved A-Star algorithm has a smaller value of Ph and takes less time. Figure 6 demonstrates that the improved A-Star algorithm can accurately adjust the attitude angle and reduce the value of RCS on the paths, and this algorithm achieves the rapid and safe path penetration planning for stealth UAV. Besides, the paths of the improved A-Star algorithm is shorter than the LRTA-Star algorithm, and the improved A-Star algorithm has higher path planning efficiency and safety, which further proves the effectiveness of the improved A-Star algorithm.


### 5.2. The Second Threat Scenario

The threat region has a range of50km×50km; the geographical location and types of threat sources are presented in Table 4. The coordinate of the starting point is (5,5) km, and the coordinate of the target point is (90,90) km. The numerical results of penetration path planning which is performed by different algorithms in the second scenario are depicted in Figure 8. The variation of heading angle and roll angle which is performed by an improved A-Star algorithm in the second scenario is depicted in Figure 9, the detection probability of netted radar in the second threat scenario is depicted in Figure 10, and the statistical result of flight in the second scenario is presented in Table 5.Table 4
The geographical location and types of threat sources in the second threat scenario.
ThreatGeographic locations (km)TypeThreat module 1(15,10)Known threatThreat module 2(15,20)Known threatThreat module 3(40,25)Known threatThreat module 4(30,35)Known threatThreat module 5T (38,15), R (23,15)Known threatThreat module 6(20,30)BB threatThreat module 7(37,40)BB threatFigure 8
Numerical results of penetration path planning in the second scenario.
(a)
Improved A-Star algorithm(b)
LRTA-Star algorithm(c)
D-Star algorithmFigure 9
The variation of attitude angle employed by improved A-Star algorithm in the second scenario.
(a)
The variation of heading angle(b)
The variation of roll angleFigure 10
The detection probability of netted radar in the second threat scenario.
(a)
Original planning(b)
ReplanningTable 5
The statistical result of flight in the second scenario.
SituationAlgorithmDistance (km)PhRun timeOriginal planningImproved A-Star65.1920.4296.053LRTA-Star68.8790.45910.504D-Star65.1520.47011.581ReplanningImproved A-Star70.4250.4469.827LRTA-Star70.1910.50216.117D-Star69.3210.53518.659Figure8 describes the penetration paths for stealth UAV which is performed by applying three different algorithms in a scenario with more complex threat sources, and from Figure 10 and Table 5. For the original path planning, the penetration paths which are safe by employing all three algorithms in the scenario. However, for the path replanning in the presence of BB threats, the path replanning has failed in the second threat scenario, because the Ph of paths by adopting the LRTA-Star algorithm and D-star algorithm are higher than 0.5. Inversely, stealth UAV can still achieve the path replanning in the second scenario with two BB threats by using an improved A-Star algorithm. Moreover, the improved A-Star algorithm takes less time and has higher path planning efficiency than the other two algorithms whether in the original path planning or replanning, and it further proves the validity of the algorithm. Figure 9 indicates that when the improved A-Star algorithm is performed to penetration path planning, the variation of heading angle and roll angle are within the preset range, which can satisfy the real flight requirements.


### 5.3. The Third Threat Scenario

The threat region has a range of50km×50km; the geographical location and types of threat sources are presented in Table 6. The coordinate of the starting point is (5,5) km, and the coordinate of the target point is (90,90) km. The numerical results of penetration path planning which is performed by different algorithms in the third scenario are depicted in Figure 11. The variation of heading angle and roll angle which is performed by an improved A-Star algorithm in the third scenario is depicted in Figure 12, the detection probability of netted radar in the third threat scenario is depicted in Figure 13, and the statistical result of flight in the third scenario is presented in Table 7.Table 6
The geographical location and types of threat sources in the third threat scenario.
ThreatGeographic locations (km)TypeThreat module 1(7,22)Known threatThreat module 2(15,30)Known threatThreat module 3(40,20)Known threatThreat module 4(26,35)Known threatThreat module 5(37,40)Known threatThreat module 6(44,30)Known threatThreat module 7(12,40)Known threatThreat module 8(42,12)Known threatThreat module 9T (33,13), R (15,13)Known threatThreat module 10(30,25)BB threatThreat module 11(20,42)BB threatThreat module 12(33,43)BB threatFigure 11
Numerical results of penetration path planning in the third scenario.
(a)
Improved A-Star algorithm(b)
LRTA-Star algorithm(c)
D-Star algorithmFigure 12
The variation of attitude angle that performed by improved A-Star algorithm in the third scenario.
(a)
The variation of heading angle(b)
The variation of roll angleFigure 13
The detection probability of netted radar in the third threat scenario.
(a)
Original planning(b)
ReplanningTable 7
The statistical result of flight in the third scenario.
SituationAlgorithmDistance (km)PhRun timeOriginal planningImproved A-Star67.3160.4579.259LRTA-Star71.9090.52116.746D-Star82.3390.63919.468ReplanningImproved A-Star69.7880.47617.362LRTA-Star73.5200.57028.732D-Star84.8250.69834.085Figure11 describes the penetration paths for stealth UAV which are performed by applying three different algorithms in a scenario with a large number of threat sources and BB threats; from Figure 13 and Table 7, we can see that the Ph paths by adopting the LRTA-Star algorithm and D-star algorithm are higher than 0.5 no matter in the original path planning or path replanning. Obviously, stealth UAV can rarely achieve the penetration path planning which is safely and rapidly by employing these two algorithms. However, the Ph of paths by adopting improved A-Star algorithms are always lower than 0.5; therefore, stealth UAV can still achieve penetration path planning and replanning by using the improved A-Star algorithm in dynamic threat scenarios. Furthermore, paths of the improved A-Star algorithm have higher efficiency and safety, which further proves the validity of the improved A-Star algorithm in the scenario with high threat density. Figure 12 indicates that stealth UAV can meet constraints of attitude angle in the field of close and low altitude combat.


## 5.1. The First Threat Scenario

The threat region has a range of50km×50km; the geographical location and types of threat sources are presented in Table 2. The coordinate of the starting point is (5,5) km, and the coordinate of the target point is (90,90) km. The numerical results of penetration path planning which is performed by different algorithms in the first scenario are depicted in Figure 5. The variation of heading angle and roll angle which is applied by the improved A-Star algorithm in the first scenario is depicted in Figure 6, the detection probability of netted radar in the first threat scenario is depicted in Figure 7, and the statistical result of flight in the first scenario is presented in Table 3.Table 2
The geographical location and types of threat sources in the first threat scenario.
ThreatGeographic locations (km)TypeThreat module 1(15,10)Known threatThreat module 2(25,25)Known threatThreat module 3(40,35)Known threatThreat module 4(30,35)BB threatFigure 5
Numerical results of penetration path planning in the first scenario.
(a)
Improved A-Star algorithm(b)
LRTA-Star algorithm(c)
D-Star algorithmFigure 6
The variation of attitude angle that performed by improved A-Star algorithm in the first scenario.
(a)
The variation of heading angle(b)
The variation of roll angleFigure 7
The detection probability of netted radar in the first threat scenario.
(a)
Original planning(b)
ReplanningTable 3
The statistical result of flight in the first scenario.
SituationAlgorithmDistance (km)PhRun timeOriginal planningImproved A-Star62.3720.4085.372LRTA-Star67.7430.4188.276D-Star61.8120.4318.749ReplanningImproved A-Star67.4150.4177.508LRTA-Star70.8330.42914.624D-Star65.1520.44213.391Figure5 demonstrates the penetration paths for stealth UAV which is performed by employing three different algorithms in a scenario. From Figure 7 and Table 3, we can infer that the stealth UAV can achieve penetration path planning by employing in the first scenario in which a static or dynamic environment. However, compared with the other two algorithms, the paths obtained by employing the improved A-Star algorithm has a smaller value of Ph and takes less time. Figure 6 demonstrates that the improved A-Star algorithm can accurately adjust the attitude angle and reduce the value of RCS on the paths, and this algorithm achieves the rapid and safe path penetration planning for stealth UAV. Besides, the paths of the improved A-Star algorithm is shorter than the LRTA-Star algorithm, and the improved A-Star algorithm has higher path planning efficiency and safety, which further proves the effectiveness of the improved A-Star algorithm.


## 5.2. The Second Threat Scenario

The threat region has a range of50km×50km; the geographical location and types of threat sources are presented in Table 4. The coordinate of the starting point is (5,5) km, and the coordinate of the target point is (90,90) km. The numerical results of penetration path planning which is performed by different algorithms in the second scenario are depicted in Figure 8. The variation of heading angle and roll angle which is performed by an improved A-Star algorithm in the second scenario is depicted in Figure 9, the detection probability of netted radar in the second threat scenario is depicted in Figure 10, and the statistical result of flight in the second scenario is presented in Table 5.Table 4
The geographical location and types of threat sources in the second threat scenario.
ThreatGeographic locations (km)TypeThreat module 1(15,10)Known threatThreat module 2(15,20)Known threatThreat module 3(40,25)Known threatThreat module 4(30,35)Known threatThreat module 5T (38,15), R (23,15)Known threatThreat module 6(20,30)BB threatThreat module 7(37,40)BB threatFigure 8
Numerical results of penetration path planning in the second scenario.
(a)
Improved A-Star algorithm(b)
LRTA-Star algorithm(c)
D-Star algorithmFigure 9
The variation of attitude angle employed by improved A-Star algorithm in the second scenario.
(a)
The variation of heading angle(b)
The variation of roll angleFigure 10
The detection probability of netted radar in the second threat scenario.
(a)
Original planning(b)
ReplanningTable 5
The statistical result of flight in the second scenario.
SituationAlgorithmDistance (km)PhRun timeOriginal planningImproved A-Star65.1920.4296.053LRTA-Star68.8790.45910.504D-Star65.1520.47011.581ReplanningImproved A-Star70.4250.4469.827LRTA-Star70.1910.50216.117D-Star69.3210.53518.659Figure8 describes the penetration paths for stealth UAV which is performed by applying three different algorithms in a scenario with more complex threat sources, and from Figure 10 and Table 5. For the original path planning, the penetration paths which are safe by employing all three algorithms in the scenario. However, for the path replanning in the presence of BB threats, the path replanning has failed in the second threat scenario, because the Ph of paths by adopting the LRTA-Star algorithm and D-star algorithm are higher than 0.5. Inversely, stealth UAV can still achieve the path replanning in the second scenario with two BB threats by using an improved A-Star algorithm. Moreover, the improved A-Star algorithm takes less time and has higher path planning efficiency than the other two algorithms whether in the original path planning or replanning, and it further proves the validity of the algorithm. Figure 9 indicates that when the improved A-Star algorithm is performed to penetration path planning, the variation of heading angle and roll angle are within the preset range, which can satisfy the real flight requirements.


## 5.3. The Third Threat Scenario

The threat region has a range of50km×50km; the geographical location and types of threat sources are presented in Table 6. The coordinate of the starting point is (5,5) km, and the coordinate of the target point is (90,90) km. The numerical results of penetration path planning which is performed by different algorithms in the third scenario are depicted in Figure 11. The variation of heading angle and roll angle which is performed by an improved A-Star algorithm in the third scenario is depicted in Figure 12, the detection probability of netted radar in the third threat scenario is depicted in Figure 13, and the statistical result of flight in the third scenario is presented in Table 7.Table 6
The geographical location and types of threat sources in the third threat scenario.
ThreatGeographic locations (km)TypeThreat module 1(7,22)Known threatThreat module 2(15,30)Known threatThreat module 3(40,20)Known threatThreat module 4(26,35)Known threatThreat module 5(37,40)Known threatThreat module 6(44,30)Known threatThreat module 7(12,40)Known threatThreat module 8(42,12)Known threatThreat module 9T (33,13), R (15,13)Known threatThreat module 10(30,25)BB threatThreat module 11(20,42)BB threatThreat module 12(33,43)BB threatFigure 11
Numerical results of penetration path planning in the third scenario.
(a)
Improved A-Star algorithm(b)
LRTA-Star algorithm(c)
D-Star algorithmFigure 12
The variation of attitude angle that performed by improved A-Star algorithm in the third scenario.
(a)
The variation of heading angle(b)
The variation of roll angleFigure 13
The detection probability of netted radar in the third threat scenario.
(a)
Original planning(b)
ReplanningTable 7
The statistical result of flight in the third scenario.
SituationAlgorithmDistance (km)PhRun timeOriginal planningImproved A-Star67.3160.4579.259LRTA-Star71.9090.52116.746D-Star82.3390.63919.468ReplanningImproved A-Star69.7880.47617.362LRTA-Star73.5200.57028.732D-Star84.8250.69834.085Figure11 describes the penetration paths for stealth UAV which are performed by applying three different algorithms in a scenario with a large number of threat sources and BB threats; from Figure 13 and Table 7, we can see that the Ph paths by adopting the LRTA-Star algorithm and D-star algorithm are higher than 0.5 no matter in the original path planning or path replanning. Obviously, stealth UAV can rarely achieve the penetration path planning which is safely and rapidly by employing these two algorithms. However, the Ph of paths by adopting improved A-Star algorithms are always lower than 0.5; therefore, stealth UAV can still achieve penetration path planning and replanning by using the improved A-Star algorithm in dynamic threat scenarios. Furthermore, paths of the improved A-Star algorithm have higher efficiency and safety, which further proves the validity of the improved A-Star algorithm in the scenario with high threat density. Figure 12 indicates that stealth UAV can meet constraints of attitude angle in the field of close and low altitude combat.


## 6. Conclusions

This paper presented a new solution for path replanning for stealth unmanned aerial vehicle in a complex radar net environment with BB threats.We are focus on analyzing the kinematics model of stealth UAV, threat source in penetration environment, and path planning algorithm. Further combined with model-based predictive control (MPC) and the LRTA-Star algorithm, an improved A-Star algorithm is proposed to achieve penetration path planning and replanning for stealth UAV. Meanwhile, the attitude angle information of stealth UAV is added to the algorithm, which demonstrates the variation characteristics of dynamic RCS. Further combined with the kinematics analysis of stealth UAV and the detection performance analysis of radar net, the original paths, and the replanning paths can satisfy the real flight requirements.Compared with the other two algorithms, the stealth UAV can rapidly achieve the penetration path planning by employing the improved A-Star algorithm in a different complex environment with BB threats, improve the survivability of stealth UAV, and the efficiency of penetration path planning.The threat scenarios are composed of single-base radar, dual-base radar, and BB threats, which are closer to the real combat environment. The model and the improved A-Star algorithm proposed in this paper can quickly generate better penetration paths in the combat area under a dynamic environment, exhibiting certain military application value.In future work, we are focus on the real-time rapid penetration path planning for stealth UAV account for a three-dimensional complex dynamic environment with the terrain and unknown motion targets in the field of close and low altitude combat.



---

*Source: 8896357-2020-08-01.xml*