Linear ODEs. Types of systems. Linear ODEs. Definition (Linear ODE) Linear ODEs. Existence of solutions to linear IVPs.

Size: px
Start display at page:

Download "Linear ODEs. Types of systems. Linear ODEs. Definition (Linear ODE) Linear ODEs. Existence of solutions to linear IVPs."

Transcription

1 Linear ODEs Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems p. 1 Linear ODEs Types of systems Definition (Linear ODE) A linear ODE is a ifferential equation taking the form x = A(t)x + B(t), t (LNH) where A(t) Mn(R) with continuous entries, B(t) R n with real value, continuous coefficients, an x R n. The associate IVP takes the form x = A(t)x + B(t) t (1) x(t0) = x0. x = A(t)x + B(t) is linear nonautonomous (A(t) epens on t) nonhomogeneous (also calle affine system). x = A(t)x is linear nonautonomous homogeneous. x = Ax + B, that is, A(t) A an B(t) B, is linear autonomous nonhomogeneous (or affine autonomous). x = Ax is linear autonomous homogeneous. If A(t + T ) = A(t) for some T > 0 an all t, then linear perioic. Linear ODEs p. 3 Linear ODEs p. 4

2 Existence an uniqueness of solutions Linear ODEs Existence of solutions to linear IVPs Resolvent matrix (Existence an Uniqueness) Solutions to (1) exist an are unique on the whole interval over which A an B are continuous. In particular, if A, B are constant, then solutions exist on R. Autonomous linear systems Existence of solutions to linear IVPs p. 6 The vector space of solutions Funamental matrix Consier the homogeneous system t x = A(t)x, with A(t) efine an continuous on an interval J. The set of solutions of (LH) forms an n-imensional vector space. (LH) Definition A set of n linearly inepenent solutions of (LH) on J, {φ1,..., φn}, is calle a funamental set of solutions of (LH) an the matrix Φ = [φ1 φ2... φn] is calle a funamental matrix of (LH). Existence of solutions to linear IVPs p. 7 Existence of solutions to linear IVPs p. 8

3 Funamental matrix solution Abel s formula Let X Mn(R) with entries [xij]. Define the erivative of X, X (or t X ) as t X (t) = [ t xij(t)]. The system of n 2 equations t X = A(t)X is calle a matrix ifferential equation. A funamental matrix Φ of (LH) satisfies the matrix equation X = A(t)X on the interval J.- If Φ is a solution of the matrix equation X = A(t)X on an interval J an τ J, then ( t ) etφ(t) = etφ(τ) exp tra(s)s τ for all t J. Existence of solutions to linear IVPs p. 9 Existence of solutions to linear IVPs p. 10 The resolvent matrix Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Definition (Resolvent matrix) Let t0 J an Φ(t) be a funamental matrix solution of (LH) on J. Since the columns of Φ are linearly inepenent, it follows that Φ(t0) is invertible. The resolvent (or state transition matrix, or principal funamental matrix) of (LH) is then efine as R(t, t0) = Φ(t)Φ(t0) 1. Autonomous linear systems Resolvent matrix p. 12

4 Proposition The resolvent matrix satisfies the Chapman-Kolmogorov ientities 1. R(t, t) = I, 2. R(t, s)r(s, u) = R(t, u), as well as the ientities 3. R(t, s) 1 = R(s, t), 4. s R(t, s) = R(t, s)a(s), 5. t R(t, s) = A(t)R(t, s). Proposition R(t, t0) is the only solution in Mn(K) of the initial value problem with M(t) Mn(K). M(t) = A(t)M(t) t M(t0) = I, Resolvent matrix p. 13 Resolvent matrix p. 14 A variation of constants formula The solution to the IVP consisting of the linear homogeneous nonautonomous system (LH) with initial conition x(t0) = x0 is given by φ(t) = R(t, t0)x0. (Variation of constants formula) Consier the IVP x = A(t)x + g(t, x) x(t0) = x0, (2a) (2b) where g : R R n R n a smooth function, an let R(t, t0) be the resolvent associate to the homogeneous system x = A(t)x, with R efine on some interval J t0. Then the solution φ of (2) is given by t φ(t) = R(t, t0)x0 + R(t, s)g(φ(s), s)s, (3) on some subinterval of J. t0 Resolvent matrix p. 15 Resolvent matrix p. 16

5 Autonomous linear systems Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Consier the autonomous affine system x = Ax + B, t (A) an the associate homogeneous autonomous system t x = Ax. (L) Autonomous linear systems p. 18 Exponential of a matrix Properties of the matrix exponential Definition (Matrix exponential) Let A Mn(K) with K = R or C. The exponential of A, enote e At, is a matrix in Mn(K), efine by e At t k = I + k! Ak, k=1 where I is the ientity matrix in Mn(K). Φ(t) = e At is a funamental matrix for (L) for t R. The resolvent for (L) is given for t J by R(t, t0) = e A(t t0) = Φ(t t0). e At1 e At2 = e A(t1+t2) for all t1, t2 R. 1 Ae At = e At A for all t R. (e At ) 1 = e At for all t R. The unique solution φ of (L) with φ(t0) = x0 is given by φ(t) = e A(t t0) x0. Autonomous linear systems p. 19 Autonomous linear systems p. 20

6 Computing the matrix exponential Let P be a nonsingular matrix in Mn(R). We transform the IVP t x = Ax x(t0) = x0 using the transformation x = Py or y = P 1 x. The ynamics of y is y = (P 1 x) = P 1 x = P 1 Ax = P 1 APy (L IVP) We have thus transforme IVP (L IVP) into t y = P 1 APy y(t0) = P 1 x0 From the earlier result, we then know that the solution of (L IVP y) is given by ψ(t) = e P 1 AP(t t0) P 1 x0, an since x = Py, the solution to (L IVP) is given by φ(t) = Pe P 1 AP(t t0) P 1 x0. So everything epens on P 1 AP. (L IVP y) The initial conition is y0 = P 1 x0. Autonomous linear systems p. 21 Autonomous linear systems p. 22 Diagonalizable case Assume P nonsingular in Mn(R) such that λ1 0 P 1 AP = 0 λn with all eigenvalues λ1,..., λn ifferent. We have k e P 1AP t k λ1 0 = I + k! k=1 0 λn Autonomous linear systems p. 23 Autonomous linear systems p. 24

7 For a (block) iagonal matrix M of the form m11 0 M = 0 mnn there hols m k 11 0 M k = 0 mnn k Therefore, e P 1AP t k λ k 1 0 = I + k! k=1 0 λ k n k=0 tk k! λk 1 0 = 0 e λ1t 0 = 0 e λnt k=0 tk k! λk n Autonomous linear systems p. 25 Autonomous linear systems p. 26 Noniagonalizable case An so the solution to (L IVP) is given by e λ1t 0 φ(t) = P P 1 x0. 0 e λnt The Joran canonical form is J0 0 P 1 AP = 0 Js so we use the same property as before (but with block matrices now), an e J0t 0 e P 1APt = 0 e Jst Autonomous linear systems p. 27 Autonomous linear systems p. 28

8 The first block in the Joran canonical form takes the form λ0 0 J0 = 0 λk an thus, as before, e λ0t 0 e J0t = 0 e λkt Other blocks Ji are written as Ji = λk+ii + Ni with I the ni ni ientity an Ni the ni ni nilpotent matrix Ni = λk+ii an Ni commute, an thus e Ji t = e λk+i t e Ni t Autonomous linear systems p. 29 Autonomous linear systems p. 30 Since Ni is nilpotent, Ni k = 0 for all k ni, an the series e Ni t terminates, an t 1 t n i 1 (ni 1)! e Ji t = e λk+i t t 0 1 n i 2 (ni 2)! 0 1 For all (t0, x0) R R n, there is a unique solution x(t) to (L IVP) efine for all t R. Each coorinate function of x(t) is a linear combination of functions of the form t k e αt cos(βt) an t k e αt sin(βt) where α + iβ is an eigenvalue of A an k is less than the algebraic multiplicity of the eigenvalue. Autonomous linear systems p. 31 Autonomous linear systems p. 32

9 Fixe points (equilibria) Orbits, limit sets Orbits an limit sets are efine as for maps. Definition A fixe point (or equilibrium point, or critical point) of an autonomous ifferential equation x = f (x) is a point p such that f (p) = 0. For a nonautonomous ifferential equation x = f (t, x), a fixe point satisfies f (t, p) = 0 for all t. A fixe point is a solution. For the equation x = f (x), the subset {x(t), t I }, where I is the maximal interval of existence of the solution, is an orbit. If the maximal solution x(t, x0) of x = f (x) is efine for all t 0, where f is Lipschitz on an open subset V of R n, then the omega limit set of x0 is the subset of V efine by ( ) ω(x0) = {x(t, x0) : t τ} V }. τ=0 Proposition A point q is in ω(x0) iff there exists a sequence {tk} such that limk tk = an limk x(tk, x0) = q V. Autonomous linear systems p. 33 Autonomous linear systems p. 34 Contracting linear equation Definition (Liapunov stable orbit) The orbit of a point p is Liapunov stable for a flow φt if, given ε > 0, there exists δ > 0 such that (x, p) < δ implies that (φt(x), φt(p)) < ε for all t 0. If p is a fixe point, then this is written (φt(x), p) < ε. Definition (Asymptotically stable orbit) The orbit of a point p is asymptotically stable (or attracting) for a flow φt if it is Liapunov stable, an there exists δ1 > 0 such that (x, p) < δ1 implies that limt (φt(x), φt(p)) = 0. If p is a fixe point, then it is asymptotically stable if it is Liapunov stable an there exists δ1 > 0 such that (x, p) < δ1 implies that ω(x) = {p}. Let A Mn(R), an consier the equation (L). Then the following conitions are equivalent. 1. There is a norm A on R n an a constant a > 0 such that for any x0 R n an all t 0, e At x0 A e at x0 A. 2. There is a norm B on R n an constants a > 0 an C 1 such that for any x0 R n an all t 0, e At x0 B Ce at x0 B. 3. All eigenvalues of A have negative real parts. In that case, the origin is a sink or attracting, the flow is a contraction (antonyms source, repelling an expansion). Autonomous linear systems p. 35 Autonomous linear systems p. 36

10 Hyperbolic linear equation Definition (Stable eigenspace) The stable eigenspace of A Mn(R) is E s = span{v : v generalize eigenvector for eigenvalue λ, with R(λ) < 0} Definition The linear ifferential equation (L) is hyperbolic if A has no eigenvalue with zero real part. Definition (Center eigenspace) The center eigenspace of A Mn(R) is E c = span{v : v generalize eigenvector for eigenvalue λ, with R(λ) = 0} Definition (Unstable eigenspace) The unstable eigenspace of A Mn(R) is Autonomous linear systems p. 37 We can write R n = E s E u +E c, an in the case that E c =, then R n = E s E u is calle a hyperbolic splitting. The symbol stans for irect sum. Definition (Direct sum) Let U, V be two subspaces of a vector space X. Then the span of U an V is efine by u + v for u U an v V. If U an V are isjoint except for 0, then the span of U an V is calle the irect sum of U an V, an is enote U V. Autonomous linear systems p. 39 E u = span{v : v generalize eigenvector for eigenvalue λ, with R(λ) > 0} Autonomous linear systems p. 38 Trichotomy Define V s = {v : there exists a > 0 an C 1 such that e At v Ce at v for t 0}. V u = {v : there exists a > 0 an C 1 such that e At v Ce a t v for t 0}. V c = {v : for all a > 0, e At v e a t 0 as t ± }. The following are true. 1. The subspaces E s, E u an E c are invariant uner the flow e At. 2. There hols that E s = V s, E u = V u an E c = V c, an thus e At E u is an exponential expansion, eat E s is an exponential contraction, an e At E c grows subexponentially as t ±. Autonomous linear systems p. 40

11 Topologically conjugate linear ODEs Definition (Topologically conjugate flows) Let φt an ψt be two flows on a space M. φt an ψt are topologically conjugate if there exists an homeomorphism h : M M such that for all x M an all t R. h φt(x) = ψt h(x), Definition (Topologically equivalent flows) Let φt an ψt be two flows on a space M. φt an ψt are topologically equivalent if there exists an homeomorphism h : M M an a function α : R M R such that Let A, B Mn(R). 1. If all eigenvalues of A an B have negative real parts, then the linear flows e At an e Bt are topologically conjugate. 2. Assume that the system is hyperbolic, an that the imension of the stable eigenspace of A is equal to the imension of the eigenspace of B. Then the linear flows e At an e Bt are topologically conjugate. h φ α(t+s,x) (x) = ψt h(x), for all x M an all t R, an where α(t, x) is monotonically increasing in t for each x an onto all of R. Autonomous linear systems p. 41 Autonomous linear systems p. 42 Let A, B Mn(R). Assume that e At an e Bt are linearly conjugate, i.e., there exists M with e Bt = Me At M 1. Then A an B have the same eigenvalues. Autonomous linear systems p. 43

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems Linear ODEs p. 1 Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Linear ODEs Definition (Linear ODE) A linear ODE is a differential equation taking the form

More information

Outline. Summary of topics relevant for the final. Outline of this part. Scalar difference equations. General theory of ODEs. Linear ODEs.

Outline. Summary of topics relevant for the final. Outline of this part. Scalar difference equations. General theory of ODEs. Linear ODEs. Outline Scalar ifference equations Summary of topics relevant for the final General theory of ODEs Linear ODEs Linear maps Analysis near fixe points (linearization) Bifurcations How to analyze a system

More information

Summary of topics relevant for the final. p. 1

Summary of topics relevant for the final. p. 1 Summary of topics relevant for the final p. 1 Outline Scalar difference equations General theory of ODEs Linear ODEs Linear maps Analysis near fixed points (linearization) Bifurcations How to analyze a

More information

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. Uniqueness for solutions of ifferential equations. We consier the system of ifferential equations given by x = v( x), () t with a given initial conition

More information

The goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T

The goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T 1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.

More information

Examples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling

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

Lecture Notes for Math 524

Lecture Notes for Math 524 Lecture Notes for Math 524 Dr Michael Y Li October 19, 2009 These notes are based on the lecture notes of Professor James S Muldowney, the books of Hale, Copple, Coddington and Levinson, and Perko They

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

MA5510 Ordinary Differential Equation, Fall, 2014

MA5510 Ordinary Differential Equation, Fall, 2014 MA551 Ordinary Differential Equation, Fall, 214 9/3/214 Linear systems of ordinary differential equations: Consider the first order linear equation for some constant a. The solution is given by such that

More information

Constant coefficients systems

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

7 Planar systems of linear ODE

7 Planar systems of linear ODE 7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

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

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

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

ODEs Cathal Ormond 1

ODEs Cathal Ormond 1 ODEs Cathal Ormond 2 1. Separable ODEs Contents 2. First Order ODEs 3. Linear ODEs 4. 5. 6. Chapter 1 Separable ODEs 1.1 Definition: An ODE An Ordinary Differential Equation (an ODE) is an equation whose

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

= A(λ, t)x. In this chapter we will focus on the case that the matrix A does not depend on time (so that the ODE is autonomous):

= A(λ, t)x. In this chapter we will focus on the case that the matrix A does not depend on time (so that the ODE is autonomous): Chapter 2 Linear autonomous ODEs 2 Linearity Linear ODEs form an important class of ODEs They are characterized by the fact that the vector field f : R m R p R R m is linear at constant value of the parameters

More information

Problem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs

Problem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable

More information

1. Find the solution of the following uncontrolled linear system. 2 α 1 1

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

Ordinary Differential Equation Theory

Ordinary Differential Equation Theory Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 1. Linear systems Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Linear systems 1.1 Differential Equations 1.2 Linear flows 1.3 Linear maps

More information

Slide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)

Slide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13) Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to

More information

Department of Mathematics IIT Guwahati

Department of Mathematics IIT Guwahati Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

Remark By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.

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

LINEAR ALGEBRA 1, 2012-I PARTIAL EXAM 3 SOLUTIONS TO PRACTICE PROBLEMS

LINEAR ALGEBRA 1, 2012-I PARTIAL EXAM 3 SOLUTIONS TO PRACTICE PROBLEMS LINEAR ALGEBRA, -I PARTIAL EXAM SOLUTIONS TO PRACTICE PROBLEMS Problem (a) For each of the two matrices below, (i) determine whether it is diagonalizable, (ii) determine whether it is orthogonally diagonalizable,

More information

5.) For each of the given sets of vectors, determine whether or not the set spans R 3. Give reasons for your answers.

5.) For each of the given sets of vectors, determine whether or not the set spans R 3. Give reasons for your answers. Linear Algebra - Test File - Spring Test # For problems - consider the following system of equations. x + y - z = x + y + 4z = x + y + 6z =.) Solve the system without using your calculator..) Find the

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

Solution of Linear State-space Systems

Solution of Linear State-space Systems Solution of Linear State-space Systems Homogeneous (u=0) LTV systems first Theorem (Peano-Baker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state

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

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

NOTES ON LINEAR ODES

NOTES ON LINEAR ODES NOTES ON LINEAR ODES JONATHAN LUK We can now use all the discussions we had on linear algebra to study linear ODEs Most of this material appears in the textbook in 21, 22, 23, 26 As always, this is a preliminary

More information

x 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7

x 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7 Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)

More information

10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections )

10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections ) c Dr. Igor Zelenko, Fall 2017 1 10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections 7.2-7.4) 1. When each of the functions F 1, F 2,..., F n in right-hand side

More information

Definition: An n x n matrix, "A", is said to be diagonalizable if there exists a nonsingular matrix "X" and a diagonal matrix "D" such that X 1 A X

Definition: An n x n matrix, A, is said to be diagonalizable if there exists a nonsingular matrix X and a diagonal matrix D such that X 1 A X DIGONLIZTION Definition: n n x n matrix, "", is said to be diagonalizable if there exists a nonsingular matrix "X" and a diagonal matrix "D" such that X X D. Theorem: n n x n matrix, "", is diagonalizable

More information

UNIVERSITY OF MANITOBA

UNIVERSITY OF MANITOBA DATE: May 8, 2015 Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones

More information

Jointly continuous distributions and the multivariate Normal

Jointly continuous distributions and the multivariate Normal Jointly continuous istributions an the multivariate Normal Márton alázs an álint Tóth October 3, 04 This little write-up is part of important founations of probability that were left out of the unit Probability

More information

Remark 1 By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.

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

Advanced Control Theory

Advanced Control Theory State Space Solution and Realization chibum@seoultech.ac.kr Outline State space solution 2 Solution of state-space equations x t = Ax t + Bu t First, recall results for scalar equation: x t = a x t + b

More 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

In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.

In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation. 1 2 Linear Systems In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation 21 Matrix ODEs Let and is a scalar A linear function satisfies Linear superposition ) Linear

More information

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the

More information

Differential Equations Practice: 2nd Order Linear: Nonhomogeneous Equations: Undetermined Coefficients Page 1

Differential Equations Practice: 2nd Order Linear: Nonhomogeneous Equations: Undetermined Coefficients Page 1 Differential Equations Practice: 2nd Order Linear: Nonhomogeneous Equations: Undetermined Coefficients Page 1 Questions Example (3.5.3) Find a general solution of the differential equation y 2y 3y = 3te

More information

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through

More information

Linear Systems Notes for CDS-140a

Linear Systems Notes for CDS-140a Linear Systems Notes for CDS-140a October 27, 2008 1 Linear Systems Before beginning our study of linear dynamical systems, it only seems fair to ask the question why study linear systems? One might hope

More information

Euler Equations: derivation, basic invariants and formulae

Euler Equations: derivation, basic invariants and formulae Euler Equations: erivation, basic invariants an formulae Mat 529, Lesson 1. 1 Derivation The incompressible Euler equations are couple with t u + u u + p = 0, (1) u = 0. (2) The unknown variable is the

More information

Linear ODE s with periodic coefficients

Linear ODE s with periodic coefficients Linear ODE s with periodic coefficients 1 Examples y = sin(t)y, solutions Ce cos t. Periodic, go to 0 as t +. y = 2 sin 2 (t)y, solutions Ce t sin(2t)/2. Not periodic, go to to 0 as t +. y = (1 + sin(t))y,

More information

We shall finally briefly discuss the generalization of the solution methods to a system of n first order differential equations.

We shall finally briefly discuss the generalization of the solution methods to a system of n first order differential equations. George Alogoskoufis, Dynamic Macroeconomic Theory, 2015 Mathematical Annex 1 Ordinary Differential Equations In this mathematical annex, we define and analyze the solution of first and second order linear

More information

Jordan normal form notes (version date: 11/21/07)

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

LECTURE NOTES ON DVORETZKY S THEOREM

LECTURE NOTES ON DVORETZKY S THEOREM LECTURE NOTES ON DVORETZKY S THEOREM STEVEN HEILMAN Abstract. We present the first half of the paper [S]. In particular, the results below, unless otherwise state, shoul be attribute to G. Schechtman.

More information

SOME LYAPUNOV TYPE POSITIVE OPERATORS ON ORDERED BANACH SPACES

SOME LYAPUNOV TYPE POSITIVE OPERATORS ON ORDERED BANACH SPACES Ann. Aca. Rom. Sci. Ser. Math. Appl. ISSN 2066-6594 Vol. 5, No. 1-2 / 2013 SOME LYAPUNOV TYPE POSITIVE OPERATORS ON ORDERED BANACH SPACES Vasile Dragan Toaer Morozan Viorica Ungureanu Abstract In this

More information

Rank, Trace, Determinant, Transpose an Inverse of a Matrix Let A be an n n square matrix: A = a11 a1 a1n a1 a an a n1 a n a nn nn where is the jth col

Rank, Trace, Determinant, Transpose an Inverse of a Matrix Let A be an n n square matrix: A = a11 a1 a1n a1 a an a n1 a n a nn nn where is the jth col Review of Linear Algebra { E18 Hanout Vectors an Their Inner Proucts Let X an Y be two vectors: an Their inner prouct is ene as X =[x1; ;x n ] T Y =[y1; ;y n ] T (X; Y ) = X T Y = x k y k k=1 where T an

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

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

ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky

ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky Equilibrium points and linearization Eigenvalue decomposition and modal form State transition matrix and matrix exponential Stability ELEC 3035 (Part

More information

University of Sydney MATH3063 LECTURE NOTES

University of Sydney MATH3063 LECTURE NOTES University of Sydney MATH3063 LECTURE NOTES Shane Leviton 6-7-206 Contents Introduction:... 6 Order of ODE:... 6 Explicit ODEs... 6 Notation:... 7 Example:... 7 Systems of ODEs... 7 Example 0: fx = 0...

More information

Linear Algebra II Lecture 13

Linear Algebra II Lecture 13 Linear Algebra II Lecture 13 Xi Chen 1 1 University of Alberta November 14, 2014 Outline 1 2 If v is an eigenvector of T : V V corresponding to λ, then v is an eigenvector of T m corresponding to λ m since

More information

Chapter 7. Canonical Forms. 7.1 Eigenvalues and Eigenvectors

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

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

Final A. Problem Points Score Total 100. Math115A Nadja Hempel 03/23/2017

Final A. Problem Points Score Total 100. Math115A Nadja Hempel 03/23/2017 Final A Math115A Nadja Hempel 03/23/2017 nadja@math.ucla.edu Name: UID: Problem Points Score 1 10 2 20 3 5 4 5 5 9 6 5 7 7 8 13 9 16 10 10 Total 100 1 2 Exercise 1. (10pt) Let T : V V be a linear transformation.

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

Darboux s theorem and symplectic geometry

Darboux s theorem and symplectic geometry Darboux s theorem an symplectic geometry Liang, Feng May 9, 2014 Abstract Symplectic geometry is a very important branch of ifferential geometry, it is a special case of poisson geometry, an coul also

More information

On the Inclined Curves in Galilean 4-Space

On the Inclined Curves in Galilean 4-Space Applie Mathematical Sciences Vol. 7 2013 no. 44 2193-2199 HIKARI Lt www.m-hikari.com On the Incline Curves in Galilean 4-Space Dae Won Yoon Department of Mathematics Eucation an RINS Gyeongsang National

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

Two Dimensional Linear Systems of ODEs

Two Dimensional Linear Systems of ODEs 34 CHAPTER 3 Two Dimensional Linear Sstems of ODEs A first-der, autonomous, homogeneous linear sstem of two ODEs has the fm x t ax + b, t cx + d where a, b, c, d are real constants The matrix fm is 31

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

Putzer s Algorithm. Norman Lebovitz. September 8, 2016

Putzer s Algorithm. Norman Lebovitz. September 8, 2016 Putzer s Algorithm Norman Lebovitz September 8, 2016 1 Putzer s algorithm The differential equation dx = Ax, (1) dt where A is an n n matrix of constants, possesses the fundamental matrix solution exp(at),

More information

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University. Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

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

model considered before, but the prey obey logistic growth in the absence of predators. In

model considered before, but the prey obey logistic growth in the absence of predators. In 5.2. First Orer Systems of Differential Equations. Phase Portraits an Linearity. Section Objective(s): Moifie Preator-Prey Moel. Graphical Representations of Solutions. Phase Portraits. Vector Fiels an

More information

Half of Final Exam Name: Practice Problems October 28, 2014

Half of Final Exam Name: Practice Problems October 28, 2014 Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half

More information

MATH 56A: STOCHASTIC PROCESSES CHAPTER 3

MATH 56A: STOCHASTIC PROCESSES CHAPTER 3 MATH 56A: STOCHASTIC PROCESSES CHAPTER 3 Plan for rest of semester (1) st week (8/31, 9/6, 9/7) Chap 0: Diff eq s an linear recursion (2) n week (9/11...) Chap 1: Finite Markov chains (3) r week (9/18...)

More information

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1 Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 PHYS 500 - Southern Illinois University September 8, 2016 PHYS 500 - Southern Illinois University Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 September 8,

More information

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

SYMPLECTIC GEOMETRY: LECTURE 3

SYMPLECTIC GEOMETRY: LECTURE 3 SYMPLECTIC GEOMETRY: LECTURE 3 LIAT KESSLER 1. Local forms Vector fiels an the Lie erivative. A vector fiel on a manifol M is a smooth assignment of a vector tangent to M at each point. We think of M as

More information

BASIC ALGORITHMS IN LINEAR ALGEBRA. Matrices and Applications of Gaussian Elimination. A 2 x. A T m x. A 1 x A T 1. A m x

BASIC ALGORITHMS IN LINEAR ALGEBRA. Matrices and Applications of Gaussian Elimination. A 2 x. A T m x. A 1 x A T 1. A m x BASIC ALGORITHMS IN LINEAR ALGEBRA STEVEN DALE CUTKOSKY Matrices and Applications of Gaussian Elimination Systems of Equations Suppose that A is an n n matrix with coefficents in a field F, and x = (x,,

More information

1. The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ = Ax + Bu is

1. The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ = Ax + Bu is ECE 55, Fall 2007 Problem Set #4 Solution The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ Ax + Bu is x(t) e A(t ) x( ) + e A(t τ) Bu(τ)dτ () This formula is extremely important

More information

Differential Topology Final Exam With Solutions

Differential Topology Final Exam With Solutions Differential Topology Final Exam With Solutions Instructor: W. D. Gillam Date: Friday, May 20, 2016, 13:00 (1) Let X be a subset of R n, Y a subset of R m. Give the definitions of... (a) smooth function

More information

Many problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx

Many problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx Math 53 Notes on turm-liouville equations Many problems in physics, engineering, an chemistry fall in a general class of equations of the form w(x)p(x) u ] + (q(x) λ) u = w(x) on an interval a, b], plus

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

Applications of the Wronskian to ordinary linear differential equations

Applications of the Wronskian to ordinary linear differential equations Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.

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

Econ 204 Differential Equations. 1 Existence and Uniqueness of Solutions

Econ 204 Differential Equations. 1 Existence and Uniqueness of Solutions Econ 4 Differential Equations In this supplement, we use the methods we have developed so far to study differential equations. 1 Existence and Uniqueness of Solutions Definition 1 A differential equation

More information

7. Symmetric Matrices and Quadratic Forms

7. Symmetric Matrices and Quadratic Forms Linear Algebra 7. Symmetric Matrices and Quadratic Forms CSIE NCU 1 7. Symmetric Matrices and Quadratic Forms 7.1 Diagonalization of symmetric matrices 2 7.2 Quadratic forms.. 9 7.4 The singular value

More information

Math 1553, Introduction to Linear Algebra

Math 1553, Introduction to Linear Algebra Learning goals articulate what students are expected to be able to do in a course that can be measured. This course has course-level learning goals that pertain to the entire course, and section-level

More information

LINEAR ALGEBRA SUMMARY SHEET.

LINEAR ALGEBRA SUMMARY SHEET. LINEAR ALGEBRA SUMMARY SHEET RADON ROSBOROUGH https://intuitiveexplanationscom/linear-algebra-summary-sheet/ This document is a concise collection of many of the important theorems of linear algebra, organized

More information

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

Autonomous system = system without inputs

Autonomous system = system without inputs Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation

More information

Nonlinear Systems and Elements of Control

Nonlinear Systems and Elements of Control Nonlinear Systems and Elements of Control Rafal Goebel April 3, 3 Brief review of basic differential equations A differential equation is an equation involving a function and one or more of the function

More information

Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems

Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 45 Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems Peter J. Hammond latest revision 2017 September

More information

Definition (T -invariant subspace) Example. Example

Definition (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 information

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R

More information

Chapter III. Stability of Linear Systems

Chapter III. Stability of Linear Systems 1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,

More information

1. What is the determinant of the following matrix? a 1 a 2 4a 3 2a 2 b 1 b 2 4b 3 2b c 1. = 4, then det

1. What is the determinant of the following matrix? a 1 a 2 4a 3 2a 2 b 1 b 2 4b 3 2b c 1. = 4, then det What is the determinant of the following matrix? 3 4 3 4 3 4 4 3 A 0 B 8 C 55 D 0 E 60 If det a a a 3 b b b 3 c c c 3 = 4, then det a a 4a 3 a b b 4b 3 b c c c 3 c = A 8 B 6 C 4 D E 3 Let A be an n n matrix

More information

1 Systems of Differential Equations

1 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 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

systems of linear di erential If the homogeneous linear di erential system is diagonalizable,

systems of linear di erential If the homogeneous linear di erential system is diagonalizable, G. NAGY ODE October, 8.. Homogeneous Linear Differential Systems Section Objective(s): Linear Di erential Systems. Diagonalizable Systems. Real Distinct Eigenvalues. Complex Eigenvalues. Repeated Eigenvalues.

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

8 Periodic Linear Di erential Equations - Floquet Theory

8 Periodic Linear Di erential Equations - Floquet Theory 8 Periodic Linear Di erential Equations - Floquet Theory The general theory of time varying linear di erential equations _x(t) = A(t)x(t) is still amazingly incomplete. Only for certain classes of functions

More information