Abstract

A nonsingular fast terminal sliding mode guidance law with an impact angle constraint is proposed to solve the problem of missile guidance accuracy and impact angle constraint for maneuvering targets. Aiming at the singularity problem of the terminal sliding mode, a fast terminal sliding mode surface with finite-time convergence and impact angle constraint is designed based on fixed-time convergence and piecewise sliding mode theory. To weaken chattering and suppress interference, a second-order sliding mode supertwisting algorithm is improved. By designing the parameter adaptive law, an adaptive smooth supertwisting algorithm is designed. This algorithm can effectively weaken chattering without knowing the upper bound information of interference, and it converges faster. Based on the proposed adaptive supertwisting algorithm and the sliding mode surface, a guidance law with the impact angle constraint is designed. The global finite-time convergence of the guidance law is proved by constructing the Lyapunov function. The simulation results verify the effectiveness of the proposed guidance law, and compared with the existing terminal sliding mode guidance laws, the proposed guidance law has higher guidance precision and angle constraint accuracy.

1. Introduction

In modern warfare, many missiles (such as some antiship missiles, antitank missiles, and air defense missiles) need to hit the target with certain impact angles to increase the damage effectiveness of the warheads. Therefore, the impact angle constraint is a problem that needs to be considered in the design of the guidance law [1].

Sliding mode control is widely used in the design of the guidance law because of its invariability to interference in the sliding mode. In [2], the line-of-sight (LOS) angular velocity and impact angle constraint have been used as the sliding surface, and the sliding mode control is applied to design the guidance law with the impact angle constraint. In [3], the adaptive exponential approach law has been used to design the sliding mode guidance law, which increases the adaptability and dynamic performance of the guidance law. However, none of these methods has finite-time convergence. For the finite-time control problem, a finite-time convergence guidance law with the impact angle constraint is designed based on the terminal sliding mode control in [46]. In [7, 8], linear terms are added to the terminal sliding mode surface to further speed up the convergence of system states. But the negative exponential term of the state quantity in the terminal sliding mode control law will cause the singular problem. For the singular problem, a nonsingular terminal sliding surface is improved to avoid the singularity problem in [911], and the corresponding guidance law is designed. However, the proposed guidance laws cannot guarantee strictly finite-time convergence of the sliding mode surface, and there are nonconvergence factors, so the convergence rate will be reduced. In [12], the problem of nonstrict convergence of the sliding surface is studied, and a nonsingular terminal sliding surface with strictly finite-time convergence is proposed. However, the sliding surface function is not smooth, the system can only converge to a bounded region, and the specific range cannot be given.

For the disturbance problems such as target maneuvering and system disturbance, there are currently three methods for processing most documents: (1) designing the disturbance observer to estimate the disturbance in real time and online [13, 14]; (2) designing the adaptive law to estimate the upper bound of disturbance [15, 16]; and (3) using the robustness of the sliding mode control to resist interference. These methods need to introduce symbolic function terms which will make the control quantity discontinuous and easy to cause the chattering phenomenon. Most literature studies smooth the symbol terms to reduce chattering. However, at the same time, they also change the inherent structure of the sliding mode control and weaken the robustness of the sliding mode control system. For the chattering problem, a second-order sliding mode supertwisting algorithm is proposed in [17]. It has the advantages of simple form, avoiding chattering, and strong robustness. However, the control law of the supertwisting algorithm is not smooth, the parameter selection needs to know the upper boundary information of system disturbance, and the convergence speed is slow when the system states are far from the equilibrium point.

In order to solve the above problems, this paper improves an adaptive smooth supertwisting algorithm, which solves the problems of slow convergence speed and the unsmooth control law of the traditional supertwisting algorithm and greatly weakens the chattering problem of the sliding mode control. At the same time, the parameter adaptive law is designed against the disturbance without knowing the upper bound information of the disturbance. Based on the idea of fixed-time convergence and piecewise sliding surface, a nonsingular fast terminal sliding surface with the impact angle constraint is designed. A nonsingular fast terminal sliding mode guidance law with the impact angle constraint is proposed based on the adaptive supertwisting algorithm. The global finite-time convergence is proved by constructing the Lyapunov function. Finally, the effectiveness and superiority of the guidance law are verified by simulation experiments.

2. Preparatory Knowledge

2.1. Relative Dynamics between the Missile and the Target

In the inertial coordinate system, the relative motion relationship between the missile and the target is established as shown in Figure 1. M and T represent the missile and the target, respectively. is the relative distance between the missile and the target, and is the LOS angle. and are velocities of the missile and the target, respectively, and and are track angles of the missile and the target, respectively.

According to the relative motion relationship of the missile and the target, the relative motion equations of the missile and the target can be obtained as follows:

Differentiating with respect to time gives

The impact angle is the angle between the missile and the target velocity vector at the time of guidance terminal, and the impact angle constraint problem can be transformed into the terminal LOS angle constraint problem [38]. Therefore, the state equation of the guidance system with the impact angle constraint can be obtained based on (2) as follows:withwhere is the desired terminal LOS angle; can be regarded as the total disturbance of the system.

2.2. Related Lemma

For the convenience of analysis and proof, the following lemmas are introduced.

Lemma 1 (see [18]). Assume that there is a smooth function defined on the neighborhood of the origin, and a1 > 0 and 0 < b1 < 1, ; then, the origin of the system is finite-time stable, and the convergence time satisfies

Lemma 2 (see [19]). Assume that Lyapunov function satisfies , and a1 > 0, a2 > 0, and 0 < b1 < 1; then, the system can converge to the origin in finite time, and the convergence time satisfies

Lemma 3 (see [20, 21]). For the nonlinear system , if a1 > 0, a2 > 0, 0 < b1 < 1, and b2 > 1, the system is stable in finite time, and the convergence time satisfies

In addition, if the system has a small disturbance, that is, and is a small positive number, the system can converge to the neighborhood of the origin in finite time, and the convergence time satisfies

3. Adaptive Fast Supertwisting Algorithm

For the following first-order system,where is the system state, is the input, and is the disturbance, the supertwisting algorithm can be expressed as follows:

The supertwisting algorithm can greatly reduce chattering and has strong robustness and high precision control performance [17]. However, the supertwisting algorithm has the following disadvantages: (1) the control law is a continuous function, but not a smooth function, which will affect the control performance; (2) the selection of control parameters needs to know the upper bound information of the system disturbances; and (3) when the system states are far from the equilibrium point, the convergence speed is slow. In view of the above shortcomings, this paper speeds up the convergence of the algorithm by adding linear terms to the algorithm. And the adaptive law does not need the information of the interference. The improved adaptive supertwisting algorithm can be expressed as follows:where .

The parameter adaptive law is designed as follows:where , , , and .

Substituting (11) into (9),

Remark 1. It can be seen from (13) that when the system state is far away from the equilibrium point, the linear term y in (11) will accelerate the convergence rate of the system. When the system is close to the equilibrium point, the nonlinear term plays an important role in accelerating the convergence rate of the system. Therefore, compared with the traditional supertwisting algorithm, the improved adaptive supertwisting algorithm (11) has a faster convergence speed.

Remark 2. Due to the measurement noise of the system, the state of the system cannot reach the equilibrium point completely. In order to avoid the parameter increasing to infinity, the term is added to the adaptive law to avoid the problem of overestimation [22].

Remark 3. It is obvious that . So, it can be seen from (12) that when , m1 and m2 will gradually increase, making the system state to converge. When the system state converges to , m1 and m2 will decrease gradually. If m1 and m2 decrease to the point where the interference cannot be eliminated, the system state will deviate from . At this point, m1 and m2 will gradually increase under the effect of the adaptive law, making the system state converge to . Repeat the previous process; m1 and m2 will gradually decrease. Therefore, m1 and m2 are globally bounded.
For the total disturbance of the system, the following assumption can be made:

Assumption 1. is bounded, , is nondifferentiable, and is differentiable; they satisfy

Theorem 1. Under Assumption 1, the existence of makes the system state converge in finite time when and .

Proof. Define a new state vector asFrom (14), it can be seen that the existence of and makes the following equation valid:Differentiating (15) with respect to time giveswithDefine the following Lyapunov function:withIt is easy to prove that P is a positive definite matrix; then, V1 is unbounded radially, i.e.,Differentiating with respect to time giveswithIf we defineit can be proved that is a positive definite matrix, and .
According to (21), we can getAccording to Lemma 2, z is finite-time convergent. The proof is complete.

Theorem 2. Under the control law (11) and the control parameters which satisfy (12), the system state can converge to in finite time.

Proof. According to (12), . If , then from Theorem 1, the system is finite-time convergent. If , define the Lyapunov function aswhere .
Differentiating with respect to time giveswithAccording to (26), we can getWhen , combining with (12) giveswithIt is known that from , soAccording to Lemma 1, can converge in finite time, y can converge to , and the convergence time satisfiesWhen , if m1 and m2 decrease to the point where the interference cannot be eliminated, the system state will deviate from . In this case, m1 and m2 will increase again under the effect of the adaptive law, making the system state converge to .

4. Guidance Law Design

The terminal sliding mode control adopts the nonlinear function as the sliding mode surface, which can make the system states converge in finite time, but the method has singular problems. In order to avoid singular problems, based on the piecewise sliding surface [12] and Lemma 3, a nonsingular fast terminal sliding surface is designed aswithwhere , , , , , and and are positive odd numbers.

Differentiating with respect to time giveswithwhere is a function of and satisfies the following conditions:(1) is a smooth function in with the same sign as (2)(3), and in

Remark 4. Condition (1) can ensure that is a continuous bounded function and eliminates singular problems, and when the system reaches the sliding surface , and are always with different signs, ensuring that the system state is convergent. Condition (2) ensures that the sliding surface s is a continuous function. Condition (3) guarantees that is bounded in , and is a continuous function, so s is a smooth function.
According to the above conditions, this paper selects function as follows:where and .
Substituting (3) into (38) yieldsThe equivalent guidance law is designed asSubstituting (42) into (41) givesIn order to counteract the disturbance, suppress chattering, and accelerate the convergence speed of the sliding surface, based on the adaptive smooth, fast supertwisting algorithm proposed in the second section, an auxiliary guidance law is designed asThe parameter adaptive law is designed asCombining with (42) and (44), we can design a nonsingular terminal sliding mode guidance law with the impact angle constraint based on the adaptive supertwisting algorithm asFor the convenience of description, the design guidance law (46) is abbreviated as ASNTSMG.

5. Simulation Analysis

In order to test the performance of the designed guidance law, ASNTSMG, this section conducts simulation analysis based on ballistic simulation in different scenarios. The initial position of the missile is (0 m, 0 m), and the initial position of the target is (1000 m, 5000 m). The missile’s velocity is  = 500 m/s, and the target’s velocity is  = 250 m/s. The acceleration of gravity is  = 9.8 m/s2, and the maximum available overload of the missile is 20 g. The parameters of ASNTSMG are set as follows: , , , , , , and .

In order to verify the superiority of the designed guidance law, this section also carries out the nonsingular fast terminal sliding mode guidance law (NFTSMG) proposed in [23] and the second-order nonsingular terminal sliding mode guidance law (SONTSMG) proposed in [24] to perform a comparative simulation. The expression of NFTSMG is

The parameters are set as follows: , , , , , and .

The expression of SONTSMG is

The parameters are set as follows: , , , , and .

The average overload (unit: g) is introduced to evaluate the energy consumption in the process of guidance, which is defined as follows:where K is the total number of simulation steps.

Case 1. Attack moving target with different impact angle constraints: set as 20°, 30°, 40°, and 50°, respectively, and θm0 = 45°. The target makes sinusoidal maneuver, and its acceleration is at = 30sin (πt/5) m/s2, and θt0 = 150°. The simulation results are shown in Figure 2.
It can be seen from Figures 2(a) and 2(b) that ASNTSMG can effectively intercept the target with different impact angle constraints. The miss distances are 0.374 m, 0.428 m, 0.408 m, and 0.479 m, respectively. This method can hit the target accurately. It can be seen from Figures 2(c)2(e) that the sliding surface and the LOS angular rate can converge to zero in finite time, and LOS angle can effectively converge to the expected impact angle. With the increase of , the convergence time increases. This is because the larger is, the larger the initial LOS angle deviation will be, and the convergence time is related to the initial value, which leads to the corresponding growth of the convergence time.
Figure 2(f) shows the overload curve of the missile, which is saturated in the early stage, and the larger θm0, the longer the saturation time, which is mainly due to the larger overload needed in the earlier stage, which makes the missile meet the requirements of angle constraint and guidance accuracy. When and approach the expected values, the overload gradually approaches zero, which ensures that the missile has sufficient overload margin to deal with other unknown disturbances in the later stage of guidance.

Case 2. Comparative simulation of ASNTSMG, NFTSMG, and SONTSMG: the relevant initial parameters are set to θm0 = 45°, , and θt0 = 180°. The movement of the target is set as follows:(1)Cosine motion: at = 30cos (πt/5) m/s2(2)Square wave motion: at = 30sgn (sin (πt/5)) m/s2The simulation results are shown in Figures 36 and Table 1.
Figure 3 shows the trajectories of the missile and the target. It can be seen that the missile can track and intercept the target under the three guidance laws. Compared with NFTSMG and SONTSMG, the trajectory of ASNTSMG is relatively smooth, indicating that its attack time is relatively short, which can be verified by Table 1. Figure 4 shows the LOS angle curve. All three guidance laws can make the LOS angle gradually converge to the expected angle. ASNTSMG can make the LOS angle converge to the expected angle more quickly. NFTSMG adopts the robustness of the sliding mode control to cancel the disturbance of target maneuver, so it can only make the system states converge to the neighborhood of origin in finite time. NFTSMG and SONTSMG adopt the traditional nonsingular terminal sliding surface, which has the nonstrict finite-time convergence problem, and the nonconvergence factor slows down the convergence speed. Figure 5 shows the LOS angular rate curve. Under the three guidance laws, the LOS angular rates converge to zero in finite time. ASNTSMG has smaller convergence error and faster convergence speed. Figure 6 shows the overload curve of the missile. Due to the finite-time convergence of LOS angle and angular rate, the missile needs a large overload in the early stage of guidance. Therefore, the overloads of three guidance laws are all saturated in the early stage. When the LOS angle and the angular rate converge, the missile overload gradually converges to zero in the later stage. And the convergence speed of ASNTSMG is faster.
Table 1 shows the simulation results of attack time, miss distance, LOS angle error, and average overload under the three guidance laws. It can be seen that compared with NFTSMG and SONTSMG, ASNTSMG has smaller attack time, miss distance, terminal LOS angle error, and average overload, so ASNTSMG has better guidance performance.
According to the analysis of the simulation results of two cases, ASNTSMG can hit the target precisely with the expected impact angle under different expected LOS angles and target maneuvering conditions. Compared with the existing guidance laws NFTSMG and SONTSMG, ASNTSMG can effectively attack the target with less impact angle error, miss distance, and energy consumption, which verifies the effectiveness and superiority of ASNTSMG.

6. Conclusion

In this paper, a nonsingular fast terminal sliding mode guidance law is proposed to solve the problem of guidance accuracy and impact angle constraint. Through theoretical analysis and simulation verification, the following conclusions can be obtained:(1)The proposed adaptive smooth supertwisting algorithm can effectively counteract the disturbance of the system and accelerate the convergence speed of the system without knowing the upper bound of the disturbance.(2)The designed nonsingular terminal sliding mode surface can realize the fast finite-time convergence of the system states and ensure the impact angle constraint and guidance accuracy requirements.(3)This guidance law can attack the target precisely under the conditions of different expected LOS angles and target maneuvers. Compared with the existing nonsingular fast terminal sliding mode guidance law and second-order nonsingular terminal sliding mode guidance law, this law has higher guidance accuracy and angle constraint accuracy and consumes less energy.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.