Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 412509, 13 pages

http://dx.doi.org/10.1155/2013/412509

Research Article

Output Tracking Control of Switched Hybrid Systems: A Fliess Functional Expansion Approach

Fenghua He,¹ Peng Zhang,¹ Yong Chen,² Long Wang,¹ Yu Yao,³ and Weishan Chen⁴

¹Control and Simulation Center, Harbin Institute of Technology, Harbin, Heilongjiang 150080, China

²Department of Mathematics, Harbin Institute of Technology, Harbin, Heilongjiang 150080, China

³Harbin Institute of Technology, Harbin, Heilongjiang 150080, China

⁴Department of Mechanical Engineering, Harbin Institute of Technology, Harbin, Heilongjiang 150080, China

Received 5 April 2013; Accepted 18 May 2013

Academic Editor: Rongni Yang

Copyright © 2013 Fenghua He et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The output tracking problem is investigated for a nonlinear affine system with multiple modes of continuous control inputs. We convert the family of nonlinear affine systems under consideration into a switched hybrid system by introducing a multiple-valued logic variable. The Fliess functional expansion is adopted to express the input and output relationship of the switched hybrid system. The optimal switching control is determined for a multiple-step output tracking performance index. The proposed approach is applied to a multitarget tracking problem for a flight vehicle aiming for one real target with several decoys flying around it in the terminal guidance course. These decoys appear as apparent targets and have to be distinguished with the approaching of the flight vehicle. The guidance problem of one flight vehicle versus multiple apparent targets should be considered if no large miss distance might be caused due to the limitation of the flight vehicle maneuverability. The target orientation at each time interval is determined. Simulation results show the effectiveness of the proposed method.

1. Introduction

In the last decade, hybrid systems have been received considerable interests from researchers in control and computer science fields; see [1–7]. One particular class of hybrid systems, switched systems, has been extensively investigated during the past several years; refer to [8–17] and the references therein. A switched system can be obtained from a hybrid system by ignoring the discrete behavior and constrain all possible switching problems form a certain class. A switched system consists of several subsystems and a switching law specifying the subsystem at each time instant. We call the switching triggered by an external event externally forced switching and the switching triggered by an internal event internally forced switching. In this work, the latter type of switching will be addressed.

For switched systems, the optimal switching control problems have attracted much attention in recent years. [18] provided a brief survey on recent results of switched systems optimal control. The available theoretical results were usually extending the classical maximum principle or the dynamic programming approach. The optimal control of discrete-time switched systems can be modelled as piecewise affine systems by partitioning the state space into polyhedral regions. The optimal control for constrained discrete-time linear hybrid systems was provided based on quadratic or linear performance criteria, and the construction of the state-feedback optimal control law was also treated by combining multiparametric programming and dynamic programming in [19]. A recent result for the optimal control of discrete-time switched linear systems is [20]; the optimal control input was given in piecewise linear state-feedback form by solving a set of difference Riccati equations and the switching sequence and time were obtained by dynamic programming. For continuous-time switched systems, a two-stage optimization method was proposed and proved when a sequence of active subsystem was prespecified; see [21–23]. In [24], a switching system under consideration was embedded into a larger family of systems and the optimization problem was then formulated, and necessary and sufficient conditions for optimality were also provided. In addition, a numerical technique was developed in [25] for solving hybrid optimal control problems proposed in [24], and the methodology presented can be extended to incorporate hybrid behavior stemming from autonomous (uncontrolled) switches that result in plant equations with piecewise smooth vector fields. The problem of optimal switching for switched linear systems was studied in [26], where the switching actions were described as multiple control inputs; Pontryagin Maximum Principle was applied for solving the optimal control problem by embedding the switched system in a larger family of systems.

Motivated by the application example that we will treat in this paper, we consider nonlinear affine systems with continuous control inputs. By introducing a multiple-valued logic variable, the systems under consideration can be converted into a switched hybrid system with candidate modes of control inputs. The output tracking problem of such a class of switched hybrid system is investigated by applying the Fliess functional expansion approach. The Fliess functional expansion was studied via numerical methods for nonlinear control system in [27, 28]. In this work, we consider the control modes of inputs are finite and constant in a sampling time duration, and then the output-input relation can be approximately expressed via the Fliess function expansion. The optimal switching control can be obtained via dynamic programming or simple brutal method if the control input modes are not too much. We apply the proposed approach in a multitarget tracking problem, which shows the effectiveness of the proposed approach.

The rest of the paper is organized as follows. Section 2 provides preliminaries of Fliess function expansion of a nonlinear affine system. A nonlinear affine system with multiple modes of continuous control inputs is converted into a switched hybrid system, and then the output tracking problem for such a class of switched hybrid system is formulated. In Section 3, the optimal switching control is obtained for one-step prediction horizon and multiple-step prediction horizon case, respectively. In Section 4, the proposed approach is applied to a multitarget tracking problem for a scenario in which one flight vehicle is aiming for a real target with several decoys flying around it, and simulation experiments are also done to illustrate the effectiveness of the proposed approach. Finally, some concluding remarks and future directions are provided in Section 5.

Before the end of this section, some remarks on notations in guidance control field are given as follows. We denote by the LOS angle between a flight vehicle and a target. is the relative distance. and are the relative velocity along the LOS and normal to the LOS, respectively. and are the acceleration of the flight vehicle along the LOS and normal to the LOS, respectively. is the navigation constant.

2. Preliminaries and Problem Formulation

In this section, we will provide preliminaries of Fliess functional expansion for a nonlinear affine system. We convert a nonlinear affine system with multiple modes of continuous control inputs into a hybrid switched system. Then the output tracking problem of such a class of switched hybrid system is formulated.

2.1. Preliminaries

To begin with, we revisit Fliess function expansion of a nonlinear affine system. Consider a nonlinear affine system described as follows:

where is the state vector, is the control input, is the output vector, and are smooth vector fields, and , , , are smooth functions.

The relative degree vector for the above nonlinear affine system (1) is defined as follows.

Definition 1. The relative degree vector at of the system (1) is

where is a neighborhood of .

The Fliess function expansions are an efficient way to express the input-output relation of a nonlinear system [29]. We denote and for simplicity. Then we can define the multi-integrals inductively as

The input-output relation of system (1) can be determined by its Fliess functional expansion. Thus, we have the following lemma. The multi-index integrations in (3) are defined recursively, via which we can define the Fliess functional expansion of system (1).

Lemma 2 (see [29]). The input-output response of system (1) can be expressed as

where and represents the Lie derivative of function about the vector field .

From Lemma 2 we can see; (4) is an infinite series functional expansion on the control inputs, where the control inputs enter the system via multi-index integrations.

2.2. Problem Formulation

Next, we will first convert an affine nonlinear system with candidate control modes of continuous inputs into a switched hybrid system and then give the formulation of the output tracking control problem for such a class of switched hybrid system.

An affine nonlinear system with candidate modes of control inputs can be described by

where is the state vector, is smooth vector field, are smooth functions, and, , is the continuous control input vector. Here, we assume that is known or has already been well designed for details refer to [21–23].

To convert the family of continues nonlinear affine systems into a switched hybrid system, we introduce a multiple-valued logic variable

where , , , and implies that the th continuous control input is activated.

To be more clear, we define a finite set , :

where , represents an identity matrix, and is the th column of identity matrix . The switching from to implies that the control mode switches from to , and correspondingly the piecewise continuous control input switches from to .

Then we can rewrite the above system (5) as follows

where ,  , , .

Choose the output of system (8) as

Then the output tracking problem under consideration can be defined as follows.

Definition 3. Consider system equations (8) and output equation (9), and let the reference output trajectory be

The output tracking problem in infinite horizon is to find an admissible switching control such that

where is defined in (10).

Next, we formulate the output tracking problem in a discrete-time form. Choose a time interval as the sampling time; assume within the sampling time the multiple-valued logic variable is constant and the possible switching instants are at

Then

where .

In the sequel, we simply denote

In (11), the cost function is defined in infinite time horizon. In this work, we concern the output tracking problem in a finite time horizon; in this case the cost function (11) becomes

where is a proper time horizon. Then purpose of the tracking problem is to find a sequence of , such that the cost function (15) is minimized.

In the cost function (11) or (15), only the tracking performance is considered. Actually, for application, the energy cost of the control inputs should be taken into account. Besides, the switching numbers also need to be restricted as well. A higher switching frequency will cause unstable operations or other problems. Thus, we will define two more comprehensive cost functions which incorporate these constraints.

Denote a column vector

and a matrix

where represents the dimension of continuous control .

Denote a control input vector

where .

Considering the control input energy constraints, we have the following cost function:

where and are symmetric positive-definite weighting matrices.

Remark 4. In order to guarantee the tracking performance, should be chosen related to . If it is too large with respect to , the tracking efficiency might be hurt.

When we concern about the switching frequency, a punishment term can be added in (19), and then the cost function with the penalty of the switching numbers is described as follows:

where , , , and are symmetric positive-definite weighting matrices.

Remark 5. Similarly like the choice of the weighting matrices in (19), the choice of affects the output tracking performance; and should be chosen related to . Generally, we choose and to be about while choosing to be 1, because the tracking performance is the most concerned.

Fliess functional expansion combined with multiple-step prediction will enhance the convergence speed. Next, we will formulate the output tracking problem for a switched hybrid system within prediction horizon.

Consider a nonlinear affine system with candidate modes of piecewise continuous control input, described as follows:

The output tracking problem formulated within a receding horizon control problem is as follows.

Definition 6. The output tracking problem within prediction horizon is, for the switched hybrid system (21) and the current state , to find a switching control sequence

such that the cost function (19) or (20) is minimized over the prediction horizon , where , , and are symmetric positive-definite weighting matrices.

As a special case, we give the output tracking problem of switched hybrid system (21) within one-step prediction horizon as follows.

Definition 7. The output tracking problem of switched hybrid system (21) within one-step prediction horizon is to find an admissible switching control

such that the cost function

or

is minimized, where , and are symmetric positive-definite weighting matrices.

3. Optimal Switching Control Design

In this section, we will give the design of the optimal switching control law via which the output tracking can be achieved. The main idea here is to express the output by using Fliess function expansion with admissible switching control. Then the cost function can be described by linear function of switching control inputs.

3.1. Optimal Switching Control in -Step Output Prediction

We consider the output tracking problem within predictive steps. Consider system (1), and assume that the continuous control inputs are constant within the sampling duration , that is,

Consider system (1), we have the following lemma about the Fliess functional expansion about , , , .

Lemma 8 (see [30]). Consider the nonlinear affine system (1), and assume the controls are constant at time interval , . Then at , corresponding to the coefficient

the integral in (4) is

where , .

Using Lemma 8, consider the nonlinear affine system (21), and assume that is the relative degree related to the th output. Then, we have the th order approximation of Fliess functional expansion being [31]

where is a linear function on , .

Remark 9. Since has finite choices, and is a linear function on , we can substitute (29) into (19) or (20). The performance index becomes a linear function of . Then the optimal switching control sequence , , can be obtained via dynamic programming.

Due to the problem of computation amount, the -step predictive output is actually not useful in application if both and are large. Next, we will discuss the optimal switching control design problem in one-step predictive case in detail.

3.2. Optimal Switching Control in One-Step Output Prediction

Choose a time duration as the sampling time, and let , . Consider system (1), and assume that the continuous control inputs are constant within the sampling duration , that is,

Then the output of system (1) can be given by the following lemma.

Lemma 10 (see [30]). The output of system (1) can be expressed via Fliess functional expansion as follows:

where is the relative degree related to the th output, , and .

From (31) we can see that is an infinite sequence about , . For convenience, we give the order approximation of the Fliess functional expansion (31) as follows [30]:

where .

The detailed expression of the and order approximation of the Fliess functional expansion (31) can be found in [30, 32].

Consider cost function (20)

Recall the switched hybrid system (21), in which . Thus, in terms of Lemma 10, using (32) to each , , we have

where , .

Define

Denote the optimal control by , recall the cost function (33), and substitute (36), (37), and (38) into (33) and then we have

Remark 11. Since , and has only finite choices at each moment, say choices, hence, the above problem (39) can be solved via dynamic programming. Considering generally is not too large, we can also use a simple brutal method to calculate the optimal control , which will not need a lot of computation amount.

4. Application in Multitarget Tracking Problem

In this section, we consider a multitarget tracking problem for a flight vehicle in the terminal guidance course. The proposed design approach is adopted.

4.1. Multitarget Tracking Problem Formulation

We consider a scenario in which one flight vehicle intercepts one real target with several decoys flying around it. These decoys appear as apparent targets before they are distinguished via the detective system. In this case, we have to consider the guidance problem of one flight vehicle versus multiple targets, otherwise a large miss distance between the flight vehicle and the real target might be caused due to the maneuverability limitation of the flight vehicle. Generally speaking, these decoys can be identified with the approaching of the flight vehicle and these targets. Such a problem is also a subproblem of multiple flight vehicles intercepting multiple targets.

First, we gave the dynamics of an engagement geometry between one flight vehicle and one target in a two-dimensional plane, shown in Figure 1. Let represent the flight vehicle and represent the target. LOS represents the line of sight. Variables associated with the flight vehicle or the target are denoted by an additional subscript or . The velocity, normal acceleration, and flight path angle are denoted by , , and , respectively.

Figure 1: Planar engagement geometry between one flight vehicle and one target.

The purpose of guidance is to stabilize the rotation of the LOS during the approaching between the flight vehicle and its target. When the approaching velocity is less than zero, the flight vehicle can intercept the target successfully once the velocity normal to the LOS converges to zero. Therefore, for a one-to-one guidance case, the guidance control for the flight vehicle is aiming to eliminate the relative velocity normal to the LOS. While considering the guidance problem of a one flight vehicle versus multiple targets, it is required to take into account each LOS rotation between the flight vehicle and each of the multiple targets in order to track multiple targets. The purpose of guidance in this case is to guarantee the relative velocity normal to the LOS between the flight vehicle and each of the multiple targets as small as possible until the real target is determined, which poses a challenging problem for the researchers in the guidance and control field.

For most existing guidance laws, the relative motion equations are usually described in the LOS coordinates. However, for the guidance problem of one flight vehicle versus multiple targets, we need to deal with the coupling relations between different LOS coordinates. Thus, it is convenient to describe these relative motion equations in a uniform framework, that is, the inertial coordinates, the detailed definition of the LOS coordinate and the inertial coordinate please refer to [33]. Label the targets through 1 to . The engagement geometry between the flight vehicle and targets in a two-dimensional inertial coordinates is shown in Figure 2.

Figure 2: One-to-many guidance engagement geometry.

Next, we will first give the relative motion equation between the flight vehicle and the th target, . Let and be the relative position between the flight and the th target and their the relative velocity, respectively. Then, the transformation of the variables can be described as follows:

where and represent the LOS angle between the flight vehicle and target and their relative distance, respectively, and represent the relative velocity along the LOS associated with the flight vehicle and target and their relative velocity normal to the LOS.

Choose the relative position and the relative velocity as the state variable

and choose the relative velocity normal to the LOS as the output of interest

Then, the relative motion equation and the output equation of the a flight vehicle and target can be described in a canonical form as follows:

where

Choose as the state variable and as the output of interest; then we can obtain the relative motion equations and output equation of the one flight vehicle versus targets problem

where

For a conventional guidance law, the orientation of the flight vehicle is fixed, that is, towards its aiming target. However, when we consider the guidance law design for one flight vehicle aiming for , targets, there are multiple relative motion information required to be dealt with. Suppose that the overload command at each time instant is a function of the relative information between the flight vehicle and some certain target; then different orientations of the flight vehicle lead to different targets indexes which will cause different overload commands and directions. Thus, in a one flight vehicle versus targets guidance problem, we need to design the functions between the overload command and the relative information associated with the flight vehicle and multiple targets, and the orientation of the flight vehicle as well. Adopting the approach proposed in the above section, we track multiple targets by switching the orientation among these targets so that these targets can be intercepted once it is determined to be the real aiming one. The switching guidance control scheme is shown in Figure 3.

Figure 3: The switching guidance scheme of one flight vehicle versus multiple targets.

Let , , be the orientation of the flight vehicle at time . If the flight vehicle chooses the target as its orientation at time , then . At each time instant, the overload command is the function of relative information between the flight vehicle and one individual target; thus the overload command can be written as

where represents the guidance law related to the relative motion information between the flight vehicle and the target whose index is chosen by the switching law. The aim of guidance command is to decrease the relative velocity normal to the LOS between the flight vehicle and chosen target. In this work, we adopt a commonly used guidance law, that is, the Proportional Navigation (PN) guidance law [33], to design . In the LOS coordinates, the Proportional Navigation guidance law for the flight vehicle intercepting target is

where is called navigation constant and in practice always chosen as a value between 3 to 5 to guarantee a good dynamic performance. In inertial coordinates, we have

Substituting (47) into (45) yields

The purpose of guidance problem for one flight vehicle versus multiple targets is to guarantee the relative velocities normal to the LOS between the flight vehicle and multiple targets are as small as possible. Meanwhile, we need to concern those targets which have higher probabilities of being real targets. Thus, we denote by , , the level of interest (LoI). The higher level of interest implies higher priority of the target or higher probability of the target being real target. Generally, the LoI can be provided by the detective system installed on the flight vehicle or can be obtained via some filter approaches; for more discussion on the LoI, please refer to [34, 35]. Thus, considering the system (50), we have the following cost function:

where is the LoI of target , is the relative velocity normal to the LOS between the flight vehicle and target , and are weighting coefficients, and

Calculate the relative degree of system (50) and we get , . In terms of (34), we have

where and are defined in (49),

Define

Denote the optimal control by , and substitute (56), (57), and (58) into (51); then we have

where is the matrix of the level of interest, described as follows:

Then can be obtained via dynamic programming. When is not too much large, we also can use brutal algorithm to obtain .

4.2. Simulation Results

Next, we will take a two-target tracking case as an example to show the effectiveness of the proposed approach, say . The setup of initial conditions is shown in Table 1.

Table 1: Initial conditions.

Initial position (m)
Initial velocity (m/s)

The flight vehicle

Target 1

Target 2

First, we consider the situation that the levels of interest of all the targets are the same, choose . The simulation time period is 10 seconds. The orientation of the flight vehicle is shown in Figure 4 and the relative velocity normal to the LOS of the flight vehicle is shown in Figure 5. From Figures 4 and 5 we can know that the flight vehicle keeps changing its target orientation to guarantee that the relative velocity normal to the corresponding LOS will not increase continuously.

Figure 4: Orientations of the flight vehicle without LoIs.

Figure 5: Relative velocity normal to the LOS without LoIs.

The flight trajectories that the flight vehicle aims to target 1 and target 2 individually respectively, and the flight trajectory that the flight vehicle switches its target orientation between target 1 and target 2 are shown in Figure 6.

Figure 6: Trajectories of the flight vehicle.

It is shown in Figure 6 that the trajectory that the flight vehicle takes into account, both of the targets, lies in the trajectories that the flight vehicle, respectively, aims for, target 1 and target 2, which implies that the flight vehicle can keep tracking both of the targets during the terminal guidance course.

Next, we consider the situation where the targets associated with different levels of interest. Suppose that target 2 has a higher priority, and choose , . The orientation of the flight vehicle is shown in Figure 7, and the relative velocity normal to the LOS is shown in Figure 8.

Figure 7: Orientations of the flight vehicle with different LoIs.

Figure 8: Relative velocity normal to the LOS with different LoIs.

From Figure 7 we can see, during the simulation period which is totally 10 seconds, there are about 3 seconds when the flight vehicles aim for target 1 and about 7 seconds when the flight vehicles aim for target 2. Correspondingly, the relation velocity normal to the LOS between the flight vehicle and the target 2 is smaller.

The trajectory of the flight vehicle is shown in Figure 9.

Figure 9: Flight trajectories of the flight vehicle considering different LoIs.

It can be seen from Figure 9 that when we consider the different priorities of the targets, the trajectory of the flight vehicle will be closer to the target with a higher level of interest, which implies that the flight vehicle can take into account both targets via the proposed method. Furthermore, the relative velocity normal to the LOS between the flight vehicle and the target with a higher level of interest has a faster convergence speed.

5. Concluding Remarks and Future Work

In this paper, the output tracking problem of a class of switched hybrid system is considered. By introducing a multiple-valued logic signal, a nonlinear affine system with multiple modes of candidate control inputs is converted into a switched hybrid system. The output tracking problem of such a class of switched hybrid system is formulated into a discrete-time form. Under the assumption that the possible switching only occurs at sampling time instant, the output of the switched hybrid system can be expressed via Fliess functional expansion. The optimal switching control can be determined for a multiple-step output tracking performance index by dynamic programming. Then we consider the guidance problem of a flight vehicle aiming to a real target with several decoys; the proposed approach is applied to a multitarget tracking problem. Simulation results show the effectiveness.

In the current work, we use an approximation of Fliess expansion for the output of the switched hybrid system. The truncated errors remain small if the time interval is chosen small enough. Actually, Fliess expansion can be finite terms for an input-output decoupled nonlinear affine system which requires in advance designing a decoupled feedback control. Future work will focus on these issues.

Acknowledgments

This research work is partially supported by the National Natural Science Foundation of China (Grant no. 61074160) and the Key Laboratory of Systems and Control, Chinese Academy of Sciences. The authors would gratefully acknowledge the support of the National Natural Science Foundation of China for this work.

References

J. A. Stiver and P. J. Antsaklis, “Modeling and analysis of hybrid control systems,” in Proceedings of the 31st IEEE Conference on Decision and Control, pp. 3748–3751, 1992.

P. Antsaklis, X. Koutsoukos, and J. Zaytoon, “On hybrid control of complex systems: a survey,” Journal Europeen des Systemes Automatises, vol. 32, no. 9-10, pp. 1023–1045, 1998.

M. S. Branicky, V. S. Borkar, and S. K. Mitter, “A unified framework for hybrid control: model and optimal control theory,” IEEE Transactions on Automatic Control, vol. 43, no. 1, pp. 31–45, 1998.

Z. Sun and S. S. Ge, “Analysis and synthesis of switched linear control systems,” Automatica, vol. 41, no. 2, pp. 181–195, 2005.

M. Lazar and W. P. M. H. Heemels, “Predictive control of hybrid systems: input-to-state stability results for sub-optimal solutions,” Automatica, vol. 45, no. 1, pp. 180–185, 2009.

M. S. Shaikh and P. E. Caines, “On the hybrid optimal control problem: theory and algorithms,” IEEE Transactions on Automatic Control, vol. 52, no. 9, pp. 1587–1603, 2007.

R. Goebel, R. G. Sanfelice, and A. R. Teel, “Hybrid dynamical systems: robust stability and control for systems that combine continuous-time and discrete-time dynamics,” IEEE Control Systems Magazine, vol. 29, no. 2, pp. 28–93, 2009.

X. Xu and P. J. Antsaklis, “Optimal control of switched systems based on parameterization of the switching instants,” IEEE Transactions on Automatic Control, vol. 49, no. 1, pp. 2–16, 2004.

L. Wu and J. Lam, “Sliding mode control of switched hybrid systems with time-varying delay,” International Journal of Adaptive Control and Signal Processing, vol. 22, no. 10, pp. 909–931, 2008.

H. Lin and P. J. Antsaklis, “Stability and stabilizability of switched linear systems: a survey of recent results,” IEEE Transactions on Automatic Control, vol. 54, no. 2, pp. 308–322, 2009.

L. Wu and W. X. Zheng, “Weighted model reduction for linear switched systems with time-varying delay,” Automatica, vol. 45, no. 1, pp. 186–193, 2009.

L. Wu, W. Zheng, and H. Gao, “Dissipativity-based sliding mode control of switched stochastic systems,” IEEE Transactions on Automatic Control, vol. 58, no. 3, pp. 785–793, 2013.

S. Ben Attia, S. Salhi, and M. Ksouri, “Static switched output feedback stabilization for linear discrete-time switched systems,” International Journal of Innovative Computing, Information and Control, vol. 8, no. 5, pp. 3203–3213, 2012.

X. Su, P. Shi, L. Wu, and Y.-D. Song, “A novel control design on discrete-time takagi-sugeno fuzzy systems with time-varying delays,” IEEE Transactions on Fuzzy Systems, article 1, no. 99, 2012.

R. Yang, P. Shi, G.-P. Liu, and H. Gao, “Network-based feedback control for systems with mixed delays based on quantization and dropout compensation,” Automatica, vol. 47, no. 12, pp. 2805–2809, 2011.

L. Wu, D. W. C. Ho, and C. W. Li, “Sliding mode control of switched hybrid systems with stochastic perturbation,” Systems & Control Letters, vol. 60, no. 8, pp. 531–539, 2011.

H. Lin and P. J. Antsaklis, “Stability and stabilizability of switched linear systems: a survey of recent results,” IEEE Transactions on Automatic Control, vol. 54, no. 2, pp. 308–322, 2009.

F. Zhu and P. J. Antsaklis, “Optimal control of switched hybrid systems: a brief survey,” Tech. Rep. ISIS-2011-003, 2011.

F. Borrelli, M. Baotić, A. Bemporad, and M. Morari, “Dynamic programming for constrained optimal control of discrete-time linear hybrid systems,” Automatica, vol. 41, no. 10, pp. 1709–1721, 2005.

D. Görges, M. Izák, and S. Liu, “Optimal control and scheduling of switched systems,” IEEE Transactions on Automatic Control, vol. 56, no. 1, pp. 135–140, 2011.

X. Xu and P. J. Antsaklis, “Optimal control of switched systems via non-linear optimization based on direct differentiations of value functions,” International Journal of Control, vol. 75, no. 16-17, pp. 1406–1426, 2002.

X. Xu and P. J. Antsaklis, “Optimal control of switched autonomous systems,” in Proceedings of the 41st IEEE Conference on Decision and Control, vol. 4, pp. 4401–4406, December 2002.

X. Xu and P. J. Antsaklis, “Results and perspectives on computational methods for optimal control of switched systems,” in Hybrid Systems: Computation and Control, pp. 540–555, Springer, Berlin, Germany, 2003.

S. C. Bengea and R. A. DeCarlo, “Optimal control of switching systems,” Automatica, vol. 41, no. 1, pp. 11–27, 2005.

S. Wei, K. Uthaichana, M. Žefran, R. A. DeCarlo, and S. Bengea, “Applications of numerical optimal control to nonlinear hybrid systems,” Nonlinear Analysis. Hybrid Systems, vol. 1, no. 2, pp. 264–279, 2007.

T. Das and R. Mukherjee, “Optimally switched linear systems,” Automatica, vol. 44, no. 5, pp. 1437–1441, 2008.

D. Nešić and L. Grüne, “Lyapunov-based continuous-time nonlinear controller redesign for sampled-data implementation,” Automatica, vol. 41, no. 7, pp. 1143–1156, 2005.

L. Grüne and P. E. Kloeden, “Higher order numerical schemes for affinely controlled nonlinear systems,” Numerische Mathematik, vol. 89, no. 4, pp. 669–690, 2001.

A. Isidori, Nonlinear Control Systems, Communications and Control Engineering Series, Springer, Berlin, Germany, 3rd edition, 1995.

Y. Yao, B. Yang, F. He, Y. Qiao, and D. Cheng, “Attitude control of missile via fliess expansion,” IEEE Transactions on Control Systems Technology, vol. 16, no. 5, pp. 959–970, 2008.

Y. Qiao, B. Yang, and D. Cheng, “Attitude control of missile via Fliess expansion and model predictive control,” in Proceedings of the 7th World Congress on Intelligent Control and Automation WCICA '08), pp. 1527–1532, June 2008.

Y. Yao, B. Yang, F. He, and D. Cheng, “Fliess expansion-based bang-bang control design and its application to attitude control of missile,” in Proceedings of the 6th World Congress on Intelligent Control and Automation (WCICA '06), pp. 188–192, June 2006.

P. Zarchan, Tactical and Strategic Missile Guidance, AIAA, San Diego, Calif, USA, 1997.

Y. Yao, P. Zhang, H. H. T. Liu, and F. He, “Optimal sweep-based persistent surveillance using multiple unmanned aerial vehicles with level of interest,” in Proceedings of the 10th World Congress on Intelligent Control and Automation, pp. 2441–2446, Beijing, China, 2012.

Y. Yao, P. Zhang, H. H. T. Liu, and F. He, “Optimal switching target-assignment based on the integral performance in cooperative tracking,” Science China Information Science, vol. 56, no. 1, pp. 1–14, 2013.