Equilibrium Jacobian Linearization: Local Linear Models

By Vegard Gjerde Based on Masterful Learning 10 min read Published
equilibrium-jacobian-linearization differential-equations math learning-strategies

Equilibrium Jacobian Linearization shifts an autonomous system to perturbation coordinates around an equilibrium, then replaces the nonlinear vector field by its Jacobian matrix at that equilibrium. It applies when FF is differentiable near x\mathbf{x}_* and F(x)=0\mathbf{F}(\mathbf{x}_*)=\mathbf{0}. The fast check is to verify the candidate is truly an equilibrium before using JF(x)J_F(\mathbf{x}_*) as the local model.

Unisium hero image titled Equilibrium Jacobian Linearization showing the principle equation and a conditions card.
Equilibrium Jacobian linearization turns the nearby perturbation dynamics into the local linear system uJF(x)u\mathbf{u}^{\prime}\approx J_F(\mathbf{x}_*)\mathbf{u}.

On this page: The Principle | Conditions | Failure Modes | EE Questions | Retrieval Practice | Practice Ground | Solve a Problem | Related Guides | FAQ


The Principle

The move: write the state as x=x+u\mathbf{x}=\mathbf{x}_*+\mathbf{u}, then replace F(x+u)\mathbf{F}(\mathbf{x}_*+\mathbf{u}) by its first-order Jacobian term at the equilibrium.

The invariant: the equilibrium becomes the origin in u\mathbf{u}-coordinates, and the linear model matches the vector field’s first-order Taylor term at that equilibrium. Higher-order nonlinear terms are deliberately discarded, so this is a local approximation, not an equivalent rewrite.

Pattern:

x=F(x),x=x+u,F(x)=0uJF(x)u\mathbf{x}^{\prime}=\mathbf{F}(\mathbf{x}),\quad \mathbf{x}=\mathbf{x}_*+\mathbf{u},\quad \mathbf{F}(\mathbf{x}_*)=\mathbf{0} \Longrightarrow \mathbf{u}^{\prime}\approx J_F(\mathbf{x}_*)\mathbf{u}
Legal routeIllegal route
F(0)=0, JF(0)=[1001]u[1001]u\mathbf{F}(\mathbf{0})=\mathbf{0},\ J_F(\mathbf{0})=\begin{bmatrix}1&0\\0&-1\end{bmatrix}\Longrightarrow \mathbf{u}^{\prime}\approx \begin{bmatrix}1&0\\0&-1\end{bmatrix}\mathbf{u}F(x)0uJF(x)u\mathbf{F}(\mathbf{x}_*)\neq\mathbf{0}\Longrightarrow \mathbf{u}^{\prime}\approx J_F(\mathbf{x}_*)\mathbf{u}

The illegal route skips the equilibrium check. Around a non-equilibrium point, the shifted equation has a constant forcing term and is not this equilibrium linearization.


Conditions of Applicability

Condition: FdifferentiablenearxF differentiable near \mathbf{x}_*; F(x)=0\mathbf{F}(\mathbf{x}_*)=0

Before applying, check: verify that the candidate point makes every component of F\mathbf{F} zero, then compute the Jacobian at that same point.

If the condition is violated: the local system may need an affine constant term, or the Jacobian may not exist, so the equilibrium linear model is not licensed.

  • The system must be autonomous in the form x=F(x)\mathbf{x}^{\prime}=\mathbf{F}(\mathbf{x}).
  • FF must be differentiable near the equilibrium so the Jacobian is meaningful as the first-order local term.
  • The same x\mathbf{x}_* must be used in the equilibrium check, coordinate shift, and Jacobian evaluation.
  • Linearization predicts nearby local behavior. It does not prove global behavior far from the equilibrium.

Want the complete framework behind this guide? Read Masterful Learning.


Common Failure Modes

Failure mode: compute JF(x)J_F(\mathbf{x}_*) before checking F(x)=0\mathbf{F}(\mathbf{x}_*)=\mathbf{0} -> the result may be a linearization about a moving point, not an equilibrium model.

Debug: plug x\mathbf{x}_* into the vector field first; if any component is nonzero, stop and do not use the equilibrium formula.

Failure mode: treat uJF(x)u\mathbf{u}^{\prime}\approx J_F(\mathbf{x}_*)\mathbf{u} as an exact nonlinear solution -> the higher-order terms disappear from the model.

Debug: say “local first-order model” whenever you write the approximation sign.


Elaborative Encoding

Use these questions to build deep understanding. (See Elaborative Encoding for the full method.)

Within the Principle

  • What does the perturbation vector u=xx\mathbf{u}=\mathbf{x}-\mathbf{x}_* measure?
  • Why does the constant term vanish only when F(x)=0\mathbf{F}(\mathbf{x}_*)=\mathbf{0}?

For the Principle

  • Before linearizing, what two checks protect the move from becoming an invalid local model?
  • Why is the Jacobian evaluated at the equilibrium rather than at a nearby moving state?

Between Principles

Generate an Example

  • Create a two-dimensional autonomous system where F(0)=0\mathbf{F}(\mathbf{0})=\mathbf{0}, then change one constant term so the same point becomes a near miss.

Retrieval Practice

Answer from memory, then click to reveal and check. (See Retrieval Practice for the full method.)

State the move in one sentence: _____Shift an autonomous system to perturbation coordinates around an equilibrium, then use the Jacobian at that equilibrium as the local linear model.
Write the canonical pattern: _____x=F(x), x=x+u, F(x)=0uJF(x)u\mathbf{x}^{\prime}=\mathbf{F}(\mathbf{x}),\ \mathbf{x}=\mathbf{x}_*+\mathbf{u},\ \mathbf{F}(\mathbf{x}_*)=\mathbf{0} \Rightarrow \mathbf{u}^{\prime}\approx J_F(\mathbf{x}_*)\mathbf{u}
State the canonical condition: _____Fdifferentiablenearx;F(x)=0F differentiable near \mathbf{x}_*; \mathbf{F}(\mathbf{x}_*)=0

Practice Ground

Use these exercises to build move-selection fluency. (See Self-Explanation for how to learn from worked examples.)

Procedure Walkthrough

Starting from the autonomous system

x=[x+y2y+xy],x=[00],\mathbf{x}^{\prime}= \begin{bmatrix} x+y^2\\ -y+xy \end{bmatrix}, \qquad \mathbf{x}_*=\begin{bmatrix}0\\0\end{bmatrix},

reach the local linear model in perturbation coordinates.

StepExpressionOperation
0F(0)=[00]\mathbf{F}(\mathbf{0})=\begin{bmatrix}0\\0\end{bmatrix}Check that the candidate is an equilibrium.
1JF(x,y)J_F(x,y) has rows (1,2y)(1,2y) and (y,1+x)(y,-1+x).Compute the Jacobian matrix.
2JF(0)=[1001]J_F(\mathbf{0})=\begin{bmatrix}1&0\\0&-1\end{bmatrix}Evaluate the Jacobian at the equilibrium.
3x=0+u\mathbf{x}=\mathbf{0}+\mathbf{u}Shift the equilibrium to the origin.
4u[1001]u\mathbf{u}^{\prime}\approx \begin{bmatrix}1&0\\0&-1\end{bmatrix}\mathbf{u}Write the local linearized system.

Drills

Micro-Chains

Reach the local linear model. Assume FF is differentiable near x\mathbf{x}_* and verify F(x)=0\mathbf{F}(\mathbf{x}_*)=\mathbf{0} before linearizing.

F(x,y)=[2x+y23y+x2],x=[00]\mathbf{F}(x,y)=\begin{bmatrix}2x+y^2\\-3y+x^2\end{bmatrix},\qquad \mathbf{x}_*=\begin{bmatrix}0\\0\end{bmatrix}
Reveal

First,

F(0)=[00].\mathbf{F}(\mathbf{0})=\begin{bmatrix}0\\0\end{bmatrix}.

Then

JF(x,y)=[22y2x3],JF(0)=[2003].J_F(x,y)=\begin{bmatrix}2&2y\\2x&-3\end{bmatrix}, \qquad J_F(\mathbf{0})=\begin{bmatrix}2&0\\0&-3\end{bmatrix}.

So the local linear model is

u[2003]u.\mathbf{u}^{\prime}\approx \begin{bmatrix}2&0\\0&-3\end{bmatrix}\mathbf{u}.

Reach the local linear model near the equilibrium. Assume differentiability where needed.

F(x,y)=[yx+x3],x=[00]\mathbf{F}(x,y)=\begin{bmatrix}y\\-x+x^3\end{bmatrix},\qquad \mathbf{x}_*=\begin{bmatrix}0\\0\end{bmatrix}
Reveal

The equilibrium check passes because F(0)=0\mathbf{F}(\mathbf{0})=\mathbf{0}.

JF(x,y)=[011+3x20],JF(0)=[0110].J_F(x,y)=\begin{bmatrix}0&1\\-1+3x^2&0\end{bmatrix}, \qquad J_F(\mathbf{0})=\begin{bmatrix}0&1\\-1&0\end{bmatrix}.

Therefore

u[0110]u.\mathbf{u}^{\prime}\approx \begin{bmatrix}0&1\\-1&0\end{bmatrix}\mathbf{u}.

Reject or linearize. Assume FF is differentiable near the candidate point.

F(x,y)=[x+1y],x=[00]\mathbf{F}(x,y)=\begin{bmatrix}x+1\\-y\end{bmatrix},\qquad \mathbf{x}_*=\begin{bmatrix}0\\0\end{bmatrix}
Reveal

Reject this as equilibrium Jacobian linearization at 0\mathbf{0}.

F(0)=[10]0.\mathbf{F}(\mathbf{0})=\begin{bmatrix}1\\0\end{bmatrix}\neq \mathbf{0}.

The point is not an equilibrium, so the condition fails before the Jacobian matters.


Reach the local linear model near the equilibrium. Assume differentiability where needed.

F(x,y)=[xy2x+y2y],x=[00]\mathbf{F}(x,y)=\begin{bmatrix}x-y\\2x+y^2-y\end{bmatrix},\qquad \mathbf{x}_*=\begin{bmatrix}0\\0\end{bmatrix}
Reveal

The equilibrium check passes:

F(0)=[00].\mathbf{F}(\mathbf{0})=\begin{bmatrix}0\\0\end{bmatrix}.

The Jacobian is

JF(x,y)=[1122y1],JF(0)=[1121].J_F(x,y)=\begin{bmatrix}1&-1\\2&2y-1\end{bmatrix}, \qquad J_F(\mathbf{0})=\begin{bmatrix}1&-1\\2&-1\end{bmatrix}.

So

u[1121]u.\mathbf{u}^{\prime}\approx \begin{bmatrix}1&-1\\2&-1\end{bmatrix}\mathbf{u}.

Action Labels

What was done between these two steps?

x=F(x),x=x+uu=F(x+u)\mathbf{x}^{\prime}=\mathbf{F}(\mathbf{x}),\quad \mathbf{x}=\mathbf{x}_*+\mathbf{u} \quad \Longrightarrow \quad \mathbf{u}^{\prime}=\mathbf{F}(\mathbf{x}_*+\mathbf{u})
Reveal

The state was shifted to perturbation coordinates around the equilibrium. Since x\mathbf{x}_* is constant, x=u\mathbf{x}^{\prime}=\mathbf{u}^{\prime}.


What was done between these two steps, assuming FF is differentiable near x\mathbf{x}_* and F(x)=0\mathbf{F}(\mathbf{x}_*)=\mathbf{0}?

u=F(x+u)uJF(x)u\mathbf{u}^{\prime}=\mathbf{F}(\mathbf{x}_*+\mathbf{u}) \quad \Longrightarrow \quad \mathbf{u}^{\prime}\approx J_F(\mathbf{x}_*)\mathbf{u}
Reveal

The nonlinear vector field was replaced by its first-order Taylor term at the equilibrium. The constant term is zero because F(x)=0\mathbf{F}(\mathbf{x}_*)=\mathbf{0}.


Which proposed steps are eligible equilibrium Jacobian linearizations?

A. F(x)=0\mathbf{F}(\mathbf{x}_*)=\mathbf{0} and FF is differentiable near x\mathbf{x}_*, so use uJF(x)u\mathbf{u}^{\prime}\approx J_F(\mathbf{x}_*)\mathbf{u}.
B. F(x)0\mathbf{F}(\mathbf{x}_*)\neq\mathbf{0}, but compute JF(x)J_F(\mathbf{x}_*) and use uJF(x)u\mathbf{u}^{\prime}\approx J_F(\mathbf{x}_*)\mathbf{u} anyway.
C. FF has a corner at x\mathbf{x}_*, so estimate slopes from one side and call that JF(x)J_F(\mathbf{x}_*).

Reveal

A is eligible.

B fails the equilibrium condition. C fails the differentiability condition, so the Jacobian linearization is not licensed.


What was done between these two steps?

JF(x,y): (2x,1),(y,xy)JF(1,0): (2,1),(0,1)J_F(x,y):\ (2x,1),(y,x-y) \quad \Longrightarrow \quad J_F(1,0):\ (2,1),(0,1)
Reveal

The Jacobian was evaluated at the chosen equilibrium point. This is only the right matrix for equilibrium linearization if the same point also satisfies F(1,0)=0\mathbf{F}(1,0)=\mathbf{0}.


Forward Steps

Apply the principle once after the equilibrium check has passed.

JF(x)=[2410]J_F(\mathbf{x}_*)=\begin{bmatrix}-2&4\\1&0\end{bmatrix}
Reveal

With x=x+u\mathbf{x}=\mathbf{x}_*+\mathbf{u} and the equilibrium condition already checked,

u[2410]u.\mathbf{u}^{\prime}\approx \begin{bmatrix}-2&4\\1&0\end{bmatrix}\mathbf{u}.

Apply the principle once after the equilibrium check has passed.

JF(x)=[0152]J_F(\mathbf{x}_*)=\begin{bmatrix}0&1\\-5&-2\end{bmatrix}
Reveal

The local linearized perturbation system is

u[0152]u.\mathbf{u}^{\prime}\approx \begin{bmatrix}0&1\\-5&-2\end{bmatrix}\mathbf{u}.

Reject or complete the step. A student writes uJF(x)u\mathbf{u}^{\prime}\approx J_F(\mathbf{x}_*)\mathbf{u} before checking F(x)=0\mathbf{F}(\mathbf{x}_*)=\mathbf{0}.

Reveal

Reject the step as incomplete. For this principle, the equilibrium check is part of the move, not a later note.


Solve a Problem

Apply what you’ve learned with Problem Solving.

Problem: Starting from the autonomous system

x=[x+y2x+yy2],x=[00],\mathbf{x}^{\prime}= \begin{bmatrix} x+y\\ -2x+y-y^2 \end{bmatrix}, \qquad \mathbf{x}_*=\begin{bmatrix}0\\0\end{bmatrix},

reach the local linearized system using Equilibrium Jacobian Linearization.

Full solution
StepExpressionMove
0F(0)=[00]\mathbf{F}(\mathbf{0})=\begin{bmatrix}0\\0\end{bmatrix}Check the equilibrium condition.
1x=0+u\mathbf{x}=\mathbf{0}+\mathbf{u}Shift the equilibrium to perturbation coordinates.
2JF(x,y)=[11212y]J_F(x,y)=\begin{bmatrix}1&1\\-2&1-2y\end{bmatrix}Compute the Jacobian.
3JF(0)=[1121]J_F(\mathbf{0})=\begin{bmatrix}1&1\\-2&1\end{bmatrix}Evaluate at the same equilibrium.
4u[1121]u\mathbf{u}^{\prime}\approx \begin{bmatrix}1&1\\-2&1\end{bmatrix}\mathbf{u}Write the local linearized system.


FAQ

What is Equilibrium Jacobian Linearization?

Equilibrium Jacobian Linearization is the move that approximates an autonomous nonlinear system near an equilibrium by a linear perturbation system. After writing x=x+u\mathbf{x}=\mathbf{x}_*+\mathbf{u}, the local model is uJF(x)u\mathbf{u}^{\prime}\approx J_F(\mathbf{x}_*)\mathbf{u}.

When is Equilibrium Jacobian Linearization valid?

It is valid when FF is differentiable near x\mathbf{x}_* and F(x)=0\mathbf{F}(\mathbf{x}_*)=\mathbf{0}. The differentiability condition gives the Jacobian, and the equilibrium condition removes the constant term from the shifted local model.

Why do we shift to perturbation coordinates?

The shift x=x+u\mathbf{x}=\mathbf{x}_*+\mathbf{u} makes u\mathbf{u} measure displacement from the equilibrium. That turns the equilibrium itself into u=0\mathbf{u}=\mathbf{0}, so the nearby motion can be compared with a linear system around the origin.

Is the linearized system an exact replacement for the nonlinear system?

No. It is a first-order local approximation. The nonlinear terms beyond the Jacobian are ignored, so conclusions from the linearized system need the usual local-stability interpretation rather than a global exact-solution interpretation.

What happens if the candidate point is not an equilibrium?

Then F(x)\mathbf{F}(\mathbf{x}_*) is not zero, and the shifted equation has a constant term. You may still approximate the vector field locally, but it is not the equilibrium Jacobian linearization described by this principle.

Once the perturbation system is uAu\mathbf{u}^{\prime}\approx A\mathbf{u} with A=JF(x)A=J_F(\mathbf{x}_*), eigenvalues of AA describe the local linear modes. If the linearized matrix has eigenvalues with zero real part, the linear model may be inconclusive for the nonlinear system’s stability. That builds on Eigenvalue-Eigenvector Solution Mode rather than replacing the equilibrium check.


How This Fits in Unisium

Within the differential equations subdomain, Unisium treats equilibrium linearization as a condition-first move: find the equilibrium, shift coordinates, compute the Jacobian, then read the local system. The Unisium Study System pairs that sequence with retrieval practice, self-explanation, and compact problem-solving chains so the equilibrium condition stays visible while the matrix work becomes fluent.

Ready to practice differential equations with structure? Check access and join the Unisium waitlist or explore the complete framework in Masterful Learning.

Masterful Learning book cover

Masterful Learning

The book behind these guides: a study system for physics, math, & programming built on retrieval, connection, explanation, and problem solving.

Ready to apply this strategy?

Unisium turns these evidence-based techniques into guided study sessions for math and physics. Unisium is currently in early access. See pricing, availability, and join the waitlist.

Check Unisium Access and Pricing Read More Guides

Already have access? Sign in