A guidance problem for impact time and angle control applicable to cooperative attack is considered based on the sliding mode control. In order to satisfy the impact angle constraint, a line-of-sight rate polynomial function is introduced with four tuning parameters. And the time-to-go derivative with respect to a downrange orientation is derived to minimize the impact time error. Then the sliding mode control surface with impact time and angle constraints is constructed using nonlinear engagement dynamics to provide an accurate solution. The proposed guidance law is easily extended to a nonmaneuvering target using the predicted interception point. Numerical simulations are performed to verify the effectiveness of the proposed guidance law for different engagement scenarios.
1. Introduction
In modern warfare, the objectives of guidance laws are not only limited to the interception of targets with zero miss distance. To be effective, guidance laws also account for terminal constraints such as the impact time and impact angle [1]. In some applications, it is desirable to control the heading angle of missiles such as antitank missiles and antiship missiles to increase warhead effectiveness and lethality. The impact time constraint is important when performing salvo attacks to enhance the survivability of the missiles against advanced self-defense systems of battleships, such as close-in weapon systems (CIWS) [2]. Typical warships have two to four CIWS to cover a certain defensive zone [3]. A simultaneous attack of multiple missiles with impact time and angle control is regarded as a cost-effective and efficient cooperative attack strategy.
Guidance laws considering terminal angle constraint have been widely studied in the past decades. One of the initial efforts in this research field was the proposal of a suboptimal guidance law for reentry vehicles with impact angle control by Kim and Grider [4]. Subsequently, Kim et al. [5] derived a biased proportional navigation guidance law which included an additional time-varying bias term to impose a desired attitude angle. A generalized form of the optimal guidance law considering the terminal angle constraint was studied in [6], in which two time-to-go calculation methods were proposed to account for path curvature. In addition to using optimal control theory, sliding mode control was also used to achieve terminal angle control [7–9], which was mostly based on nonsingular terminal sliding mode [7], second-order sliding mode [8], and finite-time convergent sliding mode control [9].
The design of impact time control guidance laws was also an active research area. Jeon et al. [10] first derived a suboptimal guidance law to control the impact time against a stationary target, which combined the proportional navigation guidance command and the feedback of the impact time error. Motivated by the idea of ITCG, Zhao and Zhou [11] designed centralized and distributed coordination algorithms for multimissile salvo attacks. Sliding mode control was also used to design impact time control guidance laws. In [12], the switching surface was defined as a combination of the impact time error and the line-of-sight (LOS) rate to meet the requirement of the impact time and zero miss distance. In [13], another nonsingular sliding mode guidance law was proposed for the impact time problem, in which a positive continuous nonlinear function of the lead angle was introduced to avoid the singularity of the guidance command. In addition, the impact time control guidance problems were also extended to three dimensions [14, 15].
Though the issues of impact time and angle control have been studied separately, only a few papers consider both constraints simultaneously. Bokyung and Youdan [16] first combined an impact angle guidance law and an impact time control guidance law to enhance high survivability and kill probabilities. Another method of accounting for both constraints was to derive an impact angle biased proportional guidance law with an additional command for the impact time [17]. Harl and Balakrishnan solved this impact time and angle control problem by shaping the LOS rate and tracking it with second-order sliding mode control [18]. The LOS rate shaping profile was established as a fourth-order polynomial with five tuning parameters to satisfy both impact time and angle constraints. However, the impact time control parameters had to be calculated using an offline complex algorithm for each engagement scenario. Recently, a new strategy was presented based on switching between an impact time control guidance law and an impact angle guidance law [19]. In this scheme, the guidance law corrected the impact time in its first step and then switched to an angle-control guidance law when the time error was reduced to an acceptable value. The potential drawback of this guidance law is that its acceleration command may jump suddenly at the switching instant. An abrupt command change is unfavorable for the missile because it may cause autopilot instability due to the actuator’s slew rate limit [20].
In this paper, a Lyapunov-based impact time and angle control guidance (ITACG) law is proposed. The guidance law is initially developed against a stationary target and then extended to a nonmaneuvering target. The switching surface is defined as a combination of the impact time error and the LOS rate error for satisfying the impact time and angle requirements. The effectiveness of the proposed guidance law is illustrated by numerical simulations under different scenarios such as the many-to-one engagement and the presence of autopilot delay. Simulations highlight the viability and robustness of the guidance law. Compared with other methods in the literatures [18, 19], the proposed guidance law does not require the offline algorithm to calculate the impact time parameter and there is no abrupt jump in the acceleration command during the engagement.
The remainder of this paper is organized as follows. Section 2 defines the problem that is considered in this work. Section 3 derives the proposed impact time and angle control guidance law, followed by numerical simulations to evaluate its performance in Section 4. Finally, some conclusions are given in Section 5.
2. Problem Statement
To define the problem, we first consider a two-dimensional homing guidance geometry for a stationary target as shown in Figure 1. The initial positions of the missile and the target are denoted as x0,y0 and xT,yT. r, θ, λ, σ denote range-to-go, heading angle, LOS angle, and leading angle, respectively. A is the missile acceleration, which is perpendicular to the constant missile velocity V.
Missile/target engagement geometry.
The equations of motion for the missile are given by(1)dxdt=Vcosθ,dydt=Vsinθ,dθdt=AVand its boundary conditions are given as(2)xt0=x0,yt0=x0,θt0=θ0,xtf=xf,ytf=yf,θtf=θd,tf=Td,where the subscripts 0 and f denote the initial and terminal instant. The boundary conditions require the missile to intercept the target at the impact time Td and the impact angle θd.
The engagement kinematics are expressed from Figure 1 as (3)dλdt=-Vsinθ-λr,drdt=-Vcosθ-λ.For the impact time and angle control problem, it is more convenient to have a set of equations with respect to downrange x as the independent variable instead of time t. This is because the final downrange always corresponds to the target’s position xT while the final time tf varies with the choice of impact time. The conversion between time t and downrange x can be performed by using some variable α as follows [18]:(4)dαdx=dαdtdtdx.Under this conversion, (1) can be rewritten as(5)t′=dtdx=1Vcosθ,y′=dydx=tanθ,θ′=dθdx=AV2cosθ.In (5), the “prime” superscript denotes that the new differentiation is performed with respect to downrange x. The same conversion process can be also applied to the boundary conditions in (2) and the engagement kinematics in (3), leading to(6)yx0=y0,θx0=θ0,tx0=0,yxT=yT,θxT=θd,txT=Td,λ′x=-sinθ-λrcosθ,r′x=-cosθ-λcosθ.
3. Impact Time and Angle Control Guidance3.1. Desired LOS Angle and Angular Rate Shaping
In this subsection, motivated by the study in [18], an approach for shaping the LOS rate profile is derived to achieve impact angle control. The desired LOS angle and angular rate can be defined as polynomial functions with respect to downrange x, which are expressed as(7)λdx=ax3xT3+bx2xT2+cxxT+d,(8)λd′x=3ax2xT3+2bxxT2+c1xT,where λdx and λd′x are the desired LOS angle and LOS rate, respectively. The symbols a, b, c, and d represent four constant coefficients, which can be solved by using boundary conditions.
The first two boundary conditions are specified by the initial conditions as(9)λdx0=λ0=tan-1yT-y0xT-x0,λd′x0=λ0′=-sinθ0-λ0cosθ0x0-xT2+y0-yT2.To guarantee a hit, the desired LOS rate should be zero at the downrange corresponding to the target location. So the third boundary condition can be expressed as(10)λd′xT=0.When the target is stationary or its speed is far less than the missile velocity, the final heading angle will be equal to the LOS angle since the desired final LOS rate is zero. Thus, another terminal boundary condition that considers the impact angle control is taken as(11)λdxT=θd,where θd is the desired terminal angle.
By applying (9), (10), and (11), we can obtain(12)ax03xT3+bx02xT2+cx0xT+d=λ0,3ax02xT3+2bx0xT2+c1xT=λ0′,3axT+2bxT+cxT=0,a+b+c+d=θd.From (12), the four tuning parameters can be calculated as(13)a=xT3xT-x02λ0′-2xT3θd-λ0xT-x03,b=3xT2x0+xTθd-λ0xT-x03-xT2x0+2xTxT-x02λ0′,c=-6x0xT2θd-λ0xT-x03+xT22x0+xTxT-x02λ0′,d=θdx023xT-x0+λ0xT2xT-3x0xT-x03-x0xT2xT-x02λ0′.The parameters in (13), which are a function of the impact angle θd, ensure that the missile intercepts the target at the desired terminal angle.
3.2. Estimation of Time-to-Go
Accurate estimation of time-to-go is vital to control the impact time. A poor estimation of time-to-go not only severely degrades guidance performance but also makes the missile off the desired trajectory. If the missile flies straight along LOS, that is, θ=λ, then the time-to-go estimation should be r/V. When θ-λ is not zero, but small, the time-to-go estimation can be approximated as [21](14)tgo=rV1+σ222N-1,where N denotes the effective navigation constant.
For the impact time and angle control, the time-to-go estimation in (14) may not be adequate because its time-to-go error only depends on the initial heading error. But the effect of the terminal angle constraint is not considered. In [6], the trajectory considering the impact angle constraint can be approximated by a third-order polynomial of range-to-go. Thus the time-to-go estimation can be calculated as the range over the average velocity, which gives(15)tgo≈rV1+σ2+σf215-σσf30,where σ and σf are defined as(16)σ=θ-λ,σf=θd-λ.By differentiating (15), the derivative of time-to-go tgo′ with respect to downrange x can be calculated as(17)tgo′=r′V1+σ2+σf215-σσf30+rV2σσ′+2σfσf′15-σ′σf+σσf′30.By differentiating (16), we can also obtain the derivatives of σ and σf as(18)σ′=θ′-λ′=AV2cosθ-λ′,σf′=-λ′.Substituting (18) into (17) yields(19)tgo′=r′V1+σ2+σf215-σσf30+rAV3cosθ4σ-σf30-rλ′σ+σf10V.Equation (19) can be rewritten in a simpler form as(20)tgo′=K1+K2A,where(21)K1=r′V21+σ2+σf215-σσf30-rλ′10Vσ+σf,K2=rV3cosθ2σ15-σf30.
3.3. Impact Time and Angle Control Guidance Law
The objective of the guidance law is to enable the missile to intercept a target with the desired impact time and impact angle. According to the above requirements, the switching surface with respect to downrange is chosen as(22)s=eλ′+etsigneλ′,where (23)eλ′=λd′x-λ′x,(24)telap+tgoet=Td-telap+tgo.In (23), λ′ is the real-time LOS rate and λd′ is the desired shaping LOS rate in (7). In (23), telap is the time elapsed after the launch of the missile. Note that the desired flight time Td should be larger than the initial time-to-go tgox0, and then et≥0.
Proposition 1.
The sliding surface in (22) can help to achieve both impact time and angle constraints.
Proof.
The dynamics of the sliding surface in (22) are given by s=0; that is,(25)eλ′+etsigneλ′=0⟹eλ′=-etsigneλ′.When the LOS rate error eλ′ is not equal to zero, we can get(26)eλ′=-et≤0.But eλ′≥0 for all eλ′, so the only term is et=0,eλ′=0. In other words, the sliding mode s=0 cannot occur until et=0 and eλ′=0. et=0 means that the estimated interception time is equal to the desired impact time. eλ′=0 means that the real-time LOS rate is equal to the designed one in (7) at the interception point; that is, λ′xT=λd′xT. According to the terminal boundary conditions in (10) and (11), we can get λ′xT=λd′xT=0 and λxT=λdxT=θd. Thus the angle constraint can be achieved.
Now consider the case et≠0 and eλ′=0 at the start of the engagement. Differentiating (23) yields(27)eλ′′=λd′′x-λ′′x=cosθ-λ+tanθsinθ-λrV2cos2θA+sin2θ-2λr2cos2θ+6axxT3+2bxT2=cosλrV2cos3θA+sin2θ-2λr2cos2θ+6axxT3+2bxT2.Note that eλ′′ is not identically equal to zero with the change of range-to-go r, lateral acceleration A, heading angle θ, and LOS angle λ. Hence, eλ′ is a time-varying value and eλ′=0 is not an attractor while s≠0. Once eλ′ is not equal to zero, we can use the sliding surface above to control et and eλ′.
To design lateral acceleration A for impact time and angle control, let us consider the downrange derivative of the switching surface.
Case 1.
When LOS rate is less than the desired LOS rate, that is, eλ′>0, then from (22)(28)s=eλ′+et.Differentiating (28) yields(29)s′x=eλ′′x+et′x=λd′′x-λ′′x-telap′+tgo′=cosλrV2cos3θA+2r′λ′r+6axxT3+2bxT2-1Vcosθ-K1+K2A.Based on (29), one can design the lateral acceleration using the Lyapunov-based nonlinear control theory as the sum of an equivalent control and a discontinuous control. Consider(30)A=aMeq+aMdis.The equivalent component aMeq is obtained by setting s′=0, which gives(31)aMeq=1/Vcosθ+K1-2r′λ′/r+6ax/xT3+2b/xT2cosλ/rV2cos3θ-K2.Once the system reaches the sliding mode, the equivalent control component keeps the system on the sliding mode. In addition, the discontinuous component guarantees that the switching surface s converges to zero, which is designed in the form of(32)aMdis=-Msignssigncosλ/rV2cos3θ-K2,where M denotes the control magnitude, M>0. By increasing the control magnitude M, the time interval to reach sliding mode is reduced. However, M should be chosen appropriately to ensure that the acceleration constraint is not violated.
Consider the Lyapunov function candidate to be VL=1/2s2. Differentiating VL and substituting (29) and (30) into the resulting expression gives(33)VL′=ss′=-MscosλrV2cos3θ-K2,which is negative definite if(34)cosλrV2cos3θ-K2≠0.Under these conditions, it is guaranteed to achieve the sliding mode. Hence the guidance law is(35)A=aMeq+aMdis=1/Vcosθ+K1-2r′λ′/r+6ax/xT3+2b/xT2cosλ/rV2cos3θ-K2+-Msignssigncosλ/rV2cos3θ-K2.Note that if (34) is not satisfied, then VL′=0, and the Lyapunov condition will not be satisfied. If (34) is violated, then (36)cosλrV2cos3θ-K2=0.After substituting (21) and neglecting smaller terms of K2, we can get(37)15V2r2=cos2θcosλσ.The left hand side (LHS) of (37) remains positive during the whole engagement. Here we consider the range of θ and λ as θ<π/2, λ<π/2. If the heading angle error σ≤0 then σcos2θ/cosλ≤0, and the condition (34) is not violated; If σ>0, that is, θ>λ, the right hand side is greater than zero. According to the terminal boundary conditions in (10) and (11), that is, λ′xT=0 and λxT=θd, the heading angle error σ converges to zero asymptotically. The highly increasing time-varying term 1/r2 on LHS forces the system cross (37) and therefore Lyapunov stability criteria is still satisfied.
Case 2.
When LOS rate is greater than the desired LOS rate, that is, eλ′<0, then from (22)(38)s=eλ′-et.Differentiating (38) yields(39)s′=eλ′′-et′=λd′′x-λ′′x+telap′+tgo′=cosλrV2cos3θA+2r′λ′r+6axxT3+2bxT2+1Vcosθ+K1+K2A.Also, the lateral acceleration of the missile is designed to be the sum of an equivalent control and a discontinuous control. Consider(40)A=aMeq+aMdis.By setting s′=0, the equivalent component aMeq is obtained as(41)aMeq=-1/Vcosθ+K1+2r′λ′/r+6ax/xT3+2b/xT2cosλ/rV2cos3θ+K2.The discontinuous controller, which drives the system towards the switching surface, is of the form(42)aMdis=-Msignssigncosλ/rV2cos3θ+K2,where M denotes the control magnitude, M>0. Now we consider the Lyapunov function candidate to be VL=1/2s2. Differentiating VL and substituting (39) and (40) into the resulting expression gives(43)VL′=ss′=-MscosλrV2cos3θ+K2≤0,which is negative definite if(44)cosλrV2cos3θ+K2≠0.Under these conditions, it is guaranteed to achieve the sliding mode. Hence the guidance law is(45)A=aMeq+aMdis=-1/Vcosθ+K1+2r′λ′/r+6ax/xT3+2b/xT2cosλ/rV2cos3θ+K2+-Msignssigncosλ/rV2cos3θ+K2.Note that if (44) is not satisfied, VL′=0, and the Lyapunov condition will not be satisfied. Similarly to the analytical method in Case 1, it can be easily shown that condition (44) is satisfied most of the time during the engagement and hence the impact time and angle criterion can be achieved.
In addition, to reduce chattering due to the application of the discontinuous control in (35) and (45), the discontinuous function signs can be approximated by the sigmoid function [19]. Consider(46)sgmfs=211+exp-s/ε-12,where ε denotes the width of the boundary layer in the sigmoid function. In (46), the function uses a continuous approximation of the discontinuous control when s≤ε. The performance of the sigmoid function as ε is varied is shown in Figure 2. It can be seen that sigmoid function comes closer to the actual signum function for lower values of ε. This method ensures that the proposed guidance law generates a continuous and smooth lateral acceleration.
Performance of the sigmoid function.
3.4. Extension to Nonmaneuvering Targets
The proposed guidance law can be extended to attack a nonmaneuvering target using the concept of predicted interception point (PIP) represented in [22]. The PIP is defined as a virtual point where the missile is expected to intercept the target. The target is assumed to have a speed VT with a heading angle θT, and the virtual position xTPyTP can be calculated online as follows:(47)xTP=xT+VTcosθTtgo,yTP=yT+VTsinθTtgo,where xTP,yTP and xT,yT are the predicted interception point and the current position of the target, respectively, and tgo is the time-to-go estimation defined in (15). Using this method, we can approximate a nonmaneuvering target as a virtual stationary target. This concept can easily be incorporated into the proposed ITACG law by replacing xT,yT with xTP,yTP.
4. Simulation Results
This section presents simulation results of the proposed ITACG law under various engagement scenarios. The proposed guidance law is first simulated in one-to-one and many-to-one scenarios for stationary targets and nonmaneuvering targets. The simulation is then extended to evaluate performance with autopilot lag. For the given engagements, the gain M of the discontinuous component in the proposed guidance law is chosen as M=50. The width of the boundary layer in the sigmoid function (46) is chosen as ε=0.2. For each simulation, it is assumed that the desired guidance time and terminal angle are specified.
4.1. Stationary Target
In this subsection, performance of the proposed guidance law considering a stationary target is given. The simulation conditions of the missile are listed in Table 1. In the first simulation, the results of the PNG and sliding mode control guidance (SMCG) [19] law are also shown for comparison with the proposed ITACG. The SMCG requires values for the tolerance parameters ε1 and ε2 in its guidance law switching logic. Once the impact time error becomes smaller than the tolerance parameters, SMCG switches to the impact angle control guidance law. Here the tolerance parameters are chosen to be ε1=0.3s and ε2=0.05s.
Simulation parameters for stationary target.
x0,y0/m
xT,yT/m
V/(m/s)
θ0/deg
θf/deg
Td/s
(0,500)
(10000,0)
250
30
-60
50
The simulation results for this scenario are presented in Figure 3. Figure 3(a) shows that the missile trajectory using ITACG is similar to that of SMCG. Both of these guidance laws use a longer route to achieve the desired impact time Td compared to PNG, which uses a more direct route but hits the target too early. The terminal time of PNG and SMCG in Figure 3(b) can be observed to be 41.41 s and 50.29 s, respectively. Even though the initial time-to-go of ITACG is about 44.34 s, its final impact time is 49.99 s, which is very close to the desired impact time. Figure 3(b) shows that the sliding surface of ITACG converges to zero quickly, while sliding surface of SMCG is minimized, but never reaches zero. This is because when SMCG switches to the impact angle control guidance law (in other words, the impact time error decreases below the tolerance parameters), it ceases to control the impact time. Figure 3(c) compares the heading angles of the three guidance laws. The heading angle for ITACG and SMCG increases at the beginning of the engagement, while it decreases monotonically for PNG. However, PNG does not achieve the desired impact angle at the finial time, while ITACG does achieve this goal. It also can be observed from Figure 3(d) that ITACG and SMCG employ higher accelerations than PNG because of the impact time and angle control. Due to the way the guidance logic is implemented in SMCG, its lateral acceleration jumps suddenly at the switching instants. In contrast, the lateral acceleration of ITACG is observed to change continuously throughout the engagement.
Simulation results compared with other guidance laws.
Missile trajectory
Time-to-go and occurrence of sliding mode
Heading angle
Missile acceleration
In order to further evaluate the performance of the proposed guidance law, simulations with the impact time of Td=45,50,55s and a fixed impact angle of θd=-60deg are performed. All other simulation conditions are the same as in the previous scenario in Table 1. Figures 4(a) and 4(b) show that the proposed guidance law generates feasible trajectories for satisfying the three impact time constraints. In order to increase the impact time, the missile needs to turn farther away from the target by increasing its acceleration, as shown in Figures 4(c) and 4(d).
Simulation results considering various impact times with a fixed impact angle.
Missile trajectory
Time-to-go and occurrence of sliding mode
Heading angle
Missile acceleration
Next simulations are conducted with a fixed impact time of 50 s and desired impact angles of −30 deg, −50 deg, and −80 deg. All other simulation conditions are the same as in the previous scenario in Table 1. Figures 5(a) and 5(c) show that three different trajectories satisfy each impact angle requirement. Figure 5(b) shows that each time-to-go estimation approaches the same desired time-to-go quickly. Figure 5(d) shows that the different angle constraints require the different lateral accelerations. For each case, the acceleration is large at the beginning of the engagement until the system reaches the sliding mode surface.
Simulation results considering various impact angles with a fixed impact time.
Missile trajectory
Time-to-go and occurrence of sliding mode
Heading angle
Missile acceleration
In the next simulation, a salvo attack engagement scenario is considered with three missiles. Three missiles are assumed to be launched from different positions with different initial heading angles and different impact angles. All three missiles are required to hit the target simultaneously. In this case, each missile has the same speed of 250 m/s, and the target position is still (10000,0) m. Other relevant parameters are given in Table 2. Figure 6(a) illustrates the salvo attack trajectories used by ITACG compared with PNG. It can be seen that the trajectories used by ITACG have a greater curvature than those with PNG. Figure 6(b) shows that the proposed guidance law reduces the impact time dispersion with negligible error. The impact angles assigned to each missile are all achieved at the final time, as shown in Figure 6(c). From Figure 6(d), missiles need larger accelerations at the beginning to enable the system to move to the sliding mode.
Simulation parameters for salvo attack.
Missile 1
Missile 2
Missile 3
Position X/m
0
3000
4000
Position Y/m
500
8000
−8000
Heading angle/deg
10
−40
30
Desired impact angle/deg
−30
−80
80
Desired impact time/s
45
45
45
Simulations results considering a salvo attack against a stationary target.
Missile trajectory
Time-to-go
Heading angle
Missile acceleration
4.2. With Autopilot Dynamics
In this subsection, first-order uncompensated autopilot dynamics are considered with a time constant τm. The autopilot dynamics lag can be expressed as(48)am=1τms+1ac,where ac and am are the acceleration command and the lateral acceleration, respectively.
Note that, in the presence of autopilot lag, the coefficients in (13) are no longer constant. They should be updated as time-varying parameters, which are given as(49)ax=xT3xT-x2λ′-2xT3θd-λxT-x3,bx=3xT2x+xTθd-λxT-x3-xT2x+2xTxT-x2λ′,cx=-6xxT2θd-λxT-x3+xT22x+xTxT-x2λ′,dx=θd3x2xT-x3+λxT3-3xT2xxT-x3-xxT2xT-x2λ′.Simulations are carried out for different initial heading angles with autopilot time constant τm=0.3 s. With the same initial position and other constraints as in Table 1, different launch angles, 0,40 and 80 deg, are considered. Figure 7 shows that the missile successfully intercepts the target at the desired impact time and angle with different launch angles. The impact time error is less than 1×10-2 s for all the cases. And the impact angle error is less than 0.1 deg. This simulation highlights the robustness of the proposed guidance law with respect to uncompensated autopilot delays.
Simulation results considering various launch angles with autopilot dynamics.
Missile trajectory
Time-to-go and occurrence of sliding mode
Heading angle
Missile acceleration
4.3. Nonmaneuvering Target
This subsection demonstrates the effectiveness of the proposed ITACG against a nonmaneuvering target. The target is assumed to move with a velocity of 20 m/s, which is a typical cruise speed for a modern warfare ship. The heading angle of the target is 90 deg, and the initial conditions of the missile and target are given in Table 3. Figure 8 indicates that the proposed ITACG successfully achieves the desired impact time and impact angle when used against a moving target. The impact time error of ITACG is 0.01 s, which is smaller than the 0.29 s error obtained using SMCG. Figure 8(d) shows that the proposed guidance law avoids the abrupt jumps of lateral acceleration.
Simulation parameters against a moving target.
x0,y0/m
θ0/deg
θf/deg
Td/s
xT,yT/m
VT/(m/s)
θT/deg
(0,500)
30
−60
50
(10000,0)
20
90
Simulation results for a nonmaneuvering target.
Missile trajectory
Time-to-go and occurrence of sliding mode
Heading angle
Missile acceleration
A further simulation is performed in a salvo attack scenario. As before, the target is moving with a velocity of 20 m/s and a heading angle of 90 deg. The other simulation conditions are the same as in the previous salvo attack scenario against a stationary target. Figure 9 shows that each missile successfully intercepts the moving target at the desired time with the correct impact angle.
Simulation results considering a salvo attack against a nonmaneuvering target.
Missile trajectory
Time-to-go
Heading angle
Missile acceleration
5. Conclusion
In this paper, a novel guidance law based on sliding mode control is presented to achieve both impact time and angle control. The sliding mode surface is formed using the LOS rate error and the time-to-go error. The LOS rate shaping method is designed to realize the terminal angle requirement and the time-to-go derivative estimation is used to minimize the impact time error. Simulations of various engagement scenarios demonstrate that the proposed guidance law successfully controls the impact time and angle against both stationary and nonmaneuvering targets. Further simulations show that the proposed guidance law also provides excellent performance considering autopilot dynamics. Future studies will consider additional constraints with regard to this guidance law, such as the seeker’s field-of-view limit.
Competing Interests
The authors declare that there are no competing interests regarding the publication of this article.
ZhangY. A.WangX. L.WuH. L.Impact time control guidance law with field of view constraintKimT.-H.LeeC.-H.JeonI.-S.TahkM.-J.Augmented polynomial guidance with impact time and angle constraintsKangS.KimH. J.Differential game missile guidance with impact angle and time constraintsKimM.GriderK. V.Terminal guidance for impact attitude angle constrained flight trajectoriesKimB. S.LeeJ. G.HanH. S.Biased PNG law for impact with angular constraintRyooC.-K.ChoH.TahkM.-J.Optimal guidance laws with terminal impact angle constraintKumarS. R.RaoS.GhoseD.Nonsingular terminal sliding mode guidance with impact angle constraintsHeS.LinD.WangJ.Continuous second-order sliding mode based impact angle guidance lawZhangY.SunM.ChenZ.Finite-time convergent guidance law with impact angle constraint based on sliding-mode controlJeonI.-S.LeeJ.-I.TahkM.-J.Impact-time-control guidance law for anti-ship missilesZhaoS.ZhouR.Cooperative guidance for multimissile salvo attackKumarS. R.GhoseD.Impact time guidance for large heading errors using sliding mode controlChoD.KimH. J.TahkM.-J.Nonsingular sliding mode guidance for impact time controlWeiX.WangY.DongS.LiuL.A three-dimensional cooperative guidance law of multimissile systemZhaoJ.ZhouR.DongZ.Three-dimensional cooperative guidance laws against stationary and maneuvering targetsBokyungJ.YoudanK.Guidance laws for anti-ship missiles using impact angle and impact timeProceedings of the AIAA Guidance, Navigation, and Control Conference and ExhibitAugust 2006Keystone, Colo, USAAmerican Institute of Aeronautics and Astronautics10.2514/6.2006-6432ZhangY. A.MaG. X.LiuA. L.Guidance law with impact time and impact angle constraintsHarlN.BalakrishnanS. N.Impact time and angle guidance with sliding mode controlKumarS. R.GhoseD.Impact time and angle control guidanceProceedings of the AIAA Guidance, Navigation, and Control ConferenceJanuary 2015Kissimmee, Fla, USAAIAALeeC.-H.KimT.-H.TahkM.-J.Interception angle control guidance using proportional navigation with error feedbackJeonI.-S.LeeJ.-I.TahkM.-J.Homing guidance law for cooperative attack of multiple missilesKumarS. R.GhoseD.Sliding mode control based guidance law with impact time constraintsProceedings of the 1st American Control Conference (ACC '13)June 2013576057652-s2.0-84883511588