ODE-ProjectODE-Project Complex Eigenvalues

Section 3.4 Complex Eigenvalues

Objectives

To understand and be able to apply Euler’s formula,

$e^{i β t} = \cos β t + i \sin β t .$
To understand that if a $2 \times 2$ matrix $A$ has two complex eigenvalues, $λ = α \pm i β,$ then the general solution will involve sines and cosines. Furthermore, the origin will be a spiral sink, a spiral source or a center.

Consider the following system,

\begin{matrix} (3.4.1) & (\begin{matrix} d x / d t \\ d y / d t \end{matrix}) = (\begin{matrix} - 3 & 1 \\ - 2 & - 1 \end{matrix}) (\begin{matrix} x \\ y \end{matrix}) \end{matrix}

The characteristic polynomial of the system (3.4.1) is

λ^{2} + 4 λ + 5 .

The roots of this polynomial are

λ_{1} = - 2 + i

and

λ_{2} = - 2 - i

with eigenvectors

v_{1} = (1, 1 + i)

and

v_{2} = (1, 1 - i),

respectively. It is clear from the phase portrait of the system that there are no straight-line solutions (Figure 3.4.1). However, we would like to have real solutions for a linear system with real coefficients.

🔗

Figure 3.4.1. A system with no straight-line solutions
🔗

Subsection 3.4.1 Complex Eigenvalues

Suppose that we have the system

🔗

(\begin{matrix} d x / d t \\ d y / d t \end{matrix}) = (\begin{matrix} 0 & β \\ - β & 0 \end{matrix}) (\begin{matrix} x \\ y \end{matrix}) = A (\begin{matrix} x \\ y \end{matrix}),

where

β \neq 0 .

The characteristic polynomial of this system is

det (A - λ I) = λ^{2} + β^{2},

and so we have imaginary eigenvalues

\pm i β .

To find the eigenvector corresponding to

λ = i β,

we must solve the system

🔗

(\begin{matrix} - i β & β \\ - β & - i β \end{matrix}) (\begin{matrix} x \\ y \end{matrix}) = (\begin{matrix} 0 \\ 0 \end{matrix});

however, this reduces to solving the equation

i β x = β y .

Thus, we can find a complex eigenvector

(1, i) .

Consequently,

🔗

x (t) = e^{i β t} (\begin{matrix} 1 \\ i \end{matrix})

is a complex solution to the system

x^{'} = A x .

The problem is that we have a real system of differential equations and would like real solutions. We can remedy the situation if we use Euler’s formula,

If you are unfamiliar with Euler’s formula, try expanding both sides as a power series to check that we do indeed have a correct identity.

🔗

e^{i β t} = \cos β t + i \sin β t .

🔗

Let us rewrite our solution as

🔗

\begin{aligned} x (t) & = e^{i β t} (\begin{array}{c} 1 \\ i \end{array}) \\ = (\begin{array}{c} \cos β t + i \sin β t \\ i (\cos β t + i \sin β t) \end{array}) \\ = (\begin{array}{c} \cos β t + i \sin β t \\ - \sin β t + i \cos β t \end{array}) \\ = (\begin{array}{c} \cos β t \\ - \sin β t \end{array}) + i (\begin{array}{c} \sin β t \\ \cos β t \end{array}) \end{aligned}

and consider the real and imaginary parts of the solution:

🔗

x_{Re} = (\begin{matrix} \cos β t \\ - \sin β t \end{matrix}) and x_{Im} = (\begin{matrix} \sin β t \\ \cos β t \end{matrix}) .

Since

🔗

\begin{aligned} x_{Re}^{'} (t) + i x_{Im}^{'} (t) & = x^{'} (t) \\ = A x (t) \\ = A (x_{Re} (t) + i x_{Im} (t)) \\ = A x_{Re} (t) + i A x_{Im} (t) . \end{aligned}

we know that

x_{Re}^{'} (t) = A x_{Re} (t)

and

x_{Im}^{'} (t) = A x_{Im} (t)

by setting the real and imaginary parts equal. Thus, both

x_{Re} (t)

and

x_{Im} (t)

are solutions to our system. Moreover, since

🔗

x_{Re} (0) = (\begin{matrix} 1 \\ 0 \end{matrix}) and x_{Im} (0) = (\begin{matrix} 0 \\ 1 \end{matrix}),

we know that the appropriate linear combination of

x_{Re} (t)

and

x_{Im} (t)

will provide a solution to any initial value problem.

🔗

We claim that

🔗

\begin{matrix} (3.4.2) & x (t) = c_{1} x_{Re} (t) + c_{2} x_{Im} (t) \end{matrix}

is a general solution to our system. That is, we must be able to write every solution of our system as a linear combination of

x_{Re} (t)

and

x_{Im} (t) .

🔗

y (t) = (\begin{matrix} u (t) \\ v (t) \end{matrix})

is another solution to the system

x^{'} = A x,

then

🔗

y^{'} (t) = (\begin{matrix} u^{'} (t) \\ v^{'} (t) \end{matrix}) = (\begin{matrix} 0 & β \\ - β & 0 \end{matrix}) (\begin{matrix} u (t) \\ v (t) \end{matrix}) = (\begin{matrix} β v (t) \\ - β u (t) \end{matrix}) .

In other words,

u^{'} (t) = β v (t)

and

v^{'} (t) = - β u (t) .

Now, define

f

🔗

f (t) = (u (t) + i v (t)) e^{i β t} .

The derivative of

f

🔗

\begin{aligned} f^{'} (t) & = (u^{'} (t) + i v^{'} (t)) e^{i β t} + i β (u (t) + i v (t)) e^{i β t} \\ = (β v (t) - i β u (t)) e^{i β t} + (i β u (t) + i^{2} β v (t)) e^{i β t} \\ = 0. \end{aligned}

Therefore,

f (t)

is a complex constant and

f (t) = (u (t) + i v (t)) e^{i β t} = a + b i .

We can now write

u (t) + i v (t) = (a + i b) e^{- i β t} .

Thus,

🔗

\begin{aligned} u (t) + i v (t) & = (a + i b) e^{- i β t} \\ = (a + i b) (\cos β t - i \sin β t) \\ = (a \cos β t + b \sin β t) + i (b \cos β t - a \sin β t) . \end{aligned}

Therefore,

🔗

\begin{aligned} u (t) & = a \cos β t + b \sin β t \\ v (t) & = b \cos β t - a \sin β t . \end{aligned}

Consequently,

🔗

\begin{aligned} (\begin{array}{c} u (t) \\ v (t) \end{array}) & = (\begin{array}{c} a \cos β t + b \sin β t \\ b \cos β t - a \sin β t \end{array}) \\ = a (\begin{array}{c} \cos β t \\ - \sin β t \end{array}) + b (\begin{array}{c} \sin β t \\ \cos β t \end{array}) \\ = a x_{Re} (t) + b x_{Im} (t) . \end{aligned}

🔗

Notice that the solutions

🔗

x (t) = c_{1} (\begin{matrix} \cos β t \\ - \sin β t \end{matrix}) + c_{2} (\begin{matrix} \sin β t \\ \cos β t \end{matrix})

are periodic with period

2 π / β .

This type of system is called a center.

🔗

Example 3.4.2.

Consider the initial value problem

\begin{aligned} \frac{d x}{d t} & = 2 y \\ \frac{d y}{d t} & = - 2 x \\ x (0) & = - 3 \\ y (0) & = 0. \end{aligned}

The eigenvalues of this system are

λ = \pm 2 i .

Therefore, the general solution to the system is

\begin{aligned} x (t) & = c_{1} \cos 2 t + c_{2} \sin 2 t \\ y (t) & = - c_{1} \sin 2 t + c_{2} \cos 2 t . \end{aligned}

Using the initial conditions to solve for

c_{1}

and

c_{2},

the solution to our initial value problem is

\begin{aligned} x (t) & = - 3 \cos 2 t \\ y (t) & = - 3 \sin 2 t . \end{aligned}

The phase portrait is a circle of radius 3 about the origin (Figure 3.4.3).

Figure 3.4.3. Phase portrait for a center

🔗

Subsection 3.4.2 Spiral Sinks and Sources

Now let us consider the system

x^{'} = A x,

where

🔗

A = (\begin{matrix} α & β \\ - β & α \end{matrix})

and

α

and

β

are nonzero real numbers. The characteristic equation of

A

🔗

λ^{2} - 2 α λ + (α^{2} + β^{2}) = 0,

so the eigenvalues are

λ = α \pm i β .

We can use the equation

🔗

(α - (α + i β)) x + β y = 0

to show that

(1, i)

is an eigenvector for

α + i β .

Therefore, we have a complex solution of the form

🔗

\begin{aligned} x (t) & = e^{(α + i β) t} (\begin{array}{c} 1 \\ i \end{array}) \\ = e^{α t} (\begin{array}{c} \cos β t \\ - \sin β t \end{array}) + i e^{α t} (\begin{array}{c} \sin β t \\ \cos β t \end{array}) \\ = x_{Re} (t) + i x_{Im} (t) . \end{aligned}

As before, both

🔗

x_{Re} (t) = e^{α t} (\begin{matrix} \cos β t \\ - \sin β t \end{matrix}) and x_{Im} (t) = e^{α t} (\begin{matrix} \sin β t \\ \cos β t \end{matrix})

are real solutions to

x^{'} = A x .

Furthermore, these solutions are linearly independent. Indeed,

x_{Re}

cannot be a multiple of

x_{Im}

for all values of

t .

Thus, we have the general solution

🔗

x (t) = c_{1} e^{α t} (\begin{matrix} \cos β t \\ - \sin β t \end{matrix}) + c_{2} e^{α t} (\begin{matrix} \sin β t \\ \cos β t \end{matrix}) .

The

e^{α t}

factor tells us that the solutions either spiral into the origin if

α < 0

or spiral out to infinity if

α > 0 .

In this case we say that the equilibrium points are spiral sinks and spiral sources, respectively.

🔗

Example 3.4.4.

Consider the initial value problem

\begin{aligned} \frac{d x}{d t} & = - x / 10 + y \\ \frac{d y}{d t} & = - x - y / 10 \\ x (0) & = 2 \\ y (0) & = 2. \end{aligned}

The matrix

(\begin{matrix} - 1 / 10 & 1 \\ - 1 & - 1 / 10 \end{matrix})

has eigenvalues

λ = - 1 / 10 \pm i .

The eigenvalue

λ = - 1 / 10 + i

has an eigenvector

v = (1, i) .

The complex solution of our system is

\begin{aligned} x (t) & = e^{(- 1 / 10 + i) t} (\begin{array}{c} 1 \\ i \end{array}) \\ = e^{- t / 10} e^{i t} (\begin{array}{c} 1 \\ i \end{array}) \\ = e^{- t / 10} (\cos t + i \sin t) (\begin{array}{c} 1 \\ i \end{array}) \\ = e^{- t / 10} (\begin{array}{c} \cos t + i \sin t \\ - \sin t + i \cos t \end{array}) \\ = e^{- t / 10} (\begin{array}{c} \cos t \\ - \sin t \end{array}) + i e^{- t / 10} (\begin{array}{c} \sin t \\ \cos t \end{array}) \end{aligned}

The real and imaginary parts of this solution are

x_{Re} (t) = e^{- t / 10} (\begin{matrix} \cos t \\ - \sin t \end{matrix}) and x_{Im} (t) = e^{- t / 10} (\begin{matrix} \sin t \\ \cos t \end{matrix}),

respectively. Thus, we have the general solution

x (t) = c_{1} e^{- t / 10} (\begin{matrix} \cos t \\ - \sin t \end{matrix}) + c_{2} e^{- t / 10} (\begin{matrix} \sin t \\ \cos t \end{matrix}) .

Applying our initial conditions, we can determine that

c_{1} = 2

and

c_{2} = 2;

hence, the solution to our initial value problem is

x (t) = 2 e^{- t / 10} (\begin{matrix} \cos t + \sin t \\ \cos t - \sin t \end{matrix}) .

The phase portrait of this solution indicates that we do indeed have a spiral sink (Figure 3.4.5).

Figure 3.4.5. Phase portrait for a spiral sink

Example 3.4.6.

The initial value problem

\begin{aligned} \frac{d x}{d t} & = x / 10 + y \\ \frac{d y}{d t} & = - x + y / 10 \\ x (0) & = 0 \\ y (0) & = 1 / 2. \end{aligned}

The matrix

(\begin{matrix} 1 / 10 & 1 \\ - 1 & 1 / 10 \end{matrix})

has an eigenvector

(1, - i)

with eigenvalue

λ = 1 / 10 - i .

Thus, the complex solution is

x (t) = e^{(1 / 10 - i) t} (\begin{matrix} 1 \\ - i \end{matrix}) .

Following the procedure that we used in the previous example, the solution to our initial value problem is

x (t) = \frac{1}{2} e^{t / 10} (\begin{matrix} \sin t \\ \cos t \end{matrix}),

and he phase portrait is a spiral source (Figure 3.4.7).

Figure 3.4.7. Phase portrait for a spiral sink

Activity 3.4.1. Systems with Complex Eigenvalues.

Consider the system

d x / d t = A x,

where

🔗

A = (\begin{matrix} 7 & 4 \\ - 4 & 7 \end{matrix})

🔗

(a)

Find the eigenvalues,

λ

and

\overset{―}{λ}

A .

🔗

(b)

Find eigenvectors,

v

and

\overset{―}{v}

for the eigenvalues

λ

and

\overset{―}{λ},

respectively.

🔗

(c)

Find the complex solution to the system

d x / d t = A x .

🔗

(d)

Find the real solution to the system

d x / d t = A x .

🔗

(e)

Is the origin a spiral source or a spiral sink? Sketch a solution curve in the

x y

-plane.

🔗

Subsection 3.4.3 Solving Systems with Complex Eigenvalues

Suppose that we have the linear system

x^{'} = A x,

where

🔗

A = (\begin{matrix} a & b \\ c & d \end{matrix}) .

The characteristic polynomial of

A

🔗

p (λ) = λ^{2} - (a + d) λ + (a d - b c) .

(a + d)^{2} - 4 (a d - b c) < 0,

then the eigenvalues of

A

are complex, and we cannot apply the strategy that we used to determine the general solution in the case of distinct real roots.

🔗

Example 3.4.8.

The system

x^{'} = A x,

where

A = (\begin{matrix} - 1 & - 2 \\ 4 & 3 \end{matrix}) .

The characteristic polynomial of

A

λ^{2} - 2 λ + 5

and so the eigenvalues are complex conjugates,

λ = 1 + 2 i

and

\overset{―}{λ} = 1 - 2 i .

It is easy to show that an eigenvector for

λ = 1 + 2 i

v = (1, - 1 - i) .

Recalling that

e^{i θ} = \cos θ + i \sin θ,

\begin{aligned} x (t) & = e^{(1 + 2 i) t} v \\ = e^{(1 + 2 i) t} (\begin{array}{c} 1 \\ - 1 - i \end{array}) \\ = e^{t} e^{2 i t} (\begin{array}{c} 1 \\ - 1 - i \end{array}) \\ = e^{t} (\cos 2 t + i \sin 2 t) (\begin{array}{c} 1 \\ - 1 - i \end{array}) \\ = e^{t} (\begin{array}{c} \cos 2 t + i \sin 2 t \\ (- 1 - i) (\cos 2 t + i \sin 2 t) \end{array}) \\ = e^{t} (\begin{array}{c} \cos 2 t + i \sin 2 t \\ (- \cos 2 t + \sin 2 t) + i (- \cos 2 t - \sin 2 t) \end{array}) \\ = e^{t} (\begin{array}{c} \cos 2 t \\ - \cos 2 t + \sin 2 t \end{array}) + i e^{t} (\begin{array}{c} \sin 2 t \\ - \cos 2 t - \sin 2 t \end{array}) \end{aligned}

is a complex solution to our system. Taking the real and imaginary parts of this solution, we obtain the general solution to our system

x (t) = c_{1} e^{t} (\begin{matrix} \cos 2 t \\ - \cos 2 t + \sin 2 t \end{matrix}) + c_{2} e^{t} (\begin{matrix} \sin 2 t \\ - \cos 2 t - \sin 2 t \end{matrix}) .

The nature of the equilibrium solution is determined by the factor

e^{α t}

in the solution. If

α < 0,

the equilibrium point is a spiral sink. If

α > 0,

the equilibrium point is a spiral source. If

α = 0,

the equilibrium point is a center.

🔗

Although we have outlined a procedure to find the general solution of

x^{'} = A x

A

has complex eigenvalues, we have not shown that this method will work in all cases. We will do so in Section 3.6.

🔗

Activity 3.4.2. Planar Systems with Complex Eigenvalues.

Consider the system

d x / d t = A x,

where

🔗

A = (\begin{matrix} 7 & - 4 \\ 10 & - 5 \end{matrix})

🔗

(a)

Find the eigenvalues,

λ

and

\overset{―}{λ}

A .

🔗

(b)

Find eigenvectors,

v

and

\overset{―}{v}

for the eigenvalues

λ

and

\overset{―}{λ},

respectively.

🔗

(c)

Find the complex solution to the system

d x / d t = A x .

🔗

(d)

Find the real solution to the system

d x / d t = A x .

🔗

(e)

Is the origin a spiral source or a spiral sink? Sketch a solution curve in the

x y

-plane.

🔗

Subsection 3.4.4 Important Lessons

If

$A = (\begin{matrix} α & β \\ - β & α \end{matrix}),$

then $A$ has two complex eigenvalues, $λ = α \pm i β .$ The general solution to the system $x^{'} = A x$ is

$x (t) = c_{1} e^{α t} (\begin{matrix} \cos β t \\ - \sin β t \end{matrix}) + c_{2} e^{α t} (\begin{matrix} \sin β t \\ \cos β t \end{matrix}) .$

If $α < 0,$ the equilibrium point is a spiral sink. If $α > 0,$ the equilibrium point is a spiral source.

🔗

Reading Questions 3.4.5 Reading Questions

1.

When are two complex numbers equal?

🔗

2.

What is Euler’s formula?

🔗

3.

For a

2 \times 2

linear system with complex eigenvalues, what are the different possibilities for solution curves?

🔗

Exercises 3.4.6 Exercises

Solving Linear Systems with Complex Eigenvalues.

Find the general solution of each of the linear systems in Exercise Group 3.4.6.1–8.

🔗

1.

\begin{aligned} x^{'} & = 2 x + 2 y \\ y^{'} & = - 4 x + 6 y \end{aligned}

🔗

2.

\begin{aligned} x^{'} & = 2 x - 5 y \\ y^{'} & = x - 2 y \end{aligned}

🔗

3.

\begin{aligned} x^{'} & = - 9 x + 26 y \\ y^{'} & = - 4 x + 11 y \end{aligned}

🔗

4.

\begin{aligned} x^{'} & = - 7 x + 26 y \\ y^{'} & = - 4 x + 13 y \end{aligned}

🔗

5.

\begin{aligned} x^{'} & = - 2 x + 13 y \\ y^{'} & = - 2 x + 8 y \end{aligned}

🔗

6.

\begin{aligned} x^{'} & = - 7 x + 13 y \\ y^{'} & = - 2 x + 3 y \end{aligned}

🔗

7.

\begin{aligned} x^{'} & = 18 x - 52 y \\ y^{'} & = 8 x - 22 y \end{aligned}

🔗

8.

\begin{aligned} x^{'} & = - 12 x + 26 y \\ y^{'} & = - 4 x + 8 y \end{aligned}

🔗

Solving Initial Value Problems.

Solve each of the following linear systems for the given initial values in Exercise Group 3.4.6.9–16.

🔗

9.

\begin{aligned} x^{'} & = 2 x + 2 y \\ y^{'} & = - 4 x + 6 y \\ x (0) & = 2 \\ y (0) & = - 3 \end{aligned}

🔗

10.

\begin{aligned} x^{'} & = 2 x - 5 y \\ y^{'} & = x - 2 y \\ x (0) & = 2 \\ y (0) & = 1 \end{aligned}

🔗

11.

\begin{aligned} x^{'} & = - 9 x + 26 y \\ y^{'} & = - 4 x + 11 y \\ x (0) & = 10 \\ y (0) & = 10 \end{aligned}

🔗

12.

\begin{aligned} x^{'} & = - 7 x + 26 y \\ y^{'} & = - 4 x + 13 y \\ x (0) & = 5 \\ y (0) & = - 5 \end{aligned}

🔗

13.

\begin{aligned} x^{'} & = - 2 x + 13 y \\ y^{'} & = - 2 x + 8 y \\ x (0) & = - 2 \\ y (0) & = - 3 \end{aligned}

🔗

14.

\begin{aligned} x^{'} & = - 7 x + 13 y \\ y^{'} & = - 2 x + 3 y \\ x (0) & = 1 \\ y (0) & = - 1 \end{aligned}

🔗

15.

\begin{aligned} x^{'} & = 18 x - 52 y \\ y^{'} & = 8 x - 22 y \\ x (0) & = 5 \\ y (0) & = 3 \end{aligned}

🔗

16.

\begin{aligned} x^{'} & = - 12 x + 26 y \\ y^{'} & = - 4 x + 8 y \\ x (0) & = 2 \\ y (0) & = - 5 \end{aligned}

🔗

17.

Consider the linear system

d x / d t = A x,

where

🔗

A = (\begin{matrix} 3 & 2 \\ - 3 & - 1 \end{matrix}) .

Suppose the initial conditions for the solution curve are

x (0) = 1

and

y (0) = 1 .

We can use the following Sage code to plot the phase portrait of this system, including a solution curve.

🔗

xxxxxxxxxx
 
x, y, t = var('x y t') #declare the variables
F = [3*x + 2*y, -3*x - y] #declare the system
# normalize the vector fields so that all of the arrows are the same length
n = sqrt(F[0]^2 + F[1]^2)
# plot the vector field
p = plot_vector_field((F[0]/n, F[1]/n), (x, -30, 30), (y, -30, 30), aspect_ratio = 1)
# solve the system for the initial condition t = 0, x = 1, y = 1
P1 = desolve_system_rk4(F, [x, y], ics=[0, 1, 1], ivar = t, end_points = 5, step = 0.01)
# grab the x and y values
S1 = [ [j, k] for i, j, k in P1]
# plot the solution
# Setting xmin, xmax, ymin, ymax will clip the window
# Try plotting without doing this to see what happens
p += line(S1, thickness = 2, axes_labels=['$x(t)$','$y(t)$'], xmin = -30, xmax = 30, ymin = -30, ymax = 30) 
p

Use Sage to graph the direction field for the system linear systems

d x / d t = A x

in Exercise Group 3.4.6.9–16. Plot the solution curve for the given initial condition.

🔗

xxxxxxxxxx
 
​

🔗

You have attempted 1 of 4 activities on this page.

🔗

Prev Top Next