Introduction to Modern Control MT 2016
|
|
- Pauline Strickland
- 5 years ago
- Views:
Transcription
1 CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate
2 Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear case and linearization procedure From continuous to discrete time: simulation of an ODE From state-space models to frequency domain: Transfer functions Alessandro Abate 1
3 Linear Ordinary Differential Equations (Linear ODE) Recall dynamical model for spring-damper system: q(t) = 1 m ( c( q(t)) kq(t) + u(t)) We ve seen it can be formulated as: { ẋ1 (t) = x 2 (t) ẋ 2 (t) = 1 m ( c(x 2(t)) kx 1 (t) + u(t)) Introduce linear approximation of nonlinear term (e.g., arctan): c(x 2 (t)) c x 2 (t) Consider output equation: y(t) = x 1 (t) Alessandro Abate 2
4 Consider state vector [ x1 (t) x 2 (t) ] Obtain state-space representation as linear ODE: d dt [ ] [ ] [ x1 (t) 0 1 x1 (t) = x 2 (t) k m c x m 2 (t) y(t) = [ 1 0 ] [ ] x 1 (t) x 2 (t) ] + [ 0 1 m ] u(t) if u(t) 0, t 0, ODE is said to be autonomous model is time-invariant: the behavior does not change under time-shifts (we thus usually pick initial time to be equal to 0) Alessandro Abate 3
5 General State-Space Model for Linear ODE Linear ODE in standard state-space form: x(t) R n, state variable u(t) R m, input variable y(t) R l, output variable d dt x(t) = Ax(t) + Bu(t), x(0) = x 0 y(t) = Cx(t) + Du(t) Dimensions of x, u, y dimensions (A, B, C, D) Alessandro Abate 4
6 General Model for Linear ODE: comments We shall often consider SISO systems (m = l = 1) Each integrator output is state variable One initial condition needed for each state variable ss structures Control System Toolbox of MATLAB Use of Simulink to represent these models Alessandro Abate 5
7 Block Diagram for Single Model (linear ODE) Linear system in state-space form 8 x 0 x(t) y(t) u(t) B + + C + + A D Each integrator output is state variable One d initial condition needed for each state variable dt x(t) = Ax(t) + Bu(t), x(0) = x 0 y(t) = Cx(t) + Du(t) Alessandro Abate 6
8 Block-Diagrams The block-diagram is obtained by interconnecting individual blocks (by The block-diagram is obtained by interconnecting individual blocks (by equating signals). Hence such diagrams (should) always correspond in equating signals). Hence such diagrams (should) always correspond in a precise fashion to algebraic equations. a precise fashion to algebraic equations. Block Diagrams for Models Interconnection Generic examples Generic examples inseries series interconnection: Series interconnection: u 1 y ẋ 1 = f(x 1,u 1 ẋ 2 2 = f(x 2,u ) ẋ 1 = f(x y 1 = h(x 1,u 1,u 1 ) ẋ 1 y 2 = f(x 2 = h(x 2,u 2,u 2 ) 2 ) u 1 y 2 y 1 = h(x 1,u 1 ) y 2 = h(x 2,u 2 ) Parallel interconnection: inparallel parallel interconnection: u u ẋ 1 y ẋ 1 1 y 1 ẋ 2 y 2 ẋ 2 y 2 = f(x 1,u) = h(x 1,u) = f(x 1,u) = h(x 1,u) = f(x 2,u) = h(x 2,u) = f(x 2,u) = h(x 2,u) + + y y 41/44 41/44 And of course in feedback (as discussed in previous lecture) Alessandro Abate 7
9 From high-order ODE to system of first order ODE Leitmotif of previous examples: exploit theory of first-order ODE Consider n-th order system: d n y dt n + a d n 1 y 1 dt n a ny = u Introduce variables: x 1 = d n 1 y/dt n 1, x 2 = d n 2 y/dt n 2,..., x n = y x 1 a 1 x 1 a 2 x 2... a n x n u d x 2 dt. x n 1 = x 1. x n x n x n 1 0 n-dim system of first-order ODE non-uniqueness of construction Alessandro Abate 8
10 More on Linearization: From Nonlinear to Linear ODE Recall model for Inverted Pendulum (Balance System, Segway) 36 CHAPTER 2. SYSTEM MODELING m θ l F M p (b) Saturn rocket (c) Cart pendulum system tems. (a) Segway Personal Transporter, (b) Saturn rocket and (c) art. Each of these examples uses forces at the bottom of the system Alessandro Abate (a) Segway (b) Saturn rock Figure 2.5: Balance systems. (a) Segway Pers inverted pendulum on a cart. Each of these exam to keep it upright. 9 a generalization of the spring mass system we saw earlier. ics for a mechanical system in the general form M(q)q + C(q, q ) + K (q) = B(q)u, Balance systems are a generalization of
11 More on Linearization: From Nonlinear to Linear ODE Overall state q = (p, θ), input u = F, output y = (p, θ) Can write dynamics as: M(q) q + K(q, q) = Bu It is often of interest to work around (equilibrium) point θ 0 sin θ θ, cos θ 1, θ 2 = o( θ) (exercise session will discuss details of linearisation) Alessandro Abate 10
12 Solution to Autonomous ODE Let us focus on the dynamical (autonomous) part of the ODE: d dt x(t) = Ax(t), x(0) = x 0 has solution where x(t) = e At x 0, (1) e At := I + At + A 2t2 2! + A3t3 3! +... Alessandro Abate 11
13 Matrix exponential: Moler and Van Loan, Nineteen Dubious Ways to Compute the Exponential of a Matrix, SIAM Review 20, 1978, pp Use of expm function in MATLAB Proof of (1): check if putative solution satisfies ODE, and i.c.: d dt x(t) = d dt ( e At x 0 ) = A(...), x(0) = x0 Assume A R n,n is diagonal, A = λ λ λ n Alessandro Abate 12
14 Then i = 1,..., n, x i (t) = e λ it x 0,i x(t) = What if A is not diagonal? e λ 1t 0 0 e λ 2t e λ nt x 0 State coordinate change via similarity transformation z = T x (T invertible n n) ż = T AT 1 z Ideal situation: T AT 1 is diagonal Procedure: study z-model, to derive properties for x-model Condition: A diagonalizable in MATLAB, command eig Alessandro Abate 13
15 Sufficient condition: A has distinct eigenvalues What if A is not diagonalizable? (e.g.: A has eigenvalues with multiplicity) For any A, existence of a Jordan form in MATLAB, command jordan Alessandro Abate 14
16 Transformation of state equations (state coordinate change) z = T x (T invertible n n) similarity transformation State space representation becomes ẋ = Ax + Bu, ż = T AT 1 z + T Bu, y = Cx + Du y = CT 1 z + Du (we will see shortly) two systems related by a similarity transformation have same I/O response CT 1 e T AT 1t T Bu(t) = Ce At Bu(t) Alessandro Abate 15
17 Eigenvalues and Dynamical Modes Importance of the eigenvalues of state matrix A: they are the solution of characteristic polynomial: det (λi A) = 0 in MATLAB poly (A) Alessandro Abate 16
18 Solution to ODE with Input The controlled ODE d dt x(t) = Ax(t) + Bu(t), x(0) = x 0 (assuming A is invertible) has solution x(t) = e At [x 0 + = e At x 0 + t 0 ] e Aλ Bu(λ)dλ t e A(t λ) Bu(λ)dλ } 0 {{} convolution of u(t) and e At B Alessandro Abate 17
19 Sum of response to initial condition and to input signal: superposition principle, which comes from model linearity Additionally, the output part: y(t) = Cx(t) + Du(t) has solution y(t) = Ce At x 0 + t 0 Ce A(t λ) Bu(λ)dλ + Du(t) Extension to non-linear, time-varying case: use of state-transition matrix Φ(t) (not in this course) Alessandro Abate 18
20 Impulse Response Assume for simplicity x 0 = 0. Impulse on (1-d) input u(t) = δ(t): y(t) = C t 0 e A(t λ) Bδ(λ)dλ + Dδ(t) = Ce At B + Dδ(t) 146 CHAPTER 5. LINEAR SYSTEMS Identification of system dynamics from data (recall end of Lec. 1): u Time t (a) Pulse and impulse functions 1 y 0.5 Pulse responses Impulse response t (b) Pulse and impulse responses Actually, LTI systems are completely characterized by impulse response Figure 5.6: Pulse response and impulse response. (a) The rectangles show pulses of width 5, 2.5 and 0.8, each with total area equal to 1. The arrow denotes an impulse δ(t) defined by equation (5.17). The corresponding pulse responses for a linear system with eigenvalues λ ={ 0.08, 0.62} are shown in (b) as dashed lines. The solid line is the true impulse response, which is well approximated by a pulse of duration 0.8. Alessandro Abate 19
21 Step Response - Transience and Steady State Step input (Heaviside function): u(t) = { 0, t = 0 1, t > 0 Assuming again x 0 = 0, we obtain y(t) = C t 0 e A(t λ) Bdλ + D = CA 1 e At B CA 1 B + D, t > 0 Presence of transient part, and steady-state term Alessandro Abate 20
22 Define controller specifications based on profile of step response Focus on dominant mode in model dynamics 5.3. INPUT/OUTPUT RESPONSE Overshoot M p Output Rise time T r Settling time T s Steady-state value y ss Time [sec] Figure 5.9: Sample step response. The rise time, overshoot, settling time and steady-state value give the key performance properties of the signal. equilibrium point in the absence of any input), then we can rewrite the solution as Alessandro Abate 21 y(t) = CA } 1 {{ e At B} + } D CA {{ 1 B }, t > 0. (5.22) transient steady-state The first term is the transient response and decays to zero as t. The second
23 Features: A Note on Non-Linear ODE non-linear dynamics and observations non-autonomous time-varying d dt x(t) = f(t, x(t), u(t)), x(t 0) = x 0 y(t) = g(t, x(t), u(t)) existence of non-standard cases (non-existence or non-uniqueness of solution) will resort to local linearization Alessandro Abate 22
24 Integration of ODE: Discrete-Time Systems Consider a general ODE of the form ẋ(t) = f(x(t)), x(0) = x 0 Introduce integration time-step h, 0 < h < Approximate time-derivative ẋ(t) at time t 0 as ẋ(t) x(t + h) x(t) (t + h) t This yields x(t + h) x(t) + hf(x(t)) Alessandro Abate 23
25 The above is a forward-euler first-order integration scheme Obtain ordinary difference equation describing discrete-time system WARNING: Property disruption through sampling Can employ more complex schemes (higher-order, variable step-size) 42 CHAPTER 2. SYSTEM MODELING 2 Position q [m] 1 0 h = 1!1 h = 0.5 h = 0.1 analytical! Time t [sec] Figure 2.9: Simulation of the forced spring mass system with different simulation time constants. The dashed line represents the analytical solution. The solid lines represent the approximate solution via the method of Euler integration, using decreasing step sizes. Integration of ODE in MATLAB performed with functions ode23, ode45 h gets smaller, the computed solution converges to the exact solution. The form L.F. Shampine, of the solution Numerical is also worth noticing: Solution after anof initial Ordinary transient, the system Differential settles into a periodic motion. The portion of the response after the transient is called the Equations. steady-state Chapman response to & thehall, input In addition to generating simulations, models can also be used to answer other types of questions. Two that are central to the methods described in this text concern Alessandro Abate the stability of an equilibrium point and the input/output frequency response. We illustrate these two computations through the examples below and return to the general computations in later chapters. Returning to the damped spring mass system, the equations of motion with no input forcing are given by 24
26 Transfer Functions and Frequency Domain Recall the block-diagram representation of dynamical models. Transfer Functions compactly represent the linear relationship between inputs and outputs of a system. As discussed, they can be easily constructed from measured data Start with a state-space representation ẋ = Ax + Bu, y = Cx + Du Consider inputs of the form u(t) = e st, s = σ + iω u(t) = e σt (cos ωt + i sin ωt) Solution is x(t) = e At x 0 + t 0 ea(t τ) Be sτ dτ = e At x 0 + e At t 0 e(si A)τ Bdτ If s λ(a) x(t) = e At x 0 + e At (si A) 1 ( e (si A)t I ) B Alessandro Abate 25
27 Thus, y(t) = Ce ( At x 0 (si A) 1 B ) + ( C(sI A) 1 B + D ) e st Consider steady-state term: let us call G yu (s) = C(sI A) 1 B + D the Transfer Function of the system: y(t) = G yu (s)u(t), u(t) = e st Notice that if s = λ i (A), G yu presents a singularity Alessandro Abate 26
28 Transfer Functions, example Consider canonical 2nd order system: q + 2ζω o q + ω 2 oq = kω 2 ou, y = q Can be re-written as: [ ] dx dt = 0 ω o x + ω o 2ζω o Assume ζ, ω o > 0 model is stable [ 0 kω o ] u, y = [1 0]x Steady-state response to u(t) = e st is C(sI A) 1 B: [1 0] ( si [ 0 ω o ω o 2ζω o ]) 1 [ 0 kω o ] = kω 2 o s 2 + 2ζω o s + ω 2 o Alessandro Abate 27
29 Transfer Functions, alternative derivation Consider integrable function f(t), f : R + R and s.t. f(t) < e s 0t Laplace transform L yields an F = Lf : C C, defined by F (s) = 0 e sτ f(τ)dτ, R(s) > s 0 Consider the system ẋ = Ax + Bu, y = Cx + Du Take Laplace transforms (with zero initial conditions): sx(s) = AX(s) + BU(s), Y (s) = CX(s) + DU(s) Substitute for X(s) to obtain: Y (s) = ( C(sI A) 1 B + D ) U(s) = G yu (s)u(s) Alessandro Abate 28
Linear 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 informationControl Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) and Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
More informationControl 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 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 informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationControl Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft Lecture 1 The concept of feedback The role of a controller What is a state? The concept of model in systems engineering: a
More informationControl Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) vs Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
More informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Outline of this module Instructors: A. Abate, P.J. Goulart, K. Margellos Teaching Assistants: A.
More informationControl Systems. Time response
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ẋ 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 informationTransfer function and linearization
Transfer function and linearization Daniele Carnevale Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata Corso di Controlli Automatici, A.A. 24-25 Testo del corso:
More informationControl 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 informationTransfer Functions. Chapter Introduction. 6.2 The Transfer Function
Chapter 6 Transfer Functions As a matter of idle curiosity, I once counted to find out what the order of the set of equations in an amplifier I had just designed would have been, if I had worked with the
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 informationLinear System Theory. Wonhee Kim Lecture 1. March 7, 2018
Linear System Theory Wonhee Kim Lecture 1 March 7, 2018 1 / 22 Overview Course Information Prerequisites Course Outline What is Control Engineering? Examples of Control Systems Structure of Control Systems
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More informationTopic # /31 Feedback Control Systems
Topic #7 16.30/31 Feedback Control Systems State-Space Systems What are the basic properties of a state-space model, and how do we analyze these? Time Domain Interpretations System Modes Fall 2010 16.30/31
More informationTopic # Feedback Control Systems
Topic #1 16.31 Feedback Control Systems Motivation Basic Linear System Response Fall 2007 16.31 1 1 16.31: Introduction r(t) e(t) d(t) y(t) G c (s) G(s) u(t) Goal: Design a controller G c (s) so that the
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 informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from
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 informationCDS 101/110: Lecture 3.1 Linear Systems
CDS /: Lecture 3. Linear Systems Goals for Today: Describe and motivate linear system models: Summarize properties, examples, and tools Joel Burdick (substituting for Richard Murray) jwb@robotics.caltech.edu,
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 informationCDS 101/110: Lecture 3.1 Linear Systems
CDS /: Lecture 3. Linear Systems Goals for Today: Revist and motivate linear time-invariant system models: Summarize properties, examples, and tools Convolution equation describing solution in response
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 informationSystems and Control Theory Lecture Notes. Laura Giarré
Systems and Control Theory Lecture Notes Laura Giarré L. Giarré 2017-2018 Lesson 7: Response of LTI systems in the transform domain Laplace Transform Transform-domain response (CT) Transfer function Zeta
More informationTime 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 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 informationDefinition of Reachability We begin by disregarding the output measurements of the system and focusing on the evolution of the state, given by
Chapter Six State Feedback Intuitively, the state may be regarded as a kind of information storage or memory or accumulation of past causes. We must, of course, demand that the set of internal states be
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 informationChapter Eight. Transfer Functions. 8.1 Frequency Domain Modeling
Chapter Eight Transfer Functions The typical regulator system can frequently be described, in essentials, by differential equations of no more than perhaps the second, third or fourth order. In contrast,
More informationControl 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 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 informationProblem Set 3: Solution Due on Mon. 7 th Oct. in class. Fall 2013
EE 56: Digital Control Systems Problem Set 3: Solution Due on Mon 7 th Oct in class Fall 23 Problem For the causal LTI system described by the difference equation y k + 2 y k = x k, () (a) By first finding
More informationEE263: Introduction to Linear Dynamical Systems Review Session 6
EE263: Introduction to Linear Dynamical Systems Review Session 6 Outline diagonalizability eigen decomposition theorem applications (modal forms, asymptotic growth rate) EE263 RS6 1 Diagonalizability consider
More informationLecture 2 ELE 301: Signals and Systems
Lecture 2 ELE 301: Signals and Systems Prof. Paul Cuff Princeton University Fall 2011-12 Cuff (Lecture 2) ELE 301: Signals and Systems Fall 2011-12 1 / 70 Models of Continuous Time Signals Today s topics:
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.6 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 5: Calculating the Laplace Transform of a Signal Introduction In this Lecture, you will learn: Laplace Transform of Simple
More informationChapter 1 Fundamental Concepts
Chapter 1 Fundamental Concepts 1 Signals A signal is a pattern of variation of a physical quantity, often as a function of time (but also space, distance, position, etc). These quantities are usually the
More informationLec 6: State Feedback, Controllability, Integral Action
Lec 6: State Feedback, Controllability, Integral Action November 22, 2017 Lund University, Department of Automatic Control Controllability and Observability Example of Kalman decomposition 1 s 1 x 10 x
More informationState Variable Analysis of Linear Dynamical Systems
Chapter 6 State Variable Analysis of Linear Dynamical Systems 6 Preliminaries In state variable approach, a system is represented completely by a set of differential equations that govern the evolution
More informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response
.. AERO 422: Active Controls for Aerospace Vehicles Dynamic Response Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. . Previous Class...........
More informationLecture 1 From Continuous-Time to Discrete-Time
Lecture From Continuous-Time to Discrete-Time Outline. Continuous and Discrete-Time Signals and Systems................. What is a signal?................................2 What is a system?.............................
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 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 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 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 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 informationAPPPHYS 217 Tuesday 6 April 2010
APPPHYS 7 Tuesday 6 April Stability and input-output performance: second-order systems Here we present a detailed example to draw connections between today s topics and our prior review of linear algebra
More informationReview: transient and steady-state response; DC gain and the FVT Today s topic: system-modeling diagrams; prototype 2nd-order system
Plan of the Lecture Review: transient and steady-state response; DC gain and the FVT Today s topic: system-modeling diagrams; prototype 2nd-order system Plan of the Lecture Review: transient and steady-state
More informationChapter 3. State Feedback - Pole Placement. Motivation
Chapter 3 State Feedback - Pole Placement Motivation Whereas classical control theory is based on output feedback, this course mainly deals with control system design by state feedback. This model-based
More informationCDS 101: Lecture 4.1 Linear Systems
CDS : Lecture 4. Linear Systems Richard M. Murray 8 October 4 Goals: Describe linear system models: properties, eamples, and tools Characterize stability and performance of linear systems in terms of eigenvalues
More informationDynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology.
Dynamic Response Assoc. Prof. Enver Tatlicioglu Department of Electrical & Electronics Engineering Izmir Institute of Technology Chapter 3 Assoc. Prof. Enver Tatlicioglu (EEE@IYTE) EE362 Feedback Control
More informationRepresent this system in terms of a block diagram consisting only of. g From Newton s law: 2 : θ sin θ 9 θ ` T
Exercise (Block diagram decomposition). Consider a system P that maps each input to the solutions of 9 4 ` 3 9 Represent this system in terms of a block diagram consisting only of integrator systems, represented
More informationLaplace Transform Part 1: Introduction (I&N Chap 13)
Laplace Transform Part 1: Introduction (I&N Chap 13) Definition of the L.T. L.T. of Singularity Functions L.T. Pairs Properties of the L.T. Inverse L.T. Convolution IVT(initial value theorem) & FVT (final
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52
1/52 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 2 Laplace Transform I Linear Time Invariant Systems A general LTI system may be described by the linear constant coefficient differential equation: a n d n
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 Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017
More informationEE/ME/AE324: Dynamical Systems. Chapter 7: Transform Solutions of Linear Models
EE/ME/AE324: Dynamical Systems Chapter 7: Transform Solutions of Linear Models The Laplace Transform Converts systems or signals from the real time domain, e.g., functions of the real variable t, to the
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 informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
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 informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries
. AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace
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 informationControls Problems for Qualifying Exam - Spring 2014
Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function
More informationLaboratory handouts, ME 340
Laboratory handouts, ME 340 This document contains summary theory, solved exercises, prelab assignments, lab instructions, and report assignments for Lab 4. 2014-2016 Harry Dankowicz, unless otherwise
More informationEE C128 / ME C134 Final Exam Fall 2014
EE C128 / ME C134 Final Exam Fall 2014 December 19, 2014 Your PRINTED FULL NAME Your STUDENT ID NUMBER Number of additional sheets 1. No computers, no tablets, no connected device (phone etc.) 2. Pocket
More informationMATH 4B Differential Equations, Fall 2016 Final Exam Study Guide
MATH 4B Differential Equations, Fall 2016 Final Exam Study Guide GENERAL INFORMATION AND FINAL EXAM RULES The exam will have a duration of 3 hours. No extra time will be given. Failing to submit your solutions
More information2 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 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 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 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 informationMath 3313: Differential Equations Second-order ordinary differential equations
Math 3313: Differential Equations Second-order ordinary differential equations Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Mass-spring & Newton s 2nd law Properties
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 informationSolutions of Spring 2008 Final Exam
Solutions of Spring 008 Final Exam 1. (a) The isocline for slope 0 is the pair of straight lines y = ±x. The direction field along these lines is flat. The isocline for slope is the hyperbola on the left
More informationAnalog Signals and Systems and their properties
Analog Signals and Systems and their properties Main Course Objective: Recall course objectives Understand the fundamentals of systems/signals interaction (know how systems can transform or filter signals)
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 18: State Feedback Tracking and State Estimation Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 18:
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
More informationState Regulator. Advanced Control. design of controllers using pole placement and LQ design rules
Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state
More informationAn Introduction for Scientists and Engineers SECOND EDITION
Feedback Systems An Introduction for Scientists and Engineers SECOND EDITION Karl Johan Åström Richard M. Murray Version v3.h (25 Sep 216) This is the electronic edition of Feedback Systems and is available
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 informationEECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators.
Name: SID: EECS C28/ ME C34 Final Wed. Dec. 5, 2 8- am Closed book. Two pages of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 2 2 6 3 4 4 5 6 6 7 8 2 Total
More informationContinuous Dynamics Solving LTI state-space equations גרא וייס המחלקה למדעי המחשב אוניברסיטת בן-גוריון
Continuous Dynamics Solving LTI state-space equations גרא וייס המחלקה למדעי המחשב אוניברסיטת בן-גוריון 2 State Space Models For a causal system with m inputs u t R m and p outputs y t R p, an nth-order
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 informationChapter Five. Linear Systems. 5.1 Basic Definitions
Feedback Systems by Astrom and Murray, v2.b http://www.cds.caltech.edu/~murray/fbswiki Chapter Five Linear Systems Few physical elements display truly linear characteristics. For example the relation between
More informationChap 4. State-Space Solutions and
Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations
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 informationPredictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado
Chapter 4 Dynamic Behavior Predictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado in Texas? Talk given by Edward Lorenz, December 972 meeting of the American Association for the
More informationChapter 3 Convolution Representation
Chapter 3 Convolution Representation DT Unit-Impulse Response Consider the DT SISO system: xn [ ] System yn [ ] xn [ ] = δ[ n] If the input signal is and the system has no energy at n = 0, the output yn
More informationChapter 1 Fundamental Concepts
Chapter 1 Fundamental Concepts Signals A signal is a pattern of variation of a physical quantity as a function of time, space, distance, position, temperature, pressure, etc. These quantities are usually
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 informationAnalysis and Design of Control Systems in the Time Domain
Chapter 6 Analysis and Design of Control Systems in the Time Domain 6. Concepts of feedback control Given a system, we can classify it as an open loop or a closed loop depends on the usage of the feedback.
More informationEE C128 / ME C134 Midterm Fall 2014
EE C128 / ME C134 Midterm Fall 2014 October 16, 2014 Your PRINTED FULL NAME Your STUDENT ID NUMBER Number of additional sheets 1. No computers, no tablets, no connected device (phone etc.) 2. Pocket calculator
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 informationACM/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 informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 8: Response Characteristics Overview In this Lecture, you will learn: Characteristics of the Response Stability Real
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 informationFEL3210 Multivariable Feedback Control
FEL3210 Multivariable Feedback Control Lecture 8: Youla parametrization, LMIs, Model Reduction and Summary [Ch. 11-12] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 8: Youla, LMIs, Model Reduction
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 information