Lecture Notes in Differential Equations - Bruce E. Shapiro

Lecture Notes in Differential Equations - Bruce E. Shapiro Lecture Notes in Differential Equations - Bruce E. Shapiro

bruce.shapiro.com
from bruce.shapiro.com More from this publisher
21.04.2015 Views

310 LESSON 31. LINEAR SYSTEMS Example 31.2. Solve the system ( y ′ 6 = Ay where y(0) = 5) and From theorem 31.4 y = e A(t−t0) y 0 = A = ( ) 1 1 4 −2 { [( 1 1 exp 4 −2 ) ]} ( 6 t 5 ) (31.65) (31.66) (31.67) The problem now is that we don’t know how to calculate the matrix exponential! We will continue this example after we have discussed some of its properties. Theorem 31.5. Properties of the Matrix Exponential. 1. If 0 is a square n × n matrix composed entirely of zeros, then e 0 = I (31.68) 2. If A and B are square n×n matrices that commute (i.e., if AB = BA), then e A+B = e A e B (31.69) 3. e A is invertible, and ( e A ) −1 = e −A (31.70) 4. If A is any n × n matrix and S is any non-singular n × n matrix, then S −1 e A S = exp ( S −1 AS ) (31.71) 5. Let D = diag(x 1 , ..., x n ) be a diagonal matrix. Then e D = diag(e x1 , ..., e xn ). (31.72) 6. If S −1 AS = D, where D is a diagonal matrix, then e A = Se D S −1 . 7. If A has n linearly independent eigenvectors x 1 , ..., x n with corresponding eigenvalues λ 1 , ..., λ n , then e A = UEU −1 , where U = (x 1 , ..., x n ) and E = diag ( e λ1 , ..., e λn ) . 8. The matrix exponential is differentiable and d dt eAt = Ae At .

311 The Jordan Form Let A be a square n × n matrix. A is diagonalizable if for some matrix U, U −1 AU = D is diagonal. A is diagonalizable if and only if it has n linearly independent eigenvectors v 1 , ..., v n , in which case U = ( v 1 · · · v n ) will diagonalize A. If λ 1 , ..., λ k are the eigenvalues of A with multiplicities m 1 , ..., m k (hence n = m 1 + · · · + m k ) then the Jordan Canonical Form of A is ⎛ ⎞ B 1 0 · · · 0 J = 0 B 2 . ⎜ ⎝ . ⎟ . .. 0 ⎠ 0 · · · 0 B k (31.73) where (a) if m i = 1, B i = λ i (a scalar or 1 × 1 matrix; and (b) if m i ≠ 1, B i is an m i ×m i submatrix with λ i repeated in every diagonal element, the number 1 in the supra-diagonal, and zeroes everywhere else, e.g., if m i = 3 then ⎛ λ i 1 0 ⎞ 0 0 λ i B i = ⎝ 0 λ i 1 ⎠ (31.74) The B i are called Jordan Blocks. For every square matrix A, there exists a matrix U such J = U −1 AU exists. If J is the Jordan form of a matrix, then Thus e A = Ue J U −1 . ⎛ ⎞ e B1 0 0 e J ⎜ = ⎝ . 0 .. ⎟ 0 ⎠ (31.75) 0 0 e B k Example 31.3. Evaluate the solution y = e A(t−t0) y 0 = { [( 1 1 exp 4 −2 ) ]} ( 6 t 5 ) (31.76) found in the previous example.

310 LESSON 31. LINEAR SYSTEMS<br />

Example 31.2. Solve the system<br />

(<br />

y ′ 6<br />

= Ay where y(0) =<br />

5)<br />

and<br />

From theorem 31.4<br />

y = e A(t−t0) y 0 =<br />

A =<br />

( ) 1 1<br />

4 −2<br />

{ [( 1 1<br />

exp<br />

4 −2<br />

) ]} ( 6<br />

t<br />

5<br />

)<br />

(31.65)<br />

(31.66)<br />

(31.67)<br />

The problem now is that we don’t know how to calculate the matrix exponential!<br />

We will cont<strong>in</strong>ue this example after we have discussed some of its<br />

properties.<br />

Theorem 31.5. Properties of the Matrix Exponential.<br />

1. If 0 is a square n × n matrix composed entirely of zeros, then<br />

e 0 = I (31.68)<br />

2. If A and B are square n×n matrices that commute (i.e., if AB = BA),<br />

then<br />

e A+B = e A e B (31.69)<br />

3. e A is <strong>in</strong>vertible, and<br />

(<br />

e<br />

A ) −1<br />

= e<br />

−A<br />

(31.70)<br />

4. If A is any n × n matrix and S is any non-s<strong>in</strong>gular n × n matrix, then<br />

S −1 e A S = exp ( S −1 AS ) (31.71)<br />

5. Let D = diag(x 1 , ..., x n ) be a diagonal matrix. Then<br />

e D = diag(e x1 , ..., e xn ). (31.72)<br />

6. If S −1 AS = D, where D is a diagonal matrix, then e A = Se D S −1 .<br />

7. If A has n l<strong>in</strong>early <strong>in</strong>dependent eigenvectors x 1 , ..., x n with correspond<strong>in</strong>g<br />

eigenvalues λ 1 , ..., λ n , then e A = UEU −1 , where U =<br />

(x 1 , ..., x n ) and E = diag ( e λ1 , ..., e λn )<br />

.<br />

8. The matrix exponential is differentiable and d dt eAt = Ae At .

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!