STATE-SPACE DYNAMIC SYSTEMS (DISCRETE-TIME)
|
|
- Gabriel Willis
- 5 years ago
- Views:
Transcription
1 ECE452/552: Multivariable Control Systems I. 3 1 STATE-SPACE DYNAMIC SYSTEMS (DISCRETE-TIME) 3.1: The transform Computer control requires analog-to-digital (A2D) and digital-toanalog (D2A) conversion. Everything shown in the dashed box below represents the internals of acomputerorembeddedsysteminadigitalcontrolscenario. r.t/ e.t/ A2D Q eœk D. / uœk D2A zoh u.t/ G.s/ w.t/ y.t/ v.t/ The A2D comprises a sampler and quantizer; the D2A performs a zero-order-hold operation. The controller is designed in discrete time, where the transform is used as a design tool rather than the Laplace transform. DEFINITION: The transform is defined for a signal x.kt/ to be X. / D x./ C x.t/ 1 C x.2t/ 2 C x.3t/ 3 C 1X D x.kt/ k kd! 1X or; X. /D xœk k : kd
2 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 2 As with the Laplace-transform variable s, the -transform variable is a complex number. Iassumethatyouhavestudiedthe transform before, and so include only a short review in these notes. EXAMPLE: Define the digital impulse (pulse) function to be ıœk D ( 1; k D I ; otherwise: Note: This is very different from an analog impulse (e.g., ıœ is defined), but plays a similar role. so; let xœk D ıœk 1X X. / D xœk k kd D xœ C xœ1 1 C xœ2 2 C D xœ D 1: This sum converges regardless of the value of, sorocd fj j >g. Selected important properties of the transform Some select properties of the transform may be useful. LINEARITY: The transform is linear. If xœk X. / and vœk V. /,then axœk C bvœk ax. / C bv. /:
3 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 3 RIGHT SHIFT IN TIME (DELAY): When shifting a signal in time, xœk 1 1 X. / C xœ 1 xœk 2 2 X. / C xœ 2 C 1 xœ 1 : xœk q q X. / C xœ q C 1 xœ q C 1 CC qc1 xœ 1 : CONVOLUTION: Discrete-time convolution is analogous to continuoustime convolution:.x v/œk D or; D kx xœi vœk i k id kx vœi xœk i : k : Then, assuming that xœk and vœk are zero for negative k, id.x v/œk X. /V. /: INITIAL-VALUE THEOREM: If xœk X. /, xœ D lim!1 X. /: FINAL-VALUE THEOREM: If a finite final value exists (all poles of X. / strictly inside the unit circle except perhaps for a single pole on the unit circle at D 1), lim k!1 xœk D lim. 1/X. /:!1 The table on the next page lists transforms for commonly encountered sampled time-domain signals:
4 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 4 Example table showing transform of sampled time-domain signals: Time function Laplace transform -Transform e.t/ E.s/ E. / 1.t/ t1.t/ 1 s t t/ 1 s 3 t k 1 1.t/ e at 1.t/ te at 1.t/ t k e at 1.t/ 1 1 s 2 T. 1/ 2.k 1/Š s k 1 sca 1.sCa/ 2.k 1/Š.sCa/ k 1 e at a s.sca/ t 1 e at a a s 2.sCa/ T 2. C1/ 2. 1/ 3 k 1 lim a!. 1/. 1/ e at h k 1 e at T e at. e at / h 2 k e at.1 e at /. 1/. e at / Œ.aT 1Ce at / C.a e at ate at / a. 1/ 2. e at / at a2 1.1 C at/e ate at s.sca/ 2 1 e at. e at / 2 e at e bt sin.at/ cos.at/ e at sin.bt/ e at cos.bt/ b a.sca/.scb/ 1 e at cos.bt/ C a b sin.bt/ a 2 Cb 2 sœ.sca/ 2 Cb 2.e at e bt /. e at /. e bt / a sin.at / s 2 Ca cos.at /C1 s. cos.at // s 2 Ca cos.at /C1 b e at sin.bt /.sca/ 2 Cb e at cos.bt /Ce 2aT sca 2 e at cos.bt /.sca/ 2 Cb e at cos.bt /Ce 2aT.A CB/. 1/. 2 2 e at cos.bt /Ce 2aT / A D 1 e at cos.bt / C a sin.bt / b 1 ab C e at a.a b/ C e bt b.b a/ 1 s.sca/.scb/ B D e 2aT C e at a sin.bt / cos.bt / b.a CB/. e at /. e bt /. 1/ A D b.1 e at / a.1 e bt / ab.b a/ B D ae at.1 e bt / be bt.1 e at / ab.b a/
5 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME : Working with the transform The following shows some correspondence between the -plane and some discrete-time unit-pulse-response signals. Unit step: hœk D 1Œk. H. / D 1 ; j j > Exponential (geometric): hœk D a k 1Œk ; jaj <1. H. / D a ; j j > jaj General cosinusoid: hœk D a k cosœ!k 1Œk ; jaj <1.. a cos!/ H. / D 2 2a.cos!/ C a 2 for j j > jaj. The radius to the two poles is a; theangletothepolesis! The zero (not at the origin) has the same real part as the two poles. If! D, H. / D a :::geometric! If! D ; a D 1, H. / D 1 :::step! Pole radius a is the geometric factor, determines settling time. 1. jaj D; finite-duration response. e.g., ıœk N N. 2. jaj >1; growing signal that will not decay. 3. jaj D1; signal with constant amplitude; either step or cosine. 4. jaj <1; decaying signal. Small a =fastdecay(seebelow).
6 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 6 a duration N Pole angle! determines number of samples per oscillation. That is, if we require cosœ!k D cosœ!.k C N/,then N D 2 ˇ D 36 ˇ :! ˇrad! N D 3 ˇdeg Solid: cst. damping ratio. Dashed: constant natural frequency! n. N D 2 I. / N D 4 N D 5 N D 8 N D 1 N D 2 Plot to right shows discrete-time unit-pulse responses versus pole locations. R. / Discrete Impulse Responses versus Pole Locations
7 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 7 Correspondence with continuous-time signals Let g.t/ D e at cos.bt/1.t/. Suppose Then, a D :3567=T b D =4 T 9 = ; T D sampling period: gœk D g.kt / D e :3567 k k cos 1Œk 4 k D :7 k cos 1Œk : 4 (This is the cosinusoid example used in the earlier example). G.s/ has poles at s1;2 D a C jb and a jb. G. / has poles at radius e at angle! D bt or at e at jbt. So, 1;2 D e s 1T and e s 2T. In general, poles convert between the s-plane and -plane via D e st. EXAMPLE: Some corresponding pole locations: j T j T s-plane -plane
8 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 8 j!-axis maps to unit circle. Constant damping ratio maps to strange spiral. When considering system response to a step input for controls purposes, the following diagrams may be helpful: Good Good Good Damping Frequency! n Settling time Higher-order systems: Pole moving toward D 1, systemslowsdown. Zero moving toward D 1, overshoot. Pole and zero moving close to each other cancel.
9 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME : Discrete-time state-space form Discrete-time systems can also be represented in state-space form. xœk C 1 D A d xœk C B d uœk yœk D C d xœk C D d uœk : The subscript d isusedheretoemphasizethat,ingeneral,the A, B, C and D matrices are different for discrete-time and continuous-time systems, even if the underlying plant is the same. Iwillusuallydropthe d andexpectyoutointerpretthesystemfrom its context. Formulating from transfer functions Discrete-time dynamics are represented as difference equations. e.g., yœkc3 Ca 1 yœkc2 Ca 2 yœkc1 Ca 3 yœk D b 1 uœkc2 Cb 2 uœkc1 Cb 3 uœk yœk Ca 1 yœk 1 Ca 2 yœk 2 Ca 3 yœk 3 D b 1 uœk 1 Cb 2 uœk 2 Cb 3 uœk 3 : This particular example has transfer function G. / D b 1 2 C b 2 C b 3 3 C a 1 2 C a 2 C a 3 D Y. / U. / : This transfer function may be converted to state-space in a very similar way to continuous-time systems. First, consider the poles: G p. / D 1 D V. / 3 C a 1 2 C a 2 C a 3 U. / vœk C 3 C a 1 vœk C 2 C a 2 vœk C 1 C a 3 vœk D uœk :
10 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 1 Choose current and advanced versions of vœk as state. h i T xœk D vœk C 2 vœk C 1 vœk : Then xœk C 1 D D We now add zeros. vœk C 3 vœk C 2 vœk C a 1 a 2 a vœk C vœk C 1 vœk 7 5 C G. / D b 1 2 C b 2 C b 3 3 C a 1 2 C a 2 C a 3 D Y. / U. / : uœk : Break up transfer function into two parts. V. / contains all of the U. / poles of Y. / U. /.Then, Y. / D b 1 2 C b 2 C b 3 V. /: Or, Then xœk C 1 D yœk D yœk D b 1 vœk C 2 C b 2 vœk C 1 C b 3 vœk : h a 1 a 2 a vœk C i h b 1 b 2 b 3 xœk C vœk C 1 vœk i uœk : 7 5 C uœk
11 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 11 Many discrete-time transfer functions are not strictly proper. Solve by polynomial long division, and setting D equal to the quotient. In MATLAB, [A,B,C,D]=tf2ss(num,den) converts a transfer function form to state-space form (same syntax as for cts.-time, but num and den must have same order, maybe zero padded). As with continuous-time systems, we have a lot of freedom when making state-space models (i.e., inchoosingcomponentsofxœk ). Canonical forms In discrete-time we have the same canonical forms: Controller, observer, controllability, observability, modal, and Jordan. They are derived in the same way, as demonstrated above for the controller form. Ablockdiagramforcontrollerformis: b 1 b2 yœk uœk x 1c x 2c x 3c b 3 a 1 a2 a 3 Time (dynamic) response Homogeneous part First, consider the scalar case xœk C 1 D axœk ; xœ :
12 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 12 Take -transform. X. / D. a/ 1 xœ : Inverse -transform. xœk D a k xœ. Similarly, the full solution (vector case) is Aside: Nilpotent systems xœk D A k xœ : A is nilpotent if some power of n exists such that A n D : A does not just decay to zero, it is exactly zero! This might be a desirable control design! (Why?) You might imagine that all the eigenvalues of A must be zero for this to work. Forced solution The full solution is: xœk D A k xœ C Xk 1 j D A k 1 j BuŒj : ƒ convolution This can be proved by induction from the equation xœk C 1 D AxŒk C BuŒk ; xœ Clearly, if yœk D CxŒk C DuŒk ; yœk D CA k xœ C ƒ initial resp: Xk 1 j D CA k 1 j BuŒk ƒ convolution C DuŒk : ƒ feedthrough
13 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME : More on discrete-time state-space models State-space to transfer function Start with the state equations -transform or and So, xœk C 1 D AxŒk C BuŒk yœk D CxŒk C DuŒk : X. / xœ D AX. / C BU. / Y. / D CX. / C DU. /. I A/X. / D BU. / C xœ X. / D. I A/ 1 BU. / C. I A/ 1 xœ Y. / D ŒC. I A/ 1 B C D ƒ transfer function of system Y. / U. / D C. I A/ 1 B C D Same form as for continuous-time systems. Poles of system are roots of detœ I A D. U. / C C. I A/ ƒ 1 xœ : response to initial conditions Transformation State-space representations are not unique. Selection of state x are quite arbitrary.
14 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 14 Analyze the transformation of xœk C 1 D AxŒk C BuŒk yœk D CxŒk C DuŒk Let xœk D TwŒk,whereT is an invertible (similarity) transformation matrix. wœk C 1 D T 1 AT ƒ AN yœk D Same as for continuous-time. wœk C ƒ T 1 B uœk NB ƒ CT wœk C ƒ D uœk NC ND so; wœk C 1 D N AwŒk C NBuŒk yœk D NCwŒk C NDuŒk : Converting plant dynamics to discrete time Combine the dynamics of the zero-order hold and the plant. uœk ZOH u.t/ A; B; C; D y.t/ The continuous-time dynamics of the plant are: Px.t/ D Ax.t/ C Bu.t/ y.t/ D Cx.t/C Du.t/: Evaluate x.t/ at discrete times. Recall x.t/ D xœk C 1 D x..k C 1/T / D Z t e A.t / Bu./d Z.kC1/T e A..kC1/T / Bu./d
15 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 15 With malice aforethought, break up the integral into two pieces. The first piece will become A d times x.kt/. Thesecondpartwillbecome B d times u.kt /. D D Z kt Z kt D e AT x.kt/ C e A..kC1/T / Bu./d C e AT e A.kT / Bu./d C Z.kC1/T kt Z.kC1/T kt Z.kC1/T kt e A..kC1/T / Bu./d: e A..kC1/T / Bu./d e A..kC1/T / Bu./d In the remaining integral, note that u./ is constant from kt to.k C 1/T; and equal to u.kt /; let D.k C 1/T I D.k C 1/T I d D d. Z T x..k C 1/T / D e AT x.kt/ C e A B d u.kt / Z T or; xœkc 1 D e AT xœk C e A B d uœk : So, we have a discrete-time state-space representation from the continuous-time representation. where A d D e AT ;B d D Similarly, xœk C 1 D A d xœk C B d uœk Z T That is, C d D C ; D d D D. e A B d. yœk D CxŒk C DuŒk : Calculating A d, B d, C d and D d Cd and D d require no calculation since C d D C and D d D D.
16 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 16 Ad is calculated via the matrix exponential A d D e AT.Thisisdifferent from taking the exponential of each element in AT. If MATLAB is handy, you can type in Ad=expm(A *T) If MATLAB is not handy, then you need to work a little harder. Recall from the previous set of notes that e At D L 1 Œ.sI A/ 1.So, e AT D L 1 Œ.sI A/ 1 ˇˇtDT ; which is probably the easiest way to work it out by hand. Or, e AT D I C AT C A2 T 2 C A3 T 3 C 2Š 3Š which is a convergent series so may be approximated with only afew terms. Now we focus on computing Bd.Recallthat if A 1 exists. B d D Z T Z T e A B d D I C A C A 22 2 C ::: B d D IT C A T 2 3 2Š C A2T 3Š C ::: B D A 1.e AT I/B D A 1.A d I/B; So, calculating Bd is easy once we have already calculated A d. Also, in MATLAB, [Ad,Bd]=c2d(A,B,T)
17 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME : Linear time-varying and nonlinear discrete-time systems Linear time-varying discrete-time systems Alineartime-varyingsystemcanbewrittenas xœk C 1 D AŒk xœk C BŒk uœk yœk D CŒk xœk C DŒk uœk : Analysis is somewhat easier than for the continuous-time counterpart. Consider first the homogeneous case xœk C 1 D AŒk xœk ; xœk D x ; k : The solution to this is xœk D ˆŒk; k x where ˆŒk; k D 8 < I; k D k : AŒk 1 AŒk 2 AŒk ; k > k : We see that the state-transition matrix can be computed readily, and does not involve any difficult integrals. Some properties of the state-transition matrix include: ˆŒk C 1; k D AŒk ˆŒk; k ˆŒk ;k D I ˆŒk; s ˆŒs; D ˆŒk; : The solution to the nonhomogeneous case can be shown to be xœk D ˆŒk; k x C Xk 1 Dk ˆŒk; C 1 BŒ uœ
18 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 18 yœk D CŒk ˆŒk;k x C Xk 1 Nonlinear discrete-time systems Dk CŒk ˆŒk; C 1 BŒ uœ C DŒk uœk : The approach to working with nonlinear discrete-time systems is similar to that used for nonlinear continuous-time systems. We can linearize around an equilibrium point or around a solution trajectory. We define a nonlinear state-space form as xœk C 1 D f.xœk ;uœk / Linearizing around an equilibrium point yœk D g.xœk ; uœk /: We can linearize around an equilibrium point if there exists an equilibrium constant solution uœk D u eq ; xœk D x eq ; yœk D y eq : Then, let the actual input signal and initial state be written as uœk D u eq C ıuœk ; xœk D x eq C ıx eq Œk for small ıuœk. Then, we can write actual state and output as xœk D x eq C ıxœk ; yœk D y eq C ıyœk : Following the same process used in the prior chapter, the linearized perturbation system is
19 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 19 ıxœk C 1 D AıxŒk C BıuŒk ıyœk D CıxŒk C DıuŒk ; where the linearized state-space matrices are df.x eq ;u eq / df.x eq ;u eq / A D I B D dx du dg.x eq ;u eq / dg.x eq ;u eq / C D I D D dx du and the overall state and output can be computed as xœk D x eq C ıxœk and yœk D y eq C ıyœk : Linearizing around a solution trajectory Alternately, suppose it is known that u sol Œk ; x sol Œk ; y sol Œk form a time-varying solution to the nonlinear dynamics of the system. This means that x sol Œk C 1 D f.x sol Œk ; u sol Œk / y sol Œk D g.x sol Œk ; u sol Œk /: Then, let general input, state, and output be uœk D u sol Œk C ıuœk xœk D x sol Œk C ıxœk yœk D y sol Œk C ıyœk : Proceeding as before, we find the perturbation system
20 ECE452/552, STATE-SPACE DYNAMIC SYSTEMS DISCRETE-TIME 3 2 df.x sol ;u sol / df.x sol ;u sol / ıxœk C 1 ıxœk C ıuœk dxœk duœk ƒ ƒ AŒk BŒk dg.x sol ;u sol / dg.x sol ;u sol / ıyœk ıxœk C ıuœk : dxœk duœk ƒ ƒ CŒk DŒk The overall state and output can be computed as xœk D x sol C ıxœk and yœk D y sol C ıyœk : Notice that even if the nonlinear system is time-invariant, the linearized system will be time-varying, in general. Where to from here? We have now seen how to model both continuous-time and discrete-time systems in state-space form. There are many commonalities between the two, which we will take advantage of in the remainder of the course. The main difference is in interpreting the result. For example, Eigenvalues of the continuous-time A matrix correspond to s-plane locations, and should be in the left-half s-plane for stability; Eigenvalues of the discrete-time A matrix correspond to -plane locations, and should be in the unit circle for stability. In fact, stability is a very important concept, which must be evaluated and ensured before we talk more about control. This will be our next topic.
OUTPUT-FEEDBACK CONTROL
ECE4520/5520: Multivariable Control Systems I. 1 OUTPUT-FEEDBACK CONTROL.1: Open-loop and closed-loop estimators Open-loop estimators State feedback is impractical since we don t know the state! But, what
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 informationSolution of Linear State-space Systems
Solution of Linear State-space Systems Homogeneous (u=0) LTV systems first Theorem (Peano-Baker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state
More informationModule 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 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 informationDiscrete Systems. Step response and pole locations. Mark Cannon. Hilary Term Lecture
Discrete Systems Mark Cannon Hilary Term 22 - Lecture 4 Step response and pole locations 4 - Review Definition of -transform: U() = Z{u k } = u k k k= Discrete transfer function: Y () U() = G() = Z{g k},
More informationModes and Roots ... mx + bx + kx = 0. (2)
A solution of the form x(t) = ce rt to the homogeneous constant coefficient linear equation a x (n) + (n 1). n a n 1 x + + a 1 x + a 0 x = 0 (1) is called a modal solution and ce rt is called a mode of
More informationTopic # /31 Feedback Control Systems
Topic #7 16.30/31 Feedback Control Systems State-Space Systems What are the basic properties of a state-space model, and how do we analyze these? Time Domain Interpretations System Modes Fall 2010 16.30/31
More 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 informationDiscrete-time linear systems
Automatic Control Discrete-time linear systems Prof. Alberto Bemporad University of Trento Academic year 2-2 Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year 2-2 / 34 Introduction
More informationControl Systems. Time response. L. Lanari
Control Systems Time response L. Lanari outline zero-state solution matrix exponential total response (sum of zero-state and zero-input responses) Dirac impulse impulse response change of coordinates (state)
More 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 informationControl Systems. Time response
Control Systems Time response L. Lanari outline zero-state solution matrix exponential total response (sum of zero-state and zero-input responses) Dirac impulse impulse response change of coordinates (state)
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationModule 09 From s-domain to time-domain From ODEs, TFs to State-Space Modern Control
Module 09 From s-domain to time-domain From ODEs, TFs to State-Space Modern Control Ahmad F. Taha EE 3413: Analysis and Desgin of Control Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/
More information20. The pole diagram and the Laplace transform
95 0. The pole diagram and the Laplace transform When working with the Laplace transform, it is best to think of the variable s in F (s) as ranging over the complex numbers. In the first section below
More informationOBSERVABILITY AND CONTROLLABILITY
ECE45/55: Multivariable Control Systems I. 5 OBSERVABILITY AND CONTROLLABILITY 5.: Continuous-time observability: Where am I? We describe dual ideas called observability and controllability. Both have
More informationLinear System Theory. Wonhee Kim Lecture 1. March 7, 2018
Linear System Theory Wonhee Kim Lecture 1 March 7, 2018 1 / 22 Overview Course Information Prerequisites Course Outline What is Control Engineering? Examples of Control Systems Structure of Control Systems
More 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 informationLinear Algebra and ODEs review
Linear Algebra and ODEs review Ania A Baetica September 9, 015 1 Linear Algebra 11 Eigenvalues and eigenvectors Consider the square matrix A R n n (v, λ are an (eigenvector, eigenvalue pair of matrix A
More informationIt is common to think and write in time domain. creating the mathematical description of the. Continuous systems- using Laplace or s-
It is common to think and write in time domain quantities, but this is not the best thing to do in creating the mathematical description of the system we are dealing with. Continuous systems- using Laplace
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 informationMath 216 Second Midterm 16 November, 2017
Math 216 Second Midterm 16 November, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material
More informationSOLVING DIFFERENTIAL EQUATIONS. Amir Asif. Department of Computer Science and Engineering York University, Toronto, ON M3J1P3
SOLVING DIFFERENTIAL EQUATIONS Amir Asif Department of Computer Science and Engineering York University, Toronto, ON M3J1P3 ABSTRACT This article reviews a direct method for solving linear, constant-coefficient
More information11.2 Basic First-order System Methods
112 Basic First-order System Methods 797 112 Basic First-order System Methods Solving 2 2 Systems It is shown here that any constant linear system u a b = A u, A = c d can be solved by one of the following
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 4884 NOVEMBER 9, 7 Summary This is an introduction to ordinary differential equations We
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 informationLinear dynamical systems with inputs & outputs
EE263 Autumn 215 S. Boyd and S. Lall Linear dynamical systems with inputs & outputs inputs & outputs: interpretations transfer function impulse and step responses examples 1 Inputs & outputs recall continuous-time
More informationDiscrete-Time State-Space Equations. M. Sami Fadali Professor of Electrical Engineering UNR
Discrete-Time State-Space Equations M. Sami Fadali Professor of Electrical Engineering UNR 1 Outline Discrete-time (DT) state equation from solution of continuous-time state equation. Expressions in terms
More information1 Steady State Error (30 pts)
Professor Fearing EECS C28/ME C34 Problem Set Fall 2 Steady State Error (3 pts) Given the following continuous time (CT) system ] ẋ = A x + B u = x + 2 7 ] u(t), y = ] x () a) Given error e(t) = r(t) y(t)
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 information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More informationECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77
1/77 ECEN 605 LINEAR SYSTEMS Lecture 7 Solution of State Equations Solution of State Space Equations Recall from the previous Lecture note, for a system: ẋ(t) = A x(t) + B u(t) y(t) = C x(t) + D u(t),
More informationCDS Solutions to the Midterm Exam
CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2
More 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 informationCourse Background. Loosely speaking, control is the process of getting something to do what you want it to do (or not do, as the case may be).
ECE4520/5520: Multivariable Control Systems I. 1 1 Course Background 1.1: From time to frequency domain Loosely speaking, control is the process of getting something to do what you want it to do (or not
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More informationDynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology.
Dynamic Response Assoc. Prof. Enver Tatlicioglu Department of Electrical & Electronics Engineering Izmir Institute of Technology Chapter 3 Assoc. Prof. Enver Tatlicioglu (EEE@IYTE) EE362 Feedback Control
More informationECE504: Lecture 8. D. Richard Brown III. Worcester Polytechnic Institute. 28-Oct-2008
ECE504: Lecture 8 D. Richard Brown III Worcester Polytechnic Institute 28-Oct-2008 Worcester Polytechnic Institute D. Richard Brown III 28-Oct-2008 1 / 30 Lecture 8 Major Topics ECE504: Lecture 8 We are
More informationControl Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 48824. JANUARY 3, 25 Summary. This is an introduction to ordinary differential equations.
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationStability of Dynamical systems
Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u
More informationLaplace 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 informationAutomatic Control Systems theory overview (discrete time systems)
Automatic Control Systems theory overview (discrete time systems) Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations
More informationDifferential Equations Practice: 2nd Order Linear: Nonhomogeneous Equations: Undetermined Coefficients Page 1
Differential Equations Practice: 2nd Order Linear: Nonhomogeneous Equations: Undetermined Coefficients Page 1 Questions Example (3.5.3) Find a general solution of the differential equation y 2y 3y = 3te
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 informationMTH 464: Computational Linear Algebra
MTH 464: Computational Linear Algebra Lecture Outlines Exam 4 Material Prof. M. Beauregard Department of Mathematics & Statistics Stephen F. Austin State University April 15, 2018 Linear Algebra (MTH 464)
More informationAPPPHYS 217 Tuesday 6 April 2010
APPPHYS 7 Tuesday 6 April Stability and input-output performance: second-order systems Here we present a detailed example to draw connections between today s topics and our prior review of linear algebra
More 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 informationEE451/551: Digital Control. Chapter 3: Modeling of Digital Control Systems
EE451/551: Digital Control Chapter 3: Modeling of Digital Control Systems Common Digital Control Configurations AsnotedinCh1 commondigitalcontrolconfigurations As noted in Ch 1, common digital control
More informationMATH 23 Exam 2 Review Solutions
MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution
More informationSolution of ODEs using Laplace Transforms. Process Dynamics and Control
Solution of ODEs using Laplace Transforms Process Dynamics and Control 1 Linear ODEs For linear ODEs, we can solve without integrating by using Laplace transforms Integrate out time and transform to Laplace
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 informationDepartment of Mathematics IIT Guwahati
Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,
More informationSelected Solutions: 3.5 (Undetermined Coefficients)
Selected Solutions: 3.5 (Undetermined Coefficients) In Exercises 1-10, we want to apply the ideas from the table to specific DEs, and solve for the coefficients as well. If you prefer, you might start
More informationMath 331 Homework Assignment Chapter 7 Page 1 of 9
Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a
More information27. The pole diagram and the Laplace transform
124 27. The pole diagram and the Laplace transform When working with the Laplace transform, it is best to think of the variable s in F (s) as ranging over the complex numbers. In the first section below
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52
1/52 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 2 Laplace Transform I Linear Time Invariant Systems A general LTI system may be described by the linear constant coefficient differential equation: a n d n
More informationOptimal Polynomial Control for Discrete-Time Systems
1 Optimal Polynomial Control for Discrete-Time Systems Prof Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning this paper should
More informationHigher Order Linear Equations
C H A P T E R 4 Higher Order Linear Equations 4.1 1. The differential equation is in standard form. Its coefficients, as well as the function g(t) = t, are continuous everywhere. Hence solutions are valid
More informationCONTROL OF DIGITAL SYSTEMS
AUTOMATIC CONTROL AND SYSTEM THEORY CONTROL OF DIGITAL SYSTEMS Gianluca Palli Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI) Università di Bologna Email: gianluca.palli@unibo.it
More informationINTRODUCTION TO TRANSFER FUNCTIONS
INTRODUCTION TO TRANSFER FUNCTIONS The transfer function is the ratio of the output Laplace Transform to the input Laplace Transform assuming zero initial conditions. Many important characteristics of
More informationwe get y 2 5y = x + e x + C: From the initial condition y(0) = 1, we get 1 5 = 0+1+C; so that C = 5. Completing the square to solve y 2 5y = x + e x 5
Math 24 Final Exam Solution 17 December 1999 1. Find the general solution to the differential equation ty 0 +2y = sin t. Solution: Rewriting the equation in the form (for t 6= 0),we find that y 0 + 2 t
More 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 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 informationCDS 101: Lecture 4.1 Linear Systems
CDS : Lecture 4. Linear Systems Richard M. Murray 8 October 4 Goals: Describe linear system models: properties, eamples, and tools Characterize stability and performance of linear systems in terms of eigenvalues
More informationLinear Systems. Chapter Basic Definitions
Chapter 5 Linear Systems Few physical elements display truly linear characteristics. For example the relation between force on a spring and displacement of the spring is always nonlinear to some degree.
More information16.30 Estimation and Control of Aerospace Systems
16.30 Estimation and Control of Aerospace Systems Topic 5 addendum: Signals and Systems Aeronautics and Astronautics Massachusetts Institute of Technology Fall 2010 (MIT) Topic 5 addendum: Signals, Systems
More informationCourse Outline. Higher Order Poles: Example. Higher Order Poles. Amme 3500 : System Dynamics & Control. State Space Design. 1 G(s) = s(s + 2)(s +10)
Amme 35 : System Dynamics Control State Space Design Course Outline Week Date Content Assignment Notes 1 1 Mar Introduction 2 8 Mar Frequency Domain Modelling 3 15 Mar Transient Performance and the s-plane
More informationC(s) R(s) 1 C(s) C(s) C(s) = s - T. Ts + 1 = 1 s - 1. s + (1 T) Taking the inverse Laplace transform of Equation (5 2), we obtain
analyses of the step response, ramp response, and impulse response of the second-order systems are presented. Section 5 4 discusses the transient-response analysis of higherorder systems. Section 5 5 gives
More informationSolutions to Homework 1, Introduction to Differential Equations, 3450: , Dr. Montero, Spring y(x) = ce 2x + e x
Solutions to Homewor 1, Introduction to Differential Equations, 3450:335-003, Dr. Montero, Spring 2009 problem 2. The problem says that the function yx = ce 2x + e x solves the ODE y + 2y = e x, and ass
More informationCDS 101/110: Lecture 3.1 Linear Systems
CDS /: Lecture 3. Linear Systems Goals for Today: Describe and motivate linear system models: Summarize properties, examples, and tools Joel Burdick (substituting for Richard Murray) jwb@robotics.caltech.edu,
More information22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes
Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration
More informationMethods for analysis and control of. Lecture 6: Introduction to digital control
Methods for analysis and of Lecture 6: to digital O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.inpg.fr www.lag.ensieg.inpg.fr/sename 6th May 2009 Outline Some interesting books:
More informationAdvanced Training Course on FPGA Design and VHDL for Hardware Simulation and Synthesis
065-3 Advanced Training Course on FPGA Design and VHDL for Hardware Simulation and Synthesis 6 October - 0 November, 009 Digital Signal Processing The z-transform Massimiliano Nolich DEEI Facolta' di Ingegneria
More informationChapter 3. State Feedback - Pole Placement. Motivation
Chapter 3 State Feedback - Pole Placement Motivation Whereas classical control theory is based on output feedback, this course mainly deals with control system design by state feedback. This model-based
More informationEE451/551: Digital Control. Chapter 7: State Space Representations
EE45/55: Digital Control Chapter 7: State Spae Representations State Variables Definition 7.: The system state is a minimal set of { variables x ( t ), i =, 2,, n needed together with the i input ut, and
More informationCopyright (c) 2006 Warren Weckesser
2.2. PLANAR LINEAR SYSTEMS 3 2.2. Planar Linear Systems We consider the linear system of two first order differential equations or equivalently, = ax + by (2.7) dy = cx + dy [ d x x = A x, where x =, and
More informationA: Brief Review of Ordinary Differential Equations
A: Brief Review of Ordinary Differential Equations Because of Principle # 1 mentioned in the Opening Remarks section, you should review your notes from your ordinary differential equations (odes) course
More informationTransfer function and linearization
Transfer function and linearization Daniele Carnevale Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata Corso di Controlli Automatici, A.A. 24-25 Testo del corso:
More 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 informationChapter 7. Digital Control Systems
Chapter 7 Digital Control Systems 1 1 Introduction In this chapter, we introduce analysis and design of stability, steady-state error, and transient response for computer-controlled systems. Transfer functions,
More informationThe Laplace transform
The Laplace transform Samy Tindel Purdue University Differential equations - MA 266 Taken from Elementary differential equations by Boyce and DiPrima Samy T. Laplace transform Differential equations 1
More informationMath 3313: Differential Equations Second-order ordinary differential equations
Math 3313: Differential Equations Second-order ordinary differential equations Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Mass-spring & Newton s 2nd law Properties
More informationDIFFERENTIAL EQUATIONS REVIEW. Here are notes to special make-up discussion 35 on November 21, in case you couldn t make it.
DIFFERENTIAL EQUATIONS REVIEW PEYAM RYAN TABRIZIAN Here are notes to special make-up discussion 35 on November 21, in case you couldn t make it. Welcome to the special Friday after-school special of That
More informationEECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators.
Name: SID: EECS C28/ ME C34 Final Wed. Dec. 5, 2 8- am Closed book. Two pages of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 2 2 6 3 4 4 5 6 6 7 8 2 Total
More informationUniversity of Toronto Department of Electrical and Computer Engineering ECE410F Control Systems Problem Set #3 Solutions = Q o = CA.
University of Toronto Department of Electrical and Computer Engineering ECE41F Control Systems Problem Set #3 Solutions 1. The observability matrix is Q o C CA 5 6 3 34. Since det(q o ), the matrix is
More informationSolutions to Math 53 Math 53 Practice Final
Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points
More informationMath 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More informationControl, Stabilization and Numerics for Partial Differential Equations
Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua
More informationME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ
ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and
More 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 informationAPPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015.
APPM 23: Final Exam :3am :pm, May, 25. ON THE FRONT OF YOUR BLUEBOOK write: ) your name, 2) your student ID number, 3) lecture section, 4) your instructor s name, and 5) a grading table for eight questions.
More informationµ = e R p(t)dt where C is an arbitrary constant. In the presence of an initial value condition
MATH 3860 REVIEW FOR FINAL EXAM The final exam will be comprehensive. It will cover materials from the following sections: 1.1-1.3; 2.1-2.2;2.4-2.6;3.1-3.7; 4.1-4.3;6.1-6.6; 7.1; 7.4-7.6; 7.8. The following
More informationMath 216 Final Exam 24 April, 2017
Math 216 Final Exam 24 April, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationESC794: Special Topics: Model Predictive Control
ESC794: Special Topics: Model Predictive Control Discrete-Time Systems Hanz Richter, Professor Mechanical Engineering Department Cleveland State University Discrete-Time vs. Sampled-Data Systems A continuous-time
More informationEE -213 BASIC CIRCUIT ANALYSIS LAB MANUAL
EE -213 BASIC CIRCUIT ANALYSIS LAB MANUAL EE 213 Fall 2009 LABORATORY #1 INTRODUCTION TO MATLAB INTRODUCTION The purpose of this laboratory is to introduce you to Matlab and to illustrate some of its circuit
More information