EE451/551: Digital Control. Chapter 7: State Space Representations
|
|
- Cynthia Dalton
- 6 years ago
- Views:
Transcription
1 EE45/55: Digital Control Chapter 7: State Spae Representations
2 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 tin the interval t, t f to uniquely determine e e the behaviour of the system in the interval ev t, t f ; where n is the order of the system The set of state variables for the state vetor is written: x [ ] T x 2 xt = x, x2,, x n = x n } )
3 State Variables The state-spae spae is the n-dimensional vetor spae where the x ( t) from x( t) represent the oordinate axes i For a seond-order system, the state spae is two- dimensional and known as the state - plane; for the speial ase where the state-variables are proportional to the derivatives of fthe output, t the state t plane is alled the phase plane and the state variables are alled phase variables Curves in state spae are known as the state trajetories, and a plot of the trajetories in the plane is the state portrait (or phase portrait for the phase plane)
4 Example 7.: Motion of a Pt. Mass Consider the motion of a point mass m driven by a fore u: u my = u y = m where y is the displaement of the point mass; defining the system states as the phase variables: x ( t) y( t) x ( t) = y ( t) = x ( t) 2 x( t) yt x ( t) = yt = 2 2 u m
5 Example 7.: Motion of a Pt. Mass If the state vetor is given by: () T x t xt () = [ x ] (), t x2() t = x 2() t then the system dynamis in state-variable form is: x2() t x x x = = = u x 2 + x u Ax+ Bu 2 m m x [ ] y = Cx+ Du x 2
6 Questions?
7 The General State Spae Forms The general form of a linear state-spaespae equation is: xt = Axt + But yt () = C xt () + D ut () where xt is a n real state vetor, ut is a m real input vetor, y () t is a p real output vetor, and the matries are defined as: A is a n n real state matrix, B is a n m real input matrix, C is a p n real output matrix, D is a p m real feedforward matrix
8 The General State Spae Forms The general form of a nonlinear state-spae spae equation is: xt = f xt, ut yt = g xt, ut where f () i and g() i are real n and p vetors of funtions that guarantee the existene and uniqueness of a solution to the system dyanmis The nonlinear system an be approximated at an equilibrium point x, u using a linear model, as speified previously, using the onept of a system's Jaobians, as defined on the next slide
9 System Jaobians System Jaobians are defined for equilibrium point x, u as: A C f f f f x x u u ( x ), u n x, u x, u m ( x, u ) =, B = f n f n fn f n x x ( x, u) n u u ( x, u) ( x, u) m ( x, u) g g g g x x u u ( x, u ) n ( x, u) ( x, u) m ( x, u ) =, D = g p g p gp g p x x ( x, u) n u u ( x, u) ( x, u) m ( x, u)
10 Linearization Example Given the nonlinear mass-spring-damper system: 2 my + by + k sin( y) = u b 2 k u y = y sin( y) + m m m The system state an be defined as: x2 y x y x = = x = = b 2 k u y x 2 y x sin( x ) 2 + m m m f = f 2 Equilibrium points our where x =, whih implies 2 yy, = x =, x = sin u k The system Jaobians an now alulated as shown:
11 Linearization Example If u = x = sin = x = f f x x ( x, u 2 ) ( x, u ) f f os( ) 2 2 m m x x ( x, u) 2 ( x, u) A = = k = k f u ( x, u) B = = f 2 m u ( x, u )
12 If Linearization Example y = x = g () i is onsidered the system output, then: g g C [ ] = = x x ( x, u ) 2 ( x, u ) g D = = [ ] u ( x, u ) Therefore, the linearized model at ( x, u ) =, is: x = x k x + u 2 m m y x = + x [ ] [ ] 2 u
13 Questions?
14 Soln. of Linear State Equations If the state equaton with ICs are given as: xt () = A xt () + B ut (), xt Then the Laplae transform yields: Solving for sx () s x( t ) = A X () s + B U () s X( s) yields: n [ ] [ U() s] = Φ () s [ x( t ) + BU() s] X() s = si A x( t ) + B A n state - transition resolvent matrix Φ s = [ si ] where the has the inverse Φ ( t t ) = L si A n = e = k = ( ) ( τ ) ( ) { [ ] } A A t t t known as the matrix exponential resulting in the state equation: t A t t A t e e ( τ ) τ Φ ( t ) + Φ ( τ ) B u ( τ ) d τ = + = Φ + Φ xt xt Bu d t xt t zero input response t zero state response t t k k! t k
15 Example of State Equation Soln. T Given the system with unit step input and x() = [,] : x x = + u x2 x y = [ ] + [ ] u x 2 The state transition matrix Φ ( t) = si A an be omputed as shown: L L {[ ] } n t s s s( s + ) e L t s + e Φ () t = = = s +
16 Transfer Funtion Matrix As shown, the state equation with zero ICs an be written in the s-domain: = = [ ] ( s) X() s si A BU() s Φ BU() s n Sub. into the s-domain expression of the output equation yields: Y() s= C Φ () s BU() s + DU() s Y () s = = U() s = Φ () t H () s C Φ B + D h () () t C Φ B + D δ () t = () s Transfer Funtion (TF) At Ce B + Dδ () t Impulse Response (IR)
17 Transfer Funtion Matrix Given the prior example, we an ompue the system TF and IR as: H () s h Y () s s s ( s+ ) = = [ ] [ ] U() s + = s s+ s + t e () t = [ ] [ ] δ () t e t + = e t Do these results agree with prior state soln., see Matlab verifiation? >> A=[,;, ];B=[;];C=[,];D=[]; >> H_s=tf(ss(A,B,C,D)) Transfer funtion: s^2 + s >> zpk(h_s) zero/pole/gain: s (s+)
18 Questions?
19 Disrete State Equations We an develop disrete state eqns. of sample data system diretly form the assoiated ontinuous state eqns. assuming that a sample and zero-order hold is plaed prior to the plant G ( s) p = Y () s U() s Given this, we an assume that t = kt and t = kt + T in the solution of the ontinuous state eqn. developed previously, to obtain: kt + T ( + ) =Φ + ( τ) τ Φ + kt xkt T T xkt kt T Bd ukt
20 Disrete State Equations Based on this, the resulting disrete e state espae eeq eqn..beomes: es: x k+ = A x( k) + B u( k) To obtain this let y k = C x ( k ) + D u ( k ) d d d kt τ = σ in the prior eqn. T where A =Φ ( T ), B = Φ ( T σ ) d σ B, C = C and D = D d d d d As noted previously, Φ () t an be found by series expansion: A T A T Φ = ! 3! T In AT, thus Bd T T B I T A A B 2! 3! d = n d is given as
21 Example Disrete State Eqn. Given the sampled data system with T =. and G ( s) =, find the disrete state spae representation p s s + t As shown previously, e Φ ; therefore, t = t e T T t T e T T + e Φ ( τ ) d τ = t = T e e.952 Thus, A d =Φ ( T ) = and T =..95 T Bd = Φ( τ) dτ B = = T =.
22 Soln. to the Disrete State Eqn. Given a disrete state eqn. of the form: xk+ = Axk + Buk the soln. an be found by reursion as: d k k k i d () d d zero input i= response zero state response xk A x A Bui = + k wh ere Φ k = A is the disrete state transiton matrix wh d d Note, the reursion is easily performed on a omputer d
23 Disrete Transfer Funtion Matrix The disrete state equation with zero ICs an be written in the z-domain: = = [ ] ( z) X( z) zi A B U( z) Φ B U( z) n d d d Sub. into the s-domain expression of the output equation yields: Y( z) = C Φ ( z) B U( z) + D U( z) d d d d Y ( z ) H( z) = = C Φ ( z) B + D U( z) hk = C Φ ( kb ) + D δ ( k ) d d d d d d d d Ak d = Ce B + Dδ ( k) d d d d Transfer Funtion (TF) Impulse Response (IR)
24 Disrete Transfer Funtion Matrix z.952 Φ = = d n d z.95 Note: ( z) [ zi A ] = ( z ) ( z )( z.95) ( z.95) k ( k) δ( k).52δ( k) +.2 ( k) Φ ( k) = d k
25 Disrete Transfer Funtion Matrix Given the prior example, we an ompute the system TF and IR as:.952 Y( z) z z z H ( z ) = = [ ] U( z) z.95 [ ] [ ] ( z ) z z z.952 (.95) ( z )( z.95 ) = = hk = k ( k ) [ ] δ k δ k + k k [] δ ( k ) k =.2 ( k) δ ( k)
26 Disrete Transfer Funtion Matrix Do these results agree with the soln. based on the methods of Chapter 3? See Matlab verifiation... >>Ad=[,.952;,.95];Bd=[.484;.952]; Cd=[,];Dd=[]; >> T=.;GZAS_z=tf(ss(Ad,Bd,Cd,Dd,T)) Transfer funtion:.484 z z^2.95 z +.95 >> zpk(gzas) Zero/pole/gain: (z+.9672) (z ) (z.948) >> Gp_s=tf(,onv([ ],[ ])) Transfer funtion: s^2 + s >> zpk(gp_s) Zero/pole/gain: s(s+) >> GZAS2_z=2d(Gp_s,T,'zoh') Transfer funtion:.4837 z z^2.95 z >> zpk(gzas2_z) Zero/pole/gain: (z+.9672) (z ) (z.948)
27 Disrete Transfer Funtion Matrix Do these results agree with the soln. based on the methods of Chapter 3? See Matlab verifiation... >> z=tf('z',t) Transfer funtion: z Sampling time:. >> [num_gzas_zoz,den_gzas_zoz]= tfdata(minreal(gzas_z/z,e 4)); >>[R,P,K] = residue(num_gzas_zoz{},den_gzas_zoz{}) R = P =..95 K = [] >> A=[,;, ];B=[;];C=[,];D=[]; >> Gp_s=tf(ss(A,B,C,D)) Transfer funtion: s^2 + s >> 2d(ss(A,B,C,D),T,'zoh') a = x x2 x.956 x2.948 b = u x.4837 x = x x2 y d = u y
28 Questions?
29 Similarity Transforms State representations are not unique, e.g., assume that for a disrete time system that state eqns. are given by: x k+ = A x( k) + B u( k) xd xd y k C x k D u k = + xd and there a new state wk exists suh that xk = Pwk where P is an invertible n n matrix P exists Sub. w into the above state representation yields: wk + = P A Pwk + P B uk xd y k = C Pw( k) + D u( k) xd A = P A P, B = P B, C = C P, D = D wd xd wd xd wd xd wd xd xd xd xd
30 Similarity Transforms It is shown in Chapter 8 that system stability is determined by the roots of the harateristi eqn. of the A matrix, defined by zi A =, and alled system eigenvalues n whih are analogous to the poles d d of H( z) = C Φ ( z) B + D It an be shown that zi A = zi A, and also n xd n wd [ ] [ ] d d d d H ( z ) = C zi A B + D = C zi A B + D xd n xd xd xd wd n wd wd wd whih implies that system stability and TF remain unhanged under the similarity transform xk = Pwk
31 In lass Assignment Show that the system:.8 xk ( + ) = xk + uk.9 [ ] yk = xk has equivalent eigenvalues and TF under the similarity transform xk = Pwk where P =
32 Matlab Verifiation >> Axd=[.8 ;.9];Bxd=[;];Cxd=[ ];Dxd=[];P=[ ; ]; >> Awd=inv(P)*Axd*P Awd = >> Bwd=inv(P)*Bxd Bxd Bwd =.5.5 >> Cwd=Cxd*P Cwd = >> Dwd=Dxd Dxd Dwd = >> eig(axd) ans =.8.9 >> eig(awd) ans =.9.8 >> Hx_z=tf(ss(Axd,Bxd,Cxd,Dxd, )) Transfer funtion: z^2.7 z +.72 Sampling time: unspeified >> zpk(hx_z) Zero/pole/gain: (z.9) (z.8) >> Hw_z=tf(ss(Awd,Bwd,Cwd,Dwd, )) Transfer funtion: z^2.7 z +.72 >> zpk(hw_z) Zero/pole/gain: (z.9) (z.8)
33 Questions?
EE451/551: Digital Control. Chapter 8: Properties of State Space Models
EE451/551: Digital Control Chapter 8: Properties of State Space Models Equilibrium State Definition 8.1: An equilibrium point or state is an initial state from which the system nevers departs unless perturbed
More informationEE451/551: Digital Control. Final Exam Review Fall 2013
EE45/55: Digital Control Final Exam Review Fall 03 Exam Overview The Final Exam will consist of four/five questions for EE45/55 students based on Chapters 7 and a bonus based on Chapters 8 9 (students
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 informationHomogeneous and particular LTI solutions
Homogeneous and particular LTI solutions Daniele Carnevale Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata Fondamenti di Automatica e Controlli Automatici A.A. 2014-2015
More informationDifferential Equations 8/24/2010
Differential Equations A Differential i Equation (DE) is an equation ontaining one or more derivatives of an unknown dependant d variable with respet to (wrt) one or more independent variables. Solution
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 informationControl Theory association of mathematics and engineering
Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology
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 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 informationLTI system response. Daniele Carnevale. Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata
LTI system response Daniele Carnevale Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata Fondamenti di Automatica e Controlli Automatici A.A. 2014-2015 1 / 15 Laplace
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 informationContinuous Dynamics Solving LTI state-space equations גרא וייס המחלקה למדעי המחשב אוניברסיטת בן-גוריון
Continuous Dynamics Solving LTI state-space equations גרא וייס המחלקה למדעי המחשב אוניברסיטת בן-גוריון 2 State Space Models For a causal system with m inputs u t R m and p outputs y t R p, an nth-order
More informationLinear System Theory
Linear System Theory Wonhee Kim Chapter 6: Controllability & Observability Chapter 7: Minimal Realizations May 2, 217 1 / 31 Recap State space equation Linear Algebra Solutions of LTI and LTV system Stability
More informationAdvanced Control Theory
State Space Solution and Realization chibum@seoultech.ac.kr Outline State space solution 2 Solution of state-space equations x t = Ax t + Bu t First, recall results for scalar equation: x t = a x t + b
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 4
Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas
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-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 informationMath 2374: Multivariable Calculus and Vector Analysis
Math 2374: Multivariable Calulus and Vetor Analysis Part 26 Fall 2012 The integrals of multivariable alulus line integral of salar-valued funtion line integral of vetor fields surfae integral of salar-valued
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 8: Solutions of State-space Models Readings: DDV, Chapters 10, 11, 12 (skip the parts on transform methods) Emilio Frazzoli Aeronautics and Astronautics Massachusetts
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 informationEE 341 Homework Chapter 2
EE 341 Homework Chapter 2 2.1 The electrical circuit shown in Fig. P2.1 consists of two resistors R1 and R2 and a capacitor C. Determine the differential equation relating the input voltage v(t) to the
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 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 informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
More informationGeometry of Transformations of Random Variables
Geometry of Transformations of Random Variables Univariate distributions We are interested in the problem of finding the distribution of Y = h(x) when the transformation h is one-to-one so that there is
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 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 informationEL 625 Lecture 11. Frequency-domain analysis of discrete-time systems
EL 625 Leture 11 1 EL 625 Leture 11 Frequeny-domain analysis of disrete-time systems Theorem: If S is a fixed linear disrete-time system, the zero-state response, S{z k } = H(z)z k (1) Proof: S{z k } =
More informationPerspective. ECE 3640 Lecture 11 State-Space Analysis. To learn about state-space analysis for continuous and discrete-time. Objective: systems
ECE 3640 Lecture State-Space Analysis Objective: systems To learn about state-space analysis for continuous and discrete-time Perspective Transfer functions provide only an input/output perspective of
More 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 informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationModule 08 Observability and State Estimator Design of Dynamical LTI Systems
Module 08 Observability and State Estimator 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 November
More informationControl engineering sample exam paper - Model answers
Question Control engineering sample exam paper - Model answers a) By a direct computation we obtain x() =, x(2) =, x(3) =, x(4) = = x(). This trajectory is sketched in Figure (left). Note that A 2 = I
More informationDigital Control & Digital Filters. Lectures 13 & 14
Digital Controls & Digital Filters Lectures 13 & 14, Professor Department of Electrical and Computer Engineering Colorado State University Spring 2017 Systems with Actual Time Delays-Application 2 Case
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 informationIntegration of the Finite Toda Lattice with Complex-Valued Initial Data
Integration of the Finite Toda Lattie with Complex-Valued Initial Data Aydin Huseynov* and Gusein Sh Guseinov** *Institute of Mathematis and Mehanis, Azerbaijan National Aademy of Sienes, AZ4 Baku, Azerbaijan
More informationSystems and Control Theory Lecture Notes. Laura Giarré
Systems and Control Theory Lecture Notes Laura Giarré L. Giarré 2017-2018 Lesson 7: Response of LTI systems in the transform domain Laplace Transform Transform-domain response (CT) Transfer function Zeta
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 informationSTATE-SPACE DYNAMIC SYSTEMS (DISCRETE-TIME)
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.
More informationDefinition of the Laplace transform. 0 x(t)e st dt
Definition of the Laplace transform Bilateral Laplace Transform: X(s) = x(t)e st dt Unilateral (or one-sided) Laplace Transform: X(s) = 0 x(t)e st dt ECE352 1 Definition of the Laplace transform (cont.)
More informationBäcklund Transformations: Some Old and New Perspectives
Bäklund Transformations: Some Old and New Perspetives C. J. Papahristou *, A. N. Magoulas ** * Department of Physial Sienes, Helleni Naval Aademy, Piraeus 18539, Greee E-mail: papahristou@snd.edu.gr **
More informationECE 3793 Matlab Project 3 Solution
ECE 3793 Matlab Project 3 Solution Spring 27 Dr. Havlicek. (a) In text problem 9.22(d), we are given X(s) = s + 2 s 2 + 7s + 2 4 < Re {s} < 3. The following Matlab statements determine the partial fraction
More informationRemark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.
Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need
More informationLecture 19 Observability and state estimation
EE263 Autumn 2007-08 Stephen Boyd Lecture 19 Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time
More informationCDS 101/110: Lecture 2.1 System Modeling. Model-Based Analysis of Feedback Systems
CDS 11/11: Leture 21 System Modeling Rihard M Murray 5 Otober 215 Goals:! Define a model and its use in answering questions about a system! Introdue the onepts of state, dynamis, inputs and outputs! Review
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 informationEE/ME/AE324: Dynamical Systems. Chapter 7: Transform Solutions of Linear Models
EE/ME/AE324: Dynamical Systems Chapter 7: Transform Solutions of Linear Models The Laplace Transform Converts systems or signals from the real time domain, e.g., functions of the real variable t, to the
More informationAn Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems
An Integrated Arhiteture of Adaptive Neural Network Control for Dynami Systems Robert L. Tokar 2 Brian D.MVey2 'Center for Nonlinear Studies, 2Applied Theoretial Physis Division Los Alamos National Laboratory,
More informationWhere as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationEXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION
Journal of Mathematial Sienes: Advanes and Appliations Volume 3, 05, Pages -3 EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION JIAN YANG, XIAOJUAN LU and SHENGQIANG TANG
More informationDiscrete Bessel functions and partial difference equations
Disrete Bessel funtions and partial differene equations Antonín Slavík Charles University, Faulty of Mathematis and Physis, Sokolovská 83, 186 75 Praha 8, Czeh Republi E-mail: slavik@karlin.mff.uni.z Abstrat
More informationBrazilian Journal of Physics, vol. 29, no. 3, September, Classical and Quantum Mechanics of a Charged Particle
Brazilian Journal of Physis, vol. 9, no. 3, September, 1999 51 Classial and Quantum Mehanis of a Charged Partile in Osillating Eletri and Magneti Fields V.L.B. de Jesus, A.P. Guimar~aes, and I.S. Oliveira
More informationEnergy Gaps in a Spacetime Crystal
Energy Gaps in a Spaetime Crystal L.P. Horwitz a,b, and E.Z. Engelberg a Shool of Physis, Tel Aviv University, Ramat Aviv 69978, Israel b Department of Physis, Ariel University Center of Samaria, Ariel
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 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 informationLinear Systems Theory
ME 3253 Linear Systems Theory Review Class Overview and Introduction 1. How to build dynamic system model for physical system? 2. How to analyze the dynamic system? -- Time domain -- Frequency domain (Laplace
More informationSolving Dynamic Equations: The State Transition Matrix II
Reading the Text Solving Dynamic Equations: The State Transition Matrix II EGR 326 February 27, 2017 Just a reminder to read the text Read through longer passages, to see what is connected to class topics.
More informationME Fall 2001, Fall 2002, Spring I/O Stability. Preliminaries: Vector and function norms
I/O Stability Preliminaries: Vector and function norms 1. Sup norms are used for vectors for simplicity: x = max i x i. Other norms are also okay 2. Induced matrix norms: let A R n n, (i stands for induced)
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 informationHILLE-KNESER TYPE CRITERIA FOR SECOND-ORDER DYNAMIC EQUATIONS ON TIME SCALES
HILLE-KNESER TYPE CRITERIA FOR SECOND-ORDER DYNAMIC EQUATIONS ON TIME SCALES L ERBE, A PETERSON AND S H SAKER Abstrat In this paper, we onsider the pair of seond-order dynami equations rt)x ) ) + pt)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 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 information). In accordance with the Lorentz transformations for the space-time coordinates of the same event, the space coordinates become
Relativity and quantum mehanis: Jorgensen 1 revisited 1. Introdution Bernhard Rothenstein, Politehnia University of Timisoara, Physis Department, Timisoara, Romania. brothenstein@gmail.om Abstrat. We first
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 state-space
More informationMA2331 Tutorial Sheet 5, Solutions December 2014 (Due 12 December 2014 in class) F = xyi+ 1 2 x2 j+k = φ (1)
MA2331 Tutorial Sheet 5, Solutions. 1 4 Deember 214 (Due 12 Deember 214 in lass) Questions 1. ompute the line integrals: (a) (dx xy + 1 2 dy x2 + dz) where is the line segment joining the origin and the
More informationSURFACE WAVES OF NON-RAYLEIGH TYPE
SURFACE WAVES OF NON-RAYLEIGH TYPE by SERGEY V. KUZNETSOV Institute for Problems in Mehanis Prosp. Vernadskogo, 0, Mosow, 75 Russia e-mail: sv@kuznetsov.msk.ru Abstrat. Existene of surfae waves of non-rayleigh
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 informationChapter 8 Hypothesis Testing
Leture 5 for BST 63: Statistial Theory II Kui Zhang, Spring Chapter 8 Hypothesis Testing Setion 8 Introdution Definition 8 A hypothesis is a statement about a population parameter Definition 8 The two
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 informationEE480.3 Digital Control Systems. Part 7. Controller Design I. - Pole Assignment Method - State Estimation
EE480.3 Digital Control Systems Part 7. Controller Design I. - Pole Assignment Method - State Estimation Kunio Takaya Electrical and Computer Engineering University of Saskatchewan February 10, 2010 **
More informationTime Domain Method of Moments
Time Domain Method of Moments Massahusetts Institute of Tehnology 6.635 leture notes 1 Introdution The Method of Moments (MoM) introdued in the previous leture is widely used for solving integral equations
More informationControl System Design
ELEC4410 Control System Design Lecture 19: Feedback from Estimated States and Discrete-Time Control Design Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science
More information9 Controller Discretization
9 Controller Discretization In most applications, a control system is implemented in a digital fashion on a computer. This implies that the measurements that are supplied to the control system must be
More informationOne-edge Diffraction. Chapter 1. Direction of incidence. Reflected ray 4. Screen y
Chapter One-edge Diffration The time domain impulse response for edge diffration is derived from Eqn.. aording to Fig... ϕ < ϕ < π is the angle between y and r. ϕ < ϕ < π is the angle between y and the
More informationME8281-Advanced Control Systems Design
ME8281 - Advanced Control Systems Design Spring 2016 Perry Y. Li Department of Mechanical Engineering University of Minnesota Spring 2016 Lecture 4 - Outline 1 Homework 1 to be posted by tonight 2 Transition
More informationChapter 2: Solution of First order ODE
0 Chapter : Solution of irst order ODE Se. Separable Equations The differential equation of the form that is is alled separable if f = h g; In order to solve it perform the following steps: Rewrite the
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 informationResearch Article Approximation of Analytic Functions by Solutions of Cauchy-Euler Equation
Funtion Spaes Volume 2016, Artile ID 7874061, 5 pages http://d.doi.org/10.1155/2016/7874061 Researh Artile Approimation of Analyti Funtions by Solutions of Cauhy-Euler Equation Soon-Mo Jung Mathematis
More informationEE451/551: Digital Control. Relationship Between s and z Planes. The Relationship Between s and z Planes 11/10/2011
/0/0 EE45/55: Digital Control Chater 6: Digital Control System Design he Relationshi Between s and Planes As noted reviously: s j e e e e r s j where r e and If an analog system has oles at: s n jn a jd
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 informationSystems and Control Theory Lecture Notes. Laura Giarré
Systems and Control Theory Lecture Notes Laura Giarré L. Giarré 2018-2019 Lesson 14: Rechability Reachability (DT) Reachability theorem (DT) Reachability properties (DT) Reachability gramian (DT) Reachability
More informationLecture Note #6 (Chap.10)
System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryook Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 7: State-space Models Readings: DDV, Chapters 7,8 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology February 25, 2011 E. Frazzoli
More informationControl Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) vs Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
More informationMath 151 Introduction to Eigenvectors
Math 151 Introdution to Eigenvetors The motivating example we used to desrie matrixes was landsape hange and vegetation suession. We hose the simple example of Bare Soil (B), eing replaed y Grasses (G)
More informationClassical Field Theory
Preprint typeset in JHEP style - HYPER VERSION Classial Field Theory Gleb Arutyunov a a Institute for Theoretial Physis and Spinoza Institute, Utreht University, 3508 TD Utreht, The Netherlands Abstrat:
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 informationEG4321/EG7040. Nonlinear Control. Dr. Matt Turner
EG4321/EG7040 Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear [System Analysis] and Control Dr. Matt
More informationEECS 120 Signals & Systems University of California, Berkeley: Fall 2005 Gastpar November 16, Solutions to Exam 2
EECS 0 Signals & Systems University of California, Berkeley: Fall 005 Gastpar November 6, 005 Solutions to Exam Last name First name SID You have hour and 45 minutes to omplete this exam. he exam is losed-book
More informationStabilisation of Linear Time-Invariant Systems subject to Output Saturation
Stabilisation of Linear Time-Invariant Systems subject to Output Saturation by Gijs Hilhorst, BSc A Thesis in Partial Fulfillment of the Requirement for the Degree of Master of Science at the Department
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 information1 sin 2 r = 1 n 2 sin 2 i
Physis 505 Fall 005 Homework Assignment #11 Solutions Textbook problems: Ch. 7: 7.3, 7.5, 7.8, 7.16 7.3 Two plane semi-infinite slabs of the same uniform, isotropi, nonpermeable, lossless dieletri with
More informationRemarks Around Lorentz Transformation
Remarks Around Lorentz Transformation Arm Boris Nima arm.boris@gmail.om Abstrat After diagonalizing the Lorentz Matrix, we find the frame where the Dira equation is one derivation and we alulate the speed
More informationControl Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) and Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
More informationGeneration of EM waves
Generation of EM waves Susan Lea Spring 015 1 The Green s funtion In Lorentz gauge, we obtained the wave equation: A 4π J 1 The orresponding Green s funtion for the problem satisfies the simpler differential
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 informationVelocity Addition in Space/Time David Barwacz 4/23/
Veloity Addition in Spae/Time 003 David arwaz 4/3/003 daveb@triton.net http://members.triton.net/daveb Abstrat Using the spae/time geometry developed in the previous paper ( Non-orthogonal Spae- Time geometry,
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 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017
More information