Historical Context & Motivation
For centuries, scientists described change using single equations — one equation tracking one variable over time. But the real world is rarely that simple. A predator population affects its prey, and the prey population in turn affects the predators. Two chemical concentrations influence each other simultaneously. To capture this interconnected behavior, mathematicians developed systems of differential equations, where multiple quantities change together. The challenge then became: how do you see what such a system does without solving every equation by hand?
The answer emerged through a powerful geometric idea. Instead of plotting each variable separately against time, mathematicians plotted the variables against each other. The resulting picture — a phase portrait — reveals the long-term behavior of a system at a glance: does it spiral inward to a stable state, explode outward, or cycle endlessly?
The central question that phase portraits address is this: given a system of differential equations, can we predict where every possible starting condition leads — without solving the equations algebraically? As you will see, the answer is a resounding yes, and the tool for doing so is both elegant and surprisingly intuitive.
Core Principles & Definitions
Before diving into diagrams, you need a handful of key ideas. A system of two first-order differential equations tracks how two quantities — call them x and y — change over time. The phase plane is the xy-plane where each point represents a specific state of the system. As time passes, the state traces out a curve called a trajectory. The collection of many trajectories from different starting points forms the phase portrait.
Phase Plane
Trajectory (Solution Curve)
Equilibrium Point
Direction Field
Stability
Visual Explanation — Anatomy of a Phase Portrait
The diagram below shows a typical phase portrait for a stable spiral (also called a spiral sink). Notice how every trajectory curves inward toward the equilibrium at the origin. The arrows on each curve indicate the direction of motion as time increases. Trajectories that start far from the origin still eventually wind their way in, illustrating asymptotic stability: the system settles toward the equilibrium state in the long run.
Several important features stand out. First, no two trajectories ever cross each other — this is a consequence of the uniqueness theorem for differential equations. If two trajectories crossed at a point, that point would have two different "futures," which is impossible for a well-defined system. Second, the winding (spiraling) pattern tells you the system involves oscillatory behavior — the state rotates as it changes. Third, the inward drift tells you the oscillations are decaying — the system is losing energy over time and settling toward the equilibrium.
Mathematical Framework
The simplest systems studied through phase portraits are linear systems with constant coefficients. These involve two first-order equations where the rate of change of each variable depends on both variables in a straightforward (linear) way. The key numbers that determine the shape of every trajectory are the trace and determinant of the system's coefficient matrix — two summary numbers you can compute directly from the four constants in the equations.
Classifying Equilibrium Points
Two numbers — the trace T and the determinant D — are all you need to classify the shape of the phase portrait. The diagram below maps out the major types based on where (T, D) falls in the plane. This classification is the single most useful tool for reading a phase portrait.
| Type | Conditions (T, D, Δ) | Trajectory Shape | Stability |
|---|---|---|---|
| Stable Node | T < 0, D > 0, Δ > 0 | Straight-line paths converging to origin | Stable |
| Unstable Node | T > 0, D > 0, Δ > 0 | Straight-line paths diverging from origin | Unstable |
| Saddle Point | D < 0 | Hyperbolic curves — some approach, some flee | Unstable |
| Stable Spiral | T < 0, D > 0, Δ < 0 | Spiraling inward toward origin | Stable |
| Unstable Spiral | T > 0, D > 0, Δ < 0 | Spiraling outward from origin | Unstable |
| Center | T = 0, D > 0 | Closed elliptical orbits | Neutrally stable |
Worked Example
Let us classify the equilibrium and sketch the phase portrait for a concrete system. Consider:
Strengths & Limitations of Phase Portrait Analysis
Phase portraits are remarkably powerful, but like any tool, they have boundaries. Understanding both their strengths and their limitations will help you decide when to use this approach and when you need something more.
| Strengths | Limitations |
|---|---|
| Reveal global behavior at a glance without solving the system algebraically. | Limited to two-variable systems for easy visualization (3D phase space is harder to interpret). |
| Classify stability and type of equilibrium quickly using just the trace and determinant. | Do not tell you the exact solution x(t) and y(t) — only the qualitative shape of solutions. |
| Work for nonlinear systems too (by approximating the system with a linear one near each equilibrium). | The linear approximation can fail at borderline cases (for example, when the trace is exactly zero in a nonlinear system). |
| The uniqueness theorem guarantees trajectories never cross, simplifying the picture. | Time information is lost — you cannot read how fast the system moves along a trajectory from the portrait alone. |
Connection to Nonlinear Systems & Advanced Theory
Everything discussed so far applies to linear systems, but the real world is overwhelmingly nonlinear. The good news is that near an equilibrium point, a nonlinear system often behaves like a simpler linear approximation — this idea is called linearization. Just as you can approximate a curved road with a straight line when you zoom in close enough, you can approximate a nonlinear system with a linear one near each equilibrium. You then compute the trace and determinant of that local linear approximation to classify the behavior — exactly the same method as before.
| Feature | Linear Systems | Nonlinear Systems |
|---|---|---|
| Coefficient matrix | Constant matrix A — same everywhere in the plane | Local linear approximation — changes at each equilibrium |
| Number of equilibria | Only one (the origin), if D ≠ 0 | Can have multiple equilibria throughout the plane |
| Global behavior | Fully determined by trace and determinant | May include limit cycles, chaos, and other complex structures |
| Classification method | Trace–determinant analysis of A | Linearize at each equilibrium, then apply trace–determinant to the local approximation |
In more advanced courses, you will encounter phenomena that have no linear counterpart: limit cycles (isolated closed trajectories that nearby solutions approach or flee), bifurcations (qualitative changes in the portrait as a parameter varies), and even chaos in three or more dimensions. The trace–determinant classification you learned here is the essential first step toward all of these topics.
Practice Problems
Summary
A phase portrait is a collection of trajectories drawn in the phase plane — a coordinate system where each axis represents one of the system's variables, not time. For a linear system dx/dt = ax + by, dy/dt = cx + dy, the behavior of every trajectory is determined by just two numbers computed from the four constants: the trace T = a + d and the determinant D = ad − bc. The discriminant Δ = T² − 4D then tells you whether trajectories spiral (Δ < 0) or move along straight-line-like paths (Δ > 0).
The six main classifications are: stable nodes (T < 0, D > 0, Δ > 0 — paths converge along straight-like curves), unstable nodes (T > 0, D > 0, Δ > 0 — paths diverge), saddle points (D < 0 — always unstable), stable spirals (T < 0, Δ < 0 — spiral inward), unstable spirals (T > 0, Δ < 0 — spiral outward), and centers (T = 0, D > 0 — closed orbits, neutrally stable). Near equilibria of nonlinear systems, you can apply the same trace–determinant method to a local linear approximation to classify each equilibrium's behavior.