Why This Matters
Phase plane analysis is where linear algebra and differential equations come together in a powerful visual framework. You're being tested on your ability to connect eigenvalue analysis, matrix operations, and solution behavior to understand how dynamical systems evolve over time. The phase plane transforms abstract equations into geometric intuition—showing you why solutions spiral, converge, or explode rather than just that they do.
This topic bridges nearly everything you've learned: eigenvalues determine stability, eigenvectors set trajectory directions, and linearization lets you apply matrix techniques to nonlinear problems. When you see a phase portrait on an exam, you should immediately think about the underlying Jacobian, its eigenvalues, and what those values tell you about long-term behavior. Don't just memorize classifications—know what mathematical features produce each type of behavior.
Foundations: Setting Up the Phase Plane
Before analyzing behavior, you need to understand how we represent dynamical systems geometrically. The phase plane converts a system of two first-order ODEs into a two-dimensional picture where every point represents a complete state of the system.
Definition of a Phase Plane
- Two-dimensional state space—each axis represents one dependent variable (like x and y, or position and velocity)
- Every point is a complete system state, meaning if you know where you are in the phase plane, you know everything about the system at that instant
- Time is implicit—trajectories show the sequence of states but don't explicitly mark when each state occurs
Direction Fields and Trajectories
- Direction fields show the vector (x˙y˙) at each point, giving the instantaneous direction of motion
- Trajectories are solution curves that follow these vectors, representing the actual path a system takes from a given initial condition
- Uniqueness theorem guarantees trajectories never cross in autonomous systems—this constraint shapes all phase portrait analysis
Compare: Direction fields vs. trajectories—direction fields show possible motion at every point, while trajectories show actual paths from specific initial conditions. FRQs often ask you to sketch trajectories consistent with a given direction field.
Finding Critical Points: Equilibria and Nullclines
Equilibrium points are where the action stops—and where analysis begins. Identifying these points and understanding the flow around them is the first step in characterizing any system.
Equilibrium Points and Their Classification
- Equilibrium occurs where x˙=0 and y˙=0 simultaneously—the system has no tendency to move
- Classifications include nodes (stable/unstable), saddle points, spirals, and centers, each with distinct trajectory patterns
- Long-term behavior of most initial conditions is determined by which equilibrium they approach (or avoid)
Nullclines and Their Significance
- x-nullcline is the curve where x˙=0; y-nullcline is where y˙=0—trajectories cross these curves horizontally or vertically
- Intersections of nullclines are equilibrium points, making nullclines a powerful tool for locating equilibria graphically
- Nullclines divide the phase plane into regions where you can determine the sign of each derivative, revealing the general flow direction
Compare: x-nullclines vs. y-nullclines—on an x-nullcline, trajectories move purely vertically; on a y-nullcline, purely horizontally. If an FRQ gives you the equations, sketch both nullclines first to find equilibria quickly.
Stability Analysis: The Eigenvalue Connection
This is where linear algebra becomes essential. The eigenvalues of the Jacobian matrix evaluated at an equilibrium point completely determine local stability for hyperbolic equilibria.
Linearization of Nonlinear Systems
- Jacobian matrix J=(∂x∂f∂x∂g∂y∂f∂y∂g) captures the local linear approximation near an equilibrium
- Linearization is valid when eigenvalues have nonzero real parts (hyperbolic equilibria)—the linear system accurately predicts nonlinear behavior nearby
- Evaluate the Jacobian at the equilibrium point to get the coefficient matrix for stability analysis
Eigenvalues and Eigenvectors in Phase Plane Analysis
- Real parts of eigenvalues determine stability: both negative means stable, both positive means unstable, opposite signs means saddle
- Imaginary parts indicate oscillation: pure imaginary gives centers, complex with nonzero real part gives spirals
- Eigenvectors determine trajectory directions—for real eigenvalues, trajectories approach or leave along eigenvector directions
Stability Analysis of Equilibrium Points
- Asymptotically stable when all eigenvalues have negative real parts—trajectories converge to equilibrium as t→∞
- Unstable when at least one eigenvalue has positive real part—small perturbations grow over time
- Trace-determinant plane offers a quick classification: τ=tr(J), Δ=det(J), with stability requiring τ<0 and Δ>0
Compare: Stable nodes vs. stable spirals—both have eigenvalues with negative real parts, but nodes have real eigenvalues (direct approach) while spirals have complex eigenvalues (oscillatory approach). The discriminant τ2−4Δ tells you which.
Beyond Equilibria: Global Behavior
Local analysis tells you what happens near equilibrium points, but phase plane analysis also reveals global structures. Limit cycles and bifurcations describe behaviors that can't be understood from linearization alone.
Limit Cycles and Periodic Solutions
- Limit cycles are isolated closed trajectories—nearby solutions spiral toward (stable) or away from (unstable) the cycle
- Represent self-sustained oscillations that persist regardless of initial conditions, unlike the centers of linear systems
- Poincaré-Bendixson theorem guarantees limit cycles exist under certain conditions in planar systems with bounded trajectories and no equilibria
Bifurcations and Their Types
- Bifurcations mark qualitative changes in system behavior as a parameter varies—equilibria can appear, disappear, or change stability
- Saddle-node bifurcation: two equilibria collide and annihilate; Hopf bifurcation: equilibrium loses stability and spawns a limit cycle
- Bifurcation diagrams plot equilibrium locations/stability against parameter values, revealing critical thresholds
Compare: Limit cycles vs. centers—both show periodic motion, but limit cycles are isolated (nearby trajectories converge to them) while centers have families of closed orbits. Limit cycles are structurally stable; centers are not.
Phase Portraits: Putting It All Together
Phase portraits synthesize everything: equilibria, stability, nullclines, and global flow. Different system types produce characteristic portrait patterns that you should recognize instantly.
Phase Portraits for Different System Types
- Linear systems produce nodes, saddles, spirals, or centers depending entirely on eigenvalues—no limit cycles possible
- Nonlinear systems can exhibit multiple equilibria, limit cycles, and complex basins of attraction
- Conservative systems (like undamped oscillators) show centers and closed orbits; dissipative systems show attracting equilibria or limit cycles
Quick Reference Table
|
| Stability from eigenvalues | Negative real parts → stable; positive → unstable; mixed → saddle |
| Oscillatory behavior | Complex eigenvalues; spirals (damped/growing) or centers (pure imaginary) |
| Finding equilibria | Nullcline intersections; solve x˙=0 and y˙=0 simultaneously |
| Linearization tool | Jacobian matrix evaluated at equilibrium point |
| Node vs. spiral distinction | Real eigenvalues → node; complex eigenvalues → spiral |
| Global periodic behavior | Limit cycles (nonlinear only); Poincaré-Bendixson theorem |
| Parameter sensitivity | Bifurcations: saddle-node, transcritical, pitchfork, Hopf |
| Quick stability check | Trace-determinant: τ<0 and Δ>0 for asymptotic stability |
Self-Check Questions
-
If a system has eigenvalues λ=−2±3i at an equilibrium point, what type of equilibrium is it, and is it stable? What would change if the real part were positive?
-
Compare and contrast nullclines and eigenvectors—both give directional information, but how do they differ in what they tell you about the system?
-
A nonlinear system has a bounded region with no equilibrium points inside. What does the Poincaré-Bendixson theorem tell you must exist, and how does this differ from what's possible in linear systems?
-
Given a Jacobian with tr(J)=0 and det(J)>0, classify the equilibrium. Why does linearization fail to fully characterize stability in this case?
-
Explain how a Hopf bifurcation differs from a saddle-node bifurcation in terms of what happens to equilibria and what new structures might appear.