An autonomous differential equation has no explicit dependence on the independent time variable: its vector field depends only on the current state. The two families of nullclines are
On the first family the horizontal velocity is zero; on the second the vertical velocity is zero. A critical point of a dynamical system is a state where both components vanish. It gives a constant solution, also called an equilibrium.
Translate the equilibrium to the origin by . If the vector field is continuously differentiable, its first-order Taylor expansion is
The Jacobian matrix gives the linearization , whose general solution is using the matrix exponential. Its eigenvalues determine the exponential factors, while its eigenvectors and, when necessary, generalized eigenvectors determine their directions.
For distinct eigenvalues and corresponding eigenvectors ,
For a real matrix with nonreal conjugate eigenvalues, choose conjugate coefficients to obtain real solutions. Equivalently, if and , a real pair of independent solutions is
For a repeated eigenvalue , there are two cases:
  • If two independent eigenvectors exist, and .
  • Otherwise choose a Jordan chain , . Then .
The factor in the second case is essential. These are exact solutions of the linearization; the remainder describes the error for the nonlinear differential equation. When no eigenvalue has zero real part, the hyperbolic equilibrium theory justifies the corresponding local qualitative classification. Zero real parts require additional nonlinear analysis.
Factor each component of the vector field. The horizontal-velocity nullclines and vertical-velocity nullclines are, respectively,
The last curve is a hyperbola. In the nonnegative quadrant its branch is for .
Intersect the two nullcline families. The axes yield and , and with yields . The possible remaining intersection would have and , requiring , whose discriminant is . Hence there are no further real equilibria.
The Jacobian matrix is
Evaluate it at each equilibrium:
All three are hyperbolic equilibria, so their classifications persist for the nonlinear system. Immediately to the right of on , the velocity is vertically downward, showing that the attracting spiral turns clockwise.
For , the signs needed for the phase portrait follow directly from the factorizations:
The coordinate axes are invariant sets. On , the logistic differential equation directs positive states toward the boundary saddle point . On , directs positive states upward. At the three equilibria, use the local directions found in part (ii); interior sample trajectories near spiral clockwise inward.
Figure 1.
Nullclines and sample trajectories for the polynomial planar system
. Blue curves have zero horizontal velocity and orange curves have zero vertical velocity. Black arrows show the vector field and green curves show forward trajectories. The origin is an unstable node, the boundary point is a saddle, and the interior point is a clockwise spiral sink. Details.
The nullclines prescribe where arrows become vertical or horizontal; they are generally not themselves trajectories. The image samples the vector field and several solutions, rather than assuming every nullcline is an invariant set.

Articles by others on the same topic (0)

There are currently no matching articles.