Theorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0).
|
|
- Octavia Hampton
- 6 years ago
- Views:
Transcription
1 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 = } and take Q H, so if Q is positive on all vectors then it s clearly positive on any subspace of V, so det(a k ) >. ( ): We know that det(a k ) > including det(a) >. By interlacing property, the matrix is positive definite.. Lyapunov functions Theorem. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < ). Idea is to transform basis in order to take a trajectory which doesn t admit the application of Lyapunov function to basis wherein the system decomposes in Jordan form, with each Jordan block admitting a quadratic form (Lyapunov function) for each direct summand; take Lyapunov function to be the direct sum of these pieces. E.g. If U = x 2 = x x, then U = ẋ x + x ẋ = x J x + x Jx, i.e. 2λ 2λ... 2λ ẋ = (λi k + N k )x () The determinant is 4λ 2, but this probably won t work, we should change the basis so that instead of having ones on the diagonal we have ɛ, which can be made arbitrarily small. Write x = y, x 2 = ay 2,..., x d a d y d, and also write, x k = λx k + x k so a k y k = λa k y k + a k 2 y k and y k = λy k + a y k. If λ / R, if A : R d R d, then P A (λ) = z d + tr(a)z d det(a); then P A (λ) = implies P A (λ) = too, which implies ( that nonreal ) complex roots always come in roots. If z λz, then a b the matrix of this operator is b a ( ) ( ) i Let A =, char A (z) = z 2 + 2z + 5 = (λ + + 2i)(λ + 2i), v = and ( ) ( ) ( ) ( ) 2i 2i 2 v 2 = ; A f =, let e + 2i = and e 2 = be the new basis. Then to calculate: ( ) ( ) ( ) ( ) = (3) /2 /2 2 (2)
2 Linear Systems Notes 2 This matrix is simply multiplication by 2i, as expected. If y = ay by 2 and ẏ 2 = by + ay 2, then (y 2 + y2 2 ). = 2ẏ y + 2ẏ 2 y 2. x Qx positive definite quadratic form such that d/dx(x Qx) = ẋ Qx + x Qẋ = x A Qx + x QAx = x (A Q + QA)x, negative definite, i.e. the derivative of a quadratic form is still quadratic. x Qx a x 2, for a > and d/dx(x Qx) a x 2. Then let take U(x(t)) = f(t), so f a x 2, then f a x 2 a 2 f, from which it follows that f(t) f()e a2t.
3 Linear Systems Notes 3 2 Lecture Theorem 2. Let ẋ = Ax. Then A is Hurwitz iff the origin is globally exponentially stable. Proof. ( ): Done. ( ): Assume that Re(λ), with λ an eigenvalue. Then we have ẋ = λx +..., ẋ 2 = λx , etc. On the other hand, if λ C \ R, then ẋ = ax bx +... and ẋ 2 = bx + ax 2. Corollary. Let ẋ = Ax+g(x) where g(x) c x 2, for small x, ẋ = f(x) = Ax+g(x) (Taylor), f() =. If A is Hurwitz, then any quadratic strict Lyapunov function for ẋ = Ax will be also strict Lyapunov for a neighborhood of x. Definition. Domain (region) of stability: If U is positive definite in Ω and U f(x) is negative definite in Ω; the region of attraction is the set of {ω Ω : x() = ω lim t x(t) }. Let ẋ = Ax, x R d. We want a quadratic form; want U = x P x, positive definit, and U Ax negative definite. Plugging in, we have d dt U(x(t)) = ẋ P x + x ẋ = x (A P + P A)x. We want A P +P A = Q with Q negative definite symmetric matrix. To solve there are d(d+) 2 equations to solve. Having solved for P, you then have to verify that P is positive definite. Proposition 2. If A is Hurwitz, then such P exists and is positive definite. Proof. Take x P x, and look at x(t) at the solution (trajectory) of ẋ = Ax, and x (t)qx(t)dt which exists. This integral is simply x e A t Qe At xdt = x e A t Qe At dtx. Claim: this integral (-) solves the Lyapunov equation. Positive definiteness of P follows immediately from Q being negative definite and its integral will still be negative definite. The only thing to check is that it solves the equation. d δ Calculate: dtp x(t) = lim e A t Qe At dt = Q, hence d/dtx P x = Q. δ δ 2. Discrete Time Systems Given ẋ = Ax, the analogous discrete time system is x[k + ] = Ax[k], and stability is given not by max Re(σ(A)) <, but that max σ(a) < (global exponential stability, or just exponential stability for nonlinear perturbations x[k + ] = Ax[k] + g(x[k])).
4 Linear Systems Notes Hamiltonian Systems Let ẋ = Ax and x R 2d, (q, ( p). ) Hamiltonian matrices are never Hurwitz. Example: q = ap + bq a b and ṗ + cp + dq, the matrix which is Hamiltonian iff tr(a) = ; can t have Hurwitz because c d R(λ + λ 2 ) < contradicts λ + λ 2 =. Nevertheless, such systems can still be Lyapunov stable (only Lyapunov stable), if det(a) >. If det(a) <, then the system is unstable.
5 Linear Systems Notes 5 3 Lecture 2 3. BIBO Stability Consider the LTV system ẋ = A(t)x + B(t)u (4) C(t)D(t)u (5) Question: if input is bounded, will output also be? Take first an example: say ẋ = v, and v = x + u, y = x; if ɛ sin(t) we have ẍ + x = ɛ sin t which has solution x(t) = A sin(t) + B cos(t) ± ɛt cos(t); hence with bounded u sup t [,T ] So to formalize this discussion, we consider bounded inputs in the L sense, namely u = u(t), if u < C R, will y < M C R as well? There is not much to say in this case, just take y(t) = t Cϕ(t, s)bu(s)ds + Du(t). For one thing, D should be bounded in the matrix norm, i.e. D(t) M D for all t; where M is as follows: For a matrix transformation M L(V, W ), where V comes equipped with a norm V and W as well with a norm W, then we define M M = sup Mv W. What remains is to demand that the integrand is bounded in v V = the operator norm. Precisely: Theorem 3. If there exists constants E, D R such that for all t, s the above bounded conditions hold, then the system is BIBO stable. Proposition 3. For a linear time invariant system (with ẋ = Ax + Bu), if A is Hurwitz then the system is BIBO stable. 3.2 Controllability Take ẋ = Ax + Bu and a finite time interval [, t ]. Define the reachable subspace R(, t ) = {x V : {u(t)} t [t,t ] s.t. x t = x t = x}. Note that R is a linear subspace in V, for if x, x 2 R(, t ), then λx + x 2 R as well because if you have controls u x and u 2 x 2, then λu + u 2 λx + x 2. ( ) Take ẋ = x. This system is asymptotically stable. If we add controls, which states 2 3 ( ) ( ) can we reach? Well if our control is say u, then if you start at the origin, since is an ( ) eigenvector of A, the trajectory x(t) will remain on that eigenvalue, i.e. traj(x(t)) span( ). ( ) If on the other hand, the control is given by u then you have control over one coordinate
6 Linear Systems Notes 6 direction of travel. This control does allow one to reach the whole space, for if you apply a control you can transverse trajectories then turn off the control and follow the natural dynamics for a period of time then turn on control again, etc. (Intuition for why R = R 2 here). On the other hand, we have controllable subspaces, defined similarly: C(, t ) = {x V : x( ) = x u t [t,t ] : x x(t ) = }. For ϕ(t, ), V t V t where C sits in the first space and R in the second. Recall ϕ(t, s) : V s V t satisfying ϕ t Define W r (, t ) = ϕ(t, s)b(s)b(s) ϕ (t, s)ds. Theorem 4. R(, t ) = Im(W R ). = Aϕ, and ϕ(t, t) = I. Proof. Define u(t) = B (t)ϕ (t, t)z, where assuming x W R z. Then simply compute t ϕ(t, s)b(s)b (s)ϕ(t, s x(t ) = W r z = x. In the other direction, if x / Im(W R ), then we want to show that X / R(, t ). But x Im(W R ) iff (c, x) = for all CW R = ; x = W R z cx = cw R z =. x / Im c row vector such that cx, cw R =, so = c t ϕ(t, s)b(s)b (s)ϕ(t, s)dsc = t (cϕ(t, s)b(s))(cϕb) ds = t cϕ(t, s)b ds which implies that all cϕ(t, s)b(s) is uniformly zero, so that x(t ) = ϕ(t, s)b(s)u(s)ds
7 Linear Systems Notes 7 4 Lecture 3 4. Lagrange Multiplier Given f : V R, with unconstrained optimization, one approach would be to find all the critical points and test whether they are local extremum. More generally, ( f x,..., f x d ), at a mininum point (assuming differentiability of f), f =. With constrained optimization, we have a set of conditions which we want to be satisfied, namely g i : V R (g : V V = R I, with (g) i = g i ), i =,..., I, with each g i = (g = ). Difficult problem. So we introduce a new vector, λ R I. We now consider the function F (x, λ) = f(x) + λg(x). Extremum points of f subject to g = corresponds exactly to extremum points of F (x, λ) on V W. In this case, F = ( F x,..., F x d, F λ,..., F λ I ) (6) The result is that f x j + λ i g x j =, with j =,..., d, g i (x) =. Optimization with Lagrange multipliers will help us in reachability. As usual we have ẋ = A(t)x(t) + B(t)u(t), with x V = R d, Take x = t ϕ(t, s)b(s)u(s)ds (7) {u} = U, where U is the space of all control functions. Given g : U V, with u ϕbuds, want g x = while minimizing the energy /2 u 2 ds (note that here what we want most, namely a solution to the dynamic system, gives us the constraint in the optimization problem). As before, λ V, form λ(g(u) x ) + 2 u uds = F (u, λ). The result t t λ( ϕ(t, s)b(s)u(s)ds x ) + /2 u (s)u(s)ds (8) equivalent to t /2u u + λϕ(t, s)b9s)u(s)ds λx (9) Now we need to take the gradient; the tricky part is the differentiating w.r.t. u, because it is a function. To address this difficulty, we consider the function at u, at u(s), some arbitrary s, denoted F u(s). F u(s) = u (s) + λϕ(t, s)b(s) () We know that if we solve this problem, then u (s) = λϕ(t, s)b(s), this is the optimal control. Alternatively u(s) = B (s)ϕ (t, s)z, with z = λ. The result: x = t ϕ(t, s)b(s)b (s)ϕ (t, s)dsz = W r z ()
8 Linear Systems Notes 8 The basis of this argument works so long as a feasible point exists. In other words, assuming that a minimal control exists, then it has this form. Controllability is similar. C(, t ) = {x : u x + ϕ(t, s)bu(s)ds = }. The controllability grammian W C = t ϕ(, s)b(s)b (s)ϕ (, s)ds. So we want W C = x V t. We have almost the same setup as before, = ϕ(t, )x + t ϕ(t, s)b(s)u(s)ds = g with u 2 ds minimum. The solution in this case follows the same process of computations as in reachability. 4.2 LTI Case Let ẋ = Ax + Bu, x V d, u U k, define the controllability matrix C = [B AB... A d B] which has kd columns. Theorem 5. R(, t ) = C(, t ) = ImC. Proof. Let x R, iff x = d t e A(t s) Bu(s)ds. By Cayley Hamilton, A d = cj A j so we can d express e A = αj A j, so t α j A j Bu(s)ds = A j B t α j (s)u(s)ds, and hence x ImC. In the other direction, let x ImC, and suppose that CR =. Want to show that c x =. R = ImW R, so = cw R = c t e A(t s) BB e A (t s) dsc = ce A(t s) B 2 ds ce A(t s) B =. Taking derivatives, of any order, will also be zero. So e.g. cb =, cae (t s) B = (when t = s, equality still holds) or cab =, etc. Hence cw R = ca k B =, but x = A k Bu k. Conclusion: if R = C = ImC = V. Theorem 6. LTI system is controllable iff there are no left eigenvectors of A in the kernel of B. u Ker(B), means Bu =, where B : U V. Similarly, B : V U, then c ker(b ) if cb =. Proof. dim R < d rank(c) < d c cc =. So cb =, i.e. c kerb. Now, we want to prove that among these c s, one of them is an eigenvector. Consider L = ker(c ) = {c : C c = } V. Claim: LA L; cb = cab =... = ca d B = (ca)b = (ca)ab = = (ca)a d B = c( α j A j )B =. Hence it s a linear operator on this space, and therefore has an eigenvector. In the other direction, if cb =, ca = λc, then cab = λcb =. Iterating, gives the same. Therefore, the rank cannot be full.
The motivating principal of linear systems is that even if a system isn t linear, it is at least locally linear.
Linear Systems Notes 1 1 Lecture 1 Introduction Consider the following linear system ẋ = Ax + Bx (1) where x R n, which describes the dynamics of x, given by two components: Ax the position of x (A L(R
More informationME 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 information1. 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 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 informationGrammians. Matthew M. Peet. Lecture 20: Grammians. Illinois Institute of Technology
Grammians Matthew M. Peet Illinois Institute of Technology Lecture 2: Grammians Lyapunov Equations Proposition 1. Suppose A is Hurwitz and Q is a square matrix. Then X = e AT s Qe As ds is the unique solution
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationLinear Control Theory
Linear Control Theory (Lecture notes) Version 0.7 Dmitry Gromov April 21, 2017 Contents Preface 4 I CONTROL SYSTEMS: ANALYSIS 5 1 Introduction 6 1.1 General notions................................ 6 1.2
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationModule 06 Stability of Dynamical Systems
Module 06 Stability of Dynamical Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha October 10, 2017 Ahmad F. Taha Module 06
More informationSolution 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 informationControl, Stabilization and Numerics for Partial Differential Equations
Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua
More informationMathematical Economics. Lecture Notes (in extracts)
Prof. Dr. Frank Werner Faculty of Mathematics Institute of Mathematical Optimization (IMO) http://math.uni-magdeburg.de/ werner/math-ec-new.html Mathematical Economics Lecture Notes (in extracts) Winter
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 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 informationControllability, Observability, Full State Feedback, Observer Based Control
Multivariable Control Lecture 4 Controllability, Observability, Full State Feedback, Observer Based Control John T. Wen September 13, 24 Ref: 3.2-3.4 of Text Controllability ẋ = Ax + Bu; x() = x. At time
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 informationCDS 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 informationExamples 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 information7.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 informationME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ
ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and
More informationThe 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 informationRobust Control 2 Controllability, Observability & Transfer Functions
Robust Control 2 Controllability, Observability & Transfer Functions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /26/24 Outline Reachable Controllability Distinguishable
More informationOPTIMAL CONTROL. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 28
OPTIMAL CONTROL Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 28 (Example from Optimal Control Theory, Kirk) Objective: To get from
More informationModule 03 Linear Systems Theory: Necessary Background
Module 03 Linear Systems Theory: Necessary Background Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html September
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 informationNonlinear 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 informationSTAT200C: Review of Linear Algebra
Stat200C Instructor: Zhaoxia Yu STAT200C: Review of Linear Algebra 1 Review of Linear Algebra 1.1 Vector Spaces, Rank, Trace, and Linear Equations 1.1.1 Rank and Vector Spaces Definition A vector whose
More informationModule 07 Controllability and Controller Design of Dynamical LTI Systems
Module 07 Controllability and Controller Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha October
More informationMath 113 Winter 2013 Prof. Church Midterm Solutions
Math 113 Winter 2013 Prof. Church Midterm Solutions Name: Student ID: Signature: Question 1 (20 points). Let V be a finite-dimensional vector space, and let T L(V, W ). Assume that v 1,..., v n is a basis
More informationEntrance 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 information1. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal
. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal 3 9 matrix D such that A = P DP, for A =. 3 4 3 (a) P = 4, D =. 3 (b) P = 4, D =. (c) P = 4 8 4, D =. 3 (d) P
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.4: Dynamic Systems Spring Homework Solutions Exercise 3. a) We are given the single input LTI system: [
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More information1 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 informationLinear System Theory
Linear System Theory Wonhee Kim Chapter 6: Controllability & Observability Chapter 7: Minimal Realizations May 2, 217 1 / 31 Recap State space equation Linear Algebra Solutions of LTI and LTV system Stability
More informationObservability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)
Observability Dynamic Systems Lecture 2 Observability Continuous time model: Discrete time model: ẋ(t) = f (x(t), u(t)), y(t) = h(x(t), u(t)) x(t + 1) = f (x(t), u(t)), y(t) = h(x(t)) Reglerteknik, ISY,
More information8 Periodic Linear Di erential Equations - Floquet Theory
8 Periodic Linear Di erential Equations - Floquet Theory The general theory of time varying linear di erential equations _x(t) = A(t)x(t) is still amazingly incomplete. Only for certain classes of functions
More informationNonlinear 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 informationHere each term has degree 2 (the sum of exponents is 2 for all summands). A quadratic form of three variables looks as
Reading [SB], Ch. 16.1-16.3, p. 375-393 1 Quadratic Forms A quadratic function f : R R has the form f(x) = a x. Generalization of this notion to two variables is the quadratic form Q(x 1, x ) = a 11 x
More informationCDS 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 informationNonlinear Control Lecture 5: Stability Analysis II
Nonlinear Control Lecture 5: Stability Analysis II Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 5 1/41
More informationControl Systems. Internal Stability - LTI systems. L. Lanari
Control Systems Internal Stability - LTI systems L. Lanari outline LTI systems: definitions conditions South stability criterion equilibrium points Nonlinear systems: equilibrium points examples stable
More informationDecay rates for partially dissipative hyperbolic systems
Outline Decay rates for partially dissipative hyperbolic systems Basque Center for Applied Mathematics Bilbao, Basque Country, Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ Numerical Methods
More informationMATH 4211/6211 Optimization Constrained Optimization
MATH 4211/6211 Optimization Constrained Optimization Xiaojing Ye Department of Mathematics & Statistics Georgia State University Xiaojing Ye, Math & Stat, Georgia State University 0 Constrained optimization
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 informationGQE ALGEBRA PROBLEMS
GQE ALGEBRA PROBLEMS JAKOB STREIPEL Contents. Eigenthings 2. Norms, Inner Products, Orthogonality, and Such 6 3. Determinants, Inverses, and Linear (In)dependence 4. (Invariant) Subspaces 3 Throughout
More informationw T 1 w T 2. w T n 0 if i j 1 if i = j
Lyapunov Operator Let A F n n be given, and define a linear operator L A : C n n C n n as L A (X) := A X + XA Suppose A is diagonalizable (what follows can be generalized even if this is not possible -
More informationNonlinear 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 informationSolution 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 informationMATH 5720: Unconstrained Optimization Hung Phan, UMass Lowell September 13, 2018
MATH 57: Unconstrained Optimization Hung Phan, UMass Lowell September 13, 18 1 Global and Local Optima Let a function f : S R be defined on a set S R n Definition 1 (minimizers and maximizers) (i) x S
More informationControl 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 informationNumerical Optimization
Constrained Optimization Computer Science and Automation Indian Institute of Science Bangalore 560 012, India. NPTEL Course on Constrained Optimization Constrained Optimization Problem: min h j (x) 0,
More informationTopic # /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 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 informationChap 3. Linear Algebra
Chap 3. Linear Algebra Outlines 1. Introduction 2. Basis, Representation, and Orthonormalization 3. Linear Algebraic Equations 4. Similarity Transformation 5. Diagonal Form and Jordan Form 6. Functions
More informationLecture 4 and 5 Controllability and Observability: Kalman decompositions
1 Lecture 4 and 5 Controllability and Observability: Kalman decompositions Spring 2013 - EE 194, Advanced Control (Prof. Khan) January 30 (Wed.) and Feb. 04 (Mon.), 2013 I. OBSERVABILITY OF DT LTI SYSTEMS
More informationIn particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with
Appendix: Matrix Estimates and the Perron-Frobenius Theorem. This Appendix will first present some well known estimates. For any m n matrix A = [a ij ] over the real or complex numbers, it will be convenient
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 informationModeling 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 / 57 Outline 1 Lecture 13: Linear System - Stability
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 informationEEE582 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 informationLINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM
LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM Unless otherwise stated, all vector spaces in this worksheet are finite dimensional and the scalar field F is R or C. Definition 1. A linear operator
More informationLecture 8 : Eigenvalues and Eigenvectors
CPS290: Algorithmic Foundations of Data Science February 24, 2017 Lecture 8 : Eigenvalues and Eigenvectors Lecturer: Kamesh Munagala Scribe: Kamesh Munagala Hermitian Matrices It is simpler to begin with
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationStabilization 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 informationj=1 x j p, if 1 p <, x i ξ : x i < ξ} 0 as p.
LINEAR ALGEBRA Fall 203 The final exam Almost all of the problems solved Exercise Let (V, ) be a normed vector space. Prove x y x y for all x, y V. Everybody knows how to do this! Exercise 2 If V is a
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 4: LMIs for State-Space Internal Stability Solving the Equations Find the output given the input State-Space:
More informationControl Systems. Internal Stability - LTI systems. L. Lanari
Control Systems Internal Stability - LTI systems L. Lanari definitions (AS) - A system S is said to be asymptotically stable if its state zeroinput response converges to the origin for any initial condition
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 informationMATH JORDAN FORM
MATH 53 JORDAN FORM Let A,, A k be square matrices of size n,, n k, respectively with entries in a field F We define the matrix A A k of size n = n + + n k as the block matrix A 0 0 0 0 A 0 0 0 0 A k It
More informationMA5510 Ordinary Differential Equation, Fall, 2014
MA551 Ordinary Differential Equation, Fall, 214 9/3/214 Linear systems of ordinary differential equations: Consider the first order linear equation for some constant a. The solution is given by such that
More informationELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky
ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky Equilibrium points and linearization Eigenvalue decomposition and modal form State transition matrix and matrix exponential Stability ELEC 3035 (Part
More information1 Continuous-time Systems
Observability Completely controllable systems can be restructured by means of state feedback to have many desirable properties. But what if the state is not available for feedback? What if only the output
More informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
More informationCriterions on periodic feedback stabilization for some evolution equations
Criterions on periodic feedback stabilization for some evolution equations School of Mathematics and Statistics, Wuhan University, P. R. China (Joint work with Yashan Xu, Fudan University) Toulouse, June,
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationNOTES ON LINEAR ODES
NOTES ON LINEAR ODES JONATHAN LUK We can now use all the discussions we had on linear algebra to study linear ODEs Most of this material appears in the textbook in 21, 22, 23, 26 As always, this is a preliminary
More informationNonlinear 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 informationSystems of Algebraic Equations and Systems of Differential Equations
Systems of Algebraic Equations and Systems of Differential Equations Topics: 2 by 2 systems of linear equations Matrix expression; Ax = b Solving 2 by 2 homogeneous systems Functions defined on matrices
More informationEigenvectors. Prop-Defn
Eigenvectors Aim lecture: The simplest T -invariant subspaces are 1-dim & these give rise to the theory of eigenvectors. To compute these we introduce the similarity invariant, the characteristic polynomial.
More informationQuadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix.
Quadratic forms 1. Symmetric matrices An n n matrix (a ij ) n ij=1 with entries on R is called symmetric if A T, that is, if a ij = a ji for all 1 i, j n. We denote by S n (R) the set of all n n symmetric
More informationMCE/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 informationNonlinear 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 information1 Overview. 2 A Characterization of Convex Functions. 2.1 First-order Taylor approximation. AM 221: Advanced Optimization Spring 2016
AM 221: Advanced Optimization Spring 2016 Prof. Yaron Singer Lecture 8 February 22nd 1 Overview In the previous lecture we saw characterizations of optimality in linear optimization, and we reviewed the
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 informationOn Controllability of Linear Systems 1
On Controllability of Linear Systems 1 M.T.Nair Department of Mathematics, IIT Madras Abstract In this article we discuss some issues related to the observability and controllability of linear systems.
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 19: Stabilization via LMIs Optimization Optimization can be posed in functional form: min x F objective function : inequality constraints
More informationEC Control Engineering Quiz II IIT Madras
EC34 - Control Engineering Quiz II IIT Madras Linear algebra Find the eigenvalues and eigenvectors of A, A, A and A + 4I Find the eigenvalues and eigenvectors of the following matrices: (a) cos θ sin θ
More informationMultivariable Control. Lecture 03. Description of Linear Time Invariant Systems. John T. Wen. September 7, 2006
Multivariable Control Lecture 3 Description of Linear Time Invariant Systems John T. Wen September 7, 26 Outline Mathematical description of LTI Systems Ref: 3.1-3.4 of text September 7, 26Copyrighted
More informationChap. 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 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 informationLyapunov 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 information1 Controllability and Observability
1 Controllability and Observability 1.1 Linear Time-Invariant (LTI) Systems State-space: Dimensions: Notation Transfer function: ẋ = Ax+Bu, x() = x, y = Cx+Du. x R n, u R m, y R p. Note that H(s) is always
More informationRobotics. 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 informationStabilization of persistently excited linear systems
Stabilization of persistently excited linear systems Yacine Chitour Laboratoire des signaux et systèmes & Université Paris-Sud, Orsay Exposé LJLL Paris, 28/9/2012 Stabilization & intermittent control Consider
More informationDesigning Information Devices and Systems II Spring 2018 J. Roychowdhury and M. Maharbiz Discussion 6B
EECS 16B Designing Information Devices and Systems II Spring 2018 J. Roychowdhury and M. Maharbiz Discussion 6B 1 Stability 1.1 Discrete time systems A discrete time system is of the form: xt + 1 A xt
More informationHalf 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 information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 8: Solutions of State-space Models Readings: DDV, Chapters 10, 11, 12 (skip the parts on transform methods) Emilio Frazzoli Aeronautics and Astronautics Massachusetts
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 informationLinear Algebra II Lecture 22
Linear Algebra II Lecture 22 Xi Chen University of Alberta March 4, 24 Outline Characteristic Polynomial, Eigenvalue, Eigenvector and Eigenvalue, Eigenvector and Let T : V V be a linear endomorphism. We
More information