M403(2012) Solutions Problem Set 1 (c) 2012, Philip D. Loewen. = ( 1 λ) 2 k, k 1 λ

Size: px
Start display at page:

Download "M403(2012) Solutions Problem Set 1 (c) 2012, Philip D. Loewen. = ( 1 λ) 2 k, k 1 λ"

Transcription

1 M43(22) Solutions Problem Set (c) 22, Philip D. Loewen. The number λ C is an eigenvalue for A exactly when λ = det(a λi) = det = ( λ) 2 k, k λ i.e., when λ = ± k. (i) When k = 2, the eigenvalues are ± 2, or approximately.77 and Corresponding eigenvectors u and v are given by λ = + 2 = 2 u u 2 = = u = ( 2, ), λ = 2 = The general solution of ( ) in this case is x(t) = C e ( +/ 2)t 2 2 v v 2 = = v = ( 2, ), ] + C 2 e ( / 2)t ] 2, C,C 2 R. () (ii) When k = 2, the eigenvalues are ± 2, or approximately.442 and Corresponding eigenvectors u and v are given by λ = + 2 = 2u u 2 = = u = (, 2), λ = 2 = 2v v 2 = = v = (, 2), The general solution of ( ) in this case is x(t) = C e ( + 2)t + C 2 2 e ( 2)t 2, C,C 2 R. (2) (iii) When k = 2, both eigenvalues of A are negative, so each trajectory in () converges to as t. By contrast, when k = 2, the first eigenvalue of A is positive. Hence each trajectory in (2) with C obeys x(t) as t. The general formula for the eigenvalues, λ = ± k, shows that transition between these situations happens when k = : for k < both eigenvalues have negative real part (they are complex when k < ), whereas whenever k >, one eigenvalue will be a positive real number. 2. (i) Here A = λi + M, where The ODE system ẋ = Mx, in components, says M =. ẋ =, ẋ 2 = x, ẋ 3 =. Therefore the general solution has x (t) = c and x 3 (t) = c 3 for some constants c and c 3, while x 2 (t) = c t + c 2 for some constant c 2. With this notation, x() = c, and it follows that x(t; ξ) = ξ + tξ ê 2. File hw, version of 7 September 22, page. Typeset at 9: September 7, 22.

2 2 UBC M43 Solutions for Problem Set # Thus we obtain e Mt = x(t;ê ) x(t;ê 2 ) x(t;ê 3 )] = t. Finally, since λi commutes with M, we can write e At = e (λi+m)t = e λit e Mt = e λt t = eλt te λt e λt. e λt (ii) Expanding the determinant along the top row, we find det(a λi) = det 3 λ 2 λ = (3 λ)(2 λ)( λ) + (2 λ) λ = (2 λ)(3 λ)( λ) + ] = (2 λ) 3. Thus A has a triple eigenvalue at 2. Since the matrix 2I commutes with any other matrix, the laws of matrix exponents give e At = e (A 2I)t+2It = e 2It e (A 2I)t. Writing M = A 2I, we have M =. Writing ẋ = Mx in components, we have ẋ = x + x 3, ẋ 2 =, ẋ 3 = x + x 2 x 3. Therefore x 2 (t) = c 2 (some constant) in general, while ẍ = ẋ + ẋ 3 = (x + x 3 ) + ( x + x 2 x 3 ) = x 2 (t) = c 2. Thus x (t) = 2 c 2t 2 + kt + c for some constants c and k, and then the first equation gives x 3 = ẋ x = (c 2 t + k) ( 2 c 2t 2 + kt + c ) = 2 c 2t 2 + (c 2 k)t + (k c ). Let s define c 3 = k c to eliminate k = c + c 3. This gives x() = c, and hence x(t; ξ) = ξ + (ξ + ξ 3 )t + 2 ξ 2t 2 ξ 2 = e Mt = x(t;ê ) x(t;ê 2 ) x(t;ê 3 )] = + t t2 /2 t. ξ 3 + (ξ 2 ξ ξ 3 )t 2 ξ 2t 2 t t 2 /2 t Consequently e At = e 2It e Mt = e 2t + t t2 /2 t. t t 2 /2 t Remark. In both parts, a power-series approach also works well: in part (i), one has M 2 = ; in part (ii), M 3 =. File hw, version of 7 September 22, page 2. Typeset at 9: September 7, 22.

3 UBC M43 Solutions for Problem Set # 3 3. In components, the ODE system ẋ = Ax becomes ẋ = x 2, ẋ 2 = ω2 R x 3, ẋ 3 = R x 2. The third equation implies ẍ 3 = Rẋ2, which the second equation reduces to ẍ 3 = R ( ) ω2 R x 3, i.e., ẍ 3 + ω2 R 2x 3 =. The general solution is x 3 (t) = αcos(ωt/r)+β sin(ωt/r) for arbitrary α,β R. Now the third equation gives x 2 (t) = Rẋ 3 (t) = αω sin(ωt/r) + βω cos(ωt/r), and the first equation gives some constant γ such that ẋ = x 2 = αω sin(ωt/r) + βω cos(ωt/r), so x (t) = αr cos(ωt/r) + βr sinωt/r) + γ. The general solution for x involves the three parameters α,β,γ: αr cos(ωt/r) + βr sin(ωt/r) + γ x = αω sin(ωt/r) + βω cos(ωt/r) α cos(ωt/r) + β sin(ωt/r) Note that x() = (αr + γ, ωβ, α). Therefore x() = e α =, β =, γ = x(t) = (,, ), x() = e 2 α =, β = /ω, γ = x(t) = ((R/ω)sin(ωt/R), cos(ωt/r), (/ω)sin(ωt/r)), x() = e 3 α =, β =, γ = R x(t) = (Rcos(ωt/R) ], ω sin(ωt/r), cos(ωt/r)), These provide the columns of the matrix e At : (R/ω)sin(ωt/R) R cos(ωt/r) R e At = cos(ωt/r) ω sin(ωt/r) (/ω) sin(ωt/r) cos(ωt/r) This reduces correctly to I when t =, and the comparison below confirms that the calculation is correct: (R/ω)sin(ωt/R) R cos(ωt/r) R cos(ωt/r) ω sin(ωt/r) d cos(ωt/r) ω sin(ωt/r) = (ω/r)sin(ωt/r) (ω 2 /R)cos(ωt/R), dt (/ω) sin(ωt/r) cos(ωt/r) (/R) cos(ωt/r) (ω/r) sin(ωt/r) (R/ω)sin(ωt/R) R cos(ωt/r) R cos(ωt/r) ω sin(ωt/r) ω 2 /R cos(ωt/r) ω sin(ωt/r) = (ω/r)sin(ωt/r) (ω 2 /R)cos(ωt/R). /R (/ω) sin(ωt/r) cos(ωt/r) (/R) cos(ωt/r) (ω/r) sin(ωt/r) 4. (a) The vector-matrix system ẋ = Ax expands to ẋ = x 2, ẋ 2 = x 2, giving ẍ = ẋ 2 = x 2 = ẋ, i.e., ( ) ẍ + ẋ =. Equation ( ) is solved by x = e st when = s 2 + s, i.e., when s = or s =, giving x (t) = α + βe, α,β R. File hw, version of 7 September 22, page 3. Typeset at 9: September 7, 22.

4 4 UBC M43 Solutions for Problem Set # Recalling x 2 = ẋ leads to the general solution of ẋ = Ax: e x(t) = α + β, α,β R. e This has x() = (α + β, β). Thus the initial condition x() = ê leads to α =, β =, and x(t;ê ) = (,). The initial condition x() = ê 2 leads to α =, β =, so x(t;ê 2 ) = ( e,e ). In summary, ] e At = x(t;ê ) e x(t;ê 2) = e. (b) Expanding the standard formula produces the component equations x(t) = e At ξ + e A(t s) Bu(s)ds e ξ e (t s) = e + ξ 2 e (t s) ξ + ξ = 2 ξ 2 e 2 e (t s) ξ 2 e + x (t) = ξ + ξ 2 ξ 2 e + x 2 (t) = ξ 2 e + e (t s) e (t s) u(s)ds. ] u(s) ds ] u(s) ds (2 e (t s) ) (c) When (ξ,ξ 2 ) = (,3), we get x() = exactly when u satisfies () = x () = 2 3e + (2) = x 2 () = 3e + (2 e ( s) ) e ( s) u(s)ds. An equivalent system is produced by adding the equations, and by multiplying e into (2): ( ) = (2 ) = 3 + e s u(s)ds. To choose a function u that solves this system of two equations, one might try to introduce just two unknowns by choosing u from some two-parameter family of functions. Endless variety exists. For example, one might try u(t) = mt + b for some constants m,b. Integration by parts leads to = 3 = m (ms + b)ds = m + b, i.e., m = 2(b + ), 2 se s ds + b e s ds = m (s )e s + b(e ) = m + b(e ). s= File hw, version of 7 September 22, page 4. Typeset at 9: September 7, 22.

5 UBC M43 Solutions for Problem Set # 5 Substituting the first equation into the second gives 3 = 2b 2 + (e )b, i.e., = b(e 3), i.e., b = 3 e ; ( ) 4 e Hence m = 2(b + ) = 2, and one possible control function is 3 e ( ) 4 e u(t) = 2 t + 3 e 3 e. Alternatively, one could use a piecewise-constant assumption like u(t) = { c, if t 2, c 2, if 2 < t. ( ) (Transition times other than 2 can also work.) Easy integrals lead to the system c + c 2 = 4, (e /2 )c + (e e /2 )c 2 = 3, so the constants that work with ( ) are c = e 4e ( e ) 2, c 2 = 4 e 7 ( e ) 2. File hw, version of 7 September 22, page 5. Typeset at 9: September 7, 22.

Section 8.2 : Homogeneous Linear Systems

Section 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 information

Matrix Solutions to Linear Systems of ODEs

Matrix 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 information

+ i. cos(t) + 2 sin(t) + c 2.

+ i. cos(t) + 2 sin(t) + c 2. MATH HOMEWORK #7 PART A SOLUTIONS Problem 7.6.. Consider the system x = 5 x. a Express the general solution of the given system of equations in terms of realvalued functions. b Draw a direction field,

More information

20D - Homework Assignment 5

20D - Homework Assignment 5 Brian Bowers TA for Hui Sun MATH D Homework Assignment 5 November 8, 3 D - Homework Assignment 5 First, I present the list of all matrix row operations. We use combinations of these steps to row reduce

More information

Modal Decomposition and the Time-Domain Response of Linear Systems 1

Modal 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 information

21 Linear State-Space Representations

21 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 information

Differential equations

Differential 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 information

Math Ordinary Differential Equations

Math 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 information

Matrices and Linear Algebra

Matrices and Linear Algebra Contents Quantitative methods for Economics and Business University of Ferrara Academic year 2017-2018 Contents 1 Basics 2 3 4 5 Contents 1 Basics 2 3 4 5 Contents 1 Basics 2 3 4 5 Contents 1 Basics 2

More information

Understand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.

Understand the existence and uniqueness theorems and what they tell you about solutions to initial value problems. Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics

More information

Linear Algebra Review (Course Notes for Math 308H - Spring 2016)

Linear 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 information

Lecture Notes for Math 251: ODE and PDE. Lecture 27: 7.8 Repeated Eigenvalues

Lecture Notes for Math 251: ODE and PDE. Lecture 27: 7.8 Repeated Eigenvalues Lecture Notes for Math 25: ODE and PDE. Lecture 27: 7.8 Repeated Eigenvalues Shawn D. Ryan Spring 22 Repeated Eigenvalues Last Time: We studied phase portraits and systems of differential equations with

More information

Math : Solutions to Assignment 10

Math : Solutions to Assignment 10 Math -3: Solutions to Assignment. There are two tanks. The first tank initially has gallons of pure water. The second tank initially has 8 gallons of a water/salt solution with oz of salt. Both tanks drain

More information

Poincaré Map, Floquet Theory, and Stability of Periodic Orbits

Poincaré Map, Floquet Theory, and Stability of Periodic Orbits Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct

More information

October 4, 2017 EIGENVALUES AND EIGENVECTORS. APPLICATIONS

October 4, 2017 EIGENVALUES AND EIGENVECTORS. APPLICATIONS October 4, 207 EIGENVALUES AND EIGENVECTORS. APPLICATIONS RODICA D. COSTIN Contents 4. Eigenvalues and Eigenvectors 3 4.. Motivation 3 4.2. Diagonal matrices 3 4.3. Example: solving linear differential

More information

State will have dimension 5. One possible choice is given by y and its derivatives up to y (4)

State 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 information

Systems of differential equations Handout

Systems of differential equations Handout Systems of differential equations Handout Peyam Tabrizian Friday, November 8th, This handout is meant to give you a couple more example of all the techniques discussed in chapter 9, to counterbalance all

More information

and let s calculate the image of some vectors under the transformation T.

and 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 information

MATH 4B Differential Equations, Fall 2016 Final Exam Study Guide

MATH 4B Differential Equations, Fall 2016 Final Exam Study Guide MATH 4B Differential Equations, Fall 2016 Final Exam Study Guide GENERAL INFORMATION AND FINAL EXAM RULES The exam will have a duration of 3 hours. No extra time will be given. Failing to submit your solutions

More information

Eigenvalues and Eigenvectors, More Direction Fields and Systems of ODEs

Eigenvalues and Eigenvectors, More Direction Fields and Systems of ODEs Eigenvalues and Eigenvectors, More Direction Fields and Systems of ODEs First let us speak a bit about eigenvalues. Defn. An eigenvalue λ of an nxn matrix A means a scalar (perhaps a complex number) such

More information

September 26, 2017 EIGENVALUES AND EIGENVECTORS. APPLICATIONS

September 26, 2017 EIGENVALUES AND EIGENVECTORS. APPLICATIONS September 26, 207 EIGENVALUES AND EIGENVECTORS. APPLICATIONS RODICA D. COSTIN Contents 4. Eigenvalues and Eigenvectors 3 4.. Motivation 3 4.2. Diagonal matrices 3 4.3. Example: solving linear differential

More information

Calculus and Differential Equations II

Calculus and Differential Equations II MATH 250 B Second order autonomous linear systems We are mostly interested with 2 2 first order autonomous systems of the form { x = a x + b y y = c x + d y where x and y are functions of t and a, b, c,

More information

CONTROL DESIGN FOR SET POINT TRACKING

CONTROL 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 information

Before you begin read these instructions carefully.

Before you begin read these instructions carefully. MATHEMATICAL TRIPOS Part IA Friday, 29 May, 2015 1:30 pm to 4:30 pm PAPER 2 Before you begin read these instructions carefully. The examination paper is divided into two sections. Each question in Section

More information

Study guide - Math 220

Study guide - Math 220 Study guide - Math 220 November 28, 2012 1 Exam I 1.1 Linear Equations An equation is linear, if in the form y + p(t)y = q(t). Introducing the integrating factor µ(t) = e p(t)dt the solutions is then in

More information

Homogeneous Linear Systems of Differential Equations with Constant Coefficients

Homogeneous 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 information

APPPHYS 217 Tuesday 6 April 2010

APPPHYS 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 information

Math 250B Final Exam Review Session Spring 2015 SOLUTIONS

Math 250B Final Exam Review Session Spring 2015 SOLUTIONS Math 5B Final Exam Review Session Spring 5 SOLUTIONS Problem Solve x x + y + 54te 3t and y x + 4y + 9e 3t λ SOLUTION: We have det(a λi) if and only if if and 4 λ only if λ 3λ This means that the eigenvalues

More information

Lyapunov Stability Theory

Lyapunov Stability Theory Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous

More information

1. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal

1. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal . Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal 3 9 matrix D such that A = P DP, for A =. 3 4 3 (a) P = 4, D =. 3 (b) P = 4, D =. (c) P = 4 8 4, D =. 3 (d) P

More information

SOLUTIONS: ASSIGNMENT Use Gaussian elimination to find the determinant of the matrix. = det. = det = 1 ( 2) 3 6 = 36. v 4.

SOLUTIONS: ASSIGNMENT Use Gaussian elimination to find the determinant of the matrix. = det. = det = 1 ( 2) 3 6 = 36. v 4. SOLUTIONS: ASSIGNMENT 9 66 Use Gaussian elimination to find the determinant of the matrix det 1 1 4 4 1 1 1 1 8 8 = det = det 0 7 9 0 0 0 6 = 1 ( ) 3 6 = 36 = det = det 0 0 6 1 0 0 0 6 61 Consider a 4

More information

Full-State Feedback Design for a Multi-Input System

Full-State Feedback Design for a Multi-Input System Full-State Feedback Design for a Multi-Input System A. Introduction The open-loop system is described by the following state space model. x(t) = Ax(t)+Bu(t), y(t) =Cx(t)+Du(t) () 4 8.5 A =, B =.5.5, C

More information

Solution via Laplace transform and matrix exponential

Solution 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 information

Linear Algebra Exercises

Linear Algebra Exercises 9. 8.03 Linear Algebra Exercises 9A. Matrix Multiplication, Rank, Echelon Form 9A-. Which of the following matrices is in row-echelon form? 2 0 0 5 0 (i) (ii) (iii) (iv) 0 0 0 (v) [ 0 ] 0 0 0 0 0 0 0 9A-2.

More information

Chapter 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. 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

Recall : Eigenvalues and Eigenvectors

Recall : Eigenvalues and Eigenvectors Recall : Eigenvalues and Eigenvectors Let A be an n n matrix. If a nonzero vector x in R n satisfies Ax λx for a scalar λ, then : The scalar λ is called an eigenvalue of A. The vector x is called an eigenvector

More information

The Rocket Car. UBC Math 403 Lecture Notes by Philip D. Loewen

The Rocket Car. UBC Math 403 Lecture Notes by Philip D. Loewen The Rocket Car UBC Math 403 Lecture Notes by Philip D. Loewen We consider this simple system with control constraints: ẋ 1 = x, ẋ = u, u [ 1, 1], Dynamics. Consider the system evolution on a time interval

More information

Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1

Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1 Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium

More information

Systems of Algebraic Equations and Systems of Differential Equations

Systems of Algebraic Equations and Systems of Differential Equations Systems of Algebraic Equations and Systems of Differential Equations Topics: 2 by 2 systems of linear equations Matrix expression; Ax = b Solving 2 by 2 homogeneous systems Functions defined on matrices

More information

Linear Algebra Practice Problems

Linear Algebra Practice Problems Linear Algebra Practice Problems Math 24 Calculus III Summer 25, Session II. Determine whether the given set is a vector space. If not, give at least one axiom that is not satisfied. Unless otherwise stated,

More information

Modeling and Simulation with ODE for MSE

Modeling and Simulation with ODE for MSE Zentrum Mathematik Technische Universität München Prof. Dr. Massimo Fornasier WS 6/7 Dr. Markus Hansen Sheet 7 Modeling and Simulation with ODE for MSE The exercises can be handed in until Wed, 4..6,.

More information

Chapter 2 Notes, Linear Algebra 5e Lay

Chapter 2 Notes, Linear Algebra 5e Lay Contents.1 Operations with Matrices..................................1.1 Addition and Subtraction.............................1. Multiplication by a scalar............................ 3.1.3 Multiplication

More information

QUANTUM MECHANICS I PHYS 516. Solutions to Problem Set # 5

QUANTUM MECHANICS I PHYS 516. Solutions to Problem Set # 5 QUANTUM MECHANICS I PHYS 56 Solutions to Problem Set # 5. Crossed E and B fields: A hydrogen atom in the N 2 level is subject to crossed electric and magnetic fields. Choose your coordinate axes to make

More information

Module 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 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 information

Solutions Problem Set 8 Math 240, Fall

Solutions Problem Set 8 Math 240, Fall Solutions Problem Set 8 Math 240, Fall 2012 5.6 T/F.2. True. If A is upper or lower diagonal, to make det(a λi) 0, we need product of the main diagonal elements of A λi to be 0, which means λ is one of

More information

Grammians. Matthew M. Peet. Lecture 20: Grammians. Illinois Institute of Technology

Grammians. Matthew M. Peet. Lecture 20: Grammians. Illinois Institute of Technology Grammians Matthew M. Peet Illinois Institute of Technology Lecture 2: Grammians Lyapunov Equations Proposition 1. Suppose A is Hurwitz and Q is a square matrix. Then X = e AT s Qe As ds is the unique solution

More information

Eigenpairs and Diagonalizability Math 401, Spring 2010, Professor David Levermore

Eigenpairs 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 information

6 Linear Equation. 6.1 Equation with constant coefficients

6 Linear Equation. 6.1 Equation with constant coefficients 6 Linear Equation 6.1 Equation with constant coefficients Consider the equation ẋ = Ax, x R n. This equating has n independent solutions. If the eigenvalues are distinct then the solutions are c k e λ

More information

Periodic Linear Systems

Periodic Linear Systems Periodic Linear Systems Lecture 17 Math 634 10/8/99 We now consider ẋ = A(t)x (1) when A is a continuous periodic n n matrix function of t; i.e., when there is a constant T>0 such that A(t + T )=A(t) for

More information

Linear Differential Equations. Problems

Linear Differential Equations. Problems Chapter 1 Linear Differential Equations. Problems 1.1 Introduction 1.1.1 Show that the function ϕ : R R, given by the expression ϕ(t) = 2e 3t for all t R, is a solution of the Initial Value Problem x =

More information

MATH 5640: Functions of Diagonalizable Matrices

MATH 5640: Functions of Diagonalizable Matrices MATH 5640: Functions of Diagonalizable Matrices Hung Phan, UMass Lowell November 27, 208 Spectral theorem for diagonalizable matrices Definition Let V = X Y Every v V is uniquely decomposed as u = x +

More information

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t

More information

Systems of Linear ODEs

Systems of Linear ODEs P a g e 1 Systems of Linear ODEs Systems of ordinary differential equations can be solved in much the same way as discrete dynamical systems if the differential equations are linear. We will focus here

More information

Math 323 Exam 2 Sample Problems Solution Guide October 31, 2013

Math 323 Exam 2 Sample Problems Solution Guide October 31, 2013 Math Exam Sample Problems Solution Guide October, Note that the following provides a guide to the solutions on the sample problems, but in some cases the complete solution would require more work or justification

More information

Solution: In standard form (i.e. y + P (t)y = Q(t)) we have y t y = cos(t)

Solution: 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 information

Math 331 Homework Assignment Chapter 7 Page 1 of 9

Math 331 Homework Assignment Chapter 7 Page 1 of 9 Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a

More information

From Lay, 5.4. If we always treat a matrix as defining a linear transformation, what role does diagonalisation play?

From Lay, 5.4. If we always treat a matrix as defining a linear transformation, what role does diagonalisation play? Overview Last week introduced the important Diagonalisation Theorem: An n n matrix A is diagonalisable if and only if there is a basis for R n consisting of eigenvectors of A. This week we ll continue

More information

CAAM 335 Matrix Analysis

CAAM 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 information

ECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77

ECEN 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 information

LS.2 Homogeneous Linear Systems with Constant Coefficients

LS.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 information

Homework 3 Solutions Math 309, Fall 2015

Homework 3 Solutions Math 309, Fall 2015 Homework 3 Solutions Math 39, Fall 25 782 One easily checks that the only eigenvalue of the coefficient matrix is λ To find the associated eigenvector, we have 4 2 v v 8 4 (up to scalar multiplication)

More information

Linearization of Differential Equation Models

Linearization of Differential Equation Models Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking

More information

Complex Dynamic Systems: Qualitative vs Quantitative analysis

Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic

More information

Dynamic interpretation of eigenvectors

Dynamic 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 information

Section 9.8 Higher Order Linear Equations

Section 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 information

Lecture 6. Eigen-analysis

Lecture 6. Eigen-analysis Lecture 6 Eigen-analysis University of British Columbia, Vancouver Yue-Xian Li March 7 6 Definition of eigenvectors and eigenvalues Def: Any n n matrix A defines a LT, A : R n R n A vector v and a scalar

More information

Cheat Sheet for MATH461

Cheat Sheet for MATH461 Cheat Sheet for MATH46 Here is the stuff you really need to remember for the exams Linear systems Ax = b Problem: We consider a linear system of m equations for n unknowns x,,x n : For a given matrix A

More information

π 1 = tr(a), π n = ( 1) n det(a). In particular, when n = 2 one has

π 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 information

A 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

A 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 information

MATH 320 INHOMOGENEOUS LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS

MATH 320 INHOMOGENEOUS LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS MATH 2 INHOMOGENEOUS LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS W To find a particular solution for a linear inhomogeneous system of differential equations x Ax = ft) or of a mechanical system with external

More information

Calculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y

Calculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y Calculus for the Life Sciences II Assignment 6 solutions Find the tangent plane to the graph of the function at the point (0, π f(x, y = 3π 3 cos 2x + 2 sin 3y Solution: The tangent plane of f at a point

More information

DIAGONALIZATION. In order to see the implications of this definition, let us consider the following example Example 1. Consider the matrix

DIAGONALIZATION. 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 information

A 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

A 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 information

MathQuest: Differential Equations

MathQuest: Differential Equations MathQuest: Differential Equations Solutions to Linear Systems. Consider the linear system given by dy dt = 4 True or False: Y e t t = is a solution. c False, but I am not very confident Y.. Consider the

More information

MATHEMATICS FOR ENGINEERS BASIC MATRIX THEORY TUTORIAL 3 - EIGENVECTORS AND EIGENVALUES

MATHEMATICS FOR ENGINEERS BASIC MATRIX THEORY TUTORIAL 3 - EIGENVECTORS AND EIGENVALUES MATHEMATICS FOR ENGINEERS BASIC MATRIX THEORY TUTORIAL 3 - EIGENVECTORS AND EIGENVALUES This is the third tutorial on matrix theory. It is entirely devoted to the subject of Eigenvectors and Eigenvalues

More information

MA Ordinary Differential Equations and Control Semester 1 Lecture Notes on ODEs and Laplace Transform (Parts I and II)

MA Ordinary Differential Equations and Control Semester 1 Lecture Notes on ODEs and Laplace Transform (Parts I and II) MA222 - Ordinary Differential Equations and Control Semester Lecture Notes on ODEs and Laplace Transform (Parts I and II Dr J D Evans October 22 Special thanks (in alphabetical order are due to Kenneth

More information

18.06 Professor Strang Quiz 3 May 2, 2008

18.06 Professor Strang Quiz 3 May 2, 2008 18.06 Professor Strang Quiz 3 May 2, 2008 Your PRINTED name is: Grading 1 2 3 Please circle your recitation: 1) M 2 2-131 A. Ritter 2-085 2-1192 afr 2) M 2 4-149 A. Tievsky 2-492 3-4093 tievsky 3) M 3

More information

Small oscillations and normal modes

Small oscillations and normal modes Chapter 4 Small oscillations and normal modes 4.1 Linear oscillations Discuss a generalization of the harmonic oscillator problem: oscillations of a system of several degrees of freedom near the position

More information

CS 246 Review of Linear Algebra 01/17/19

CS 246 Review of Linear Algebra 01/17/19 1 Linear algebra In this section we will discuss vectors and matrices. We denote the (i, j)th entry of a matrix A as A ij, and the ith entry of a vector as v i. 1.1 Vectors and vector operations A vector

More information

EE263: Introduction to Linear Dynamical Systems Review Session 5

EE263: 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 information

The Contraction Mapping Theorem and the Implicit and Inverse Function Theorems

The Contraction Mapping Theorem and the Implicit and Inverse Function Theorems The Contraction Mapping Theorem and the Implicit and Inverse Function Theorems The Contraction Mapping Theorem Theorem The Contraction Mapping Theorem) Let B a = { x IR d x < a } denote the open ball of

More information

Math 215 HW #9 Solutions

Math 215 HW #9 Solutions Math 5 HW #9 Solutions. Problem 4.4.. If A is a 5 by 5 matrix with all a ij, then det A. Volumes or the big formula or pivots should give some upper bound on the determinant. Answer: Let v i be the ith

More information

LMI Methods in Optimal and Robust Control

LMI Methods in Optimal and Robust Control LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 4: LMIs for State-Space Internal Stability Solving the Equations Find the output given the input State-Space:

More information

Math Bootcamp An p-dimensional vector is p numbers put together. Written as. x 1 x =. x p

Math Bootcamp An p-dimensional vector is p numbers put together. Written as. x 1 x =. x p Math Bootcamp 2012 1 Review of matrix algebra 1.1 Vectors and rules of operations An p-dimensional vector is p numbers put together. Written as x 1 x =. x p. When p = 1, this represents a point in the

More information

SOLUTIONS TO PRACTICE EXAM 3, SPRING 2004

SOLUTIONS TO PRACTICE EXAM 3, SPRING 2004 8034 SOLUTIONS TO PRACTICE EXAM 3, SPRING 004 Problem Let A be a real matrix and consider the linear system of first order differential equations, y y (t) Ay(t), y(t) (t) y (t) Let α be a real number,

More information

Theorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0).

Theorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0). Linear Systems Notes Lecture Proposition. A M n (R) is positive definite iff all nested minors are greater than or equal to zero. n Proof. ( ): Positive definite iff λ i >. Let det(a) = λj and H = {x D

More information

EN221 - Fall HW # 2 Solutions

EN221 - Fall HW # 2 Solutions EN1 - Fall8 - HW # Solutions Prof. Vivek Shenoy 1) Let A and B be arbitrary tensors, A and B the adjugates of A and B, and α and β arbitrary scalars. Show that det(αa + βb) = deta + β tr (B T A ) + αβ

More information

Math 250B Midterm III Information Fall 2018 SOLUTIONS TO PRACTICE PROBLEMS

Math 250B Midterm III Information Fall 2018 SOLUTIONS TO PRACTICE PROBLEMS Math 25B Midterm III Information Fall 28 SOLUTIONS TO PRACTICE PROBLEMS Problem Determine whether the following matrix is diagonalizable or not If it is, find an invertible matrix S and a diagonal matrix

More information

Diagonalization of Matrices

Diagonalization of Matrices LECTURE 4 Diagonalization of Matrices Recall that a diagonal matrix is a square n n matrix with non-zero entries only along the diagonal from the upper left to the lower right (the main diagonal) Diagonal

More information

Topic # /31 Feedback Control Systems

Topic # /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

PH.D. PRELIMINARY EXAMINATION MATHEMATICS

PH.D. PRELIMINARY EXAMINATION MATHEMATICS UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem

More information

MATH 315 Linear Algebra Homework #1 Assigned: August 20, 2018

MATH 315 Linear Algebra Homework #1 Assigned: August 20, 2018 Homework #1 Assigned: August 20, 2018 Review the following subjects involving systems of equations and matrices from Calculus II. Linear systems of equations Converting systems to matrix form Pivot entry

More information

STEP Support Programme. STEP 2 Matrices Topic Notes

STEP Support Programme. STEP 2 Matrices Topic Notes STEP Support Programme STEP 2 Matrices Topic Notes Definitions............................................. 2 Manipulating Matrices...................................... 3 Transformations.........................................

More information

LINEAR EQUATIONS OF HIGHER ORDER. EXAMPLES. General framework

LINEAR EQUATIONS OF HIGHER ORDER. EXAMPLES. General framework Differential Equations Grinshpan LINEAR EQUATIONS OF HIGHER ORDER. EXAMPLES. We consider linear ODE of order n: General framework (1) x (n) (t) + P n 1 (t)x (n 1) (t) + + P 1 (t)x (t) + P 0 (t)x(t) = 0

More information

3. Identify and find the general solution of each of the following first order differential equations.

3. Identify and find the general solution of each of the following first order differential equations. Final Exam MATH 33, Sample Questions. Fall 6. y = Cx 3 3 is the general solution of a differential equation. Find the equation. Answer: y = 3y + 9 xy. y = C x + C is the general solution of a differential

More information

Robust Control 2 Controllability, Observability & Transfer Functions

Robust 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 information

Matrix-Exponentials. September 7, dx dt = ax. x(t) = e at x(0)

Matrix-Exponentials. September 7, dx dt = ax. x(t) = e at x(0) Matrix-Exponentials September 7, 207 In [4]: using PyPlot INFO: Recompiling stale cache file /Users/stevenj/.julia/lib/v0.5/LaTeXStrings.ji for module LaTeXString Review: Solving ODEs via eigenvectors

More information

Lecture Notes: Eigenvalues and Eigenvectors. 1 Definitions. 2 Finding All Eigenvalues

Lecture Notes: Eigenvalues and Eigenvectors. 1 Definitions. 2 Finding All Eigenvalues Lecture Notes: Eigenvalues and Eigenvectors Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk 1 Definitions Let A be an n n matrix. If there

More information

EE16B - Spring 17 - Lecture 12A Notes 1

EE16B - Spring 17 - Lecture 12A Notes 1 EE6B - Spring 7 - Lecture 2A Notes Murat Arcak April 27 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4. International License. Sampling and Discrete Time Signals Discrete-Time

More information

Linear maps. Matthew Macauley. Department of Mathematical Sciences Clemson University Math 8530, Spring 2017

Linear maps. Matthew Macauley. Department of Mathematical Sciences Clemson University  Math 8530, Spring 2017 Linear maps Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 8530, Spring 2017 M. Macauley (Clemson) Linear maps Math 8530, Spring 2017

More information