21.04.2015 Views

Lecture Notes in Differential Equations - Bruce E. Shapiro

Lecture Notes in Differential Equations - Bruce E. Shapiro

Lecture Notes in Differential Equations - Bruce E. Shapiro

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

308 LESSON 31. LINEAR SYSTEMS<br />

We can f<strong>in</strong>d v by substitution <strong>in</strong>t (31.36):<br />

u = 29 5 e2t + 1 5 e−3t (31.43)<br />

v = u ′ − u (31.44)<br />

= 58 5 e2t − 3 5 e−3t − 29<br />

5 e2t − 1 5 e−3t (31.45)<br />

= 29 5 e2t − 4 5 e−3t (31.46)<br />

The Matrix Exponential<br />

Def<strong>in</strong>ition 31.3. Let M be any square matrix. Then we def<strong>in</strong>e the exponential<br />

of the matrix as<br />

exp(M) = e M = I + M + 1 2 M2 + 1 3! M3 + 1 4! M4 + · · · (31.47)<br />

assum<strong>in</strong>g that the series converges.<br />

Theorem 31.4. The solution of y ′ = Ay with <strong>in</strong>itial conditions y(t 0 ) = y 0<br />

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

Proof. Use Picard iteration:<br />

Then<br />

and <strong>in</strong> general<br />

Φ 0 (t) = y 0 (31.48)<br />

Φ k (t) = y 0 +<br />

∫ t<br />

t 0<br />

AΦ k−1 (s)ds (31.49)<br />

∫ t<br />

Φ 1 = y 0 + AΦ 0 (s)ds<br />

t 0<br />

(31.50)<br />

⎡<br />

Φ k = ⎣I +<br />

∫ t<br />

= y 0 + Ay 0 ds (31.51)<br />

t 0<br />

= y 0 + Ay 0 (t − t 0 ) (31.52)<br />

= [I + A(t − t 0 )]y 0 (31.53)<br />

k∑<br />

j=1<br />

1<br />

k! Ak (t − t 0 ) k ⎤<br />

⎦ y 0 (31.54)

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

Saved successfully!

Ooh no, something went wrong!