Lecture Note #6 (Chap.10)
|
|
- Franklin Garrett
- 6 years ago
- Views:
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)
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 informationAnalysis 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 informationCBE507 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 informationAnalysis 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 informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
More informationEEE582 Homework Problems
EEE582 Homework Problems HW. Write a state-space realization of the linearized model for the cruise control system around speeds v = 4 (Section.3, http://tsakalis.faculty.asu.edu/notes/models.pdf). Use
More informationD(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 informationControl 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 informationSome 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 informationTime 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 informationCHBE320 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 informationModel 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 informationECE504: 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 information4F3 - 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 informationLecture 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 informationRobust 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 informationState will have dimension 5. One possible choice is given by y and its derivatives up to y (4)
A Exercise State will have dimension 5. One possible choice is given by y and its derivatives up to y (4 x T (t [ y(t y ( (t y (2 (t y (3 (t y (4 (t ] T With this choice we obtain A B C [ ] D 2 3 4 To
More informationSystem 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 informationControl Systems. Laplace domain analysis
Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output
More informationAutomatic 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 informationSTABILITY 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 informationModel 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 informationIntro. Computer Control Systems: F8
Intro. Computer Control Systems: F8 Properties of state-space descriptions and feedback Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 22 dave.zachariah@it.uu.se F7: Quiz! 2
More information4F3 - 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 informationECEEN 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 informationOutput 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 information4F3 - 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 informationBasic 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 informationIngegneria 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 information6.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 informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response
.. AERO 422: Active Controls for Aerospace Vehicles Dynamic Response Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. . Previous Class...........
More informationThe 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 informationModule 03 Linear Systems Theory: Necessary Background
Module 03 Linear Systems Theory: Necessary Background Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html September
More information3 Gramians and Balanced Realizations
3 Gramians and Balanced Realizations In this lecture, we use an optimization approach to find suitable realizations for truncation and singular perturbation of G. It turns out that the recommended realizations
More informationNeed 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 informationStabilization and Passivity-Based Control
DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationChap 4. State-Space Solutions and
Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations
More informationTransfer Functions. Chapter Introduction. 6.2 The Transfer Function
Chapter 6 Transfer Functions As a matter of idle curiosity, I once counted to find out what the order of the set of equations in an amplifier I had just designed would have been, if I had worked with the
More informationFall 線性系統 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 informationChemical 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 informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from
More informationECE557 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 informationControls Problems for Qualifying Exam - Spring 2014
Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function
More informationPrü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 informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationModule 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 informationKalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and
Kalman Decomposition Controllable / uncontrollable decomposition Suppose that the controllability matrix C R n n of a system has rank n 1
More informationẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty
More informationEE 380. Linear Control Systems. Lecture 10
EE 380 Linear Control Systems Lecture 10 Professor Jeffrey Schiano Department of Electrical Engineering Lecture 10. 1 Lecture 10 Topics Stability Definitions Methods for Determining Stability Lecture 10.
More informationControl Systems. Frequency domain analysis. L. Lanari
Control Systems m i l e r p r a in r e v y n is o Frequency domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic
More informationTime Response of Systems
Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =
More informationControl Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli
Control Systems I Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 13, 2017 E. Frazzoli (ETH)
More informationIntroduction 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 informationSTABILITY 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 informationControl 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 informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
More informationConformal 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 informationLinear 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 informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationAdvanced 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 informationControl 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 informationRobust Control 2 Controllability, Observability & Transfer Functions
Robust Control 2 Controllability, Observability & Transfer Functions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /26/24 Outline Reachable Controllability Distinguishable
More informationLecture 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 informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationLecture 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 informationEECS 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 informationModelling 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 informationAutomatique. 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 informationSampling 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 informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries
. AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace
More informationMath 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!
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 informationControl 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 informationKernel 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 information1 (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 informationProperties 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 informationPower 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 information9.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 informationMA 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 informationMCE693/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!
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 informationFrequency 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 informationCourse 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 informationExam 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 informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationSchool 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 informationUnit 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 informationChapter 1 Fundamental Concepts
Chapter 1 Fundamental Concepts 1 Signals A signal is a pattern of variation of a physical quantity, often as a function of time (but also space, distance, position, etc). These quantities are usually the
More informationCME 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 informationEquilibrium 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 informationIterative 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 informationMAE143a: 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 informationELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky
ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky Equilibrium points and linearization Eigenvalue decomposition and modal form State transition matrix and matrix exponential Stability ELEC 3035 (Part
More informationMODELING 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 information1+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 informationLecture 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 informationAutonomous system = system without inputs
Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation
More information