ODE-ProjectODE-Project Bifurcations

Section 1.7 Bifurcations

Objectives

To understand that a one-parameter family of differential equations

$\frac{d x}{d t} = f_{λ} (x)$

has a bifurcation at $λ = λ_{0}$ if a change in the number or type of equilibrium solutions occurs.
To understand that bifurcation diagrams are an effective way of representing the nature of the solutions of a one-parameter family of differential equations.

Many of the equations that we have examined have a parameter, which means that we actually have a family of differential equations. For example,

🔗

\frac{d x}{d t} = k x

has the growth rate parameter

k .

The logistic equation

🔗

\frac{d P}{d t} = k P (1 - \frac{P}{N})

has two parameters,

k

and

N .

In this section we will investigate how the solutions of a differential equation vary as we change the value of a parameter.

🔗

Subsection 1.7.1 The Logistic Model with Harvesting Revisited

Recall how we modeled logistic growth in a trout pond in Example 1.3.9 with the equation

🔗

\frac{d P}{d t} = P (1 - \frac{P}{200}) .

If we allowed fishing in our pond at a rate of 32 fish per year, then the equation became

🔗

\frac{d P}{d t} = P (1 - \frac{P}{200}) - 32.

There are two equilibrium solutions for this equation,

P_{1} = 160

(a sink) and

P_{2} = 40

(a source). If the population of the pond falls below 40, then the fish will die out unless the pond is restocked or fishing is banned (Figure 1.7.1).

🔗

Figure 1.7.1. Harvesting with $H = 32$
🔗

Now let us see what happens when we allow more fishing in our pond, say

H = 100 .

Our differential equation now becomes

🔗

\frac{d P}{d t} = P (1 - \frac{P}{200}) - 100.

To determine the equilibrium solutions, we must solve

🔗

\begin{matrix} (1.7.1) & \frac{d P}{d t} = P (1 - \frac{P}{200}) - 100 = 0 \end{matrix}

for

P .

This last equation can be rewritten as

P^{2} - 200 P + 20000 = 0 .

Thus,

🔗

P = \frac{200 \pm \sqrt{200^{2} - 4 \cdot 20000}}{2} = 100 \pm \sqrt{- 10000},

which means that equation (1.7.1) has no real solutions and that we have no equilibrium solutions. Furthermore,

d P / d t < 0

for all values of

P .

This means that no matter how many fish are in the pond initially, the trout population will eventually die out due to overfishing (Figure 1.7.2).

🔗

Figure 1.7.2. Harvesting with $H = 100$
🔗

Finally, we will let

H = 50 .

In this case, we must solve

🔗

\frac{d P}{d t} = P (1 - \frac{P}{200}) - 50 = 0

in order to determine any equilibrium solutions. We now obtain a single equilibrium solution at

P = 100 .

In fact,

P = 100

will be a node. For values of

P > 100

as well as values of

P < 100,

we have

d P / d t < 0,

and the number of fish in the pond will decrease (Figure 1.7.3).

🔗

Figure 1.7.3. Harvesting with $H = 50$
🔗

To better understand what is happening, we will generalize our model. Suppose that a population with a limited carrying capacity

N

is modeled with the logistic equation

🔗

\frac{d P}{d t} = k P (1 - \frac{P}{N}) .

If we allow harvesting at a constant rate

H,

our model now becomes

🔗

\frac{d P}{d t} = k P (1 - \frac{P}{N}) - H .

To analyze our model, we will first find the equilibrium solutions. If we will let

🔗

f_{H} (P) = k P (1 - \frac{P}{N}) - H,

each equilibrium solution must satisfy

f_{H} (P) = 0

🔗

- k P^{2} + k N P - H N = 0.

Therefore, our equilibrium solutions are given by

🔗

P = \frac{- k N \pm \sqrt{k^{2} N^{2} - 4 k H N}}{- 2 k} = \frac{N}{2} \pm \sqrt{\frac{N^{2}}{4} - \frac{H N}{k}} .

🔗

The explanation of how our model behaves lies in the discriminant,

🔗

\frac{N^{2}}{4} - \frac{H N}{k} .

🔗

\frac{N^{2}}{4} - \frac{H N}{k} < 0

or, equivalently if

H > k N / 4,

there are no equilibrium solutions and

🔗

\frac{d P}{d t} = f_{H} (P) < 0

for all values of

P .

In particular, all solutions of

d P / d t = f_{H} (P)

tend towards negative infinity as

t \to \infty .

In this case, the population is doomed to extinction no matter how large the initial population is. Since negative populations do not make sense, we say that the population is extinct when

P = 0 .

🔗

On the other hand, if

H < k N / 4,

we have equilibrium solutions at

🔗

P_{1} = \frac{N}{2} + \sqrt{\frac{N^{2}}{4} - \frac{H N}{k}}

and

🔗

P_{2} = \frac{N}{2} - \sqrt{\frac{N^{2}}{4} - \frac{H N}{k}} .

The first equilibrium solution,

P_{1}

is a sink, while the second,

P_{2}

is a source.

🔗

Finally, if

H = k N / 4,

then we will have exactly one equilibrium solution at

P = N / 2 .

Although

d P / d t < 0

for all

P \neq N / 2,

we see that

P \to N / 2

t \to \infty

for all initial values of

P

greater than

N / 2 .

For initial values of

P

less than

N / 2,

solutions tend towards

- \infty

t \to \infty .

Thus, the initial population of fish must be at least

k N / 4;

otherwise, the fish will go extinct.

🔗

In our example, we have a family of differential equations—one for each value of

H,

🔗

\begin{matrix} (1.7.2) & \frac{d P}{d t} = P (1 - \frac{P}{200}) - H . \end{matrix}

A small change in

H

can have a dramatic effect on how the solutions of the differential equation behave. Changing the value of

H

from

50

50.1

will doom the population of fish to extinction no matter what the initial population is. As we increase the value of

H,

the number of equilibrium solutions changes from two to one and then to none. This change occurs exactly at

H = 50 .

We say that a bifurcation occurs at

H = 50

for equation (1.7.2).

🔗

Activity 1.7.1. Upland Bird Populations.

The chukar partridge, or simply chukar, is a upland gamebird in the pheasant family. Originally native to Asia and ranging from the eastern Mediterranean to Himalayas, the chukar has been widely introduced as an upland game bird with populations now established in the United States, Canada, Chile, Argentina, New Zealand and Hawaii. One particularly good area for hunting chukar is the western Great Basin area of the U.S. (eastern Oregon and Washington and western Idaho).

🔗

(a)

Suppose that the population of chukar on a private game ranch in eastern Oregon grows logisitically. Estimates tell us that the one hundred square mile ranch and that the ranch can support at most

120

birds per square mile. The growth rate of the chukar population is estimated to be

1.5

birds per year. Model the growth of the chukar population with an initial value problem.

🔗

(b)

Suppose that hunting on the ranch is restricted to guests and the average guest harvest

10

chukars per visit. Modify the model in part (a) to take into account the effect that hunting has on the chukar population.

🔗

(c)

What is the maximum number of guests that the ranch can accommodate and still maintain a healthy population of game birds? How many chukar per square mile would be needed to allow this many guests?

🔗

Subsection 1.7.2 One-Parameter Families

Let us consider the equation

🔗

\begin{matrix} (1.7.3) & \frac{d x}{d t} = x^{2} - 4 x + λ \end{matrix}

as a family of differential equations indexed by the parameter

λ .

If we let

f_{λ} (x) = x^{2} - 4 x + λ,

then

🔗

\frac{d x}{d t} = f_{λ} (x)

is a called one-parameter family of differential equations. For each value of

λ,

we obtain an autonomous differential equation, and for each value of

λ,

we have a different phase line to examine.

🔗

For

λ = 0,

the differential equation

🔗

\frac{d x}{d t} = f_{0} (x) = x^{2} - 4 x = x (x - 4),

there is a sink at

x = 0

and a source at

x = 4

(Figure 1.7.4).

🔗

For

λ = 4,

the differential equation

🔗

\frac{d x}{d t} = f_{4} (x) = x^{2} - 4 x + 4 = (x - 2)^{2},

we have exactly one equilibrium solution, a node at

x = 2

(Figure 1.7.5).

🔗

λ = 8,

then the differential equation

🔗

\frac{d x}{d t} = f_{8} (x) = x^{2} - 4 x + 8

has no equilibrium solutions (Figure 1.7.6).

🔗

Figure 1.7.6. $x^{'} = x^{2} - 4 x + λ$ for $λ = 8$
🔗

In fact, the number of equilibrium solutions for (1.7.3) changes at

λ = 4 .

We say that

λ = 4

is a bifurcation value for the differential equation

🔗

\begin{matrix} (1.7.4) & \frac{d x}{d t} = f_{λ} (x) = x^{2} - 4 x + λ . \end{matrix}

For

λ < 4,

we have two equilibrium solutions.

🔗

x = 2 \pm \sqrt{4 - λ} .

For values of

λ > 4,

there are no equilibrium solutions. We can record all of the information for the various values in a graph called the bifurcation diagram. The horizontal axis is

λ

and the vertical axis is

x .

Over each value of

λ,

we will plot the corresponding phase line. The curve in the graph represents the various equilibrium solutions for the different values of

λ .

The bifurcation diagram for equation (1.7.4) is a parabola (Figure 1.7.7). We have a phase line for each value of

λ .

🔗

Figure 1.7.7. Bifurcation diagram for $x^{'} = x^{2} - 4 x + λ$
🔗

Bifurcations for a one-parameter family of differential equations

d x / d t = f_{λ} (x)

are, in fact, rare. Let us consider a bifurcation where a sink changes to a source as we vary the parameter

λ .

Suppose that for

λ = λ_{0},

we have a sink at

x_{0} .

Then

🔗

\frac{d x}{d t} = f_{λ_{0}} (x_{0}) = 0.

Furthermore, the graph of

f_{λ_{0}} (x)

must be decreasing for

x

near

x_{0},

since

f_{λ_{0}} (x)

must be postive for values of

x < x_{0}

and negative for values of

x > x_{0} .

In other words,

f_{λ_{0}}^{'} (x) < 0

for

x

near

x_{0}

with

f_{λ_{0}}^{'} (x_{0}) < 0,

then for all

λ_{1}

sufficiently close to

λ_{0},

the differential equation

🔗

\frac{d x}{d t} = f_{λ_{1}} (x)

must have sink at a point

x = x_{1}

very close to

x_{0} .

A similar situation holds if

x_{0}

is a source and

f_{λ_{0}}^{'} (x_{0}) > 0 .

Thus, bifurcations can only occur when

f_{λ_{0}} (x_{0}) = 0

and

f_{λ_{0}}^{'} (x_{0}) = 0 .

🔗

Example 1.7.8.

Now consider the one-parameter family

\frac{d y}{d t} = y^{3} - α y = y (y^{2} - α) .

We will have an equilibrium solution at zero for all values of

α

and two additional equilibrium solutions at

\pm \sqrt{α}

for

α > 0 .

This type of bifurcation is a pitch fork bifurcation (Figure 1.7.9).

Figure 1.7.9. The bifurcation diagram for $y^{'} = y^{3} - α y$

Activity 1.7.2. Bifurcations.

For each of the following parametrized family of differential equations, plot phase lines for

c = - 2, - 1, 0, 1, 2,

find any bifurcation values, and sketch the bifurcation diagram.

🔗

(a)

y^{'} = (1 - y) y + c

🔗

(b)

y^{'} = (c - y^{2}) y

🔗

Example 1.7.10.

Let us find the bifurcation values of the one-parameter family

\begin{matrix} (1.7.5) & \frac{d y}{d t} = y (y - 2)^{2} + λ . \end{matrix}

g_{λ} (y) = y (y - 2)^{2} + λ,

then

g_{λ}^{'} (y) = 3 y^{2} - 8 y + 4 .

The roots of

g_{λ}^{'} (y) = 0

are

y = 2

and

y = 2 / 3 .

In order for

λ

to be a bifurcation value, we must have

g_{λ} (2) = λ = 0

g_{λ} (2 / 3) = \frac{32}{27} + λ = 0

Thus, equation (1.7.5) has two bifurcation values,

λ = - 32 / 27

and

λ = 0 .

The bifurcation diagram for this one-parameter family is given in Figure 1.7.11.

Figure 1.7.11. The bifurcation diagram for $y^{'} = y (1 - y)^{2} + λ$

🔗

Subsection 1.7.3 Important Lessons

A one-parameter family of differential equations

$\frac{d x}{d t} = f_{λ} (x)$

has a bifurcation at $λ = λ_{0}$ if a change in the number of equilibrium solutions occurs.
Bifurcation diagrams are an effective way of representing the nature of the solutions of a one-parameter family of differential equations.
Bifurcations for a one-parameter family of differential equations $d x / d t = f_{λ} (x)$ are rare. Bifurcations occur when $f_{λ_{0}} (x_{0}) = 0$ and $f_{λ_{0}}^{'} (x_{0}) = 0 .$