Math 677. Fall 2009. Homework #2 Solutions.: 1, 2 1 1 3, 4 2 2 1, 2 1, 2 1, 2 1 T 1 T 2 T 2 T
Math 677. Fall 2009. Homework #2 Solutions.: 1, 2 1 1 3, 4 2 2 1, 2 1, 2 1, 2 1 T 1 T 2 T 2 T
Math 677. Fall 2009. Homework #2 Solutions.: 1, 2 1 1 3, 4 2 2 1, 2 1, 2 1, 2 1 T 1 T 2 T 2 T
Homework #2 Solutions.
Part I. Exercises are taken from ”Diffential Equations and Dynamical Sys-
tems” by Perko, 3rd edition.
Problem Set 6: # 4
The system ẋ = Ax with
−1 −1 0 0
1 −1 0 0
A=
0
0 0 −2
0 0 1 2
has eigenvalues λ1,2 = a1 ± ib1 = −1 ± i, λ3,4 = a2 ± ib2 = 1 ± i. Eigenvectors
can be chosen for example as w1,2 = u1,2 ± iv1,2 with u1 = (1, 0, 0, 0)T , v1 =
(0, −1, 0, 0)T , u2 = (0, 0, −2, 1)T , v2 = (0, 0, 0, 1)T , so that
0 1 0 0
−1 0 0 0
P = [v1 u1 v2 u2 ] =
0 0 0 −2
0 0 1 1
The solution of the IVP is then provided by
aj t cos bj t − sin bj t
x(t) = P diag{e }P −1 x0
sin bj t cos bj t
where j = 1, 2.
1
w1 = (0, 1, 0, 0)T , w2 = (0, −4, 1, 0)T , w3 = (0, 12, −3, 1)T are the generalized
eigenvectors produced by the procedure (A − λI)wi = wi−1 , with i = 1, 2, 3
and w0 = v. In the basis determined by P = [v, w1 , w2 , w3 ], the Jordan form
of A reduces to
2 1 0 0
0 2 1 0
J = P −1 AP = 0 0 2 1
0 0 0 2
Problem Set 9: # 3
0 2 0
ẋ = −2 0 0 x
2 0 6
Eigenvalues are λ1 = 6, λ2,3 = ±2i, with eigenvectors w1 = (0, 0, 1)T and
w2,3 = u ± iv = (3, 1, −1)T ± i(−1, 3, 0)T , for instance. Notice that w1
spans x3 -axis and u, v span x1 , x2 -plane. Since λ1 > 0 and Re(λ2,3 ) = 0,
E u = x3 -axis, E c = x1 , x2 -plane. The phase portrait is a helix around x3 -
axis, diverging away from the x1 , x2 -plane, where it forms concentric circles
upon projection.
Problem Set 9: # 6
See Notes, lecture 8.
Problem Set 10: # 2
1 1 t
ẋ = x+
0 −1 1
with x(0) = [1, 0]T , is a linear non-homogeneous system, which is solved via
Variations of Constants formula:
Z t
At At
x(t) = e x0 + e e−Aτ b(τ )dτ
0
Since t
−t
Atet e −e
e = 2 ,
0 e−t
t t et −e−t Z t −τ e−τ −eτ
e e e τ
x(t) = + 2 2
τ dτ
0 0 e −t
0 0 e 1