Math 211 Math 211 Lecture #26 Solutions of a Planar System March - - PowerPoint PPT Presentation

math 211 math 211
SMART_READER_LITE
LIVE PREVIEW

Math 211 Math 211 Lecture #26 Solutions of a Planar System March - - PowerPoint PPT Presentation

1 Math 211 Math 211 Lecture #26 Solutions of a Planar System March 25, 2001 2 Polar Representation Polar Representation z = x + iy = r [cos + i sin ] . is the argument of z : tan = y/x. r = | z | . Eulers


slide-1
SLIDE 1

1

Math 211 Math 211

Lecture #26 Solutions of a Planar System March 25, 2001

slide-2
SLIDE 2

Return

2

Polar Representation Polar Representation

  • z = x + iy = r[cos θ + i sin θ].

⋄ θ is the argument of z: tan θ = y/x. ⋄ r = |z|.

  • Euler’s formula: eiθ = cos θ + i sin θ.

⋄ z = |z|eiθ. ⋄ z = |z|e−iθ.

slide-3
SLIDE 3

Return

3

Multiplication Multiplication

  • Two complex numbers

z = |z|eiθ and w = |w|eiφ

  • The product is

zw = |z|eiθ · |w|eiφ = |z||w|ei(θ+φ).

  • |zw| = |z||w|.
  • The argument of zw is the sum of the

arguments of z and w.

slide-4
SLIDE 4

Return

4

Complex Exponential Complex Exponential

Definition: For z = x + iy we define ez = ex+iy = ex · eiy = ex[cos y + i sin y]. Properties:

  • ez+w = ez · ew;

ez−w = ez · e−w = ez/ew

  • ez = ez
  • |ez| = ex = eRez
  • If λ is a complex number, then d

dteλt = λeλt

slide-5
SLIDE 5

Return

5

Complex Matrices Complex Matrices

Matrices (or vectors) with complex entries inherit many of the properties of complex numbers.

  • M = A + iB where A = ReM and

B = ImM are real matrices.

  • M = A − iB;

M = M ⇔ M is real.

  • ReM = 1

2(M + M);

ImM = 1

2i(M − M)

  • M + N = M + N
  • Mz = Mz
slide-6
SLIDE 6

Return

6

Procedure to Solve x′ = Ax Procedure to Solve x′ = Ax

  • Find the eigenvalues of A

⋄ the roots of p(λ) = det(A − λI) = 0

  • For each eigenvalue λ find the eigenspace

⋄ = null(A − λI)

  • If λ is an eigenvalue and v is an associated

eigenvector, x(t) = eλtv is a solution.

  • Show that n of these are linearly

independent.

slide-7
SLIDE 7

Return

7

Cases Cases

  • Distinct real eigenvalues.

⋄ In this case the method works as described.

  • Complex eigenvalues.

⋄ The method yields complex solutions.

  • Repeated eigenvalues.

⋄ The method does not always give enough solutions.

slide-8
SLIDE 8

Return Complex Matrices

8

Complex Eigenpairs Complex Eigenpairs

A a real matrix

  • λ a complex eigenvalue with associated

eigenvector w, so Aw = λw. Aw = Aw = Aw λw = λw

  • Aw = λw ⇒ Aw = λw ⇒ Aw = λw
  • ⇒ λ is an eigenvalue of A with associated

eigenvector w

slide-9
SLIDE 9

Return Complex Matrices

9

  • Thus complex eigenvalues come in

conjugate pairs λ and λ.

  • The associated eigenvectors also come in

conjugate pairs w and w.

  • λ = λ

⇒ w and w are linearly independent.

slide-10
SLIDE 10

Return Complex Matrices Complex eigenpairs

10

  • Complex exponential solutions

z(t) = eλtw and z(t) = eλtw.

  • z and z are linearly independent complex

valued solutions to x′ = Ax.

slide-11
SLIDE 11

Return Complex Matrices Complex eigenpairs Complex solutions

11

z(t) = x(t) + iy(t) & z(t) = x(t) − iy(t) x(t) = Re(z(t)) = z(t) + z(t) 2 y(t) = Im(z(t)) = z(t) − z(t) 2i

  • x(t) and y(t) are real valued solutions.
  • x(t) and y(t) are linearly independent.
slide-12
SLIDE 12

Return

12

Planar System x′ = Ax Planar System x′ = Ax

A = a11 a12 a21 a22

  • and

x(t) = x1(t) x2(t)

  • Characteristic polynomial:

p(λ) = λ2 − Tλ + D. ⋄ T = tr A = a11 + a22; D = det A ⋄ The eigenvalues of A are the roots of p.

slide-13
SLIDE 13

Return

13

Eigenvalues of A Eigenvalues of A

  • Roots of p(λ) = λ2 − Tλ + D = 0.

λ = T ± √ T 2 − 4D 2 .

  • Three cases:

⋄ 2 distinct real roots if T 2 − 4D > 0 ⋄ 2 complex conjugate roots if T 2 − 4D < 0 ⋄ Double real root if T 2 − 4D = 0

slide-14
SLIDE 14

Return

14

Example Example

x′ = Ax where A = −5 20 −2 7

  • p(λ) = λ2 − 2λ + 5.
  • Eigenvalues: λ = 1 + 2i

and λ = 1 − 2i

slide-15
SLIDE 15

Return Example

15

λ = 1 + 2i

  • A − λI =

−6 − 2i 20 −2 6 − 2i

  • .
  • Eigenvector: w =

3 − i 1

slide-16
SLIDE 16

Return Example λ

16

  • Complex Solutions

z(t) = eλtw = e(1+2i)t

  • 3 − i

1

  • z(t) = eλtw = e(1−2i)t
  • 3 + i

1

  • Real Solutions

x(t) = Re(z(t)) = et

  • 3 cos 2t + sin 2t

cos 2t

  • y(t) = Im(z(t)) = et
  • 3 sin 2t − cos 2t

sin 2t

slide-17
SLIDE 17

Return Fundamental set of solutions

17

Initial Value Problem Initial Value Problem

Solve x′ = Ax where A =

  • −5

20 −2 7

  • with the intial condition

x(0) =

  • 5

3

  • .
slide-18
SLIDE 18

Return Real solutions

18

Initial Value Problem Initial Value Problem

Solution is u(t) = 3et

  • 3 cos 2t + sin 2t

cos 2t

  • + 4et
  • 3 sin 2t − cos 2t

sin 2t

  • = et
  • 5 cos 2t + 15 sin 2t

3 cos 2t + 4 sin 2t

slide-19
SLIDE 19

Return

19

Summary — Complex Eigenvalues Summary — Complex Eigenvalues

Suppose A is a real 2 × 2 matrix with

  • complex conjugate eigenvalues λ and λ, and
  • associated nonzero eigenvectors w and w.

Then

  • z(t) = eλtw and z(t) = eλtw form a complex

valued fundamental set of solutions, and

  • x(t) = Re(z(t)) and y(t) = Im(z(t)) form a real

valued fundamental set of solutions.