Section 5.1 Linearization
Objectives

To understand that a nonlinear system
\begin{align*} x'(t) \amp = f(x,y)\\ y'(t) \amp = g(x, y) \end{align*}can be approximated near each equilibrium point \((x_0, y_0)\) by a linear system by using a Taylor series approximation for \(f\) and \(g\text{.}\)

To understand that the Jacobian matrix,
\begin{equation*} J = \begin{pmatrix} f_x(x_0, y_0) \amp f_y(x_0, y_0) \\ g_x(x_0, y_0) \amp g_y(x_0, y_0) \end{pmatrix}, \end{equation*}of a nonlinear system
\begin{align*} x'(t) \amp = f(x,y)\\ y'(t) \amp = g(x, y) \end{align*}is useful for classifying equilibrium solutions of the system as stable or unstable.
To understand that linearization can fail for some systems.
In ChapterÂ 3, we explored how to solve linear systems. However, if we are asked to solve a nonlinear system such as
we are faced with a much more difficult task. In general, it may not be possible to find solutions for a nonlinear system in terms of elementary functions. However, for a given modeling problem, we can ask many questions that may be answered without finding an explicit solution for the associated system of differential equations. For example, what are the equilibrium points? Are the equilibrium points stable? Do closed solution curves exist in the phase plane?
Subsection 5.1.1 Equilibrium Solutions
An equilibrium solution for a firstorder system of differential equations
is a point \((x_0, y_0)\) such that
Notice that neither \(x\) nor \(y\) is changing at this point. If we have the initial conditions, \(x(0) = x_0\) and \(y(0) = y_0\text{,}\) then the solution to the system is \(x(t) = x_0\) and \(y(t) = y_0\text{.}\) Of course the interesting question is what happens if our initial conditions are close to an equilibrium solution. Do solutions tend toward the equilibrium solution, away from the equilibrium solution, or is there a combination of the two?
One of the most useful methods of determining the nature of an equilibrium solution for a given nonlinear system is to approximate the nonlinear system with a linear system. More specifically, an equilibrium solution occurs for the linear system
when an \(x\)nullcline intersects a \(y\)nullcline. That is, when a curve defined by \(dx/dt = f(x,y) = 0\) intersects a curve defined \(dy/dt = g(x,y) = 0\text{,}\) we have an equilibrium solution. Since the \(x\) and \(y\)nullclines are simply curves in the \(xy\)plane, we can approximate them locally by intersecting straight lines. Translating the pair of intersection lines to the origin, we obtain a linear system, and we can apply everything that we learned about such systems in ChapterÂ 3.
Example 5.1.1.
Consider the system
From
we can quickly conclude that the only equilibrium solution to the system is \((0, 0)\text{.}\) The phase portrait for this system is given in FigureÂ 5.1.2. If we have the initial conditions \(x(0) = 1\) and \(y(0) = 1\text{,}\) we can see that the solution tends toward the origin as \(t \to \infty\text{.}\) However, it is unclear from the phase portrait if the solution curves of all initial value problems with initial conditions near the origin tend towards the equilibrium solution as \(t \to \infty\text{.}\)
Since the nonlinear terms of the system do not contribute much towards \(dx/dt\) and \(dy/dt\) for values of \(x\) and \(y\) near zero, we can determine the nature of the equilibrium solution by examining the system consisting of only linear terms on the righthand side of the equation,
The matrix for this linear system,
has eigenvalues \(\lambda = 1\) and \(\mu = 2\) with eigenvectors \(\mathbf u = (3,2)\) and \(\mathbf v = (1,1)\text{,}\) respectively. The general solution for this linear system is
This indeed suggests that solutions near the origin tend towards the origin as \(t \to \infty\text{.}\) In this case, we say that the equilibrium solution is stable. Of course, if we are given an initial condition such as \(x(0) = 0.5\) and \(y(0) = 1\text{.}\)
Example 5.1.3. A Competing Species Model.
Suppose that \(x\) and \(y\) are the population of two distinct species that compete for the same resources. For example, two species of fish may compete for the same food in a lake or sheep and cattle competing for the same grazing land. Recall from SectionÂ 2.2 that we can model two competing species using the following system of firstorder differential equations,
The first term in each equation is the logistic growth of each species. The second term in each equation tells how each species is affected by interacting with the competing species.
More specifically, let us consider the following system,
It is easy to see that the four equilibrium solutions: \((0,0)\text{,}\) \((0,3)\text{,}\) \((2,0)\text{,}\) and \((1, 1)\text{.}\) We can view the direction field, the phase plane, and some solution curves for this system in FigureÂ 5.1.4.
Let us analyze what happens at the equilibrium solution \((1, 1)\text{.}\) If we decide on an appropriate change of variables, we can translate the entire system so that this equilibrium solution is at the origin. If we let
then
EquationsÂ (5.1.1) and (5.1.2) now become
As before we consider only the linear part of equationsÂ (5.1.4) and (5.1.4) are
The idea is that the linear part is a good local approximation to the original equations much like a tangent line is a good local approximation to a smooth function in calculus. We can determine the local nature of the equilibrium solution by examining the eigenvalues of the matrix
The eigenvalues of \(A\) are \(\lambda = 1 + \sqrt{2} = 0.4142\) and \(\mu = 1  \sqrt{2} \approx 2.4142\) with eigenvectors \(\mathbf u = (1,\sqrt{2}\,) \approx (1, 1.4142)\) and \(\mathbf v = (1,\sqrt{2}\,) \approx (1, 1.4142)\text{,}\) respectively. The solution to the linear system is now
or
If we have initial conditions \(x(0) = 3\) and \(y(0) = 4\text{,}\) one can quickly determine that \(c_1 = 3/2 + \sqrt{2} \approx 2.9142\) and \(c_2 = 3/2  \sqrt{2} \approx 0.0859\text{,}\) and equationsÂ (5.1.7) and (5.1.8) become
As \(t \to \infty\) notice that both \(x(t)\) and \(y(t)\) become very large and tend away from the origin.
We can conclude that the equilibrium solution \((1,1)\) is unstable. That is, tend away from the equilibrium solution as \(t \to \infty\) if one population begins with a slight advantage over the other. If neither popluation has an initial advantage over the other, then the solution curve will approach the equilibrium solution as \(t \to \infty\text{.}\)
Example 5.1.5. A Nonpolynomial Example.
If systems such as the following can be approximated by linear systems,
Certainly, this sytems has an equilibrium solution at \((0,0)\) (FigureÂ 5.1.6). We can expand \(\sin x\) into a power series,
Thus, this system can be approximated by the linear system
See FigureÂ 5.1.7.
Example 5.1.8.
For example, consider the system
The \(x\) and \(y\)nullclines of this system are the curves \(y = x^2  4x + 4\) and \(y = x^3\text{,}\) respectively. Since the two nullclines intersect only at \((1, 1)\text{,}\) we have a single equilibrium solution. From the phase plane, it appears that \((1, 1)\) is a stable equilibrium solution. That is, all solution curves starting sufficiently close to \((1, 1)\) will approach the equilibrium solution as \(t \to \infty\) (FigureÂ 5.1.9).
Making the substitution \(u = x  1\) and \(v = y  1\text{,}\) simply translates the entire system to the origin, and we obtain the system
Our new system has \(u\) and \(v\)nullclines \(v = u(u  2)\) and \(v = u(u^2 + 3u + 3)\text{,}\) respectively. Notice that we have simply moved the phase portrait of the original system so that our equilibrium solution is now at the origin. Furthermore, we can approximate the \(u\) and \(v\)nullclines by their tangent lines \(v = 2u\) and \(v = 3u\text{,}\) respectively. From FigureÂ 5.1.9, it appears that we are approximating our original system with a linear system.
To determine the general case, we must approximate a nonlinear system
near each equilibrium point \((x_0, y_0)\) with a linear system. The idea is to move the equilibrium solution to the origin with a change of coordinates and then approximate the nonlinear system with a linear system. In order to move the equilibrium solution to the origin, we will introduce new variables
If \((x,y)\) is close to the equilibrium solution \((x_0, y_0)\text{,}\) then \((u, v)\) will be close to the origin. Under this change of coordinates,
Similarly, \(dv/dt = dy/dt = g(x_0 + u, y_0 + v)\text{.}\)
In order to approximate the nonlinear system with a linear system, we will use a Taylor series in two variables. That is, we can write
If we only use the linear terms of the Taylor series, we obtain a fairly accurate approximation of \(f\) provided we are near the equilibrium solution. Of course, \(f\) and its linear approximation may be quite different for values far away from the equilibrium solution. The linearization of the system of equationsÂ (5.1.15)â€“(5.1.16) now becomes
Since \((x_0, y_0)\) is an equilibrium solution, the constant terms vanish in each equation and
If we write our linearization in matrix form, we have
The matrix
is the Jacobian matrix of the system. We can now classify equilibrium solutions of nonlinear systems by examining the eigenvalues of the Jacobian matrix of the system or by using the tracedeterminant plane. For example, if \(D =\det(J) \gt 0\) and \(T =\trace(J) \lt 0\text{,}\) we have a sink. We can determine the sink to be spiral or nodal by examining whether \((T, D)\) lies above or below the parabola \(T^2 = 4D\) in the tracedeterminant plane.
It is important to note that linearization only tells us the local story. A solution curve might behave quite differently if it is far away from the equilibrium solution.
Example 5.1.11.
In ExampleÂ 5.1.8, we considered the system
The Jacobian matrix of this system is
For the equilibrium solution \((1, 1)\text{,}\)
Since \(J\) has eigenvalues \(\lambda = (3 \pm i \sqrt{11})/2\text{,}\) the equilibrium solution will act as a spiral sink for initial values close to \((1,1)\text{.}\)
Activity 5.1.1. Classifying Equilibrium Solutions for Nonlinear Systems.
Find all of the nullclines and equilbrium solutions for each of the following systems. Classify each equilibrium solution as stable or unstable. Plot the nullclines and a direction field for each system.^{â€‰1â€‰}
(a)
(b)
(c)
(d)
Subsection 5.1.2 When Linearization Fails
There are at least two cases when linearization does not give us the information that we seek. First, it might well be the case that the linear terms vanish in the nonlinear system. For example, the system
has an equilibrium solution at the origin, but this linearization of this system vanishes.
A more subtle example is the system
The origin is an equilibrium solution, and the linearization of this system is
The eigenvalue of the matrix corresponding to the linear system,
are \(\lambda = \pm i\text{.}\) Thus, the solution curves of the linear system are simply circles about the origin. However, the nonlinear system acts quite differently. In the nonlinear system, solutions spiral out slowly from the origin (FigureÂ 5.1.12). This system has no periodic solutions.
Subsection 5.1.3 Important Lessons
We can approximate a nonlinear system
\begin{align*} x'(t) \amp = f(x,y)\\ y'(t) \amp = g(x, y) \end{align*}near each equilibrium point \((x_0, y_0)\) with a linear system by using a Taylor series approximation for \(f\) and \(g\text{.}\) The matrix of our linear approximation,\begin{equation*} J = \begin{pmatrix} f_x(x_0, y_0) \amp f_y(x_0, y_0) \\ g_x(x_0, y_0) \amp g_y(x_0, y_0) \end{pmatrix}, \end{equation*}is the Jacobian matrix of the system. We can classify nonlinear systems by examining the Jacobian matrix of the system and using the tracedeterminant plane.Linearization only tells us how solutions behave near the equilibrium point. A solution curve might behave quite differently if it is far away from the equilibrium solution.
In some cases, linearization can fail.
Reading Questions 5.1.4 Reading Questions
1.
What does it mean to linearize a nonlinear system?
2.
How might the linearization of a nonlinear system fail?
Exercises 5.1.5 Exercises
Linearization.
For each of the systems in Exercise GroupÂ 5.1.5.1â€“6,
Plot and label the nullclines for equation in the system.
Find all of the equilibrium solutions for the system.
Use the Jacobian to classify each equilibrium solution (spiral source, nodal sink, etc.).
1.
2.
3.
4.
5.
6.
7.
Consider the following three systems
 \begin{align*} x' & = 3 \sin x + y\\ y' & = 4x + \cos y  1 \end{align*}
 \begin{align*} x' & = 3 \sin x + y\\ y' & = 4x + \cos y  1 \end{align*}
 \begin{align*} x' & = 3 \sin x + y\\ y' & = 4x + 3\cos y  3 \end{align*}
All three systems have an equilibrium solution at \((0, 0)\text{.}\) Which two systems have phase portraits with the same â€ślocal pictureâ€ť near \((0, 0)\text{?}\) Justify your answer.
8.
Let us consider an epidemic model for a city. We make the following additional assumptions about our model.
The city has a constant birth rate rate of \(\alpha\) persons per day. All new born babies are susceptible to the disease.
Infected people either recover or die after a certain number of days. If an individual recovers, he or she is immune.
If we let \(S(t)\) be the number of susceptible individuals at time \(t\) and \(I(t)\) be the number of infected individuals at time \(t\text{,}\) our assumptions give rise to the following system of differential equations,
The constant \(\alpha\) is determined by the probability of an interaction between a susceptible individual and an infected individual, and \(\gamma\) is the rate at which infected individuals are removed from the population. If
then both the susceptible and infected populations do not change. This will occur at
We are interested in the behavior of solutions near \((S_0, I_0)\text{.}\) If solutions approach this equilbrium point, then the disease will become endemic to the population.
9.
Consider the predatorprey system modeled by the following equations,
Find the equilibrium solutions of this system.
What are the eigenvalues of the Jacobian for each equilibrium solution?
What, if anything, can be said about the nature of the equilibrium solutions of the system?