Modeling and Analysis of Dynamic Systems

Size: px
Start display at page:

Download "Modeling and Analysis of Dynamic Systems"

Transcription

1 Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57

2 Outline 1 Lecture 13: Linear System - Stability Analysis 2 Definitions Equilibrium Sets 3 G. Ducard c 2 / 57

3 Outline 1 Lecture 13: Linear System - Stability Analysis 2 Definitions Equilibrium Sets 3 G. Ducard c 3 / 57

4 Zero Dynamics The dynamic behavior of linear system described as d dtx(t) = ẋ(t) = Ax(t)+Bu(t) y(t) = Cx(t)+Du(t) can be studied through its poles (eigenvalues of A) for the stability of the state vector x. See previous lecture. Let s consider the dynamic behavior when the output equation is considered y(t) through the output matrix C. In the Laplace domain, the relationship between input and output can be represented by a transfer function matrix: P(s) = C [si A] 1 B +D G. Ducard c 4 / 57

5 SISO Case: Transfer Function In the SISO case, the transfer function matrix: P(s) = C [si A] 1 B +D is a scalar rational transfer function which can always be written in the following form P(s) = Y(s) U(s) = k s n r +b n r 1 s n r b 1 s+b 0 s n +a n 1 s n 1 +a n 2 s n a 2 s 2 +a 1 s+a 0 G. Ducard c 5 / 57

6 SISO Case: Transfer Function P(s) = Y(s) U(s) = k s n r +b n r 1 s n r b 1 s+b 0 s n +a n 1 s n 1 +a n 2 s n a 2 s 2 +a 1 s+a 0 Discussions: The order of the highest power of s is n. Input gain: k The relative degree r: difference between highest power of s at denominator and the highest power of s at numerator. r plays an important role in the discussion of system zeros. A dynamic system can possess - not only poles - but also zeros. Question: What is the influence of the zeros? G. Ducard c 6 / 57

7 There is an equivalence between the transfer function P(s) = Y(s) U(s) = k s n r +b n r 1 s n r b 1 s+b 0 s n +a n 1 s n 1 +a n 2 s n a 2 s 2 +a 1 s+a 0 and its state-space representation (controller canonical form) d dt x(t) = x(t)+ 0 0 a 0 a 1 a 2... a n 1 k u(t) y(t) = [ b 0... b n r ] x(t) = Cx(t) The state vector is x(t) = [ x(t) ẋ(t) ẍ(t)... x (n) (t) ] Controller canonical form with gain k (min. number of parameters) The terms involved in the numerator are those of the C output vector ( transmission zeros ). G. Ducard c 7 / 57

8 Zero Dynamics - Definition The Zero Dynamics of a system: corresponds to its behavior for those special non-zero inputs u (t) and initial conditions x for which its output y(t) is identical to zero for a finite interval. 1 Study of the influence of the zeros on the dynamic properties of the system. 2 Study of the internal dynamics : analyze the stability of the system states, which are not directly controlled by the input u(t). G. Ducard c 8 / 57

9 Zero Dynamics - Problem In all the reference-tracking control problems, the controller tries to force the error to zero. ǫ = y ref y If y ref = 0, then y(t) is to be zero for all times, all its derivatives are to be zero as well. If a plant has internal dynamics, which are unstable, but not visible at the system s output, problems are to be expected. Let s see the form of the derivatives. G. Ducard c 9 / 57

10 The relative degree r is the number of differentiations needed to have the input u(t) explicitly appear in the output y (r) (t) y(t) = Cx(t), ẏ(t) = Cẋ(t) = CAx(t)+CBu(t) = CAx(t),. y (r 1) (t) = d dt y(r 2) (t) = CA r 1 x(t)+ca r 2 Bu(t) = CA r 1 x(t), y (r) (t) = d dt y(r 1) (t) = CA r x(t)+ca r 1 Bu(t) = CA r x(t)+ku(t) where k 0 and r n. G. Ducard c 10 / 57

11 The relative degree r z 1 = y = Cx = [b 0 x 1 +b 1 x b n r 1 x n r +x n r+1 ] z 2 = ẏ = CAx = [b 0 x 2 +b 1 x b n r 1 x n r+1 +x n r+2 ].. z r = y (r 1) = CA r 1 x = [b 0 x r +b 1 x r b n r 1 x n 1 +x n ] y (r) = CA r x+ku = [b 0 x r+1 +b 1 x r b n r x n +ẋ n ] ẋ n is found from the state-space representation: d dt x(t) = a 0 a 1 a 2... a n 1 x(t) k u(t) G. Ducard c 11 / 57

12 Outline 1 Lecture 13: Linear System - Stability Analysis 2 Definitions Equilibrium Sets 3 G. Ducard c 12 / 57

13 Formulation The following coordinate transformation z = Φx is introduced z 1 = y = Cx = [b 0 x 1 +b 1 x b n r 1 x n r +x n r+1 ] z 2 = ẏ = CAx = [b 0 x 2 +b 1 x b n r 1 x n r+1 +x n r+2 ]. z r = y (r 1) = CA r 1 x = [b 0 x r +b 1 x r b n r 1 x n 1 +x n ] G. Ducard c 13 / 57

14 Formulation The remaining n r coordinates are chosen such that the transformation Φ is regular and such that their derivatives also do not depend on the input u. Obviously the choice z r+1 = x 1 z r+2 = x 2... z n = x n r satisfies both requirements. To simplify notation the vector z is partitioned into two subvectors [ ] ξ z =, ξ = z 1..., η = z r+1... η z r z n G. Ducard c 14 / 57

15 Formulation In the new coordinates the system has the form [ ξ η ] = r T p T s T q T [ ξ η ] k u and obviously y = z 1 = ξ(1). In order to have an identically vanishing output it is therefore necessary and sufficient to choose the following control and initial conditions ξ = 0, u (t) = 1 k st η (t) where the initial condition η 0 0 can be chosen arbitrarilyġ. Ducard c 15 / 57

16 Formulation The internal states (zero dynamics states) evolve according to the equation d dt η(t) = Minimum phase system q T η (t) = Qη (t), η (0) = η 0 If the matrix Q is asymptotically stable (all eigenvalues with negative real part) then the system is minimum phase. Equivalence: a minimum phase system, is a system whose zeros have all negative real parts. These two definitions are consistent: see definition of the vector q T q T = [ b 0, b 1,..., b n r 2, b n r 1 ] G. Ducard c 16 / 57

17 Unstable Zero dynamics As soon as there is a zero with positive real part: system is non-minimum phase, system has its zero dynamics unstable, its internal states η can diverge without y(t) being affected. Consequences: the input u(t) may not be chosen such that the output y(t) is (almost) zero before the states η associated with the zero dynamics are (almost) zero. Feedback control more difficult to design. This imposes a constraint of the bandwidth of the closed-loop system: significantly slower than the slowest non-minimum phase zero. G. Ducard c 17 / 57

18 Formulation This equation can be derived using the definition of z n, the original system equation, the definition of z 1 (1), again the coordinate transformation, ż n = ẋ n r = x n r+1 = z 1 b 0 x 1... b n r 1 x n r = z 1 b 0 z r+1... b n r 1 z n = z 1 +q T η Therefore, the eigenvalues of Q coincide with the transmission zeros of the original system and with the roots of the numerator of its transfer function. G. Ducard c 18 / 57

19 Formulation u η n r 1 η s T 1 k ξ r ξ r 1 ξ 1 = q T r T p T G. Ducard c 19 / 57

20 Outline 1 Lecture 13: Linear System - Stability Analysis 2 Definitions Equilibrium Sets 3 G. Ducard c 20 / 57

21 Zero Dynamics Analysis on a Small SISO System P(s) = Y(s) U(s) = k s 2 +b 1 s+b 0 s 4 +a 3 s 3 +a 2 s 2 +a 1 s+a 0 Step 1: Convert the plant s transfer function into a state-space controller canonical form number of states n = 4, relative degree r = 2. d dt x(t) = a 0 a 1 a 2 a 3 x(t) k u(t) y(t) = [ b 0 b ] x(t) + [ 0 ] u(t) G. Ducard c 21 / 57

22 Step 2: Coordinate transformation Relative degree r = 2, therefore y(t) = b 0 x 1 (t)+b 1 x 2 (t)+x 3 (t) ẏ(t) = b 0 x 2 (t)+b 1 x 3 (t)+x 4 (t) ÿ(t) = a 0 x 1 (t) a 1 x 2 (t)+(b 0 a 2 )x 3 (t)+(b 1 a 3 )x 4 (t)+ku(t) The coordinate transformation z = Φ 1 x has the form z 1 = y = b 0 x 1 +b 1 x 2 +x 3 z 2 = ẏ = b 0 x 2 +b 1 x 3 +x 4 z 3 = x 1 z 4 = x 2 G. Ducard c 22 / 57

23 Step 3: Find the transformation matrices Φ 1, such that z = Φ 1 x and then compute Φ b 0 b b 0 b 1 1 Φ 1 = Φ = 1 0 b 0 b 1 and b 1 1 b 0 b 1 b 2 1 b 0 Remark: Notice that, by construction, det(φ) = det(φ 1 ) = 1 G. Ducard c 23 / 57

24 [ ξ Step 4: Build a new state-space representation in z = η ] ξ = [ z1 in the new coordinates the system z 2 ], η = [ z3 z 4 ] d dt d dt z(t) = Φ 1 AΦ z(t)+φ 1 Bu(t), y(t) = C Φ z (t) ξ 1 (t) ξ 2 (t) η 1 (t) = r 1 r 2 s 1 s ξ 1 (t) ξ 2 (t) η 1 (t) + 0 k 0 u(t) η 2 (t) 1 0 b 0 b 1 η 2 (t) 0 G. Ducard c 24 / 57

25 The coefficients r 1,r 2,s 1,s 2 are listed below r 1 = b 0 a 2 b 1 (b 1 a 3 ) r 2 = b 1 a 3 s 1 = b 0 b 1 (b 1 a 3 ) a 0 b 0 (b 0 a 2 ) s 2 = (b 1 a 3 )(b 2 1 b 0) a 1 (b 0 a 2 )b 1 G. Ducard c 25 / 57

26 u s 2 k η 2 η 1 s 1 ξ 2 ξ 1 b 1 r 2 b 0 r 1 Figure: System structure of the example s zero dynamics. G. Ducard c 26 / 57

27 Step 5: Study the submatrix Q of à = Φ 1 AΦ corresponding to the zero-dynamics vector η Choosing the following initial conditions ξ1 (0) = ξ 2 (0) = 0 and control signal u (t) = 1 k [s 1η1(t)+s 2 η2(t)] yields a zero output y(t) = 0 for all t 0. The initial conditions η1 (0) 0 and η 2 (0) 0 may be chosen arbitrarily. The trajectories of state variables η 1 (t) and η 2 (t), in this case, are defined by the equations [ d dt η (t) = 0 1 b 0 b 1 ] η (t) = Q η (t) G. Ducard c 27 / 57

28 Step 6: Conclude on the conditions to have Q asymptotically stable [ d dt η (t) = 0 1 b 0 b 1 ] η (t) = Q η (t) G. Ducard c 28 / 57

29 Definitions Equilibrium Sets Outline 1 Lecture 13: Linear System - Stability Analysis 2 Definitions Equilibrium Sets 3 G. Ducard c 29 / 57

30 Definitions Equilibrium Sets Lyapunov Stability Nonlinear differential equation d dt x(t) = f(x(t),u(t),t), x(t 0) = x 0 0 (1) Equilibrium x e is such that f(x e,t) = 0, t. G. Ducard c 30 / 57

31 Definitions Equilibrium Sets Stability of Equilibrium Points The point x e is Uniformly Lyapunov Stable (ULS) if for each scalar R > 0, there is a scalar r(r) > 0, such that if the initial condition x 0 satisfies x 0 < r, (2) the corresponding solution of (1) will satisfy the condition for all times t greater than t 0. x(t) < R (3) Remark: The notion uniformly indicates that neither r not R depend on t 0. G. Ducard c 31 / 57

32 Definitions Equilibrium Sets Lyapunov Stability For isolated equilibrium points: Lyapunov stability is the most useful stability concept. always connected to an equilibrium point x e of a system. The same point x e is asymptotically stable if it is ULS and attractive, i.e. lim t x(t) = x e G. Ducard c 32 / 57

33 Definitions Equilibrium Sets Lyapunov Stability A system is exponentially asymptotically stable if there exist constant scalars a > 0 and b > 0 such that Remarks: x(t) a e b (t t 0) x 0 1 In general, only an exponentially asymptotically stable equilibrium is acceptable for technical applications, 2 this form of stability is robust with respect to modeling errors. G. Ducard c 33 / 57

34 Definitions Equilibrium Sets Outline 1 Lecture 13: Linear System - Stability Analysis 2 Definitions Equilibrium Sets 3 G. Ducard c 34 / 57

35 Definitions Equilibrium Sets Equilibrium Sets For linear systems an equilibrium set can be: one isolated point; entire subspaces; periodic orbits (same frequency but arbitrary amplitude). G. Ducard c 35 / 57

36 Definitions Equilibrium Sets Equilibrium Sets Nonlinear systems can have (infinitely) many isolated equilibrium points. An equilibrium point can: have a finite region of attraction; be non-exponentially asymptotically stable; and be unstable the state of the system can escape to infinity in finite time. G. Ducard c 36 / 57

37 Outline 1 Lecture 13: Linear System - Stability Analysis 2 Definitions Equilibrium Sets 3 G. Ducard c 37 / 57

38 Stability of Nonlinear First-Order Systems First-order nonlinear systems are easy d x(t) = f(x(t),u(t)), x,u R dt because they can be separated for u(t) = 0 dx f(x) = dt = t+c The constant c is needed to satisfy the initial condition x(t 0 ) = x 0. G. Ducard c 38 / 57

39 Example: Non-Exponential Asymptotic Stability System d dt x(t) = x3 (t)+u(t), x(0) = x 0 Explicit solution for u(t) = 0 dx x 3 = and, therefore dt 1 2x 2 = t+c x(t) = x 0 (2 t x ) 1/2 Therefore, the system is asymptotically stable. G. Ducard c 39 / 57

40 Example: Non-Exponential Asymptotic Stability But the solution approaches the equilibrium slower than exponentially. To see this use the inequality x(t) = x 0 (2 t x ) 1/2 a e bt x 0 which yields 1 a e bt (2 t x )1/2 For limt, the right-hand side of this inequality tends to 0 for all possible 0 < a,b <. Contradiction! No exponential asymptotic stability. G. Ducard c 40 / 57

41 Outline 1 Lecture 13: Linear System - Stability Analysis 2 Definitions Equilibrium Sets 3 G. Ducard c 41 / 57

42 Stability of Second-Order Systems Two variables involved: x 1 and x 2 d dt x 1(t) = f 1 (x 1,x 2 ), x 1 (0) = x 1,0 (4) d dt x 2(t) = f 2 (x 1,x 2 ), x 2 (0) = x 2,0 G. Ducard c 42 / 57

43 Taylor expansion d dt x 1(t) = a 11 x 1 (t)+a 12 x 2 (t)+ f 1 (x 1,x 2 ), x 1 (0) = x 1,0 d dt x 2(t) = a 21 x 1 (t)+a 22 x 2 (t)+ f 2 (x 1,x 2 ), x 2 (0) = x 2,0 a i,j = f i f i (x 1,x 2 ), i,j = 1,2 lim = 0 x j x 0 x [ ] a11 a 12 A 2 2 = a 21 a 22 Linear system is given by: d dt δx(t) = A δx(t), δx(t) = [δx 1(t),δx 2 (t)] T (5) G. Ducard c 43 / 57

44 Outline 1 Lecture 13: Linear System - Stability Analysis 2 Definitions Equilibrium Sets 3 G. Ducard c 44 / 57

45 2nd-Order Systems Excluding the case where the matrix A has two eigenvalues with zero real part, the local behavior of the nonlinear system (4) and the behavior of the linear system (5) are topologically equivalent (i.e., they have the same geometric characteristics ). eigenvalues linearized system nonlinear system λ 1 C,λ 2 C Stable Focus Stable Focus λ 1 R,λ 2 R Stable Node Stable Node λ 1 R +,λ 2 R Saddle Saddle λ 1 R +,λ 2 R + Unstable Node Unstable Node λ 1 C +,λ 2 C + Unstable Focus Unstable Focus Re(λ 1,2 ) = 0 Center??? G. Ducard c 45 / 57

46 1 Stable Focus Stable Node Saddle x2 0 x2 0 x x x x Unstable Node Unstable Focus Center x2 0 x2 0 x x x x 1 G. Ducard c 46 / 57

47 Example Critical Nonlinear System d dt x 1 = x 1 +x 2 d dt x 2 = x 3 2 System has only one isolated equilibrium at x e,1 = x e,2 = 0. Linearization has one eigenvalue at 1 and one at 0. Observations: The solution of the linear system is stable (although not asymptotically stable), while the nonlinear system is unstable (it has even finite escape times). x 2,0 x 2 (t) = 1 2tx 2 2,0 With limt 1/(2x 2 2,0 ) x 2(t) escapes to infinity. G. Ducard c 47 / 57

48 General Systems 1 The is valid for all finite-order systems: as long as the linearized system has no eigenvalues on the imaginary axis, 2 The local stability properties of an arbitrary-order nonlinear system are fully understood once the eigenvalues of the linearization are known. 3 Particularly, if the linearization of a nonlinear system around an isolated equilibrium point x e is asymptotically stable (unstable resp.), then this equilibrium is an asymptotically stable (unstable resp.) equilibrium of the nonlinear system as well. G. Ducard c 48 / 57

49 General Systems Remarks: 1 Nothing is said about the large signal behavior. 2 In the critical case (Re{λ i } = 0), the nonlinear system s behavior is not defined by its first-order approximation, but by the higher-order terms f(x). G. Ducard c 49 / 57

50 Outline 1 Lecture 13: Linear System - Stability Analysis 2 Definitions Equilibrium Sets 3 G. Ducard c 50 / 57

51 Local stability properties of equilibrium x = 0 of system fully described by ẋ(t) = f(x(t)), x(0) 0 (6) A = f x x=0 provided A has no eigenvalues with zero real part. If this is not true or if one is interested in non-local results for more general systems, then use Lyapunov s direct method. G. Ducard c 51 / 57

52 Definitions A scalar function α(p) with α : R + R + is a nondecreasing function if α(0) = 0 and α(p) α(q) p > q. A function V : R n+1 R is a candidate global Lyapunov function if the function is strictly positive, i.e., V(x,t) > 0 x 0, t and V(0) = 0; and there are two nondecreasing functions α and β that satisfy the inequalities β( x ) V(x, t) α( x ). Obviously, if these conditions are met only in a neighborhood of the equilibrium point x = 0 only local assertions can be made. G. Ducard c 52 / 57

53 Theorem 1 The system ẋ(t) = f(x(t),t), x(t 0 ) = x 0 0 is globally/locally stable in the sense of Lyapunov if there is a global/local Lyapunov function candidate V(x, t) for which the following inequality holds true x(t) 0 and t Theorem 2 V(x(t),t) = V(x,t) t + V(x,t) f(x(t),t) 0 x The system (6) is globally/locally asymptotically stable if there is a global/local Lyapunov function candidate V(x, t) such that V(x(t),t) satisfies all conditions of a global/local Lyapunov function candidate. G. Ducard c 53 / 57

54 In general, it is very difficult to find a suitable function V(x,t) Hint: use physical insight, Lyapunov functions can be seen as generalized energy functions. Linear systems are simple V(x) = x T Px where P = P T > 0 is the solution of the Lyapunov equation PA+A T P = Q For arbitrary Q = Q T > 0, a solution to this equation exists iff A is Hurwitz matrix. Remark: Lyapunov theorems provide sufficient but not necessary conditions. Many extension proposed (LaSalle, Hahn, etc.). G. Ducard c 54 / 57

55 Example System d dt x 1 = x 1 (x 2 1 +x2 2 1) x 2 (7) d dt x 2 = x 1 +x 2 (x 2 1 +x2 2 1) One isolated equilibrium at x 1 = x 2 = 0. Linearizing the system around this equilibrium yields A = [ Eigenvalues are λ 1,2 = 1±j, system is locally asymptotically stable (Lyapunov principle). ] G. Ducard c 55 / 57

56 Example System has a finite region of attraction. Use a Lyapunov analysis to obtain a conservative estimation of the region of attraction. Candidate Lyapunov function V(x) = x 2 1 +x 2 2 Time derivative of V along a trajectory of (7) V(t) = 2(x 2 1 +x2 2 )(x2 1 +x2 2 1) Therefore, at least the region x 2 = x 2 1 +x2 2 < 1 must be part of the region of attraction. G. Ducard c 56 / 57

57 Vector field 0.1 f(x), x=end point x x 1 G. Ducard c 57 / 57

Modeling and Analysis of Dynamic Systems

Modeling and Analysis of Dynamic Systems Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 54 Outline 1 G. Ducard c 2 / 54 Outline 1 G. Ducard

More information

1 The Observability Canonical Form

1 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 information

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Time-varying Systems ẋ = f(t,x) f(t,x) is piecewise continuous in t and locally Lipschitz in x for all t 0 and all x D, (0 D). The origin

More information

Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1

Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1 Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium

More information

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.

Lecture 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 information

1. Find the solution of the following uncontrolled linear system. 2 α 1 1

1. 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 information

CONTROL DESIGN FOR SET POINT TRACKING

CONTROL DESIGN FOR SET POINT TRACKING Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded

More information

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points Nonlinear Control Lecture # 3 Stability of Equilibrium Points The Invariance Principle Definitions Let x(t) be a solution of ẋ = f(x) A point p is a positive limit point of x(t) if there is a sequence

More information

EN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015

EN 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 information

Prof. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait

Prof. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,

More information

1 Lyapunov theory of stability

1 Lyapunov theory of stability M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability

More information

Nonlinear Control. Nonlinear Control Lecture # 2 Stability of Equilibrium Points

Nonlinear Control. Nonlinear Control Lecture # 2 Stability of Equilibrium Points Nonlinear Control Lecture # 2 Stability of Equilibrium Points Basic Concepts ẋ = f(x) f is locally Lipschitz over a domain D R n Suppose x D is an equilibrium point; that is, f( x) = 0 Characterize and

More information

Lyapunov Stability Theory

Lyapunov 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 information

Zeros and zero dynamics

Zeros 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 information

Solution of Linear State-space Systems

Solution of Linear State-space Systems Solution of Linear State-space Systems Homogeneous (u=0) LTV systems first Theorem (Peano-Baker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state

More information

Department of Mathematics IIT Guwahati

Department of Mathematics IIT Guwahati Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,

More information

Nonlinear Control. Nonlinear Control Lecture # 6 Passivity and Input-Output Stability

Nonlinear Control. Nonlinear Control Lecture # 6 Passivity and Input-Output Stability Nonlinear Control Lecture # 6 Passivity and Input-Output Stability Passivity: Memoryless Functions y y y u u u (a) (b) (c) Passive Passive Not passive y = h(t,u), h [0, ] Vector case: y = h(t,u), h T =

More information

Topic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis

Topic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis Topic # 16.30/31 Feedback Control Systems Analysis of Nonlinear Systems Lyapunov Stability Analysis Fall 010 16.30/31 Lyapunov Stability Analysis Very general method to prove (or disprove) stability of

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A

More information

ECE504: 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 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 information

Nonlinear systems. Lyapunov stability theory. G. Ferrari Trecate

Nonlinear systems. Lyapunov stability theory. G. Ferrari Trecate Nonlinear systems Lyapunov stability theory G. Ferrari Trecate Dipartimento di Ingegneria Industriale e dell Informazione Università degli Studi di Pavia Advanced automation and control Ferrari Trecate

More information

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems

Introduction 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 information

L2 gains and system approximation quality 1

L2 gains and system approximation quality 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION L2 gains and system approximation quality 1 This lecture discusses the utility

More information

Some solutions of the written exam of January 27th, 2014

Some solutions of the written exam of January 27th, 2014 TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following

More information

CDS Solutions to Final Exam

CDS Solutions to Final Exam CDS 22 - Solutions to Final Exam Instructor: Danielle C Tarraf Fall 27 Problem (a) We will compute the H 2 norm of G using state-space methods (see Section 26 in DFT) We begin by finding a minimal state-space

More information

High-Gain Observers in Nonlinear Feedback Control. Lecture # 2 Separation Principle

High-Gain Observers in Nonlinear Feedback Control. Lecture # 2 Separation Principle High-Gain Observers in Nonlinear Feedback Control Lecture # 2 Separation Principle High-Gain ObserversinNonlinear Feedback ControlLecture # 2Separation Principle p. 1/4 The Class of Systems ẋ = Ax + Bφ(x,

More information

BIBO STABILITY AND ASYMPTOTIC STABILITY

BIBO STABILITY AND ASYMPTOTIC STABILITY BIBO STABILITY AND ASYMPTOTIC STABILITY FRANCESCO NORI Abstract. In this report with discuss the concepts of bounded-input boundedoutput stability (BIBO) and of Lyapunov stability. Examples are given to

More information

1. The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ = Ax + Bu is

1. The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ = Ax + Bu is ECE 55, Fall 2007 Problem Set #4 Solution The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ Ax + Bu is x(t) e A(t ) x( ) + e A(t τ) Bu(τ)dτ () This formula is extremely important

More information

Solution of Additional Exercises for Chapter 4

Solution of Additional Exercises for Chapter 4 1 1. (1) Try V (x) = 1 (x 1 + x ). Solution of Additional Exercises for Chapter 4 V (x) = x 1 ( x 1 + x ) x = x 1 x + x 1 x In the neighborhood of the origin, the term (x 1 + x ) dominates. Hence, the

More information

MCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory

MCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory MCE/EEC 647/747: Robot Dynamics and Control Lecture 8: Basic Lyapunov Stability Theory Reading: SHV Appendix Mechanical Engineering Hanz Richter, PhD MCE503 p.1/17 Stability in the sense of Lyapunov A

More information

Nonlinear Control Systems

Nonlinear Control Systems Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 7. Feedback Linearization IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs1/ 1 1 Feedback Linearization Given a nonlinear

More information

Pole placement control: state space and polynomial approaches Lecture 2

Pole placement control: state space and polynomial approaches Lecture 2 : state space and polynomial approaches Lecture 2 : a state O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.fr www.gipsa-lab.fr/ o.sename -based November 21, 2017 Outline : a state

More information

Convergence Rate of Nonlinear Switched Systems

Convergence 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 information

5. Observer-based Controller Design

5. Observer-based Controller Design EE635 - Control System Theory 5. Observer-based Controller Design Jitkomut Songsiri state feedback pole-placement design regulation and tracking state observer feedback observer design LQR and LQG 5-1

More information

MCE693/793: Analysis and Control of Nonlinear Systems

MCE693/793: Analysis and Control of Nonlinear Systems MCE693/793: Analysis and Control of Nonlinear Systems Lyapunov Stability - I Hanz Richter Mechanical Engineering Department Cleveland State University Definition of Stability - Lyapunov Sense Lyapunov

More information

Chap. 3. Controlled Systems, Controllability

Chap. 3. Controlled Systems, Controllability Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 1. Linear systems Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Linear systems 1.1 Differential Equations 1.2 Linear flows 1.3 Linear maps

More information

Time Response of Systems

Time Response of Systems Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =

More information

Stabilization and Passivity-Based Control

Stabilization and Passivity-Based Control DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive

More information

2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ

2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ 1 2 Lyapunov Stability Whereas I/O stability is concerned with the effect of inputs on outputs, Lyapunov stability deals with unforced systems: ẋ = f(x, t) (1) where x R n, t R +, and f : R n R + R n.

More information

LMI Methods in Optimal and Robust Control

LMI Methods in Optimal and Robust Control LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear

More information

Stability theory is a fundamental topic in mathematics and engineering, that include every

Stability theory is a fundamental topic in mathematics and engineering, that include every Stability Theory Stability theory is a fundamental topic in mathematics and engineering, that include every branches of control theory. For a control system, the least requirement is that the system is

More information

Control Systems. Laplace domain analysis

Control Systems. Laplace domain analysis Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output

More information

Lyapunov stability ORDINARY DIFFERENTIAL EQUATIONS

Lyapunov stability ORDINARY DIFFERENTIAL EQUATIONS Lyapunov stability ORDINARY DIFFERENTIAL EQUATIONS An ordinary differential equation is a mathematical model of a continuous state continuous time system: X = < n state space f: < n! < n vector field (assigns

More information

Nonlinear Control Lecture # 14 Input-Output Stability. Nonlinear Control

Nonlinear Control Lecture # 14 Input-Output Stability. Nonlinear Control Nonlinear Control Lecture # 14 Input-Output Stability L Stability Input-Output Models: y = Hu u(t) is a piecewise continuous function of t and belongs to a linear space of signals The space of bounded

More information

10 Transfer Matrix Models

10 Transfer Matrix Models MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important

More information

CDS Solutions to the Midterm Exam

CDS Solutions to the Midterm Exam CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2

More information

Lecture 8. Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control. Eugenio Schuster.

Lecture 8. Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control. Eugenio Schuster. Lecture 8 Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture

More information

STABILITY ANALYSIS OF DYNAMIC SYSTEMS

STABILITY ANALYSIS OF DYNAMIC SYSTEMS C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY STABILITY ANALYSIS OF DYNAMIC SYSTEMS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e

More information

Hybrid Systems Course Lyapunov stability

Hybrid 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 information

EE363 homework 7 solutions

EE363 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 information

Stability and Robustness Analysis of Nonlinear Systems via Contraction Metrics and SOS Programming

Stability and Robustness Analysis of Nonlinear Systems via Contraction Metrics and SOS Programming arxiv:math/0603313v1 [math.oc 13 Mar 2006 Stability and Robustness Analysis of Nonlinear Systems via Contraction Metrics and SOS Programming Erin M. Aylward 1 Pablo A. Parrilo 1 Jean-Jacques E. Slotine

More information

Robotics. Control Theory. Marc Toussaint U Stuttgart

Robotics. Control Theory. Marc Toussaint U Stuttgart Robotics Control Theory Topics in control theory, optimal control, HJB equation, infinite horizon case, Linear-Quadratic optimal control, Riccati equations (differential, algebraic, discrete-time), controllability,

More information

Half of Final Exam Name: Practice Problems October 28, 2014

Half of Final Exam Name: Practice Problems October 28, 2014 Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half

More information

Nonlinear Control Systems

Nonlinear Control Systems Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 3. Fundamental properties IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Example Consider the system ẋ = f

More information

Fall 線性系統 Linear Systems. Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian. NTU-EE Sep07 Jan08

Fall 線性系統 Linear Systems. Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian. NTU-EE Sep07 Jan08 Fall 2007 線性系統 Linear Systems Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian NTU-EE Sep07 Jan08 Materials used in these lecture notes are adopted from Linear System Theory & Design, 3rd.

More information

ME Fall 2001, Fall 2002, Spring I/O Stability. Preliminaries: Vector and function norms

ME Fall 2001, Fall 2002, Spring I/O Stability. Preliminaries: Vector and function norms I/O Stability Preliminaries: Vector and function norms 1. Sup norms are used for vectors for simplicity: x = max i x i. Other norms are also okay 2. Induced matrix norms: let A R n n, (i stands for induced)

More information

Lecture 4: Numerical solution of ordinary differential equations

Lecture 4: Numerical solution of ordinary differential equations Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor

More information

7 Planar systems of linear ODE

7 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 information

Examples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling

Examples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling 1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples

More information

Intro. Computer Control Systems: F8

Intro. Computer Control Systems: F8 Intro. Computer Control Systems: F8 Properties of state-space descriptions and feedback Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 22 dave.zachariah@it.uu.se F7: Quiz! 2

More information

Control Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli

Control Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli Control Systems I Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 13, 2017 E. Frazzoli (ETH)

More information

1 Relative degree and local normal forms

1 Relative degree and local normal forms THE ZERO DYNAMICS OF A NONLINEAR SYSTEM 1 Relative degree and local normal orms The purpose o this Section is to show how single-input single-output nonlinear systems can be locally given, by means o a

More information

Lecture 10: Singular Perturbations and Averaging 1

Lecture 10: Singular Perturbations and Averaging 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.243j (Fall 2003): DYNAMICS OF NONLINEAR SYSTEMS by A. Megretski Lecture 10: Singular Perturbations and

More information

Stability lectures. Stability of Linear Systems. Stability of Linear Systems. Stability of Continuous Systems. EECE 571M/491M, Spring 2008 Lecture 5

Stability 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 information

Equilibrium points: continuous-time systems

Equilibrium points: continuous-time systems Capitolo 0 INTRODUCTION 81 Equilibrium points: continuous-time systems Let us consider the following continuous-time linear system ẋ(t) Ax(t)+Bu(t) y(t) Cx(t)+Du(t) The equilibrium points x 0 of the system

More information

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Time-varying Systems ẋ = f(t,x) f(t,x) is piecewise continuous in t and locally Lipschitz in x for all t 0 and all x D, (0 D). The origin

More information

Modeling and Analysis of Dynamic Systems

Modeling and Analysis of Dynamic Systems Modeling and Analysis of Dynamic Systems by Dr. Guillaume Ducard Fall 2016 Institute for Dynamic Systems and Control ETH Zurich, Switzerland based on script from: Prof. Dr. Lino Guzzella 1/33 Outline 1

More information

System Types in Feedback Control with Saturating Actuators

System Types in Feedback Control with Saturating Actuators System Types in Feedback Control with Saturating Actuators Yongsoon Eun, Pierre T. Kabamba, and Semyon M. Meerkov Department of Electrical Engineering and Computer Science, University of Michigan, Ann

More information

Dynamical 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 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 information

High-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation

High-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation High-Gain Observers in Nonlinear Feedback Control Lecture # 3 Regulation High-Gain ObserversinNonlinear Feedback ControlLecture # 3Regulation p. 1/5 Internal Model Principle d r Servo- Stabilizing u y

More information

Stability of Dynamical systems

Stability of Dynamical systems Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u

More information

Stability in the sense of Lyapunov

Stability in the sense of Lyapunov CHAPTER 5 Stability in the sense of Lyapunov Stability is one of the most important properties characterizing a system s qualitative behavior. There are a number of stability concepts used in the study

More information

Stability of Parameter Adaptation Algorithms. Big picture

Stability of Parameter Adaptation Algorithms. Big picture ME5895, UConn, Fall 215 Prof. Xu Chen Big picture For ˆθ (k + 1) = ˆθ (k) + [correction term] we haven t talked about whether ˆθ(k) will converge to the true value θ if k. We haven t even talked about

More information

Nonlinear Autonomous Systems of Differential

Nonlinear Autonomous Systems of Differential Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such

More information

The goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T

The goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T 1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.

More information

Feedback Linearization

Feedback Linearization Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)

More information

Nonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México

Nonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre

More information

Massachusetts Institute of Technology. Department of Electrical Engineering and Computer Science : MULTIVARIABLE CONTROL SYSTEMS by A.

Massachusetts Institute of Technology. Department of Electrical Engineering and Computer Science : MULTIVARIABLE CONTROL SYSTEMS by A. Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Q-Parameterization 1 This lecture introduces the so-called

More information

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R

More information

Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.

Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4. Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the

More information

Complex Dynamic Systems: Qualitative vs Quantitative analysis

Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic

More information

7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system

7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system 7 Stability 7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system ẋ(t) = A x(t), x(0) = x 0, A R n n, x 0 R n. (14) The origin x = 0 is a globally asymptotically

More information

CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 101/110a: Lecture 2.1 Dynamic Behavior CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium

More information

ACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016

ACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016 ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe

More information

EEE582 Homework Problems

EEE582 Homework Problems EEE582 Homework Problems HW. Write a state-space realization of the linearized model for the cruise control system around speeds v = 4 (Section.3, http://tsakalis.faculty.asu.edu/notes/models.pdf). Use

More information

Control Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

Control Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich Control Systems I Lecture 5: Transfer Functions Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 20, 2017 E. Frazzoli (ETH) Lecture 5: Control Systems I 20/10/2017

More information

A conjecture on sustained oscillations for a closed-loop heat equation

A 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 information

Chapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems

Chapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.

More information

Control Systems. Frequency domain analysis. L. Lanari

Control Systems. Frequency domain analysis. L. Lanari Control Systems m i l e r p r a in r e v y n is o Frequency domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic

More information

1.7. Stability and attractors. Consider the autonomous differential equation. (7.1) ẋ = f(x),

1.7. Stability and attractors. Consider the autonomous differential equation. (7.1) ẋ = f(x), 1.7. Stability and attractors. Consider the autonomous differential equation (7.1) ẋ = f(x), where f C r (lr d, lr d ), r 1. For notation, for any x lr d, c lr, we let B(x, c) = { ξ lr d : ξ x < c }. Suppose

More information

Control Systems. Time response. L. Lanari

Control Systems. Time response. L. Lanari Control Systems Time response L. Lanari outline zero-state solution matrix exponential total response (sum of zero-state and zero-input responses) Dirac impulse impulse response change of coordinates (state)

More information

Control Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

Control Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich Control Systems I Lecture 7: Feedback and the Root Locus method Readings: Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 2, 2018 J. Tani, E. Frazzoli (ETH) Lecture 7:

More information

Theorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0).

Theorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0). Linear Systems Notes Lecture Proposition. A M n (R) is positive definite iff all nested minors are greater than or equal to zero. n Proof. ( ): Positive definite iff λ i >. Let det(a) = λj and H = {x D

More information

Balanced Truncation 1

Balanced 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 information

Output Feedback and State Feedback. EL2620 Nonlinear Control. Nonlinear Observers. Nonlinear Controllers. ẋ = f(x,u), y = h(x)

Output Feedback and State Feedback. EL2620 Nonlinear Control. Nonlinear Observers. Nonlinear Controllers. ẋ = f(x,u), y = h(x) Output Feedback and State Feedback EL2620 Nonlinear Control Lecture 10 Exact feedback linearization Input-output linearization Lyapunov-based control design methods ẋ = f(x,u) y = h(x) Output feedback:

More information

ECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27

ECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27 1/27 ECEN 605 LINEAR SYSTEMS Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability Feedback System Consider the feedback system u + G ol (s) y Figure 1: A unity feedback system

More information

Lyapunov Stability Analysis: Open Loop

Lyapunov Stability Analysis: Open Loop Copyright F.L. Lewis 008 All rights reserved Updated: hursday, August 8, 008 Lyapunov Stability Analysis: Open Loop We know that the stability of linear time-invariant (LI) dynamical systems can be determined

More information