This paper deals with missile autopilot design for agile turn phase in airtoair engagement scenarios. To attain a fast response,angleofattack (AOA) is adopted for an autopilot command structure. Since a high operational AOA is generally required during the agile turn phase, dealing with the aerodynamic uncertainties can be a challenge for autopilot design. As a remedy, a new controller design method based on robust nonlinear control methodology such as time delay control is proposed in this paper.Nonlinear observer is also proposed to estimate the AOA in the presence of the model uncertainties. The performance of the proposed controller with variation of the aerodynamic coefficients is investigated through numerical simulations.
1. Introduction
In recent years, the maneuverability and survivability of aircraft against airtoair missiles during closedincombat scenarios has grown. For airtoair missile systems, an agile turn capability is recommended to effectively shoot down an enemy's aircraft. To improve the missile's maneuverability and achieve agile turn capability, a so called high angleofattack (AOA) maneuver is required. When the missile is operating at high AOA regions, aerodynamic nonlinearity and uncertainty can be a problem. Accordingly, an autopilot with robustness and a fast response is needed. As a remedy to this,some previous works have been carried out as follows.
Agile missile autopilots have been designed based on the classical control approach with gainscheduling method(Ryu et al., 2010; Wise and Broy, 1998). In Mehrabian’s work,(Mehrabian and Roshanian, 2006) a linear parameter varying scheme has been applied to the missile autopilot design. A pitch controller with aerodynamic fin and side thrust has also been designed by using dynamic inversion and the extendedmean assignment by Choi et al. (2003). As shown by Menon and Yousefpor (1996), the authors proposed the autopilot design method by combining the feedback linearization technique and linear control theory. The sliding mode control has been also applied to the autopilot problem for agile missile design by Thukral and Innocenti (1998). In the work of Kim et al. (2010), adaptive backstepping methodology has been used for pitch autopilot design. Theorists have also devised sliding mode controllers and time delay controllers in longitudinal missile motion to facilitate agile turn (Kim et al., 2010; Lee et al., 2009c).
In this paper, the authors build on previous works as discussed above (Lee et al., 2009b, c). The AOA is adopted for an appropriate autopilot command structure to achieve a fast response. Then, the robust pitch autopilot is designed based on the time delay control (TDC) methodology (YoucefToumi and Wu, 1991). The advantage of the TDC approach is that if an input gain is partially known, a controller can be designed even though the plant dynamics are unknown. For this reason, the TDC can be applicable in the case of agile turn under the presence of model uncertainties. To design an autopilot, we should derive the approximated aerodynamic model, which is given by a polynomial function of AOA, from the experimental aerodynamic data. Next, the inputoutput linearization technique is used to attain an appropriate form of dynamics for the AOA controller. By using this equation that represents an explicit relationship between tracking state and input, we apply the TDC scheme to pitch autopilot design. In order to implement the proposed AOA controller,feedback of AOA information is needed. Unlike normal acceleration, it is hard to measure the AOA in flight. Instead,since the AOA is directly related to the normal acceleration and the pitch rate, we construct a nonlinear observer in order to estimate AOA that uses the normal acceleration and the pitch rate measurement. Finally, we develop the autopilot design scheme by combining the robust AOA controller and the nonlinear observer.
This article is composed of four sections. In Section 2, the missile model used in the autopilot design is discussed. The proposed autopilot design method is provided in Section 3.To show the validity of the proposed controller, numerical simulations are performed and its results are given in Section 4. Finally, conclusions of this study are offered in Section 5.
2. Missile Model
In this section, the nonlinear missile model considered hereafter is explained. The missile motion can be decoupled into two perpendicular channels in the skidtoturn cruciform type missile. A longitudinal motion with a high AOA capability is considered. Under the assumption that the missile has rigid body dynamics and the gravity force is ignored, the missile dynamics including the firstorder actuator dynamics can be expressed as:
The approximation of C_{N0}(left) and C_{Nδ}(right).
where the notation of α, q and δ describe AOA, pitch rate and fin deflection angle, respectively. The missile velocity,mass, moment of inertia with respect to pitch axis, dynamic pressure, reference area, and reference length are expressed by V, m, Iyy, Q, S and l, respectively. τ represents time constant of the actuator and M is the Mach number.
In Eqs. (1) and (2), the missile dynamics include the aerodynamic coefficients, CN0, C
_{Nδ}
, C
_{m0}
, C
_{mδ}
, C
_{mq}
, which are dependent upon the variation of AOA, Mach number, and fin deflection. In general, these parameters are not directly given by the function of AOA, Mach number, and fin deflection.Instead, these are computed by interpolation of experimental data tables. For controller design purposes, the aerodynamic coefficients must be given by a function form. As shown in
Figs. 1
and
2
, for a fixed Mach number, these parameters can be represented as 3rd and 4th order polynomial functions of AOA. Accordingly, under Mach number
M
holds on constant, we can approximate aerodynamic parameters as the polynomial function forms as shown by Choi et al.(2003).
where
a_{i}
,
b_{i}
,
c_{i}
,
d_{i}
denote the coefficients of each polynomial function, which are depending on Mach number. These parameters can be computed using a lookup table at each Mach number. It is noted that fin deflection and AOA are coupled to each other. Substituting Eqs.(4) (7) into Eqs.(1) and (2), we have equations as follows:
where
K_{α}
,
K_{q}
and
C_{q}
can be regarded as the AOA gain,the pitch rate gain, and the aerodynamic pitch rate gain,respectively. These variables are defined as:
3. Autopilot Design
In this section, we discuss the missile autopilot design based on the TDC with nonlinear observer (YoucefToumi and Wu, 1991). The proposed controller configuration is shown in Fig. 3 In order to design the AOA controller, we first derive the dynamics which represent an explicit relationship between
y
and
u
by using the inputoutput linearization technique. In the autopilot design step, we take aerodynamic uncertainties into consideration.
Fig. 3 The configuration of the proposed control method.
 3.1 Inputoutput linearization
In order to derive an appropriate equation for applying the TDC, the inputoutput linearization technique (Slotine and Li, 1991) is first applied to the above missile model as given in Eqs.(1) through (3). For the control objective, we choose AOA as the output state. The output is not directly related to control input in the missile model so that we take a timederivative of output to find the direct relationship between the system output y and control input
u
:
The approximation of C_{m0}(left) and C_{mδ}(right).
At the second derivative of AOA, control input
u
is appeared. In Eq. (13), the timederivative of aerodynamic coefficients are given as:
For simplicity's sake, let the new functions be defined as:
Substituting Eqs.(14) through (17) into Eq.(13), the second timederivative of AOA can be rewritten as:
Equation (18) represents an explicit relationship between
y
and
u
. In this equation, the relative degree is two. This means that there is one internal dynamic in this system.To design the AOA controller by using Eq. (18), the internal dynamics must be stable. In advance of applying the TDC methodology, we need to check the stability of the internal dynamics. It can be proven by the following theorem:
Theorem 1. When angle of attack is chosen for output state in inputoutput linearization of the system as given in Eq.(1)through (3), the internal dynamics is stable.
Proof . From Eq. (18), the control input can be formulated as:
To obtain the internal dynamics, we substitute the obtained control input as shown in Eq.(19) into Eq.(3). Then,we get the internal dynamic as below:
It is noted that the internal dynamic is considerably complicated to check its stability. Subsequently, we indirectly consider the zero dynamic instead of the internal dynamics.The zero dynamics is given by the internal dynamics when the output is kept at zero by the control input as α = α?= α¨ =0.Applying this condition to the internal dynamics, we have the following zero dynamics as:
where C
_{1}
and C
_{2}
are defined as:
According to approximated aerodynamic coefficients, C
_{1}
and C
_{2}
are given by positive values. In addition, d
_{2}
is positive and K
_{α}
is negative. Consequently, C
_{3}
is a positive value. From Eq.(21), it can be easily shown that the zero dynamics is asymptotically stable.
 3.2 Time delay controller design
The dynamics of AOA as shown in Eq. (18) can be rewritten in a simplified form,
where
f
and
b
are determined as:
To design the AOA controller based on the TDC, Eq. (23)is reformulated as:
In this equation, it is assumed that the dynamics
f
and the control input gain
b
are not exactly known because of aerodynamic uncertainty. Instead, it is assumed that only the input gain □ can be estimated from the partially known model information. Then, the right side of Eq. (26) can be represented by using a partially known term and an unknown term:
where the unknown term
h
can be expressed as:
Next, we construct the control input
u
as following equation:
where
v
represents the pseudocontrol input and is given by:
where e
_{d}
=
^{Δ}
α
_{c}
α, which denotes tracking error between command and response. Substituting Eq. (29) and (30) into Eq. (27), the error dynamics can be obtained as:
where
K_{p}
and
K_{v}
denote the gains of the error dynamics that can be chosen to satisfy the desired closedloop performance. In order to implement the control law as shown in Eq. (29), the unknown term
h
should be estimated.In TDC methodology, the following estimation method is used to determine
h
. If a time delay
L
is sufficiently small, we can use the following approximation as:
Substituting Eq. (32) into Eq. (29), we finally have the TDC law as follows:
According to Eq. (33), it is noted that the control law does not require the plant model. However, to estimate an unknown term, the previous control input and the second derivative of AOA measurement are additionally needed.
 3.3 Nonlinear observer design
In order to implement the proposed controller as given in Eq. (33), the AOA measurements are required, but it is difficult to obtain those measurements from conventional devices. Accordingly, in this paper, the authors propose a nonlinear observer to estimate AOA as studied suggested by Devaud et al. (1998). In their work, the nonlinear observer was designed based on the known plant model, while we assume that the model of the nonlinear observer also has uncertainty. If the model of the nonlinear observer is exactly known, the normal acceleration feedback is enough to estimate AOA information. Conversely, we use the normal acceleration feedback and the pitch rate to improve the performance of the nonlinear observer against the model uncertainty as follows:
where â
_{n}
and qˆ denote the estimation of normal acceleration and the estimation of pitch rate, respectively. ân and q? can be determined as:
where
K_{n}
=
^{Δ}
Qs / m
. Let the AOA estimation error be defined as
e_{d}
Δ α
_{c}
 α, then we choose the Lyapunov candidate function as follows:
Taking the timederivative of
V_{2}
, we have the following equation:
where:
To make V˙
_{2}
< 0 and to ensure the exponential stability of the nonlinear observer, for a positive real number ε and some real number
K_{1}
and
K_{2}
are chosen as follows:
In order to select proper values of
K_{1}
and
K_{2}
, a trial and error procedure may be required. In practice,
K_{1}
and
K_{2}
are designed in such a way as follows:
From Eq. (43), the first term of the lefthand side is independent of
K_{1}
and
K_{2}
. Therefore, Eq. (43) can be rewritten as follows:
In Eq. (44), it can be easily observed that the value of the lefthand side has two degrees of freedom related to
K_{1}
and
K_{2}
. Therefore, some proper combination of
K_{1}
and
K_{2}
can satisfy the condition of Eq. (44). A detailed explanation of the closedstability issue including these types of estimators can be found the work of Devaud et al. (1998).
4. Numerical Analysis
In this section, to investigate the performance of the proposed controller, numerical simulations are carried out in conjunction with the nonlinear observer. The missile model used in simulation is given by the following parameters. The mass and pitch axis moment of inertia are
m
= 56.3
kg
and
I_{yy}
= 44.36
kg

m^{2}
, respectively. We assume that the missile is launched with Mach. 0.85. The reference area and length are
S
= 0.01824
m^{2}
and
l
= 0.1524
m
. The natural frequency of the firstorder actuator model is 200
rad
/
s
. The missile's aerodynamic model parameters are given by
Table. 1
.The time delay controller are designed by the following parameters:
Aerodynamic model parameters
Aerodynamic model parameters
We assume that initial estimation of angleof attack is zero, α? (0) = 0˚. The parameters for AOA estimator are chosen as:
The remainder of this section is concerned with three kinds of nonlinear simulation results. In a high operational AOA range as the agile missile, the aerodynamic parameters could not be obtained accurately. Therefore, to test the robustness of the proposed controller, we consider 30%aerodynamic uncertainties in all simulations. In order to implement the proposed controller, partial information of the plant model is required so that the proposed controller may have robustness against model uncertainties.
In the first place, the proposed controller with known and unknown AOA measurement is applied to the longitudinal missile model. In the known AOA measurement case, we assume that AOA can be directly measured. In the unknown AOA measurement case, it is assumed that the AOA is estimated from the plant model which contains model uncertainties; it is the same as the estimation of the AOA using the nonlinear observer as shown in Eq. (34) with
K_{1}
=0 and
K_{2}
= 0. Second, the proposed controller conjunction with the proposed nonlinear observer is performed. In this simulation, we compare the estimation result and controller performance between the established observer (use of
K_{1}
only) and the proposed observer (use of
K_{1}
and
K_{2}
). Finally,the proposed controller is applied to an agile turn engagement scenario. In this case, the AOA command is obtained from the optimization results as studied by previous theorists(Kim et al., 2009, 2010; Lee et al., 2009a).
Angleofattack and control input with known and unknown measurement.
Angleofattack response and estimation error with the nonlinear observer.
Angleofattack and control input in 180° turn scenario.
Figure 4
provides the step input response of AOA and control input in cases of known and unknown AOA measurement. If we use the true AOA measurement, the proposed controller has a good tracking performance in spite of the absence of the plant model information and the uncertainty of input gain. When the unknown AOA measurement is applied to the proposed controller, the tracking error increases as shown in
Fig. 4
For this reason, we need the nonlinear observer to compensate the tracking error under the absence of the true AOA measurement.
Figure 5
provides the AOA response and the estimation error under the proposed controller with the nonlinear observer. Two simulation results as shown in
Fig. 5
provide a good tracking performance in the presence of model uncertainties. However, in the nonlinear observer as studied in Devaud et al. (1998), which uses
K_{1}
only, the estimation error increases, while the proposed nonlinear observer can compensate for this divergence.
Figure 6
provides the results of the controller in case of 180 degree heading reversal. In such an agile turn scenario, the proposed controller shows excellent performances under the aerodynamic uncertainty.
5. Conclusions
In this paper, based on the TDC methodology, we discuss the pitch autopilot design for missiles operating at high AOA regions in the presence of aerodynamic nonlinearity and uncertainty. The authors first derive the appropriate nonlinear dynamics for AOA control based on inputoutput linearization techniques and prove its internal dynamics stability. From the results, the design method of the robust AOA controller against model uncertainty using the time delay estimation techniques is provided. In order to implement the proposed controllers, a nonlinear observer is developed to estimate AOA under model uncertainty.Numerical simulations show that the proposed method can provide a good tracking performance against model uncertainty and can be successfully applied to the agile turn phase of airtoair missiles.
Acknowledgements
This research was supported by Agency for Defense Development in Korea under the contract ADD09010303.
Choi Y. S
,
Lee H. C
,
Choi J. W
2003
Autopilot design for agile missile with aerodynamic fin and side thruster
SICE 2003 Annual Conference
Fukui Japan
1476 
1481
Devaud E
,
Siguerdidjane H
,
Font S
Nonlinear dynamic autopilot design for the non minimum phase missile
Proceedings of the IEEE Conference on Decision and Control
Tampa FL
1998
4691 
4696
Kim K. U
,
Kang S. H
,
Kim H. J
,
Lee C. H
,
Tahk M. J
Realtime agileturn guidance and control for an airtoair missile
AIAA Guidance Navigation and Control Conference
Toronto Canada
2010
29 
31
Kim Y. H
,
Park J. W
,
Lee C. H
,
Tahk M. J
Control input discretization and optimization method for halfturn maneuver analysis of short range airtoair missile
Proceeding of KSAS Spring Conference
Pyeongchang Korea
2009
373 
376
Lee C. H
,
Kim T. H
,
Ryu S. M
,
Tahk M. J
Analysis of turning maneuver for agile missile using trajectory optimization
Proceeding of KSAS Spring Conference
Pyeongchang Korea
2009
369 
372
Lee C. H
,
Kim T. H
,
Tahk M. J
Autopilot design for agile turn using timedelay control with observer
Proceeding of KSAS Fall Conference
Gyeongju Korea
2009
822 
825
Lee C. H
,
Kim T. H
,
Tahk M. J
Design of missile autopilot for agile turn using nonlinear control
Proceeding of KSAS Fall Conference
Gyeongju Korea
2009
667 
670
Mehrabian A. R
,
Roshanian J
Design of gainscheduled autopilot for a highlyagile missile
Proceedings of the 1st International Symposium on Systems and Control in Aerospace and Astronautics
Harbin China
2006
144 
149
DOI : 10.1109/ISSCAA.2006.1627599
Menon P. K
,
Yousefpor M
Design of nonlinear autopilots for high angle of attack missiles
AIAA Guidance Navigation and Control Conference
San Diego CA.
1996
29 
31
Ryu S. M
,
Won D. Y
,
Lee C. H
,
Tahk M. J
High angle of attack missile autopilot design by pole placement approach
Proceedings of the 3rd International Symposium on Systems and Control in Aeronautics and Astronautics
Harbin China
2010
535 
539
DOI : 10.1109/ISSCAA.2010.5632233
Slotine J. J. E
,
Li W
1991
Applied Nonlinear Control
Prentice Hall.
Englewood Cliffs
Thukral A
,
Innocenti M
1998
A sliding mode missile pitch autopilot synthesis for high angle of attack maneuvering
IEEE Transactions on Control Systems Technology
6
359 
371
DOI : 10.1109/87.668037
Wise K. A
,
Broy D. J
1998
Agile missile dynamics and control.
Journal of Guidance Control and Dynamics
21
441 
449
DOI : 10.2514/2.4256
YoucefToumi K
,
Wu S. T
Input/output linearization using time delay control
Proceedings of the American Control Conference
Boston MA
1991
2601 
2606