Graduate → Differential Equations → Ordinary Differential Equations ↓
Stability Analysis
Stability analysis is an important part of understanding ordinary differential equations (ODEs) in mathematics, especially at the undergraduate level. It involves determining how solutions to ODEs behave over time, especially in response to small disturbances or changes in initial conditions.
1. Introduction to sustainability
When dealing with differential equations, it is essential to understand whether the solutions remain close to a particular state over time. This concept is known as stability. A stable system will return to its equilibrium state after a small disturbance, while an unstable system will drift away from equilibrium.
2. Equilibrium point
Before considering stability, we need to identify the equilibrium points of the system. The equilibrium point of an ordinary differential equation is the point where the derivative is zero.
Consider the simple ODE: dx/dt = f(x) An equilibrium point x0 satisfies f(x0) = 0.
These points are like the resting states of the dynamic system. If you imagine a ball in a bowl, the bottom of the bowl will be the equilibrium point.
Visual example
3. Types of stability
3.1 Asymptotic stability
An equilibrium point is asymptotically stable if solutions starting near the equilibrium point not only remain close to equilibrium but also tend toward equilibrium as time approaches infinity.
3.2 Lyapunov stability
An equilibrium point is Lyapunov stable if whenever solutions start close enough to the equilibrium point, they remain close, but as time goes on they may not necessarily converge to the equilibrium point.
3.3 Volatility
An equilibrium point is unstable if solutions that start close together do not remain close; rather, they drift apart from the equilibrium point.
4. Linear stability analysis
For linear systems, determining stability is relatively simple. Consider this system:
dx/dt = Ax
where A is a matrix of constants. The stability of this system can be determined by looking at the eigenvalues of the matrix A.
- If the real parts of all eigenvalues are negative, the equilibrium is asymptotically stable.
- If the real part of any eigenvalue is positive, the equilibrium is unstable.
- If the real part of the eigenvalue is zero, further analysis is needed to determine stability.
Example
Consider the linear system: dx/dt = [[-2, 0], [0, -3]] x The eigenvalues are -2 and -3, both negative. Thus, this system is asymptotically stable.
5. Non-linear stability analysis
For non-linear systems, stability analysis is more complicated. The linearization method allows us to study a non-linear system by approximating it with a linear point near an equilibrium point.
If you have a non-linear system:
dx/dt = f(x)
The system can be linearized around the equilibrium point x0 by computing the Jacobian matrix J of the partial derivatives at x0:
J = [ [∂f₁/∂x₁, ..., ∂f₁/∂xₙ], ..., [∂fₙ/∂x₁, ..., ∂fₙ/∂xₙ] ]
Example
Consider:
dx/dt = y dy/dt = -x + y - x²
Balance point: (0, 0)
Jacobian matrix at (0, 0):
J = [ [0, 1], [-1, 1] ]
The eigenvalues of J determine the stability of the origin.
6. Lyapunov's direct method
Another way to study stability is Lyapunov's direct method. This method does not require solving a differential equation.
Find a scalar function V(x), called the Lyapunov function, that satisfies:
- V(x) > 0, and V(0) = 0 for x ≠ 0
- dV/dt ≤ 0 along the trajectory of the system
If such a function exists, then the point is Lyapunov stable. If dV/dt is actually <, then the point is asymptotically stable.
Example
Consider the system:
dx/dt = -x
We choose V(x) = x². Then:
dV/dt = 2x(dx/dt) = 2x(-x) = -2x² ≤ 0
This system is asymptotically stable at x = 0.
7. Visualizing stability with phasor diagrams
The phase diagram is a graphical representation of the trajectories of a dynamic system in the phase plane. Let's look at a simple case:
Example
In the above phase diagram, the trajectories converge to a point, indicating an asymptotically stable system.
8. Summary and conclusion
Stability analysis of ordinary differential equations provides insight into the long-term behavior of dynamic systems. By examining equilibrium points and understanding concepts such as asymptotic stability, we can predict whether small perturbations will cause a system to recover or fall apart.
Various methods including linearization, Lyapunov methods, and phase delineation provide us with tools to assess stability for both linear and nonlinear systems.
In practice, these methods are invaluable in fields such as engineering, physics, and economics, where they help design stable systems and understand complex dynamic phenomena.