---
title: Discrete-Time High-Order CBFs
url: https://www.emergentmind.com/topics/discrete-time-high-order-control-barrier-functions-dhcbfs
type: topic
---

# Discrete-Time High-Order CBFs

Searching arXiv for the cited DHCBF literature to ground the article.
First, I’ll look up the 2022 DHOCBF-MPC paper and a few adjacent discrete-time CBF papers.
Discrete-Time High-Order Control Barrier Functions (DHCBFs), also termed Discrete-Time High Order Control Barrier Functions (DHOCBFs), are barrier-based safety certificates for discrete-time dynamical systems whose safety outputs have arbitrary relative degree. In the formulation developed for the discrete-time nonlinear system \(x_{t+1}=l(x_t,u_t)\), a scalar safety function \(h:\mathbb{R}^n\to\mathbb{R}\) defines the safe set \(\mathcal{C}=\{x\in\mathbb{R}^n\mid h(x)\ge 0\}\), and DHCBFs generalize relative-degree-one discrete-time control barrier functions by constructing a hierarchy of auxiliary functions whose highest-order inequality guarantees forward invariance of the intersection of associated safe sets [2210.04361]. The framework is motivated by safety-critical control and model predictive control (MPC), where many discrete-time systems and constraints exhibit higher relative degree, so enforcing only a one-step barrier condition on \(h\) is either inadequate or overly conservative [2210.04361].

## 1. Formal setting and core definition

The canonical discrete-time model is
\[
x_{t+1}=l(x_t,u_t),
\]
where \(x_t\in\mathcal{X}\subset\mathbb{R}^n\), \(u_t\in\mathcal{U}\subset\mathbb{R}^q\), and \(l:\mathcal{X}\times\mathcal{U}\to\mathbb{R}^n\) is locally Lipschitz [2210.04361]. Safety is encoded by a scalar function \(h:\mathbb{R}^n\to\mathbb{R}\), with safe set
\[
\mathcal{C}\coloneqq \{x\in\mathbb{R}^n\mid h(x)\ge 0\}.
\]
The discrete-time safety requirement is \(x_t\in\mathcal{C}\) for all \(t\in\mathbb{N}\), equivalently \(h(x_t)\ge 0\) for all \(t\in\mathbb{N}\) [2210.04361].

For relative degree \(m\), DHCBFs introduce the sequence
\[
\psi_0(x_t)\coloneqq h(x_t),
\]
and, for \(i=1,\dots,m\),
\[
\psi_i(x_t)\coloneqq \psi_{i-1}(x_{t+1})-\psi_{i-1}(x_t)+\gamma_i\psi_{i-1}(x_t),\quad 0<\gamma_i\le 1,
\]
with associated sets
\[
\mathcal{C}_i\coloneqq \{x_t\in\mathbb{R}^n\mid \psi_{i-1}(x_t)\ge 0\},\quad i=1,\dots,m
\]
[2210.04361]. Thus \(\mathcal{C}_1\) is the original safety set, while higher \(\mathcal{C}_i\) encode multi-step conditions on the evolution of \(h\).

The formal DHOCBF definition is: \(h\) is a discrete-time High Order Control Barrier Function with relative degree \(m\) for \(x_{t+1}=l(x_t,u_t)\) if, for all \(x_t\in \mathcal{C}_1\cap\dots\cap\mathcal{C}_m\),
\[
\psi_{m-1}(x_{t+1})\ge (1-\gamma_m)\psi_{m-1}(x_t),\quad 0<\gamma_m\le 1
\]
[2210.04361]. This is an exponential-type condition on the highest-order auxiliary function.

A closely related discrete-time higher-relative-degree construction appears in distributed safety-critical MPC for nonlinear multi-agent systems, where, for each agent, one defines
\[
\psi_0(x_i(t)):=h(x_i(t)),
\]
\[
\psi_l(x_i(t)):=\Delta \psi_{l-1}(x_i(t),u_i(t))+\alpha_l(\psi_{l-1}(x_i(t))),\quad l=1,\dots,m,
\]
with
\[
\Delta \psi_{l-1}(x_i(t),u_i(t)):=\psi_{l-1}(x_i(t+1))-\psi_{l-1}(x_i(t)),
\]
and \(\alpha_l(\psi_{l-1}(x_i(t)))=\phi\,\psi_{l-1}(x_i(t))\), \(\phi\in(0,1]\) [2508.19678]. In that formulation, \(h\) is a DHCBF if \(\psi_m(x_i(t))\ge 0\) on \(\bigcap_{l=0}^{m-1}\mathcal{C}^l\), where \(\mathcal{C}^l=\{x_i\mid \psi_l(x_i)\ge 0\}\) [2508.19678].

## 2. Forward invariance and relative degree

The defining property of DHCBFs is forward invariance. In the 2022 DHOCBF framework, if \(x_0\in\mathcal{C}_1\cap\dots\cap\mathcal{C}_m\) and the control satisfies
\[
\psi_{m-1}(x_{t+1})\ge (1-\gamma_m)\psi_{m-1}(x_t)
\]
for all \(t\in\mathbb{N}\), then \(\mathcal{C}_1\cap\dots\cap\mathcal{C}_m\) is forward invariant for \(x_{t+1}=l(x_t,u_t)\) [2210.04361]. The recursive construction propagates nonnegativity backward through the hierarchy of \(\psi_i\), so nonnegativity of the highest-order term yields nonnegativity of all lower-order terms [2210.04361].

The relative-degree interpretation is central. In discrete time, the relative degree of \(h(x_t)\) with respect to \(x_{t+1}=l(x_t,u_t)\) is informally the smallest number of steps after which \(u_t\) directly affects the evolution of a function of \(h\) [2210.04361]. The multi-agent DSMPC formulation states this through derivatives of composed maps: the relative degree of the output \(z_i(t)=g(x_i(t))\) with respect to \(u_i(t)\) is \(m\) if
\[
\frac{\partial}{\partial u_i(t)}\big[g\circ \bar f^{l}(f(x_i(t),u_i(t)))\big]=\mathbf{0}_q,\quad 0\le l\le m-2,
\]
and
\[
\frac{\partial}{\partial u_i(t)}\big[g\circ \bar f^{m-1}(f(x_i(t),u_i(t)))\big]\neq \mathbf{0}_q
\]
[2508.19678]. This is the discrete-time analogue of high relative degree in continuous-time HOCBF theory.

The framework also supports a weaker enforcement principle. The 2022 paper remarks that satisfying the \(i\)-th order CBF constraint is sufficient to render \(\mathcal{C}_1\cap\dots\cap\mathcal{C}_i\) forward invariant, so it is not strictly necessary to go all the way to order \(m\) in an optimal control problem; one can choose a suitable \(m_{\mathrm{CBF}\le m}\) to reduce computational complexity [2210.04361]. This suggests a design trade-off between barrier order and solver burden.

## 3. Relation to standard discrete-time CBFs, sampled-data CBFs, and continuous-time HOCBFs

The relative-degree-one discrete-time CBF condition used in prior MPC-DCBF and NMPC-DCBF formulations is
\[
h(x_{t+1})\ge \omega_t(1-\gamma)\,h(x_t),\quad 0<\gamma\le 1,
\]
with slack variable \(\omega_t\) [2210.04361]. When \(m=1\), the DHOCBF recursion collapses to \(\psi_0=h\), and the highest-order condition applies directly to \(h\), so the relative-degree-one case is recovered [2210.04361]. The distinction is therefore not categorical but structural: DHCBFs lift the safety condition into auxiliary functions so that systems with multi-step input influence can still be constrained through an exponential-type inequality.

This differs from the discrete-time CBF notion that is necessary and sufficient for controlled invariance in a one-step sense,
\[
\sup_{u\in U} h(f(x)+\tilde g(x)u)\ge 0,\quad \forall x\in\mathcal{S},
\]
with safe input set
\[
K_{\mathcal{S}}(x)=\{u\in U: h(f(x)+\tilde g(x)u)\ge 0\}
\]
[2004.01858]. That formulation is explicitly least restrictive among several discrete-time barrier inequalities discussed there, because prior relaxed forms imply it but not conversely [2004.01858]. By contrast, DHCBFs deliberately introduce higher-order auxiliary conditions to account for relative degree rather than to relax one-step invariance.

Sampled-data CBF work provides a parallel discrete-time viewpoint even when the plant evolves continuously between samples. In that setting, a sampled-data barrier condition constrains the increment of \(s_h\) over one sample,
\[
s_h(F_h(x,u))-s_h(x)\ge -h\,\alpha(s_h(x)),
\]
and practical safety is established for the exact sampled-data system when the approximate discrete-time model is one-step consistent with the exact flow [2203.11470]. Although that framework does not explicitly define DHCBFs, it provides a direct bridge from sampled implementation to discrete-time barrier design and shows how approximate discrete-time models and Runge–Kutta maps can support convex synthesis for higher-relative-degree structures [2203.11470].

Continuous-time HOCBFs instead use derivatives and Lie derivatives,
\[
\zeta_i(x)=\dot\zeta_{i-1}(x)+\alpha_i(\zeta_{i-1}(x)),
\]
and enforce a highest-order inequality involving \(L_f^r h\) and \(L_gL_f^{r-1}h\) [2403.09573]. Discrete-time DHCBFs replace derivatives with finite differences and use inequalities like
\[
\psi_i(x_t)=\psi_{i-1}(x_{t+1})-\psi_{i-1}(x_t)+\gamma_i\psi_{i-1}(x_t)
\]
or the equivalent multi-step Taylor-based conditions discussed below [2210.04361; 2503.15014].

## 4. Embedding in model predictive control

A principal application of DHCBFs is safety-critical MPC. The baseline nonlinear MPC with a relative-degree-one discrete-time CBF minimizes
\[
p(x_{t+N})+\sum_{k=0}^{N-1} q(x_{t+k},u_{t+k},\omega_{t+k})
\]
subject to dynamics, input/state constraints, and
\[
h(x_{t+k+1})\ge \omega_{t+k}(1-\gamma)\,h(x_{t+k})
\]
over the horizon [2210.04361]. The DHOCBF extension preserves the MPC structure but replaces the one-step barrier with high-order recursion on a candidate barrier \(\psi_0\) and its auxiliary functions \(\psi_i\) [2210.04361].

In the iterative MPC-DCBF formulation, the dynamics and barrier functions are linearized at each time step. The nonlinear dynamics are approximated by
\[
x_{t+1}^{j+1}-x_{t+1}^{j}=A_j(x_t^{j+1}-x_t^j)+B_j(u_t^{j+1}-u_t^j),
\]
with
\[
A_j=D_x l(x_t^j,u_t^j),\quad B_j=D_u l(x_t^j,u_t^j)
\]
[2210.04361]. The nonlinear safety function is replaced by a linearized obstacle boundary \(h_e(x_t\mid x_{s,t})\), obtained from the tangent of the obstacle boundary at an intersection point \(x_{s,t}\), and this linearized function becomes \(\psi_0\) in the high-order recursion [2210.04361].

To improve feasibility, per-order, per-time slack variables \(\omega_{i,t}^j\) are introduced:
\[
\psi_{i-1}(x_{t+1}^j)\ge \omega_{i,t}^j(1-\gamma_i)\psi_{i-1}(x_t^j),\quad \omega_{i,t}^j\in\mathbb{R}
\]
[2210.04361]. Because this form can remain nonconvex, the paper derives a linear convex reparameterization,
\[
\psi_{i-1}(x_t^j)+\sum_{v=1}^i Z_{v,i}(1-\gamma_i)^t\psi_0(x_v^j)\ge \omega_{i,t}^j Z_{0,i}(1-\gamma_i)^t\psi_0(x_0^j),
\]
where \(Z_{0,i},Z_{v,i}\) depend only on \(\gamma_1,\dots,\gamma_i\) and \(Z_{i,i}=0\) [2210.04361]. This yields linear and convex DHOCBF constraints over the MPC horizon.

A distributed variant for nonlinear multi-agent systems uses DHCBFs inside distributed safety-critical MPC alongside discrete-time control Lyapunov functions (DCLFs). Each agent enforces
\[
\psi_m(x_i(k|t),u_i(k|t))\ge 0
\]
as a hard MPC constraint, with neighbor states replaced by estimated neighbor trajectories to enable distributed implementation [2508.19678]. The resulting DSMPC formulation includes dynamics, box constraints, DHCBF constraints, compatibility bounds on estimation errors, and DCLF-based terminal constraints [2508.19678].

## 5. Computational structure, convexity, and performance

The computational challenge of DHCBFs is explicit in the literature: nonlinear dynamics and nonlinear barrier functions render MPC nonconvex, while higher-order barriers introduce more auxiliary functions and more constraints across larger horizons [2210.04361]. The iterative convexification strategy addresses this by solving a sequence of convex subproblems based on local linearization [2210.04361].

The reported numerical benchmark compares NMPC-DCBF and iMPC-DCBF for relative degree \(m_{\mathrm{CBF}}=1\) and \(2\), with horizons \(N=4,8,12,16,20,24\) [2210.04361]. The key observations are summarized below.

| Aspect | Observation | Context |
|---|---|---|
| Computation time | iMPC-DCBF achieves about 10x faster computation times | Compared with NMPC-DCBF |
| Infeasibility rate | iMPC-DCBF has about 30% more infeasibility | Trade-off with speed |
| Relative degree handling | Both \(m_{\mathrm{CBF}}=1\) and \(2\) are evaluated | High-order handled in iMPC |

These results indicate a speed–feasibility trade-off rather than a universal dominance relation [2210.04361]. In distributed MPC, simulation results similarly report improved performance and reduced computation time relative to MPC with distance constraints, while retaining obstacle avoidance and formation control feasibility under a mild assumption [2508.19678]. This suggests that DHCBFs can reduce the need for longer horizons by imposing more informative safety structure at each step.

A different computational route appears in sampled-data discrete-time CBF synthesis using approximate discrete-time models. For systems with block-integrator structure, if the barrier depends on the first \(q\) blocks and is concave in its last argument, then choosing Runge–Kutta order \(p=\gamma-q+1\) makes the function
\[
\phi_h(x,u):=-s_h(F_h^{a,p}(x,u))+s_h(x)-h\,\alpha(s_h(x))
\]
convex in \(u\) [2203.11470]. This result is not stated as a DHCBF theorem, but it supplies a convexity template for discrete-time higher-relative-degree barriers.

## 6. Alternative formulations, robustness, and verification

Several adjacent directions broaden the DHCBF landscape.

A truncated Taylor formulation proposes a discrete-time high-order condition by approximating the one-step barrier increment with derivatives up to the relative degree \(r\). For high-order safety constraints, the Truncated Taylor CBF (TTCBF) enforces
\[
\sum_{i=1}^{r-1}\frac{(\Delta T)^i}{i!}h^{(i)}(\bm x(t_k),t_k)+\frac{(\Delta T)^r}{r!}h^{(r)}(\bm x(t_k),\bm u(t_k),t_k)+\alpha(h(\bm x(t_k),t_k))+\inf_{\bm u(t_k)\in\mathcal U}R_{\mathrm T}\ge 0
\]
and requires only one class-\(\mathcal K\) function, independent of the relative degree [2601.15196]. The earlier continuous-time-to-discrete-time exposition of TTCBF presents the same design motivation: standard HOCBFs require multiple class-\(\mathcal K\) functions whose tuning scales with relative degree, whereas the truncated Taylor approach reduces design complexity while ensuring safety in sampled-data settings [2503.15014]. This suggests a distinct DHCBF interpretation: rather than constructing auxiliary forward-difference sequences, one can enforce an approximate \(r\)-step discrete-time CBF condition using truncated Taylor expansions.

Robust sampled-data CBF methods with zero-order hold and input delay offer complementary tools. They derive CBF constraints affine in the input by evaluating barrier inequalities at predicted delayed states and robustifying them over reachable sets or uncertainty sets [2005.06418]. The paper explicitly notes that it does not define DHCBFs, but it provides a robust sampled-data barrier framework and delay-handling methodology directly useful for discrete-time high-order barrier design under ZOH, delays, and uncertainty [2005.06418].

Verification has also become a distinct topic. An optimization-based branch-and-bound method verifies or falsifies candidate discrete-time CBFs by solving global nonconvex optimization problems over the state set and admissible inputs [2409.15448]. Although it targets first-order DTCBFs, the paper states that the same machinery extends conceptually to DHCBFs by replacing the single-step expression \(h(f(x,u))-h(x)+\gamma(h(x))\) with multi-step barrier expressions involving iterated dynamics and control sequences [2409.15448]. This suggests that high-order barrier verification is an optimization problem over more complicated but structurally similar inequalities.

Robust discrete-time CBF synthesis under bounded disturbances uses a sufficient condition
\[
h(f(x,u))-h(x)+\gamma(h(x))\ge L_h w_{\max}
\]
to replace infinitely many disturbance constraints by a single inequality involving a Lipschitz constant and a disturbance bound [2506.13011]. That work does not explicitly define DHCBFs, but it provides a counterexample-guided synthesis and verification loop for robust discrete-time barrier functions, which a plausible extension would carry over to multi-step or higher-order barrier conditions.

Sample-efficient certification of discrete-time CBFs via Lipschitz arguments similarly focuses on first-order conditions,
\[
r(x):=h(f(x))-h(x)\le -\bar\alpha h(x),
\]
and certifies them on the whole zero-sublevel set using \(\epsilon\)-nets and level-wise segmentation [2509.03899]. The paper explicitly notes that the method can be generalized to multi-step maps \(f^m\) and higher-order discrete-time barrier conditions involving forward differences and multi-step compositions [2509.03899]. This suggests a verification pipeline for DHCBFs based on multi-step Lipschitz constants.

## 7. Related constructions and open directions

A distinct branch concerns synthesis rather than online enforcement. For control-affine polynomial discrete-time systems with input constraints and semi-algebraic safe sets, an alternating-descent sum-of-squares framework synthesizes quadratic DTCBFs and polynomial control policies, then extends to higher-degree polynomial DTCBFs [2504.19330]. That work does not use “high-order” in the relative-degree sense, but it directly addresses higher-degree polynomial barrier construction in discrete time and is relevant whenever DHCBF implementations need expressive polynomial certificates [2504.19330].

Discrete-time CBF compositions also matter. A tractable formulation for partially control-affine discrete-time systems introduces nonlinear DT-CBFs that remain affine in control and develops mixed-integer constraints for conjunction, disjunction, implication, exclusive OR, equivalence, and piecewise barrier functions [2004.01858]. This provides a route to expressing mode-dependent safety logic, which a plausible implication is especially useful when DHCBFs must switch between obstacle-avoidance modes or region-specific higher-order constraints.

Beyond model-based constructions, a continuous-time model-free approach to higher-order CBFs for systems of arbitrary relative degree uses funnel control to construct vector-valued auxiliary variables \(\beta_k\) and a first-order zeroing CBF
\[
b(t,\xi):=\tfrac12(1-\|\beta_r(t,\xi)\|^2)
\]
for the highest-order auxiliary variable [2606.04790]. The paper explicitly demonstrates that a straightforward extension of relative-degree-one model-free CBF design to higher-order systems fails, and replaces scalar HOCBF recursion with a funnel-based recursive structure [2606.04790]. While the treatment is continuous-time, this suggests that discrete-time higher-order barrier design need not always mirror the standard HOCBF recursion; alternative recursive transforms may be viable.

Across the literature, the recurring issues are clear. DHCBFs are introduced to handle arbitrary relative degree in discrete time [2210.04361; 2508.19678], but this raises computational complexity, tuning complexity, and feasibility concerns [2210.04361; 2503.15014]. Sampled-data formulations emphasize the gap between continuous-time theory and digital implementation [2005.06418; 2203.11470]. Verification and synthesis remain nontrivial global optimization problems [2409.15448; 2504.19330; 2506.13011; 2509.03899].

This suggests three broad research directions already visible in the cited works. One is simplification of high-order design, exemplified by TTCBFs requiring only one class-\(\mathcal K\) function [2503.15014; 2601.15196]. A second is tractable synthesis and verification, through convexification, branch-and-bound, or sum-of-squares methods [2210.04361; 2409.15448; 2504.19330]. A third is robustness to sampled implementation artifacts, delays, uncertainty, and distributed coupling [2005.06418; 2203.11470; 2508.19678].

In that sense, DHCBFs are best understood not as a single formula but as a family of discrete-time safety constructions for higher-relative-degree constraints. Their common objective is the same: to make the safe set, or an auxiliary hierarchy of safe sets, forward invariant under discrete-time dynamics when the control acts on the safety output only through multi-step evolution.

Source: https://www.emergentmind.com/topics/discrete-time-high-order-control-barrier-functions-dhcbfs