Discrete-Time State-Space Equations. M. Sami Fadali Professor of Electrical Engineering UNR
|
|
- Georgina Flowers
- 5 years ago
- Views:
Transcription
1 Discrete-Time State-Space Equations M. Sami Fadali Professor of Electrical Engineering UNR 1
2 Outline Discrete-time (DT) state equation from solution of continuous-time state equation. Expressions in terms of constituent matrices. Solution of DT state equation. Example. 2
3 Solution of State Equation x t f = e A t f t 0 x t 0 t f + න e A t f τ Bu τ dτ t 0 Obtain difference equation from the solution of the analog state, over a sampling period T. Initial time t 0 = kt, and final time t f = k + 1 T k+1 T x k + 1 = e AT x k + න kt e A k+1 T τ Bu τ dτ x k = x kt = state vector at time kt 3
4 Digital Control Piecewise constant inputs over a sampling period. u t = u k, Move input outside the integral x k + 1 = e AT x k + න kt k+1 T kt t k + 1 T e A k+1 T τ Bdτ u k 4
5 Simplify Integral Change of variable of integration λ = k + 1 T τ, λ = ቊ 0, τ = k + 1 T T, τ = kt dλ = dτ න kt k+1 T e A k+1 T τ Bdτ = න 0 T T = න 0 e Aλ B dλ e Aλ Bdλ 5
6 Discrete State-Space Equations x k + 1 = A d x k + B d u k, k = 0,1, T A d = e AT, B d = න e Aτ Bdτ 0 y k = Cx k + Du k Obtain A d using state-transition matrix of A Integrate to obtain B d 6
7 State & Input Matrices A d = n n discrete state matrix B d = n m discrete input matrix (same orders as their continuous counterparts). Discrete state matrix = state transition matrix of the analog system evaluated at the sampling period T. Properties of the matrix exponential: integral of the matrix exponential for invertible matrix A e At dt = A 1 e At = e At A 1 B d = A 1 e AT I n B = e AT I n A 1 B 7
8 Constituent Matrices Use expansion of the matrix exponential in terms of the constituent matrices. Eigenvalues of discrete state matrix related to those of the analog system. n n A d = Z i e λ it = Z i e λ it i=1 t=t i=1 8
9 Input Matrix T T B d = න e Aτ Bdτ = න 0 0 n B d = i=1 n Z i e λ iτ i=1 T Z i B න e λiτ dτ 0 Bdτ Scalar integrands: easily evaluate integral. B d = n i=1 Z i B 1 eλ it λ i, λ i 0 n Z 1 BT + i=2 Z i B 1 eλ it λ i, λ 1 = 0 Assume distinct eigenvalues (only one zero eigenvalue) 9
10 Discrete-time State-space Representation Discrete state & output equation. Discrete-time state equation: approximately valid for a general input vector u t provided that the sampling period T is sufficiently short. x k + 1 = A d x k + B d u k Output equation evaluated at time kt y k = Cx k + Du k 10
11 ሶ Example 7.15 Obtain the DT state equations for the system of Example 7.7 x x 1 0 x2 ሶ = x u xሶ x 3 10 for a sampling period T=0.01 s. Solution: From Example 7.7, the state-transition matrix is e At = e e t e 10t 90 11
12 Discrete state matrix A d = e A e = e e Simplifies to A d =
13 Discrete-time Input matrix B d = Z 1 B Z 2 B 1 e Z 3 B 1 e = e e Simplifies to B d =
14 MATLAB Obtain A d, B d, C, D from A, B, C, D» A=[1,2;3,4]; B=[0;1];C=[1,1]; D=0;» p=ss(a, B, C, D); % State-space matrices» pd = c2d(p,0.01) % piecewise constant input (default) Alternatively, use the MATLAB commands» ad = expm(a * 0.05)» bd = a\ (ad-eye(3) )* b 14
15 a = MATLAB Output x1 x2 x x b = u1 x x c = x1 x2 y1 1 1 d = u1 y1 0 Sampling time: 0.01 Discrete-time model. 15
16 Solution of DT State-Space Equation DT State Equation: state at time k in terms of the initial condition vector x 0 and the input sequence u k, k = 0, 1,, k 1. x k + 1 = A d x k + B d u k At k = 0, 1, we have x 1 = A d x 0 + B d u 0 x 2 = A d x 1 + B d u 1 = A d 2 x 0 + A d B d u 0 + B d u 1 16
17 Solution by Induction 2 1 x 2 = A 2 d x 0 + A 2 i 1 d B d u i i=0 k 1 x k = A k d x 0 + A k i 1 d B d u i i=0 k 1 k k x k = A 0 d x k 0 + A k i 1 d B d u i Solution i=k 0 =zero-input response+ zero-state response 17
18 State-transition Matrix State-transition matrix for the DT system A d k. State-transition matrix for time-varying DT system: not a matrix power dependent on both time k and initial time k 0 18
19 Output Solution Substitute in output equation y k = Cx k + Du k k 1 = C A k d x 0 + A k i 1 d B d u i i=0 + Du k 19
20 Z-Transform Solution of DT State Equation z-transform the discrete-time state equation zx z zx 0 = A d X z + B d U z X z = zi n A 1 d zx 0 + B d U z zi n A 1 d z = I n 1 z A d = I n + A d z 1 + A 2 d z A i d z i
21 Inverse z-transform Inverse z-transform zi n A d 1 z zi n A d 1 z = I n + A d z A d i z i + Z 1 zi n A d 1 z = I n, A d, A d 2,, A d i, = A d k k=0 Analogous to the scalar transform pair z z a d Z a d k k=0 21
22 Matrix Inversion Evaluate using the Leverrier algorithm. Partial fraction expansion then multiply by z. zi A d 1 z = z adj zi A d det zi A d = P 0z + P 1 z P n 1 z n a 0 + a 1 z + + a n 1 z n 1 + z n n = i=1 z z λ i Z i, n A k k d = Z i λ i i=1 22
23 DT State Matrix n A d = i=1 n Z i λ i = i=1 Z i A e λ i A T Parentheses: pertaining to the CT state matrix A. Equality for any sampling period T and any matrix A Z i = Z i A λ i = e λ i A T Same constituent matrices for DT state matrix & CT state matrix A DT eigenvalues are exponential functions of the CT eigenvalues times the sampling period. 23
24 Zero-state Response x ZS t = zi n A d 1 z z 1 B d U z Known inverse transform for. term. Multiplication by z 1 : delay by T. Convolution theorem: inverse of product is the convolution summation k 1 x zs k = A k i 1 d B d u i, A k d = Z j e λ j A kt i=0 j=1 k 1 n x zs k = Z j e λ j A k i 1 T B d u i i=0 j=1 n 24
25 Alternative Expression k 1 n x zs k = Z j e λ j A k i 1 T B d u i i=0 j=1 Change the order of summation n x zs k = j=1 k 1 Z j B d e λ j A k 1 T i=0 u i e λ j A it Useful when the i summation has a closed form. 25
26 Step Inputs n k 1 x zs k = Z j B d e λ j A k 1 T u i e λ j A it j=1 i=0 u k = 1 = 1 1 T k 1 a i = 1 ak 1 a, a 1 i=1 n x zs k = j=1 Z j B d 1e λ j A k 1 T 1 e λ j A kt 1 e λ j A T 26
27 ሶ Example 7.16 x = x u (a) Solve the state equation for a unit step input and the initial condition vector x 0 = 1 0 T (b) Use the solution to obtain the discrete-time state equations for a sampling period of 0.1s. (c) Solve the discrete-time state equations with the same initial conditions and input as in (a) and verify that the solution is the same as that of (a) evaluated at multiples of the sampling period T. 27
28 Solution (a)the resolvent matrix s Φ s = s s = 2 s + 3 s 2 + 3s s = s + 1 s + 2 Partial fraction expansions s + 1 s + 2 = s s + 1 s + 2 = s s s s
29 State-transition Matrix s Φ s = s + 1 s = s = s + 1 s + 2 φ t = e t s e 2t 29
30 Zero-input Response x ZI t = e At x 0 = e t + = 2 2 e t e 2t e 2t
31 Zero-state Response x ZS t = e t e 2t t = 1 1 e t 1 t e 2t 1 t x ZS t = = 1/ (1 e t ) (1 e 2t )/2 1 1 e t 1/2 1 e 2t 31
32 Total Response x t = x ZI t + x ZS t x t = 2 2 e t e 2t + 1/2 0 = 1/ e t 1/ e t + 1/2 1 At the sampling points: e 2t e 2t t = multiples of 0.1s 32
33 (b) Discrete-time state equations A d = φ 0.1 = e e 0.2 = B d = e e = 1/2 0 e e = x 1 = A d x 0 + B d u 0 CT system response to a step input of duration one sampling period=response of CT system due to a piecewise constant input over one period 33
34 Zero-input Response A d k = Φ 0.1k = e 0.1k + = e 0.1k + x ZI k = e 0.2k x ZI k = Φ 0.1k x e 0.2k e 0.1k e 0.2k Same as the zero-input response of the continuous-time system x ZI t = 2 2 e t e 2t at all sampling points k = 0, 1, 2, 34
35 z-transform of the zero-state Φ z = X ZS z = Φ z z 1 B d U z response z z e z z e 0.2 = 2 1 z z e = z z e z z e 0.1 z 1 z z e 0.2 z 1 z 1 z z 1 35
36 Partial Fractions z z e 0.1 z 1 = 1 z 1 e 0.1 z z z e 0.1 = z z z z z z e 0.2 z 1 = 1 z 1 e 0.2 z z z e 0.1 = z z z z
37 Expand zero-state response X ZS z 1 z = 2 z 1 0 x ZS k = 1/ z z e e 0.1k 1/2 1 z z e 0.2 e 0.2k Same as zero-state response for the CT system x ZS t = 1/ e t 1/2 e 2t 1 at time t = 0.1 k, k = 0,1,2, 37
38 Zero-state response k 1 n x zs k = Z j e λ j A k i 1 T B d u i i=0 j=1 n k 1 x zs k = Z j B d e λ j A k 1 T u i e λ j A it j=1 i=0 k 1 u k = 1 = 1 1 T, a i = 1 ak 1 a, a 1 n x zs k = j=1 i=1 Z j B d 1e λ j A k 1 T 1 eλ j A kt 1 e λ j A T 38
39 z-transfer Function x k + 1 = A d x k + B d u k y k = Cx k + Du k Y z = CX s + DU z For zero initial conditions zx z = A d X z + B d U z X z = zi n A d 1 B d U z Substitute Y z = C zi n A d 1 B d U z + DU z = C zi n A d 1 B d + D U z = G z U z 39
40 Impulse Response & Modes G z = C zi n A d 1 B d + D Inverse transform of the transfer function G z = C zi n A d 1 B d + D Z G k = ቊ CA d k 1 B d, k 1 D, k = 0 Substitute in terms of constituent matrices n G z = i=1 CZ i B d 1 z λ i + D Z G k = n i=1 CZ i B d λ i k, k 1 D, k = 0 40
41 Poles and Stability Poles= eigenvalues of discrete-time state matrix A d = exponential functions of λ i A (continuous-time state matrix A). For stable A, λ i A have negative real parts and have magnitude less than unity. Discretization yields a stable DT system for a stable CT system. 41
42 Minimal Realizations Product CZ j B can vanish & eliminate eigenvalues from the transfer function: if CZ j = 0, Z j B = 0, or both. If cancellation occurs, the system is said to have an CZ j = 0 : output-decoupling zero at λ j Z j B = 0: an input-decoupling zero at λ j CZ j = 0, & Z j B = 0: an input-output-decoupling zero at λ j Poles of the reduced transfer function are a subset of the eigenvalues of the state matrix A d. A state-space realization that leads to pole-zero cancellation is said to be reducible or nonminimal. If no cancellation occurs, the realization is said to be irreducible or minimal. 42
43 Decoupling Modes Output-decoupling zero at λ j : the forced system response does not include the mode λ j k. Input-decoupling zero: the mode is decoupled from or unaffected by the input. Input-output-decoupling zero: the mode is decoupled both from the input and the output. These properties are related to the concepts of controllability and observability discussed later in this chapter. 43
44 Example 7.17 Obtain the z-transfer function for the position control system of Example 7.16 (a) With x 1 as output. (b) With x 1 + x 2 as output. The resolvent matrix and input matrix were obtained in Example
45 Solution Φ z z 1 = z e z e 0.2 B d = 1/ e e = (a) Output x 1 : C = 1 0, D = G z = z e z e = z e 0.1 z e
46 (c) Output x 1 +x z 2 1 G z = z e z z e = + z e 0.1 z e Output-decoupling zero at e 0.1 since CZ 1 = [0] 2. System response to any input does not include the decoupling term. 46
47 Step Response Y z = z e 0.2 z z z = 1 e 0.2 z z z e 0.2 z = 0.5 z z z e 0.2 Z-transform inverse y k = e 0.2k 47
48 z-transfer Function: MATLAB Let T = 0.05s» P = ss(ad, Bd, C, D, 0.05)» g = tf(p) % Obtain z-domain transfer function» zpk(g) % Obtain transfer function poles and zeros The command reveals that the system has a zero at and poles at (0.9048, ) with a gain of With pole-zero cancellation, the transfer function is the same as that of Example 7.17(b).» minreal(g) % Cancel poles and zeros 48
Digital Control System Models. M. Sami Fadali Professor of Electrical Engineering University of Nevada
Digital Control System Models M. Sami Fadali Professor of Electrical Engineering University of Nevada 1 Outline Model of ADC. Model of DAC. Model of ADC, analog subsystem and DAC. Systems with transport
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 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 informationLinear Algebra Review (Course Notes for Math 308H - Spring 2016)
Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,
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 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 informationStochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno
Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.
More informationCayley-Hamilton Theorem
Cayley-Hamilton Theorem Massoud Malek In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n Let A be an n n matrix Although det (λ I n A
More informationConvolution and Linear Systems
CS 450: Introduction to Digital Signal and Image Processing Bryan Morse BYU Computer Science Introduction Analyzing Systems Goal: analyze a device that turns one signal into another. Notation: f (t) g(t)
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 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 information(V.C) Complex eigenvalues and other topics
V.C Complex eigenvalues and other topics matrix Let s take a closer look at the characteristic polynomial for a 2 2 namely A = a c f A λ = detλi A = λ aλ d bc b d = λ 2 a + dλ + ad bc = λ 2 traλ + det
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 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 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 informationMath 1060 Linear Algebra Homework Exercises 1 1. Find the complete solutions (if any!) to each of the following systems of simultaneous equations:
Homework Exercises 1 1 Find the complete solutions (if any!) to each of the following systems of simultaneous equations: (i) x 4y + 3z = 2 3x 11y + 13z = 3 2x 9y + 2z = 7 x 2y + 6z = 2 (ii) x 4y + 3z =
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 informationZ - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.
Z - Transform The z-transform is a very important tool in describing and analyzing digital systems. It offers the techniques for digital filter design and frequency analysis of digital signals. Definition
More informationECEEN 5448 Fall 2011 Homework #4 Solutions
ECEEN 5448 Fall 2 Homework #4 Solutions Professor David G. Meyer Novemeber 29, 2. The state-space realization is A = [ [ ; b = ; c = [ which describes, of course, a free mass (in normalized units) with
More informationLecture 2 and 3: Controllability of DT-LTI systems
1 Lecture 2 and 3: Controllability of DT-LTI systems Spring 2013 - EE 194, Advanced Control (Prof Khan) January 23 (Wed) and 28 (Mon), 2013 I LTI SYSTEMS Recall that continuous-time LTI systems can be
More informationWe use the overhead arrow to denote a column vector, i.e., a number with a direction. For example, in three-space, we write
1 MATH FACTS 11 Vectors 111 Definition We use the overhead arrow to denote a column vector, ie, a number with a direction For example, in three-space, we write The elements of a vector have a graphical
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 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 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 informationChap 2. Discrete-Time Signals and Systems
Digital Signal Processing Chap 2. Discrete-Time Signals and Systems Chang-Su Kim Discrete-Time Signals CT Signal DT Signal Representation 0 4 1 1 1 2 3 Functional representation 1, n 1,3 x[ n] 4, n 2 0,
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2016 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationMULTIVARIABLE ZEROS OF STATE-SPACE SYSTEMS
Copyright F.L. Lewis All rights reserved Updated: Monday, September 9, 8 MULIVARIABLE ZEROS OF SAE-SPACE SYSEMS If a system has more than one input or output, it is called multi-input/multi-output (MIMO)
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 information1.4 Unit Step & Unit Impulse Functions
1.4 Unit Step & Unit Impulse Functions 1.4.1 The Discrete-Time Unit Impulse and Unit-Step Sequences Unit Impulse Function: δ n = ቊ 0, 1, n 0 n = 0 Figure 1.28: Discrete-time Unit Impulse (sample) 1 [n]
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 informationEE451/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 information06/12/ rws/jMc- modif SuFY10 (MPF) - Textbook Section IX 1
IV. Continuous-Time Signals & LTI Systems [p. 3] Analog signal definition [p. 4] Periodic signal [p. 5] One-sided signal [p. 6] Finite length signal [p. 7] Impulse function [p. 9] Sampling property [p.11]
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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.4: Dynamic Systems Spring Homework Solutions Exercise 3. a) We are given the single input LTI system: [
More informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2018 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationUsing MATLAB with the Convolution Method
ECE 350 Linear Systems I MATLAB Tutorial #5 Using MATLAB with the Convolution Method A linear system with input, x(t), and output, y(t), can be described in terms of its impulse response, h(t). x(t) h(t)
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 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 informationEigenvalues and Eigenvectors
5 Eigenvalues and Eigenvectors 5.2 THE CHARACTERISTIC EQUATION DETERMINANATS nn Let A be an matrix, let U be any echelon form obtained from A by row replacements and row interchanges (without scaling),
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 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 informationHomework Set #8 Solutions
Exercises.2 (p. 19) Homework Set #8 Solutions Assignment: Do #6, 8, 12, 14, 2, 24, 26, 29, 0, 2, 4, 5, 6, 9, 40, 42 6. Reducing the matrix to echelon form: 1 5 2 1 R2 R2 R1 1 5 0 18 12 2 1 R R 2R1 1 5
More informationEIGENVALUES AND EIGENVECTORS 3
EIGENVALUES AND EIGENVECTORS 3 1. Motivation 1.1. Diagonal matrices. Perhaps the simplest type of linear transformations are those whose matrix is diagonal (in some basis). Consider for example the matrices
More informationLinear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0.
Matrices Operations Linear Algebra Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0 The rectangular array 1 2 1 4 3 4 2 6 1 3 2 1 in which the
More informationKalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and
Kalman Decomposition Controllable / uncontrollable decomposition Suppose that the controllability matrix C R n n of a system has rank n 1
More 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 informationChapter 3 Convolution Representation
Chapter 3 Convolution Representation DT Unit-Impulse 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 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 informationSIGNALS AND SYSTEMS. Unit IV. Analysis of DT signals
SIGNALS AND SYSTEMS Unit IV Analysis of DT signals Contents: 4.1 Discrete Time Fourier Transform 4.2 Discrete Fourier Transform 4.3 Z Transform 4.4 Properties of Z Transform 4.5 Relationship between Z
More informationLecture 20 Aspects of Control
CS 460/560 Introduction to Computational Robotics Fall 2017, Rutgers University Lecture 20 Aspects of Control Instructor: Jingjin Yu Outline Feedback (closed-loop) control Mathematical models of dynamical
More information2. CONVOLUTION. Convolution sum. Response of d.t. LTI systems at a certain input signal
2. CONVOLUTION Convolution sum. Response of d.t. LTI systems at a certain input signal Any signal multiplied by the unit impulse = the unit impulse weighted by the value of the signal in 0: xn [ ] δ [
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 informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
More informationDiscrete State Space Models
Discrete State Space Models In this lecture we will discuss the relation between transfer function and state space model for a discrete time system and various standard or canonical state variable models.
More informationDigital Signal Processing Lecture 5
Remote Sensing Laboratory Dept. of Information Engineering and Computer Science University of Trento Via Sommarive, 14, I-38123 Povo, Trento, Italy Digital Signal Processing Lecture 5 Begüm Demir E-mail:
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 informationStudy Guide for Linear Algebra Exam 2
Study Guide for Linear Algebra Exam 2 Term Vector Space Definition A Vector Space is a nonempty set V of objects, on which are defined two operations, called addition and multiplication by scalars (real
More informationA = 3 1. We conclude that the algebraic multiplicity of the eigenvalues are both one, that is,
65 Diagonalizable Matrices It is useful to introduce few more concepts, that are common in the literature Definition 65 The characteristic polynomial of an n n matrix A is the function p(λ) det(a λi) Example
More informationCHAPTER 5: LAPLACE TRANSFORMS
CHAPTER 5: LAPLACE TRANSFORMS SAMANTHA RAMIREZ PREVIEW QUESTIONS What are some commonly recurring functions in dynamic systems and their Laplace transforms? How can Laplace transforms be used to solve
More informationThe Laplace Transform
The Laplace Transform Generalizing the Fourier Transform The CTFT expresses a time-domain signal as a linear combination of complex sinusoids of the form e jωt. In the generalization of the CTFT to the
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 informationSteady-state DKF. M. Sami Fadali Professor EE UNR
Steady-state DKF M. Sami Fadali Professor EE UNR 1 Outline Stability of linear estimators. Lyapunov equation. Uniform exponential stability. Steady-state behavior of Lyapunov equation. Riccati equation.
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 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 informationRobust Control 2 Controllability, Observability & Transfer Functions
Robust Control 2 Controllability, Observability & Transfer Functions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /26/24 Outline Reachable Controllability Distinguishable
More informationL2 gains and system approximation quality 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION L2 gains and system approximation quality 1 This lecture discusses the utility
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 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 informationMATH 312 Section 7.1: Definition of a Laplace Transform
MATH 312 Section 7.1: Definition of a Laplace Transform Prof. Jonathan Duncan Walla Walla University Spring Quarter, 2008 Outline 1 The Laplace Transform 2 The Theory of Laplace Transforms 3 Conclusions
More informationDigital Control Engineering Analysis and Design
Digital Control Engineering Analysis and Design M. Sami Fadali Antonio Visioli AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY TOKYO Academic Press is
More informationMATRICES. a m,1 a m,n A =
MATRICES Matrices are rectangular arrays of real or complex numbers With them, we define arithmetic operations that are generalizations of those for real and complex numbers The general form a matrix of
More informationSpring, 2012 CIS 515. Fundamentals of Linear Algebra and Optimization Jean Gallier
Spring 0 CIS 55 Fundamentals of Linear Algebra and Optimization Jean Gallier Homework 5 & 6 + Project 3 & 4 Note: Problems B and B6 are for extra credit April 7 0; Due May 7 0 Problem B (0 pts) Let A be
More information(b) If a multiple of one row of A is added to another row to produce B then det(b) =det(a).
.(5pts) Let B = 5 5. Compute det(b). (a) (b) (c) 6 (d) (e) 6.(5pts) Determine which statement is not always true for n n matrices A and B. (a) If two rows of A are interchanged to produce B, then det(b)
More informationDIGITAL STABILIZATION OF LINEAR CONTINUOUS-TIME PERIODIC PROCESSES WITH PURE DELAY
DIGITAL STABILIZATION OF LINEAR CONTINUOUS-TIME PERIODIC PROCESSES WITH PURE DELAY B.P. Lampe E.N. Rosenwasser University of Rostock, Department of Computer Science and Electrical Engineering, D-18051
More informationBasics. A VECTOR is a quantity with a specified magnitude and direction. A MATRIX is a rectangular array of quantities
Some Linear Algebra Basics A VECTOR is a quantity with a specified magnitude and direction Vectors can exist in multidimensional space, with each element of the vector representing a quantity in a different
More information1.6: 16, 20, 24, 27, 28
.6: 6, 2, 24, 27, 28 6) If A is positive definite, then A is positive definite. The proof of the above statement can easily be shown for the following 2 2 matrix, a b A = b c If that matrix is positive
More informationMATH 2050 Assignment 8 Fall [10] 1. Find the determinant by reducing to triangular form for the following matrices.
MATH 2050 Assignment 8 Fall 2016 [10] 1. Find the determinant by reducing to triangular form for the following matrices. 0 1 2 (a) A = 2 1 4. ANS: We perform the Gaussian Elimination on A by the following
More informationUNIT 3 INTEGRATION 3.0 INTRODUCTION 3.1 OBJECTIVES. Structure
Calculus UNIT 3 INTEGRATION Structure 3.0 Introduction 3.1 Objectives 3.2 Basic Integration Rules 3.3 Integration by Substitution 3.4 Integration of Rational Functions 3.5 Integration by Parts 3.6 Answers
More informationECE 301 Division 1 Exam 1 Solutions, 10/6/2011, 8-9:45pm in ME 1061.
ECE 301 Division 1 Exam 1 Solutions, 10/6/011, 8-9:45pm in ME 1061. Your ID will be checked during the exam. Please bring a No. pencil to fill out the answer sheet. This is a closed-book exam. No calculators
More information6 The Fourier transform
6 The Fourier transform In this presentation we assume that the reader is already familiar with the Fourier transform. This means that we will not make a complete overview of its properties and applications.
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 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 informationChapter 6 Controllability and Obervability
Chapter 6 Controllability and Obervability Controllability: whether or not the state-space equation can be controlled from input. Observability: whether or not the initial state can be observed from output.
More informationNonlinear Control Lecture # 14 Input-Output Stability. Nonlinear Control
Nonlinear Control Lecture # 14 Input-Output Stability L Stability Input-Output Models: y = Hu u(t) is a piecewise continuous function of t and belongs to a linear space of signals The space of bounded
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 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 information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 12: I/O Stability Readings: DDV, Chapters 15, 16 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology March 14, 2011 E. Frazzoli
More informationAn Energy Balance Model for Arctic Sea Ice
An Energy Balance Model for Arctic Sea Ice Kaitlin Hill October 1, 217 Outline Motivation: Arctic sea ice decline Energy balance model Bifurcation analysis Results Arctic sea ice age Sept 1986 Sept 216
More information21.4. Engineering Applications of z-transforms. Introduction. Prerequisites. Learning Outcomes
Engineering Applications of z-transforms 21.4 Introduction In this Section we shall apply the basic theory of z-transforms to help us to obtain the response or output sequence for a discrete system. This
More informationDO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.
EE 537 Homewors Friedland Text Updated: Wednesday November 8 Some homewor assignments refer to Friedland s text For full credit show all wor. Some problems require hand calculations. In those cases do
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 informationA = 3 B = A 1 1 matrix is the same as a number or scalar, 3 = [3].
Appendix : A Very Brief Linear ALgebra Review Introduction Linear Algebra, also known as matrix theory, is an important element of all branches of mathematics Very often in this course we study the shapes
More informationFundamentals of Engineering Analysis (650163)
Philadelphia University Faculty of Engineering Communications and Electronics Engineering Fundamentals of Engineering Analysis (6563) Part Dr. Omar R Daoud Matrices: Introduction DEFINITION A matrix is
More informationA VERY BRIEF LINEAR ALGEBRA REVIEW for MAP 5485 Introduction to Mathematical Biophysics Fall 2010
A VERY BRIEF LINEAR ALGEBRA REVIEW for MAP 5485 Introduction to Mathematical Biophysics Fall 00 Introduction Linear Algebra, also known as matrix theory, is an important element of all branches of mathematics
More informationLinear Algebra Practice Problems
Math 7, Professor Ramras Linear Algebra Practice Problems () Consider the following system of linear equations in the variables x, y, and z, in which the constants a and b are real numbers. x y + z = a
More informationName: Final Exam MATH 3320
Name: Final Exam MATH 3320 Directions: Make sure to show all necessary work to receive full credit. If you need extra space please use the back of the sheet with appropriate labeling. (1) State the following
More informationEigenvalues and Eigenvectors. Review: Invertibility. Eigenvalues and Eigenvectors. The Finite Dimensional Case. January 18, 2018
January 18, 2018 Contents 1 2 3 4 Review 1 We looked at general determinant functions proved that they are all multiples of a special one, called det f (A) = f (I n ) det A. Review 1 We looked at general
More informationECE 546 Lecture 15 Circuit Synthesis
ECE 546 Lecture 15 Circuit Synthesis Spring 018 Jose E. Schutt-Aine Electrical & Computer Engineering University of Illinois jesa@illinois.edu ECE 546 Jose Schutt Aine 1 MOR via Vector Fitting Rational
More informationSAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015
FACULTY OF ENGINEERING AND SCIENCE SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015 Lecturer: Michael Ruderman Problem 1: Frequency-domain analysis and control design (15 pt) Given is a
More information