NOTES ON LINEAR NON-AUTONOMOUS SYSTEMS

Similar documents
6 Linear Equation. 6.1 Equation with constant coefficients

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

Chapter 1: Systems of Linear Equations

Math Ordinary Differential Equations

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

1 Determinants. 1.1 Determinant

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

Linear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0.

Math 60. Rumbos Spring Solutions to Assignment #17

Section 9.2: Matrices.. a m1 a m2 a mn

Determine whether the following system has a trivial solution or non-trivial solution:

More chapter 3...linear dependence and independence... vectors

Chapter 1: Systems of linear equations and matrices. Section 1.1: Introduction to systems of linear equations

Chapter 1. Vectors, Matrices, and Linear Spaces

Chapter 7. Linear Algebra: Matrices, Vectors,

MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2

Chapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form

Lecture 18: Section 4.3

1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )

Linear Systems and Matrices

Matrices and RRE Form

Periodic Linear Systems

Lecture Notes in Linear Algebra

Section 9.2: Matrices. Definition: A matrix A consists of a rectangular array of numbers, or elements, arranged in m rows and n columns.

Equality: Two matrices A and B are equal, i.e., A = B if A and B have the same order and the entries of A and B are the same.

13. Systems of Linear Equations 1

Cayley-Hamilton Theorem

Chapter 1 Matrices and Systems of Equations

II. Determinant Functions

MAT 2037 LINEAR ALGEBRA I web:

Algebraic Properties of Solutions of Linear Systems

LS.1 Review of Linear Algebra

Ma 227 Review for Systems of DEs

AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES

Chapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form

Fundamentals of Engineering Analysis (650163)

Linear Equations in Linear Algebra

Linear Algebra: Lecture notes from Kolman and Hill 9th edition.

MATH 2331 Linear Algebra. Section 1.1 Systems of Linear Equations. Finding the solution to a set of two equations in two variables: Example 1: Solve:

Characteristic Polynomial

Chapter 5. Linear Algebra. Sections A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form

Midterm 1 Review. Written by Victoria Kala SH 6432u Office Hours: R 12:30 1:30 pm Last updated 10/10/2015

6.4 BASIS AND DIMENSION (Review) DEF 1 Vectors v 1, v 2,, v k in a vector space V are said to form a basis for V if. (a) v 1,, v k span V and

LU Factorization. A m x n matrix A admits an LU factorization if it can be written in the form of A = LU

Glossary of Linear Algebra Terms. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Properties of the Determinant Function

Chapter 5. Linear Algebra. Sections A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form

Solutions to Homework 5 - Math 3410

MATRICES. knowledge on matrices Knowledge on matrix operations. Matrix as a tool of solving linear equations with two or three unknowns.

Gaussian Elimination -(3.1) b 1. b 2., b. b n

System of Linear Equations

Ch 10.1: Two Point Boundary Value Problems

Systems of Linear Equations and Matrices

Linear Algebra Primer

CHAPTER 3. Matrix Eigenvalue Problems

Systems of Linear Equations and Matrices

Linear Independence x

UNIT 1 DETERMINANTS 1.0 INTRODUCTION 1.1 OBJECTIVES. Structure

The Solution of Linear Systems AX = B

ELEMENTARY LINEAR ALGEBRA WITH APPLICATIONS. 1. Linear Equations and Matrices

Lecture 6: Spanning Set & Linear Independency

LINEAR ALGEBRA SUMMARY SHEET.

Introduction to Determinants

Section Gaussian Elimination

Linear Algebra: Matrix Eigenvalue Problems

Basic Concepts in Linear Algebra

Math113: Linear Algebra. Beifang Chen

MAT Linear Algebra Collection of sample exams

Review of Basic Concepts in Linear Algebra

Linear Algebra. Preliminary Lecture Notes

3 Matrix Algebra. 3.1 Operations on matrices

Introduction to Matrices

Elementary Matrices. MATH 322, Linear Algebra I. J. Robert Buchanan. Spring Department of Mathematics

Lecture Summaries for Linear Algebra M51A

Linear Equations in Linear Algebra

Determinants. Samy Tindel. Purdue University. Differential equations and linear algebra - MA 262

SOLUTION OF GENERALIZED LINEAR VECTOR EQUATIONS IN IDEMPOTENT ALGEBRA

ECON 186 Class Notes: Linear Algebra

Matrices: 2.1 Operations with Matrices

MTH 464: Computational Linear Algebra

c c c c c c c c c c a 3x3 matrix C= has a determinant determined by

Linear Independence. Linear Algebra MATH Linear Algebra LI or LD Chapter 1, Section 7 1 / 1

Row Space, Column Space, and Nullspace

IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET

Practical Linear Algebra: A Geometry Toolbox

Sec. 7.4: Basic Theory of Systems of First Order Linear Equations

Math 240 Calculus III

EIGENVALUES AND EIGENVECTORS 3

Linear Algebra. Preliminary Lecture Notes

Review of Linear Algebra

Chapter 2:Determinants. Section 2.1: Determinants by cofactor expansion

Review of Vectors and Matrices

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

Math 4377/6308 Advanced Linear Algebra

ORIE 6300 Mathematical Programming I August 25, Recitation 1

SPRING OF 2008 D. DETERMINANTS

ELEMENTARY LINEAR ALGEBRA

IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET

MATH 323 Linear Algebra Lecture 6: Matrix algebra (continued). Determinants.

Evaluating Determinants by Row Reduction

Transcription:

NOTES ON LINEAR NON-AUTONOMOUS SYSTEMS XIAO-BIAO LIN 1 General Linear Systems Consider the linear nonhomogeneous system (11) y = A(t)y g(t) where A(t) and g(t) are continuous on an interval I Theorem 11 If A(t), g(t) are continuous on some interval a t b, if a t b, and if η R n, then the system (11) has a unique solution φ(t) satisfying the initial condition φ(t ) = η and existing on the interval a t b 11 Linear Homogeneous Systems Consider the liner homogeneous system associated with (11) (12) y = A(t)y For a homogeneous system, φ(t) = is the only solution that satisfies φ(t ) = Moreover, if φ 1 and φ 2 are any solutions of (12) on an interval I, and c 1 and c 2 are any constants, then c 1 φ 1 c 2 φ 2 is again a solution of (12) Definition 11 A set of vectors v 1, v 2,, v k is linearly dependent if there exist scalers c 1, c 2,, c k, not all zero, such that the linear combination c 1 v 1 c 2 v 2 c k v k = A set of vectors v 1, v 2,, v k is linearly independent if it is not linearly dependent A set S of vectors is said to form a basis of a vector space V if it is linearly independent and if every vector in V can be expressed as a linear combination of vectors in S We can define the dimension of a particular vector space V to be the number of elements in any basis of V A vector space is called finite-dimensional if is has a finite basis Theorem 12 if the complex n n matrix A(t) is continuous on an interval I, then the solutions of the system (12) on I form a vector space of dimension n over the complex numbers We say that the linearly independent solutions φ 1, φ 2,, φ n form a fundamental set of solutions There are infinitely many different fundamental sets of solutions of (12) A matrix of n rows whose columns are solutions (12) is called a solution matrix An n n solution matrix whose columns form a fundamental set of solutions is called a fundamental matrix for (12) on I Denote the fundamental matrix formed from the solutions φ 1,, φ n by Φ The statement that every solution φ of (12) is the 1

2 XIAO-BIAO LIN linear combination of φ 1,, φ n for some unique choice of the constants c 1,, c n is simply that (13) φ(t) = Φ(t)c where Φ is the fundamental matrix solution and c is the column vector with the components c 1,, c n Theorem 13 If Φ is a solution matrix of (12) on I and if t is any point of I, then (detφ) = ( a jj (t))detφ, j=1 detφ(t) = detφ(t )exp[ t j=1 a jj (s)ds], for every t in I It follows that either detφ(t) for each t I or detφ(t) = for every t I Proof Because the column vectors of Φ is a solution of (12), we have (14) φ ij = a ik φ kj, i, j = 1,, n k=1 Therefore, φ 11 φ 12 φ 1n (detφ) φ 21 φ 22 φ 2n = φ n1 φ n2 φ nn φ 11 φ 12 φ 1n φ 21 φ 22 φ 2n φ n1 φ n2 φ nn φ 11 φ 12 φ 1n φ 21 φ 22 φ 2n φ n1 φ n2 φ nn

NOTES ON LINEAR NON-AUTONOMOUS SYSTEMS 3 Using (14), we have n k=1 a n 1kφ k1 k=1 a n 1kφ k2 k=1 a 1kφ kn (detφ) φ 21 φ 22 φ 2n = φ n1 φ n2 φ nn φ 11 φ 12 φ 1n n k=1 a n 2kφ k1 k=1 a n 2kφ k2 k=1 a 2kφ kn φ n1 φ n2 φ nn φ 11 φ 12 φ 1n φ 21 φ 22 φ 2n n k=1 a n nkφ k1 k=1 a n nkφ k2 k=1 a nkφ kn Using elementary row elimination, we find that the first row of the first determinant simplifies to a 11 φ 11 a 11 φ 12 a 11 φ 1n Similarly for the ith row of the ith determinant Thus (detφ) = a 11 detφ a 22 detφ a nn detφ for every t I This proves the first part of the theorem The rest of the proof follows by solving the scalar equation for detφ Theorem 14 A solution matrix Φ of (12) on an interval I is a fundamental matrix of (12) if and only if detφ(t) for every t I Theorem 15 If Φ is a fundamental matrix for (12) on I and C is a nonsingular constant matrix, then ΦC is also a fundamental matrix for (12) Every fundamental matrix of (12) is of this form for some nonsingular matrix C 12 Linear Nonhomogeneous Systems Consider the linear nonhomogeneous system as in (11) (11) y = A(t)y g(t) Suppose φ 1 and φ 2 are any two solutions of (11) on I Then φ 1 φ 2 is a solution of the associated homogeneous systems (12) on I By the remark following Theorem 12, there exists a constant vector c such that The general solutions for (11) are φ 1 φ 2 = Φc (15) ψ = Φc φ, where Φ is a fundamental matrix solution of (12), c is an arbitrary constant vector and φ is a particular solution of (11) If the initial condition ψ(t ) = η is given for t I, then, the constant vector c can be solved from the given vector η Note that the matrix Φ(t ) is nonsingular

4 XIAO-BIAO LIN Theorem 16 If Φ is a fundamental matrix of (12) on I, then the function ψ(t) = Φ(t) Φ 1 (s)g(s)ds t is the (unique) solution of (11) valid on I and satisfying the initial condition ψ(t ) = If the initial condition ψ(t ) = η is given, then φ(t) = φ h (t) ψ(t), where ψ is the solution given in Theorem 16, and ψ h is the solution of the homogeneous system (12) satisfying the initial condition φ h (t ) = η, the same η as the initial condition for φ φ h (t) = Φ(t)Φ 1 (t )η We have φ(t) = Φ(t)Φ 1 (t )η Φ(t) Φ 1 (s)g(s)ds t Let T (t, s) = Φ(t)Φ 1 (s) Then T (t, s) is a special fundamental matrix solution that satisfies T (s, s) = I, and is called the principal matrix solution We can write the solution to (11) as φ(t) = T (t, t )η T (t, s)g(s)ds t 2 Linear Systems with Periodic Coefficients 21 Homogeneous systems with Periodic Coefficients Consider the homogeneous system (12) y = A(t)y, where A(t) is a continuous periodic n n matrix of period ω, ie, A(t ω) = A(t), t R If Φ(t) is a fundamental matrix for (12), then Φ(t ω) is also a fundamental matrix of (12) Therefore there exists a nonsingular constant matrix C such that Φ(t ω) = Φ(t)C, < t < In the rest of this paper, we assume that Φ(t) satisfies the initial condition Φ() = I Then we have Φ(t ω) = Φ(t)Φ(ω) proof? The nonsingular matrix C = Φ(ω) is called the monodromy matrix Since C is nonsingular, there exists a matrix R such that C = exp(ωr) Note that (21) Φ(ω) = exp(ωr) Theorem 21 Let A(t) be a continuous periodic matrix of period ω and let Φ(t) be any fundamental matrix of (12) Then there exists a periodic nonsingular matrix P (t) of period ω and a constant matrix R such that Φ(t) = P (t)exp(tr) Remark 21 (1) Note that Φ(t) may not be periodic (2) If A(t) is a constant matrix, hence periodic of any period ω, then P (t) = I and R = A

NOTES ON LINEAR NON-AUTONOMOUS SYSTEMS 5 Proof Let R be as in (21) Define P (t) is nonsingular and P (t) = Φ(t)exp( tr), < t < P (t ω) = Φ(t ω)exp( (t ω)r) = Φ(t)Cexp( (t ω)r) = Φ(t)Cexp( ωr)exp( tr) = Φ(t)exp( tr) = P (t) The Floquet theorem can be used to transform (12) to a linear system with constant coefficients Let (22) y = P (t)u Since y is a solutioin and Φ(t) is a fundamental matrix to (12), there exits a constant c such that y = Φ(t)c = P (t)exp(tr)c = P (t)u, it follows that Therefore u = exp(tr)c (23) u = Ru Corollary 22 The change of variable y = P (t)u transforms the periodic system (12) to the system (23) with constant coefficients Corollary 23 A nontrivial solution φ(t) of (12) has the property φ(t ω) = kφ(t), < t <, where k is a constant, if and only if k is an eigenvalue of Φ(ω) = exp(ωr) The eigenvalues ρ i of the nonsingular matrix Φ(ω) = exp(ωr) is called the characteristic multipliers of the system (12), and the eigenvalues λ i of the matrix R is called the characteristic exponents of (12) One can show that ρ j = e ωλj, λ j = 1 ω log ρ j, (mod 2πi ), j = 1,, n ω Π n j=1ρ j = exp( tra(s)ds), λ j = 1 tra(s)ds, (mod 2πi ω ω ) j=1 Corollary 24 If the characteristic exponents of (12) have negative real parts (or equivalently, it the multipliers (12) have magnitude strictly less than 1), then all solutions of (12) approach zero as t

6 XIAO-BIAO LIN 22 Nonhomogeneous Systems with Periodic Coefficients We consider the nonhomogeneous system (11) y = A(t)y g(t), where we assume that A(t) and g(t) are continuous and periodic in t of the same period ω Theorem 25 A solution φ(t) of (11) is periodic of period ω if and only if φ(ω) = φ() Proof The condition is clearly necessary To prove the sufficiency, let ψ(t) = φ(t ω) Then φ and ψ are both solutions of (11) and ψ() = φ(ω) = φ() By the uniqueness theorem, φ(t) = ψ(t) = φ(t ω) for all t Theorem 26 The system (11) has a periodic solution of period ω for any periodic forcing vector g of period ω if and only if the homogeneous system (12) has no periodic solution of period ω except for the trivial solution y = Proof Since Φ(t) be a fundamental matrix of (12) with Φ() = I A solution of (11) has the form ψ(t) = Φ(t)ψ() Φ(t) Φ 1 (s)g(s)ds The solution ψ is periodic of and only if ψ() = ψ(ω) But ψ(ω) = Φ(ω)ψ() Φ(ω) Therefore, ψ() = ψ(ω) if and only if [I Φ(ω)]ψ() = Φ(ω) Φ 1 (s)g(s)ds Φ 1 (s)g(s)ds The homogeneous algebraic system must be solved for every periodic forcing term g This is possible if and only if det(i Φ(ω)) On the other hand, (12) has a periodic solution of period ω if and only if 1 is a multiplier of y = A(t)y Thus, det(i Φ(ω)) if and only if (12) has no nontrivial periodic solution of period ω Theorem 27 (Fredholm s alternative) If A and g are continuous and periodic in t with period ω, then (11) has a period ω solution if and only if z(t)g(t)dt =, for all period ω solutions z of the adjoint equation ż = za(t) Note that z(t) is a row vector in the adoint equation