Math 536 - Mathematical Modeling Fall Semester, 2000 © 1999, All Rights Reserved, SDSU & Joseph M. Mahaffy San Diego State University -- This page last updated 20-Sep-00

Discrete Logistic Growth Model

In the previous section, we studied the Discrete Malthusian growth model, which showed exponential growth. This model is appropriate for early phases of population growth for most animal populations. However, as a population grows, it encounters crowding pressure due to many factors such as toxic build up or space and resource limitation.

The Malthusian growth model is a linear model, so the next obvious addition to the updating function is the addition of a quadratic term, which should be negative to reflect a decrease in the growth of a population due to crowding effects. This is the Logistic Growth model and can be written: This equation is the Malthusian growth model with the additional term -rPn2/M. The parameter M is called the carrying capacity of the population.

The behavior of the Logistic growth model is substantially more complicated than that of the Malthusian growth model. There is no exact solution to this discrete dynamical system. The ecologist Robert May (1974) studied this equation for populations and discovered that it could produce very complicated dynamics. In its simplest form the Logistic growth model can be written: where the parameter m varies between 0 and 4. For a good description of this model complete with cobwebbing and Java applet simulations see the website of B. Fraser.

Yeast Model

In 1913, Carlson  studied a growing culture of yeast. Below is a table of the population for these yeast at one hour intervals. We would like to develop a mathematical model to describe the growth of this culture.

 Time Population Time Population Time Population 1 9.6 7 174.6 13 594.8 2 18.3 8 257.3 14 629.4 3 29.0 9 350.7 15 640.8 4 47.2 10 441.0 16 651.1 5 71.1 11 513.3 17 655.9 6 119.1 12 559.7 18 659.6

The general discrete dynamical system is given by the

pn+1 = f(pn).

This is an iterative map where the population at the (n + 1)st generation depends on the population at the nth generation. The function f(p) is called the updating function as it produces the next population in the iterative scheme. A graph of the updating function has the (n + 1)st generation on the vertical axis, while the nth generation is on the horizontal axis.

It is clear that the population of yeast in the table above does not satisfy a Malthusian growth, which has a linear updating function and grows exponentially without bound. The next obvious addition to the updating function is the addition of a quadratic term, which should be negative to reflect a decrease in the growth of a population due to crowding effects. This is the Logistic Growth model and can be written: This equation has the Malthusian growth model seen in the previous section with the additional term -rPn2/M. The parameter M is called the carrying capacity of the population.

The behavior of the Logistic growth model is substantially more complicated than that of the Malthusian growth model. There is no exact solution to this discrete dynamical system. In its simplest form the Logistic growth model can be written: where the parameter m varies between 0 and 4. For a good description of this model complete with Java applet simulations see the website of B. Fraser.

The graph below shows a plot of pn+1 vs. pn from the table above. (This is accomplished by plotting the population from one time against the population from the previous time. For example, the first two points are (9.6, 18.3) and (18.3, 29.0).) The graph of the data is fit applying Excel's polynomial fit with trendline, using a quadratic passing through the origin, and is shown below. This equation becomes the updating function. Our discrete logistic growth model for the yeast experiment above is given by

pn+1 = 1.56pn - 0.000861pn2.

Below is a simulation showing both the data and the model. As we can see, the model does a fairly reasonable job of simulating the data from this fairly simplistic model. Qualitatively, we see the same initial roughly exponential growth, then both models seem to level off at approximately the same value. This is the carrying capacity of the population. The equation above shows r = 0.56 and r/M = 0.000861, so M = 650.4. This is clearly a little low based on the original table.

Cobwebbing

The website of B. Fraser shows a process of geometrically viewing the updating function, which is called cobwebbing. The updating function, f(pn), is graphed with the identity map, pn+1 = pn on a single graph with the vertical axis being pn+1 and the horizontal axis being pn. The idea is that you start at any p0, then go vertically to p1 = f(p0). Next you go horizontally to the identity map to locate p1 on the horizontal axis. From there you find p2 by going vertically to p2 = f(p1). This process is repeated to generate the cobweb of points by this discrete dynamical system. The sequence of points on the horizontal axis form the solution set generated by the discrete dynamical system. The graphical representation allows you make some projections of the behavior of the system. Below is a diagram showing several steps in the cobwebbing scheme for the quadratic map

pn+1 = 3pn(1 - pn ). U. S. Census with Logistic Growth Model

How well does this model work in actual applications? We return to the example using the U. S. census data in the previous section. We saw that the Malthusian growth model was very limited in its ability to project the population into the future with the exponential rate of growth causing the model to eventually over predict the population.

When you learned about cobwebbing, you found that one way to graphically represent the updating function was to plot Pn+1 versus Pn. This is easily done with the U. S. census data. So it remains a question of how to best fit an updating function of the quadratic form given by the Logistic growth model through this data. It turns out that if you simply perform the least squares best fit of a quadratic function (without the constant term) through this data, it gives too much weight to the later years of the model and grossly underestimates the rate of growth r. Based on the assumptions of the Logistic growth model that increasing population should decrease the growth rate, I chose r = 0.351 (matching the initial growth rate between 1790 and 1800), then performed a least squares best fit to the data. This gave the Logistic growth model: The graph of this updating function with the data is seen below: The simulation of this discrete Logistic growth model compared to the two models of the previous section and with the data is seen below: It is clear that the discrete Logistic growth model follows the data better than the Malthusian growth model (even the Malthusian model with r = 0.349), but not as well as the non-autonomous growth model of the previous section. The Logistic growth model shown above predicts the population relatively well up until 1900. However, it even has the wrong concavity after 1930, suggesting a premature leveling off of the population of the U. S. Thus, a time varying growth rate seems the best of these models. This shows that human populations have a more complicated dynamics than these models can predict with both time-varying and crowding factors. (Sociological and technical factors are especially difficult to incorporate into mathematical models.) Your ecology courses should help explain more details underlying the assumptions in these models, so explain a little better when the models are applicable and why they fail in other predictions.

Equilibria

Later we would like to study the qualitative behavior of discrete dynamical equations in more detail. Consider the general discrete dynamical equation:

Pn+1 = f(Pn).

The first step in any analysis is finding equilibria, which is simply an algebraic equation (and the area that we are reviewing now). An equilibrium point of a discrete dynamical system is a point where there is no change in the variable from one iteration to the next. Mathematically, this occurs whenever there is a solution to

Pe = f(Pe).

Graphically, this is when f(P) crosses the line y = P, which is one reason why this line was shown above.

Consider the original discrete Logistic equation listed above with r > 0. The equilibria are found by solving: Thus, the equilibria for the Logistic growth model are either the trivial solution 0 (no population) or the carrying capacity M.

If we examine the Logistic growth model proposed above for the U. S. census data, then we find the equilibria satisfying: Thus, Pe = 0 or 2.849*108. Note that this later population will be reached in less than 20 years, while actual projections have the U. S. population rising to over 400 million by the middle of the next century. It is unlikely that human population can continue its current course, but what will be the actual scenario? Mathematical modeling can provide reasonable estimates for short term growth and allows one to predict several different possibilities for longer term depending on the assumptions that are entered into the model.

There are Worked Examples available to help you understand the material.

Other Behavior of the Logistic Growth Model

Below we present an applet that shows simulations of the logistic growth model for various choices of r.

Robert May (1974) demonstrated that the discrete Logistic growth model could display very complicated dynamics.Watch what happens in the applet above as we choose different values of r. For example, try the values r = 0.5, 1.8, 2.3, and 2.65. (Note that the solution of the discrete Logistic equation only gives solutions at the integer values of n, so the connecting lines are only drawn to help visualize the behavior of the system.) The first value of r shows the curve smoothly ascending to carrying capacity of 1000. The second value of r has the population ascend and actually overshoot 1000, then oscillates about 1000, getting closer to the carrying capacity as n increases. In both of these cases, the equilibrium population of 1000 is said to be stable. When r = 2.3, the solution oscillates about 1000, taking on the values of approximately 690 and 1180. This solution is said to have a period of 2. The last case shows the population oscillating almost randomly about 1000. This last situation could have either a very high period of oscillation or actually be chaotic. Later in this course we will study more of the mathematics behind this sequence of changes, known as bifurcations. We will find that derivatives play an important role in the study of the behavior of these systems.

Stability of the Logistic Growth Model

When studying the growth of a culture of yeast, we found that the discrete logistic growth model simulated the data quite well. However, the applet for the discrete logistic growth model showed that the behavior of this model changes very dramatically as the parameter r varies. The derivative of the function f(P) is a valuable tool for determining the behavior of the discrete dynamical system near an equilibrium point.

The discrete logistic growth model, as given by the equation above, has two equilibria,

Pe = 0 and Pe = M.

It is easy to see that the derivative of f(P) is given by

f '(P) = 1 + r - 2 rP/M.

At Pe = 0, the derivative satisfies

f '(0) = 1 + r,

which for r positive always results in solutions growing away from this equilibrium. (If r is a negative number between -1 and 0, then the solution of the discrete logistic growth model decays to 0 or the population goes to extinction.)

The more interesting behavior occurs around the second equilibrium, Pe = M. Below is a summary of the types of behavior that are observed for a discrete dynamical system near an equilibrium, Pe. (In all of the descriptions below, we are assuming that the simulation begins "near" the equilibrium value.

Behavior of the Discrete Dynamical Model near an Equilibrium

1. If f '(Pe) > 1, then the solutions of the discrete dynamical model grow away from the equilibrium (monotonically). Thus, the equilibrium is unstable.
2. If 0 < f '(Pe) < 1, then the solutions of the discrete dynamical model approach the equilibrium (monotonically). Thus, the equilibrium is stable.
3. If -1 < f '(Pe) < 0, then the solutions of the discrete dynamical model oscillate about the equilibrium and again the equilibrium is stable.
4. If f '(Pe) < -1, then the solutions of the discrete dynamical model oscillate but move away from the equilibrium. Again, the equilibrium is unstable.

Returning to the logistic growth model, we can evaluate the derivative at the larger equilibrium, Pe = M. From the formula for the derivative, it is easy to see that

f '(M) = 1 - r.

From our list of behaviors above, it is follows that

1. If 0 < r < 1, then the solution of the discrete logistic model monotonically approaches the equilibrium, Pe = M, which was the case observed for the experiment with the yeast.
2. If 1 < r < 2, then the solution of the discrete logistic model oscillates about the equilibrium, Pe = M, but the solution still asymptotically approaches this equilibrium.
3. If 2 < r < 3, then the solution of the discrete logistic model oscillates about the equilibrium, Pe = M, but the solution grows away from this equilibrium. ( r > 3 results in negative solutions.

References

 T. Carlson Über Geschwindigkeit und Grösse der Hefevermehrung in Würze. Biochem. Z. 57: 313-334, 1913.

 Statistical Abstracts of the United States (1993) 113 th ed., U. S. Department of Commerce, Bureau of the Census, Washington, DC.