M403(2012) Solutions Problem Set 1 (c) 2012, Philip D. Loewen. = ( 1 λ) 2 k, k 1 λ
|
|
- Thomas Owen
- 6 years ago
- Views:
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 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 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 information+ 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 information20D - 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 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 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 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 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 informationMatrices 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 informationUnderstand 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 informationLinear Algebra Review (Course Notes for Math 308H - Spring 2016)
Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,
More informationLecture 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 informationMath : 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 informationPoincaré 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 informationOctober 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 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 informationSystems 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 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 informationMATH 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 informationEigenvalues 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 informationSeptember 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 informationCalculus 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 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 informationBefore 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 informationStudy 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 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 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 informationMath 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 informationLyapunov 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 information1. 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 informationSOLUTIONS: 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 informationFull-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 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 informationLinear 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 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 informationRecall : 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 informationThe 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 informationNonlinear 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 informationSystems 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 informationLinear 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 informationModeling 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 informationChapter 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 informationQUANTUM 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 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 informationSolutions 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 informationGrammians. 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 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 information6 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 informationPeriodic 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 informationLinear 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 informationMATH 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 informationIntroduction 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 informationSystems 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 informationMath 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 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 informationMath 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 informationFrom 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 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 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 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 informationHomework 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 informationLinearization 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 informationComplex 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 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 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 informationLecture 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 informationCheat 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
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 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 informationMATH 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 informationCalculus 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 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 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 informationMathQuest: 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 informationMATHEMATICS 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 informationMA 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 information18.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 informationSmall 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 informationCS 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 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 informationThe 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 informationMath 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 informationLMI 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 informationMath 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 informationSOLUTIONS 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 informationTheorem 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 informationEN221 - 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 informationMath 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 informationDiagonalization 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 informationTopic # /31 Feedback Control Systems
Topic #7 16.30/31 Feedback Control Systems State-Space Systems What are the basic properties of a state-space model, and how do we analyze these? Time Domain Interpretations System Modes Fall 2010 16.30/31
More informationPH.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 informationMATH 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 informationSTEP Support Programme. STEP 2 Matrices Topic Notes
STEP Support Programme STEP 2 Matrices Topic Notes Definitions............................................. 2 Manipulating Matrices...................................... 3 Transformations.........................................
More informationLINEAR 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 information3. 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 informationRobust Control 2 Controllability, Observability & Transfer Functions
Robust Control 2 Controllability, Observability & Transfer Functions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /26/24 Outline Reachable Controllability Distinguishable
More informationMatrix-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 informationLecture 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 informationEE16B - 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 informationLinear 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