Control Systems. Laplace domain analysis

Size: px
Start display at page:

Download "Control Systems. Laplace domain analysis"

Transcription

1 Control Systems Laplace domain analysis L. Lanari

2 outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output representation of the system through the transfer function exploring the structure of the transfer function solve a realization problem Lanari: CS - Laplace domain analysis

3 Laplace transform real valued function in the real variable t f(t) L F (s) complex valued function in the complex variable s t R s C F (s) = Z f(t)e st dt Laplace (unilateral) transform for f(t) with no impulse and no discontinuities at t = F (s) =L[f(t)] Im Region of existence: for all s with real part greater equal to an abscissa of convergence ¾ example f(t) =e at, a >, = a -a region of convergence Re Lanari: CS - Laplace domain analysis

4 inverse Laplace transform f(t) L F (s) L f(t) =L [F (s)] = j Z j j F (s)e st ds unique for one-sided functions f(t) defined for t or f(t) = for t < (one-to-one) all give the same F(s) with this assumption the inverse Laplace transform is unique F(s)? f(t) t Lanari: CS - Laplace domain analysis 4

5 Laplace transform Laplace transform of the impulse L[ (t)] = Z F (s) = Z (t)e st dt = f(t)e st dt Z (t)e st dt = e s = correct definition of the Laplace transform Linearity L[ f(t)+ g(t)] = L[f(t)] + L[g(t)] Derivative property L apple df (t) dt = sl[f(t)] f() can be applied iteratively L h f(t) i = s L[f(t)] sf() f() very useful in model derivation Lanari: CS - Laplace domain analysis 5

6 LTI systems also useful here ẋ(t) =Ax(t)+Bu(t) x() = x (proof): linearity + derivative property algebraic solution X(s) =(si A) x +(si A) BU(s) and therefore Y (s) =C(sI A) x + C(sI A) B + D U(s) Lanari: CS - Laplace domain analysis 6

7 LTI systems state ZIR transform state ZSR transform X(s) =(si A) x +(si A) BU(s) {z } {z } x(t) =e At x + Z t e A(t ) Bu( )d and therefore, comparing L e At =(si A) Heaviside step function L e at = s a L [ (t)] = s (t) for t for t< Lanari: CS - Laplace domain analysis 7 t

8 LTI systems y(t) =Ce At x + Z t w(t )u( )d with w(t) =Ce At B + D (t) transform W (s) =C(sI A) B + D Convolution theorem Y (s) = C(sI A) x + C(sI A) B + D U(s) = C(sI A) x + W (s) U(s) L applez t w(t )u( )d = W (s) U(s) being L[ (t)] = the output response transform corresponding to u(t) = ±(t) is indeed W(s) same for H(s) Lanari: CS - Laplace domain analysis 8

9 transfer function Input/Output behavior x = (ZSR) y ZS (t) = Z t w(t )u( )d Y ZS (s) =W (s) U(s) Transfer function W (s) = C(sI A) B + D = Y ZS(s) U(s) = L [w(t)] Input/Output behavior independent from state choice? independent from state space representation? Lanari: CS - Laplace domain analysis 9

10 transfer function (A, B, C, D) ẋ = Ax + Bu y = Cx + Du z = Tx det(t ) 6= ( e A, e B, e C, e D) ż = e Az + e Bu y = e Cz + e Du ea = TAT e B = TB e C = CT e D = D W (s) =C(sI A) B + D = e C(sI e A) e B + e D for a given system, the transfer function is unique Lanari: CS - Laplace domain analysis

11 shape of the transfer function W (s) =C(sI A) B + D inverse through the adjoint (transpose of the cofactor matrix) C (adjoint of (si A)) B + D det(si A) cofactor(i, j) =( ) i+j minor(i, j) polynomial of order n - polynomial of order n - det(si A) C (cofactor of (si A))T B + D polynomial of order n rational function (strictly proper) strictly proper rational function: proper rational function: degree of numerator < degree of denominator degree of numerator = degree of denominator Lanari: CS - Laplace domain analysis

12 shape of the transfer function W (s) = strictly proper rational function + D proper rational function W (s) = N(s) D(s) D = D 6= strictly proper rational function proper rational function W (s) = N(s) D(s) roots zeros poles (for coprime N(s) & D(s)) i.e. no common roots from previous analysis the poles are a subset of the eigenvalues of A {poles} {eigenvalues} more on this later Lanari: CS - Laplace domain analysis

13 partial fraction expansion (distinct roots case) Let F (s) = N(s) D(s) be a strictly proper rational function with coprime N(s) and D(s) and distinct roots of D(s) i.e. D(s) =a n n Y i= (s p i ) ) F (s) = N(s) a n Q n i= (s p i) then F(s) can be expanded as F (s) = nx i= s R i p i with the residues Ri computed as R i =[(s p i )F (s)] s=p i this result will be used for the transfer function W(s) output zero-state response Y(s) Lanari: CS - Laplace domain analysis

14 ( poles & eigenvalues W (s) = N(s) D(s)? from def W (s) = det(si A) N (s) characteristic polynomial if det(si A) and N (s) coprime {poles} = {eigenvalues} if det(si A) and N (s) not coprime {poles} subset of {eigenvalues} Input/Output representation Transfer function vs Input/State/Output representation State space here visible we need to understand when & why this happens (so to understand when we can consider the transfer function equivalent to a state space representation) Lanari: CS - Laplace domain analysis 4

15 poles & eigenvalues (distinct eigenvalues of A case) n = state space dimension = dimension of A = number of eigenvalues np = number of poles in W(s) W (s) = N(s) D(s) = N(s) Q np i= (s p i) = partial fraction expansion n p X R i p i s i= nx W (s) =L[w(t)] = L[Ce At B]=L[C@ e jt u j vj T j= A B] = done the same analysis in t for n = nx j= Cu j v T j B s j if v T j B = and/or Cu j = spectral form the eigenvalue j does not appear as a pole we have a hidden mode associated to the eigenvalue j (see structural properties) NB distinct eigenvalues is different from diagonalizable Lanari: CS - Laplace domain analysis 5

16 poles & eigenvalues (distinct eigenvalues of A case) fact: all the natural modes appear in the state ZIR for some generic initial condition e At If for an eigenvalue j we have that v T j B = implies the corresponding mode will not appear in the state impulsive response e At B = H(t) the corresponding mode (or eigenvalue) is said to be uncontrollable Cu j = implies the corresponding mode will not appear in the output transition matrix Ce At = (t) the corresponding mode (or eigenvalue) is said to be unobservable Lanari: CS - Laplace domain analysis 6

17 poles & eigenvalues (distinct eigenvalues of A case) Theorem Every pole is an eigenvalue. An eigenvalue i becomes a pole if and only if it is both controllable and observable vi T B 6= and Cu i 6= or equivalently the following two PBH rank tests are both verified rank A ii B = n and rank A C ii! = n Popov-Belevitch-Hautus controllability test Popov-Belevitch-Hautus observability test (the PBH test could be tested for a generic but matrix A - I loses rank only for = i ) Lanari: CS - Laplace domain analysis 7

18 poles & eigenvalues (distinct eigenvalues of A case) Where does the PBH test comes from? Observability (sketch): rank A C ii! <n means that the rectangular (n+) x n matrix has not full column rank and therefore it has a non-zero nullspace, that is there exists a n vector ui such that A C ii! u i = ( ( (A ii) u i = Au i = i u i Cu i = Cu i = that is there exists an eigenvector which belongs to the nullspace of C (or the corresponding mode is unobservable) Lanari: CS - Laplace domain analysis 8

19 example (si A) = with s+ A = (s+)(s ) s B = = M s + + M M = (s + )(si A) / = s= M = (s )(si A) / = s= C = s fractional decomposition works also for rational matrices a different way to compute the matrix exponential e At = M e t + M e t both natural modes appear (as it should be) in the state transition matrix Lanari: CS - Laplace domain analysis 9

20 example (si A) B = s+ mode corresponding to does not appear e At B W (s) = C(sI A) = s mode corresponding to does not appear Ce At no poles equivalently w(t) =Ce At B = forced response will always be zero independently from the input applied (look at the first order ODE) u / (A I)u = u / (A I)u = apple apple u = u = apple apple u = u = v T = / v T = / v T B 6= v T B = Cu = Cu 6= Lanari: CS - Laplace domain analysis

21 example equivalently with the PBH rank test controllability test mode corresponding to is controllable rank ==n mode corresponding to is uncontrollable rank =<n observability test rank C A =<n rank C A ==n mode corresponding to is unobservable mode corresponding to is observable Lanari: CS - Laplace domain analysis

22 poles & eigenvalues (general case) Theorem Every pole is an eigenvalue. An eigenvalue i becomes a pole with the multiplicity ma (algebraic multiplicity) if and only if both PBH rank tests are verified A ii rank A ii B = n rank C controllability observability! = n NB - If one of the two conditions is not verified then the eigenvalue i will appear as a pole with multiplicity strictly less than the algebraic multiplicity, possibly even (in this case we will have a hidden eigenvalue). In particular the eigenvalue will appear at most as a pole with multiplicity equal to its index (dimension of the largest Jordan block). NB - If for an eigenvalue i the geometric mg( i) > then there exists a hidden dynamics. Lanari: CS - Laplace domain analysis

23 example A = 4 5 (si A ) = (s ) PBH rank test verified for B = (s ) (s ) 6 4 (s ) 7 (s ) 5 (s ) C = 6= 8 8 6= easily seen from A I = 4 5 B = B = 4 5 C = F (s) = (s ) 4 5 C = C = F (s) = (s ) (s ) = index of s B = B = 4 5 C = F (s) = Lanari: CS - Laplace domain analysis

24 example A 4 = 4 5 since A 4 I = 4 5 the PBH rank test will never be satisfied independently from B and C. At most the eigenvalue will appear as a pole with multiplicity = index of = (si A 4 ) = (s ) (s ) (s ) 4 (s ) 5 = (s ) (s ) (s ) 4 (s ) 5 (s ) B 4 = 4 5 C 4 = F 4 (s) = (s ) B 5 = 4 5 C 5 = C = F 5 (s) = s (s ) = s B 6 = B = 4 5 C 6 = C = F 6 (s) = Lanari: CS - Laplace domain analysis 4

25 example A 7 = 4 5 (si A 7 ) = s 4 5 etc... the PBH rank test will never be satisfied independently from B and C. At most the eigenvalue will appear as a pole with multiplicity = index of = Lanari: CS - Laplace domain analysis 5

26 Laplace transform tables (t) sin t s + = /j s j /j s + j (t) s cos t s s + = / s j + / s + j e at s a sin( t + ) s sin + cos s + t k k! s k+ e at sin t (s a) + t k k! eat (s a) k+ e at cos t (s a) (s a) + Lanari: CS - Laplace domain analysis 6

27 realizations ok (A, B, C, D) W (s) =C(sI A) B + D how? W (s) =C(sI A) B + D realization infinite solutions (A, B, C, D) state dimension? how can we easily find one state space representation (A, B, C, D)? may be complicated for MIMO systems (here SISO) we see only one, obtainable directly from the coefficients of the transfer function (others are obtainable by simple similarity transformations) Lanari: CS - Laplace domain analysis 7

28 realizations W (s) = N(s) D(s) coprime N(s) & D(s) - find D W (s) = b n s n + b n s n + + b s + b s n + a n s n + a n s n + + a s + a + D - state with dimension n - one possible realization for A, B, C A c = 6 4 a a a a n C c = b b b b n 7 5 B c = controller canonical form (useful for eigenvalue assignment) Lanari: CS - Laplace domain analysis 8

Control Systems. Frequency domain analysis. L. Lanari

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

More information

Control Systems. System response. L. Lanari

Control Systems. System response. L. Lanari Control Systems m i l e r p r a in r e v y n is o System response L. Lanari Outline What we are going to see: how to compute in the s-domain the forced response (zero-state response) using the transfer

More information

Control Systems. Dynamic response in the time domain. L. Lanari

Control Systems. Dynamic response in the time domain. L. Lanari Control Systems Dynamic response in the time domain L. Lanari outline A diagonalizable - real eigenvalues (aperiodic natural modes) - complex conjugate eigenvalues (pseudoperiodic natural modes) - phase

More information

MULTIVARIABLE ZEROS OF STATE-SPACE SYSTEMS

MULTIVARIABLE ZEROS OF STATE-SPACE SYSTEMS Copyright F.L. Lewis All rights reserved Updated: Monday, September 9, 8 MULIVARIABLE ZEROS OF SAE-SPACE SYSEMS If a system has more than one input or output, it is called multi-input/multi-output (MIMO)

More information

Control Systems. Time response

Control 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

9.5 The Transfer Function

9.5 The Transfer Function Lecture Notes on Control Systems/D. Ghose/2012 0 9.5 The Transfer Function Consider the n-th order linear, time-invariant dynamical system. dy a 0 y + a 1 dt + a d 2 y 2 dt + + a d n y 2 n dt b du 0u +

More information

Control Systems. Time response. L. Lanari

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

More information

State will have dimension 5. One possible choice is given by y and its derivatives up to y (4)

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

Control Systems. Internal Stability - LTI systems. L. Lanari

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

Multivariable Control. Lecture 03. Description of Linear Time Invariant Systems. John T. Wen. September 7, 2006

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

Linear Systems. Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output

Linear Systems. Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output Linear Systems Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output Our interest is in dynamic systems Dynamic system means a system with memory of course including

More information

1 Continuous-time Systems

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

Control Systems. Linear Algebra topics. L. Lanari

Control Systems. Linear Algebra topics. L. Lanari Control Systems Linear Algebra topics L Lanari outline basic facts about matrices eigenvalues - eigenvectors - characteristic polynomial - algebraic multiplicity eigenvalues invariance under similarity

More information

ECE557 Systems Control

ECE557 Systems Control ECE557 Systems Control Bruce Francis Course notes, Version.0, September 008 Preface This is the second Engineering Science course on control. It assumes ECE56 as a prerequisite. If you didn t take ECE56,

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)

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

Identification Methods for Structural Systems

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

2.4 REALIZATION AND CANONICAL FORMS

2.4 REALIZATION AND CANONICAL FORMS Copyright F.L. Lewis 2008 All rights reserved Updated: Saturday, October 12, 2013 These notes are taken from F.L. Lewis, Applied Optimal Control and Estimation: Digital Design and Implementation, Prentice-Hall,

More information

Chap 4. State-Space Solutions and

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

Solution of Linear State-space Systems

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

More information

Perspective. ECE 3640 Lecture 11 State-Space Analysis. To learn about state-space analysis for continuous and discrete-time. Objective: systems

Perspective. ECE 3640 Lecture 11 State-Space Analysis. To learn about state-space analysis for continuous and discrete-time. Objective: systems ECE 3640 Lecture State-Space Analysis Objective: systems To learn about state-space analysis for continuous and discrete-time Perspective Transfer functions provide only an input/output perspective of

More information

Control Systems. Internal Stability - LTI systems. L. Lanari

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

Introduction & Laplace Transforms Lectures 1 & 2

Introduction & Laplace Transforms Lectures 1 & 2 Introduction & Lectures 1 & 2, Professor Department of Electrical and Computer Engineering Colorado State University Fall 2016 Control System Definition of a Control System Group of components that collectively

More information

Chapter 6: The Laplace Transform. Chih-Wei Liu

Chapter 6: The Laplace Transform. Chih-Wei Liu Chapter 6: The Laplace Transform Chih-Wei Liu Outline Introduction The Laplace Transform The Unilateral Laplace Transform Properties of the Unilateral Laplace Transform Inversion of the Unilateral Laplace

More information

Introduction to Modern Control MT 2016

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

Robust Control 2 Controllability, Observability & Transfer Functions

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

Observability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)

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

EE263: Introduction to Linear Dynamical Systems Review Session 6

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

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

João P. Hespanha. January 16, 2009

João P. Hespanha. January 16, 2009 LINEAR SYSTEMS THEORY João P. Hespanha January 16, 2009 Disclaimer: This is a draft and probably contains a few typos. Comments and information about typos are welcome. Please contact the author at hespanha@ece.ucsb.edu.

More information

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

Laplace Transforms Chapter 3

Laplace Transforms Chapter 3 Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important

More information

Intro. Computer Control Systems: F8

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

More information

Unit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace

Unit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace Unit : Modeling in the Frequency Domain Part : Engineering 81: Control Systems I Faculty of Engineering & Applied Science Memorial University of Newfoundland January 1, 010 1 Pair Table Unit, Part : Unit,

More information

ECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77

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

EE Experiment 11 The Laplace Transform and Control System Characteristics

EE Experiment 11 The Laplace Transform and Control System Characteristics EE216:11 1 EE 216 - Experiment 11 The Laplace Transform and Control System Characteristics Objectives: To illustrate computer usage in determining inverse Laplace transforms. Also to determine useful signal

More information

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

= e t sin 2t. s 2 2s + 5 (s 1) Solution: Using the derivative of LT formula we have

= e t sin 2t. s 2 2s + 5 (s 1) Solution: Using the derivative of LT formula we have Math 090 Midterm Exam Spring 07 S o l u t i o n s. Results of this problem will be used in other problems. Therefore do all calculations carefully and double check them. Find the inverse Laplace transform

More information

CONTROL DESIGN FOR SET POINT TRACKING

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

More information

Kalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and

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

Laplace Transforms and use in Automatic Control

Laplace Transforms and use in Automatic Control Laplace Transforms and use in Automatic Control P.S. Gandhi Mechanical Engineering IIT Bombay Acknowledgements: P.Santosh Krishna, SYSCON Recap Fourier series Fourier transform: aperiodic Convolution integral

More information

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

e st f (t) dt = e st tf(t) dt = L {t f(t)} s

e st f (t) dt = e st tf(t) dt = L {t f(t)} s Additional operational properties How to find the Laplace transform of a function f (t) that is multiplied by a monomial t n, the transform of a special type of integral, and the transform of a periodic

More information

Transfer function and linearization

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

Controllability, Observability, Full State Feedback, Observer Based Control

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

Basic Procedures for Common Problems

Basic Procedures for Common Problems Basic Procedures for Common Problems ECHE 550, Fall 2002 Steady State Multivariable Modeling and Control 1 Determine what variables are available to manipulate (inputs, u) and what variables are available

More information

ECE504: Lecture 9. D. Richard Brown III. Worcester Polytechnic Institute. 04-Nov-2008

ECE504: Lecture 9. D. Richard Brown III. Worcester Polytechnic Institute. 04-Nov-2008 ECE504: Lecture 9 D. Richard Brown III Worcester Polytechnic Institute 04-Nov-2008 Worcester Polytechnic Institute D. Richard Brown III 04-Nov-2008 1 / 38 Lecture 9 Major Topics ECE504: Lecture 9 We are

More information

Contents. 1 State-Space Linear Systems 5. 2 Linearization Causality, Time Invariance, and Linearity 31

Contents. 1 State-Space Linear Systems 5. 2 Linearization Causality, Time Invariance, and Linearity 31 Contents Preamble xiii Linear Systems I Basic Concepts 1 I System Representation 3 1 State-Space Linear Systems 5 1.1 State-Space Linear Systems 5 1.2 Block Diagrams 7 1.3 Exercises 11 2 Linearization

More information

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition

More information

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

1. Find the solution of the following uncontrolled linear system. 2 α 1 1 Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +

More information

EE/ME/AE324: Dynamical Systems. Chapter 7: Transform Solutions of Linear Models

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

MODELING OF CONTROL SYSTEMS

MODELING OF CONTROL SYSTEMS 1 MODELING OF CONTROL SYSTEMS Feb-15 Dr. Mohammed Morsy Outline Introduction Differential equations and Linearization of nonlinear mathematical models Transfer function and impulse response function Laplace

More information

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

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

More information

MATHEMATICAL MODELING OF CONTROL SYSTEMS

MATHEMATICAL MODELING OF CONTROL SYSTEMS 1 MATHEMATICAL MODELING OF CONTROL SYSTEMS Sep-14 Dr. Mohammed Morsy Outline Introduction Transfer function and impulse response function Laplace Transform Review Automatic control systems Signal Flow

More information

Lecture 3. Chapter 4: Elements of Linear System Theory. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.

Lecture 3. Chapter 4: Elements of Linear System Theory. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University. Lecture 3 Chapter 4: Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 3 p. 1/77 3.1 System Descriptions [4.1] Let f(u) be a liner operator, u 1 and u

More information

2.161 Signal Processing: Continuous and Discrete Fall 2008

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

Definition of the Laplace transform. 0 x(t)e st dt

Definition of the Laplace transform. 0 x(t)e st dt Definition of the Laplace transform Bilateral Laplace Transform: X(s) = x(t)e st dt Unilateral (or one-sided) Laplace Transform: X(s) = 0 x(t)e st dt ECE352 1 Definition of the Laplace transform (cont.)

More information

EE Control Systems LECTURE 9

EE 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

Multivariable Control. Lecture 05. Multivariable Poles and Zeros. John T. Wen. September 14, 2006

Multivariable Control. Lecture 05. Multivariable Poles and Zeros. John T. Wen. September 14, 2006 Multivariable Control Lecture 05 Multivariable Poles and Zeros John T. Wen September 4, 2006 SISO poles/zeros SISO transfer function: G(s) = n(s) d(s) (no common factors between n(s) and d(s)). Poles:

More information

ENGIN 211, Engineering Math. Laplace Transforms

ENGIN 211, Engineering Math. Laplace Transforms ENGIN 211, Engineering Math Laplace Transforms 1 Why Laplace Transform? Laplace transform converts a function in the time domain to its frequency domain. It is a powerful, systematic method in solving

More information

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52

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

MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction

MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 11/2/214 Outline Solving State Equations Variation

More information

EE 380. Linear Control Systems. Lecture 10

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

Linear Systems Theory

Linear Systems Theory ME 3253 Linear Systems Theory Review Class Overview and Introduction 1. How to build dynamic system model for physical system? 2. How to analyze the dynamic system? -- Time domain -- Frequency domain (Laplace

More information

Linear System Theory

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

Laplace Transform Part 1: Introduction (I&N Chap 13)

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

10 Transfer Matrix Models

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

More information

Name: Solutions Final Exam

Name: Solutions Final Exam Instructions. Answer each of the questions on your own paper. Put your name on each page of your paper. Be sure to show your work so that partial credit can be adequately assessed. Credit will not be given

More information

Lecture 7: Laplace Transform and Its Applications Dr.-Ing. Sudchai Boonto

Lecture 7: Laplace Transform and Its Applications Dr.-Ing. Sudchai Boonto Dr-Ing Sudchai Boonto Department of Control System and Instrumentation Engineering King Mongkut s Unniversity of Technology Thonburi Thailand Outline Motivation The Laplace Transform The Laplace Transform

More information

Laplace Transforms. Chapter 3. Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France

Laplace Transforms. Chapter 3. Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France Laplace Transforms Dr. M. A. A. Shoukat Choudhury 1 Laplace Transforms Important analytical

More information

Time Response Analysis (Part II)

Time Response Analysis (Part II) Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary

More information

EEE582 Homework Problems

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

More information

Dynamical system. The set of functions (signals) w : T W from T to W is denoted by W T. W variable space. T R time axis. W T trajectory space

Dynamical system. The set of functions (signals) w : T W from T to W is denoted by W T. W variable space. T R time axis. W T trajectory space Dynamical system The set of functions (signals) w : T W from T to W is denoted by W T. W variable space T R time axis W T trajectory space A dynamical system B W T is a set of trajectories (a behaviour).

More information

Observability and state estimation

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

I Laplace transform. I Transfer function. I Conversion between systems in time-, frequency-domain, and transfer

I Laplace transform. I Transfer function. I Conversion between systems in time-, frequency-domain, and transfer EE C128 / ME C134 Feedback Control Systems Lecture Chapter 2 Modeling in the Frequency Domain Alexandre Bayen Department of Electrical Engineering & Computer Science University of California Berkeley Lecture

More information

Stability. X(s) Y(s) = (s + 2) 2 (s 2) System has 2 poles: points where Y(s) -> at s = +2 and s = -2. Y(s) 8X(s) G 1 G 2

Stability. X(s) Y(s) = (s + 2) 2 (s 2) System has 2 poles: points where Y(s) -> at s = +2 and s = -2. Y(s) 8X(s) G 1 G 2 Stability 8X(s) X(s) Y(s) = (s 2) 2 (s 2) System has 2 poles: points where Y(s) -> at s = 2 and s = -2 If all poles are in region where s < 0, system is stable in Fourier language s = jω G 0 - x3 x7 Y(s)

More information

Module 4. Related web links and videos. 1. FT and ZT

Module 4. Related web links and videos. 1.  FT and ZT Module 4 Laplace transforms, ROC, rational systems, Z transform, properties of LT and ZT, rational functions, system properties from ROC, inverse transforms Related web links and videos Sl no Web link

More information

CHEE 319 Tutorial 3 Solutions. 1. Using partial fraction expansions, find the causal function f whose Laplace transform. F (s) F (s) = C 1 s + C 2

CHEE 319 Tutorial 3 Solutions. 1. Using partial fraction expansions, find the causal function f whose Laplace transform. F (s) F (s) = C 1 s + C 2 CHEE 39 Tutorial 3 Solutions. Using partial fraction expansions, find the causal function f whose Laplace transform is given by: F (s) 0 f(t)e st dt (.) F (s) = s(s+) ; Solution: Note that the polynomial

More information

Advanced Control Theory

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

The Laplace Transform

The Laplace Transform The Laplace Transform Introduction There are two common approaches to the developing and understanding the Laplace transform It can be viewed as a generalization of the CTFT to include some signals with

More information

Raktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response

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

Lecture 19 Observability and state estimation

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

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

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67

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

Transform Solutions to LTI Systems Part 3

Transform Solutions to LTI Systems Part 3 Transform Solutions to LTI Systems Part 3 Example of second order system solution: Same example with increased damping: k=5 N/m, b=6 Ns/m, F=2 N, m=1 Kg Given x(0) = 0, x (0) = 0, find x(t). The revised

More information

5. Observer-based Controller Design

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

More information

Design of Positive Linear Observers for Positive Systems via Coordinates Transformation and Positive Realization

Design of Positive Linear Observers for Positive Systems via Coordinates Transformation and Positive Realization Design of Positive Linear Observers for Positive Systems via Coordinates Transformation and Positive Realization Juhoon Back 1 Alessandro Astolfi 1 1 Department of Electrical and Electronic Engineering

More information

x 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7

x 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7 Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)

More information

f(t)e st dt. (4.1) Note that the integral defining the Laplace transform converges for s s 0 provided f(t) Ke s 0t for some constant K.

f(t)e st dt. (4.1) Note that the integral defining the Laplace transform converges for s s 0 provided f(t) Ke s 0t for some constant K. 4 Laplace transforms 4. Definition and basic properties The Laplace transform is a useful tool for solving differential equations, in particular initial value problems. It also provides an example of integral

More information

Lecture 7 and 8. Fall EE 105, Feedback Control Systems (Prof. Khan) September 30 and October 05, 2015

Lecture 7 and 8. Fall EE 105, Feedback Control Systems (Prof. Khan) September 30 and October 05, 2015 1 Lecture 7 and 8 Fall 2015 - EE 105, Feedback Control Systems (Prof Khan) September 30 and October 05, 2015 I CONTROLLABILITY OF AN DT-LTI SYSTEM IN k TIME-STEPS The DT-LTI system is given by the following

More information

L2 gains and system approximation quality 1

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

More information

Intro. Computer Control Systems: F9

Intro. Computer Control Systems: F9 Intro. Computer Control Systems: F9 State-feedback control and observers Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 21 dave.zachariah@it.uu.se F8: Quiz! 2 / 21 dave.zachariah@it.uu.se

More information

12/20/2017. Lectures on Signals & systems Engineering. Designed and Presented by Dr. Ayman Elshenawy Elsefy

12/20/2017. Lectures on Signals & systems Engineering. Designed and Presented by Dr. Ayman Elshenawy Elsefy //7 ectures on Signals & systems Engineering Designed and Presented by Dr. Ayman Elshenawy Elsefy Dept. of Systems & Computer Eng. Al-Azhar University Email : eaymanelshenawy@yahoo.com aplace Transform

More information

Topic # /31 Feedback Control Systems

Topic # /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 information

Systems Analysis and Control

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

Reachability, Observability and Minimality for a Class of 2D Continuous-Discrete Systems

Reachability, Observability and Minimality for a Class of 2D Continuous-Discrete Systems Proceedings of the 7th WSEAS International Conference on Systems Theory and Scientific Computation, Athens, Greece, August 24-26, 27 Reachability, Observability and Minimality for a Class of 2D Continuous-Discrete

More information

Modeling and Analysis of Dynamic Systems

Modeling and Analysis of Dynamic Systems Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability

More information

CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang

CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 5-1 Road Map of the Lecture V Laplace Transform and Transfer

More information

MA 527 first midterm review problems Hopefully final version as of October 2nd

MA 527 first midterm review problems Hopefully final version as of October 2nd MA 57 first midterm review problems Hopefully final version as of October nd The first midterm will be on Wednesday, October 4th, from 8 to 9 pm, in MTHW 0. It will cover all the material from the classes

More information

CLTI System Response (4A) Young Won Lim 4/11/15

CLTI System Response (4A) Young Won Lim 4/11/15 CLTI System Response (4A) Copyright (c) 2011-2015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2

More information

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

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

More information