Journal of Applied Mathematics

Volume 2014 (2014), Article ID 391617, 9 pages

http://dx.doi.org/10.1155/2014/391617

Research Article

Sign Stability for Switched Linear Systems and Its Application in Flight Control

Qing Wang,¹ Tong Wang,¹ Yanze Hou,² and Chaoyang Dong³

¹Department of Automation Science and Electrical Engineering, Beijing University of Aeronautics and Astronautics, Beijing 100191, China

²Institute of Manned Space System Engineering, China Academy of Space Technology, Beijing 100094, China

³Department of Aeronautic Science and Engineering, Beijing University of Aeronautics and Astronautics, Beijing 100191, China

Received 8 December 2013; Revised 21 April 2014; Accepted 23 April 2014; Published 13 May 2014

Academic Editor: Yuxin Zhao

Copyright © 2014 Qing Wang 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 sign stability concept in ecological systems is introduced into the analysis and synthesis of switched linear system to explore new control design technique. The necessary and sufficient condition for sign stability of a switched linear system under arbitrary switching is achieved via the notion of complete isogenous sign stable set (CISSS). A new approach for the stabilization of switched system is presented. Although the controllers are devised for each subsystem, respectively, the switched system is sign stabilized by the constitution of CISSS. The provided method has natural robustness and more design freedoms than the familiar Lyapunov function method, which bears relative conservativeness as the requirement of solving LMIs. The presented technique is validated by an example of flight control within a large-scale flight envelop. Simulation results indicate that the proposed method can stabilize the flight attitude under large variations of system parameters and external perturbations.

1. Introduction

The sign stability (or qualitative stability) concept is first proposed in bionomics and utilized to analyze the interactions of different species in a large-scale ecosystem which lacks exact model but presents high robustness under various perturbations [1, 2]. The sign stability approach is also applied to population biology and economics in respect that these systems are also short of quantitative mathematical models. This qualitative analysis technique provides an avenue to research the linear system stability by the Jacobian matrix with only signs, and still attracts growing attention in system science community [3, 4].

For a given matrix, the signs of its elements are taken to make up a new matrix named the sign-pattern. A matrix is called sign stable if arbitrary matrix which has the same sign-pattern is Hurwitz stable, regardless of the elements’ magnitudes. Hence, the sign stability of a matrix equates to that of the corresponding sign-pattern. Jeffries developed necessary and sufficient conditions for sign stability of matrix and proposed an approach named the color test to verify an arbitrary matrix in ecological terms [5, 6]. Yedavalli translated the color test conditions in matrix theory notation and devised a programmable set of conditions for the color test in an irreducible matrix [7]. It is worth noticing that the conclusions in this work are based on the conditions and criterions for the matrix sign stability presented in the aforementioned literatures.

Based on theoretical researches, the sign stability technique is applied in control engineering, especially the robust controllers design for aerospace flight control system. A three-axis attitude stabilization controller for an axisymmetric satellite is provided in [7]. The closed-loop system matrix is designed to possess the specific sign-pattern and sign stability property which will bring the system robust stability under arbitrarily large variations in the spin angular velocity. In [8], the effect of the elements’ signs on the matrix properties such as eigenvalues and condition number is shown. Efforts are made to identify target closed-loop systems that incorporate the desirable features of ecological systems, and an algorithm for the design of controller is given. The control design procedure is illustrated with two applications in the aerospace field: satellite attitude control and aircraft lateral dynamics control.

However, the existing applications of sign stability only consider linear time-invariant systems which cannot approximate the flight dynamics within the full flight envelope. One of the alternative solutions is the switched linear system which is a good approximation of complicated system characteristics. Flight-control-oriented analysis and synthesis of switched systems are studied in the prior literatures and show large potential in engineering practice [9–12]. For example, the gain schedule control that is widely applied in flight control systems can be abstracted as a switching control law if variations of system parameters are regarded as a series of switching with the transition of the schedule variable. Hence, the aim of this paper is to generalize the sign stability concept to switched linear system and develop a new and practicable control scheme.

Compared to the prior references, the significant contribution of this work is the conclusion on sign stability analysis of switched linear system, which is presented by a necessary and sufficient condition. The system model treated with sign stability approach is extended from LTI system to a hybrid system. The mix of continuous linear systems and discrete signals make switched system suitable to describe a wide range of engineering systems such as power systems, automotive engine control, flight control, and networked control systems. As a result, the application areas of sign stability approach are enlarged remarkably.

Another important contribution is that a new stabilization technique of switched system is presented. As for the asymptotic stability of switched system under arbitrary switching, the common Lyapunov function is a necessary and sufficient condition and a usual method. In [13, 14], the existing condition of common quadratic Lyapunov function for switched linear system is discussed. Furthermore, multiple Lyapunov functions method [15], dwell-time method [16], and average dwell-time method [17] are proposed to adapt different design situations. In recent years, Zhao and his coworkers propose the mode-dependent average dwell-time method in order to decrease conservativeness [18, 19]. To utilize the above methods, the LMI technique is widely used, whereas the feasible solution of LMI is only a sufficient condition for the existence of common Lyapunov function, and an expectant controller may not exist or may not be found even if one does exist [20]. That leads to uncertainty of the controller solvability, even though modern mathematic tools are utilized. In contrast, the sign stability theory offers a fire-new approach with which designers can configure target sign-patterns of closed-loop switched subsystems via state feedback and guarantee the stability under arbitrary switching. The primary advantage of the sign stability technique is that there is no need to solve LMIs. That provides more design freedoms than Lyapunov function approach. Besides, the provided control synthesis technique possesses natural robustness because the sign stability is independent of the elements’ magnitudes. This property allows the sign stabilization controllers to fit discretionarily large sign-preserving parameter perturbations, especially in the flight control system.

With this consideration, the paper is organized as follows. In Section 2, the basic concepts and conditions for the sign stability of matrix are reviewed. The definition of sign stability for switched linear system is presented to formulate the main issue of this paper. In Section 3, the main result of this paper is given in the form of necessary and sufficient condition for the sign stability of switched system along with a new approach for the stabilization problem. In Section 4, the proposed technique is illustrated by an application to flight control of the HiMAT vehicle. Finally, Section 5 provides a detailed conclusion of this work.

2. Problem Formulation

In this section, preliminary knowledge about the sign stability is reviewed, and then the definition of sign stability for switched linear system is proposed to formulate the main problem of this paper.

2.1. Preliminary Knowledge

A switched linear autonomous system is usually described as

in which is the switching signal and is the set of switching subsystem matrices. Then the aim of this section is to establish the concept of sign stability for the switched system (1). Above all, basic knowledge and results are needed.

For the real matrix and , if , for all , it is defined that the two matrices have the same sign-pattern. The sign-pattern can be represented as a matrix, of which all the entries simply consist of signs , or 0. The is defined as the sign-pattern matrix of . Similarly, is defined as the sign-pattern set of the matrix set ; that is, . For two sets of matrices and , if for all and for all , it is defined that the two matrix sets have the same sign-pattern set.

Definition 1. The real matrix is sign stable if each matrix of the same sign-pattern as is Hurwitz stable.

The necessary and sufficient condition for the sign stability of matrix is provided in referenced researches [5, 21–23]. The main necessary conditions are listed in Lemma 2 as the foundation of main results in this paper.

Lemma 2. The following are necessary conditions for sign stability of matrix .

(1)For all .(2) for at least one .(3)For all ,   .(4) for any sequences of three or more distinct indices .(5).

It is worth announcing that the color test, given by Jeffries in [5], is confirmed as the criterion for the sign stability of any matrix or sign-pattern in this work.

2.2. Problem Statement

Compared with the definition of sign stability for matrix, sign stability of switched linear system is defined as below.

Definition 3. The switched linear system (1) is sign stable under arbitrary switching if each switched system which has the same sign-pattern set as is asymptotically stable under arbitrary switching.

From Definition 3 it is seen that, for a given switched system, that is, sign stable, it is obviously asymptotically stable under arbitrary switching. Besides, the determinant of sign stability is the sign-pattern set, and the concept of sign stability for the sign-pattern set is similar to that for the switched system. We can also know that a sign stable switched system will be achieved by assigning arbitrary values for elements of sign-patterns in a sign-pattern set which is sign stable. In addition, the number of subsystems will not affect the stability of the switched system.

Remark 4. In a sign stable switched system, there may be several subsystems of which the matrices have the same sign-pattern. In other words, the switched system derived from a given sign stable sign-pattern set is stable even if more than one matrices are built according to the same sign-pattern.

Remark 5. The subset of an arbitrary sign stable switched system (or sign-pattern set) is also sign stable. The conclusion is obtained under arbitrary switching laws.

With the concept given in Definition 3, we are now interested in the properties of a sign stable switched linear system, and it is also important to develop a technique to constitute a sign stable switched system. Towards the above targets, the main results are provided in the next section. To accomplish the proof of the main theorem, a primary conclusion given by [24] is used and presented here as a lemma.

Lemma 6. The following statements are equivalent:

(1)the switched linear system

 where is asymptotically stable under arbitrary switching;(2)the linear time-variant system

is asymptotically stable.

3. Main Results

3.1. Necessary Condition

From Definition 3 and Lemma 6, the sign stability of the switched linear system (2) under arbitrary switching implies that each sign-pattern in is sign stable. The necessary conditions of a sign stable switched linear system are given in the following theorem, where denotes the element of the subsystem matrix of switched system (2).

Theorem 7. The following are necessary conditions for sign stability of switched linear system (2) under arbitrary switching.

(1)For all , which is the subsystem matrix of switched system (2), the necessary conditions in Lemma 2 are satisfied.(2)For all , for all , one of the following three conditions is satisfied: (i) ; (ii) ; (iii) .

Proof. For (1), noting that each subsystem of switched system (2) is sign stable, the conclusion is obvious.

For (2), with no loss of generality, it is supposed that there are 2 subsystems in switched system (2). From Lemma 6, the linear time-variant system (3) is asymptotically stable, and the system matrix can be written as

It is easily understood that is sign stable because the elements of and may take arbitrary values. Define as the element of , and it holds that

Due to condition (3) of Lemma 2, the following inequality holds for all :

(i) For the situation that , since and may take arbitrary values, it must hold that or . That is, for all , for all , or is satisfied.

(ii) For the situation that , with no loss of generality, supposing gives

That is,

Considering that all the elements permit arbitrary sign-preserving variations, let

which leads to

with denoting the sign of .

From condition (3) of Lemma 2, it is known that and have the opposite signs (including 0), so do and . It implies that and have the same sign. That is, . Then (8) can be rewritten as

That is,

To ensure (12), and must have the opposite signs. Therefore, and have the same sign (including 0). That is, for all , for all , .

Synthesizing situations (i) and (ii), condition (2) is proved.

It can be concluded by Theorem 7 that, for all the subsystem matrices of a sign stable switched system, the elements in the same position must satisfy specific conditions. According to this, the definition of an isogenous sign-pattern set and correlative corollary are given as below.

Definition 8. A sign-pattern set is called an isogenous set if the elements in the same position of all the sign-patterns satisfy one of the following terms: (i) they have the same sign (including 0); (ii) they are all 0; (iii) the elements in the opposite position are all 0.

Corollary 9. The sign-pattern set of a sign stable switched system is an isogenous set.

Remark 10. We also say that two sign-patterns are isogenous with each other if they conform to the condition of Definition 8. Apparently, if arbitrary two sign-patterns of a known set are isogenous, this set is an isogenous set. Corollary 9 can be acquired from Theorem 7 and Definition 8 easily. However, it should be noticed that the isogenous set is only a necessary condition for sign stable switched system. For example, an isogenous set by is in (13). However, the set is not sign stable as the first sign-pattern is not sign stable obviously:

Considering an arbitrary isogenous sign-pattern set, we define the joining operation of all the sign-patterns to develop a special sign-pattern named the original sign-pattern of the set. The joining operation should keep the following rules.

(1)If elements in a certain position of all the sign-patterns have the same nonzero sign or 0, then the corresponding element of the original sign-pattern is the same nonzero sign.(2)If elements in a certain position of all the sign-patterns are all 0, then the corresponding element of the original sign-pattern is 0.(3)If elements in a certain position of all the sign-patterns have different signs, then the corresponding element of the original sign-pattern is written as to denote arbitrary signs.

The original sign-pattern obtained from a known isogenous set may be a new sign-pattern or an existent one which is already in the set. The original sign-pattern of set (13) is given in as follows as an illustration:

Furthermore, it is obvious that the new original sign-pattern and the known isogenous set can be joined to build a new set which is still isogenous. Actually, new sign-patterns can be derived to extend the isogenous set by replacing some nonzero elements of the original sign-pattern by 0. To demonstrate the aforementioned operation, the set which is developed from (13) is given below with the first matrix as the original sign-pattern and the last two as accessorial ones:

Of course, we can continue to add new isogenous sign-patterns of (15) on the basis of the original sign-pattern and list all the other ones. However, among all these sign-patterns, the sign stable ones are what we focused on and can be selected out. That leads to another important definition in this paper.

Definition 11. A sign-pattern set is named a complete isogenous sign stable set (CISSS) if it contains a sign stable original sign-pattern and all of the derivative isogenous sign-patterns that are sign stable.

Remark 12. Definition 11 indicates that we can obtain a CISSS from an isogenous set, for example, (15), if the original sign-pattern and some other ones are sign stable. Unfortunately, as the original sign-pattern of (15) is not sign stable, the subset of (15) cannot be a CISSS. It can be seen that a sign stable original sign-pattern is necessary for a CISSS. This notion can be interpreted by another example. Set (16) is an isogenous set, but not a CISSS even though all the sign-patterns are sign stable. Actually, the original sign-pattern given by (17) is not included in (16), and it is more important that the original sign-pattern is not sign stable:

Remark 13. It also should be emphasized that the selection of the uncertain element does not affect the sign stability, which is related to the products of the off-diagonal pairs. Since the sign must appear in the form of an off-diagonal pair “0 and ,” leading the product to 0 at all time, the above pair can still exist in CISSS.

3.2. Necessary and Sufficient Condition

In what follows, the main result of this work is proposed.

Theorem 14. The necessary and sufficient condition for sign stability of switched linear system (2) under arbitrary switching is that the sign-pattern set of (2) is a subset of a CISSS.

Proof. (1) The necessity can be easily proved by Corollary 9 and Definition 11.

(2) If the sign-pattern set of (2) is a subset of a CISSS, then all the sign-patterns of are isogenous and sign stable. For another arbitrary matrix set that has the same sign-pattern set and , let

According to Definition 11, it is certain that belongs to the foregoing CISSS. Hence, is sign stable and thus Hurwitz stable. By Lemma 6, the switched system is asymptotically stable. At last, by Definition 3, the sufficiency is proved.

With Theorem 14, it is easy to judge whether a given switched system is sign stable. To construct a sign stable switched system, the original sign-pattern and the CISSS are required first. Furthermore, a sign stable original sign-pattern with the least element “0” is the basis to constitute a CISSS.

3.3. Simple Examples

Two simple examples are given to illuminate the constitution of CISSS.

Example 1. For matrices, the sets given by (19) are CISSS with the first sign-pattern as the original sign-pattern of each set, respectively:

It has been summarized that the original sign-pattern is needed to have as few “0” elements as possible. Hence, all the diagonal elements are chosen to be “−,” and the off-diagonal elements are selected to possess as many pairs of opposite signs as possible. For matrices, it is not necessary for us to configure any “0” element in the original sign-pattern as shown in Example 1, whereas, for matrix , according to condition of Lemma 2, there must be at least two “0” elements in the off-diagonal positions to guarantee and . It then can be concluded that, in the matrix below, there must be at least one 0 in the positions represented by and one 0 in the positions represented by , respectively:

Example 2. Consider the sign stable original sign-pattern given in (21), where each presents an uncertain sign that may be , or 0. The aim is to discover the corresponding CISSS:

Utilizing the criterion in [5] to select sign stable sign-patterns in the isogenous set of (21), the CISSS is acquired and shown in (22), with each denoting an arbitrary sign. As each sign-pattern in (22) contains sign-patterns, there is a total of sign-patterns in the CISSS:

In fact, in the above process of finding CISSS, the two pairs of “0 and ” do not need to be changed. The first term of the diagonal is also unchangeable. As a result, the only way is to transform the block on the lower-right. Indeed, it can be seen that the lower-right blocks are the same as the set of Example 1.

4. Example

The controller design problem for the highly maneuverable technology (HiMAT) vehicle is presented to demonstrate the main results in this paper. As declared in [10], a switched system can be applied to describe flight conditions with fast varying parameters. Controllers for the subsystems can be synthesized via a switching law to adapt the parameter variations among different operating points. In this section, the sign stability approach is utilized to design a switching control scheme for the unstable longitudinal dynamics. The considered flight envelope and 14 operating points within it are depicted in Figure 1. We suppose that the dynamics in the contiguous region of each operating point can be approximated by the corresponding linear subsystem. These linear models are given in [25].

Figure 1: Partial flight envelope and operating points.

The linearization model for each operating point is regarded as a subsystem of the switched system. The switched system is then given by

where is the state vector consisting of three state variables (angle of attack), (pitch rate), and (pitch angle), and is the control vector consisting of (elevator), (elevon), (canard), and (symmetric aileron).

The system matrices and control matrices for several operating points are shown to illustrate the sign-patterns:

Furthermore, the sign-pattern set of the open-loop system matrices is given by the following:

As explained in Section 3, a subset of a CISSS is required for the closed-loop switched system matrices. According to the structure of , the original sign-pattern is chosen as follows:

Although the complete form of the CISSS derived from (26) has more sign-patterns, we can regard (26) as the target sign-pattern of closed-loop subsystems directly. Compared to the sign-patterns in (25), it is only needed to change two signs in each subsystem matrix. A state feedback control gain is designed for each operating point in the envelope, and the results for several points are given as follows:

It can be seen that the closed-loop system matrices , , and have different sign-patterns. However, these sign-patterns (also including those of other operating points) belong to the CISSS derived from the original sign-pattern (26). Based on Theorem 14, the closed-loop switched system is sign stable under arbitrary switching.

The simulation across different regions of the flight envelope is taken to validate the proposed technique. The flight trajectory travels through vicinities of 7 operating points (; see Figure 1). The switching law is depicted in Figure 2.

Figure 2: Switching law of the flight trajectory.

To validate the stability within the full flight envelope, gust disturbance is considered in the form of perturbations acted on the angle of attack, which is a group of pulse signals with a period of 10 seconds, a width of 1 second, and a magnitude of 3 degrees. The initial values of the angle of attack, pitch rate, and pitch angle are set to 3 degrees, 3 degrees per second, and 5 degrees, respectively. To illustrate the robustness of the proposed approach, multiplicative perturbations on the closed-loop subsystem matrices are introduced. Each element in the matrices is multiplied by an arbitrary positive value, which will preserve all the signs of the products. The perturbation matrices are given by the following:

The resulted closed-loop subsystem matrices are obtained by multiplying each element of the original matrix with the corresponding element in the perturbation matrix. Thus, it leads to 3 closed-loop switched systems with different multiplicative perturbations. The time histories of the system states and control surface deflections are depicted in Figure 3.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 3: Time histories of system states and control inputs under perturbations: (a) angle of attack, (b) pitch rate, (c) pitch angle, (d) elevator, (e) elevon, (f) canard, and (g) symmetric aileron.

As depicted in Figure 3, the system states converge quickly under gust disturbances and keep stable while the switching occurs along the flight trajectory. There are vibrations in the deflections of the control surfaces at the time of the first switching, nevertheless the amplitudes are acceptable. It is concluded that the devised approach is robust for the sign-preserving parameter perturbations on the closed-loop systems.

5. Conclusion

A new concept of sign stable switched linear system is established to develop novel techniques for flight control. The main result is provided in the form of necessary and sufficient condition for the sign stability of switched system under arbitrary switching. Hence, the application areas of sign stability approach are enlarged remarkably. A new robust stabilization approach for switched linear system is proposed via the notion of CISSS. The sign stabilization controllers are devised for each subsystem, respectively; therefore the design process possesses more freedoms. Compared with the traditional Lyapunov function method, the proposed technique has natural robustness and decreased conservativeness. The aforementioned points are verified by a stabilization control problem of aircraft with a large-scale flight envelope.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work is supported by the National Natural Science Foundation of China (61273083).

References

R. M. May, “Will a large complex system be stable?” Nature, vol. 238, no. 5364, pp. 413–414, 1972.

L. Edelstein-Keshet, Mathematical Models in Biology, McGraw-Hill, New York, NY, USA, 1988.

D. O. Logofet, “Stronger-than-Lyapunov notions of matrix stability, or how “flowers” help solve problems in mathematical ecology,” Linear Algebra and its Applications, vol. 398, pp. 75–100, 2005.

D. A. Grundy, D. D. Olesky, and P. van den Driessche, “Constructions for potentially stable sign patterns,” Linear Algebra and its Applications, vol. 436, no. 12, pp. 4473–4488, 2012.

C. Jeffries, V. Klee, and P. van den Driessche, “When is a matrix sign stable?” Canadian Journal of Mathematics, vol. 29, no. 2, pp. 315–326, 1977.

C. Jeffries, “Qualitative stability and digraphs in model ecosystems,” Ecology, vol. 55, no. 6, pp. 1415–1419, 1974.

R. K. Yedavalli, “Robust control design for linear systems using an ecological sign-stability approach,” Journal of Guidance, Control, and Dynamics, vol. 32, no. 1, pp. 348–352, 2009.

R. K. Yedavalli and N. Devarakonda, “Sign-stability concept of ecology for control design with aerospace applications,” Journal of Guidance, Control, and Dynamics, vol. 33, no. 2, pp. 333–346, 2010.

Y. Hou, Q. Wang, and C. Dong, “Gain scheduled control: switched polytopic system approach,” Journal of Guidance, Control, and Dynamics, vol. 34, no. 2, pp. 623–628, 2011.

Y. Hou, C. Dong, and Q. Wang, “Stability analysis of switched linear systems with locally overlapped switching law,” Journal of Guidance, Control, and Dynamics, vol. 33, no. 2, pp. 396–403, 2010.

T. Shimomura, “Hybrid control of gain-scheduling and switching: a design example of aircraft control,” in Proceedings of the American Control Conference, pp. 4639–4644, Denver, Colo, USA, June 2003.

B. Lu, F. Wu, and S. Kim, “Switching LPV control of an F-16 aircraft via controller state reset,” IEEE Transactions on Control Systems Technology, vol. 14, no. 2, pp. 267–277, 2006.

K. S. Narendra and J. Balakrishnan, “A common Lyapunov function for stable LTI systems with commuting -matrices,” Institute of Electrical and Electronics Engineers. Transactions on Automatic Control, vol. 39, no. 12, pp. 2469–2471, 1994.

D. Liberzon, J. P. Hespanha, and A. S. Morse, “Stability of switched systems: a Lie-algebraic condition,” Systems & Control Letters, vol. 37, no. 3, pp. 117–122, 1999.

M. S. Branicky, “Multiple Lyapunov functions and other analysis tools for switched and hybrid systems,” Institute of Electrical and Electronics Engineers. Transactions on Automatic Control, vol. 43, no. 4, pp. 475–482, 1998.

A. S. Morse, “Supervisory control of families of linear set-point controllers. I. Exact matching,” Institute of Electrical and Electronics Engineers. Transactions on Automatic Control, vol. 41, no. 10, pp. 1413–1431, 1996.

J. P. Hespanha and A. S. Morse, “Stability of switched systems with average dwell-time,” in Proceedings of the 38th IEEE Conference on Decision and Control (CDC), pp. 2655–2660, Phoenix, Ariz, USA, December 1999.

X. Zhao, L. Zhang, P. Shi, and M. Liu, “Stability and stabilization of switched linear systems with mode-dependent average dwell time,” Institute of Electrical and Electronics Engineers. Transactions on Automatic Control, vol. 57, no. 7, pp. 1809–1815, 2012.

Q. Lu, L. Zhang, H. R. Karimi, and Y. Shi, “ control for asynchronously switched linear parameter-varying systems with mode-dependent average dwell time,” IET Control Theory & Applications, vol. 7, no. 5, pp. 673–683, 2013.

O. Mason, R. N. Shorten, and D. Leith, “Evaluation of piecewise quadratic methods applied to the wind-turbine benchmark example,” Tech. Rep. NUIM-SS-2001-11, Signals and Systems Group, NUI Maynooth, 2001.

V. Klee and P. van den Driessche, “Linear algorithms for testing the sign stability of a matrix and for finding -maximum matchings in acyclic graphs,” Numerische Mathematik, vol. 28, no. 3, pp. 273–285, 1977.

J. Quirk and R. Ruppert, “Qualitative economics and the stability of equilibrium,” The Review of Economic Studies, vol. 32, no. 4, pp. 311–326, 1965.

R. May, Stability and Complexity in Model Ecosystems, Princeton University Press, Princeton, NJ, USA, 1973.

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

G. L. Hartmann, M. F. Barrett, and C. S. Greene, “Control design for an unstable vehicle,” Tech. Rep. NASA CR-170393, 1979.