Integrated Guidance and Control Based Air-to-Air Autonomous Attack Occupation of UCAV

An approach of air-to-air autonomous attack occupation for Unmanned Combat Aerial Vehicles (UCAVs) is proposed to improve attack precision and combat effectiveness. According to the shortage of UCAV in the task of attack occupation, kinematic and dynamicmodels ofUCAVandmissile loaded on it are formed.Then, attack zone andno-escape zone are calculated by pattern search algorithm, and the optimum attack position is indicated. To arrive at the optimum attack position accurately with restriction of gesture, a novel adaptive sliding mode control method is suggested to design the integrated guidance and control system of UCAV in the process of autonomous attack occupation. Key parameters of the control system are adaptively regulated, which further economize control energy at the same time.The simulation results show that compared with traditional methods our approach can guide the UCAV to the optimum attack position with stable gesture and economize nearly 25% control energy.


Introduction
In most UCAV systems (e.g., Predator.MQ-1B), operator on the ground is a key component of the whole control loop.However, the connatural limitation of operator brings drawback of the UCAV system by limiting the precision of estimation and time of response.In fact, the UCAV and operator are in different frame of time-space.To solve this problem, technology of autonomous attack for UCAV has become a hot topic in recent years [1][2][3][4][5].Yet, many achievements of researchers are aimed at attacking ground fixed targets [2][3][4].Technologies of air-to-air autonomous attack occupation for UCAV still remain as a novel field [5].
Air-to-air autonomous attack occupation of UCAV is mainly comprised of two relative processes.Firstly, after having detected target in the air, calculating the optimal position that UCAV should maneuver to is crucial for the follow-up processes.As to aerial maneuvering target, Zhang et al. [6] calculate the real-time attack zone with golden section method.However, the maneuver capability of UCAV is relatively poor for its high-aspect-ratio, and its attack zone is relatively small.Compared with the timeliness, the precision of attack occupation is more important [7].In [8], a dynamic attack zone is calculated, which can be reliably and accurately used during any flight antagonism time.But, for autonomous attack occupation of UCAV, more factors should be taken into consideration, such as the no-escape zone and optimum attack position [1,9].Secondly, after having attained the optimum attack position, UCAV should maneuver to it accurately and rapidly.Integrated guidance and control (IGC) can eliminate the time-lag that inevitably exists in the conventional guidance and control design that separates the guidance loop from the control loop [10][11][12].In addition, IGC improves the precision of the terminal phase of guidance [13][14][15].Inspired by this, we propose combining IGC with the maneuver process of attack occupation of UCAV.Many researches have confirmed the effectiveness of IGC in various guidance and control scenarios [16][17][18][19][20].In [19], high order sliding mode control (SMC) approach is incorporated with IGC to design the interception of maneuvering hypersonic vehicle.Guo and Liang [20] design an IGC system based on the block backstepping sliding mode and extended state observer for the near-space interceptor.However, most researches above are based on the scenario 2 Mathematical Problems in Engineering in which the maneuver of target could be estimated exactly.For the scenario of autonomous attack occupation of UCAV, more attention should be placed on automaticity and adaptability.
Therefore, motivated by the aforementioned considerations we compute the attack zone and no-escape zone of UCAV with pattern search algorithm, which is famous for its precision.Then, by combining the practicality of air combat to the no-escape zone computed before, we attain the optimum attack position.Furthermore, a novel adaptive SMC method is used to design the IGC system of UCAV during the process of maneuvering to the optimum attack position.The proposed control system satisfies the gesture constraints of UCAV at the attack moment.In addition, the key parameters of the control system are timely and adaptively regulated in the attack occupying process, which brings better performance and less consumption of control energy in the simulation scenario.
Our main contribution is that we design the air-toair autonomous attack occupation of UCAV based on the key factors of air combat.After analyzing the process of autonomous attack occupation theoretically and practically, we compute the optimum attack position.Then, we combine a novel adaptive SMC method to the IGC system of UCAV, so that the UCAV can maneuver to the optimum attack position rapidly and accurately and satisfy the gesture constraints at the attack moment.To our knowledge, there have not been any attempts of applying IGC approach to UCAV in the airto-air autonomous attack occupation.
The rest of the paper is organized as follows.In Section 2, the task and kinematic model of air-to-air autonomous attack are derived.Section 3 introduces the computation of attack zone, no-escape zone, and optimum attack position, which provides the target position that the UCAV should maneuver to with our proposed integrated guidance and control approach explained in Section 4. Numerical simulations are demonstrated in Section 5 and we conclude the paper in Section 6.

Air-to-Air Autonomous Attack
Occupation of UCAV  1.In the phase of takeoff and cruise, UCAV fulfills path planning and real-time preplanning according to the battlefield environment and mission.When airborne detection system finds aerial target, the tracking system provides target information to the integrated fire and flight (IFF) computer.The computer estimates target state and computes the attack zone.Then according to the condition of the guided missile, the range of attack zone, and the threat from target, UCAV is guided to the optimum attack position.When the weapon management system meets the need of launch conditions, UCAV launches the guided missile.

Kinematic Modeling of Air-to-Air Attack Task.
UCAV fulfills air-to-air autonomous attack in three-dimensional space.The relative movement of UCAV and the target can be separately studied in the vertical and lateral plane, as is shown in Figure 2.
and   denote the line-of-sight (LOS) angle in vertical and lateral planes, respectively.
In the vertical plane, the relative movement of UCAV and the target is shown in Figure 3.   is the relative distance between them.  ,   , and   are the velocity, normal acceleration, and flight-path angle of UCAV, respectively.The velocity, normal acceleration, and path angle of target are denoted by   ,   , and   respectively.
The relative movement of UCAV and the target shown in Figure 3 can be described by Similar to (1), the relative movement in lateral plane is given as follows: In this paper, the guided missile is axisymmetric and irrotational.Because the calculation of attack zone is a process of generating trajectory substantially, the guided missile is modeled in three degrees of freedom: where (  ,   , ℎ  ) denotes the missile's position.V  ,   ,   denote missile's velocity, path angle, and azimuth angle, respectively.  is the thrust of the missile, and   is the drag force of it.  is the mass of the missile.  and   are the -load of the missile in pitching and yawing direction, respectively.

The Optimum Attack Position
3.1.Calculation of Attack Zone.During the implementation of autonomous attack task, UCAV faces up to corresponding threat from the target.Concerning the properties of UCAV, such as abundant useable -load, relatively small payload, and small attack zone, the precision of attack zone becomes essential.We compute the attack zone with the pattern search algorithm for higher precision [21][22][23][24].Pattern search algorithm includes two types of moves: (i) exploratory move, which is made on current point by investigation along each independent search direction, and (ii) pattern move, which finds the new point by jumping from present base point along with pattern direction [25].The two types of moves are alternated to find the boundary of attack distance.Pattern search algorithm provides global optimal results without derivatives of objective functions [26], and it is further illustrated in Figure 4.
In Figure 4, search directions along the axes at points  () and  (+1) indicate exploratory move, which are drawn as solid lines with arrowheads.On the other hand, pattern direction along the broken line from  () to  (+1) indicates pattern move.Furthermore, we introduce a contraction parameter
Algorithm 1: Pattern search algorithm with variable step size. ∈ (0, 1) to the traditional pattern search algorithm.The contraction parameter can change the size of search step to speed up search process at the initial phase as well as improve the search accuracy at the terminal phase.Suppose () is the distance between UCAV and the target in accordance with all the constraints.Set the iterative number  = 0, and we do as shown in Algorithm 1.
The pattern search algorithm above is used to calculate the near boundary of attack distance.The far boundary of attack distance can be obtained by replacing the object function () with −().
The calculation flow of attack distance can be further shown in Figure 5.
The process of calculating attack zone involves the inner loop and outer loop.The inner loop is about to search the maximum or minimum distance in a certain entry angle with respect to the target.The outer loop searches through the entry angles, whose bound is [0, 2].Then, attack zone can be calculated by means of interpolating or fitting method.
The outer loop of attack zone calculation is presented in Figure 6.
For the missile loaded on UCAV, assume that its offaxis angle is 0, initial weight is 22 kg, initial velocity V 0 = 300 m/s, consumption of fuel per second   = 1.5 kg/s, working time of engine  0 = 4 s, average thrust  = 1000 N, delay time of actuator  = 0.02 s, maximum LOS angle  max = 60 ∘ , maximum tracking rate of LOS angle  max = 20 ∘ /s, available -load   = 20 g, unlock time of fuze   = 0.5 s, and maximum available distance of fuze is 50 m.At the same time, assume that the target flies with no lateral -load and its velocity   = 300 m/s.Considering the combat height  = 6 km, the attack zone calculated by pattern search algorithm is illustrated in Figure 7.
As shown in Figure 7, the far boundary of attack zone becomes larger with the increase of entry angle , and it reaches the maximum when entry angle  = 180 ∘ .Then, it gradually shrinks.This is consistent with our intuition that the far boundary of attack zone is mainly affected by the approaching speed of the missile and target.The missile can reach the target in a relatively far distance when it is launched in front of the target.On the contrary, the attack distance decreases when the missile is launched behind the target.On the other hand, the near boundary of attack zone is mainly affected by the approaching speed of them within the work time of fuze on the missile.The variation trend of near boundary is the same as that of far boundary.

Calculation of No-Escape Zone.
The practice of air combat shows that launching missile within the attack zone cannot guarantee destroying the target when it implements large overload maneuver.Thus, after calculating the attack zone, we further investigate the attack zone in two other cases, in which the aerial target maneuvers in the left direction and right direction with maximum -load.Then the noescape zone is attained by choosing the minimum of the three far boundaries of attack zone as the far boundary and the maximum of the three near boundaries of attack zone as the near boundary:

𝑅
no-escape min Assuming the maximum -load of target | max | = 3, the calculated no-escape zone is shown in Figure 8.
Compared with the attack zone in Figure 7, the far boundary of no-escape zone in Figure 8 shrinks.In the case of launching missile in front of the target, the maximum attack distance decreases about 1 km.On the other hand, the near boundary of no-escape zone becomes a little larger than that of attack zone in Figure 7.

Optimum Attack Position for UCAV.
Combining the practicality of air combat to the implementation of air-toair autonomous attack, UCAV should satisfy the following qualifications: (1) Keep away from the detection range of fire-control radar that is loaded on the target.
(2) Considering the target's maneuver ability, the probability of destroying it by guided missile should be as high as possible.
(3) Distance between UCAV and the target should be as far as possible, and, before being discovered by the target, attack first.
Assume that the range of the azimuth angle of target's fire-control radar is ±45 ∘ .Considering the calculating result of attack zone and no-escape zone, as well as the specified qualifications above, the optimum attack position of UCAV can be illustrated in Figure 9.
Figure 9 shows two optimum attack positions of the UCAV relative to the left and right direction of the target.Both attack positions are on the far boundary of the no-escape zone.When UCAV launches missile at these two positions, no matter how the target maneuvers at the combat height the launched missile can destroy it.In addition, these two points are both on the boundary of the range that the target's firecontrol radar can reach.The probability of being detected by target is relatively low in this case.Finally, the choice of the specific position depends on the UCAV's azimuth relative to the target when UCAV detects it.

Autonomous Attack Occupation Based on Integrated Guidance and Control
Before autonomous attack occupation, UCAV usually makes full use of its stealth characteristics and gets close to the enemy aircraft with negative height.Then it maneuvers to the optimum attack position and launches the guided missile.Figure 10 describes the whole process.To achieve the process of autonomous attack occupation rapidly and accurately, integrated guidance and control system is designed, respectively, for pitch channel, yaw channel, and roll channel of the UCAV.

Gesture Constraints at the Attack Moment.
When launching the guided missile, the constraints of UCAV's gesture are substantially about constraints of pitch angle  and yaw angle .In the terminal phase of guidance, the angle of attack  ≈ 0, the angle of sideslip  ≈ 0,  =   +, and  =   +.Define   and   as the pitch angle and yaw angle of the UCAV at the attack moment.The gesture constraints can be described as According to conclusion in Section 3.3, which is illustrated in Figure 9, we can obtain where   is the azimuth angle of UCAV with respect to the enemy aircraft at the attack moment.
During the process of autonomous attack occupation, the UCAV should meet the requirements that the -load of it may be as small as possible and the maneuver path of it may be as straight as possible.That is to say, we hope that the UCAV can be guided with its LOS angular rate approaching 0 as soon as possible, and it can be described as follows: Then, the terminal LOS angle can be described as Equations ( 8) and ( 9), respectively, describe the gesture constraints on condition that the UCAV implements autonomous attack occupation from the left and right directions with relativity to the target.

Modeling and Standardizing the IGC System of UCAV 4.2.1. Pitch Channel.
According to the relative movement between UCAV and the target described by (1), And in the terminal phase of attack occupation So λ  can be further described as where    = cos(  −   )  ≈   .What is more, because the target's acceleration is small before it has detected the UCAV, we could consider its acceleration as an unknown bounded variable.Let   1 = cos(  −   )  , and ( 12) becomes Because of the high-aspect-ratio structure, lateral force caused by sideslip angle has little effect on the UCAV, and hence where , , , and  are wing area, engine thrust, roll angle, and mass of the UCAV, respectively. denotes the atmospheric density and    is an aerodynamic constant and    denotes the approximate error of lateral acceleration in the vertical plane.In this case, (13) can be substituted by Mathematical Problems in Engineering In (15),    denotes an unknown bounded variable, which contains the approximate error of lateral acceleration and the unknown information of target's maneuver.
According to the gesture constraints analyzed in Section 4.1, UCAV's angle of attack is small when it arrives at the optimum attack position.
So, at this moment Then, the state variables, control variable, and output variables can be described as follows: Combining ( 17) with (15) and considering the coupling between the pitch channel and other channels as unknown bounded variable, we can obtain the IGC model of pitch channel as follows: where To design the control law with adaptive SMC method, IGC model should be simplified into a standard type, and it can be presented as follows.
Consider x = [ 1 ,  2 ,  3 ,  4 ]  and   =   , and we can obtain where In the terminal phase of attack occupation,   ≈   .Considering the target's acceleration as an unknown bounded variable, we obtain where    = cos(  −  )  ≈   .Furthermore, according to the dynamic analysis of the UCAV, its roll angle is always restricted in small scope in the process of autonomous attack occupation.Hence, let sin  ≈  and    can be attained as In ( 23)    is the approximate error of lateral acceleration in the lateral plane.
Substituting (23) into (22) yields where    denotes an unknown bounded variable, which contains the approximate error of lateral acceleration and the unknown information of target's maneuver.
According to the gesture constraints of UCAV in lateral plane, which have been analyzed in Section 4.1, we can obtain Define the state variables, control variable, and output variables as follows: Combining ( 26) with (24) and considering the coupling between the roll channel and other channels as unknown bounded variable, we can obtain the IGC model of roll channel as follows: where The IGC model of roll channel has the same form as that of the pitch channel, and its standard form can be got by using the same method in the analysis of pitch channel.

Yaw Channel.
As to yaw channel, define the state variables, control variable, and output variables as follows Seeing the coupling between the yaw channel and other channels as unknown bounded variable, we can obtain where  11 = ( 2    − 2)/2,  21 =  2    /2  , and  2 =  2     /2  .In order to standardize the control system, we further set the state variables x = [ 1 ,  2 ]  and control variable   =   .Then, the state representation can be described as follows: where (31)

A SMC Law for the IGC System of UCAV.
By taking external disturbances and internal parameter perturbation into consideration, a sliding mode control law is designed to make the IGC system of UCAV fleetly reach the sliding surface with a stable sliding mode.As to the process of autonomous attack occupation, the purpose of the SMC law is to make the LOS angular rate between UCAV and the optimum attack position approach 0 in the condition that the UCAV meets the gesture constraints at the attack moment.SMC is an effective robust control algorithm since it is insensitive to model uncertainties, external disturbances, and parameter variations [27].In SMC, the state of the control system is first guided to a designed surface (i.e., the sliding mode surface) and then it is forced there with a shifting law [28].The control strategy of SMC is determined by the sum of two control laws: switching and equivalent control laws.The switching control law is discontinuous, and it enforces the system state to the predetermined sliding mode surface in a finite time.When the sliding mode occurs, the closedloop system becomes insensitive to parameter variations and robust to matched disturbances by means of the equivalent control law which is usually linear [29].

Design of the SMC Law for IGC System in Pitch and
Roll Channels.As to the IGC system denoted by (19), set   =  1 =   −   and   =   and the sliding surface can be designed as where   is a positive constant called sliding mode coefficient and / is the differential operator with respect to   .Then, the control law can be designed as follows: where In (34),   is a designed positive constant called controller coefficient and d ≥ |3 2   2 + 3   3 +  4 | is the bound of unknown uncertainties.
Furthermore, the existence, accessibility, and stability of the sliding mode are analyzed as follows.
Consider the following Lyapunov function candidate: By differentiating (32) with respect to time, it implies (37) Therefore, the system is stable and it can reach the sliding surface in finite time.
The design of control law for roll channel is similar to that of pitch channel in format.However, its specific parameters are decided by the IGC model of roll channel.

Design of the SMC Law for IGC Model in Yaw Channel.
As to the IGC system denoted by (30), set   =  1 =  and   =   and the sliding surface can be designed as follows: where   is a positive constant called sliding mode coefficient and / is the differential operator with respect to   .Then, the control law yields where In ( 40),   is a designed positive constant which is called controller coefficient and d ≥ |    + ḋ  | is the bound of unknown uncertainties.
Consider the following Lyapunov function candidate: By differentiating (38) with respect to time, it implies Substitute ( 42) into (41), and we obtain (43) Therefore, the system is stable and it can reach the sliding surface in finite time.
In order to further smooth the control process, we replace the switching function sgn() with a saturation function sat(/), in which  > 0 denotes the width of boundary layer.Then the designed control laws of pitch, yaw, and roll channels are shown as follows: where   ,   , and   are the angle of elevator, rudder, and aileron of the UCAV, respectively.

Regulation of Parameters in the Designed Control System of UCAV.
Taking a close look into the designed control system, we find that the sliding mode, response time, and control error of it are mainly affected by the system parameters  1 ,  2 , , and .(1) Larger values of  1 and  2 shorten the time that control system spends to reach the sliding mode surface.However, they bring obvious chattering phenomenon and more consumption of control energy.So we prefer larger values of  1 and  2 when the system state is far from sliding mode surface and smaller values of  1 and  2 are more suitable on the contrary.(2) After reaching the sliding mode state, the control system tends to be stable in negative exponential form with parameter , and we prefer a larger value of  to save the convergence time.But a large value of  may trigger the unmodeled high frequency state when the control system is not in the sliding mode state, and a smaller  is preferred in this case.(3) The width of boundary layer  can indirectly handle the steady-state error of our control system.A smaller  brings smaller steady-state error but higher frequency of switching control and more consumption of control energy.Therefore, we consider a larger  when the system state is far from sliding mode surface to speed up the reaching process.On the contrary, a smaller  is preferred to reduce the steady-state error when our control system is reaching the sliding mode state.Moreover, we find that all of the parameters have a close relationship with the sliding surface .Therefore, we design a law in which main parameters of the designed control system are regulated adaptively with the change of sliding surface  to meet better control performance [30].Figure 11 shows the regulating law of the parameters in the designed control system.
In this paper, a negative exponential function is used, within which || is an independent variable.The variable || can adaptively regulate the main parameters in the designed control system.According to the regulating law shown in Figure 11, the detailed description of  1 ,  2 , , and  is as follows: Considering the regulating law described above, the control input in the pitch channel is   = (1/ 4 )(û+ũ), where the linear control input part û and the discontinuous control input part ũ are as follows: where   ,  1 , and  2 are in the format described in (45).In this case, we further set  2min ≥ d +   , where   is the controller coefficient and d ≥ |3 In addition, from the regulating law, ε  can be presented as Substitute ( 46) and ( 48) into (47), and we obtain where Thus, the control law given by ( 46) and adaptive regulating law in (45) guarantee the reaching and sustaining of the sliding mode.
We further replace sgn(  ) in (46) by the saturation function sat(  /  ), where the width of boundary layer   obeys the adaptive regulating law   =  max + ( min −  max ) − 2  /2 .By regulating   adaptively, the control system brings smaller steady-state error.At the same time, there is less chattering in the control process.
In the yaw channel and roll channel, we can draw the same conclusion as we discussed above.

Simulation Scenario and Initial Conditions.
Assume that the enemy aircraft implements "S" maneuver on a constant height   = 6 km.Its velocity   = 300 m/s, and its lateral acceleration   = 10sgn(sin(   + /2)) within which   = /5.According to the conclusion of Section 3, the target position, namely, the optimum attack position that UCAV should maneuver to, is point  or point  in Figure 12.In addition, assume the UCAV's velocity   = 300 m/s, and it approaches the enemy aircraft with negative height Δ = −1 km.The initial horizontal distance between them is 20 km, and they are in the same vertical plane.The specific geometry relationship of them is illustrated in Figure 12.
The structure characteristics of the UCAV are shown in Table 1.The initial states of the UCAV and enemy aircraft are shown in Table 2. Table 3 sets the key parameters of the designed control system.In addition, we set the amplitude of elevator, rudder, and aileron of the UCAV as ±20 ∘ .

Simulation Results and Analysis.
In this section, G&C represents the traditional design that separates the guidance loop from the control loop, IGC represents the design of integrated guidance and control, and AIGC represents the IGC design with our proposed adaptive SMC method, which regulates key parameters of the IGC system timely.The simulation results are shown in Figures 13-24.The simulation results of pitch channel show that, no matter what kind of designed control system is applied, when the UCAV arrives at the attack position, its attack angle and pitch rate converge to 0. However, in the case of IGC and AIGC, the values of attack angle and pitch rate are relatively smaller in the whole process of attack occupation.What is more, they change more smoothly compared with that in the case of G&C. Figure 15 further shows that the angle of UCAV's elevator is relatively small in the terminal phase of attack occupation in the case of IGC and AIGC.On the contrary, the angle of elevator is in saturated state for a long time in the case of G&C.Because the response from autopilot lags behind the guidance command in the case of G&C, therefore, from   In the yaw channel, both the sideslip angle and yaw rate of UCAV change in a narrow range in the terminal phase of attack occupation.In addition, the range of sideslip angle is [−1.5 ∘ , 0.5 ∘ ], which further provides proof for the assumption that "lateral force caused by sideslip angle has little effect on the UCAV" in Section 4.2. Figure 19 further implies that during the process of autonomous attack occupation the angle of UCAV's rudder approaches 0 within 1 second, and it is always in the range of amplitude.

Pitch
Figures 21 and 22 indicate that in the roll channel the roll angle and roll rate of UCAV change tempestuously during the process of autonomous attack occupation.This results from the constraints of UCAV's gesture at the attack moment and the limitation of aileron amplitude.On the other hand, both the roll angle and roll rate finally converge to 0 at the attack moment in the case of IGC and AIGC.However, in the case of G&C, UCAV is not in a stable state at the attack moment and cannot fulfill the attack occupation task.As the results in Figures 16, 20, and 24, in the case of AIGC, the sliding mode coefficient  converges to its maximum value.On the contrary, the reaching law coefficient  1 , control gain coefficient  2 , and width of boundary layer  converge to their minimum value.These results are all in accordance with the regulation that we design for the key parameters of AIGC system.In addition, as shown in Figures 15, 19, and 23, in the case of AIGC the angle of elevator, rudder, and aileron of UCAV is relatively smaller than that when we apply G&C or IGC.This is because the key parameters of AIGC system are time-varying during the autonomous attack occupation and the UCAV can maneuver according to the relative movement of enemy aircraft accurately and timely.

Mathematical Problems in Engineering
Let () = (  (),   (),   ())  , and we can obtain the equivalent consumption of control energy: ( According to the simulation results, the parameters of UCAV's performance in the case of three designed control systems are compared in Table 4.In Table 4, compared with G&C and IGC, AIGC system can lead UCAV to the optimum attack position more stably and drive it to meet the gesture constraints more accurately.Namely, the kinematic parameters of UCAV are in a relatively smaller range in the process of autonomous attack occupation, and the attack angle, sideslip angle, roll angle, pitch rate, yaw rate, and roll rate of UCAV converge to 0 rapidly and accurately.What is more, AIGC system can economize control energy and weaken the chattering in the process of autonomous attack occupation.In our simulation, it economizes nearly 25% of the equivalent consumption of control energy  compared with the G&C and IGC system.
In the case of AIGC, the trajectory of autonomous attack occupation in three-dimensional space is shown in Figure 25.

Conclusions
In this paper, we have presented an approach for air-to-air autonomous attack occupation of UCAV, which calculates the optimum attack position and guides UCAV to it accurately with gesture constraints.The two main conclusions of this work are the following: (1) considering the requirement of precision, attack zone and no-escape zone can be calculated based on pattern search algorithm, and the optimum attack position can be designed by combining the calculated noescape zone with the practicality of air combat; (2) a novel adaptive SMC method is used to design the IGC system of UCAV, which can guide the UCAV to the optimum attack position accurately with gesture constraints at the attack moment.In the simulation, we demonstrate that compared with traditional methods the proposed approach can drive Mathematical Problems in Engineering the UCAV to fulfill the task of air-to-air autonomous attack occupation effectively with less consumption of control energy.
We believe that the approach of air-to-air autonomous attack occupation for UCAV in this work can inspire more new control systems of UCAV for the task of autonomous attack.For future work, we can explore not only suitable control systems for single UCAV in the task of autonomous attack but also new control and decision systems for multi-UCAVs.

Figure 3 :
Figure 3: Kinematic relation of UCAV and the target in vertical plane.

Figure 4 :
Figure 4: Diagram of pattern search algorithm.

Figure 5 :Figure 6 :
Figure 5: Flowchart of pattern search algorithm for calculating attack distance.

Figure 9 :
Figure 9: The optimum attack position of UCAV.

Figure 10 :
Figure 10: Diagram of guide process of UCAV.

Figure 11 :
Figure 11: Regulating law of main parameters in the designed control system.

Figure 16 :
Figure 16: Key parameters in pitch channel.

Table 1 :
Structure characteristics of the UCAV.
the simulation results of pitch channel we can see that our proposed AIGC system can overcome the shortage of timelag in G&C system, and it performs better than the simple IGC system.

Table 2 :
Initial state of UCAV and the enemy aircraft.

Table 3 :
Controller parameters of the UCAV.
51)   ,   ,   ,   ,   , and   , respectively, denote attack angle, sideslip angle, roll angle, pitch rate, yaw rate, and roll rate of the UCAV at the moment that it arrives at the optimum attack position.|| max , || max , || max , |  | max , |  | max , and |  | max , respectively, denote the maximum absolute value of UCAV's attack angle, sideslip angle, roll angle, pitch rate, yaw rate, and roll rate during the process of autonomous attack occupation.