MISSILE INTERCEPTOR GUIDANCE AND CONTROL USING SECOND ORDER SLIDING

MODES

1
Yuri B. Shtessel, 2Ilya A. Shkolnikov and 3Arie Levant

1
Department of Electrical and Computer Engineering, The University of Alabama in
Huntsville, Huntsville, AL 35899 USA, phone: (256) 824-6164, fax: (256) 824-6803,
[email protected]
2
Z/I Imaging Corporation, Intergraph Company, 230 Business Park Blvd., Madison, AL
35757, USA, [email protected]
3
Applied Mathematics Department, School of Mathematical Sciences, Tel-Aviv
University, Tel-Aviv 69978, Israel, [email protected]

Abstract: Abstract: The missile interceptors integrated guidance and control technology
for achieving the hit-to-kill accuracy against targets performing evasive maneuvers
including spiraling motion is developed on the basis of high (second) order sliding mode
control (HOSM). The integration of guidance and flight control systems is achieved in a
two-loop guidance and flight control system designed in the combined state space of
engagement kinematics and vehicle dynamics. The designed guidance-control system
performance is verified via computer simulations using miniature hypervelocity kinetic
energy endo-atmospheric interceptor planar model. Copyright © 2005 IFAC

Keywords: sliding mode control, missile guidance and control

1. INTRODUCTION follows: Section 2 is dedicated to intercept strategy

The miss distance as an ultimate performance and a guidance law development based on a smooth
criterion of a homing interceptor is crucially second order sliding mode. In Section 3 integration of
dependent on guidance, navigation and control (see HOSM-based guidance and flight control system is
Bezick, Rusnak, & Gray, 1995; Garnell, & East, presented. Simulations are discussed in Section 4.
1977; and Zarchan, 1998) working together in a Conclusions are summarized in Section 5.
closed loop. This situation calls for the integrated
design of all interceptor modules: sensor information 2. INTERCEPT STRATEGY AND SLIDING MODE
processing and necessary data estimation, a homing GUIDANCE
guidance law, and flight control (autopilot). The 2.1 Introduction to Sliding Mode Guidance and
primary goal of this work is to develop the enabling Problem Formulation
guidance and control technology for missile
interceptors to achieve the hit-to-kill accuracy It’s known (see Zarchan, 1998) that proportional
against targets performing evasive maneuvers navigation (PN) homing missile guidance can be
including spiralling motion, using high order sliding derived as the optimal solution in mean-square sense,
mode control (HOSM) developed in works of when there is no target maneuver, and only initial
Floquet, Barbot, & Perruquetti (2003), Fridman heading error should be nullified. In presence of a
(2001, 2003), Levant (1993,1998, 2001, 2003), and step-constant target maneuver augmented
Orlov, Aguilar, & Cadiou (2003). proportional navigation (APN) guidance is optimal
for zero-lag guidance system. However, in reality the
The main concern of the paper is to use minimum
unmodeled series dynamics can cause instability of
possible information in order to achieve the goal the homing loop that yields large miss distance.
(target intercept) in presence of uncertainties and Optimal guidance law (OGL) effectively
disturbances acting in the homing loop. In this work, compensates for a single-lag frame dynamics in
an integrated two-loop HOSM-based guidance and presence of a step-constant target maneuver
flight control system is designed to robustly enforce modifying APN by using information of time-to-go
the target intercept in presence of target maneuvers, and effective time-constant of the guidance system.
atmospheric disturbances, and dynamic uncertainty However, OGL is difficult to derive and implement
of airframe-actuator. The paper is organized as (see Zarchan, 1998).Recent application of sliding
mode control (SMC) (see Utkin 1998) to the homing rD = V r ,
missile guidance (see Brierly and Longchamp, 1990; D
V r = Vλ / r + AT r − sin(λ − γ M )n L ,
2
Babu, Sarma, and Swamy, 1994; Zhou, Mu, and Xu,
D (3)
λ = V λ / r ,
1999; Moon and Kim, 2000; Shkolnikov, Shtessel,
et. al., 2000), resulted in a series of very effective D
algorithms in terms of smaller acceleration advantage Vλ = −V r Vλ / r + AT λ − cos(λ − γ M )n L ,
required for intercept of weaving targets as compared where Vλ = rω λ is a transversal component of
to PN and APN guidance. Though SMC guidance relative velocity in the reference frame rotating with
doesn’t use direct estimate of target acceleration, it’s LOS. It is well known that for a direct hit, it’s
comparable to OGL. Since the guidance law must be
necessary to keep Vr < 0 . It was shown in
smooth, all SMC guidance algorithms include
smoothing procedures that yields partial loss of Shkolnikov, Shtessel, et. al., 2000 that direct hit can
robustness to disturbances and uncertainties. In this be achieved if ω λ = 0 or Vλ = 0 . Another less
work the problem is addresses in the frame of HOSM aggressive hit-to-kill guidance strategy is known
(see Levant 1993,1998, 2001, 2003; Shtessel, ω λ = c 0 / r or Vλ = c0 r , (4)
Shkolnikov, and Brown, 2003). A smooth asymptotic
where c 0 is some constant.
second order sliding mode control algorithm (SOSM)
(see Shtessel, Shkolnikov, and Brown, 2003) is Now, the following guidance task can be formulated:
reassessed and is proved to be SOSM control. Next, Stabilize the system (2) or (3) on the manifold
the guidance algorithm based on a smooth SOSM σ 1 = ω λ = 0 , or σ 1 = Vλ = 0 (5)
control and SOSM controller/autopilot based on or
nonlinear dynamic sliding manifold (see Krupp,
Shkolnikov, & Shtessel, 2000) are used for designing σ 2 = ω λ − c 0 / r . = 0 , or σ 2 = Vλ − c0 r = 0 (6)
an integrated guidance-flight-control system. The by means of the normal acceleration (guidance)
main idea developed in this work is to enforce the command n Lc . This command is usually followed by
suitable closed-loop missile-target engagement means of corresponding aerodynamic surface
kinematics using the input voltage signal to the deflections that is treated as a control function
actuator as control. In this case engagement compensating for the interceptor dynamics.
kinematics and missile dynamics are integrated into
one state space, and the intercept problem is In the work of Shkolnikov, Shtessel, et. al., (2000) a
transformed into the output regulation problem, different approach was proposed. After stabilization
which is addressed using SOSM-based control for of (5) or (6) at zero the resulted compensated
the missile aerodynamic surface actuator. kinematics is derived. On the basis of this
compensated kinematics a corresponding sliding
surface is designed and stabilized by controlling the
2.2 Intercept Strategy: Geometric Approach
aerodynamic surface deflection. Doing this the
One can derive planar missile-target engagement guidance law in terms of n Lc is dissolved in the
kinematics without account for gravity as sliding surface designed on the basis compensated
RD = V , V = V −V ,
R R T M kinematics of the interceptor. Unlike in the work
D (1) Shkolnikov, Shtessel, et. al., (2000), where various
VR = A R , A R = A T − A M .
continuous approximations of SMC were studied for
such that R, VR , A R , VT , VM , A T , A M ∈ ℜ 2 ; the normal acceleration command, n Lc , design, in
R, VR , A R are range to the target and its first and this work we propose to use HOSM-based control.
second time-derivatives respectively; VT , A T are The expected advantages are in increasing robustness
target velocity and acceleration, VM , A M are missile and accuracy of hit-to-kill intercept.
velocity and acceleration. The guidance command is
2.3 Prescribed Sliding Variable Dynamics
missile normal acceleration, n L . In polar coordinate
system the missile-target relative position is Consider a SISO sliding variable dynamics given as
presented by R = (r , λ ) , where r = range along Line- σD = f ( x, t ) + u , (7)
Of-Site (LOS), and λ = LOS angle. The following where σ ( x, t ) ∈ ℜ is the sliding quantity, such that
state model (see Shkolnikov, Shtessel, et. al., 2000) σ = 0 defines the system motion on the sliding
of missile-target engagement process is obtained
surface, u ∈ ℜ1 is a control input that is supposed to
rD = V r , be smooth, and f ( x, t ) is an uncertain differentiable
D
V r = rω λ + AT r − sin(λ − γ M )n L ,
2
(smooth) nonlinear time-varying function, and
D (2)
λ = ω λ , x ∈ ℜ n is a state vector of the system the sliding
ωD = (− 2V ω + A − cos(λ − γ )n ) / r , variable dynamics (7) is considered for. The drift
 λ r λ Tλ M L
term f ( x, t ) is cancelled by means of a special
where we consider ω λ as a commanded output,
observer to be developed further.
missile normal acceleration as a control input, and
projections of target acceleration along and The following compensated σ-dynamics in (7) is
orthogonal to LOS, AT r , AT λ , are considered as considered
unknown bounded disturbances. The system (2) can
be written also as
  xD = −α x (r −1) / r sign( x ) + x , T. The right-hand sides of (8) are uniformly
 1 1 1 1 2
continuous in D0, therefore, due to the continuous
 (r − 2 ) / r
(8)
 xD 2 = −α 2 x1 sign( x1 ), dependence on the right-hand side, with any τ
where x1 = σ . smaller than some sufficiently small τ0 the
trajectories of (9), (10) concentrate in some small
Definition. We call a system finite-time stable
closed disk D1’ ⊂ D0. In their turn the trajectories of
(Bacciotti, & Rosier, 2001), if it is asymptotically
(8) starting from D1’ terminate at O in time T. With
stable with a finite settling time for any solution and
sufficiently small τ1 ≤ τ0, τ ≤ τ1, the trajectories
initial conditions.
terminate in some other small disk D1’’ in time T,
Lemma 1. Let r ≥ 2, α1, α2 > 0. Then the system (8) is D1’’ ⊂ D1’. Let D1 be some disk containing all the
finite-time stable, and the settling time being a trajectories’ segments starting from D1’ with t
continuous and vanishing at the origin function of
varying in the range [0, T], τ ≤ τ1, D1 ⊂ D0.
the initial conditions.
Obviously, with τ ≤ τ1 any trajectory, which starts
Proof. Consider the following Lyapunov function
candidate: from D0, enters D1 in the time T to stay there forever.
Note that system (9), (10) is invariant with respect to
x1
x2 ( r − 2) / r the transformation
V = 2 + α2 z
2 ∫
0
sign( z )dz = r
Gκ: (t, τ, x1, x2) v ( κt, κτ, κ x1, κ x2),
r-1

x 22 for any κ > 0. There exists such κ > 1 that GκD1 ⊂

α 2 x1(2 r − 2 ) / r
r
+ D0. Thus, with τ ≤ κτ1 all trajectories starting from
2 2r − 2
Its derivative is GκD0 enter GκD1 in the time not longer than κT and
∂V then proceed into D1 in time not longer than T to stay
VD = ⋅ xD = there forever. In particular, it is true also for τ ≤ τ1.
∂x 2
Thus a set sequence D1 ⊂ D0 ⊂ GκD0 ⊂ Gκ D0 ⊂ ...
[ ]
 x − α x ( r −1) / r sign x 
( r − 2) / r
sign x1 , x 2 ⋅  =
2 1 1 1
α 2 x1 covering the whole plane is constructed. Hence,
− α x ( r − 2 ) / r
sign x1  with τ ≤ τ1 all trajectories enter D1 in finite time to
 2 1
(r − 2 ) / r + (r −1) / r stay there forever.
α 2 x 2 x1(r − 2 ) / r sign( x1 ) − α 1α 2 x1 −
Let D1 ⊂ { x1, x2 | |x1| ≤ c1, |x2| ≤ c2}, and let τ be
(2 r −3 ) / r
α 2 x 2 x1(r − 2 ) / r sign( x1 ) = −α 1α 2 x1 some arbitrary sampling step. Applying Gκ with κ =
Apply the La Salle theorem. The set x : VD ( x ) = 0 { } τ/τ2 achieve that with the sampling step τ any
r r
trajectory enters the region {x1, x2| |x1| ≤ (c1/τ2 ) τ ,
consists of the axis x1 = 0 . It is easy to see that the r-1 r-1
||x2| ≤ (c2/τ2 ) τ } to stay there forever. Lemma 2 is
only invariant set inside x1 = 0 is the origin proven.

x1 = x 2 = 0 . Thus, the asymptotic convergence of Remark. Since r > 2 can be any real number, Lemma
x1 and x 2 to zero is assured. It is easy to see that 2 provides the first published example of a non-
system (8) is homogeneous with the dilation integer real-sliding order.
r r-1
dκ: (x1, x2) v (κ x1, κ x2) and the negative
2.4 Nonlinear disturbance observer
homogeneity degree -1. Following the results of
Bacciotti, & Rosier, (2001), the asymptotical stability The sliding variable dynamics (7) is sensitive to the
implies here the finite-time stability and the unknown bounded term f (σ ( x(t ), t ), t ) = g (t ) , which
continuity of the settling time function.
is assumed to be smooth with | IgII | ≤ L. In order to
When the variables are sampled at discrete times ti estimate g (t ) the modified 2nd order differentiator
with some sampling step τ > 0 the system (8) take on
(Levant, 2003) is used
the form
sign(z 0 − f ) + z1
2/ 3
zI0 = v0 , v0 = −2 ⋅ L1/ 3 z 0 − σ
(r - 1)/ r
xD1 = - α1 |x1(ti)| sign x1(ti) + x2(ti) (9)
sign(z1 − v0 ) + z 2 (12)
(r - 2)/ r 1/ 2
xD 2 = - α2 |x1(ti)| sign x1(ti) (10) zI1 = v1 , v1 = −1.5 ⋅ L1/ 2 z1 − v0
where the current time satisfies the condition zI 2 = 1.1⋅ L ⋅ sign(z 2 − v1 )
t i ≤ t < t i +1 = t i + τ (11)
and g (t ) = z1 (t ) in a finite time, if the sliding
Lemma 2. Let r ≥ 2, and α1, α2 > 0. Then with τ > 0 variable σ and control u are measured without noise.
all trajectories of the discrete-sampling system (9), Suppose that σ and u are measured with some
(10) converge in finite time into some vicinity of the
r Lebesgue-measurable noises bounded by ε ≥ 0 , then
origin defined by the inequalities | x1| < γ1 τ , |x2| < γ2 it is can be shown that there exists µ > 0 such that
r-1
τ with γ1, γ2 > 0. (Levant, 2003)
Proof. It is proven in Lemma 1 that with τ = 0 the z1 (t ) − g (t ) ≤ µε 2 / 3 (13)
system is finite-time stable. In particular all
trajectories starting in some closed disk D0 centred at 2.5 Smooth High Order Sliding Mode Control
the origin converge to the origin in some finite time
Consider a sliding-variable dynamics (7) with In this paper, we employ smooth second order sliding
smooth f (σ , t ) = g (t ) . The problem is to design mode control (14) with r = 3 . From (6) the σ-
control u providing for the compensated dynamics is identified as (we omit the subscript)
σ − dynamics in (8). The result is given in the σD = −Vr Vλ / r + AT ,λ − c0Vr /(2 r ) − cos(λ − γ M )nL (15)
following Theorem. so the commanded acceleration, n Lc , for the
Theorem. Let r = 3, α1, α2 > 0, and | DgDD | ≤ L. Then interceptor normal acceleration n L is selected to be
the closed-loop system (7) controlled by the smooth 1 2/ 3
control function nLc = (α1 σ signσ +
(r −1) / r (r −2) / r
cos(λ − γ M )
u = z1 − α1 σ ∫
sign(σ ) − α0 σ sign(σ )dτ (14)
VV c V
(16)
∫ sign(σ )dτ − r λ − 0 r + Aˆsmooth T ,λ ,
1/ 3
α0 σ
where z1 is the output of the disturbance observer r 2 r
(11) is finite-time stable and
• If σ (t ) and u(t) are measured with some where Aˆ smoothT ,λ is estimated using (12).
Lebesgue-measurable noises bounded by ε > 0 , The homing (when a seeker is activated and the
designed guidance law is executed) close loop
then inequality σ (t ) ≤ γε is established in interceptor kinematics (3), (6), (14) is derived as
finite time for some constant γ > 0 depending follows:
exclusively on the parameters of the controller. rD = Vr ,
• Let τ > 0 be the sampling interval and the noises D
Vr = c0 + AT r − sin(λ − γ M )nL ,
2
be absent, then the inequality σ(t ) ≤ µτ3 is D (17)
established in finite time for some constant µ > 0 λ = c0 / r ,

depending exclusively on the parameters of the Vλ = c0 r ,
controller. Since the initial condition Vr (0) = M << 0 , it is easy
The proof follows from the overall homogeneity of to identify c 0 such that Vr (t ) < 0 ∀t ≤ t go , where t go
the closed system (7), (12), (14) and is similar to the is time required for the interceptor to hit the target.
proof of Lemma 2. Let w0=z0-σ, w1=z1-g, w1=z2- gD . So, the guidance law (15) based on asymptotic
The following transformation is used: smooth second order sliding mode control inevitably
Gκ: (t,τ,ε,x1,x2w0,w1,w2) guarantees hit-to-kill intercept.
3 3 2 3 2
v ( κt, κτ,κ ε,κ x1,κ x2,κ w0,κ w1,κw2). 3. INTEGRATION OF GUIDANCE AND FLIGHT
CONTROL SYSTEMS: 2-LOOP SECOND ORDER
2.6 Smooth Second Order Sliding Mode Guidance SMC APPROACH
At this moment, the chosen intercept strategy 3.1 Missile Interceptor Dynamics
transforms the intercept problem to an output
The following pitch plane model of kinetic energy
regulation problem, which can be considered as a
kill vehicle (KEKV) without account for gravity (see
fundamental problem in geometric methods of
Shkolnikov, Shtessel, et. al., 2000) is considered
control synthesis. The control goal of the output
assuming the model to be affine in controls, i.e. in
regulation problem is to stabilize σ i , i = 1,2 and to
case of only aerodynamic control ( δ e is a virtual
satisfy the condition ∀t ≥ t o , Vr < 0 provided given
deflection, which is allocated by the fin’s mixing logic)
limits on control input, n L ≤ n L max , and known we have
bounds of uncertainties variations. We wish to solve  X Z 
− mV sin α + mV cos α + q   0 
this problem without explicit estimation of α     
uncertainties but to apply robust control methods,    X Z   0 
where control action will automatically counteract V  =  m cos α + m sin α + 0  ⋅ δ e ,(18)
θ     
the influence, which steers the system away from    q   q Sl δ e 
q   q Sl   I C m (⋅) 
specified path. A number of solutions to this problem
using SMC theory is developed in the last decade.
 ( C mo (⋅) + ∆C m (t ) )
  yy 
 I yy 
Each of them can be used to create a reference
qS
closed-loop model to be enforced in presence of the nL = C L (α ,...) = − Ax sin α + Az cos α = γ M ⋅ V (19)
m
unmodeled at this stage vehicle dynamics of a
homing interceptor. The important issue here is to where γ M = θ − α is the flight path angle, n L is a
design a continuous/smooth control, n Lc (t ) , to be X Z
normal acceleration, Ax = , Az = , q = 12 ρV 2 ,
enforced for the actual missile normal acceleration m m
n L (t ) . That’s why absolute robustness to X = q S (C x (⋅) + ∆C x (t ) ) , Z = q S (C z (⋅) + ∆C z (t ) ) ,
uncertainties and finite time convergence are traded C m (α , αD , Mach, q, δ e ,...) = C mo (⋅) + C m
δe
(⋅) ⋅ δ e , and
for approximation of an ideal discontinuous SMC
∆C x (t ), ∆C z (t ), ∆C m (t ) are the external disturbances.
law by a linear saturation function (see Moon, &
Kim, 2000; Shkolnikov, Shtessel, et. al., 2000). In It’s also assumed that the model (18) is of minimum
this case the output is robustly regulated to a small phase. In case of aerodynamic control, it means that
domain of convergence around zero. ∂C L ∂C D
= 0, = 0 , which is true, if the extra small
∂δ e ∂δ e
fins are employed. The actuators have the highly Although, the inner loop second order sliding mode
nonlinear dynamics with the relative degree equal to control can be designed to track nD L , it is more
one, which, after mixing control commands to convenient to track qc , a command for the pitch rate,
individual actuators actuator, can be presented as
since the pitch rate is measuring accessible. From the
δDe = − a(δ e , t )(δ e − u ) , (20) interceptor dynamics (18) we identify
where the nonlinear uncertain term a (δ e , t ) is allowed 1 
to be non-smooth but bounded.
q=
V ∫
 Tα nD L + nD L dτ  ,

(26)
where V is an interceptor speed (magnitude), and
3.2 Model Behavior of the Homing Loop
Tα is known as the turning rate coefficient
Now, the problem is to design the control input to the −1
α  m  ∂C L 
actuator, u (t ) , in order to enforce the given closed Tα = = V   (⋅)  (27)
γD M  q S  ∂α 
loop performance of a homing loop robustly to
uncertainty of vehicle dynamics. Under the second Considering Tα as a known slowly varying quantity
order smooth asymptotic SMC guidance (15), (16) (that is true for the majority of intercept scenarios),
the closed-loop Vλ -dynamics are we finally obtain the following profile for the
interceptor pitch rate to follow
VDλ = c 0V r /(2 r ) (21)
1 2/3
It was shown in the work of Shkolnikov, Shtessel, et. qc (t ) = [Tα (α1 σ 0 sign(σ 0 ) +
V cos(λ − γ M )
al., 2000) by simulations that even in presence of
∫
1/ 3
uncompensated AT ,λ the closed-loop Vλ -dynamics α0 σ 0 sign(σ 0 )dτ + ϕˆsmooth)dτ +
enforced by the traditional SMC guidance can
∫ (α σ ∫
2/3 1/ 3
provide for the reliable target intercept. The model 1 0 signσ 0 + α0 σ 0 sign(σ 0 )dτ )dτ ]
behavior (21) should be robustly enforced by control Implementing the pitch rate command it necessary to
u in presence of vehicle dynamics uncertainties and
disturbances and target maneuvers represented by compute VD , since σ = VD − c V / 2 r . It can be
λ 0 λ 0 r ( )
target acceleration components AT ,λ and AT ,r . This done using the first order exact differentiator (see
Levant 1998, 2003), or the first order robust-to noise
problem is addressed via SOSM approach (see
differentiator based on nonlinear dynamic sliding
Levant 1993, 1998, 2001, 2003; Fridman 2001,
manifold (NDSM) (see Krupp, Shkolnikov, &
2003) and a two-loop controller structure (see
Shtessel, 2000), where NDSM is developed to
Shtessel, Buffington, and Banda, 1999).
incorporate a low-pass nonlinear filter in the observer
3.3 Outer Loop SMC Design structure producing a robust to noise differentiator. A
robust differentiator of a given signal, x (t ) , polluted
Now we consider the composite state space of the
systems (3) and (18). If we hold with noise, has the form
σ 0 = VDλ − c 0V / 2 r ( ) (22)  DxˆD = ρo sign( J o ),
 J o = χ + e,
at zero, then even under small residual perturbation,  χD = b e − a ( χ + e) , (28)
 e = x(t ) − xˆ,
| σ 0 (t ) |≤ N , an acceptable closed-loop performance e
0.5
χ +e
0.5

that leads to the hit-to-kill intercept can be achieved
where the parameters should be appropriately
(see Shkolnikov, Shtessel, et. al., 2000).
selected in accordance with the upper boundary for
Following the SMC approach, we call the constraint DxD(t ) ≤ L . Apparently, when the 2-sliding mode,
(22) to be kept at zero as the sliding manifold in the
outer loop, where σ 0 is the sliding quantity. To J o = JD o = 0 , is established, we have the estimation
stabilize σ 0 to zero, its dynamics is identified 0.5
error dynamics eD = −b e sign(e) , which converges

σD0 = −
(VDrVλ + VrVDλ )r −Vr2Vλ + AD to zero in a finite time. Thus, a command on missile
T ,λ −
r2 maneuver is obtained in terms of pitch rate command
, (23)
 Vλ nL  q c . To obtain this command we had to consider the
sin(λ − γ M ) − nL − cos(λ − γ M )nDL composite state space of engagement kinematics and
r V
vehicle dynamics. The inner loop control is designed
One can rearrange the terms in (23) and write it in
next to robustly enforce q c in presence of
the short cut notation as
σD 0 = ϕ − cos(λ − γ M )nD L , (24) uncertainties and disturbances.
where 3.4 Inner Loop Second Order SMC Design
ϕ = ϕ (VDr ,Vr ,VDλ ,Vλ , r ,V , ADT ,λ , λ , γ M , nL ) = The regulated output in the inner loop is the pitch
(VD V + VrVDλ )r − Vr2Vλ + AD
− r λ
 Vλ nL  tracking error
T ,λ − sin(λ − γ M ) −  nL e q = q c (t ) − q . (29)
r2  r V 
Designing a smooth second order SMC (13) in terms From (18) and (20) we determine that the relative
of nD Lc (“jerk” command), we obtain the following degree of the input-output dynamics for e q

eDDq = ϕ1 (⋅) − α (Tα )−1ω n2 (⋅)a (.)u ,

command on missile normal jerk nD L c T
(30)
nDLc =
1
cos(λ − γ M )
( 2/ 3 1/ 3
∫
α1 σ0 signσ0 +α0 σ0 sign(σ0)dτ , (25) ) V
is equal to two, where ϕ1 (.) is a nonlinear, time- IEEE Transaction on Aerospace and Electronic Systems,
No.2, 306-325.
varying uncertain bounded term. Thus, one can apply
Floquet, T., Barbot, J.-P., & Perruquetti, W. (2003).
the SOSM approach (see Levant 1993, 1998, 2003) Higher-order sliding mode stabilization for a class of
to provide for finite time convergence of the tracking nonholonomic perturbed systems. Automatica, 39, 1077 –
error e q to zero robustly to time-varying additive 1083.
uncertainty ϕ1 (⋅) and multiplicative uncertainties of Fridman, L., (2001). An averaging approach to chattering.
IEEE Transactions of Automatic Control, 46, 1260 –
q Sl δ e 1265.
airframe dynamics, ω n2 = ω n2 (⋅) = Cm (⋅) ( ω n is
I yy Fridman, L. (2003). Chattering analysis in sliding mode
systems with inertial sensors, International Journal of
an equivalent undamped natural frequency of Control, 76 (9/10), 906-912.
airframe) and actuator dynamics a = a (δ e , t ) Garnell, P., and East, D. J. (1977). Guided Weapon Control
(actuator bandwidth). The only requirement is to Systems, Oxford: Pergamon Press.
know the limit of ϕ1 (⋅) variations, and the nominal Krupp D. R., Shkolnikov I. A., and Shtessel Y. B., (2000).
2-Sliding Mode Control for Nonlinear Plants with
1 2
value and sign of the term ω n (⋅)a (⋅) . For instance, Parametric and Dynamic Uncertainties. Proceedings of
V AIAA Guidance, Navigation and Control Conference,
SOSM that provides convergence of eq to zero in Denver, CO, Paper AIAA-2000-3968.
Levant A. ( Levantovsky L.V.) (1993). Sliding order and
finite time in presence of bounded uncertainties can sliding accuracy in sliding mode control. International
be designed in a format (see Levant 1993, 2001) Journal of Control, Vol. 58, No. 6, 1247-1263,
V 1/ 2 Levant, A. (1998). Robust exact differentiation via sliding
u = 2 α sign(eDq + eq sign(eq )) (31) mode technique. Automatica, (34), 3, 379-384.
ωn Levant A., (2001). Universal SISO sliding-mode
In order to avoid differentiation of e q control (31) controllers with finite-time convergence. IEEE
Transactions on Automatic Control, 46(9), 1447-1451.
can be also designed based on NDSM (see Krupp, Levant A., (2003). Higher-order sliding modes,
Shkolnikov, & Shtessel, 2000) as follows: differentiation and output-feedback control. International
u = 2 ρ sign( J ), J = e q + χ
V Journal of Control, 76 (9/10).
ωn a Moon, J., Kim, Y., (2000). Design of missile guidance law
(32) via variable structure control. Proceedings of AIAA
0.5 0.5
χD = β 1 e q sign(e q ) − β 2 J sign( J ) Guidance, Navigation and Control Conference, Denver,
CO, Paper AIAA-2000-4068.
where J is the nonlinear dynamic sliding quantity, Orlov, Y., Aguilar, L., and Cadiou, J. C. (2003). Switched
the nominal values for uncertainties ω n2 (⋅) and a (⋅) chattering control vs. backlash/friction phenomena in
electrical servo-motors. International Journal of Control,
are used, and the coefficient ρ is selected according 76 (9/10), 959 - 967.
to the upper absolute limit for additive uncertainty Shkolnikov, I. A., Shtessel, Y. B., Lianos, D. P., and Thies,
ϕ1 (⋅) . Thus, using only output feedback e q , the A. T., (2000), Robust Missile Autopilot Design via High-
Order Sliding Mode Control. Proceedings of AIAA
control voltage (32) to the actuator (20) provides for Guidance, Navigation, and Control Conference, AIAA
the output e q convergence to zero in a finite time. Paper 2000-3968.
4. SIMULATION EXAMPLE Y. Shtessel, I. Shkolnikov and M. Brown, (2003). An
Asymptotic Second-Order Smooth Sliding Mode
Simulation results of KEKV using smooth SOSM Control,” Asian Journal of Control, 4, No. 5, 498-504.
guidance and SOSM/NDSM-based autopilot are Shtessel, Y., Buffington, J., and Banda, S., (1999).
available and are omitted for brevity. Multiple Time Scale Flight Control Using Re-
configurable Sliding Modes, AIAA Journal on Guidance,
5. CONCLUSIONS Control, and Dynamics, (22), 6, 873-883.
An integrated two-loop guidance and flight control Utkin, V., Guldner, J., and Shi, J. (1999). Sliding Modes in
system is designed based on second order sliding Electromechanical Systems, Taylor and Francis, London.
mode control to robustly enforce hit-to-kill guidance Zarchan P., (1998). Tactical and Strategic Missile
strategy in presence of target maneuvers, Guidance. Vol. 176, Progress in Astronautics and
Aeronautics, AIAA Publications.
atmospheric disturbances, and dynamic uncertainty
Zhou, D., Mu, C. & Xu, W. (1999). Adaptive Sliding-
of airframe-actuator. The zero miss target intercept is Mode Guidance of a Homing Missile. Journal of
achieved with the reduced acceleration ratio Guidance, Control and Dynamics, (22), 4, 589-594.
requirements due to prolonged stability of the
homing loop.
REFERENCES
Bacciotti A. & Rosier L., (2001). Liapunov Functions and
Stability in Control Theory, Lecture notes in control and
information sciences, 267, Springer-Verlag, New-York.
Babu, K., Sarma, I., and Swamy, K., (1994). Switched Bias
Proportional Navigation for Homing Guidance Against
Highly Maneuvering Target. Journal of Guidance,
Control and Dynamics, 17(6), 1357-1363.
Brierly, S.D., & R. Longchamp, (1990). Application of
Sliding Mode Control to Air-Air Interception Problem.