Discrete Time Systems:
|
|
- Gwendoline Harris
- 6 years ago
- Views:
Transcription
1 Discrete Time Systems: Finding the control sequence {u(k): k =, 1, n} that will drive a discrete time system from some initial condition x() to some terminal condition x(k). Given the discrete time system (DTS) whose state equation is given by: x(k+1) = Ax(k) + Bu(k): (1) with the solution given as: k 1 k k j 1 xk ( ) = Ax( ) + A Buj ( ) j= (2) we wish to determine how to find the sequence {u(l): l =, 1, k} that drives the system from an arbitrary initial condition x() to a terminal condition x(k). Equation (2) contains both the transient and forced response of the system. If we define an initial condition x() and a terminal condition x(k) and expand the convolution summation for some final value of k, the equation may be written in vector matrix form as: [x(k) - A k x()] = [B, AB, A 2 B,. A K-1 B] [u(k-1) u(k-2) u(1) u()] T () Copyright Dr. Thomas Chmielewski 27 Page 1 of 1
2 The input or forcing function sequence {u(k)} may be found by solving vector matrix equation (). A key question to be answered is how big does k have to be to drive the system from some initial condition to a terminal condition. Alternatively we may ask if there are any limits to the size of k. At this point it may be obvious that the answer lies in the controllability test matrix : C T = [B, AB, A 2 B, A k-1 B] (4) Note answering the first question also allows us to understand controllability in a better sense. Rewrite Eq. (2) as follows: [B, AB, A 2 B,. A K-1 B] [u(k-1); u(k-2); u(1) u()] = [x(k) - A k x()] (2a) consider Eq. (2a) to be in the form C T U = Y (5) where: C T : n by m U: m by 1 Y: n by 1 Copyright Dr. Thomas Chmielewski 27 Page 2 of 1
3 From basic linear algebra, there exists a solution for U (possibly more than one), if and only if Y can be expressed as a linear combination of the columns of C T. The equivalent rank calculation for the existence of a solution of U is given as Rank [C T Y] = Rank [C T ] (6) that is by forming an augmented system, the addition of the column y does not change the rank. Since in our situation Y corresponds to an n by 1 vector (since x() and x(k) are n by 1 and a must be square and of order n) then a solution for U only exists if the partitioned matrix has rank n. Alternatively there must be n linearly independent columns out of C T. This allows us to put a bound on the size of C T and therefore the size of the vector U. C T contains n rows, based on the number of states, but could have as many columns as desired (by Eq. 2 and ). By the Cayley Hamilton theorem powers of a matrix A whose order is n by n can be expressed as A m = f(i, A, A 2, A n-1 ) where m n. Thus the additional partitions BA N provide no new linearly independent columns. In summary, if the rank of Equation (4) equals the order of the system, n, then the system is completely controllable and a control sequence {u(k): k =, 1, n-1} exists that will drive the system from some arbitrary initial condition to the origin or any other desired terminal state. Copyright Dr. Thomas Chmielewski 27 Page of 1
4 Example 1 : driving a system to the origin in k = n steps Given the system described by the matrices A = [ 1; 2] B = [ ;] x2 = A = B = 1 2 with initial condition: x= [6;] x = 6 and terminal condition x2 = [ ;] The solution is computed from the relationship: [x(2) - A 2 x()] = [B AB] [u(1); u()] yielding [u(1); u()] = [B AB] -1 [x(2) - A 2 x()] Using Matlab to evaluate the expression yields: u = inv([b, A*B])*(x2 - A^2*x) u = -2 Copyright Dr. Thomas Chmielewski 27 Page 4 of 1
5 Note that u() = -2 and u(1) = (this is the reverse order from the natural indexing for the vector u) The best way to validate the result is by recursive substitution at k = x1 = A x + Bu() x1 = A*x + B*(-2) x1 = at k = 1 x2 = A x1 + Bu(1) x2 = A*x1 + B*() x2 = which is the desired result. Copyright Dr. Thomas Chmielewski 27 Page 5 of 1
6 Example 2 : driving a system to an arbitrary location in k = n steps Now let us consider driving the system to the terminal condition x2 = [1;1] x2 = 1 1 The solution is the same with the new numerical value for x2 u = inv([b, A*B])*(x2 - A^2*x) u = Remember to modify the vector u so that the indices match the discrete time equation, Validating the solution: at k = x1 = A x + Bu() x1 = A*x + B*1. x1 = at k = 1 x2 = A x1 + Bu(1) x2 = A*x1 + B*(-.) x2 = Copyright Dr. Thomas Chmielewski 27 Page 6 of 1
7 which is the desired result within the precision of the numbers used 12 It is interesting to plot the trajectory in state space. We define a time history of the states in a trajectory matrix whose columns correspond to the state at a given time. traj = [x, x1, x2] traj = now we will plot the second state (i.e., second row of traj) versus the first state, using arrows to show the direction of the trajectory. plot(traj(1, : ), traj(2, : ), '->') axis([, 12,, 12]) The trajectory plot shows the two transitions necessary to drive the initial state to the terminal state Copyright Dr. Thomas Chmielewski 27 Page 7 of 1
8 Using more than n steps to drive the system from x() to x(k) Consider the case for k > n steps, the equation to find {u(k)} i is the same as before: [x(k) - A k x()] = [B, AB, A 2 B,. A K-1 B] [u(k-1); u(k-2); u(1) u()] (2) here the controllability test matrix C T = [B, AB, A 2 B,. A K-1 B] is no longer terminated at a power of n-1 but contains the number of partitions equals that of the desired number of steps. By the Cayley Hamilton Theorem, the extra columns are linearly dependent on the first n columns With these additional columns, C T is no longer square and the inverse is undefined. This is called the underspecified case. The solution of Eq. (5) and utilizes a least squares solution. C T U = Y (5) Since there are more unknowns k than equations, n, there are many solutions, we will specify the solution, U, that has the minimum norm. The minimum norm solution uses the pseudo inverse defined for a general non square matrix as: [C T T (C T C T T ) -1 ]. Hence U = [C T T (C T C T T ) -1 ] Y (7) The sizes of the matrices in Eq. (7), these are given in the ordinal position of the matrices of Eq. (7) as:(k by 1) = [(k by n) [(n by k)(k by n)] (n by 1) Copyright Dr. Thomas Chmielewski 27 Page 8 of 1
9 Example : driving a system to the origin in k = steps Using the system of Example 1, we wish to drive the system to the origin in three steps. Computing the pseudo inverse yieids: CT = [B, A*B, A^2*B] A = [ 1; 2] B = [ ; ] x= [6;] x = [ ;] A = B = x = x = CT = P_INV = CT'*inv(CT*CT') P_INV = Solving for vector U yields : U = P_INV*(x - A^*x) U = Copyright Dr. Thomas Chmielewski 27 Page 9 of 1
10 Once again we will verify via recursive substitution: x1 = A*x + B*(-1.6) x1 =. 1.2 x2 = A*x1 + B*(-.8) x2 = traj = [x, x1, x2, x] traj = Plotting the trajectory as before: plot(traj(1, : ), traj(2, : ), '-*') axis([-1, 7, -1, 7]) x = A*x2 + B*() x = 1.e-14 * x is essentially zero x = [;] x = Copyright Dr. Thomas Chmielewski 27 Page 1 of 1
11 Example 4 : driving a system to the origin in k = 5 steps Here we find 5 inputs that will drive the system from the initial condition to the origin CT = [B, A*B, A^2*B, A^*B, A^4*B] CT = x5 = [;] x5 = P_INV = CT' * inv(ct * CT') P_INV = Copyright Dr. Thomas Chmielewski 27 Page 11 of 1
12 U = P_INV*(x5 - A^5*x) U = x1 = A*x + B*(-1.559) x1 = x2 = A*x1 + B*(-.7529) x2 = x = A*x2 + B*(-.765) x = x4 = A*x + B*(-.1882) x4 =.282. x5 = A*x4 + B*() x5 = 1.e-14 * so again x5 is essentially zero, the small numbers are due to floating point implementation used in Matlab traj = [x, x1, x2, x, x4, x5] traj = Copyright Dr. Thomas Chmielewski 27 Page 12 of 1
13 plot(traj(1,: ), traj(2, : ), '-*') axis([-1, 7, -1, 7]) Copyright Dr. Thomas Chmielewski 27 Page 1 of 1
Solving 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 informationModule 9: State Feedback Control Design Lecture Note 1
Module 9: State Feedback Control Design Lecture Note 1 The design techniques described in the preceding lectures are based on the transfer function of a system. In this lecture we would discuss the state
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 information7.5 Operations with Matrices. Copyright Cengage Learning. All rights reserved.
7.5 Operations with Matrices Copyright Cengage Learning. All rights reserved. What You Should Learn Decide whether two matrices are equal. Add and subtract matrices and multiply matrices by scalars. Multiply
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationLecture 4 and 5 Controllability and Observability: Kalman decompositions
1 Lecture 4 and 5 Controllability and Observability: Kalman decompositions Spring 2013 - EE 194, Advanced Control (Prof. Khan) January 30 (Wed.) and Feb. 04 (Mon.), 2013 I. OBSERVABILITY OF DT LTI SYSTEMS
More informationMULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.
est Review-Linear Algebra Name MULIPLE CHOICE Choose the one alternative that best completes the statement or answers the question Solve the system of equations ) 7x + 7 + x + + 9x + + 9 9 (-,, ) (, -,
More informationChap 3. Linear Algebra
Chap 3. Linear Algebra Outlines 1. Introduction 2. Basis, Representation, and Orthonormalization 3. Linear Algebraic Equations 4. Similarity Transformation 5. Diagonal Form and Jordan Form 6. Functions
More informationRural/Urban Migration: The Dynamics of Eigenvectors
* Analysis of the Dynamic Structure of a System * Rural/Urban Migration: The Dynamics of Eigenvectors EGR 326 April 11, 2019 1. Develop the system model and create the Matlab/Simulink model 2. Plot and
More informationMATH 2331 Linear Algebra. Section 2.1 Matrix Operations. Definition: A : m n, B : n p. Example: Compute AB, if possible.
MATH 2331 Linear Algebra Section 2.1 Matrix Operations Definition: A : m n, B : n p ( 1 2 p ) ( 1 2 p ) AB = A b b b = Ab Ab Ab Example: Compute AB, if possible. 1 Row-column rule: i-j-th entry of AB:
More informationLecture 7 and 8. Fall EE 105, Feedback Control Systems (Prof. Khan) September 30 and October 05, 2015
1 Lecture 7 and 8 Fall 2015 - EE 105, Feedback Control Systems (Prof Khan) September 30 and October 05, 2015 I CONTROLLABILITY OF AN DT-LTI SYSTEM IN k TIME-STEPS The DT-LTI system is given by the following
More informationMath 4377/6308 Advanced Linear Algebra
2. Linear Transformations Math 4377/638 Advanced Linear Algebra 2. Linear Transformations, Null Spaces and Ranges Jiwen He Department of Mathematics, University of Houston jiwenhe@math.uh.edu math.uh.edu/
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationSolution for Homework 5
Solution for Homework 5 ME243A/ECE23A Fall 27 Exercise 1 The computation of the reachable subspace in continuous time can be handled easily introducing the concepts of inner product, orthogonal complement
More informationCase Study: Space Flight and Control Systems
Case Study: Space Flight and Control Systems In this case study, we examine how concepts in Chapter 4 may be used in the design of engineering control systems, such as the one in Figure on page of your
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 informationChapter 8: Linear Algebraic Equations
Chapter 8: Linear Algebraic Equations Matrix Methods for Linear Equations Uniqueness and Existence of Solutions Under-Determined Systems Over-Determined Systems Linear Algebraic Equations For 2 Equations
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 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 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 informationMATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem
MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem Pulemotov, September 12, 2012 Unit Outline Goal 1: Outline linear
More informationMATH 2331 Linear Algebra. Section 1.1 Systems of Linear Equations. Finding the solution to a set of two equations in two variables: Example 1: Solve:
MATH 2331 Linear Algebra Section 1.1 Systems of Linear Equations Finding the solution to a set of two equations in two variables: Example 1: Solve: x x = 3 1 2 2x + 4x = 12 1 2 Geometric meaning: Do these
More informationModelling and Mathematical Methods in Process and Chemical Engineering
FS 07 February 0, 07 Modelling and Mathematical Methods in Process and Chemical Engineering Solution Series. Systems of linear algebraic equations: adjoint, determinant, inverse The adjoint of a (square)
More informationReview Questions REVIEW QUESTIONS 71
REVIEW QUESTIONS 71 MATLAB, is [42]. For a comprehensive treatment of error analysis and perturbation theory for linear systems and many other problems in linear algebra, see [126, 241]. An overview of
More informationSolutions of Linear system, vector and matrix equation
Goals: Solutions of Linear system, vector and matrix equation Solutions of linear system. Vectors, vector equation. Matrix equation. Math 112, Week 2 Suggested Textbook Readings: Sections 1.3, 1.4, 1.5
More information6-2 Matrix Multiplication, Inverses and Determinants
Find AB and BA, if possible. 1. A = A = ; A is a 1 2 matrix and B is a 2 2 matrix. Because the number of columns of A is equal to the number of rows of B, AB exists. To find the first entry of AB, find
More informationMTH 5102 Linear Algebra Practice Final Exam April 26, 2016
Name (Last name, First name): MTH 5 Linear Algebra Practice Final Exam April 6, 6 Exam Instructions: You have hours to complete the exam. There are a total of 9 problems. You must show your work and write
More informationANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3
ANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3 ISSUED 24 FEBRUARY 2018 1 Gaussian elimination Let A be an (m n)-matrix Consider the following row operations on A (1) Swap the positions any
More informationMA 242 LINEAR ALGEBRA C1, Solutions to First Midterm Exam
MA 242 LINEAR ALGEBRA C Solutions to First Midterm Exam Prof Nikola Popovic October 2 9:am - :am Problem ( points) Determine h and k such that the solution set of x + = k 4x + h = 8 (a) is empty (b) contains
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 informationMath 3191 Applied Linear Algebra
Math 9 Applied Linear Algebra Lecture : Null and Column Spaces Stephen Billups University of Colorado at Denver Math 9Applied Linear Algebra p./8 Announcements Study Guide posted HWK posted Math 9Applied
More informationMATH10212 Linear Algebra B Homework Week 4
MATH22 Linear Algebra B Homework Week 4 Students are strongly advised to acquire a copy of the Textbook: D. C. Lay Linear Algebra and its Applications. Pearson, 26. ISBN -52-2873-4. Normally, homework
More informationMATH10212 Linear Algebra B Homework Week 3. Be prepared to answer the following oral questions if asked in the supervision class
MATH10212 Linear Algebra B Homework Week Students are strongly advised to acquire a copy of the Textbook: D. C. Lay Linear Algebra its Applications. Pearson, 2006. ISBN 0-521-2871-4. Normally, homework
More informationMatrix operations Linear Algebra with Computer Science Application
Linear Algebra with Computer Science Application February 14, 2018 1 Matrix operations 11 Matrix operations If A is an m n matrix that is, a matrix with m rows and n columns then the scalar entry in the
More informationLinear Equation: a 1 x 1 + a 2 x a n x n = b. x 1, x 2,..., x n : variables or unknowns
Linear Equation: a x + a 2 x 2 +... + a n x n = b. x, x 2,..., x n : variables or unknowns a, a 2,..., a n : coefficients b: constant term Examples: x + 4 2 y + (2 5)z = is linear. x 2 + y + yz = 2 is
More informationNOTES ON MATRICES OF FULL COLUMN (ROW) RANK. Shayle R. Searle ABSTRACT
NOTES ON MATRICES OF FULL COLUMN (ROW) RANK Shayle R. Searle Biometrics Unit, Cornell University, Ithaca, N.Y. 14853 BU-1361-M August 1996 ABSTRACT A useful left (right) inverse of a full column (row)
More informationSolution: By inspection, the standard matrix of T is: A = Where, Ae 1 = 3. , and Ae 3 = 4. , Ae 2 =
This is a typical assignment, but you may not be familiar with the material. You should also be aware that many schools only give two exams, but also collect homework which is usually worth a small part
More informationLecture 2 Discrete-Time LTI Systems: Introduction
Lecture 2 Discrete-Time LTI Systems: Introduction Outline 2.1 Classification of Systems.............................. 1 2.1.1 Memoryless................................. 1 2.1.2 Causal....................................
More information4F3 - Predictive Control
4F3 Predictive Control - Lecture 3 p 1/21 4F3 - Predictive Control Lecture 3 - Predictive Control with Constraints Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 3 p 2/21 Constraints on
More informationJim Lambers MAT 610 Summer Session Lecture 1 Notes
Jim Lambers MAT 60 Summer Session 2009-0 Lecture Notes Introduction This course is about numerical linear algebra, which is the study of the approximate solution of fundamental problems from linear algebra
More informationDesigning Information Devices and Systems II Fall 2015 Note 22
EE 16B Designing Information Devices and Systems II Fall 2015 Note 22 Notes taken by John Noonan (11/12) Graphing of the State Solutions Open loop x(k + 1) = Ax(k) + Bu(k) y(k) = Cx(k) Closed loop x(k
More informationSpan & Linear Independence (Pop Quiz)
Span & Linear Independence (Pop Quiz). Consider the following vectors: v = 2, v 2 = 4 5, v 3 = 3 2, v 4 = Is the set of vectors S = {v, v 2, v 3, v 4 } linearly independent? Solution: Notice that the number
More information4F3 - Predictive Control
4F3 Predictive Control - Lecture 2 p 1/23 4F3 - Predictive Control Lecture 2 - Unconstrained Predictive Control Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 2 p 2/23 References Predictive
More informationOutline. 1 Linear Quadratic Problem. 2 Constraints. 3 Dynamic Programming Solution. 4 The Infinite Horizon LQ Problem.
Model Predictive Control Short Course Regulation James B. Rawlings Michael J. Risbeck Nishith R. Patel Department of Chemical and Biological Engineering Copyright c 217 by James B. Rawlings Outline 1 Linear
More informationMATH 2030: ASSIGNMENT 4 SOLUTIONS
MATH 23: ASSIGNMENT 4 SOLUTIONS More on the LU factorization Q.: pg 96, q 24. Find the P t LU factorization of the matrix 2 A = 3 2 2 A.. By interchanging row and row 4 we get a matrix that may be easily
More informationIf A is a 4 6 matrix and B is a 6 3 matrix then the dimension of AB is A. 4 6 B. 6 6 C. 4 3 D. 3 4 E. Undefined
Question 1 If A is a 4 6 matrix and B is a 6 3 matrix then the dimension of AB is A. 4 6 B. 6 6 C. 4 3 D. 3 4 E. Undefined Quang T. Bach Math 18 October 18, 2017 1 / 17 Question 2 1 2 Let A = 3 4 1 2 3
More informationSection 4.5. Matrix Inverses
Section 4.5 Matrix Inverses The Definition of Inverse Recall: The multiplicative inverse (or reciprocal) of a nonzero number a is the number b such that ab = 1. We define the inverse of a matrix in almost
More informationLinear Algebra. Solving SLEs with Matlab. Matrix Inversion. Solving SLE s by Matlab - Inverse. Solving Simultaneous Linear Equations in MATLAB
Linear Algebra Solving Simultaneous Linear Equations in MATLAB Solving SLEs with Matlab Matlab can solve some numerical SLE s A b Five techniques available:. method. method. method 4. method. method Matri
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 informationCSL361 Problem set 4: Basic linear algebra
CSL361 Problem set 4: Basic linear algebra February 21, 2017 [Note:] If the numerical matrix computations turn out to be tedious, you may use the function rref in Matlab. 1 Row-reduced echelon matrices
More informationLet x be an approximate solution for Ax = b, e.g., obtained by Gaussian elimination. Let x denote the exact solution. Call. r := b A x.
ESTIMATION OF ERROR Let x be an approximate solution for Ax = b, e.g., obtained by Gaussian elimination. Let x denote the exact solution. Call the residual for x. Then r := b A x r = b A x = Ax A x = A
More informationSYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER
SYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER Exercise 54 Consider the system: ẍ aẋ bx u where u is the input and x the output signal (a): Determine a state space realization (b): Is the
More informationEE363 Review Session 1: LQR, Controllability and Observability
EE363 Review Session : LQR, Controllability and Observability In this review session we ll work through a variation on LQR in which we add an input smoothness cost, in addition to the usual penalties on
More informationJune 2011 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations
June 20 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations The topics covered in this exam can be found in An introduction to differential equations
More information3.2 Systems of Two First Order Linear DE s
Agenda Section 3.2 Reminders Lab 1 write-up due 9/26 or 9/28 Lab 2 prelab due 9/26 or 9/28 WebHW due 9/29 Office hours Tues, Thurs 1-2 pm (5852 East Hall) MathLab office hour Sun 7-8 pm (MathLab) 3.2 Systems
More informationMODULE 8 Topics: Null space, range, column space, row space and rank of a matrix
MODULE 8 Topics: Null space, range, column space, row space and rank of a matrix Definition: Let L : V 1 V 2 be a linear operator. The null space N (L) of L is the subspace of V 1 defined by N (L) = {x
More informationSAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear Algebra
1.1. Introduction SAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear algebra is a specific branch of mathematics dealing with the study of vectors, vector spaces with functions that
More informationMath Camp II. Basic Linear Algebra. Yiqing Xu. Aug 26, 2014 MIT
Math Camp II Basic Linear Algebra Yiqing Xu MIT Aug 26, 2014 1 Solving Systems of Linear Equations 2 Vectors and Vector Spaces 3 Matrices 4 Least Squares Systems of Linear Equations Definition A linear
More information2018 Fall 2210Q Section 013 Midterm Exam I Solution
8 Fall Q Section 3 Midterm Exam I Solution True or False questions ( points = points) () An example of a linear combination of vectors v, v is the vector v. True. We can write v as v + v. () If two matrices
More informationElementary maths for GMT
Elementary maths for GMT Linear Algebra Part 2: Matrices, Elimination and Determinant m n matrices The system of m linear equations in n variables x 1, x 2,, x n a 11 x 1 + a 12 x 2 + + a 1n x n = b 1
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 informationMethods for Solving Linear Systems Part 2
Methods for Solving Linear Systems Part 2 We have studied the properties of matrices and found out that there are more ways that we can solve Linear Systems. In Section 7.3, we learned that we can use
More informationObservability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)
Observability Dynamic Systems Lecture 2 Observability Continuous time model: Discrete time model: ẋ(t) = f (x(t), u(t)), y(t) = h(x(t), u(t)) x(t + 1) = f (x(t), u(t)), y(t) = h(x(t)) Reglerteknik, ISY,
More informationAdvanced Engineering Mathematics Prof. Pratima Panigrahi Department of Mathematics Indian Institute of Technology, Kharagpur
Advanced Engineering Mathematics Prof. Pratima Panigrahi Department of Mathematics Indian Institute of Technology, Kharagpur Lecture No. #07 Jordan Canonical Form Cayley Hamilton Theorem (Refer Slide Time:
More informationSuppose that we have a specific single stage dynamic system governed by the following equation:
Dynamic Optimisation Discrete Dynamic Systems A single stage example Suppose that we have a specific single stage dynamic system governed by the following equation: x 1 = ax 0 + bu 0, x 0 = x i (1) where
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 Combination. v = a 1 v 1 + a 2 v a k v k
Linear Combination Definition 1 Given a set of vectors {v 1, v 2,..., v k } in a vector space V, any vector of the form v = a 1 v 1 + a 2 v 2 +... + a k v k for some scalars a 1, a 2,..., a k, is called
More informationControllability. Chapter Reachable States. This chapter develops the fundamental results about controllability and pole assignment.
Chapter Controllability This chapter develops the fundamental results about controllability and pole assignment Reachable States We study the linear system ẋ = Ax + Bu, t, where x(t) R n and u(t) R m Thus
More informationElementary Matrices. MATH 322, Linear Algebra I. J. Robert Buchanan. Spring Department of Mathematics
Elementary Matrices MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Outline Today s discussion will focus on: elementary matrices and their properties, using elementary
More informationChapter 5. Linear Algebra. Sections A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra Sections 5.1 5.3 A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are
More informationFor GLM y = Xβ + e (1) where X is a N k design matrix and p(e) = N(0, σ 2 I N ), we can estimate the coefficients from the normal equations
1 Generalised Inverse For GLM y = Xβ + e (1) where X is a N k design matrix and p(e) = N(0, σ 2 I N ), we can estimate the coefficients from the normal equations (X T X)β = X T y (2) If rank of X, denoted
More informationMTH501- Linear Algebra MCQS MIDTERM EXAMINATION ~ LIBRIANSMINE ~
MTH501- Linear Algebra MCQS MIDTERM EXAMINATION ~ LIBRIANSMINE ~ Question No: 1 (Marks: 1) If for a linear transformation the equation T(x) =0 has only the trivial solution then T is One-to-one Onto Question
More informationSAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear Algebra
SAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to 1.1. Introduction Linear algebra is a specific branch of mathematics dealing with the study of vectors, vector spaces with functions that
More informationOn the Equivalence of OKID and Time Series Identification for Markov-Parameter Estimation
On the Equivalence of OKID and Time Series Identification for Markov-Parameter Estimation P V Albuquerque, M Holzel, and D S Bernstein April 5, 2009 Abstract We show the equivalence of Observer/Kalman
More informationChapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are real numbers. 1
More information(c)
1. Find the reduced echelon form of the matrix 1 1 5 1 8 5. 1 1 1 (a) 3 1 3 0 1 3 1 (b) 0 0 1 (c) 3 0 0 1 0 (d) 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 (e) 1 0 5 0 0 1 3 0 0 0 0 Solution. 1 1 1 1 1 1 1 1
More informationSpan and Linear Independence
Span and Linear Independence It is common to confuse span and linear independence, because although they are different concepts, they are related. To see their relationship, let s revisit the previous
More informationMath 4A Notes. Written by Victoria Kala Last updated June 11, 2017
Math 4A Notes Written by Victoria Kala vtkala@math.ucsb.edu Last updated June 11, 2017 Systems of Linear Equations A linear equation is an equation that can be written in the form a 1 x 1 + a 2 x 2 +...
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 informationChapter 8 Linear Algebraic Equations
PowerPoint to accompany Introduction to MATLAB for Engineers, Third Edition William J. Palm III Chapter 8 Linear Algebraic Equations Copyright 2010. The McGraw-Hill Companies, Inc. This work is only for
More information1 9/5 Matrices, vectors, and their applications
1 9/5 Matrices, vectors, and their applications Algebra: study of objects and operations on them. Linear algebra: object: matrices and vectors. operations: addition, multiplication etc. Algorithms/Geometric
More informationLinear Algebra. Session 8
Linear Algebra. Session 8 Dr. Marco A Roque Sol 08/01/2017 Abstract Linear Algebra Range and kernel Let V, W be vector spaces and L : V W, be a linear mapping. Definition. The range (or image of L is the
More informationMATH 33A LECTURE 2 SOLUTIONS 1ST MIDTERM
MATH 33A LECTURE 2 SOLUTIONS ST MIDTERM MATH 33A LECTURE 2 SOLUTIONS ST MIDTERM 2 Problem. (True/False, pt each) Mark your answers by filling in the appropriate box next to each question. 2 3 7 (a T F
More informationMath 308 Discussion Problems #4 Chapter 4 (after 4.3)
Math 38 Discussion Problems #4 Chapter 4 (after 4.3) () (after 4.) Let S be a plane in R 3 passing through the origin, so that S is a two-dimensional subspace of R 3. Say that a linear transformation T
More informationSingular Value Decomposition
Singular Value Decomposition Motivatation The diagonalization theorem play a part in many interesting applications. Unfortunately not all matrices can be factored as A = PDP However a factorization A =
More informationObservers for Linear Systems with Unknown Inputs
Chapter 3 Observers for Linear Systems with Unknown Inputs As discussed in the previous chapters, it is often the case that a dynamic system can be modeled as having unknown inputs (e.g., representing
More informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
More informationCDS Solutions to Final Exam
CDS 22 - Solutions to Final Exam Instructor: Danielle C Tarraf Fall 27 Problem (a) We will compute the H 2 norm of G using state-space methods (see Section 26 in DFT) We begin by finding a minimal state-space
More informationHamilton-Jacobi-Bellman Equation Feb 25, 2008
Hamilton-Jacobi-Bellman Equation Feb 25, 2008 What is it? The Hamilton-Jacobi-Bellman (HJB) equation is the continuous-time analog to the discrete deterministic dynamic programming algorithm Discrete VS
More informationReachability and Controllability
Capitolo. INTRODUCTION 4. Reachability and Controllability Reachability. The reachability problem is to find the set of all the final states x(t ) reachable starting from a given initial state x(t ) :
More informationECE 275A Homework #3 Solutions
ECE 75A Homework #3 Solutions. Proof of (a). Obviously Ax = 0 y, Ax = 0 for all y. To show sufficiency, note that if y, Ax = 0 for all y, then it must certainly be true for the particular value of y =
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 informationChapter 3. Directions: For questions 1-11 mark each statement True or False. Justify each answer.
Chapter 3 Directions: For questions 1-11 mark each statement True or False. Justify each answer. 1. (True False) Asking whether the linear system corresponding to an augmented matrix [ a 1 a 2 a 3 b ]
More informationCS 143 Linear Algebra Review
CS 143 Linear Algebra Review Stefan Roth September 29, 2003 Introductory Remarks This review does not aim at mathematical rigor very much, but instead at ease of understanding and conciseness. Please see
More informationMATH 300, Second Exam REVIEW SOLUTIONS. NOTE: You may use a calculator for this exam- You only need something that will perform basic arithmetic.
MATH 300, Second Exam REVIEW SOLUTIONS NOTE: You may use a calculator for this exam- You only need something that will perform basic arithmetic. [ ] [ ] 2 2. Let u = and v =, Let S be the parallelegram
More informationChapter 5. Linear Algebra. Sections A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra Sections 5.1 5.3 A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are
More informationAN ITERATION. In part as motivation, we consider an iteration method for solving a system of linear equations which has the form x Ax = b
AN ITERATION In part as motivation, we consider an iteration method for solving a system of linear equations which has the form x Ax = b In this, A is an n n matrix and b R n.systemsof this form arise
More informationLinear Systems. Manfred Morari Melanie Zeilinger. Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich
Linear Systems Manfred Morari Melanie Zeilinger Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich Spring Semester 2016 Linear Systems M. Morari, M. Zeilinger - Spring
More informationMath 2233 Homework Set 7
Math 33 Homework Set 7 1. Find the general solution to the following differential equations. If initial conditions are specified, also determine the solution satisfying those initial conditions. a y 4
More informationMATH 152 Exam 1-Solutions 135 pts. Write your answers on separate paper. You do not need to copy the questions. Show your work!!!
MATH Exam -Solutions pts Write your answers on separate paper. You do not need to copy the questions. Show your work!!!. ( pts) Find the reduced row echelon form of the matrix Solution : 4 4 6 4 4 R R
More information