Control Systems. Time response. L. Lanari


 Rosalind Hicks
 4 years ago
 Views:
Transcription
1 Control Systems Time response L. Lanari
2 outline zerostate solution matrix exponential total response (sum of zerostate and zeroinput responses) Dirac impulse impulse response change of coordinates (state) Lanari: CS  Time response 2
3 system x() = x u(t) ẋ(t) =Ax(t)+Bu(t) y(t) =Cx(t)+Du(t) y(t) input x(t) state output Linear Time Invariant (LTI) dynamical system (Continuous Time) x IR n u IR m y IR p Lanari: CS  Time response 3
4 system representation ẋ(t) =Ax(t)+Bu(t) x() = x implicit representation v = F m v() = v example x(t) =... explicit representation (solution) v(t) =v + 1 m F ( )d we want to study the solution of the set of differential equations in order to have a qualitative knowledge of the system motion we need the general expression of the solution  first we look at the solution when no input is applied (zeroinput response)  then we look at the contribution due to the input (zerostate response) Lanari: CS  Time response 4
5 solution: zeroinput response zeroinput response (i.e. with u(t) = ) ẋ(t) =Ax(t)+Bu(t) scalar case ẋ = ax x() = x x(t) =e at x matrix case ẋ = Ax x() = x x(t) =e At x dimensions n 1 n 1 n n definition e At = I + At + A 2 t2 2! + = 1 X k= A k tk k! check ẋ = d dt e At x = = Ax Lanari: CS  Time response 5
6 matrix exponential definition e At = I + At + A 2 t2 2! + = 1 X k= A k tk k! properties e At = I t= e At1 e At 2 = e A(t 1+t 2 ) e A 1t e A2t = e (A 1+A 2 )t e At 1 = e At consistency composition = () A1A2 = A2A1 equality holds iff d dt eat = Ae At Lanari: CS  Time response 6
7 solution: zeroinput response the exponential matrix propagates the initial condition into state at time t x e At propagation x(t) the exponential matrix propagates the state t seconds forward in time x( ) since x( + t) =e At x( ) x( + t) =e A( +t) x() = e At e A x() = e At x( ) curiosity: Euler approximation exact solution x( + t) x( )+tẋ( )=(I + At)x( ) x( + t) =e At x( )=(I + At + A 2 t 2 /2! +...)x( ) Lanari: CS  Time response 7
8 solution: zeroinput response phase plane x = x1 x 2 it is possible to represent the vector field ẋ = Ax ẋ(t) =Ax(t) x(t) vector field ẋ(t) = apple x(t) and build a phase portrait (geometric representation of the trajectories) Lanari: CS  Time response 8
9 solution on phase plane: pendulum phase portrait Pendulum with no damping (nonlinear differential equation) #(t) solution #(t)+ g` sin #(t) = (#(), #()) = (, ) unstable equilibrium state stable equilibrium state vector field # (#(), #()) = (, ) #(t) #() (#(), #()) = (.45,.8) Lanari: CS  Time response 9
10 solution: total response scalar case ẋ = ax + bu solution is x(t) =e at x + e a(t ) bu( )d (check) Leibniz integral rule d dt ẋ = ae at x + f(t, )d = d dt f(t, )d + f(t, ) =t ae a(t ) bu( )d + bu(t) =ax + bu example: see the point mass v(t) =v + 1 m F ( )d Lanari: CS  Time response 1
11 solution: total response matrix case state x(t) =e At x + e A(t ) Bu( )d zeroinput response (ZIR) zerostate response (ZSR) output y(t) =Ce At x + Ce A(t ) Bu( )d + Du(t) e At Ce At statetransition matrix outputtransition matrix N.B. product is not commutative for matrices Lanari: CS  Time response 11
12 solution: total response total response = zeroinput response + zerostate response x() 6= x() = u(t) = u(t) 6= two distinct contributions to the motion of a linear system a nonzero initial condition causes motion as well as a nonzero input alternative names zeroinput response = free response/evolution zerostate response = forced response/evolution Lanari: CS  Time response 12
13 superposition principle (linearity) from general solution x(t) =e At x + y(t) =Ce At x + e A(t ) Bu( )d Ce A(t ) Bu( )d + Du(t) if (x a,u a (t)) generates x a (t) and y a (t) if (x b,u b (t)) generates x b (t) and y b (t) then ( x a + x b, u a (t)+ u b (t)) generates x a (t)+ x b (t) only for same linear combination and y a (t)+ y b (t) special case ZSR (forced response) x() = u a (t) ZSRa if then u a (t) + u a (t) ZSRa + ZSRb u b (t) ZSRb Lanari: CS  Time response 13
14 superposition principle: special case  example x() = disturbance d(t) control input u(t) ẋ(t) =Ax(t)+B 1 u(t)+b 2 d(t) y(t) =Cx(t)+D 1 u(t)+d 2 d(t) output y(t) total output ZSR (forced response) is made up of two contributions a first due only to the control input (setting the disturbance to zero) a second due only to the disturbance (setting the control input to zero) y(t) =y ci (t)+y d (t) u(t) 6= d(t) = u(t) = d(t) 6= property frequently used in control design Lanari: CS  Time response 14
15 Dirac s delta (impulse) generalized function (t) = if t 6= Z +1 1 properties ( )d =1 f(t) (t ) =f( ) (t ) impulse centered in t = (or centered in = t) 1 " approximation (as a function of t) f( ) (t ) =f(t) (t ) (as a function of ) (t ) if " t "! t Z +1 1 f( ) (t )d = f(t) sifting (sampling) property (as a function of ) Lanari: CS  Time response 15 t f( )Z
16 Dirac s delta: application zerostate output response Ce A(t ) Bu( )d + Du(t) = rewritten as sifting property h i Ce A(t ) B + D (t ) u( )d W (t )u( )d with W (t) =Ce At B + D (t) for now just a more compact way to rewrite the ZSR, but it can be also given an interesting physical interpretation W (t )u( )d if u(t) = ±(t) W(t) defined as the (output) impulse response i.e. the response to a specific input, the Dirac impulse W (t ) ( )d = W (t) Lanari: CS  Time response 16
17 impulse response for any input u(t), the zerostate output response is a convolution integral of the impulse response W(t) with the input u(t) output ZSR = W (t )u( )d to the input u(t) the knowledge of the sole (output) impulse response W(t) allows us to predict the zerostate (output) response to any input u(t) a unique experiment aimed at the determination of the impulse response W(t) is, theoretically, sufficient to characterize any zerostate response Lanari: CS  Time response 17
18 impulse response  example mass + friction (with coefficient µ), measuring velocity v = µ m v + 1 m F x = v u = F ẋ = Ax + Bu A = µ m B = 1 m C = 1 D = impulse response W (t) =Ce At B = 1 m e µ m t F(t) t F a m F W(t) 1 m output impulse response physical interpretation time larger friction coefficient gives faster impulse response decay larger mass gives slower impulse response decay Lanari: CS  Time response 18
19 impulse response  example mass + friction (with coefficient µ), measuring velocity v = µ m v + 1 m F impulse response W (t) =CeAt B = 1 m e µ m t for a different input, for example u(t) = sin!t, we can compute explicitly the zerostate response from the knowledge of the impulse response y(t) = = 1 m W (t )u( )d = µ m 1 2 +! 2 output 1 h µ m sin!t! cos!t +!e µ m ti zerostate response to sin(5t) time Lanari: CS  Time response 19 y(t) u(t)
20 Dirac s delta: application similarly for the state zerostate response define H(t) as H(t) =e At B the zerostate response (of the state) can be written as the following convolution integral interpretation H( ) (t e A(t ) Bu( )d = H(t )u( )d )d = H(t) sifting property Lanari: CS  Time response 2 if u(t) = ±(t) H(t )u( )d H(t H(t) state impulse response e At B t ) ( )d = =! = t = t! =
21 state impulse response H(t) state impulse response e At B other interpretation (= ZIR from x ) H(t) =e At B formally looks like e At x with x = B ẋ = x u 2 state impulse response 1.5 the impulse is transferring the state instantaneously from to B, and then evolves with no input as a free evolution from x = B x x 1 Lanari: CS  Time response 21
22 impulse response impact hammer Impulse Response Amplitude Time (sec) impulse response (here not experimental) this beam is on the space telescope here Hubble Space Telescope Lanari: CS  Time response 22
23 impulse response  experimental determination Modal testing for vibration analysis: the impulse, which has an infinitely small duration, is the ideal testing impact to a structure: all vibration modes will be excited with the same amount of energy (more on this in the frequency analysis section) the impact hammer should be able to replicate this ideal impulse but in reality the strike cannot have an infinitesimal small duration the finite duration of the real impact influences the frequency content of the applied force: the longer is the duration the smaller bandwidth more in the Mechanical Vibrations course Lanari: CS  Time response 23
24 solution: total response more compact notation x(t) = (t)x + y(t) = (t)x + H(t )u( )d W (t )u( )d state (t) =e At H(t) =e At B output (t) =Ce At W (t) =Ce At B + D (t) transition matrix impulse response Dirac impulse Lanari: CS  Time response 24
25 change of coordinates (A, B, C, D) ( A, e B, e C, e D) e state x z = Tx det(t ) 6= state z ẋ = Ax + Bu y = Cx + Du state space representation change of coordinates same system input u & output y do not change, only state is chosen differently ż = e Az + e Bu y = e Cz + e Du the matrices of the two equivalent system representations are related as ea = TAT 1 e B = TB e C = CT 1 e D = D equivalent system representation (proof) Lanari: CS  Time response 25
26 impulse response 1 experiment 1 response ẋ = Ax + Bu y = Cx + Du same system (different representation) ż = e Az + e Bu y = e Cz + e Du W (t) =Ce At B + D (t) f W (t) = e Ce e At e B + e D (t) (proof) W (t) = f W (t) same impulse response i.e. independent from the chosen set of coordinates (state) Lanari: CS  Time response 26
27 general solution (recap) x(t) = (t)x + y(t) = (t)x + H(t )u( )d W (t )u( )d (t) =e At H(t) =e At B (t) =Ce At W (t) =Ce At B + D (t) the matrix exponential appears everywhere do we need to compute the exponential using its definition? e At = I + At + A 2 t2 2! + = 1 X k= A k tk k! Lanari: CS  Time response 27
28 matrix exponential SIAM REVIEW c 23 Society for Industrial and Applied Mathematics Vol. 45, No. 1, pp. 3 Nineteen Dubious Ways to Compute the Exponential of a Matrix, TwentyFive Years Later Cleve Moler Charles Van Loan there are many different ways to compute the matrix exponential... Abstract. In principle, the exponential of a matrix could be computed in many ways. Methods involving approximation theory, differential equations, the matrix eigenvalues, and the matrix characteristic polynomial have been proposed. In practice, consideration of computational stability and efficiency indicates that some of the methods are preferable to others, but that none are completely satisfactory. Most of this paper was originally published in An update, with a separate bibliography, describes a few recent developments. Key words. matrix, exponential, roundoff error, truncation error, condition AMS subject classifications. 15A15, 65F15, 65F3, 65L99 PII. S Introduction. Mathematical models of many physical, biological, and economic processes involve systems of linear, constant coefficient ordinary differential equations ẋ(t) =Ax(t). Here A is a given, fixed, real or complex nbyn matrix. A solution vector x(t) is sought which satisfies an initial condition x() = x. In control theory, A is known as the state companion matrix and x(t) isthesystem response. In principle, the solution is given by x(t) =e ta x where e ta can be formally defined by the convergent power series e ta = I + ta + t2 A ! Published electronically February 3, 23. A portion of this paper originally appeared in SIAM Review, Volume2,Number4,1978,pages The MathWorks, Inc., 3 Apple Hill Drive, Natick, MA Department of Computer Science, Cornell University, 413 Upson Hall, Ithaca, NY Lanari: CS  Time response 28
29 matrix exponential we will use changes of coordinates (state) z = T x ẋ = Ax + Bu y = Cx + Du e At = T 1 e At e T easier to compute if 9 T : z = Tx det(t ) 6= then ż = Az e + Bu e y = Cz e + Du e with e e At easier to compute e At e = e TAT 1t = Te At T 1 (proof) Lanari: CS  Time response 29
30 since it will be used frequently, we prove that if A = T 1 AT e then e At = T 1 e At e T from ea = TAT 1 we obtain A = T 1 AT e e At = e T 1 e AT t = = 1X k= = T 1 " 1 X Lanari: CS  Time response 3 1X k= t k k! T 1 ( e A) k T k= = T 1 e e At T t k k! ( e A) k t k k! (T 1 e AT ) k # T
31 vocabulary English phase plane ZIR/ZSR vector field (state) impulse response convolution integral transition matrix sampling property superposition principle Italiano piano delle fasi evoluzione libera/forzata campo vettoriale risposta impulsiva (nello stato) integrale di convoluzione matrice di transizione proprietà di campionamento principio di sovrapposizione Lanari: CS  Time response 31
Control Systems. Time response
Control Systems Time response L. Lanari outline zerostate solution matrix exponential total response (sum of zerostate and zeroinput responses) Dirac impulse impulse response change of coordinates (state)
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 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 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 informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 Firstorder ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationControl Systems. Internal Stability  LTI systems. L. Lanari
Control Systems Internal Stability  LTI systems L. Lanari outline LTI systems: definitions conditions South stability criterion equilibrium points Nonlinear systems: equilibrium points examples stable
More 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 informationObservability and state estimation
EE263 Autumn 2015 S Boyd and S Lall Observability and state estimation state estimation discretetime observability observability controllability duality observers for noiseless case continuoustime observability
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 informationLecture 19 Observability and state estimation
EE263 Autumn 200708 Stephen Boyd Lecture 19 Observability and state estimation state estimation discretetime observability observability controllability duality observers for noiseless case continuoustime
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 informationSolution of Linear Statespace Systems
Solution of Linear Statespace Systems Homogeneous (u=0) LTV systems first Theorem (PeanoBaker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state
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 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 DMAVT ETH Zürich October 13, 2017 E. Frazzoli (ETH)
More informationMath Ordinary Differential Equations
Math 411  Ordinary Differential Equations Review Notes  1 1  Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationLinear 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 continuoustime
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 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 information06/12/ rws/jMc modif SuFY10 (MPF)  Textbook Section IX 1
IV. ContinuousTime Signals & LTI Systems [p. 3] Analog signal definition [p. 4] Periodic signal [p. 5] Onesided signal [p. 6] Finite length signal [p. 7] Impulse function [p. 9] Sampling property [p.11]
More informationChap 4. StateSpace Solutions and
Chap 4. StateSpace Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear TimeVarying (LTV) Equations
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 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 informationEE Control Systems LECTURE 9
Updated: Sunday, February, 999 EE  Control Systems LECTURE 9 Copyright FL Lewis 998 All rights reserved STABILITY OF LINEAR SYSTEMS We discuss the stability of input/output systems and of statespace
More informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear timeinvariant systems Katalin Hangos University of Pannonia Faculty
More information1 Continuoustime Systems
Observability Completely controllable systems can be restructured by means of state feedback to have many desirable properties. But what if the state is not available for feedback? What if only the output
More informationCyberPhysical Systems Modeling and Simulation of Continuous Systems
CyberPhysical Systems Modeling and Simulation of Continuous Systems Matthias Althoff TU München 29. May 2015 Matthias Althoff Modeling and Simulation of Cont. Systems 29. May 2015 1 / 38 Ordinary Differential
More information1 (30 pts) Dominant Pole
EECS C8/ME C34 Fall Problem Set 9 Solutions (3 pts) Dominant Pole For the following transfer function: Y (s) U(s) = (s + )(s + ) a) Give state space description of the system in parallel form (ẋ = Ax +
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 409. March, 2016 Fundamentals Overview Multiple DOF Systems Statespace Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
More informationChapter 6: The Laplace Transform. ChihWei Liu
Chapter 6: The Laplace Transform ChihWei Liu Outline Introduction The Laplace Transform The Unilateral Laplace Transform Properties of the Unilateral Laplace Transform Inversion of the Unilateral Laplace
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 informationControl Systems I. Lecture 2: Modeling and Linearization. Suggested Readings: Åström & Murray Ch Jacopo Tani
Control Systems I Lecture 2: Modeling and Linearization Suggested Readings: Åström & Murray Ch. 23 Jacopo Tani Institute for Dynamic Systems and Control DMAVT ETH Zürich September 28, 2018 J. Tani, E.
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 informationAnalog Signals and Systems and their properties
Analog Signals and Systems and their properties Main Course Objective: Recall course objectives Understand the fundamentals of systems/signals interaction (know how systems can transform or filter signals)
More 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 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 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 informationFinal Exam Sample Problems, Math 246, Spring 2018
Final Exam Sample Problems, Math 246, Spring 2018 1) Consider the differential equation dy dt = 9 y2 )y 2. a) Find all of its stationary points and classify their stability. b) Sketch its phaseline portrait
More informationLecture 2. Linear Systems
Lecture 2. Linear Systems Ivan Papusha CDS270 2: Mathematical Methods in Control and System Engineering April 6, 2015 1 / 31 Logistics hw1 due this Wed, Apr 8 hw due every Wed in class, or my mailbox on
More informationContinuous Dynamics Solving LTI statespace equations גרא וייס המחלקה למדעי המחשב אוניברסיטת בןגוריון
Continuous Dynamics Solving LTI statespace equations גרא וייס המחלקה למדעי המחשב אוניברסיטת בןגוריון 2 State Space Models For a causal system with m inputs u t R m and p outputs y t R p, an nthorder
More informationSolving Dynamic Equations: The State Transition Matrix
Overview Solving Dynamic Equations: The State Transition Matrix EGR 326 February 24, 2017 Solutions to coupled dynamic equations Solutions to dynamic circuits from EGR 220 The state transition matrix Discrete
More informationControl Systems. Internal Stability  LTI systems. L. Lanari
Control Systems Internal Stability  LTI systems L. Lanari definitions (AS)  A system S is said to be asymptotically stable if its state zeroinput response converges to the origin for any initial condition
More informationDifferential Equations and Modeling
Differential Equations and Modeling Preliminary Lecture Notes Adolfo J. Rumbos c Draft date: March 22, 2018 March 22, 2018 2 Contents 1 Preface 5 2 Introduction to Modeling 7 2.1 Constructing Models.........................
More informationLecture 1 From ContinuousTime to DiscreteTime
Lecture From ContinuousTime to DiscreteTime Outline. Continuous and DiscreteTime Signals and Systems................. What is a signal?................................2 What is a system?.............................
More informationSolutions of Spring 2008 Final Exam
Solutions of Spring 008 Final Exam 1. (a) The isocline for slope 0 is the pair of straight lines y = ±x. The direction field along these lines is flat. The isocline for slope is the hyperbola on the left
More informationLecture Notes of EE 714
Lecture Notes of EE 714 Lecture 1 Motivation Systems theory that we have studied so far deals with the notion of specified input and output spaces. But there are systems which do not have a clear demarcation
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 informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the HodgkinHuxley 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 informationSTATESPACE DYNAMIC SYSTEMS (DISCRETETIME)
ECE452/552: Multivariable Control Systems I. 3 1 STATESPACE DYNAMIC SYSTEMS (DISCRETETIME) 3.1: The transform Computer control requires analogtodigital (A2D) and digitaltoanalog (D2A) conversion.
More informationChapter 3 Convolution Representation
Chapter 3 Convolution Representation DT UnitImpulse Response Consider the DT SISO system: xn [ ] System yn [ ] xn [ ] = δ[ n] If the input signal is and the system has no energy at n = 0, the output yn
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informatione st f (t) dt = e st tf(t) dt = L {t f(t)} s
Additional operational properties How to find the Laplace transform of a function f (t) that is multiplied by a monomial t n, the transform of a special type of integral, and the transform of a periodic
More informationMath 215/255 Final Exam (Dec 2005)
Exam (Dec 2005) Last Student #: First name: Signature: Circle your section #: Burggraf=0, Peterson=02, Khadra=03, Burghelea=04, Li=05 I have read and understood the instructions below: Please sign: Instructions:.
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 5: Calculating the Laplace Transform of a Signal Introduction In this Lecture, you will learn: Laplace Transform of Simple
More informationIntroduction ODEs and Linear Systems
BENG 221 Mathematical Methods in Bioengineering ODEs and Linear Systems Gert Cauwenberghs Department of Bioengineering UC San Diego 1.1 Course Objectives 1. Acquire methods for quantitative analysis and
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 informationEE263: Introduction to Linear Dynamical Systems Review Session 6
EE263: Introduction to Linear Dynamical Systems Review Session 6 Outline diagonalizability eigen decomposition theorem applications (modal forms, asymptotic growth rate) EE263 RS6 1 Diagonalizability consider
More information信號與系統 Signals and Systems
Spring 2015 信號與系統 Signals and Systems Chapter SS2 Linear TimeInvariant Systems FengLi Lian NTUEE Feb15 Jun15 Figures and images used in these lecture notes are adopted from Signals & Systems by Alan
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 informationDiscretetime linear systems
Automatic Control Discretetime linear systems Prof. Alberto Bemporad University of Trento Academic year 22 Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year 22 / 34 Introduction
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More information21 Linear StateSpace Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear StateSpace Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More informationSection 9.8 Higher Order Linear Equations
Section 9.8 Higher Order Linear Equations Key Terms: Higher order linear equations Equivalent linear systems for higher order equations Companion matrix Characteristic polynomial and equation A linear
More informationOrdinary differential equations
Class 11 We will address the following topics Convolution of functions Consider the following question: Suppose that u(t) has Laplace transform U(s), v(t) has Laplace transform V(s), what is the inverse
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 informationSignals and Systems Chapter 2
Signals and Systems Chapter 2 ContinuousTime Systems Prof. Yasser Mostafa Kadah Overview of Chapter 2 Systems and their classification Linear timeinvariant systems System Concept Mathematical transformation
More informationAPPM 2360: Midterm exam 1 February 15, 2017
APPM 36: Midterm exam 1 February 15, 17 On the front of your Bluebook write: (1) your name, () your instructor s name, (3) your recitation section number and () a grading table. Text books, class notes,
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 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 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 informationEquilibrium points: continuoustime systems
Capitolo 0 INTRODUCTION 81 Equilibrium points: continuoustime systems Let us consider the following continuoustime linear system ẋ(t) Ax(t)+Bu(t) y(t) Cx(t)+Du(t) The equilibrium points x 0 of the system
More information5. Observerbased Controller Design
EE635  Control System Theory 5. Observerbased Controller Design Jitkomut Songsiri state feedback poleplacement design regulation and tracking state observer feedback observer design LQR and LQG 51
More informationApplications of CAS to analyze the step response of a system with parameters
Applications of CAS to analyze the step response of a system with parameters Takuya Kitamoto kitamoto@yamaguchiu.ac.jp Faculty of Education Yamaguchi University 7538513 Japan Abstract Recently, Computer
More informationControl Systems I. Lecture 2: Modeling. Suggested Readings: Åström & Murray Ch. 23, Guzzella Ch Emilio Frazzoli
Control Systems I Lecture 2: Modeling Suggested Readings: Åström & Murray Ch. 23, Guzzella Ch. 23 Emilio Frazzoli Institute for Dynamic Systems and Control DMAVT ETH Zürich September 29, 2017 E. Frazzoli
More informationECE504: Lecture 8. D. Richard Brown III. Worcester Polytechnic Institute. 28Oct2008
ECE504: Lecture 8 D. Richard Brown III Worcester Polytechnic Institute 28Oct2008 Worcester Polytechnic Institute D. Richard Brown III 28Oct2008 1 / 30 Lecture 8 Major Topics ECE504: Lecture 8 We are
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu.6 Signal Processing: Continuous and Discrete Fall 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
More informationToday s goals So far Today 2.004
Today s goals So far Feedback as a means for specifying the dynamic response of a system Root Locus: from the openloop poles/zeros to the closedloop poles Moving the closedloop poles around Today Proportional
More informationDynamical systems: basic concepts
Dynamical systems: basic concepts Daniele Carnevale Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata Fondamenti di Automatica e Controlli Automatici A.A. 20142015
More information11.2 Basic Firstorder System Methods
112 Basic Firstorder System Methods 797 112 Basic Firstorder 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 informationA Model of Evolutionary Dynamics with Quasiperiodic Forcing
paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 25, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth
More informationEEE582 Homework Problems
EEE582 Homework Problems HW. Write a statespace 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 informationEN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015
EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions
More informationObservability. It was the property in Lyapunov stability which allowed us to resolve that
Observability We have seen observability twice already It was the property which permitted us to retrieve the initial state from the initial data {u(0),y(0),u(1),y(1),...,u(n 1),y(n 1)} It was the property
More informationNonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México
Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October
More information信號與系統 Signals and Systems
Spring 2010 信號與系統 Signals and Systems Chapter SS2 Linear TimeInvariant Systems FengLi Lian NTUEE Feb10 Jun10 Figures and images used in these lecture notes are adopted from Signals & Systems by Alan
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 informationModule 1: Signals & System
Module 1: Signals & System Lecture 6: Basic Signals in Detail Basic Signals in detail We now introduce formally some of the basic signals namely 1) The Unit Impulse function. 2) The Unit Step function
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationCLTI System Response (4A) Young Won Lim 4/11/15
CLTI System Response (4A) Copyright (c) 20112015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2
More informationENSC327 Communications Systems 2: Fourier Representations. Jie Liang School of Engineering Science Simon Fraser University
ENSC327 Communications Systems 2: Fourier Representations Jie Liang School of Engineering Science Simon Fraser University 1 Outline Chap 2.1 2.5: Signal Classifications Fourier Transform Dirac Delta Function
More information2.004 Dynamics and Control II Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 2.004 Dynamics and Control II Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 8 I * * Massachusetts
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationAC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI)
More informationLinear Systems. Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output
Linear Systems Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output Our interest is in dynamic systems Dynamic system means a system with memory of course including
More 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 information3.2 Complex Sinusoids and Frequency Response of LTI Systems
3. Introduction. A signal can be represented as a weighted superposition of complex sinusoids. x(t) or x[n]. LTI system: LTI System Output = A weighted superposition of the system response to each complex
More information16. Local theory of regular singular points and applications
16. Local theory of regular singular points and applications 265 16. Local theory of regular singular points and applications In this section we consider linear systems defined by the germs of meromorphic
More informationSourceFree RC Circuit
First Order Circuits SourceFree RC Circuit Initial charge on capacitor q = Cv(0) so that voltage at time 0 is v(0). What is v(t)? Prof Carruthers (ECE @ BU) EK307 Notes Summer 2018 150 / 264 First Order
More informationA conjecture on sustained oscillations for a closedloop heat equation
A conjecture on sustained oscillations for a closedloop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping
More informationLecture Notes for Math 251: ODE and PDE. Lecture 22: 6.5 Dirac Delta and Laplace Transforms
Lecture Notes for Math : ODE and PDE. Lecture : 6. Dirac Delta and Laplace Transforms Shawn D. Ryan Spring 0 Dirac Delta and the Laplace Transform Last Time: We studied differential equations with discontinuous
More information