ECEEN 5448 Fall 2011 Homework #5 Solutions
|
|
- Rebecca Potter
- 6 years ago
- Views:
Transcription
1 ECEEN 5448 Fall 211 Homework #5 Solutions Professor David G. Meyer December 8, Consider the 1-dimensional time-varying linear system ẋ t (u x) (a) Find the state-transition matrix, Φ(t, τ). Here A(t) t. Since the A-matrix is a scalar, of course A(t) and e commute. So in this case a fundamental matrix is Ψ(t) e A(θ) dθ e 2t/2 and the state transition matrix is thus Φ(t, τ) Ψ(t)Ψ 1 (τ) e 2t/2 e 2τ/2 e 2 (τ 2/ t 2/ ) You can also get this result by solving the IVP A(θ) dθ by separation of variables: dx dt tx; x(τ) x i dx dt tx dx x tdt log(x) 2t/2 + C x Ce 2t/2 and x(τ) x i tells us and so x i Ce 2τ/2 C x i e 2τ/2 x(t) e 2τ/2 e 2t/2 x i which leads to the same result for Φ(t, τ). (b) Given x() 1 and u(t) 1fort, find x(t) for t. The solution is x(t) Φ(t, )x()+ Φ(t, τ)b(τ)u(τ) dτ e 2t/2 + and so the solution when x() 1 and u(t) 1is x(t) 1 e 2 (τ /2 t /2 ) τdτ e 2t/2 +1 e 22t/2
2 2 2. Consider the nonlinear system where α is a real parameter. q αq q (a) Find all the equilibrium points. The equilibrium points are determined by the solutions to αq q q ( α q 2) q or q 2 α and so the equilibrium points are { when α q e, ± α when α> (b) Using linearization, classify, as much as possible, each equilibrium point found in (2a). so we get when α< when α when α> when α> when α> d ( αq q ) α q 2 dq q e is stable no information available. q e is critical case q e is unstable q e α is stable q e α is stable. We have the nonlinear system ẋ F (x, u) and F (, ). (a) If A F F x (, ) and b u (, ) and det(si A) s(s 1)(s + 2), then what can we conclude about the asymptotic stability of the equilibrium point x,u? A has an eigenvalue at s 1which is in the RHP. This means when u, the equilibrium point at x is unstable (b) Suppose that there is a k such that det(si A + bk) (s +7) 2 (s + 1). What can you say about the nonlinear system By the chain rule ẋ F (x, kx)? F(x, kx) (, ) F F (, ) + x x u A bk (, ) ( kx) (, ) x And A bk has all eigenvalues in the (open) LHP so for the nonlinear system ẋ F (x, kx) the equilibrium point at x is asymptotically stable. This little calculation is the basis for designing feedback laws for nonlinear system based on a linear approximation.
3 4. Let us say that a time-varying square matrix A(t) is constantly diagonable if, and only if, there is a constant invertible matrix, T, such that T 1 A(t)T D(t) for all t and D(t) is diagonal for all t (a) If A(t) is constantly diagonable, prove that A(t) and A(θ) dθ commute. Assume A(t) is constantly diagonable. First note that, since T does not depend on time, A(θ) dθ T 1 t D(θ) dθt then A(t) A(θ) dθ A(θ) dθa(t) T 1 D(t)TT 1 D(θ) dθt T 1 D(θ) dθt T 1 D(t)T T 1 D(t) D(θ) dθt T 1 D(θ) dθd(t)t but D(t) is diagonal, so D(t) D(θ) dθ D(θ) dθd(t), hence: T 1 D(θ) dθd(t)t T 1 D(θ) dθd(t)t (b) Find the state-transition matrix for the parametrically chirped harmonic oscillator system: 1 ẋ(t) 1+t 1 1 x(t)+ u(t) 1+t g(t) g(t) cos(g(t)) sin(g(t)) Hint: e which is a rotation matrix. And sin(g(t)) cos(g(t)) 1 cos(θ1 ) sin(θ 1 ) cos(θ2 ) sin(θ 2 ) cos(θ1 θ 2 ) sin(θ 1 θ 2 ) sin(θ 1 ) cos(θ 1 ) sin(θ 2 ) cos(θ 2 ) sin(θ 1 θ 2 ) cos(θ 1 θ 2 ) One can easily check that the A(t) for the parametrically chirped harmonic oscillator is constantly diagonable. Hence the fundamental matrix is Ψ(t) e A(θ) dθ and we just need to do some computations. First log(1 + t) A(θ) dθ log(1 + t) From the hint, then cos (log(1 + t)) sin (log(1 + t)) Ψ(t) sin (log(1 + t)) cos (log(1 + t)) Notice Ψ() I as expected. Now, then we find Φ(t, τ) from Φ(t, τ)ψ(t)ψ 1 cos (log(1 + t) log(1 + τ)) sin (log(1 + t) log(1 + τ)) (τ) sin (log(1 + t) log(1 + τ)) cos (log(1 + t) log(1 + τ)) Since log(a/b) log(a) log(b) you can simplify this somewhat: ( ) ( ) Φ(t, τ) cos log 1+t sin log 1+t ( ) ( ) sin log 1+t cos log 1+t
4 4 (c) Is the system asymptotically stable? Prove your answer. NO. Since cos (log(1 + t)) sin (log(1 + t)) lim Φ(t, ) lim t t sin (log(1 + t)) cos (log(1 + t)) and we see this limit is NOT zero (in fact, the limit does not even exist!). (d) If we add the output equation y(t) 1 x(t) (1) show we have an observable system on any interval, T 2 with T 2 >. On,T 2 the observability gramian is 2 2 W o,t 2 Φ T (t, )c T (t)c(t)φ(t, ) dt Φ T 1 (t, ) 1 Φ(t, ) dt In order for W o,t 2 to be singular, we d need to find a vector v such that cφ(t, )v for all t between and T 2 but Φ(t, ) is rotation through log(1 + t) radians. And since log(1), near t, Φ(t, ) rotates through as small an angle as we choose. Clearly one vector v cannot be orthogonal to both c and c rotated by a suitably small amount. Hence W o,t 2 is non-singular. Notice that the impulse responses for the system are ( ) 1+t h(t, τ)c(t)φ(t, τ)b(τ) cos log So a check for BIBO (Bounded-Input, Bounded-Output) stability is ( ) sup t 1+t log cos dz < 1+z 5. Deducing asymptotic stability through the output. (a) Consider the observable time-invariant linear system ẋ Ax x() x y cx Prove the following: If we find y(t) ast for all x then in fact e At as t. An easy way to proceed is by contradiction. Suppose e At doesn t go to zero as t. We know this means A has at least one eigenvalue with non-negative real part. Let λ be an eigenvalue of A with non-negative real part and let v be a corresponding eigenvector. Then if x v the output is y(t) ce At v cv e λ t Since {A, c} is observable, by PBH we know cv. Since Real (λ ) we know e λt does not go to zero as t. Thus y(t) does not go to zero as t which is a contradiction.
5 (b) Would this theorem generalize to time-varying systems? Consider ẋ A(t)x x() x y c(t)x and suppose we find y(t) ast for all x. Can we conclude Φ(t, τ) as t for each τ? Either prove it or find a counter-example. Since our proof in the time-invariant case relied on relating the asymptotic behavior of e At to the eigenvalues of A, we are immediately skeptical that this theorem will hold true in the time-varying case, where the aymptotic behavior of Φ(t, ) cannot be related to the eigenvalues of A(t). Consider the 1-dimensional system with A(t) 1and c(t) e 2t. The statetransition matrix is clearly { e t τ for t τ Φ(t, τ) otherwise and so here clearly lim t Φ(t, ). On the other hand, however, y(t) e 2t x(t) e 2t e t x e t x as t We just need to show that this {A(t),c(t)} pair is observable and we have our counterexample. The t o,t f Observability Gramian is W o t o,t f f t o Φ T (t, t o )c T (t)c(t)φ(t, t o ) dt f t o e 2t e 2(t to) e 4t dt f t o e 2t dt e 2t 2 f t o ( e 2t e 2t f ) which is clearly PD if t f >t o so the system is observable. e t to e 2t e 2t e t to dt
6 6 6. Given a T>and a function h(, ):R R R satisfying sup h 2 (t, τ) dτ < τ consider the operator H defined by Hu(t) h(t, τ)u(τ) dτ If we think of H as mapping L 2,TtoL 2,T. Compute H (the adjoint of H). H will be a map from L 2,T back to L 2,T. The defining property is y(t),hu( )(t) H y( )(t),u(t) The inner product is the standard L 2 inner product, so we require y(t) y(t)hu(t) dt h(t, τ)u(τ) dτ dt y(t)h(t, τ)u(τ) dτdt {y(t)h(t, τ) dt} u(τ) dτ H y(t)u(t) dt H y(t)u(t) dt H y(t)u(t) dt H y(t)u(t) dt (2) If we look at equation (6) realizing that t and τ are just dummy variables we can rewrite it as (y(t)h(t, z) dt) u(z) dz H y(z)u(z) dz and from this it is clear that H y(t) h(τ,t)y(τ) dτ Notice that this is much like finite-dimensional linear maps. For a matrix A ij A ji. If we think of h(t, τ) as a matrix with uncountably many rows and columns, the adjoint is found by transposition as well: h (t, τ) h(τ,t)!
7 7. A time-varying linear system is a so-called jump-parameter system and its A(t), b(t), and c(t) are given by b(t) 1 ; c(t) 1 1 ; A(t) for t<2 for t 2. (a) Find the state-transition matrix, Φ(t, τ). If t, τ < 2 then the system looks time invariant and Φ(t, τ)e A1(t τ) where A 1 1. Likewise, if t, τ 2 we have Φ(t, τ) e 1 A2(t τ) where A 2. It 1 is easy to calculate that 1 e t L 1 (si A) 1 1 e t 1 1 and e t 1 e t 1 e t Now, when τ 2 but t>2 we cross the jump boundary. The system state is propagated first by A 1 for 2 τ and then by A 2 for t 2 so Φ(t, τ)e A2(t 2) e A1(2 τ) 1 e t e t 2 e 2 τ Notice the ORDER here!! So, finally 1 e t τ 1 e t τ e 2 τ e t τ Φ(t, τ) 1 e τ t e 2 t 1 1 e e τ t 2 τ 1 e t τ 1 e t τ 1 e t τ for t, τ < 2 for τ 2 and t 2 for t<2 and τ 2 for t, τ 2 (b) Give an exhaustive enumeration of the regions in (t o,t f )-space where the Controllability Gramian, W c (t o,t f ), of the system is positive definite. The controllability Gramian will be PD on those intervals t o,t f where we can hit any desired end state, x(t f ) we want starting from any initial state x(t o ) we want. For t<2 the time-invariant system {A 1,b} is uncontrollable. For t 2 the TI system {A 2,b} is controllable. Thus W c (t o,t f ) is { not PD when tf < 2 PD whenever t f 2 (c) Repeat the above for the Observability Gramian, W o (t o,t f ). Similar reasoning shows that { PD when t W o (t o,t f ) is o < 2 not PD whenever t o 2
8 8. Show that if P 1 > and P 2 > (both are symmetric) then the product P 1 P 2, while neither symmetric nor PD, actually has real, positive eigenvalues. Since P 2 >, a symmetric and positive definite square root, P 1/2 2, exists. Since P 1/2 2 is PD, it is invertible so P 1/2 2 exists and is also symmetric and PD. Since similarity transform preserves eigenvalues, the eigenvalues of P 1 P 2 are the same as those of P 1/2 2 (P 1 P 2 )P 1/2 2 but P 1/2 2 (P 1 P 2 )P 1/2 2 P 1/2 2 P 1 P 1/2 2 and this matrix is symmetric. So the eigenvalues of P 1 P 2 are real. Further since P 1/2 2 is symmetric and P 1 is PD x T P 1/2 2 P 1 P 1/2 2 x (P 1/2 2 x) T (P 1 )(P 1/2 2 x) > if P 1/2 2 x but since P 1/2 2 is invertible, P 1/2 2 x x and hence P 1/2 2 P 1 P 1/2 2 is PD. So the eigenvalues of P 1 P 2 are positive as well. That completes the proof.
ECEEN 5448 Fall 2011 Homework #4 Solutions
ECEEN 5448 Fall 2 Homework #4 Solutions Professor David G. Meyer Novemeber 29, 2. The state-space realization is A = [ [ ; b = ; c = [ which describes, of course, a free mass (in normalized units) with
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 informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
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 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 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 information16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1
16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1 Charles P. Coleman October 31, 2005 1 / 40 : Controllability Tests Observability Tests LEARNING OUTCOMES: Perform controllability tests Perform
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 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 informationEE 380. Linear Control Systems. Lecture 10
EE 380 Linear Control Systems Lecture 10 Professor Jeffrey Schiano Department of Electrical Engineering Lecture 10. 1 Lecture 10 Topics Stability Definitions Methods for Determining Stability Lecture 10.
More information06/12/ rws/jMc- modif SuFY10 (MPF) - Textbook Section IX 1
IV. Continuous-Time Signals & LTI Systems [p. 3] Analog signal definition [p. 4] Periodic signal [p. 5] One-sided signal [p. 6] Finite length signal [p. 7] Impulse function [p. 9] Sampling property [p.11]
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 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 information6 Linear Equation. 6.1 Equation with constant coefficients
6 Linear Equation 6.1 Equation with constant coefficients Consider the equation ẋ = Ax, x R n. This equating has n independent solutions. If the eigenvalues are distinct then the solutions are c k e λ
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems
. AERO 632: Design of Advance Flight Control System Norms for Signals and. Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Norms for Signals ...
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 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 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 informationSupplementary chapters
The Essentials of Linear State-Space Systems Supplementary chapters J. Dwight Aplevich This document is copyright 26 2 J. D. Aplevich, and supplements the book The Ussentials of Linear State-Space Systems,
More information2. Review of Linear Algebra
2. Review of Linear Algebra ECE 83, Spring 217 In this course we will represent signals as vectors and operators (e.g., filters, transforms, etc) as matrices. This lecture reviews basic concepts from linear
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 information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 12: I/O Stability Readings: DDV, Chapters 15, 16 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology March 14, 2011 E. Frazzoli
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 informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 5. Input-Output Stability DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Input-Output Stability y = Hu H denotes
More informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty
More informationẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationModule 02 CPS Background: Linear Systems Preliminaries
Module 02 CPS Background: Linear Systems Preliminaries Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html August
More informationDepartment 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 informationLinear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems
Linear ODEs p. 1 Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Linear ODEs Definition (Linear ODE) A linear ODE is a differential equation taking the form
More informationApplied 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 informationCDS 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 informationProblem Set 5 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel
More informationAdvanced Mechatronics Engineering
Advanced Mechatronics Engineering German University in Cairo 21 December, 2013 Outline Necessary conditions for optimal input Example Linear regulator problem Example Necessary conditions for optimal input
More informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
More informationLinear dynamical systems with inputs & outputs
EE263 Autumn 215 S. Boyd and S. Lall Linear dynamical systems with inputs & outputs inputs & outputs: interpretations transfer function impulse and step responses examples 1 Inputs & outputs recall continuous-time
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 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 information10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections )
c Dr. Igor Zelenko, Fall 2017 1 10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections 7.2-7.4) 1. When each of the functions F 1, F 2,..., F n in right-hand side
More informationx(t) = t[u(t 1) u(t 2)] + 1[u(t 2) u(t 3)]
ECE30 Summer II, 2006 Exam, Blue Version July 2, 2006 Name: Solution Score: 00/00 You must show all of your work for full credit. Calculators may NOT be used.. (5 points) x(t) = tu(t ) + ( t)u(t 2) u(t
More informationObservability and state estimation
EE263 Autumn 2015 S Boyd and S Lall Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time observability
More informationEE221A Linear System Theory Final Exam
EE221A Linear System Theory Final Exam Professor C. Tomlin Department of Electrical Engineering and Computer Sciences, UC Berkeley Fall 2016 12/16/16, 8-11am Your answers must be supported by analysis,
More informationNonlinear 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 informationMATH 221: SOLUTIONS TO SELECTED HOMEWORK PROBLEMS
MATH 221: SOLUTIONS TO SELECTED HOMEWORK PROBLEMS 1. HW 1: Due September 4 1.1.21. Suppose v, w R n and c is a scalar. Prove that Span(v + cw, w) = Span(v, w). We must prove two things: that every element
More informationSmall oscillations and normal modes
Chapter 4 Small oscillations and normal modes 4.1 Linear oscillations Discuss a generalization of the harmonic oscillator problem: oscillations of a system of several degrees of freedom near the position
More informationCONTROL 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 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 informationProblem 2 (Gaussian Elimination, Fundamental Spaces, Least Squares, Minimum Norm) Consider the following linear algebraic system of equations:
EEE58 Exam, Fall 6 AA Rodriguez Rules: Closed notes/books, No calculators permitted, open minds GWC 35, 965-37 Problem (Dynamic Augmentation: State Space Representation) Consider a dynamical system consisting
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationSample Solutions of Assignment 9 for MAT3270B
Sample Solutions of Assignment 9 for MAT370B. For the following ODEs, find the eigenvalues and eigenvectors, and classify the critical point 0,0 type and determine whether it is stable, asymptotically
More informationKalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and
Kalman Decomposition Controllable / uncontrollable decomposition Suppose that the controllability matrix C R n n of a system has rank n 1
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 informationProf. 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 informationECE 388 Automatic Control
Controllability and State Feedback Control Associate Prof. Dr. of Mechatronics Engineeering Çankaya University Compulsory Course in Electronic and Communication Engineering Credits (2/2/3) Course Webpage:
More informationControl 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 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 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 informationEmpirical Gramians and Balanced Truncation for Model Reduction of Nonlinear Systems
Empirical Gramians and Balanced Truncation for Model Reduction of Nonlinear Systems Antoni Ras Departament de Matemàtica Aplicada 4 Universitat Politècnica de Catalunya Lecture goals To review the basic
More informationLecture 3: Review of Linear Algebra
ECE 83 Fall 2 Statistical Signal Processing instructor: R Nowak Lecture 3: Review of Linear Algebra Very often in this course we will represent signals as vectors and operators (eg, filters, transforms,
More informationState will have dimension 5. One possible choice is given by y and its derivatives up to y (4)
A Exercise State will have dimension 5. One possible choice is given by y and its derivatives up to y (4 x T (t [ y(t y ( (t y (2 (t y (3 (t y (4 (t ] T With this choice we obtain A B C [ ] D 2 3 4 To
More informationLecture Note #6 (Chap.10)
System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryook Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of
More informationLecture 3: Review of Linear Algebra
ECE 83 Fall 2 Statistical Signal Processing instructor: R Nowak, scribe: R Nowak Lecture 3: Review of Linear Algebra Very often in this course we will represent signals as vectors and operators (eg, filters,
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More informationThird In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 3 December 2009 (1) [6] Given that 2 is an eigenvalue of the matrix
Third In-Class Exam Solutions Math 26, Professor David Levermore Thursday, December 2009 ) [6] Given that 2 is an eigenvalue of the matrix A 2, 0 find all the eigenvectors of A associated with 2. Solution.
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationwe get y 2 5y = x + e x + C: From the initial condition y(0) = 1, we get 1 5 = 0+1+C; so that C = 5. Completing the square to solve y 2 5y = x + e x 5
Math 24 Final Exam Solution 17 December 1999 1. Find the general solution to the differential equation ty 0 +2y = sin t. Solution: Rewriting the equation in the form (for t 6= 0),we find that y 0 + 2 t
More informationTime Domain Analysis of Linear Systems Ch2. University of Central Oklahoma Dr. Mohamed Bingabr
Time Domain Analysis of Linear Systems Ch2 University of Central Oklahoma Dr. Mohamed Bingabr Outline Zero-input Response Impulse Response h(t) Convolution Zero-State Response System Stability System Response
More informationLeast Squares Based Self-Tuning Control Systems: Supplementary Notes
Least Squares Based Self-Tuning Control Systems: Supplementary Notes S. Garatti Dip. di Elettronica ed Informazione Politecnico di Milano, piazza L. da Vinci 32, 2133, Milan, Italy. Email: simone.garatti@polimi.it
More informationIntro. 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 informationMATH 23 Exam 2 Review Solutions
MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution
More informationLOCALLY POSITIVE NONLINEAR SYSTEMS
Int. J. Appl. Math. Comput. Sci. 3 Vol. 3 No. 4 55 59 LOCALLY POSITIVE NONLINEAR SYSTEMS TADEUSZ KACZOREK Institute of Control Industrial Electronics Warsaw University of Technology ul. Koszykowa 75 66
More informationLinear Differential Equations. Problems
Chapter 1 Linear Differential Equations. Problems 1.1 Introduction 1.1.1 Show that the function ϕ : R R, given by the expression ϕ(t) = 2e 3t for all t R, is a solution of the Initial Value Problem x =
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 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 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 informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationSTABILITY ANALYSIS. Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated using cones: Stable Neutral Unstable
ECE4510/5510: Feedback Control Systems. 5 1 STABILITY ANALYSIS 5.1: Bounded-input bounded-output (BIBO) stability Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated
More information3 Gramians and Balanced Realizations
3 Gramians and Balanced Realizations In this lecture, we use an optimization approach to find suitable realizations for truncation and singular perturbation of G. It turns out that the recommended realizations
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 informationLinear Systems. Chapter Basic Definitions
Chapter 5 Linear Systems Few physical elements display truly linear characteristics. For example the relation between force on a spring and displacement of the spring is always nonlinear to some degree.
More informationECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77
1/77 ECEN 605 LINEAR SYSTEMS Lecture 7 Solution of State Equations Solution of State Space Equations Recall from the previous Lecture note, for a system: ẋ(t) = A x(t) + B u(t) y(t) = C x(t) + D u(t),
More informationAdvanced Control Theory
State Space Solution and Realization chibum@seoultech.ac.kr Outline State space solution 2 Solution of state-space equations x t = Ax t + Bu t First, recall results for scalar equation: x t = a x t + b
More informationProfessor Fearing EECS120/Problem Set 2 v 1.01 Fall 2016 Due at 4 pm, Fri. Sep. 9 in HW box under stairs (1st floor Cory) Reading: O&W Ch 1, Ch2.
Professor Fearing EECS120/Problem Set 2 v 1.01 Fall 20 Due at 4 pm, Fri. Sep. 9 in HW box under stairs (1st floor Cory) Reading: O&W Ch 1, Ch2. Note: Π(t) = u(t + 1) u(t 1 ), and r(t) = tu(t) where u(t)
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 informationSolution for Homework 5
Solution for Homework 5 ME243A/ECE23A Fall 27 Exercise 1 The computation of the reachable subspace in continuous time can be handled easily introducing the concepts of inner product, orthogonal complement
More informationModule 09 From s-domain to time-domain From ODEs, TFs to State-Space Modern Control
Module 09 From s-domain to time-domain From ODEs, TFs to State-Space Modern Control Ahmad F. Taha EE 3413: Analysis and Desgin of Control Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/
More informationApplied Differential Equation. November 30, 2012
Applied Differential Equation November 3, Contents 5 System of First Order Linear Equations 5 Introduction and Review of matrices 5 Systems of Linear Algebraic Equations, Linear Independence, Eigenvalues,
More informationEE363 homework 2 solutions
EE363 Prof. S. Boyd EE363 homework 2 solutions. Derivative of matrix inverse. Suppose that X : R R n n, and that X(t is invertible. Show that ( d d dt X(t = X(t dt X(t X(t. Hint: differentiate X(tX(t =
More informationLecture 19 Observability and state estimation
EE263 Autumn 2007-08 Stephen Boyd Lecture 19 Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time
More informationL2 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 informationObservability. It was the property in Lyapunov stability which allowed us to resolve that
Observability We have seen observability twice already It was the property which permitted us to retrieve the initial state from the initial data {u(0),y(0),u(1),y(1),...,u(n 1),y(n 1)} It was the property
More informationAutonomous system = system without inputs
Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation
More informationDifferential Equations and Modeling
Differential Equations and Modeling Preliminary Lecture Notes Adolfo J. Rumbos c Draft date: March 22, 2018 March 22, 2018 2 Contents 1 Preface 5 2 Introduction to Modeling 7 2.1 Constructing Models.........................
More informationProperties of Transformations
6. - 6.4 Properties of Transformations P. Danziger Transformations from R n R m. General Transformations A general transformation maps vectors in R n to vectors in R m. We write T : R n R m to indicate
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 informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More informationMath 308 Practice Final Exam Page and vector y =
Math 308 Practice Final Exam Page Problem : Solving a linear equation 2 0 2 5 Given matrix A = 3 7 0 0 and vector y = 8. 4 0 0 9 (a) Solve Ax = y (if the equation is consistent) and write the general solution
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 informationNewtonian Mechanics. Chapter Classical space-time
Chapter 1 Newtonian Mechanics In these notes classical mechanics will be viewed as a mathematical model for the description of physical systems consisting of a certain (generally finite) number of particles
More informationEE Control Systems LECTURE 9
Updated: Sunday, February, 999 EE - Control Systems LECTURE 9 Copyright FL Lewis 998 All rights reserved STABILITY OF LINEAR SYSTEMS We discuss the stability of input/output systems and of state-space
More information