Linear ODEs. Types of systems. Linear ODEs. Definition (Linear ODE) Linear ODEs. Existence of solutions to linear IVPs.
|
|
- Geraldine Payne
- 5 years ago
- Views:
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 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 informationOutline. 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 informationSummary 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 informationSYSTEMS 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 informationThe 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 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 informationLecture 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 informationDiagonalization 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 informationMA5510 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 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 information7 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 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 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 informationEuler 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 informationODEs 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 information6 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):
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 informationProblem 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 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 informationOrdinary 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 informationB5.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 informationSlide10 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 informationDepartment 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 informationSturm-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 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 informationLINEAR 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 information5.) 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 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 informationSolution 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 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 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 informationNOTES 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 informationx 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 information10. 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 informationDefinition: 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 informationUNIVERSITY 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 informationJointly 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 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 informationAdvanced 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 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 informationIn 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 informationCHAPTER 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 informationDifferential 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 informationDesigning 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 informationLinear 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 informationEuler 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 informationLinear 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 informationWe 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 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 informationLECTURE 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 informationSOME 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 informationRank, 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 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 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 informationELEC 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 informationUniversity 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 informationLinear 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 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 information19 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 informationFinal 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 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 informationDarboux 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 informationOn 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 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 informationTwo 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 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 informationPutzer 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 informationLecture 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 informationu 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 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 informationmodel 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 informationHalf 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 informationMATH 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 informationTensors, 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 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 informationSYMPLECTIC 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 informationBASIC 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 information1. 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 informationDifferential 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 informationMany 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 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 informationApplications 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 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 informationEcon 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 information7. 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 informationMath 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 informationLINEAR 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 informationSection 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 informationAutonomous 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 informationNonlinear 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 informationLecture 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 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 informationApplied 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 informationChapter 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 information1. 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 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 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 informationsystems 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 informationTable 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 information8 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