Lecture Note #6 (Chap.10)

Size: px
Start display at page:

Download "Lecture Note #6 (Chap.10)"

Transcription

1 System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryook Yang

2 Chap. Model Approximation Model approximation Simplification, approximation and order reduction of models Model reduction Model order reduction in a linear system Balanced realization Pade approximation Moment matching Continued fraction approximation Model approximation of a nonlinear system by a linear one Linearization Describing function analysis Approximation of the nonlinear system by ignoring higher-order harmonics

3 Linearization Nonlinear differential equation: x= f( xu, ) Linearization: x f + f ( x, u )( x x ) + f ( x, u )( u u ) Discretization Sampling time: h xt () = Fx() t + Gu() t x + = Φx + Γu yt () = Cx() t yk = Cxk where h Fh Fs Φ e and Γ e Gds = = Some heuristic model reduction methods (For linear system only) H( z) =. z /(.7z.8 z ) =. z /(.8 z )(+. z ) Polynomial truncation H ( z ).3 z /(.7 z = ) Method of dominating poles H ( z ). z /(.8 z = ) Pole-zero cancellations 3 ( ). /(. ) x u k k k H z = z + z Cancel (-.8z - ) with (-.78z - )

4 Balanced Realization and Model Reduction The Reachability Gramian xk+ = Φxk + Γuk n S : xk R yk = Cxk What states can be reached with a given input energy assuming that x =? N T Given finite input energy: Ju ( ) = e = uu E uu k= k k States: x = N Φ N k Γu [ Φ N Φ N ][ u u N-] T U k = Γ = ψ = T T Input Sequence for desired x N : U = ( ψ ψ ) ψ x N k N N N T T T k= k k N N N N N Ju ( ) = uu = U U = x P x E T N k T T k N N N k = N N N N N where P = ψ ψ = Φ ΓΓ ( Φ ) Bound on the reachable states at time N P N : Reachability Gramian T T T T P = N ΦPNΦ + ΓΓ + ΦPΦ P + ΓΓ = Lyapunov equation

5 The Observability Gramian What states energy is necessary for all u k = in order to obtain a specified output energy? Specified output energy: J T ( y ) = e = y y = E T J( y) = y y = xq x = E k = k k T yy k= k k k yk = Cxk = CΦxk = = CΦ x ( uk =, k) where ( T ) k T k Q = Φ CCΦ k= Q: Observability Gramian T T T T Q = Φ QΦ + CC Φ QΦ Q+ CC= Lyapunov equation ( Φ T QΦ + CC T = ( Φ T ) k+ CC T Φ k+ + CC T = Φ ) k = Reachability Gramian decides the minimum input energy to drive a system from x to x N. Observability Gramian decides the maximum output energy that the initial state x can generate.

6 Balanced Realization Reachability and observability Gramians P and Q define matrices that describe the sensitivity of the input-output map in different directions of state space. Consider a state-space transformation z k =Tx k, then the system becomes T zk+ = TΦT zk + TΓuk Pz = TPT S : T yk = CT zk Qz = T QT Different state-space realization may result in different Gramians Balanced realization: P z =Q z P = Q = Σ = diag( σ, σ ) (where σ = λ ( PQ)) z z i i Algorithm using Cholesky factorization T Q = QQ QPQ = UΣ U ( UU = I) Σ T T T = Σ Σ T Hankelsingular value T = Σ U Q T If some σ i s are relatively small, the Σ can be reduced and the resulting transformation produces a balanced model resuction.

7 Example. (Balanced model reduction) S : H( z) =. z /(.7z.8 z ) Controllable canonical realization.7.8 xk+ = xk + uk yk = Balanced realization.57.7 P = Q = Q =.3 T = Σ = z = z + u k+ k k (. ) y Σ k x k From Lyapunov equation.743 =.3 ( ) = z k

8 Reduced model Since first singular value is dominant, second variable can be reduced as if it has no dynamics (z k+ = z k )..55 Σ z = k+ Φ Φ zk Γ.69 = u + z Φ k Φ z Γ + k z = Φ z + Φ z + Γ u z = ( I Φ ) ( Φ z + Γ u ) k k k k k k k z = Φ z + Φ ( I Φ ) ( Φ z + Γ u ) + Γ u k+ k k k k = ( Φ + Φ ( I Φ ) Φ ) z + ( Γ + Φ ( I Φ ) Γ ) u k y = Cz + C ( I Φ ) ( Φ z + Γ u ) k k k k = ( C + C ( I Φ ) Φ ) z + C ( I Φ ) Γ u k z =.7976z u k k k y =.4478z +.95u k k k k k k

9 Balanced model reduction for continuous systems x A A x B x u = + y = ( C C) + Du x A A x B x x is assumed to have no dynamics x = and x = A A x A B u This method preserves the essential low-frequency properties (gain). The reduced-order state-space model x = ( A A A A ) x + ( B A A B ) u y = ( C CA A) x + ( D CA B) u Other choices x, x = Ax + Bu y = Cx + Du x α x x = ( A A ( αi A ) A ) x + ( B A ( αi A ) B ) u, Comparison Direct transmission term (non-zero even if D=) y = ( C C ( αi A ) A ) x + ( D C ( αi A ) B ) u x Gred () = G() x G ( ) ( ) x red = G αx Gred ( α) = G( α)

10 Example. (Rohr s system) x = x+.338 u y = x ( ) Σ =.9 T = /(s+) Reduced second-order model x = x u y = ( ) x+.44u 9 Gs ( ) = Reduced first-order model ( s+ ) s + 3s+ 9 x =.9686x+.46u.8955s G() s = +.44 y =.46x.37u s s+.353 G( s) = /( s+ )

11 Continued Fraction Approximation Routh array For an asymptotically stable system Bs () Gs () = R() s R() s As () c + c3 + c5 + c c4 c6 + R () s + R() s + R3() s s s s With m coefficients, the approximating transfer function is m-th order. Bm() s Gm() s = (with Rm ) Am () s Calculation of coefficients ( ) ( ) As () = a + as + = P () s Bs () = b + bs + = P () s a a a c a/ b b b b c b/ p () () () p p p c 3 p / p = () i ( i ) ( i ) ( i ) ( i ) p k = p k p p k / p + + () i () i () i i i p c i+ p / p p p

12 Interpretation of continued fraction approximation The model reduction presuppose that the inner most TF block may be eliminated. R () m s The condition for good approximation: c / s R () s m m s= For discrete-time system Y( z) = c z U ( z) + Y ( z) i i i+ i+ U ( z) = c Y( z) + U ( z) i+ i i i Yi+ ( z) + c Y( ) i i i i z c z c z Ui+ ( z) = c ( ) i Ui z Yi( z) c ( ) iz Yi+ z Ui( z) = c U i ( z) i + c i c iz + (Used in lattice algorithms) Unimodular (det=)

13 Example Gs () = s + 3s + 59s+ 9 = +.3s+.354s +.44s 3 3 Routh array c c c c.9 c.44 c Reduced-order models G () s = =.5 + s/ s G() s =.5+ /(.7683/ s+ /( s/.473)) /(s+).65s = +.985s +.986s =.5 =.7683 = 4.74 =.47 = 34.9 =.659

14 Moment Matching Transfer function as infinite series H( z ) = h + hz + hz + Matching of the reduced-order model, B m /A m H m m = m m m Am( z ) = + az + + amz M k = nh k n = ( z ) Moment: k n= Moment matching M M M = H m z = dh = dz m B ( z ) b + bz + + bz z= dhm dh m = + dz dz z = z= (,,3, )

15 Example.4 (Moment matching).z 3 H( z ) =.z.54z.54z.7z.8z = Reduced-order model Bm ( z ) b bz Hm ( z ) = Am ( z ) = + + az Moment matching b + b Hm = = z = + a m z= ( + a) dh ab + b = = 4.99 dz dh ( ab ab) = = dz m 3 ( + a ) z = Reduced first-order model H ( z ) = z.7993z

16 Pade Approximation Truncated Taylor series expansion Bs () Gs () = Gm() s = g+ gs + + gm s As () Pade approximation B () s = G () sa () s m m m Example.5 First order: m Ys () 458 Gs ().6s.876s.898s Us () s + 3s + 59s+ 9 3 = = B ()/ s A() s = b /( + as) b/( + as ) =.6 s B()/ s A() s = /(+.3 s) Second order: B ()/ s A () s = ( b + bs)/( + as+ as ) ( b + bs)/( + as+ as ) =.6s+.876s.898s 3 B s A s = s + s+ s ()/ () (.645 )/( )

17 Example.6 (Unstable Pade approximation) B () s (+ ) s A () s ( + s+ s ) b ( ) 3 = = + s s + s + ( g gs ) b g ( ag g) s ags + as = + + = B() s b b = g = B() s = = A() s ( + as ) a = g/ g = A() s ( s) Unstable The Pade approximation can be very poor.

18 Describing Function Analysis Analysis of Linear systems Laplace and Fourier transforms Describing function analysis Extension of frequency analysis to nonlinear systems Based on harmonic analysis Assumptions Assumed periodic solution is sufficiently close to sinusoidal oscillation y(t)=csin(ωt). The nonlinearity can be represented by u=n(x,x ). Fourier series expansion f ( x ikx ) = cos( ) sin( ) k k k k k k a + a kx + a kx ce = = = = / / where a T T k = ()cos( ), k ()sin( ) T f t kt dt b = f t kt dt T / T T/ ( +ϕ ) c = a + b, ϕ = arctan( a / b ) k k k k k k k

19 Description of periodic oscillation with input amplitude C c( C) iϕ b + ia NC ( ) = e = Describing function C C T / where ak = nc ( sin tc, cos t)cos( kt) dt, T ω ω ω T / T / bk = nc ( sin ωtc, ωcos ωt)sin( kt) dt T T / Describing function N(C) is obtained as the amplitude dependent gain N(C) and its phase shift ϕ (ω) for the nonlinear elements. Harmonics caused by nonlinearities are ignored in describing function analysis Characteristic equation Bs () G () s = Apxt ( ) () + Bpnxx ( ) (, ) =, p d/ dt As () ( + G ( pnc ) ( )) xt () = G ( p) = / NC ( ) If the Nyquist curves of G (iω) and /N(C) cross, it indicates the possible existence of a limit cycle with amplitude C and frequency ω.

20 Example.7 (A rate-limited servo) (tan S S / ), + S C C> S NC ( ) = π C C, C S Saturating amplifier For small C, /N(C) (max.) Phase angle of /N(C) is π G (.44 i ) =.667 The Nyquist contour and the describing function cross when K>6. This implies that the limit cycle appears when K>6. K= K=6 K=

21 Condition of application Assumption: the input to nonlinearity is a sinusoid Linear element is usually low-pass: higher harmonics will be attenuated The conditions G ( ikω) G ( iω) for k =,3, and G ( ikω) as k G (s) must not have any imaginary poles s=ik. The function n(x,x ) should have finite partial derivatives w.r.t. x and x and must not be an explicit function of time. The zero-order coefficient c = of n(x,x ). If the prerequisites are not satisfied, the describing function analysis may predict oscillations that do not exist and may fail to predict periodic solutions that indeed do exist.

22 Example.8 (Pulse-width modulation, PWM) PWM leads to asymmetric inputs It requires more complicated analysis with a nonzero constant term c of the describing function. Consider a PWM system with switching nonlinearity The pulse-width modulating signal z is often a sawtooth-shaped signal of high frequency. The pulse-width modulated signal is a square wave with nonzero mean. ut () = nvv (, ) (/ z) xt () + nz If the modulating frequency ω z is high compare to G (s), the high frequency components deriving from z are efficiently absorbed by the low-pass link G (s). ut () (/ z )() xt (linearizedform) (n is the ordinary describing function terms)

23 Balanced Model Reduction in Identification Usage of balanced model reduction Further reduction from a linearized model Identify a reduced-order linear model Identification assisted with model reduction Estimation of high-order model (with sufficient excitation) Model reduction with elimination of less important states Example.9 (Model reduction in identification) S: yk+ =.9yk+. uk + d ( d = ) Identification as first-order model: biased yk+ =.9838yk +.638uk Estimate -th order model, then reduce to first-order model xk+ =.9466xk +.9uk.9z +.9 Y( z) = U( z) y =.9x.9u z.9466 k k k

24 Example. (An impulse response test) Impulse response data t t gt () =. (.5) +. (.75) State-space model for impulse response, g xk+ = xk + uk Balanced model reduction with n= t gt () =.38 (.53) +.38 (.73) t ( ) y = g g g g x k n n k Impulse response coefficients

Lecture Note #6 (Chap.10)

Lecture Note #6 (Chap.10) System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryoo Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of

More information

Analysis of Discrete-Time Systems

Analysis of Discrete-Time Systems TU Berlin Discrete-Time Control Systems 1 Analysis of Discrete-Time Systems Overview Stability Sensitivity and Robustness Controllability, Reachability, Observability, and Detectabiliy TU Berlin Discrete-Time

More information

CBE507 LECTURE III Controller Design Using State-space Methods. Professor Dae Ryook Yang

CBE507 LECTURE III Controller Design Using State-space Methods. Professor Dae Ryook Yang CBE507 LECTURE III Controller Design Using State-space Methods Professor Dae Ryook Yang Fall 2013 Dept. of Chemical and Biological Engineering Korea University Korea University III -1 Overview States What

More information

Analysis of Discrete-Time Systems

Analysis of Discrete-Time Systems TU Berlin Discrete-Time Control Systems TU Berlin Discrete-Time Control Systems 2 Stability Definitions We define stability first with respect to changes in the initial conditions Analysis of Discrete-Time

More information

GATE EE Topic wise Questions SIGNALS & SYSTEMS

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

D(s) G(s) A control system design definition

D(s) G(s) A control system design definition R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure

More information

Control Design. Lecture 9: State Feedback and Observers. Two Classes of Control Problems. State Feedback: Problem Formulation

Control Design. Lecture 9: State Feedback and Observers. Two Classes of Control Problems. State Feedback: Problem Formulation Lecture 9: State Feedback and s [IFAC PB Ch 9] State Feedback s Disturbance Estimation & Integral Action Control Design Many factors to consider, for example: Attenuation of load disturbances Reduction

More information

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

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

More information

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

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

Model Reduction for Unstable Systems

Model Reduction for Unstable Systems Model Reduction for Unstable Systems Klajdi Sinani Virginia Tech klajdi@vt.edu Advisor: Serkan Gugercin October 22, 2015 (VT) SIAM October 22, 2015 1 / 26 Overview 1 Introduction 2 Interpolatory Model

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

4F3 - Predictive Control

4F3 - Predictive Control 4F3 Predictive Control - Lecture 3 p 1/21 4F3 - Predictive Control Lecture 3 - Predictive Control with Constraints Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 3 p 2/21 Constraints on

More information

( ) ( = ) = ( ) ( ) ( )

( ) ( = ) = ( ) ( ) ( ) ( ) Vρ C st s T t 0 wc Ti s T s Q s (8) K T ( s) Q ( s) + Ti ( s) (0) τs+ τs+ V ρ K and τ wc w T (s)g (s)q (s) + G (s)t(s) i G and G are transfer functions and independent of the inputs, Q and T i. Note

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

Robust Multivariable Control

Robust Multivariable Control Lecture 2 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Today s topics Today s topics Norms Today s topics Norms Representation of dynamic systems Today s topics Norms

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

System Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang

System Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang System Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang 1-1 Course Description Emphases Delivering concepts and Practice Programming Identification Methods using Matlab Class

More information

Control Systems. Laplace domain analysis

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

More information

Automatic Control 2. Model reduction. Prof. Alberto Bemporad. University of Trento. Academic year

Automatic Control 2. Model reduction. Prof. Alberto Bemporad. University of Trento. Academic year Lecture: Automatic Control 2 Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 17 Lecture:

More information

STABILITY ANALYSIS TECHNIQUES

STABILITY ANALYSIS TECHNIQUES ECE4540/5540: Digital Control Systems 4 1 STABILITY ANALYSIS TECHNIQUES 41: Bilinear transformation Three main aspects to control-system design: 1 Stability, 2 Steady-state response, 3 Transient response

More information

Model reduction for linear systems by balancing

Model reduction for linear systems by balancing Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,

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

4F3 - Predictive Control

4F3 - Predictive Control 4F3 Predictive Control - Discrete-time systems p. 1/30 4F3 - Predictive Control Discrete-time State Space Control Theory For reference only Jan Maciejowski jmm@eng.cam.ac.uk 4F3 Predictive Control - Discrete-time

More information

ECEEN 5448 Fall 2011 Homework #5 Solutions

ECEEN 5448 Fall 2011 Homework #5 Solutions ECEEN 5448 Fall 211 Homework #5 Solutions Professor David G. Meyer December 8, 211 1. Consider the 1-dimensional time-varying linear system ẋ t (u x) (a) Find the state-transition matrix, Φ(t, τ). Here

More information

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

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

More information

4F3 - Predictive Control

4F3 - Predictive Control 4F3 Predictive Control - Lecture 2 p 1/23 4F3 - Predictive Control Lecture 2 - Unconstrained Predictive Control Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 2 p 2/23 References Predictive

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

Ingegneria dell Automazione - Sistemi in Tempo Reale p.1/28

Ingegneria dell Automazione - Sistemi in Tempo Reale p.1/28 Ingegneria dell Automazione - Sistemi in Tempo Reale Selected topics on discrete-time and sampled-data systems Luigi Palopoli palopoli@sssup.it - Tel. 050/883444 Ingegneria dell Automazione - Sistemi in

More information

6.241 Dynamic Systems and Control

6.241 Dynamic Systems and Control 6.241 Dynamic Systems and Control Lecture 22: Balanced Realization Readings: DDV, Chapter 26 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology April 27, 2011 E. Frazzoli

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

The first order quasi-linear PDEs

The first order quasi-linear PDEs Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.

More information

Module 03 Linear Systems Theory: Necessary Background

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

3 Gramians and Balanced Realizations

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

Need for transformation?

Need for transformation? Z-TRANSFORM In today s class Z-transform Unilateral Z-transform Bilateral Z-transform Region of Convergence Inverse Z-transform Power Series method Partial Fraction method Solution of difference equations

More information

Stabilization and Passivity-Based Control

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

More information

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

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

Transfer Functions. Chapter Introduction. 6.2 The Transfer Function

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

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

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

More information

Chemical Process Dynamics and Control. Aisha Osman Mohamed Ahmed Department of Chemical Engineering Faculty of Engineering, Red Sea University

Chemical Process Dynamics and Control. Aisha Osman Mohamed Ahmed Department of Chemical Engineering Faculty of Engineering, Red Sea University Chemical Process Dynamics and Control Aisha Osman Mohamed Ahmed Department of Chemical Engineering Faculty of Engineering, Red Sea University 1 Chapter 4 System Stability 2 Chapter Objectives End of this

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

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

Controls Problems for Qualifying Exam - Spring 2014

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

Prüfung Regelungstechnik I (Control Systems I) Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Prüfung Regelungstechnik I (Control Systems I) Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 29. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid

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

Module 07 Controllability and Controller Design of Dynamical LTI Systems

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

ẋ 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

Discrete and continuous dynamic systems

Discrete 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

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

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

Time Response of Systems

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

More information

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

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

More information

Introduction to Nonlinear Control Lecture # 4 Passivity

Introduction to Nonlinear Control Lecture # 4 Passivity p. 1/6 Introduction to Nonlinear Control Lecture # 4 Passivity È p. 2/6 Memoryless Functions ¹ y È Ý Ù È È È È u (b) µ power inflow = uy Resistor is passive if uy 0 p. 3/6 y y y u u u (a) (b) (c) Passive

More information

STABILITY OF CLOSED-LOOP CONTOL SYSTEMS

STABILITY OF CLOSED-LOOP CONTOL SYSTEMS CHBE320 LECTURE X STABILITY OF CLOSED-LOOP CONTOL SYSTEMS Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 10-1 Road Map of the Lecture X Stability of closed-loop control

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

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

Conformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.

Conformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder. Lent 29 COMPLEX METHODS G. Taylor A star means optional and not necessarily harder. Conformal maps. (i) Let f(z) = az + b, with ad bc. Where in C is f conformal? cz + d (ii) Let f(z) = z +. What are the

More information

Linear second-order differential equations with constant coefficients and nonzero right-hand side

Linear second-order differential equations with constant coefficients and nonzero right-hand side Linear second-order differential equations with constant coefficients and nonzero right-hand side We return to the damped, driven simple harmonic oscillator d 2 y dy + 2b dt2 dt + ω2 0y = F sin ωt We note

More information

Practice Problems For Test 3

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

Advanced Process Control Tutorial Problem Set 2 Development of Control Relevant Models through System Identification

Advanced Process Control Tutorial Problem Set 2 Development of Control Relevant Models through System Identification Advanced Process Control Tutorial Problem Set 2 Development of Control Relevant Models through System Identification 1. Consider the time series x(k) = β 1 + β 2 k + w(k) where β 1 and β 2 are known constants

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

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

Lecture 3 - Design of Digital Filters

Lecture 3 - Design of Digital Filters Lecture 3 - Design of Digital Filters 3.1 Simple filters In the previous lecture we considered the polynomial fit as a case example of designing a smoothing filter. The approximation to an ideal LPF can

More information

Practice Problems For Test 3

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

Lecture Note #7 (Chap.11)

Lecture Note #7 (Chap.11) System Modeling and Identification Lecture Note #7 (Chap.) CBE 702 Korea University Prof. Dae Ryoo Yang Chap. Real-time Identification Real-time identification Supervision and tracing of time varying parameters

More information

EECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators.

EECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators. Name: SID: EECS C28/ ME C34 Final Thu. May 4, 25 5-8 pm Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 4 2 4 3 6 4 8 5 3

More information

Modelling and Control of Dynamic Systems. Stability of Linear Systems. Sven Laur University of Tartu

Modelling and Control of Dynamic Systems. Stability of Linear Systems. Sven Laur University of Tartu Modelling and Control of Dynamic Systems Stability of Linear Systems Sven Laur University of Tartu Motivating Example Naive open-loop control r[k] Controller Ĉ[z] u[k] ε 1 [k] System Ĝ[z] y[k] ε 2 [k]

More information

Automatique. A. Hably 1. Commande d un robot mobile. Automatique. A.Hably. Digital implementation

Automatique. A. Hably 1. Commande d un robot mobile. Automatique. A.Hably. Digital implementation A. Hably 1 1 Gipsa-lab, Grenoble-INP ahmad.hably@grenoble-inp.fr Commande d un robot mobile (Gipsa-lab (DA)) ASI 1 / 25 Outline 1 2 (Gipsa-lab (DA)) ASI 2 / 25 of controllers Signals must be sampled and

More information

Sampling of Linear Systems

Sampling of Linear Systems Sampling of Linear Systems Real-Time Systems, Lecture 6 Karl-Erik Årzén January 26, 217 Lund University, Department of Automatic Control Lecture 6: Sampling of Linear Systems [IFAC PB Ch. 1, Ch. 2, and

More information

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries

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

Math 211. Substitute Lecture. November 20, 2000

Math 211. Substitute Lecture. November 20, 2000 1 Math 211 Substitute Lecture November 20, 2000 2 Solutions to y + py + qy =0. Look for exponential solutions y(t) =e λt. Characteristic equation: λ 2 + pλ + q =0. Characteristic polynomial: λ 2 + pλ +

More information

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 9. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -

More information

Control Systems Lab - SC4070 Control techniques

Control Systems Lab - SC4070 Control techniques Control Systems Lab - SC4070 Control techniques Dr. Manuel Mazo Jr. Delft Center for Systems and Control (TU Delft) m.mazo@tudelft.nl Tel.:015-2788131 TU Delft, February 16, 2015 (slides modified from

More information

Kernel Method: Data Analysis with Positive Definite Kernels

Kernel Method: Data Analysis with Positive Definite Kernels Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University

More information

1 (2n)! (-1)n (θ) 2n

1 (2n)! (-1)n (θ) 2n Complex Numbers and Algebra The real numbers are complete for the operations addition, subtraction, multiplication, and division, or more suggestively, for the operations of addition and multiplication

More information

Properties of Open-Loop Controllers

Properties of Open-Loop Controllers Properties of Open-Loop Controllers Sven Laur University of Tarty 1 Basics of Open-Loop Controller Design Two most common tasks in controller design is regulation and signal tracking. Regulating controllers

More information

Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0

Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0 Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can

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

MA 266 Review Topics - Exam # 2 (updated)

MA 266 Review Topics - Exam # 2 (updated) MA 66 Reiew Topics - Exam # updated Spring First Order Differential Equations Separable, st Order Linear, Homogeneous, Exact Second Order Linear Homogeneous with Equations Constant Coefficients The differential

More information

MCE693/793: Analysis and Control of Nonlinear Systems

MCE693/793: Analysis and Control of Nonlinear Systems MCE693/793: Analysis and Control of Nonlinear Systems Introduction to Describing Functions Hanz Richter Mechanical Engineering Department Cleveland State University Introduction Frequency domain methods

More information

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3.. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -

More information

Frequency Response of Linear Time Invariant Systems

Frequency Response of Linear Time Invariant Systems ME 328, Spring 203, Prof. Rajamani, University of Minnesota Frequency Response of Linear Time Invariant Systems Complex Numbers: Recall that every complex number has a magnitude and a phase. Example: z

More information

Course roadmap. ME451: Control Systems. Example of Laplace transform. Lecture 2 Laplace transform. Laplace transform

Course roadmap. ME451: Control Systems. Example of Laplace transform. Lecture 2 Laplace transform. Laplace transform ME45: Control Systems Lecture 2 Prof. Jongeun Choi Department of Mechanical Engineering Michigan State University Modeling Transfer function Models for systems electrical mechanical electromechanical Block

More information

Exam in Systems Engineering/Process Control

Exam in Systems Engineering/Process Control Department of AUTOMATIC CONTROL Exam in Systems Engineering/Process Control 27-6-2 Points and grading All answers must include a clear motivation. Answers may be given in English or Swedish. The total

More information

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

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

More information

School of Mechanical Engineering Purdue University

School of Mechanical Engineering Purdue University Case Study ME375 Frequency Response - 1 Case Study SUPPORT POWER WIRE DROPPERS Electric train derives power through a pantograph, which contacts the power wire, which is suspended from a catenary. During

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

Chapter 1 Fundamental Concepts

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

CME 345: MODEL REDUCTION

CME 345: MODEL REDUCTION CME 345: MODEL REDUCTION Balanced Truncation Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu These slides are based on the recommended textbook: A.C. Antoulas, Approximation of

More information

Equilibrium points: continuous-time systems

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

More information

Iterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations

Iterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations Iterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations Klajdi Sinani Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University

More information

MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions

MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions Question 1. SIGNALS: Design of a noise-cancelling headphone system. 1a. Based on the low-pass filter given, design a high-pass filter,

More information

ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky

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

1+t 2 (l) y = 2xy 3 (m) x = 2tx + 1 (n) x = 2tx + t (o) y = 1 + y (p) y = ty (q) y =

1+t 2 (l) y = 2xy 3 (m) x = 2tx + 1 (n) x = 2tx + t (o) y = 1 + y (p) y = ty (q) y = DIFFERENTIAL EQUATIONS. Solved exercises.. Find the set of all solutions of the following first order differential equations: (a) x = t (b) y = xy (c) x = x (d) x = (e) x = t (f) x = x t (g) x = x log

More information

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

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

More information

Autonomous system = system without inputs

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