4.6.22Ordinary Differential Equations

Phase plane analysis — trajectories, critical points

1,868 words8 min readdifficulty · medium1 backlinks

1. What is a trajectory? (WHAT)

WHY autonomous matters: because tt is absent, the velocity vector (x˙,y˙)(\dot x,\dot y) depends only on where you are, not when. So through each point there is a unique direction → trajectories never cross (except at critical points).


2. Critical points (WHAT & WHY)

HOW to find them: solve the simultaneous algebraic equations f=0, g=0f=0,\ g=0. These are the only points where the direction (x˙,y˙)(\dot x,\dot y) is undefined as a "where to go" arrow, hence the only places trajectories may meet.


3. Classifying critical points — linearization (the 80/20 core)

Derivation from scratch. Shift origin: let u=xx0, v=yy0u=x-x_0,\ v=y-y_0. Taylor expand: u˙=fxu+fyv+(higher order),v˙=gxu+gyv+\dot u = f_x u + f_y v + \text{(higher order)},\qquad \dot v = g_x u + g_y v + \cdots all partials evaluated at the critical point. Drop higher-order terms: (u˙v˙)=J(uv),J=(fxfygxgy).\begin{pmatrix}\dot u\\ \dot v\end{pmatrix}=J\begin{pmatrix}u\\ v\end{pmatrix},\qquad J=\begin{pmatrix} f_x & f_y\\ g_x & g_y\end{pmatrix}. JJ is the Jacobian at the critical point.

WHY eλte^{\lambda t}: try u=veλt\mathbf u = \mathbf v\,e^{\lambda t}. Then λveλt=JveλtJv=λv\lambda \mathbf v e^{\lambda t} = J\mathbf v e^{\lambda t}\Rightarrow J\mathbf v = \lambda\mathbf v. So λ\lambda must be an eigenvalue. The sign of Reλ\operatorname{Re}\lambda says grow (unstable) or decay (stable); imaginary part says rotate.

Classification table

Eigenvalues λ1,2\lambda_{1,2} Type Stability
Real, same sign, both <0<0 Node (sink) Stable
Real, same sign, both >0>0 Node (source) Unstable
Real, opposite signs Saddle Unstable
Complex, Re<0\operatorname{Re}<0 Spiral (focus) Stable
Complex, Re>0\operatorname{Re}>0 Spiral Unstable
Pure imaginary (Re=0\operatorname{Re}=0) Center (neutral)
Figure — Phase plane analysis — trajectories, critical points

4. Worked Example A — a linear system (full derivation)

System: x˙=x+y,y˙=2x+4y\dot x = x + y,\quad \dot y = -2x + 4y? Let's instead use clean one: x˙=x+3y,y˙=3xy?\dot x = -x + 3y,\qquad \dot y = -3x - y?

Let me do a transparent saddle: x˙=x+2y,y˙=2x+y.\dot x = x + 2y,\qquad \dot y = 2x + y.

Step 1 — critical point. Set x+2y=0, 2x+y=0x+2y=0,\ 2x+y=0. Why? equilibrium needs both rates zero. Subtract: 3x+3y=03x+3y=0\Rightarrow... solve: only solution (0,0)(0,0).

Step 2 — Jacobian. Here f,gf,g are already linear so J=(1221)J=\begin{pmatrix}1&2\\2&1\end{pmatrix}. Why? fx=1,fy=2,gx=2,gy=1f_x=1,f_y=2,g_x=2,g_y=1.

Step 3 — eigenvalues. τ=2, Δ=1122=3\tau=2,\ \Delta=1\cdot1-2\cdot2=-3. λ22λ3=0(λ3)(λ+1)=0λ=3,1.\lambda^2 -2\lambda -3=0\Rightarrow(\lambda-3)(\lambda+1)=0\Rightarrow \lambda=3,-1. Step 4 — classify. Real, opposite signs \Rightarrow saddle, unstable. Why this step? opposite signs mean one direction grows, one shrinks — the hallmark of a saddle.


5. Worked Example B — a spiral (nonlinear, linearized)

x˙=y+x(1x2y2),y˙=x+y(1x2y2).\dot x = -y + x(1-x^2-y^2),\qquad \dot y = x + y(1-x^2-y^2).

Step 1 — critical point. At origin the cubic terms vanish: x˙=y, y˙=x\dot x=-y,\ \dot y=x, both zero only at (0,0)(0,0). Why drop cubics? near origin x2+y2x^2+y^2 is tiny.

Step 2 — Jacobian at origin. Linear part is x˙=y+x, y˙=x+y\dot x=-y+x,\ \dot y=x+y (since 1x2y211-x^2-y^2\to1). So J=(1111).J=\begin{pmatrix}1&-1\\1&1\end{pmatrix}. Why this step? differentiate including the x1x\cdot1 term from x(1)x(1-\dots).

Step 3 — eigenvalues. τ=2, Δ=1+1=2\tau=2,\ \Delta=1+1=2, discriminant τ24Δ=48=4<0\tau^2-4\Delta=4-8=-4<0. λ=1±i.\lambda = 1\pm i. Step 4 — classify. Complex with Reλ=1>0\operatorname{Re}\lambda=1>0\Rightarrow unstable spiral outward. (In fact trajectories spiral out toward a circular limit cycle r=1r=1 — that's the full nonlinear truth.)


6. Common mistakes (Steel-man + fix)


Recall Feynman: explain to a 12-year-old

Imagine a windy field where a leaf is blown by wind that's the same at the same spot every day. Drop the leaf anywhere and it follows a path — that's a trajectory. Some spots have no wind at all; a leaf there just sits — those are critical points. Some no-wind spots suck leaves in (sinks), some blow them away (sources), some swirl them (spirals), and some are passes between hills where leaves come close then get flung off (saddles). By only studying the no-wind spots and the wind right around them, you can guess where every leaf ends up — without watching each leaf the whole day.


Flashcards

What is the phase plane?
The xyxy-plane in which solution curves (x(t),y(t))(x(t),y(t)) of an autonomous 2D system are drawn.
Define a trajectory/orbit.
The path traced by a solution (x(t),y(t))(x(t),y(t)) as tt varies.
How is a critical point defined?
A point where f(x0,y0)=0f(x_0,y_0)=0 and g(x0,y0)=0g(x_0,y_0)=0 (system at rest).
Trajectory slope formula (eliminating tt)?
dydx=g(x,y)f(x,y)\dfrac{dy}{dx}=\dfrac{g(x,y)}{f(x,y)}.
Why do trajectories never cross except at critical points?
Autonomous field gives a unique direction at each ordinary point (uniqueness of solutions).
What matrix governs local behavior near a critical point?
The Jacobian J=fxfygxgyJ=\begin{smallmatrix}f_x&f_y\\g_x&g_y\end{smallmatrix} evaluated there.
Characteristic equation for eigenvalues of JJ?
λ2(trJ)λ+detJ=0\lambda^2-(\operatorname{tr}J)\lambda+\det J=0.
Real eigenvalues, opposite signs → which type?
Saddle (unstable).
Real eigenvalues, same negative sign → ?
Stable node (sink).
Complex eigenvalues with negative real part → ?
Stable spiral (focus).
Pure imaginary eigenvalues → ?
Center (linearization inconclusive for nonlinear systems).
What does detJ<0\det J<0 always imply?
A saddle point.
Which sign determines stability when detJ>0\det J>0?
Sign of trJ\operatorname{tr}J: negative ⇒ stable, positive ⇒ unstable.
Discriminant separating nodes from spirals?
τ24Δ\tau^2-4\Delta: 0\ge0 node, <0<0 spiral (τ=trJ,Δ=detJ\tau=\operatorname{tr}J,\Delta=\det J).

Connections

Concept Map

too hard to solve

solutions traced as

collection forms

t absent so velocity fixed

implies

find where nothing moves

solve f=0 and g=0

except at

Taylor expand near point

first-order terms give

eigenvalues lambda

divide to drop t

Coupled autonomous ODEs

Phase plane analysis

Trajectories

Phase portrait

Unique direction per point

Trajectories never cross

Critical points

Equilibrium locations

Linearization

Jacobian J

Classify stability

Slope dy/dx = g/f

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, kabhi-kabhi do ODE aapas mein itne coupled hote hain ki x(t)x(t) aur y(t)y(t) ka exact formula nikalna namumkin ho jata hai. Tab hum phase plane ka trick use karte hain: tt ko bhool jao, sirf xyxy-plane mein curve ka shape dekho. Har solution ek trajectory banata hai, jaise hawa mein patta udta hai. Important baat: autonomous system mein har point pe direction fix hai, isliye do trajectory kabhi cross nahi karte (sirf critical point pe milte hain).

Critical point wo jagah hai jahan f=0f=0 aur g=0g=0 — yaani system ekdum rest pe, kuch move nahi karta. Yeh points pure flow ko organize karte hain, jaise floor pe drain aur fountain paani ka pattern decide karte hain. Inhe dhoondhne ke liye bas f=0,g=0f=0,g=0 algebraically solve karo.

Ab in points ke around behaviour janne ke liye hum linearize karte hain — Taylor expand karke Jacobian JJ banate hain. JJ ke eigenvalues sab kuch bata dete hain: real opposite signs = saddle (unstable), dono negative = stable node, complex with negative real part = stable spiral, pure imaginary = center. Shortcut: detJ<0\det J<0 hamesha saddle; agar detJ>0\det J>0 to trace ka sign stability decide karta hai (negative = stable). Bas itna yaad rakho aur 80% questions ho jayenge — baaki sirf discriminant τ24Δ\tau^2-4\Delta se node vs spiral confirm kar lo.

Go deeper — visual, from zero

Test yourself — Ordinary Differential Equations

Connections