x 0 := (d CA 1 Bκ) P 1
|
|
- Valentine Porter
- 5 years ago
- Views:
Transcription
1 Global Asymptotic tability of the Limit Cycle in Piecewise Linear versions of the Goodwin Oscillator Adrián A. alinas-varela, Guy-Bart tan, Jorge Gonçalves Department of Engineering, University of Cambridge. {aas47, gvs, Abstract: Conditions in the form of linear matrix inequalities (LMIs) are used in this paper to guarantee the global asymptotic stability of a limit cycle oscillation for a class of piecewise linear (PWL) systems defined as the feedback interconnection of a saturation controller with a single input, single output (IO) linear time-invariant (LTI) system. The proposed methodology extends previous results on impact maps and surface Lyapunov functions to the case when the sets of expected switching times are arbitrarily large. The results are illustrated on a PWL version of the Goodwin oscillator. 1. INTRODUCTION Motivated by the pervasiveness of dynamical systems that exhibit stable limit cycle oscillations (see Goldbeter [1996], Mosekilde [1997], and trogatz [003] for numerous examples), this paper focuses on the global analysis of limit cycle oscillations of piecewise linear (PWL) systems. uch systems are defined by a set of affine linear systems ẋ = A σ x + B σ (1) where x R n is the state, together with a piecewise constant rule to switch among them σ(x) {1,...,M} () that depends on present values (and possibly also on past values) of x. This kind of systems are typically used to model processes which encompass several modes of operation with a different (but linear) dynamical behaviour in each mode. Unfortunately, the majority of the analytical tools devised specifically for this kind of systems - whether based on extensions of Lyapunov s theory (e.g. Branicky [1998], Johansson and Rantzer [1998], Petterson and Lennartson [00], Prajna and Papachristodoulou [003]) or not (e.g Margaliot [006], Iwatani and Hara [006], Xu and Antsaklis [1999], Boscain [00]) - is incompatible with the study of limit cycle oscillations. This apparent limitation was addressed by Gonçalves [000]. In the latter reference, the stability analysis of a PWL system s limit cycle is achieved via the use of linear matrix inequalities (LMIs) to construct Lyapunov functions on the switching surfaces (the subsets of the state space in which the system is allowed to switch from one mode of operation to another). Unlike the phase plane methods and Poincaré-Bendixson theorem, this methodology generalises easily to high dimensional systems. It has been successfully applied to the global analysis of oscillations in relay feedback systems [Gonçalves, 000, chapter 5] and to the characterisation of regions of stability of limit cycles for more general PWL systems in Gonçalves [005]. Nevertheless, the construction of the Lyapunov functions depends upon the impact maps (maps from one switching surface to another), which are in turn characterised by their associated switching times (the time between switches). Consequently, it is impracticable when such switching times are arbitrarily large. This paper is devoted to the analysis of a specific class of systems in which such a situation arises. This work was supported in part by the Mexican Council of cience and Technology (CONACyT), the Cambridge Overseas Trust and the Overseas Research tudentships Programme. This paper is organised as follows. ection introduces the class of systems that we intend to analyse. ection 3 starts by stating the conditions which guarantee the asymptotic stability of the system s limit cycle in the entirety of its state space with the exception of its stable manifold. ince the presence of this stable manifold implies the existence of arbitrarily large switching times, conditions to overcome such obstacle are then presented. A PWL version of the Goodwin oscillator is used in ection 4 to illustrate our results. ection 5 concludes and gives directions for future research.. PROBLEM DEFINITION Consider the feedback interconnection of a IO LTI system satisfying the following linear dynamic equations ẋ = Ax + Bu (3) y = Cx where x R n, with a saturation controller defined as { κ if y(t) < d u(t) = κ/d y(t) if y(t) d (4) κ if y(t) > d where d > 0 and κ > 0. By a solution of (3) - (4) we mean functions (x(t),y(t),u(t)) satisfying (3) - (4). As discussed in [Gonçalves, 000, chapter 7.], this system is symmetric around the origin and has a unique solution for any initial state. In the state space, the saturation controller introduces two switching surfaces consisting of hyperplanes of dimension n 1: := {x R n Cx = d}, := {x R n Cx = d} On one side of the switching surface (Cx > d), the system is governed by ẋ = Ax Bκ. In between the two switching surfaces ( Cx d), the system is given by ẋ = (A (κ/d)bc) x. Finally, on the other side of (Cx < d), the system is governed by ẋ = Ax+Bκ. Define the subsets + and of as follows: + := {x C (A (κ/d)bc) x 0} := {x C (A (κ/d)bc) x 0} As shown in Figure 1(b), + ( ) is the set of points in that can be reached by trajectories of (3) - (4) when governed by the subsystem ẋ = (A (κ/d)bc) x (ẋ = Ax Bκ). Define also + := + and :=, where x s denotes x s {x x }..1 Class of systems under consideration. ystems of the form (3) - (4) are nonlinear and as such, they can exhibit extremely complex behaviours. ome may
2 _ u IO LTI y (a) aturation system Fig. 1. Problem definition ẋ = Ax B κ + ẋ = (A (κ/d)bc) x 0 (b) The sets +, C x > d d < C x < d be chaotic; others may have have several isolated equilibria; others might have limit cycles, or even some combination of all these behaviours. Here we consider a class of systems which satisfies the following conditions: There exists a periodic solution with 4 switches per cycle, i.e. there exists a set of time instants {t 1,t,t 3,t 4} which complies with [Gonçalves, 000, Proposition 3.]. The periodic solution is locally stable, i.e. the conditions in [Gonçalves, 000, Proposition 3.3] are satisfied. The IO LTI system (3) is of dimension 3. The matrix A in (3) is Hurwitz. The inequality κca 1 B d holds, ensuring that the origin is the only equilibrium point. The matrix (A (κ/d)bc) is assumed to have a real negative eigenvalue and a pair of complex conjugate ones with positive real parts. The last condition implies that the unstable equilibrium x = 0 possesses a stable manifold, which we will denote as M s. It is readily seen that solutions of (3) - (4) belonging to M s will approach the origin asymptotically. But what happens when this is not the case? Will all other solutions eventually converge to the limit cycle described by the first two conditions? This is what we understand as global asymptotic stability of the limit cycle and constitutes the question that we propose to answer.. tructural description. The following discussion examines the implications entailed by the assumption regarding the eigenvalues of (A (κ/d)bc). This will play an important role when verifying a set of sufficient conditions which guarantee that the limit cycle is asymptotically stable in R 3 \ M s. Let B be a basis of R 3. uppose the linear map f : B B is represented by the matrix (A (κ/d)bc) and denote its eigenvalues as λ 1, = α ± jβ and λ 3 = γ. ince (A (κ/d)bc) is real, there exists a real nonsingular matrix T such that A 1 := T 1 (A (κ/d)bc) T is in real Jordan form [Horn and Johnson, 1985, section 3.4]. In other words, α β 0 A 1 := T 1 (A (κ/d)bc) T = β α 0 (5) 0 0 γ Note that the three columns of the matrix T, which will be denoted as T a, T b, and T c, are linearly independent vectors and hence constitute a basis for R 3 which will be denoted as B. Without loss of generality, from this point onwards we will assume that x = [x 1, x, x 3 ] R 3 is a coordinate vector relative to the basis B. For the system ẋ = A 1 x, an initial condition in the subspace defined by V u := span{t a,t b } will only excite the complex conjugate mode. Hence, a necessary condition for (3) - (4) to have a limit cycle is that the subspace V u is not parallel to the switching surface. Otherwise, any trajectory of the system ẋ = A 1 x starting in would asymptotically approach the subspace V u and spiral away from the origin without reaching neither nor. imilarly, any trajectory with an initial condition in the subspace V s := T c only excites the stable real mode. If V s is not parallel to the switching surface (i.e. T c is not orthogonal to the vector C ), then there exists a point x V s such that the trajectory of (3) - (4) starting in x will converge asymptotically towards the origin. Due to the continuity of solutions in terms of the initial states [Khalil, 000, Theorem 3.5], a simple argument shows that for points in arbitrarily close to V s, the associated switching times become arbitrarily large. This fact, as we will see in the following section, plays a very important role when analysing the stability of the system s limit cycle. 3. MAIN REULT This section starts by providing a sufficient condition which guarantees that the system (3) - (4) possesses an asymptotically stable limit cycle in R 3 \M s. The following subsections show how to overcome certain technicalities which might arise when one verifies the aforementioned condition. 3.1 tability analysis of the limit cycle A sufficient condition for the system (3) - (4) to have a globally asymptotically stable limit cycle in R 3 \ M s is given by an adaptation of the results presented in [Gonçalves, 000, ection 7.3] to the class of systems considered herein. These results are based on the discovery that impact maps can be represented as linear transformations parametrised by the switching times. The stability of the limit cycle is then ensured by proving that the impact maps are contractive around the points in + and which belong to the limit cycle. These points are denoted as x 0 and x 1, respectively. There are three important aspects to consider during the aforementioned adaptation. The first one concerns the selection of the points x 0 and x 1 such that [Gonçalves, 000, Proposition 7.1] holds. Given that A is Hurwitz, let P A > 0 satisfy P A A + A P A < 0 and define the point x 0 as follows: x 0 := (d CA 1 Bκ) P 1 A C CP 1 + A 1 Bκ A C If V s, let x 1 := x; otherwise, x 1 := x 1. Our choice of x 0 guarantees that Cx 0(t) d for all t > 0, as shown in [Gonçalves, 000, ection 7.]. In a similar fashion, if V s, then Cx 1(t) d and Cx 1(t) d for all t > 0. Otherwise, the functions H a (t) and H b (t) defined in [Gonçalves, 000, Proposition 7.1] can be defined via continuation for all t > 0 such that Cx 1(t) = d. In the second place, the Lyapunov function V 1 (V ) introduced in [Gonçalves, 000, ection 7.3] has to be centred around the point in + ( ) which belongs to the limit cycle. To do so, let P i R be symmetric, positive definite matrices, g i R, δ R and α i R such that V i (δ) = δ P i δ δ g i + α i (6) for i = 1,. Let Π C (C stands for the orthogonal complement of C). Defining d 0 = Π x 0 Π x 0 (d 1 = Π x 1 Π x 1) and fixing g 1 = P 1 d 0 (g = P d 1 ) and α 1 = d 0 P 1 d 0 (α = d 1 P d 1 ) centres the ellipse V 1 (δ) = ς (V (δ) = ς) around x 0 (x 1). In the third place, the conditions presented in [Gonçalves, 000, ection 6.7.] must also be modified to incorporate our particular choice of g i and α i for i = 1,. This is done using standard algebra and is left to the reader. The numerical verification of the conditions sketched herein is straightforward when lower and upper bounds on the set of expected switching times can be found. This is done by obtaining a finite sequence t 0 < t 1 < < t k
3 whose elements belong to the set of expected switching times and verifying the conditions presented in [Gonçalves, 000, ection 7.3] on t = {t i }, i = 0,1,...,k. For large enough k, it can be shown that such conditions are also satisfied for all the set of expected switching times, as the latter reference explains. Nevertheless, it was seen in ection. that one cannot always find an upper bound for all the sets of expected switching times associated with the class of systems under consideration. A methodology to address this situation is presented next. 3. Verifying the conditions for stability of the limit cycle. From this point onwards, it will be assumed that the set V s is not empty. Define T a (T b ) as the set of expected switching times for the impact map a (b) which takes points from ( ) and maps them in + ( + ). As discussed in section, the sets T a and T b are not bounded, so a numerical verification of the conditions mentioned in ection 3.1 is futile. A methodology to overcome this problem is presented herein. 1.- Enforcing boundedness of T a and T b. We start by characterising a neighbourhood E ε of points in centred around V s. We also provide an estimation for a bounded set I l that any trajectory of (3) - (4) will eventually enter. If contraction of the impact maps around the limit cycle is guaranteed for points in this bounded set, one can simply wait until the considered trajectory of the system gets there; from that moment onwards, the trajectory will converge to the limit cycle. Using geometrical arguments, we provide an upper bound for the switching time of points in which belong to I l but do not belong to E ε. The impact maps for all such points can then shown to be contractive using the conditions mentioned in ection 3.1. Verifying the contraction of the impact maps for points which belong to E ε is deferred to the following subsection. tart by fixing ε > 0. The set E ε := { x x 1 + x ε } defines an ε-neighbourhood of the subspace V s in the switching surface. For reasons that will become transparent later, the fixed value of ε must be such that E ε V u = (see Theorem 3), E ε {x CA 1 x = 0} = (see Lemma 4), and x 1 / E ε (see the remark after Theorem 7). Now let x be the point which belongs to both V s and. This point fulfills C x = d and x = [0, 0, x 3 ] for a given x 3 > 0. The following lemma shows how to coalesce the definitions of E ε and x. Lemma 1. Let 0 i j be an i j matrix of zeros and let I k be the k k identity matrix. An equivalent description for the set E ε is given by E ε := { x 1 + Πδ 1 δ1 P E δ 1 1 } where P E := 1 I 0 1 ε Π Π is positive definite. Proof: The proof follows from straightforward algebra, the parametrisation of points in as x 1 + Πδ 1 where x 1 = x, and the fact that the coordinates x 1 and x relative to the basis B are both zero for x. The positive definiteness of P E is a direct consequence of C T c 0. ince the matrix A is Hurwitz and u κ is a bounded input, there is a bounded set I l that any trajectory of (3) - (4) will eventually enter, as the following proposition demonstrates. Proposition. Recall the definition of Π as a matrix belonging to the orthogonal complement of C. Let Π 1 Eε Tc Fig.. The value of ε > 0 defines the set E ε. and Π be the columns of the matrix Π. Let m be an integer and pick ψ i for i {1,...,m} such that 0 ψ 1 <... < ψ m < π. Define the row vector F i := cos(ψ i )Π 1 + sin(ψ i )Π. For the system (3) - (4), F 1 x κ F 1 e At B L1 I l := x. (7) F m x κ F m e At B L1 defines a bounded set that any trajectory will eventually enter. Proof: ee [Gonçalves, 000, Proposition 7.]. The definitions introduced below will be needed to show that the switching times for all points in which belong to I l but do not belong to E ε are bounded. Let x 3p be a fixed nonzero real number between 0 and x 3. Consider a plane P parallel to the subspace V u as follows: P := { x } [0, 0, 1]x = x3p For > 0, define a circle in V u relative to the basis B by C 1 ( ) = { x V u x 1 + x = }. A circle in P relative to the basis B can be defined analogously by C ( ) = { x P x 1 + x = }. Let 1 ( ; ) be the minimum value of such that C 1 ( ) (C ( ) ; C ( ) ) (see Figure 3). Now define x 3max as the maximum distance between a point in I l and the subspace V u. P ) C ( C ( ) Fig. 3. Planes V u and P; circles C ( ) and C ( ) As the following theorem demonstrates, defining the sets P and E ε permits the computation of an upper bound for the switching time of all points in the switching surface which belong to the set I l but do not belong to the neighbourhood E ε. Theorem 3. Fix x 3p 0 between 0 and x 3. Pick a value of ε such that ε > 0. Define the time instants t 1 := (ln 1 lnε) /α, t := π/β, t 3 := (lnx 3p lnx 3max ) /γ, t := max(t 1,t 3 ) + t, t 1 := ( ln lnε ) /α, t := π/β, t 3 := (lnx 3p lnx 3max ) /γ, t := max(t 1,t 3) + t for the system ẋ = A 1 x, x(0) ( I inv ) \ E ε (8) Then there exists at least one time instant t s ( t s ) which fulfills x(t s ) (x( t s ) ) and 0 < t s t (0 < t s t). Ta Tc Tb Ta Vu Tb Vu
4 Proof: Define r := x 1 + x and θ := arctan(x /x 1 ). Hence, ẋ = A 1 x can be expressed via the uncoupled differential equations ṙ = α r, θ = β, and w3 = γ w 3 with initial conditions r(0) ε, θ(0) = θ 0, and w 3 (0) w 3max. By assumption, α > 0, β > 0, and γ < 0, so straightforward calculations show that r ( t ) > r(t 1 ) 1, that w 3p w 3 (t 3 ) > w 3 ( t ) > 0, and that θ ( t ) = θ 0 + β max(t 1,t 3 ) + π. From the second set of inequalities it can be seen that the trajectory of (8) remains between the planes P and V u for all time t t 3. In addition, after the time t max(t 1,t 3 ) has elapsed, the distance between the subspace V s and the trajectory of (8) has grown to be greater than or equal to 1 (this is given by the first set of inequalities). Finally, the expression for θ ( t ) shows that the additional time interval t allows a complete revolution around V s to occur, thus completing the proof. As an aside, observe that the choice of ε guarantees that E ε V u =. The proof for t s is similar and is thus left to the reader. Notice that in this case only half a revolution around V s is needed; this is a direct consequence of the fact that > 1 >. The above theorem permits the conditions mentioned in ection 3.1 to be verified for all points in I l which do not belong to E ε. As the points in E ε cannot be neglected, the next section will show how to deal with them..- Investigation of the set E ε. This subsection shows how to verify that trajectories starting in the set E ε \ V s converge to the limit cycle of system (3) - (4). Although the approach used here involves showing that the impact maps are contracting around the limit cycle, its novelty resides in the fact that it does not recur to the switching times associated with points in E ε \ V s. We start by proving that, under some assumptions, a trajectory of ẋ = A 1 x which starts in a point belonging to E ε \ V s will take at least some time t E,a (t E,b ) to reach the switching surface (). Lemma 4. Define w a (t) := w a (0) := CA1Π CA 1x. Let 1 d a (t) = 1 P E 1 w a(t), d a (t) = CeA 1 t Π d Cx 1 (t) for t > 0 and d a (t) ( d a (t)p E da (t) ) 1/ for all t 0. Assume that the value of ε chosen to characterise the set E ε is such that (w a (0)d a (0) 1)( w a (0)d a (0) 1) > 0 Let t E,a be the smallest value of t T a \ {0} for which the expression (w a (t)d a (t) 1) ( w a (t)d a (t) 1) (9) is less than or equal to zero. Then, all points in E ε \ V s which switch in have an associated switching time greater than or equal to t E,a. It can also be shown that all points in E ε \ V s which switch in have an associated switching time greater than or equal to t E,b by defining w b (t) = t T b and using this expression instead of w a (t) in the abovementioned result. Proof: Define ta ( tb ) as the set of initial conditions x 1a (x 1b ) such that d Cx(t) d on [0,t a ] and Cx(t a ) = d ( d Cx(t) d on [0,t b ] and Cx(t b ) = d). Corollary 4.1 in Gonçalves [000] shows CeA 1 t Π d Cx 1 (t) for that, for all t a in T a \ {0}, the set ta is a subset of the linear manifold ta := {x 1 + Πδ 1 w a (t a )δ 1 = 1}. For t a = 0, the set ta=0 is given by {x 1 + Πδ 1 w a,0 δ 1 = 1} since w a,0 δ 1 = lim t 0 w a (t)δ 1. It is then obvious that if the set E ε has no points in common with ta for all t a [0,t E,a ), then the switching time for all points in E ε must be equal to or greater than t E,a. As a shorthand, we will be using w at for w a (t a ) and d at for d a (t a ). Define the sets E,t := { δ 1 R δ at P E δ 1 = 1 } w,t := { δ 1 R w at δ 1 = 1 } traightforward manipulations show that δat P E = cw at for c = ( ) 1/, wat P E w at so the sets w,t, E,t and E,t represent lines in R parallel to each other. Furthermore, E,t and E,t are tangent to the ellipse δ 1 P E δ 1 = 1 at the points δ 1 = δ at and δ 1 = δ at. It then follows that ta E ε if and only if (w at d at 1)( w at d at 1) > 0. This equivalence is explained by the fact that (w at d at 1) and (w at ( d at ) 1) have the same sign if and only if w,t is not located between the lines E,t and E,t. This means that no δ 1 can satisfy w at δ 1 = 1 and δ1 P E δ 1 1 simultaneously, as Figure 4 shows. In this case, the sets E ε and ta have no points in common, so our claim follows. The proof for t E,b follows the same line of reasoning and is thus omitted. Πδat Πδat x 1 Eε Fig. 4. ta E ε {. The dashed lines } represent the sets x 1 + Πδ 1 δ 1 E,t and { } x 1 + Πδ 1 δ1 E,t. Using time t E,a (t E,b ), it can be determined how close to the subspace V u the trajectory has to be when it reaches (). Lemma 5. Let d ( \ V s ) ( d ( \ V s )) be the set of points that will eventually switch in (). Define t E,a and t E,b as in Lemma 4. Let x 3E+ be the maximum distance between a point in E ε and the subspace V u. Define x 3Ea := x 3E+ exp(γ t E,a ) and x 3Eb := x 3E+ exp(γ t E,b ). Then any trajectory of ẋ = A 1 x starting in d E ε ( d E ε ) will not be away from the subspace V u by a distance greater than x 3Ea (x 3Eb ) when it switches at (). Proof: The proof follows from the definition of t E,a, t E,b and x 3E+ by recalling that the stable behaviour of ẋ = A 1 x is uncoupled from the unstable one. An additional geometrical argument then characterises subsets of and that this trajectory will reach. Lemma 6. Assume that the value of ε used to characterise the set E ε complies with the conditions given in Theorem 3 and Lemma 4. Define := d + c ( ) 3 x 3Ea πα Ea ( c1 ) + ( c ), q := exp (10) Ea β where [ c 1 c c 3 ] := C, and x 3Ea is defined as in Lemma 5. Let the set η,a be defined by Πw at ta
5 { CA1 x C x = d, x 1 + x q, 0 x 3 x 3Ea } and denote its infimum and supremum values as η a1 and η a, respectively. If η a1 is negative, then re-define it as η a1 := 0. Then all trajectories of ẋ = A 1 x which start in the set E ε \V s and eventually switch in will do so in the set Ea defined by Ea := {x 0 x 3 x 3Ea, η a1 CA 1 x η a } Now define q b via the substitution of x 3Ea by x 3Eb in (10). Define the set η,b by { CA1 x C x = d, x 1 + x q b, 0 x 3 x 3Eb } and denote its infimum and supremum values by η b1 and η b, respectively. If η b is positive, then re-define it as η b := 0. Then all trajectories of ẋ = A 1 x which start in the set E ε \ V s and eventually switch in will do so in the set Eb defined by {x 0 x 3 x 3Eb, η b1 CA 1 x η b } Proof: We begin by noticing that x 3Ea can be used, as in ection 3..1, to define a plane P Ea which is parallel to the subspace V u. The value of associated with P Ea is denoted as and is greater than ε due to the assumptions placed on the latter. Now recall the definitions Ea r := x 1 + x, θ := arctan(x /x 1 ), and the fact that the solutions at time t = π/β of the differential equations ṙ = α r; θ = β with initial conditions r(0) = Ea ; θ(0) = θ 0 are equal to q; θ 0 + π. This shows that once any trajectory of ẋ = A 1 x starting in d E ε has reached the set { x R } 3 x 1 + x = Ea, it must undergo another complete revolution around V s before reaching { x R 3 x 1 + x q }. Given the geometrical interpretation of Ea, this cannot be done without reaching the switching surfaces and first. To prove that a point in d E ε cannot eventually switch outside Ea, we proceed by contradiction. uppose that a point x d E ε switches outside Ea. Notice that switching in the set {x x 3 < 0} is impossible due to the dynamics of the system. Furthermore, Lemma 5 shows that switching in the set {x x 3Ea < x 3 } is impossible too. Hence, the point in d E ε is restricted to switch in the set := 1 where 1 := {x 0 x 3 x 3Ea, η 1 > CA 1 x} := {x 0 x 3 x 3Ea, η < CA 1 x} ince the definitions of η 1 and η ensure that the set { x R 3 Cx = d, 0 x3 x 3Ea, x 1 + x q } is contained in the set Ea, it is readily seen that x 1 + x > q for all points x. This contradicts the information presented in the previous paragraph. In addition, all trajectories of ẋ = A 1 x starting in are restricted to switch either in the set + or in the set +. Hence, our claim follows. The proof for Eb follows the same line of reasoning and is thus omitted. We now present a condition which guarantees that the mapping of points in E ε \V s to the aforementioned subsets of and is contracting around the limit cycle. Define E := Ea and E := Eb. We furthermore consider the following notation: if is a set, then define the set x s as x s := {x x s x } To verify that trajectories starting in the set E ε \ V s converge to the limit cycle of system (3) - (4), it is sufficient to show that the impact maps which take points from E ε \ V s and map them in and are contracting. A sufficient condition which guarantees such contraction is { } { } max V1 (δ a ) s.t. Πδ a E x < min V (δ 1 ) 0 s.t. Πδ 1 E ε x 1 { } { } max V1 (δ b ) s.t. Πδ b E x < min V (δ 1 ) 0 s.t. Πδ 1 E ε x 1 (11) We now present a set of linear matrix inequalities which, if fulfilled, guarantee that conditions (11) hold. Theorem 7. Define Γ := [0, 0, 1] and let M be a 3 3 matrix with rows given by Γ, C, and CA 1, respectively. Let ϑ a,1 := [ 0, d, η a1 ] ϑ a, := [ x 3Ea, d, η a1 ] ϑ a,3 := [ 0, d, η a ] ϑ a,4 := [ x 3Ea, d, η a ] ϑ b,1 := [ 0, d, η b1 ] ϑ b, := [ x 3Eb, d, η b1 ] ϑ b,3 := [ 0, d, η b ] ϑ b,4 := [ x 3Eb, d, η b ] where w 3Ea, w 3Eb, η a1, η a, η b1, and η b define the sets E and E. ( Define σ a,i = V { 1 Π ( M 1 ϑ a,i x0}) and σb,i = V { 1 Π M 1 ϑ b,i x0}) for i {1,,3,4}. If there exist nonnegative scalars τ a,i and τ b,i such that the matrices P + τ a,i P E g g α σ a,i τ a,i P + τ b,i P E g g α σ b,i τ b,i are positive definite for all i {1,,3,4}, then conditions (11) hold. Proof: The matrix M can be shown to be non-singular by considering that A 1 ξi 3 together with the fact that the vector Γ / V u. Let i {1,,3,4}. The fact that M 1 exists shows that E ( E ) can be visualised as the set of points in enclosed within the parallelogram with vertices M 1 ϑ a,i (M 1 ϑ b,i ). ince the function V 1 is convex by construction (see section 3.1), maximising its value over the sets E x 0 and E x 0 can be done by computing its value at the points Π {M 1 ϑ a,i x 0} and Π {M 1 ϑ b,i x 0}; the biggest of such values is the desired result. Hence, enforcing conditions (11) is equivalent to asking the inequalities σ a,i < V (δ) = δ P δ δ g + α σ b,i < V (δ) = δ P δ δ g + α to hold for all δ 0 such that δ P E δ 1. It is obvious that if there exist nonnegative scalars τ a,1 and τ b,1 such that δ P δ δ g + α σ a,i τ a,i (1 δ P E δ) > 0 δ P δ δ g + α σ b,i τ b,i (1 δ P E δ) > 0 hold, then the aforementioned inequalities hold too. Realising that one can write this as δ1 P + τ a,i P E g δ1 g > 0 α σ a,i τ a,i δ1 P + τ b,i P E g δ1 g > 0 α σ b,i τ b,i completes the proof. Remark: For Theorem 7 to be useful, the set E ε must be defined in such a way that it does not contain the point x 1. If it did, min V (δ 1 ) s.t. Πδ 1 E ε x 1
6 would be equal to zero and, given that V 1 (δ) is nonnegative for all values of δ, the inequalities in (11) would never be fulfilled. 4. EXAMPLE. The following example is a PWL version of the third order dimensionless Goodwin oscillator model whose parameters are congruent with those presented in tan et al. [007]. The IO LTI block obeys ẋ = [ ] [ 10 x + 0 ] u y = [ ]x while the saturation block is described by { 1/ if y(t) < 1/9 u(t) = 9/ y(t) if y(t) 1/9 1/ if y(t) > 1/9 The feedback system has a locally stable limit cycle that switches four times per cycle with period t ( ). The intersection of the limit cycle with the switching surface occurs at x 0 [0.065, 0.634, ] and x 1 [ , 0.035, ]. The set I l was constructed according to Proposition by selecting m = 18 and ψ i = (m 1)(π/18). The choice of x 3p and ε = yields T a = [0,15] and T b = [0.47,18.43], (refer to section 3..1 and [Gonçalves, 000, Proposition 7.5]). For this system, the Lyapunov functions on are defined by the matrices P 1 ( ) ( and P ) To show the validity of these functions, the system was simulated on the time interval [0, 30] starting from the initial condition x(0) = [0.6, 0.6, 1/9]. The most relevant results of such simulation are shown in table 1; as expected, the value of the Lyapunov functions decreases every time the trajectory reaches a switching surface. t x(t) V i 0 [ 0.6, 0.6, ] V 1 = [ , -0.7, ] V = [ , , ] V 1 = [ 0.647, , ] V = [ 0.193, 0.39, ] V 1 = [ , , ] V = [ , , ] V 1 = [ 0.546, , ] V = [ , 0.833, ] V 1 = [ , -0.05, ] V = [ , , ] V 1 = [ 0.501, , ] V = [ , 0.700, ] V 1 = [ , , ] V = [ , -0.67, ] V 1 = [ , , ] V = Table 1. Results of the system s simulation over the time interval [0, 30] with initial condition x(0) = [0.6, 0.6, 1/9]. The first and second columns show the time instants in which the system s trajectory reaches a switching surface and the points in which this happens, while the third one gives the value of the appropriate Lyapunov function at such points. 5. CONCLUION. This paper presented a methodology to analyse the global stability of the limit cycle of a PWL saturation system which cannot be directly analysed using the methodology presented in tan [005]. It is based on the construction of quadratic Lyapunov functions on the system s switching surfaces and covers the case in which the expected switching times associated with the system s impact maps are arbitrarily large. The first step involves defining a subset of the system s switching surfaces in which the switching times are bounded; the methodology developed by Gonçalves [000] is then applied to the study of such subset. The second step recurs to a geometrical argument in order to guarantee that all the trajectories which do not start in the aforementioned subset converge either to the origin or to the limit cycle. Taken together, both steps achieve the desired objective. Future work involves the extension of the present results to higher-dimensional systems, thus opening the door to the analysis of models other than the third-order Goodwin oscillator presented herein. REFERENCE U. Boscain. tability of planar switched systems: the linear single input case. IAM Journal on Control and Optimization, 41(1):89 11, 00. M.. Branicky. Multiple Lyapunov functions and other analysis tools for switched and hybrid systems. IEEE Transactions on Automatic Control, 43(4):475 48, A. Goldbeter. Biochemical oscillations and cellular rhytms: the molecular bases of periodic and chaotic behaviour. Cambridge University Press, UK, J. M. Gonçalves. Constructive global analysis of hybrid systems. PhD thesis, Massachusetts Institute of Technology, eptember 000. URL /papers/cued_control_893.pdf. J. M. Gonçalves. Regions of stability for limit cycle oscillations in piecewise linear systems. IEEE Transactions on Automatic Control, 50(11): , Nov R. A. Horn and C. R. Johnson. Matrix analysis. Cambridge University Press, UK, first edition, Y. Iwatani and. Hara. tability tests and stabilization for piecewise linear systems based on poles and zeros of subsystems. Automatica, 4(10): , 006. M. Johansson and A. Rantzer. Computation of piecewise quadratic Lyapunov functions for hybrid systems. IEEE Transactions on Automatic Control, 43(4): , April H. Khalil. Nonlinear systems. Prentice Hall, UA, third edition, 000. M. Margaliot. tability analysis of switched systems using variational principles: an introduction. Automatica, 4 (1): , Dec E. Mosekilde. Topics in nonlinear dynamics; applications to physics, biology, and economic systems. World cientific Press, ingapore, Petterson and B. Lennartson. Hybrid system stability and robustness verification using linear matrix inequalities. International Journal of Control, 75: , 00.. Prajna and A. Papachristodoulou. Analysis of switched and hybrid systems - beyond piecewise quadratic methods. In American Control Conference (ACC), Denver, CO, UA, 003. G.-B. tan. Global analysis and synthesis of oscillations: a dissipativity approach. PhD thesis, Université de Liège, 005. G.-B. tan, A. Hamadeh, R. epulchre, and J. M. Gonçalves. Output synchronization in networks of cyclic biochemical oscillators. In American Control Conference (ACC), Times quare, New York City, UA, July H. trogatz. ynch: the emerging science of spontaneous order. Hyperion, 003. X. Xu and P. J. Antsaklis. tabilization of second-order LTI switched systems. In IEEE Conference on Decision and Control (CDC), pages , Phoenix, Arizona, UA, 1999.
GLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS
GLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS Jorge M. Gonçalves, Alexandre Megretski y, Munther A. Dahleh y California Institute of Technology
More informationGlobal Analysis of Piecewise Linear Systems Using Impact Maps and Quadratic Surface Lyapunov Functions
Global Analysis of Piecewise Linear Systems Using Impact Maps and Quadratic Surface Lyapunov Functions Jorge M. Gonçalves, Alexandre Megretski, Munther A. Dahleh Department of EECS, Room 35-41 MIT, Cambridge,
More informationGlobal Analysis of Piecewise Linear Systems Using Impact Maps and Surface Lyapunov Functions
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 48, NO 12, DECEMBER 2003 2089 Global Analysis of Piecewise Linear Systems Using Impact Maps and Surface Lyapunov Functions Jorge M Gonçalves, Alexandre Megretski,
More informationThe servo problem for piecewise linear systems
The servo problem for piecewise linear systems Stefan Solyom and Anders Rantzer Department of Automatic Control Lund Institute of Technology Box 8, S-22 Lund Sweden {stefan rantzer}@control.lth.se Abstract
More informationTrajectory Tracking Control of Bimodal Piecewise Affine Systems
25 American Control Conference June 8-1, 25. Portland, OR, USA ThB17.4 Trajectory Tracking Control of Bimodal Piecewise Affine Systems Kazunori Sakurama, Toshiharu Sugie and Kazushi Nakano Abstract This
More informationStabilization of Second-Order LTI Switched Systems
F'roceedings of the 38' Conference on Decision & Control Phoenix, Arizona USA - December 1999 Stabilization of Second-Order LTI Switched Systems Xuping Xu' and Panos J. Antsaklis Department of Electrical
More informationOn Linear Copositive Lyapunov Functions and the Stability of Switched Positive Linear Systems
1 On Linear Copositive Lyapunov Functions and the Stability of Switched Positive Linear Systems O. Mason and R. Shorten Abstract We consider the problem of common linear copositive function existence for
More informationPassivity-based Stabilization of Non-Compact Sets
Passivity-based Stabilization of Non-Compact Sets Mohamed I. El-Hawwary and Manfredi Maggiore Abstract We investigate the stabilization of closed sets for passive nonlinear systems which are contained
More informationSTABILITY OF PLANAR NONLINEAR SWITCHED SYSTEMS
LABORATOIRE INORMATIQUE, SINAUX ET SYSTÈMES DE SOPHIA ANTIPOLIS UMR 6070 STABILITY O PLANAR NONLINEAR SWITCHED SYSTEMS Ugo Boscain, régoire Charlot Projet TOpModel Rapport de recherche ISRN I3S/RR 2004-07
More informationHybrid Systems - Lecture n. 3 Lyapunov stability
OUTLINE Focus: stability of equilibrium point Hybrid Systems - Lecture n. 3 Lyapunov stability Maria Prandini DEI - Politecnico di Milano E-mail: prandini@elet.polimi.it continuous systems decribed by
More informationx= Ax+Bd S + Cx>d d<cx<d x* 1 x Px=k x= Ax Bd
Quadratic urface Lyapunov Functions in Global tability Analysis of aturation ystems Jorge M. Gonοcalves Λ Department of EEC, Room 35-41 MIT, Cambridge, MA jmg@mit.edu eptember 3, 2 Abstract This paper
More informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More informationStability Analysis and Synthesis for Scalar Linear Systems With a Quantized Feedback
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 48, NO 9, SEPTEMBER 2003 1569 Stability Analysis and Synthesis for Scalar Linear Systems With a Quantized Feedback Fabio Fagnani and Sandro Zampieri Abstract
More informationGramians based model reduction for hybrid switched systems
Gramians based model reduction for hybrid switched systems Y. Chahlaoui Younes.Chahlaoui@manchester.ac.uk Centre for Interdisciplinary Computational and Dynamical Analysis (CICADA) School of Mathematics
More informationOn Dwell Time Minimization for Switched Delay Systems: Free-Weighting Matrices Method
On Dwell Time Minimization for Switched Delay Systems: Free-Weighting Matrices Method Ahmet Taha Koru Akın Delibaşı and Hitay Özbay Abstract In this paper we present a quasi-convex minimization method
More informationLeast Squares Based Self-Tuning Control Systems: Supplementary Notes
Least Squares Based Self-Tuning Control Systems: Supplementary Notes S. Garatti Dip. di Elettronica ed Informazione Politecnico di Milano, piazza L. da Vinci 32, 2133, Milan, Italy. Email: simone.garatti@polimi.it
More informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More informationThe ϵ-capacity of a gain matrix and tolerable disturbances: Discrete-time perturbed linear systems
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 11, Issue 3 Ver. IV (May - Jun. 2015), PP 52-62 www.iosrjournals.org The ϵ-capacity of a gain matrix and tolerable disturbances:
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationQuadratic and Copositive Lyapunov Functions and the Stability of Positive Switched Linear Systems
Proceedings of the 2007 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 2007 WeA20.1 Quadratic and Copositive Lyapunov Functions and the Stability of
More informationZeros and zero dynamics
CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)
More informationEE363 homework 7 solutions
EE363 Prof. S. Boyd EE363 homework 7 solutions 1. Gain margin for a linear quadratic regulator. Let K be the optimal state feedback gain for the LQR problem with system ẋ = Ax + Bu, state cost matrix Q,
More informationEvent-Triggered Decentralized Dynamic Output Feedback Control for LTI Systems
Event-Triggered Decentralized Dynamic Output Feedback Control for LTI Systems Pavankumar Tallapragada Nikhil Chopra Department of Mechanical Engineering, University of Maryland, College Park, 2742 MD,
More informationDisturbance Attenuation Properties for Discrete-Time Uncertain Switched Linear Systems
Disturbance Attenuation Properties for Discrete-Time Uncertain Switched Linear Systems Hai Lin Department of Electrical Engineering University of Notre Dame Notre Dame, IN 46556, USA Panos J. Antsaklis
More informationPractical Stabilization of Integrator Switched Systems
Practical Stabilization of Integrator Switched Systems Xuping Xu and Panos J. Antsaklis 2 Department of Electrical and Computer Engineering Penn State Erie, Erie, PA 6563, USA. E-mail: Xuping-Xu@psu.edu
More informationLimit Cycles in High-Resolution Quantized Feedback Systems
Limit Cycles in High-Resolution Quantized Feedback Systems Li Hong Idris Lim School of Engineering University of Glasgow Glasgow, United Kingdom LiHonIdris.Lim@glasgow.ac.uk Ai Poh Loh Department of Electrical
More informationNetwork Analysis of Biochemical Reactions in Complex Environments
1 Introduction 1 Network Analysis of Biochemical Reactions in Complex Environments Elias August 1 and Mauricio Barahona, Department of Bioengineering, Imperial College London, South Kensington Campus,
More informationHybrid Systems Course Lyapunov stability
Hybrid Systems Course Lyapunov stability OUTLINE Focus: stability of an equilibrium point continuous systems decribed by ordinary differential equations (brief review) hybrid automata OUTLINE Focus: stability
More informationObserver-based quantized output feedback control of nonlinear systems
Proceedings of the 17th World Congress The International Federation of Automatic Control Observer-based quantized output feedback control of nonlinear systems Daniel Liberzon Coordinated Science Laboratory,
More informationMinimum-Phase Property of Nonlinear Systems in Terms of a Dissipation Inequality
Minimum-Phase Property of Nonlinear Systems in Terms of a Dissipation Inequality Christian Ebenbauer Institute for Systems Theory in Engineering, University of Stuttgart, 70550 Stuttgart, Germany ce@ist.uni-stuttgart.de
More informationOn the simultaneous diagonal stability of a pair of positive linear systems
On the simultaneous diagonal stability of a pair of positive linear systems Oliver Mason Hamilton Institute NUI Maynooth Ireland Robert Shorten Hamilton Institute NUI Maynooth Ireland Abstract In this
More informationResearch Article Convex Polyhedron Method to Stability of Continuous Systems with Two Additive Time-Varying Delay Components
Applied Mathematics Volume 202, Article ID 689820, 3 pages doi:0.55/202/689820 Research Article Convex Polyhedron Method to Stability of Continuous Systems with Two Additive Time-Varying Delay Components
More informationAnalysis and design of switched normal systems
Nonlinear Analysis 65 (2006) 2248 2259 www.elsevier.com/locate/na Analysis and design of switched normal systems Guisheng Zhai a,, Xuping Xu b, Hai Lin c, Anthony N. Michel c a Department of Mechanical
More informationEN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015
EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions
More information2.10 Saddles, Nodes, Foci and Centers
2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationConverse Lyapunov theorem and Input-to-State Stability
Converse Lyapunov theorem and Input-to-State Stability April 6, 2014 1 Converse Lyapunov theorem In the previous lecture, we have discussed few examples of nonlinear control systems and stability concepts
More informationNonlinear Control Design for Linear Differential Inclusions via Convex Hull Quadratic Lyapunov Functions
Nonlinear Control Design for Linear Differential Inclusions via Convex Hull Quadratic Lyapunov Functions Tingshu Hu Abstract This paper presents a nonlinear control design method for robust stabilization
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More information550 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 46, NO. 4, APRIL Global Stability of Relay Feedback Systems
550 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 46, NO. 4, APRIL 2001 Global Stability of Relay Feedback Systems Jorge M. Gonçalves, Alexre Megretski, Munther A. Dahleh Abstract For a large class of relay
More informationObserver design for a general class of triangular systems
1st International Symposium on Mathematical Theory of Networks and Systems July 7-11, 014. Observer design for a general class of triangular systems Dimitris Boskos 1 John Tsinias Abstract The paper deals
More informationEE222 - Spring 16 - Lecture 2 Notes 1
EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued
More informationA Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems
53rd IEEE Conference on Decision and Control December 15-17, 2014. Los Angeles, California, USA A Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems Seyed Hossein Mousavi 1,
More informationPARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT
PARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT Hans Norlander Systems and Control, Department of Information Technology Uppsala University P O Box 337 SE 75105 UPPSALA, Sweden HansNorlander@ituuse
More informationIMPROVED MPC DESIGN BASED ON SATURATING CONTROL LAWS
IMPROVED MPC DESIGN BASED ON SATURATING CONTROL LAWS D. Limon, J.M. Gomes da Silva Jr., T. Alamo and E.F. Camacho Dpto. de Ingenieria de Sistemas y Automática. Universidad de Sevilla Camino de los Descubrimientos
More informationECE504: Lecture 8. D. Richard Brown III. Worcester Polytechnic Institute. 28-Oct-2008
ECE504: Lecture 8 D. Richard Brown III Worcester Polytechnic Institute 28-Oct-2008 Worcester Polytechnic Institute D. Richard Brown III 28-Oct-2008 1 / 30 Lecture 8 Major Topics ECE504: Lecture 8 We are
More informationState and Parameter Estimation Based on Filtered Transformation for a Class of Second-Order Systems
State and Parameter Estimation Based on Filtered Transformation for a Class of Second-Order Systems Mehdi Tavan, Kamel Sabahi, and Saeid Hoseinzadeh Abstract This paper addresses the problem of state and
More informationNecessary and Sufficient Conditions for Reachability on a Simplex
Necessary and Sufficient Conditions for Reachability on a Simplex Bartek Roszak a, Mireille E. Broucke a a Edward S. Rogers Sr. Department of Electrical and Computer Engineering, University of Toronto,
More informationStability Analysis of a Proportional with Intermittent Integral Control System
American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July, ThB4. Stability Analysis of a Proportional with Intermittent Integral Control System Jin Lu and Lyndon J. Brown Abstract
More informationCopyright (c) 2006 Warren Weckesser
2.2. PLANAR LINEAR SYSTEMS 3 2.2. Planar Linear Systems We consider the linear system of two first order differential equations or equivalently, = ax + by (2.7) dy = cx + dy [ d x x = A x, where x =, and
More informationStability lectures. Stability of Linear Systems. Stability of Linear Systems. Stability of Continuous Systems. EECE 571M/491M, Spring 2008 Lecture 5
EECE 571M/491M, Spring 2008 Lecture 5 Stability of Continuous Systems http://courses.ece.ubc.ca/491m moishi@ece.ubc.ca Dr. Meeko Oishi Electrical and Computer Engineering University of British Columbia,
More informationStability of Hybrid Control Systems Based on Time-State Control Forms
Stability of Hybrid Control Systems Based on Time-State Control Forms Yoshikatsu HOSHI, Mitsuji SAMPEI, Shigeki NAKAURA Department of Mechanical and Control Engineering Tokyo Institute of Technology 2
More informationOn the Kalman-Yacubovich-Popov lemma and common Lyapunov solutions for matrices with regular inertia
On the Kalman-Yacubovich-Popov lemma and common Lyapunov solutions for matrices with regular inertia Oliver Mason, Robert Shorten and Selim Solmaz Abstract In this paper we extend the classical Lefschetz
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More information4 Second-Order Systems
4 Second-Order Systems Second-order autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization
More informationOn common quadratic Lyapunov functions for stable discrete-time LTI systems
On common quadratic Lyapunov functions for stable discrete-time LTI systems Oliver Mason, Hamilton Institute, NUI Maynooth, Maynooth, Co. Kildare, Ireland. (oliver.mason@may.ie) Robert Shorten 1, Hamilton
More information1 The Observability Canonical Form
NONLINEAR OBSERVERS AND SEPARATION PRINCIPLE 1 The Observability Canonical Form In this Chapter we discuss the design of observers for nonlinear systems modelled by equations of the form ẋ = f(x, u) (1)
More informationA Characterization of the Hurwitz Stability of Metzler Matrices
29 American Control Conference Hyatt Regency Riverfront, St Louis, MO, USA June -2, 29 WeC52 A Characterization of the Hurwitz Stability of Metzler Matrices Kumpati S Narendra and Robert Shorten 2 Abstract
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More informationEnergy-based Swing-up of the Acrobot and Time-optimal Motion
Energy-based Swing-up of the Acrobot and Time-optimal Motion Ravi N. Banavar Systems and Control Engineering Indian Institute of Technology, Bombay Mumbai-476, India Email: banavar@ee.iitb.ac.in Telephone:(91)-(22)
More informationAQUANTIZER is a device that converts a real-valued
830 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 57, NO 4, APRIL 2012 Input to State Stabilizing Controller for Systems With Coarse Quantization Yoav Sharon, Member, IEEE, Daniel Liberzon, Senior Member,
More informationL 2 -induced Gains of Switched Systems and Classes of Switching Signals
L 2 -induced Gains of Switched Systems and Classes of Switching Signals Kenji Hirata and João P. Hespanha Abstract This paper addresses the L 2-induced gain analysis for switched linear systems. We exploit
More informationIMPULSIVE CONTROL OF DISCRETE-TIME NETWORKED SYSTEMS WITH COMMUNICATION DELAYS. Shumei Mu, Tianguang Chu, and Long Wang
IMPULSIVE CONTROL OF DISCRETE-TIME NETWORKED SYSTEMS WITH COMMUNICATION DELAYS Shumei Mu Tianguang Chu and Long Wang Intelligent Control Laboratory Center for Systems and Control Department of Mechanics
More informationAn asymptotic ratio characterization of input-to-state stability
1 An asymptotic ratio characterization of input-to-state stability Daniel Liberzon and Hyungbo Shim Abstract For continuous-time nonlinear systems with inputs, we introduce the notion of an asymptotic
More informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
More informationProblem List MATH 5173 Spring, 2014
Problem List MATH 5173 Spring, 2014 The notation p/n means the problem with number n on page p of Perko. 1. 5/3 [Due Wednesday, January 15] 2. 6/5 and describe the relationship of the phase portraits [Due
More informationA LaSalle version of Matrosov theorem
5th IEEE Conference on Decision Control European Control Conference (CDC-ECC) Orlo, FL, USA, December -5, A LaSalle version of Matrosov theorem Alessro Astolfi Laurent Praly Abstract A weak version of
More informationControl for Coordination of Linear Systems
Control for Coordination of Linear Systems André C.M. Ran Department of Mathematics, Vrije Universiteit De Boelelaan 181a, 181 HV Amsterdam, The Netherlands Jan H. van Schuppen CWI, P.O. Box 9479, 19 GB
More informationSingular perturbation analysis of an additive increase multiplicative decrease control algorithm under time-varying buffering delays.
Singular perturbation analysis of an additive increase multiplicative decrease control algorithm under time-varying buffering delays. V. Guffens 1 and G. Bastin 2 Intelligent Systems and Networks Research
More informationDelay-dependent Stability Analysis for Markovian Jump Systems with Interval Time-varying-delays
International Journal of Automation and Computing 7(2), May 2010, 224-229 DOI: 10.1007/s11633-010-0224-2 Delay-dependent Stability Analysis for Markovian Jump Systems with Interval Time-varying-delays
More informationarzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.1 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS STABILITY ANALYSIS Didier HENRION www.laas.fr/ henrion henrion@laas.fr Denis ARZELIER www.laas.fr/
More informationHYBRID LIMIT CYCLES AND HYBRID POINCARÉ-BENDIXSON. Slobodan N. Simić
HYBRID LIMIT CYCLES AND HYBRID POINCARÉ-BENDIXSON Slobodan N. Simić Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, CA 94720-1774, U.S.A. Email: simic@eecs.berkeley.edu
More informationOn the Stabilization of Neutrally Stable Linear Discrete Time Systems
TWCCC Texas Wisconsin California Control Consortium Technical report number 2017 01 On the Stabilization of Neutrally Stable Linear Discrete Time Systems Travis J. Arnold and James B. Rawlings Department
More information3118 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 57, NO. 12, DECEMBER 2012
3118 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 57, NO 12, DECEMBER 2012 ASimultaneousBalancedTruncationApproachto Model Reduction of Switched Linear Systems Nima Monshizadeh, Harry L Trentelman, SeniorMember,IEEE,
More informationOutput Input Stability and Minimum-Phase Nonlinear Systems
422 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 47, NO. 3, MARCH 2002 Output Input Stability and Minimum-Phase Nonlinear Systems Daniel Liberzon, Member, IEEE, A. Stephen Morse, Fellow, IEEE, and Eduardo
More informationFINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES. Danlei Chu, Tongwen Chen, Horacio J. Marquez
FINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES Danlei Chu Tongwen Chen Horacio J Marquez Department of Electrical and Computer Engineering University of Alberta Edmonton
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More information2nd Symposium on System, Structure and Control, Oaxaca, 2004
263 2nd Symposium on System, Structure and Control, Oaxaca, 2004 A PROJECTIVE ALGORITHM FOR STATIC OUTPUT FEEDBACK STABILIZATION Kaiyang Yang, Robert Orsi and John B. Moore Department of Systems Engineering,
More informationDOMAIN OF ATTRACTION: ESTIMATES FOR NON-POLYNOMIAL SYSTEMS VIA LMIS. Graziano Chesi
DOMAIN OF ATTRACTION: ESTIMATES FOR NON-POLYNOMIAL SYSTEMS VIA LMIS Graziano Chesi Dipartimento di Ingegneria dell Informazione Università di Siena Email: chesi@dii.unisi.it Abstract: Estimating the Domain
More informationTangent spaces, normals and extrema
Chapter 3 Tangent spaces, normals and extrema If S is a surface in 3-space, with a point a S where S looks smooth, i.e., without any fold or cusp or self-crossing, we can intuitively define the tangent
More informationConverse Results on Existence of Sum of Squares Lyapunov Functions
2011 50th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC) Orlando, FL, USA, December 12-15, 2011 Converse Results on Existence of Sum of Squares Lyapunov Functions Amir
More informationOn Piecewise Quadratic Control-Lyapunov Functions for Switched Linear Systems
On Piecewise Quadratic Control-Lyapunov Functions for Switched Linear Systems Wei Zhang, Alessandro Abate, Michael P. Vitus and Jianghai Hu Abstract In this paper, we prove that a discrete-time switched
More informationNonlinear Tracking Control of Underactuated Surface Vessel
American Control Conference June -. Portland OR USA FrB. Nonlinear Tracking Control of Underactuated Surface Vessel Wenjie Dong and Yi Guo Abstract We consider in this paper the tracking control problem
More informationENERGY DECAY ESTIMATES FOR LIENARD S EQUATION WITH QUADRATIC VISCOUS FEEDBACK
Electronic Journal of Differential Equations, Vol. 00(00, No. 70, pp. 1 1. ISSN: 107-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp ENERGY DECAY ESTIMATES
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More informationIn these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.
1 2 Linear Systems In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation 21 Matrix ODEs Let and is a scalar A linear function satisfies Linear superposition ) Linear
More informationLinear Algebra. Paul Yiu. 6D: 2-planes in R 4. Department of Mathematics Florida Atlantic University. Fall 2011
Linear Algebra Paul Yiu Department of Mathematics Florida Atlantic University Fall 2011 6D: 2-planes in R 4 The angle between a vector and a plane The angle between a vector v R n and a subspace V is the
More informationReview of Controllability Results of Dynamical System
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 13, Issue 4 Ver. II (Jul. Aug. 2017), PP 01-05 www.iosrjournals.org Review of Controllability Results of Dynamical System
More informationRow and Column Representatives in Qualitative Analysis of Arbitrary Switching Positive Systems
ROMANIAN JOURNAL OF INFORMATION SCIENCE AND TECHNOLOGY Volume 19, Numbers 1-2, 2016, 127 136 Row and Column Representatives in Qualitative Analysis of Arbitrary Switching Positive Systems Octavian PASTRAVANU,
More informationThe norms can also be characterized in terms of Riccati inequalities.
9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements
More informationStability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay
1 Stability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay Wassim M. Haddad and VijaySekhar Chellaboina School of Aerospace Engineering, Georgia Institute of Technology, Atlanta,
More informationChapter III. Stability of Linear Systems
1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,
More informationA conjecture on sustained oscillations for a closed-loop heat equation
A conjecture on sustained oscillations for a closed-loop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping
More informationResearch Article Stabilization Analysis and Synthesis of Discrete-Time Descriptor Markov Jump Systems with Partially Unknown Transition Probabilities
Research Journal of Applied Sciences, Engineering and Technology 7(4): 728-734, 214 DOI:1.1926/rjaset.7.39 ISSN: 24-7459; e-issn: 24-7467 214 Maxwell Scientific Publication Corp. Submitted: February 25,
More informationRobust Connectivity Analysis for Multi-Agent Systems
Robust Connectivity Analysis for Multi-Agent Systems Dimitris Boskos and Dimos V. Dimarogonas Abstract In this paper we provide a decentralized robust control approach, which guarantees that connectivity
More informationOperators with numerical range in a closed halfplane
Operators with numerical range in a closed halfplane Wai-Shun Cheung 1 Department of Mathematics, University of Hong Kong, Hong Kong, P. R. China. wshun@graduate.hku.hk Chi-Kwong Li 2 Department of Mathematics,
More informationUNDERGROUND LECTURE NOTES 1: Optimality Conditions for Constrained Optimization Problems
UNDERGROUND LECTURE NOTES 1: Optimality Conditions for Constrained Optimization Problems Robert M. Freund February 2016 c 2016 Massachusetts Institute of Technology. All rights reserved. 1 1 Introduction
More informationGlobal output regulation through singularities
Global output regulation through singularities Yuh Yamashita Nara Institute of Science and Techbology Graduate School of Information Science Takayama 8916-5, Ikoma, Nara 63-11, JAPAN yamas@isaist-naraacjp
More information