Systems of Linear Differential Equations Chapter 7
|
|
- Darren Tucker
- 5 years ago
- Views:
Transcription
1 Systems of Linear Differential Equations Chapter 7 Doreen De Leon Department of Mathematics, California State University, Fresno June 22, 25 Motivating Examples: Applications of Systems of First Order Differential Equations Taken from an example in A First Course in Differential Equations with Modeling Applications, 8th ed, by Dennis G Zill Radioactive Series In our earlier discussion of radioactive decay, we assumed that the rate of decay was proportional to the number y(t of nuclei of the substance present at time t When a substance decays by radioactivity, it usually does not just transmute into one stable substance and then the proces stops Rather, the first substance decays into another radioactive substance, this substance in turn decays into a third radioactive substance, and so on This process, called a radioactive decay series, continues until a stable element is reached For example, the uranium decay series is U-238 Th-234 Pb-26, where Pb-26 is a stable isotope of lead The half-lives of the various elements in a radioactive series can range from billions of years (eg, 45 9 years for U-238 to a fraction of a second Suppose a radioactive series is described schematically by X λ Y λ 2 Z, where k λ < and k 2 λ 2 < are the decay constants for substances X and Y, respectively, and Z is a stable element Suppose that x(t, y(t, and z(t denote the amounts of substances X, Y, and Z, respectively, remaining at time t The decay of element X is described by dx dt λ x However, since X decays into Y (so Y is gaining atoms and Y is, itself, decaying (thereby losing atoms, the rate at which element Y decays is the net rate dy dt λ x λ 2 y
2 Since Z is a stable element, it is simply gaining atoms from the decay of element Y, so dz dt λ 2y In other words, a model of the radioactive decay series for the three elements is the linear system of first-order differential equations, 2 Electrical Networks dx dt λ x dy dt λ x λ 2 y dz dt λ 2y An electrical network having more than one loop also gives rise to a system of differential equations For example, in the figure below the current i (t splits in the directions shown at point B, called a branch point of the network By Kirchhoff s first law, we have i (t i 2 (t + i 3 (t ( Applying Kirchhoff s second law to each loop gives the following For loop A B B 2 A 2 A, summing the voltage drops across each part of the loop gives Similarly, for loop A B C C 2 B 2 A 2 A, we see that E(t i R + L di 2 dt + i 2R 2 (2 E(t i R + L 2 di 3 dt (3 2
3 We can use Equation ( to eliminate i from Equations (2 and (3 to obtain the following system of equations: L di 2 dt + (R + R 2 i 2 + R 3 i 3 E(t (4 L 2 di 3 dt + R i 2 + R i 3 E(t (5 2 First Order Systems of Differential Equations Section 7 We can write a higher order linear differential equation as a system of first order equations Note: The number of equations in the system equals the order of the differential equation Example: Write as a system of first order equations: x + p(tx + q(tx f(t So, the solution is Let x x x x x 2 x 2 x x 2 x q(tx p(tx f(t x 2 q(tx p(tx 2 f(t x x 2 x 2 q(tx p(tx 2 + f(t Our goal in this chapter will be to solve systems of first order constant coefficient differential equations using eigenvalue-eigenvector theory First, we must write the systems of differential equations in matrix form 3 Matrices and Systems of Linear Differential Equations Section 72 We can write a system of first order differential equations in matrix form Examples: ( Write the system below in matrix form x 2x 3x 2 x 2 4x + 2x 2 3
4 (2 Write the system below in matrix form ( x 2 3 x x x 2 dx ( 2 3 dt x 4 2 ( x x 2 x x 2 x 2 q(tx p(tx 2 + f(t x + q(t p(t x 2 ( f(t Consider the homogeneous differential equation dx dt P (tx(t (6 Theorem (Superposition Principle Given x (t, x 2 (t,, x n (t n solutions to (6 on the open interval I, if c, c 2,, c n are constants, then is also a solution of (6 on I Proof Since x (t, x 2 (t,, x n (t solve (6, So, x(t c x (t + c 2 x 2 (t + + c n x n (t x P (tx, x 2 P (tx 2,, x n P (tx n x c x + c 2 x c n x n c P (tx + c 2 P (tx c n P (tx n P (t (c x + c 2 x c n x n P (tx We can check the linear independence of solutions using the Wronskian If x, x 2,, x n are solutions of (6, their Wronskian is the determinant of the matrix whose columns are the solution vectors Theorem 2 Suppose x, x 2,, x n are solutions of the n n system (6 on an open interval I Suppose P (t is continuous on I Let W W (t Then x, x 2,, x n are linearly independent if and only if W for all t I They are linearly dependent if W on I 4
5 The dimension of the solution space of (6 equals the number of variables If x, x 2,, x n are n linearly independent solutions of the n n system (6, then the general solution is x(t c x (t + c 2 x 2 (t + + c n x n (t Initial Value Problems Example: Given the initial value problem 3 2 x x, x( : e 3t 2e 2t a Verify that x and x 2 are solutions 3e 3t e 2t b Show that x, x 2 are linearly independent c Solve the initial value problem Solution: a 3e x 3t 3 2 e 3t 3e 3t 9e 3t 3 4 3e 3t 9e 3t x 4e x 2t 3 2 2e 2t 4e 2t 2 2e 2t 3 4 e 2t x 2 2e 2t b W (t c e3t 3e 3t 2e 2t e 2t et 6e t 5e t, so x and x 2 are linearly independent e 3t 2e 2t x c x + c 2 x 2 c 3e 3t + c 2 e 2t ( 2 x( c 5 + c 3 2 c + 2c 2 3c + c 2 5 c 2 c 2 So, ( e 3t x 2 3e 3t ( ( 2e 2t 2e 3t 2e 2t 6e 3t e 2t e 2t 5
6 4 The Eigenvalue Method for Systems of First Order Linear Differential Equations Section 73 Goal: Solve x Ax, < t < (A a constant matrix Idea: Search for solutions of the form x e λt v Why? Consider the separable equation y λy The solution to this equation is y ce λt Plugging x into the system gives λe λt v Ae λt v λv Av Therefore, λ and v are the eigenvalues and corresponding eigenvectors of the matrix A Possibilities: Distinct eigenvalues Section 73 2 Complex conjugate eigenvalues Section 73 3 Repeated eigenvalues Section 75 4 Distinct Eigenvalues Procedure: Find the eigenvalues of A Solve det(a λi 2 For each eigenvalue, find corresponding eigenvector(s Find a basis for the solution space of (A λiv 3 Linear combinations yield general solutions Examples: ( Find a general solution of ( 4 x x 3 8 det(a λi λ λ (λ 2(λ 6 λ 2, λ 2 6 λ(8 λ + 2 λ 2 8λ + 2 6
7 λ 2 Find a basis for the solution space of (A 2Iv In matrix form, we have 2 4 v 3 6 v r r So, we have 2v + 4v 2 v 2v 2 ( r 2 r r ( 2 4 v 2 is the only arbitrary variable, so we may choose a value for it Let v 2 v 2, and ( 2 v Then, one solution is ( 2 x e λt v e 2t λ 2 6 Find a basis for the solution space of (A 6Iu In matrix form, we have ( 6 4 u 3 2 u 2 ( 6 4 r 2 r 3 2 r 2 r 2 +2r So, we have 3u + 2u 2 u 2 3 u 2 u 2 is the only arbitrary variable, so let u 2 3 u 2, and ( 2 u 3 Then, a second solution is The general solution is x c x + c 2 x 2 ( 2 x 2 e λ2t u e 6t Therefore, x c e 2t + c 2 e 6t 3 7
8 (2 Solve the initial value problem x x, x( 6 General solution: 3 λ 2 det(a λi λ 6 λ (3 λ( λ(6 λ λ,λ 2 3, λ 3 6 λ Find a basis for the solution space of (A Iv In matrix form, we have 2 2 v v 2 5 v r 3 r 3 5r 2 5 So we have v 3 2v + 2v 2 + v 3 v v 2 v 2 arbitrary, so let v 2 v and v Then, one solution is x e t λ 2 3 Find a basis for the solution space of (A 3Iu In matrix form, we have 2 u 2 u 2 3 u 3 8
9 2 2 3 So we have r 3 r r r 2 r 2 +r u 3 2u 2 + u 3 u 2 u is arbitrary, so let u u Then, a second solution is λ 3 6 Solve (A 6Iw x 2 e 3t In matrix form, solve 3 2 w 5 w 2 w 3 The coefficient matrix is already in echelon form, so we have 5w 2 + w 3 w 2 5 w 3 3w + 2w 2 + w 3 w 2 3 w w 3 w 3 is arbitrary, so since w 3 will be divided by both 3 and 5, let w 3 5 (eliminates fractions 7 w 2 3 and w 2(3 + (5 7 Therefore, w Then, a third solution is 7 x 3 e 6t 3 5 General solution: x c x + c 2 x 2 + c 3 x 3 7 x c e t + c 2 e 3t + c 3 e 6t 3 5 9
10 Find c, c 2, c 3 : So, Therefore, or 2 7 c + c 2 + c c + c 2 + 7c 3 2 c 2 3 c + 3c 3 c 5c 3 c 3 x e t + 3e 3t, e t + 3e 3t x e t 42 Complex Eigenvalues Goal: Obtain real-valued solutions Suppose λ a + bi, λ a bi are two eigenvalues corresponding to λ, then If v u + iw is an eigenvector (A λiv (A λiv (A λiv In other words, v is an eigenvector corresponding to λ The solution associated with λ and v is x ve λt (u + iwe (a+bit (u + iwe at (cos(bt + i sin(bt e at (u cos(bt w sin(bt +i e at (w cos(bt + u sin(bt }{{}}{{} x x 2
11 So, x and x 2 are also solutions to the system Therefore, we have two linearly independent real-valued solutions, x e at (u cos(bt w sin(bt and x 2 e at (w cos(bt + u sin(bt We get the same solutions by looking at ve λt So, we need only find one eigenvector corresponding to one of the complex conjugate eigenvalue pair in order to obtain a general solution Example: Solve ( 4 x x 5 det(a λi λ 4 5 λ ( λ( 5 λ + 4 λ 2 + 6λ + 45 λ 6 ± 36 4( ± 2i 2 3 ± 6i λ 3 + 6i Find a basis for the solution space of (A λiv In matrix form, we have ( ( 3 + 6i 4 ( v 5 ( 3 + 6i v 2 ( 2 6i 4 v 2 6i We can show that the augmented matrix can be reduced to ( 2 6i 4 Therefore, we have v 2 (2 6iv + 4v 2 ( v 2 is the only arbitrary variable, so let v 2 2 6i v 4 Thus, ( 4 v 2 6i
12 ( ( 4 4 e λt v e ( 3+6it e 3t (cos(6t + i sin(6t 2 6i 2 6i ( e 3t 4 cos(6t + i( 4 sin(6t 2 cos(6t + 6 sin(6t + i(2 sin(6t 6 cos(6t ( ( e 3t 4 cos(6t +i e 3t 4 sin(6t General solution: x c x + c 2 x 2, or 2 cos(6t + 6 sin(6t }{{} x 2 sin(6t 6 cos(6t }{{} x 2 ( ( x c e 3t 4 cos(6t + c 2 cos(6t + 6 sin(6t 2 e 3t 4 sin(6t 2 sin(6y 6 cos(6t 43 Repeated Eigenvalues Section 75 Definitions: Algebraic multiplicity: the order of the root Example:(λ 3 λ has algebraic multiplicity 3 Geometric multiplicity: the dimension of the eigenspace associated with an eigenvalue Case : algebraic multiplicity geometric multiplicty In this case, we obtain an eigenvector for each instance of the eigenvalue same solution strategy as for distinct eigenvalues Example: Find a general solution of x ( 3 x 3 det(a λi 3 λ 3 λ (3 λ 2 λ 3 (multiplicity 2 λ 3 Find a basis for the solution space of (A 3Iv In matrix form, we have v v 2 ( 2
13 So, both v and v 2 are arbitrary we have two linearly independent eigenvectors Let v r, v 2 s Then ( r v s r + s r + s The basis for the eigenspace gives us the two eigenvectors we seek: {}, Therefore, the general solution is x c e 3t + c 2 e 3t Case 2: algebraic multiplicity > geometric multiplicity For the case of a double root λ, this means that we will only obtain one eigenvector v, giving us only one solution, x e λt v We need to find a second solution, x 2 To do this, we will try a product of an exponential with a polynomial Let x 2 e λt (tu + w, where u and w are determined so that x 2 is a solution of x Ax Plug x 2 into the system: Divide both sides by e λt to get: Equating coefficients of like terms gives: x 2 e λt tu + e λt w x 2 λe λt tu + e λt u + λe λt w x 2 Ax 2 λe λt tu + e λt u + λe λt w A ( e λt tu + e λt w e λt tau + e λt Aw e λt (λtu + (u + λw e λt (tau + Aw λtu + (u + λw tau + Aw Au λu Aw λw + u (A λiw u 3
14 The first equation tells us that u is an eigenvector corresponding to λ, so let u v To summarize, we do the following: Determine the eigenvalues 2 Given an eigenvalue λ of multiplicity 2: (a If the dimension of the associated eigenspace is : i Find an eigenvector, v, giving one solution, x e λt v ii Find a second solution x 2 e λt (tv + w, where w solves (A λiw v (7 (b If the dimension of the associated eigenspace is 2, follow the standard procedure 3 For eigenvalues of multiplicity, follow the standard procedure Example: Solve x ( 4 x, x( 4 8 General solution: det(a λi λ λ λ(8 λ + 6 λ 2 8λ + 6 ( 2 (λ 4 2 λ 4 (multiplicity 2 λ 4 Find a basis for the solution space of (A 4Iv In matrix form, we have 4 4 v 4 4 This gives v 2 4v + 4v 2 v v 2 ( d v 2 is arbitrary, so let v 2 v So, ( ( v and x e 4t is one solution Next, find a second solution, x 2 e 4t (tv w, 4
15 where w solves (A 4Iw v, or 4 4 w 4 4 Solving, ( So, w 2 r 2 r 2 +r 4w 4w 2 w 2 is arbitrary, so let w 2 w 4 So, and General solution: x c x + c 2 x 2, so Find c and c 2 : x( ( w 4, ( x 2 e 4t (tv + w ( ( e (t 4t + 4 ( t + e 4t 4 t x c e 4t ( ( 2 So, c 2 and c + 4 c 2 c 2 4c 8 So, x 2e 4t ( ( t + + 8e 4t 4 t ( 4 4 ( t + + c 2 e 4t 4 t ( c 2 + c 2 4 5
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 informationGeneralized Eigenvectors and Jordan Form
Generalized Eigenvectors and Jordan Form We have seen that an n n matrix A is diagonalizable precisely when the dimensions of its eigenspaces sum to n. So if A is not diagonalizable, there is at least
More 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 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 informationAPPM 2360: Midterm exam 3 April 19, 2017
APPM 36: Midterm exam 3 April 19, 17 On the front of your Bluebook write: (1) your name, () your instructor s name, (3) your lecture section number and (4) a grading table. Text books, class notes, cell
More informationChapter 7. Canonical Forms. 7.1 Eigenvalues and Eigenvectors
Chapter 7 Canonical Forms 7.1 Eigenvalues and Eigenvectors Definition 7.1.1. Let V be a vector space over the field F and let T be a linear operator on V. An eigenvalue of T is a scalar λ F such that there
More information1. In this problem, if the statement is always true, circle T; otherwise, circle F.
Math 1553, Extra Practice for Midterm 3 (sections 45-65) Solutions 1 In this problem, if the statement is always true, circle T; otherwise, circle F a) T F If A is a square matrix and the homogeneous equation
More informationConstant coefficients systems
5.3. 2 2 Constant coefficients systems Section Objective(s): Diagonalizable systems. Real Distinct Eigenvalues. Complex Eigenvalues. Non-Diagonalizable systems. 5.3.. Diagonalizable Systems. Remark: We
More information4. Linear Systems. 4A. Review of Matrices ) , show that AB BA (= A A A).
4A-. Verify that 4A-2. If A = 2 3 2 2 4A-3. Calculate A if A = and A A = I. 4. Linear Systems 4A. Review of Matrices 2 2 and B = = 3 2 6, show that AB BA. 2 2 2, and check your answer by showing that AA
More informationMath Matrix Algebra
Math 44 - Matrix Algebra Review notes - 4 (Alberto Bressan, Spring 27) Review of complex numbers In this chapter we shall need to work with complex numbers z C These can be written in the form z = a+ib,
More informationMATH 304 Linear Algebra Lecture 23: Diagonalization. Review for Test 2.
MATH 304 Linear Algebra Lecture 23: Diagonalization. Review for Test 2. Diagonalization Let L be a linear operator on a finite-dimensional vector space V. Then the following conditions are equivalent:
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 informationFind the general solution of the system y = Ay, where
Math Homework # March, 9..3. Find the general solution of the system y = Ay, where 5 Answer: The matrix A has characteristic polynomial p(λ = λ + 7λ + = λ + 3(λ +. Hence the eigenvalues are λ = 3and λ
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 information4. Linear Systems. 4A. Review of Matrices A-1. Verify that =
4. Linear Systems 4A. Review of Matrices 0 2 0 0 0 4A-. Verify that 0 2 =. 2 3 2 6 0 2 2 0 4A-2. If A = and B =, show that AB BA. 3 2 4A-3. Calculate A 2 2 if A =, and check your answer by showing that
More informationLecture 3 Eigenvalues and Eigenvectors
Lecture 3 Eigenvalues and Eigenvectors Eivind Eriksen BI Norwegian School of Management Department of Economics September 10, 2010 Eivind Eriksen (BI Dept of Economics) Lecture 3 Eigenvalues and Eigenvectors
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 informationThird In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 3 December 2009 (1) [6] Given that 2 is an eigenvalue of the matrix
Third In-Class Exam Solutions Math 26, Professor David Levermore Thursday, December 2009 ) [6] Given that 2 is an eigenvalue of the matrix A 2, 0 find all the eigenvectors of A associated with 2. Solution.
More informationApplied Differential Equation. November 30, 2012
Applied Differential Equation November 3, Contents 5 System of First Order Linear Equations 5 Introduction and Review of matrices 5 Systems of Linear Algebraic Equations, Linear Independence, Eigenvalues,
More informationEigenvalues and Eigenvectors
Eigenvalues and Eigenvectors week -2 Fall 26 Eigenvalues and eigenvectors The most simple linear transformation from R n to R n may be the transformation of the form: T (x,,, x n ) (λ x, λ 2,, λ n x n
More informationSymmetric and anti symmetric matrices
Symmetric and anti symmetric matrices In linear algebra, a symmetric matrix is a square matrix that is equal to its transpose. Formally, matrix A is symmetric if. A = A Because equal matrices have equal
More informationRemark By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.
Sec 6 Eigenvalues and Eigenvectors Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called an eigenvalue of A if there is a nontrivial
More informationJordan normal form notes (version date: 11/21/07)
Jordan normal form notes (version date: /2/7) If A has an eigenbasis {u,, u n }, ie a basis made up of eigenvectors, so that Au j = λ j u j, then A is diagonal with respect to that basis To see this, let
More informationChapter 5. Eigenvalues and Eigenvectors
Chapter 5 Eigenvalues and Eigenvectors Section 5. Eigenvectors and Eigenvalues Motivation: Difference equations A Biology Question How to predict a population of rabbits with given dynamics:. half of the
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 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 informationDifferential Equations and Modeling
Differential Equations and Modeling Preliminary Lecture Notes Adolfo J. Rumbos c Draft date: March 22, 2018 March 22, 2018 2 Contents 1 Preface 5 2 Introduction to Modeling 7 2.1 Constructing Models.........................
More informationMTH 464: Computational Linear Algebra
MTH 464: Computational Linear Algebra Lecture Outlines Exam 4 Material Prof. M. Beauregard Department of Mathematics & Statistics Stephen F. Austin State University April 15, 2018 Linear Algebra (MTH 464)
More informationLinear Algebra and ODEs review
Linear Algebra and ODEs review Ania A Baetica September 9, 015 1 Linear Algebra 11 Eigenvalues and eigenvectors Consider the square matrix A R n n (v, λ are an (eigenvector, eigenvalue pair of matrix A
More informationHomogeneous Constant Matrix Systems, Part II
4 Homogeneous Constant Matrix Systems, Part II Let us now expand our discussions begun in the previous chapter, and consider homogeneous constant matrix systems whose matrices either have complex eigenvalues
More informationMath 240 Calculus III
Generalized Calculus III Summer 2015, Session II Thursday, July 23, 2015 Agenda 1. 2. 3. 4. Motivation Defective matrices cannot be diagonalized because they do not possess enough eigenvectors to make
More informationExam in TMA4110 Calculus 3, June 2013 Solution
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 8 Exam in TMA4 Calculus 3, June 3 Solution Problem Let T : R 3 R 3 be a linear transformation such that T = 4,
More 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 informationRemark 1 By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.
Sec 5 Eigenvectors and Eigenvalues In this chapter, vector means column vector Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called
More informationMATH 213 Linear Algebra and Ordinary Differential Equations Spring 2015 Study Sheet for Final Exam. Topics
MATH 213 Linear Algebra and Ordinary Differential Equations Spring 2015 Study Sheet for Final Exam This study sheet will not be allowed during the test. Books and notes will not be allowed during the test.
More informationCalculating determinants for larger matrices
Day 26 Calculating determinants for larger matrices We now proceed to define det A for n n matrices A As before, we are looking for a function of A that satisfies the product formula det(ab) = det A det
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 informationDiagonalization. P. Danziger. u B = A 1. B u S.
7., 8., 8.2 Diagonalization P. Danziger Change of Basis Given a basis of R n, B {v,..., v n }, we have seen that the matrix whose columns consist of these vectors can be thought of as a change of basis
More information1 Systems of Differential Equations
March, 20 7- Systems of Differential Equations Let U e an open suset of R n, I e an open interval in R and : I R n R n e a function from I R n to R n The equation ẋ = ft, x is called a first order ordinary
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 informationMath 310 Final Exam Solutions
Math 3 Final Exam Solutions. ( pts) Consider the system of equations Ax = b where: A, b (a) Compute deta. Is A singular or nonsingular? (b) Compute A, if possible. (c) Write the row reduced echelon form
More informationEigenspaces and Diagonalizable Transformations
Chapter 2 Eigenspaces and Diagonalizable Transformations As we explored how heat states evolve under the action of a diffusion transformation E, we found that some heat states will only change in amplitude.
More informationLU Factorization. A m x n matrix A admits an LU factorization if it can be written in the form of A = LU
LU Factorization A m n matri A admits an LU factorization if it can be written in the form of Where, A = LU L : is a m m lower triangular matri with s on the diagonal. The matri L is invertible and is
More information4. Linear transformations as a vector space 17
4 Linear transformations as a vector space 17 d) 1 2 0 0 1 2 0 0 1 0 0 0 1 2 3 4 32 Let a linear transformation in R 2 be the reflection in the line = x 2 Find its matrix 33 For each linear transformation
More informationSolutions of the Sample Problems for the Third In-Class Exam Math 246, Fall 2017, Professor David Levermore
Solutions of the Sample Problems for the Third In-Class Exam Math 6 Fall 07 Professor David Levermore Compute the Laplace transform of ft t e t ut from its definition Solution The definition of the Laplace
More informationINHOMOGENEOUS LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS
INHOMOGENEOUS LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS Definitions and a general fact If A is an n n matrix and f(t) is some given vector function, then the system of differential equations () x (t) Ax(t)
More informationReview. To solve ay + by + cy = 0 we consider the characteristic equation. aλ 2 + bλ + c = 0.
Review To solve ay + by + cy = 0 we consider the characteristic equation aλ 2 + bλ + c = 0. If λ is a solution of the characteristic equation then e λt is a solution of the differential equation. if there
More information(a) II and III (b) I (c) I and III (d) I and II and III (e) None are true.
1 Which of the following statements is always true? I The null space of an m n matrix is a subspace of R m II If the set B = {v 1,, v n } spans a vector space V and dimv = n, then B is a basis for V III
More informationMath 3313: Differential Equations Second-order ordinary differential equations
Math 3313: Differential Equations Second-order ordinary differential equations Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Mass-spring & Newton s 2nd law Properties
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 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 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 informationDefinition (T -invariant subspace) Example. Example
Eigenvalues, Eigenvectors, Similarity, and Diagonalization We now turn our attention to linear transformations of the form T : V V. To better understand the effect of T on the vector space V, we begin
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 informationAPPM 2360 Section Exam 3 Wednesday November 19, 7:00pm 8:30pm, 2014
APPM 2360 Section Exam 3 Wednesday November 9, 7:00pm 8:30pm, 204 ON THE FRONT OF YOUR BLUEBOOK write: () your name, (2) your student ID number, (3) lecture section, (4) your instructor s name, and (5)
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 informationThis is a closed book exam. No notes or calculators are permitted. We will drop your lowest scoring question for you.
Math 54 Fall 2017 Practice Final Exam Exam date: 12/14/17 Time Limit: 170 Minutes Name: Student ID: GSI or Section: This exam contains 9 pages (including this cover page) and 10 problems. Problems are
More informationLast lecture: Recurrence relations and differential equations. The solution to the differential equation dx
Last lecture: Recurrence relations and differential equations The solution to the differential equation dx = ax is x(t) = ce ax, where c = x() is determined by the initial conditions x(t) Let X(t) = and
More informationSolving a system by back-substitution, checking consistency of a system (no rows of the form
MATH 520 LEARNING OBJECTIVES SPRING 2017 BROWN UNIVERSITY SAMUEL S. WATSON Week 1 (23 Jan through 27 Jan) Definition of a system of linear equations, definition of a solution of a linear system, elementary
More informationEigenvalues, Eigenvectors, and Diagonalization
Math 240 TA: Shuyi Weng Winter 207 February 23, 207 Eigenvalues, Eigenvectors, and Diagonalization The concepts of eigenvalues, eigenvectors, and diagonalization are best studied with examples. We will
More informationMath 113 Solutions: Homework 8. November 28, 2007
Math 113 Solutions: Homework 8 November 28, 27 3) Define c j = ja j, d j = 1 j b j, for j = 1, 2,, n Let c = c 1 c 2 c n and d = vectors in R n Then the Cauchy-Schwarz inequality on Euclidean n-space gives
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 informationMATH 1553 PRACTICE MIDTERM 3 (VERSION B)
MATH 1553 PRACTICE MIDTERM 3 (VERSION B) Name Section 1 2 3 4 5 Total Please read all instructions carefully before beginning. Each problem is worth 10 points. The maximum score on this exam is 50 points.
More informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More information235 Final exam review questions
5 Final exam review questions Paul Hacking December 4, 0 () Let A be an n n matrix and T : R n R n, T (x) = Ax the linear transformation with matrix A. What does it mean to say that a vector v R n is an
More informationIMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET
IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET This is a (not quite comprehensive) list of definitions and theorems given in Math 1553. Pay particular attention to the ones in red. Study Tip For each
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 informationMATH 304 Linear Algebra Lecture 34: Review for Test 2.
MATH 304 Linear Algebra Lecture 34: Review for Test 2. Topics for Test 2 Linear transformations (Leon 4.1 4.3) Matrix transformations Matrix of a linear mapping Similar matrices Orthogonality (Leon 5.1
More informationJordan Canonical Form Homework Solutions
Jordan Canonical Form Homework Solutions For each of the following, put the matrix in Jordan canonical form and find the matrix S such that S AS = J. [ ]. A = A λi = λ λ = ( λ) = λ λ = λ =, Since we have
More 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 informationLecture 31. Basic Theory of First Order Linear Systems
Math 245 - Mathematics of Physics and Engineering I Lecture 31. Basic Theory of First Order Linear Systems April 4, 2012 Konstantin Zuev (USC) Math 245, Lecture 31 April 4, 2012 1 / 10 Agenda Existence
More information3. Identify and find the general solution of each of the following first order differential equations.
Final Exam MATH 33, Sample Questions. Fall 7. y = Cx 3 3 is the general solution of a differential equation. Find the equation. Answer: y = 3y + 9 xy. y = C x + C x is the general solution of a differential
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 informationMath 4B Notes. Written by Victoria Kala SH 6432u Office Hours: T 12:45 1:45pm Last updated 7/24/2016
Math 4B Notes Written by Victoria Kala vtkala@math.ucsb.edu SH 6432u Office Hours: T 2:45 :45pm Last updated 7/24/206 Classification of Differential Equations The order of a differential equation is the
More information12. Homogeneous Linear Systems of ODEs
2. Homogeneous Linear Systems of ODEs A system of first order differential equations is linear, homogeneous, and has constant coefficients if it has the specific form x = a x + a 2 x 2 + +a n x n, x 2
More informationMath 308 Final Exam Practice Problems
Math 308 Final Exam Practice Problems This review should not be used as your sole source for preparation for the exam You should also re-work all examples given in lecture and all suggested homework problems
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 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 informationIMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET
IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET This is a (not quite comprehensive) list of definitions and theorems given in Math 1553. Pay particular attention to the ones in red. Study Tip For each
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 informationMath 21b Final Exam Thursday, May 15, 2003 Solutions
Math 2b Final Exam Thursday, May 5, 2003 Solutions. (20 points) True or False. No justification is necessary, simply circle T or F for each statement. T F (a) If W is a subspace of R n and x is not in
More informationSolutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 28, 2011
Solutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 8, Solutions to Problems 5 are omitted since they involve topics no longer covered on the Comprehensive
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More informationDIAGONALIZABLE LINEAR SYSTEMS AND STABILITY. 1. Algebraic facts. We first recall two descriptions of matrix multiplication.
DIAGONALIZABLE LINEAR SYSTEMS AND STABILITY. 1. Algebraic facts. We first recall two descriptions of matrix multiplication. Let A be n n, P be n r, given by its columns: P = [v 1 v 2... v r ], where the
More informationMATH 423 Linear Algebra II Lecture 20: Geometry of linear transformations. Eigenvalues and eigenvectors. Characteristic polynomial.
MATH 423 Linear Algebra II Lecture 20: Geometry of linear transformations. Eigenvalues and eigenvectors. Characteristic polynomial. Geometric properties of determinants 2 2 determinants and plane geometry
More informationChapter 4 & 5: Vector Spaces & Linear Transformations
Chapter 4 & 5: Vector Spaces & Linear Transformations Philip Gressman University of Pennsylvania Philip Gressman Math 240 002 2014C: Chapters 4 & 5 1 / 40 Objective The purpose of Chapter 4 is to think
More informationEven-Numbered Homework Solutions
-6 Even-Numbered Homework Solutions Suppose that the matric B has λ = + 5i as an eigenvalue with eigenvector Y 0 = solution to dy = BY Using Euler s formula, we can write the complex-valued solution Y
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 informationTutorial 6 (week 6) Solutions
THE UNIVERSITY OF SYDNEY PURE MATHEMATICS Linear Mathematics 9 Tutorial 6 week 6 s Suppose that A and P are defined as follows: A and P Define a sequence of numbers { u n n } by u, u, u and for all n,
More information1 Last time: least-squares problems
MATH Linear algebra (Fall 07) Lecture Last time: least-squares problems Definition. If A is an m n matrix and b R m, then a least-squares solution to the linear system Ax = b is a vector x R n such that
More informationMath 2331 Linear Algebra
5. Eigenvectors & Eigenvalues Math 233 Linear Algebra 5. Eigenvectors & Eigenvalues Shang-Huan Chiu Department of Mathematics, University of Houston schiu@math.uh.edu math.uh.edu/ schiu/ Shang-Huan Chiu,
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationSolutions to Final Exam 2011 (Total: 100 pts)
Page of 5 Introduction to Linear Algebra November 7, Solutions to Final Exam (Total: pts). Let T : R 3 R 3 be a linear transformation defined by: (5 pts) T (x, x, x 3 ) = (x + 3x + x 3, x x x 3, x + 3x
More informationWhat is on this week. 1 Vector spaces (continued) 1.1 Null space and Column Space of a matrix
Professor Joana Amorim, jamorim@bu.edu What is on this week Vector spaces (continued). Null space and Column Space of a matrix............................. Null Space...........................................2
More informationDiagonalization of Matrix
of Matrix King Saud University August 29, 2018 of Matrix Table of contents 1 2 of Matrix Definition If A M n (R) and λ R. We say that λ is an eigenvalue of the matrix A if there is X R n \ {0} such that
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 informationAnnouncements Wednesday, November 7
Announcements Wednesday, November 7 The third midterm is on Friday, November 6 That is one week from this Friday The exam covers 45, 5, 52 53, 6, 62, 64, 65 (through today s material) WeBWorK 6, 62 are
More informationMATH 1553-C MIDTERM EXAMINATION 3
MATH 553-C MIDTERM EXAMINATION 3 Name GT Email @gatech.edu Please read all instructions carefully before beginning. Please leave your GT ID card on your desk until your TA scans your exam. Each problem
More informationdx n a 1(x) dy
HIGHER ORDER DIFFERENTIAL EQUATIONS Theory of linear equations Initial-value and boundary-value problem nth-order initial value problem is Solve: a n (x) dn y dx n + a n 1(x) dn 1 y dx n 1 +... + a 1(x)
More informationA = 3 1. We conclude that the algebraic multiplicity of the eigenvalues are both one, that is,
65 Diagonalizable Matrices It is useful to introduce few more concepts, that are common in the literature Definition 65 The characteristic polynomial of an n n matrix A is the function p(λ) det(a λi) Example
More information