Two Theorems on Stability

In the module on Cobweb Diagrams we looked at linear discrete dynamical systems -- that is, models of the form

pn = f(pn-1)

where the function f(p) on the right side of the equation is a linear function

f(p) = m p + b

We observed that if m is not one then there is a unique equilibrium point

missing equation

and that sometimes this equilibrium point is attracting -- that is, over the long term the terms pn approach p* -- and sometimes it is repelling -- that is, the terms go off to plus or minus infinity.

In that module you experimented with this phenomenon and were asked to find and prove a theorem that would determine when the equilibrium point was attracting and when it was repelling. We also looked at some nonlinear models. Now, in this module, we formulate and prove two theorems that determine the character of equilibrium points for linear and nonlinear models.

The Classification Theorem for Linear Discrete Dynamical Systems

Consider the linear dynamical system

missing equation

and if m is not one let

missing equation

denote the unique equilibrium point. Then

Proofs

We present three different proofs. This theorem is very accessible to high school students. The different proofs can be used at different times in the high school curriculum. The first proof uses geometric series.

A proof using geometric series

This proof applies only to the first clause of the conclusion. We begin by computing p2,

missing equation

then p3,

missing equation

and then p4.

missing equation

Now we can see a pattern emerging (the proof is a nice exercise in mathematical induction).

missing equation

Thus, we see that

missing equation

If |m| < 1 the first term on the right side of the equation above goes to zero as n goes to infinity and the last term on the right side of the equation above involves a geometric series whose sum is given by

missing equation

Thus, we see that

Missing equation

The next two proofs are both based on the same idea. We compare the distance

Missing equation

from the n-th term to the equilibrium point to the distance

Missing equation

from the preceding term to the equilibrium point. That is, we want to compute

Missing equation

If this fraction is less than one than the terms are getting closer to the equilibrium point and the equilibrium point is attracting -- that is,

missing equation

but, if the fraction is bigger than one the terms are going away from the equilibrium point and the equilibrium point is repelling.

A geometric proof

Looking at the cobweb diagram

Missing fihure

we see that

Missing equation

is just the absolute value of the slope of the line f(p) = m p + b, or m. Thus, if |m| < 1 the equilibrium point is attracting and if |m| > 1 th equilibrium point is repelling.

An algebraic proof

Our algebraic proof is based on the following computation

Missing equation

So, once again, we see that if |m| < 1 the terms of the sequence are approaching p* and if |m| > 1 they are getting further away.

Nonlinear Discrete Dynamical Systems

For nonlinear discrete dynamical systems -- that is, systems of the form

pn = f(pn-1)

where the function f(p) is nonlinear -- the situation is different. The first thing to notice is that a nonlinear discrete dynamical system can have many equlibrium points. The graph below shows one example of a nonlinear dynamical system. The beginning of a cobweb diagram is shown on the green graph paper and the beginning of a graph of the sequence is shown on the blue graph paper. This figure is a live Java applet and is very similar to the Java applets in the preceding two modules. You can change the initial condition by clicking along the bottom edge of the green graph and you can graph the resulting sequence and cobweb diagram one term at a time by clicking repeatedly in the upper two-thirds of the applet.

The phenomena that we observe in the example above can be used in two different ways in the high school curriculum.

The Classification Theorem for Nonlinear Discrete Dynamical Systems

Consider the dynamical system

pn = f(pn-1)

and let p* be an equilibrium point -- that is, suppose

f(p*) = p*.

Then, if |f'(p*)| < 1, the equilibrium point is attracting in the sense that if the initial condition p1 is sufficiently close to p* then

Missing equation

We will not prove this theorem here.

An equilibrium point that is attracting in the sense described above is sometimes called stable because if a system moves slightly away from the equilibrium because of some transient problem then the system will return to the equilibrium.


Consider the logistic models

pn = f(pn-1)

where

Missing equation

In the module on logistic models we experimented with these models and looked at how the constant a affected their long term behavior. Now, using the Classification Theorem for Nonlinear Discrete Dynamical Systems, we can learn more about these logistic models.

[Next section -- Age Dependent Population Models]


Copyright c 1997 by Frank Wattenberg, Department of Mathematics, Montana State University, Bozeman, MT 59717