Topic # /31 Feedback Control Systems
|
|
- Stephany Crawford
- 6 years ago
- Views:
Transcription
1 Topic # /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
2 Fall / Time Response Can develop a lot of insight into the system response and how it is modeled by computing the time response x(t) Homogeneous part Forced solution Homogeneous Part Take Laplace transform so that But can show ẋ = Ax, x(0) known X(s) = (si A) 1 x(0) x(t) = L 1 (si A) 1 x(0) I A A 2 (si A) 1 = s s2 s3 1 so L 1 (si A) 1 = I + At + (At) ! At = e x(t) = e At x(0) e At is a special matrix that we will use many times in this course Transition matrix or Matrix Exponential Calculate in MATLAB using expm.m and not exp.m 1 Note that e (A+B)t = e At e Bt iff AB = BA 1 MATLAB is a trademark of the Mathworks Inc.
3 Fall / Example: ẋ = Ax, with A = (si A) 1 = = = At e = [ s 2 1 s + 3 ] 1 s s (s + 2)(s + 1) 2 s s s s s s s s + 2 2e t e 2t e t e 2t 2e t + 2e 2t e t + 2e 2t We will say more about e At when we have said more about A (eigenvalues and eigenvectors) Computation of e At = L 1 {(si A) 1 } straightforward for a 2-state system More complex for a larger system, see this paper
4 Fall / Forced Solution Consider a scalar case: SS: Forced Solution ẋ = ax + bu, x(0) given t x(t) = e at x(0) + e a(t τ) bu(τ)dτ 0 where did this come from? 1. ẋ ax = bu 2. e at [ẋ ax] = d (e at x(t)) = e at dt bu(t) t 3. d e aτ x(τ)dτ = e at x(t) x(0) = t e aτ bu(τ)dτ 0 dτ 0 Forced Solution Matrix case: ẋ = Ax + Bu where x is an n-vector and u is a m-vector Just follow the same steps as above to get t x(t) = e At x(0) + e A(t τ) Bu(τ )dτ and if y = Cx + Du, then y(t) = Ce At x(0) + 0 t 0 Ce A(t τ) Bu(τ)dτ + Du(t) Ce At x(0) is the initial response Ce A(t) B is the impulse response of the system.
5 Fall / Have seen the key role of e At in the solution for x(t) Determines the system time response But would like to get more insight! Consider what happens if the matrix A is diagonalizable, i.e. there exists a T such that λ 1 T 1 AT = Λ which is diagonal Λ =... Then e At = T e Λt T 1 where e λ 1 t Λt e =... e λ n t λ n Follows since e At = I + At + 1 (At) and that A = T ΛT 1 2!, so we can show that 1 e At = I + At + (At) ! 1 = I + T ΛT 1 t + (T ΛT 1 t) ! = T e Λt T 1 This is a simpler way to get the matrix exponential, but how find T and λ? Eigenvalues and Eigenvectors
6 Fall / Eigenvalues and Eigenvectors Recall that the eigenvalues of A are the same as the roots of the characteristic equation (page 6 7) λ is an eigenvalue of A if det(λi A) = 0 which is true iff there exists a nonzero v (eigenvector) for which (λi A)v = 0 Av = λv Repeat the process to find all of the eigenvectors. Assuming that the n eigenvectors are linearly independent Av i = λ i v i i = 1,..., n λ 1 A v 1 v n = v 1 v n... AT = T Λ T 1 AT = Λ λ n
7 Fall / Jordan Form One word of caution: Not all matrices are diagonalizable A = det(si A) = s 0 0 only one eigenvalue s = 0 (repeated twice). The eigenvectors solve 0 1 r 1 = r 2 r eigenvectors are of the form 1, r 1 = 0 would only be one. 0 Need Jordan Form to handle the case with repeated roots 2 Jordan form of matrix A R n n is block diagonal: A λ j A 0 λ j 1 0 A = 2 0 with A j = A k λ j Observation: any matrix can be transformed into Jordan form with the eigenvalues of A determining the blocks A j. The matrix exponential of a Jordan form matrix is then given by 1 t t 2 t /2! n 1 (n 1)! e A 1t 0 0 t n 2 0 e A 2t t (n 2)! At e = with e A j t.. = e.. A k t... t 0 0 e λ j t 2 see book by Strang, Matlab command Jordan, or here
8 Fall / EV Mechanics 1 1 Consider A = 8 5 s (si A) = 8 s 5 det(si A) = (s + 1)(s 5) + 8 = s 2 4s + 3 = 0 so the eigenvalues are s 1 = 1 and s 2 = 3 Eigenvectors (si A)v = 0 s v (s 1 I A)v 1 = 11 = 0 8 s 5 v 21 s=1 2 1 v 11 = 0 2v 8 v 11 v 21 = 0, v 21 = 2v v 11 is then arbitrary (= 0), so set v 11 = 1 1 v 1 = 2 (s 2 I A)v 2 = 4 1 v 12 = 0 4v 12 v 22 = 0, v 22 = 4v v 22 Confirm that Av i = λ i v i 1 v 2 = 4
9 Fall / Dynamic Interpretation Since A = T ΛT 1, then e λ 1t T w 1 e At = T e Λt T 1 = v 1 v n.... e λ nt w T where we have written T w 1 T 1 =. w which is a column of rows. Multiply this expression out and we get that T n n e = e v i w i At λ i t T i=1 n Assume A diagonalizable, then ẋ = Ax, x(0) given, has solution x(t) = e At x(0) = T e Λt T 1 x(0) n = e λit T v i {w i x(0)} i=1 n λ i t = e v i β i i=1 State solution is linear combination of the system modes v i e λ it e λ it Determines nature of the time response v i Determines how each state contributes to that mode β i Determines extent to which initial condition excites the mode
10 Fall / Note that the v i give the relative sizing of the response of each part of the state vector to the response. 1 v 1 (t) = e t mode v 2 (t) = e 3t mode Clearly e λ it gives the time modulation λ i real growing/decaying exponential response λ i complex growing/decaying exponential damped sinusoidal Bottom line: The locations of the eigenvalues determine the pole locations for the system, thus: They determine the stability and/or performance & transient behavior of the system. It is their locations that we will want to modify when we start the control work
11 Fall / Diagonalization with Complex Roots If A has complex conjugate eigenvalues, the process is similar but a little more complicated. Consider a 2x2 case with A having eigenvalues a ± bi and associated eigenvectors e 1, e 2, with e 2 = ē 1. Then a + bi 0 1 A = e 1 e 2 e 1 e 2 0 a bi a + bi 0 1 = e 1 ē 1 T DT 1 ē 1 e 1 0 a bi Now use the transformation matrix 1 i M M = 0.5 = 1 i i i Then it follows that A = T DT 1 = (T M)(M 1 DM)(M 1 T 1 ) = (T M)(M 1 DM)(T M) 1 which has the nice structure: [ ] a b [ A = Re(e 1 ) Im(e 1 ) b a where all the matrices are real. Re(e 1 ) Im(e 1 ) ] 1 With complex roots, diagonalization is to a block diagonal form.
12 Fall / For this case we have that [ ] cos(bt) sin(bt) [ e At = Re(e 1 ) Im(e 1 ) e at sin(bt) cos(bt) ] Re(e 1 ) Im(e 1 ) 1 Note that [ ] Re(e 1 ) Im(e 1 ) Re(e 1 ) Im(e 1 ) = 0 1 So for an initial condition to excite just this mode, can pick x(0) = [Re(e 1 )], or x(0) = [Im(e 1 )] or a linear combination. Example x(0) = [Re(e 1 )] [ ] cos(bt) sin(bt) x(t) = e At x(0) = Re(e 1 ) Im(e 1 ) e at sin(bt) cos(bt) Re(e Im(e 1 ) 1 1 ) [Re(e 1 )] [ Im(e 1 ) ] e at cos(bt) sin(bt) 1 = Re(e 1 ) sin(bt) cos(bt) 0 [ ] cos(bt) = e at Re(e 1 ) Im(e 1 ) sin(bt) = e at (Re(e 1 ) cos(bt) Im(e 1 ) sin(bt)) which would ensure that only this mode is excited in the response
13 Fall / Example: Spring Mass System Classic example: spring mass system consider simple case first: m i = 1, and k i = 1 k 5 Wall k 1 k m 2 k 1 m 3 k 3 m 4 2 Wall x A K = = = [ z1 z 2 z 3 ż 1 ż 2 ż 3 ] [ ] 0 I M 1 M = diag(m K 0 i ) k 1 + k 2 + k 5 k 5 k 2 k 5 k 3 + k 4 + k 5 k 3 k 2 k 3 k 2 + k 3 Eigenvalues and eigenvectors of the undamped system λ 1 = ±0.77i λ 2 = ±1.85i λ 3 = ±2.00i v 1 v 2 v ±0.77i ±0.77i ±1.08i ±1.85i ±1.85i 2.61i ±2.00i 2.00i 0.00
14 Fall / Initial conditions to excite just the three modes: x i (0) = α 1 Re(v i ) + α 2 Im(v i ) α j R Simulation using α 1 = 1, α 2 = 0 Visualization important for correct physical interpretation Mode 1 λ 1 = ±0.77i m 1 m 3 m 2 Lowest frequency mode, all masses move in same direction Middle mass has higher amplitude motions z 3, motions all in phase
15 Fall / Mode 2 λ 2 = ±1.85i m 1 m 3 m 2 Middle frequency mode has middle mass moving in opposition to two end masses Again middle mass has higher amplitude motions z 3
16 Fall / Mode 3 λ 3 = ±2.00i m 1 m 3 m 2 Highest frequency mode, has middle mass stationary, and other two masses in opposition Eigenvectors that correspond to more constrained motion of the system are associated with higher frequency eigenvalues
17 Fall / Result if we use a random input is a combination of all three modes
18 Fall / Code: Simulation of Spring Mass System Modes % Simulate modal response for a spring mass system % Jonathan How, MIT % Fall 2009 alp1=1; % weighting choice on the IC m=eye(3); % mass k=[3 1 1; 1 3 1; 1 1 2]; % stiffness a=[m*0 eye(3); inv(m)*k m*0]; [v,d]=eig(a); t=[0:.01:20]; l1=1:25:length(t); G=ss(a,zeros(6,1),zeros(1,6),0); % use the following to cll the function above close all set(0, 'DefaultAxesFontSize', 14, 'DefaultAxesFontWeight','demi') set(0, 'DefaultTextFontSize', 14, 'DefaultTextFontWeight','demi') set(0,'defaultaxesfontname','arial'); set(0,'defaulttextfontname','arial');set(0,'defaultlinemarkersize',10) figure(1);clf x0=alp1*real(v(:,1))+(1 alp1)*imag(v(:,1)) [y,t,x]=lsim(g,0*t,t,x0); plot(t,x(:,1),' ','LineWidth',2);hold on plot(t(l1),x(l1,2),'rx','linewidth',2) plot(t,x(:,3),'m ','LineWidth',2);hold off xlabel('time');ylabel('displacement') title(['mode with \lambda=', num2str(imag(d(1,1))),' i']) legend('z1','z2','z3') print dpng r300 v1.png figure(2);clf x0=alp1*real(v(:,5))+(1 alp1)*imag(v(:,5)) [y,t,x]=lsim(g,0*t,t,x0); plot(t,x(:,1),' ','LineWidth',2);hold on plot(t(l1),x(l1,2),'rx','linewidth',2) plot(t,x(:,3),'m ','LineWidth',2);hold off xlabel('time');ylabel('displacement') title(['mode with \lambda=', num2str(imag(d(5,5))),' i']) legend('z1','z2','z3') print dpng r300 v3.png figure(3);clf x0=alp1*real(v(:,3))+(1 alp1)*imag(v(:,3)) [y,t,x]=lsim(g,0*t,t,x0); plot(t,x(:,1),' ','LineWidth',2);hold on plot(t,x(:,2),'r.','linewidth',2) plot(t,x(:,3),'m ','LineWidth',2);hold off xlabel('time');ylabel('displacement') title(['mode with \lambda=', num2str(imag(d(3,3))),' i']) legend('z1','z2','z3') print dpng r300 v2.png figure(4);clf x0=[ ]'; [y,t,x]=lsim(g,0*t,t,x0); plot(t,x(:,1),' ','LineWidth',2) hold on plot(t,x(:,2),'r.','linewidth',2) plot(t,x(:,3),'m ','LineWidth',2) hold off xlabel('time');ylabel('displacement') title(['random Initial Condition']) legend('z1','z2','z3') print dpng r300 v4.png
19 MIT OpenCourseWare / Feedback Control Systems Fall 2010 For information about citing these materials or our Terms of Use, visit:
Control Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli
Control Systems I Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 13, 2017 E. Frazzoli (ETH)
More 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 informationMTH 464: Computational Linear Algebra
MTH 464: Computational Linear Algebra Lecture Outlines Exam 4 Material Prof. M. Beauregard Department of Mathematics & Statistics Stephen F. Austin State University April 15, 2018 Linear Algebra (MTH 464)
More informationModal Decomposition and the Time-Domain Response of Linear Systems 1
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING.151 Advanced System Dynamics and Control Modal Decomposition and the Time-Domain Response of Linear Systems 1 In a previous handout
More informationChapter 7. Canonical Forms. 7.1 Eigenvalues and Eigenvectors
Chapter 7 Canonical Forms 7.1 Eigenvalues and Eigenvectors Definition 7.1.1. Let V be a vector space over the field F and let T be a linear operator on V. An eigenvalue of T is a scalar λ F such that there
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 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 informationEE263: Introduction to Linear Dynamical Systems Review Session 5
EE263: Introduction to Linear Dynamical Systems Review Session 5 Outline eigenvalues and eigenvectors diagonalization matrix exponential EE263 RS5 1 Eigenvalues and eigenvectors we say that λ C is an eigenvalue
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 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 informationControl Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 5: Transfer Functions Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 20, 2017 E. Frazzoli (ETH) Lecture 5: Control Systems I 20/10/2017
More informationẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More 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 information22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes
Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration
More 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 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 informationand let s calculate the image of some vectors under the transformation T.
Chapter 5 Eigenvalues and Eigenvectors 5. Eigenvalues and Eigenvectors Let T : R n R n be a linear transformation. Then T can be represented by a matrix (the standard matrix), and we can write T ( v) =
More informationEigenvalues and Eigenvectors
Eigenvalues and Eigenvectors week -2 Fall 26 Eigenvalues and eigenvectors The most simple linear transformation from R n to R n may be the transformation of the form: T (x,,, x n ) (λ x, λ 2,, λ n x n
More informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More 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 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 informationDiagonalization. P. Danziger. u B = A 1. B u S.
7., 8., 8.2 Diagonalization P. Danziger Change of Basis Given a basis of R n, B {v,..., v n }, we have seen that the matrix whose columns consist of these vectors can be thought of as a change of basis
More informationModes and Roots ... mx + bx + kx = 0. (2)
A solution of the form x(t) = ce rt to the homogeneous constant coefficient linear equation a x (n) + (n 1). n a n 1 x + + a 1 x + a 0 x = 0 (1) is called a modal solution and ce rt is called a mode of
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More information1. The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ = Ax + Bu is
ECE 55, Fall 2007 Problem Set #4 Solution The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ Ax + Bu is x(t) e A(t ) x( ) + e A(t τ) Bu(τ)dτ () This formula is extremely important
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 informationSection 8.2 : Homogeneous Linear Systems
Section 8.2 : Homogeneous Linear Systems Review: Eigenvalues and Eigenvectors Let A be an n n matrix with constant real components a ij. An eigenvector of A is a nonzero n 1 column vector v such that Av
More informationEigenvalues, Eigenvectors. Eigenvalues and eigenvector will be fundamentally related to the nature of the solutions of state space systems.
Chapter 3 Linear Algebra In this Chapter we provide a review of some basic concepts from Linear Algebra which will be required in order to compute solutions of LTI systems in state space form, discuss
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 interpretation of eigenvectors
EE263 Autumn 2015 S. Boyd and S. Lall Dynamic interpretation of eigenvectors invariant sets complex eigenvectors & invariant planes left eigenvectors modal form discrete-time stability 1 Dynamic interpretation
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 information9.6: Matrix Exponential, Repeated Eigenvalues. Ex.: A = x 1 (t) = e t 2 F.M.: If we set
9.6: Matrix Exponential, Repeated Eigenvalues x Ax, A : n n (1) Def.: If x 1 (t),...,x n (t) is a fundamental set of solutions (F.S.S.) of (1), then X(t) x 1 (t),...,x n (t) (n n) is called a fundamental
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 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 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 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 informationGeneralized Eigenvectors and Jordan Form
Generalized Eigenvectors and Jordan Form We have seen that an n n matrix A is diagonalizable precisely when the dimensions of its eigenspaces sum to n. So if A is not diagonalizable, there is at least
More informationMatrix Solutions to Linear Systems of ODEs
Matrix Solutions to Linear Systems of ODEs James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 3, 216 Outline 1 Symmetric Systems of
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 information1 Continuous-time 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 informationProblem 2 (Gaussian Elimination, Fundamental Spaces, Least Squares, Minimum Norm) Consider the following linear algebraic system of equations:
EEE58 Exam, Fall 6 AA Rodriguez Rules: Closed notes/books, No calculators permitted, open minds GWC 35, 965-37 Problem (Dynamic Augmentation: State Space Representation) Consider a dynamical system consisting
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 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 informationTopic # Feedback Control
Topic #11 16.31 Feedback Control State-Space Systems State-space model features Observability Controllability Minimal Realizations Copyright 21 by Jonathan How. 1 Fall 21 16.31 11 1 State-Space Model Features
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 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 informationSystems of Linear Differential Equations Chapter 7
Systems of Linear Differential Equations Chapter 7 Doreen De Leon Department of Mathematics, California State University, Fresno June 22, 25 Motivating Examples: Applications of Systems of First Order
More informationSOLVING DIFFERENTIAL EQUATIONS. Amir Asif. Department of Computer Science and Engineering York University, Toronto, ON M3J1P3
SOLVING DIFFERENTIAL EQUATIONS Amir Asif Department of Computer Science and Engineering York University, Toronto, ON M3J1P3 ABSTRACT This article reviews a direct method for solving linear, constant-coefficient
More informationLinear Algebra and ODEs review
Linear Algebra and ODEs review Ania A Baetica September 9, 015 1 Linear Algebra 11 Eigenvalues and eigenvectors Consider the square matrix A R n n (v, λ are an (eigenvector, eigenvalue pair of matrix A
More informationπ 1 = tr(a), π n = ( 1) n det(a). In particular, when n = 2 one has
Eigen Methods Math 246, Spring 2009, Professor David Levermore Eigenpairs Let A be a real n n matrix A number λ possibly complex is an eigenvalue of A if there exists a nonzero vector v possibly complex
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 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 information16.30/31, Fall 2010 Recitation # 13
16.30/31, Fall 2010 Recitation # 13 Brandon Luders December 6, 2010 In this recitation, we tie the ideas of Lyapunov stability analysis (LSA) back to previous ways we have demonstrated stability - but
More informationRemark By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.
Sec 6 Eigenvalues and Eigenvectors Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called an eigenvalue of A if there is a nontrivial
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 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 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 informationJordan Canonical Form Homework Solutions
Jordan Canonical Form Homework Solutions For each of the following, put the matrix in Jordan canonical form and find the matrix S such that S AS = J. [ ]. A = A λi = λ λ = ( λ) = λ λ = λ =, Since we have
More informationExam in TMA4110 Calculus 3, June 2013 Solution
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 8 Exam in TMA4 Calculus 3, June 3 Solution Problem Let T : R 3 R 3 be a linear transformation such that T = 4,
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 informationDifferential equations
Differential equations Math 7 Spring Practice problems for April Exam Problem Use the method of elimination to find the x-component of the general solution of x y = 6x 9x + y = x 6y 9y Soln: The system
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 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 informationNOTES ON LINEAR ODES
NOTES ON LINEAR ODES JONATHAN LUK We can now use all the discussions we had on linear algebra to study linear ODEs Most of this material appears in the textbook in 21, 22, 23, 26 As always, this is a preliminary
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.24: Dynamic Systems Spring 20 Homework 9 Solutions Exercise 2. We can use additive perturbation model with
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 information16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1
16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1 Charles P. Coleman October 31, 2005 1 / 40 : Controllability Tests Observability Tests LEARNING OUTCOMES: Perform controllability tests Perform
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 informationLinear System Theory
Linear System Theory Wonhee Kim Lecture 4 Apr. 4, 2018 1 / 40 Recap Vector space, linear space, linear vector space Subspace Linearly independence and dependence Dimension, Basis, Change of Basis 2 / 40
More informationEigenpairs and Diagonalizability Math 401, Spring 2010, Professor David Levermore
Eigenpairs and Diagonalizability Math 40, Spring 200, Professor David Levermore Eigenpairs Let A be an n n matrix A number λ possibly complex even when A is real is an eigenvalue of A if there exists a
More informationLS.2 Homogeneous Linear Systems with Constant Coefficients
LS2 Homogeneous Linear Systems with Constant Coefficients Using matrices to solve linear systems The naive way to solve a linear system of ODE s with constant coefficients is by eliminating variables,
More informationLinear Algebra- Final Exam Review
Linear Algebra- Final Exam Review. Let A be invertible. Show that, if v, v, v 3 are linearly independent vectors, so are Av, Av, Av 3. NOTE: It should be clear from your answer that you know the definition.
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 informationReduction to the associated homogeneous system via a particular solution
June PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 5) Linear Algebra This study guide describes briefly the course materials to be covered in MA 5. In order to be qualified for the credit, one
More informationMatrices and Linear transformations
Matrices and Linear transformations We have been thinking of matrices in connection with solutions to linear systems of equations like Ax = b. It is time to broaden our horizons a bit and start thinking
More informationCAAM 335 Matrix Analysis
CAAM 335 Matrix Analysis Solutions to Homework 8 Problem (5+5+5=5 points The partial fraction expansion of the resolvent for the matrix B = is given by (si B = s } {{ } =P + s + } {{ } =P + (s (5 points
More informationPart II Problems and Solutions
Problem 1: [Complex and repeated eigenvalues] (a) The population of long-tailed weasels and meadow voles on Nantucket Island has been studied by biologists They measure the populations relative to a baseline,
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 informationUnderstanding the Matrix Exponential
Transformations Understanding the Matrix Exponential Lecture 8 Math 634 9/17/99 Now that we have a representation of the solution of constant-coefficient initial-value problems, we should ask ourselves:
More informationSpring 2019 Exam 2 3/27/19 Time Limit: / Problem Points Score. Total: 280
Math 307 Spring 2019 Exam 2 3/27/19 Time Limit: / Name (Print): Problem Points Score 1 15 2 20 3 35 4 30 5 10 6 20 7 20 8 20 9 20 10 20 11 10 12 10 13 10 14 10 15 10 16 10 17 10 Total: 280 Math 307 Exam
More informationSolution: In standard form (i.e. y + P (t)y = Q(t)) we have y t y = cos(t)
Math 380 Practice Final Solutions This is longer than the actual exam, which will be 8 to 0 questions (some might be multiple choice). You are allowed up to two sheets of notes (both sides) and a calculator,
More informationEigenvalues, Eigenvectors, and Diagonalization
Math 240 TA: Shuyi Weng Winter 207 February 23, 207 Eigenvalues, Eigenvectors, and Diagonalization The concepts of eigenvalues, eigenvectors, and diagonalization are best studied with examples. We will
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More informationUniversity of Toronto Department of Electrical and Computer Engineering ECE410F Control Systems Problem Set #3 Solutions = Q o = CA.
University of Toronto Department of Electrical and Computer Engineering ECE41F Control Systems Problem Set #3 Solutions 1. The observability matrix is Q o C CA 5 6 3 34. Since det(q o ), the matrix is
More informationHomogeneous Linear Systems of Differential Equations with Constant Coefficients
Objective: Solve Homogeneous Linear Systems of Differential Equations with Constant Coefficients dx a x + a 2 x 2 + + a n x n, dx 2 a 2x + a 22 x 2 + + a 2n x n,. dx n = a n x + a n2 x 2 + + a nn x n.
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 informationDIAGONALIZATION. In order to see the implications of this definition, let us consider the following example Example 1. Consider the matrix
DIAGONALIZATION Definition We say that a matrix A of size n n is diagonalizable if there is a basis of R n consisting of eigenvectors of A ie if there are n linearly independent vectors v v n such that
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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering 2.04A Systems and Controls Spring 2013
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering 2.04A Systems and Controls Spring 2013 Problem Set #4 Posted: Thursday, Mar. 7, 13 Due: Thursday, Mar. 14, 13 1. Sketch the Root
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 informationSymmetric and anti symmetric matrices
Symmetric and anti symmetric matrices In linear algebra, a symmetric matrix is a square matrix that is equal to its transpose. Formally, matrix A is symmetric if. A = A Because equal matrices have equal
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
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 information1. In this problem, if the statement is always true, circle T; otherwise, circle F.
Math 1553, Extra Practice for Midterm 3 (sections 45-65) Solutions 1 In this problem, if the statement is always true, circle T; otherwise, circle F a) T F If A is a square matrix and the homogeneous equation
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 informationRemark 1 By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.
Sec 5 Eigenvectors and Eigenvalues In this chapter, vector means column vector Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called
More informationAPPPHYS 217 Tuesday 6 April 2010
APPPHYS 7 Tuesday 6 April Stability and input-output performance: second-order systems Here we present a detailed example to draw connections between today s topics and our prior review of linear algebra
More 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 informationMath 3191 Applied Linear Algebra
Math 9 Applied Linear Algebra Lecture 9: Diagonalization Stephen Billups University of Colorado at Denver Math 9Applied Linear Algebra p./9 Section. Diagonalization The goal here is to develop a useful
More information