---
title: Singular Perturbation Lyapunov Function
url: https://www.emergentmind.com/topics/singular-perturbation-lyapunov-function
type: topic
---

# Singular Perturbation Lyapunov Function

A Singular Perturbation Lyapunov Function is a composite Lyapunov candidate used to verify stability, robustness, and convergence properties in dynamical systems exhibiting multiple distinct time scales, typically parameterized by a small positive parameter $\varepsilon$. Such systems naturally arise in various engineering, physics, and control contexts and are formally expressed as two (or more) coupled subsystems: the so-called slow and fast (boundary-layer) dynamics. The singular perturbation Lyapunov framework leverages appropriately constructed Lyapunov functions for each time scale and combines them—via convex combinations or structured quadratic forms—to obtain system-level stability guarantees, including asymptotic, fixed-time, or input-to-state stability, even under disturbances, stochastic effects, or hybrid (continuous/discrete) switching.

## 1. Mathematical Structure of Singularly Perturbed Systems

Singularly perturbed models partition the state vector $y=(x, z) \in \mathbb{R}^{n_1} \times \mathbb{R}^{n_2}$ into “slow” states $x$ and “fast” states $z$, and feature dynamics of the prototypical form
\[
\begin{aligned}
\dot{x} &= f(x, z, u),\\
\varepsilon\,\dot{z} &= g(x, z, u),
\end{aligned}
\]
where $u$ is an input and $\varepsilon > 0$ captures the time-scale separation. As $\varepsilon \to 0$, $z$ converges rapidly to a quasi–steady-state manifold $z = h(x)$, yielding two limiting subsystems:
- The reduced (slow) system: $\dot{x} = f(x, h(x), u)$,
- The boundary-layer (fast) system: $\frac{dy}{d\tau} = g(x, y+h(x), u)$ with $y = z - h(x)$ and rescaled time $\tau = t/\varepsilon$.

Extensions to hybrid and stochastic settings include flow and jump sets, outer-semicontinuous and convex-valued right-hand sides, and random or measure-driven events, as formalized in [2310.09712]. The composite flow map is frequently written as $F_\varepsilon(x, z) = F_x(x, z) \times \varepsilon^{-1} F_z(x, z)$.

## 2. Construction of Composite Lyapunov Functions

Standard methodology requires establishing individual Lyapunov (or Foster–Lyapunov) functions for the slow and fast subsystems:

- Fast subsystem: $W(x, z)$, positive definite in $z$ (relative to the manifold $\mathcal{M}(x)$), satisfying
  \[
  \max_{\nu \in \partial_z W} \langle \nu, f_z \rangle \le -k_z \varphi_z^2(\lvert z \rvert_{\mathcal{M}(x)})
  \]
  for some $k_z > 0$, and appropriate positive-definiteness and decrease-along-jumps conditions.

- Slow subsystem: $V(x)$, positive definite with respect to a compact set $\mathcal{A}$, satisfying
  \[
  \max_{\nu \in \partial V} \langle \nu, \tilde{f}_x \rangle \le -k_x \varphi_x^2(x) + \mu_F\,\mathcal{O}(x),
  \]
  where $\tilde{f}_x$ is the reduced slow dynamics and $\mathcal{O}$ a target or neighborhood set.

The composite Lyapunov function is then constructed as
\[
E_{\theta}(x, z) = (1 - \theta)\,V(x) + \theta\,W(x, z), \quad \theta \in (0, 1),
\]
or, equivalently, as a convex combination, with $\theta$ determined according to the coupling strengths between subsystems. In the context of fixed-time or input-to-state stability, additional “mixed-power” or “half-power” transforms (e.g., $\tilde{V}(x) = V(x)^{a_1/2} + V(x)^{a_2/2}$) are employed to facilitate bounding the interconnection terms [2412.16797, 2408.16905].

## 3. Sufficient Conditions and Main Theorems

Stability and convergence of the overall system via the composite Lyapunov approach hinge on several key sufficient conditions:

- Positive-definiteness and dissipativity of $V$ and $W$ individually (with possible power-type or squared dissipation, depending on the setting).
- Explicit interconnection bounds on cross-terms resulting from the slow–fast coupling, typically formulated as:
  \[
  I_1(x, y) \le \chi_1 \tilde{V} \tilde{W} + \delta_1 (\tilde{V})^2 + c_1 (\tilde{W})^2,\\
  I_2(x, y) \le \chi_2 \tilde{V} \tilde{W} + \delta_2 (\tilde{V})^2 + c_2 (\tilde{W})^2,
  \]
  with a requirement that at least one of the self-dissipation parameters (e.g., $\delta_1$ or $\delta_2$) be sufficiently negative.

- A “smallness” condition for $\varepsilon$: For fixed convex weight $\theta$, there exists $\varepsilon^* > 0$ such that, for all $\varepsilon \in (0, \varepsilon^*)$, the weighted quadratic form in the composite Lyapunov time-derivative is negative definite. For stochastic hybrid systems, analogous average or probabilistic decrease conditions on jumps are imposed [2310.09712].

- In ISS or fixed-time ISS settings, the composite Lyapunov function $\Psi_\zeta(x, y) = \zeta V_r(x) + (1 - \zeta) W(x, y)$ must satisfy
  \[
  \dot \Psi_\zeta \le -c_1 \Psi_\zeta^{\gamma_1} - c_2 \Psi_\zeta^{\gamma_2} + \rho(\|u\|),
  \]
  with exponents $\gamma_1 \in (\max\{a_1, b_1\}, 1)$, $\gamma_2 \in (1, \min\{a_2, b_2\})$, and $c_1, c_2 > 0$ [2412.16797, 2408.16905]. This structure guarantees uniform global fixed-time convergence to the equilibrium (or, in the ISS case, to a disturbance-dependent ball).

## 4. Applications: Stochastic Hybrid, Fixed-Time, and Input-to-State Stability

The singular perturbation Lyapunov method has been extended beyond classical deterministic ODEs to address:

- **Stochastic hybrid dynamical systems**, with both continuous flows (constrained differential inclusions) and jumps (constrained difference inclusions subject to random inputs), as in [2310.09712]. Here, composite Foster–Lyapunov and Lagrange–Foster functions yield criteria for uniform global asymptotic stability in probability (UGASp) of the compact set $\widetilde{\mathcal{A}} = \{ (x, z): x \in \mathcal{A}, z \in \mathcal{M}(x) \}$.

- **Fixed-time and input-to-state stability (ISS)** for nonlinear singularly perturbed systems, where solutions must converge in a time independent of initial state for vanishing disturbance, with explicit uniform upper bounds for settling time given by
  \[
  T^* \le \frac{1}{c_1(1-\gamma_1)} + \frac{1}{c_2(\gamma_2-1)}
  \]
  [2408.16905, 2412.16797]. This is achieved by imposing two-term power-type dissipation inequalities on $V$ and $W$ and explicitly controlling the coupling via the matrix $P(\varepsilon, \theta)$ in the time-derivative of $\Psi$.

- **Hybrid specializations** such as Markov jump-linear systems and deterministic switched systems—with the theory recovering classical results (e.g., Kokotović–O’Malley, Saberi & Kokotović) for ODEs and extending to nonsmooth cases [2310.09712].

## 5. Verification Procedure and Illustrative Examples

Application of the singular perturbation Lyapunov framework typically involves the following steps (see [2412.16797, 2408.16905]):

1. **Model Reduction**: Write the original system, identify the quasi-steady-state manifold $h(x)$, and split into reduced and boundary-layer subsystems.
2. **Lyapunov Construction**: Find candidate $V(x)$, $W(x, y)$ satisfying regularity and dissipation inequalities (with or without explicit ISS/fixed-time extensions).
3. **Interconnection Bounds**: Bound $I_1$, $I_2$ in quadratic form and determine the corresponding constants $\nu_i, \omega_i, \rho_i$.
4. **Parameter Selection**: Choose convex weight $\theta$ (or $\zeta$) to satisfy self-dissipation requirements; compute $\varepsilon^*$ to guarantee negative definiteness.
5. **Composite Lyapunov Derivative and Stability Assessment**: Check that the required decrease condition (possibly stochastic/ISS/fixed-time) is met.
6. **Explicit Convergence Bound**: If applicable, read off explicit settling time from the derived dissipation rates.

Illustrative examples provided in [2412.16797, 2408.16905] include scalar systems with nonlinearities, quadratic cost optimization with singular perturbations, and high-order interconnected systems, demonstrating uniform convergence times and quantitative robustness to bounded disturbance.

## 6. Specializations, Broader Connections, and Methodological Guidelines

The singular perturbation Lyapunov approach unifies several stability frameworks:

- **Quadratic and nonsmooth Lyapunov theory**: Classical quadratic Lyapunov functions are a special case for ODEs with linear/quadratic dynamics [2310.09712].
- **Composite method and power-type dissipation**: Extension to fixed-time and ISS contexts utilizes two-term dissipation and convex combinations, requiring careful cross-term management [2408.16905, 2412.16797].
- **Guidelines**: Researchers are advised to establish strong self-dissipation for individual subsystems, control cross-coupling via quadratic forms, judiciously select $\theta$ and $\varepsilon$, and apply the composite Lyapunov dissipation to deduce global—or probabilistic—stability with explicit convergence rate or time bounds.

A plausible implication is that advances in singular perturbation Lyapunov theory continue to broaden the admissible classes of systems (including stochastic, hybrid, nonlinear, and fixed-time stable) and enable rigorous, rate-optimal stabilization and control techniques in multi-scale dynamical settings.

Source: https://www.emergentmind.com/topics/singular-perturbation-lyapunov-function