## Characterization of the Dynamical Properties of Safety Filters for Linear Planar Systems

Yiting Chen ∗ Pol Mestres ∗ Emiliano Dall'Anese Jorge Cort´ es

Abstract -This paper studies the dynamical properties of closed-loop systems obtained from control barrier functionbased safety filters. We provide a sufficient and necessary condition for the existence of undesirable equilibria and show that the Jacobian matrix of the closed-loop system evaluated at an undesirable equilibrium always has a nonpositive eigenvalue. In the special case of linear planar systems and ellipsoidal obstacles, we give a complete characterization of the dynamical properties of the corresponding closed-loop system. We show that for underactuated systems, the safety filter always introduces a single undesirable equilibrium, which is a saddlepoint. We prove that all trajectories outside the global stable manifold of such equilibrium converge to the origin. In the fully actuated case, we discuss how the choice of nominal controller affects the stability properties of the closed-loop system. Various simulations illustrate our results.

## I. INTRODUCTION

Modern autonomous systems and cyber-physical systems - from self-driving vehicles and robotic systems to critical infrastructures - must provide safety guarantees while performing complex operational tasks [1]. A popular approach to promote safety, where the term 'safety' here refers to the ability to render a predefined set of states forward invariant, relies on the so-called safety filters ; these filters take a potentially unsafe nominal controller, designed to provide stability or optimality guarantees, and minimally modify it to account for safety constraints. While the filtered controller ensures safety, it may not preserve the stability or optimality properties of the nominal controller. This challenge is the main motivation for this work.

Literature Review: One of the main approaches for rendering a given set forward invariant is via Control Barrier Functions (CBFs) [2]-[5]. Given a nominal controller with desirable properties such as asymptotic stability of an equilibrium, CBFs acts on top of the nominal controller to ensure safety. This technique is often referred to as a safety filter [6]. The main research question here is whether the closed-loop system with safety filters retains the stability guarantees of the nominal controller. This was studied in, e.g., [7], which provides an estimate of the region of attraction of the equilibrium. However, it is unclear how conservative such estimate may be for general systems. The seminal works in [8]-[12] show that designs similar to safety filters can introduce undesirable equilibria that may be stable or unstable.

∗ Equal contribution of the authors. Y. Chen and E. Dall'Anese are with the Department of Electrical and Computer Engineering at Boston University; P. Mestres and J. Cort´ es are with the Department of Mechanical and Aerospace Engineering at the University of California San Diego.

This work was supported by the AFOSR Award FA9550-23-1-0740.

Statement of Contributions: The goal of this paper is to advance the understanding of the dynamical properties of closed-loop systems obtained from CBF-based safety filters. The main contribution of the paper is two-fold:

(i) Our first contribution is to characterize the undesirable equilibria that emerge in the closed-loop system formed by a control-affine dynamical system, a stabilizing nominal controller, and a CBF-based safety filter. General obstacles are considered (this is the subject of Section III).

(ii) Next, we focus our attention to linear time-invariant (LTI) planar systems (Section IV). We show that, for these systems, the dynamical properties of systems with ellipsoidal obstacles are equivalent to those of systems with circular obstacles. For underactuated LTI planar systems, we give a complete characterization of the trajectories of the closedloop system. We show that such systems always have a single undesirable equilibrium. Moreover, we show that such undesirable equilibrium is a saddle point and show that all trajectories that lie outside the global stable manifold of this equilibrium converge to the origin. For fully actuated LTI planar systems, we show that the closed-loop system can have up to three undesirable equilibria, and characterize their stability properties.

Additionally, we show that in the fully actuated case there always exists a nominal controller (which can be explicitly computed) that makes the closed-loop system have a single undesirable saddle point equilibrium. Therefore, our findings can be used to inform the design of the nominal controller. For reasons of space, proofs are included in an extended version [13].

## II. PRELIMINARIES AND PROBLEM STATEMENT

Notation . We denote by N &gt; 0 and R the set of positive integers, real, and nonnegative numbers. We use bold symbols to represent vectors and non-bold symbols to represent scalar quantities; 0 n represents the n -dimensional zero vector. Given x ∈ R n , ∥ x ∥ denotes its Euclidean norm. Given a matrix G ∈ R n × n , ∥ x ∥ G = √ x T G x . A function β : R → R is of extended class K ∞ if β (0) = 0 , β is strictly increasing and lim s →±∞ β s ( ) = ±∞ . Given a set S ⊂ R n , we denote by Int ( S ) and ∂S the interior and boundary of S , respectively. For a continuously differentiable function h : R n → R , ∇ h ( x ) denotes its gradient at x .

Consider the system ˙ x = f ( x ) , with f : R n → R n locally Lipschitz. Then, for any initial condition x 0 ∈ R n at time t 0 , there exists a maximal interval of existence [ t 0 , t 1 ) such that x ( ; t x 0 ) is the unique solution to ˙ x = f ( x ) on [ t 0 , t 1 ) , cf. [14]. For f continuously differentiable and x ∗ an equilibrium point of f (i.e., f ( x ∗ ) = 0 n ), x ∗ is degenerate if the Jacobian of f evaluated at x ∗ has at least one eigenvalue with real part equal to zero (otherwise, we refer to x ∗ as hyperbolic ). Given a hyperbolic equilibrium point with k ∈ Z &gt; 0 eigenvalues with negative real part, the Stable Manifold Theorem [15, Section 2.7] ensures that there exists an invariant k -dimensional manifold S for which all trajectories with initial conditions lying on S converge to x ∗ . The global stable manifold at x ∗ is defined as W s ( x ∗ ) =

⋃ { ≤ t 0 , x 0 ∈ } S x ( ; t x 0 ) . Given a complex number z ∈ C , Re ( z ) denotes its real part.

## A. Control barrier functions and safety filters

Consider a control-affine dynamical system of the form

$$\dot { x } = f ( { \mathbf x } ) + g ( { \mathbf x } ) { \mathbf u }, \quad \quad \ ( 1 ) \quad \text{are}$$

where f : R n → R n and g : R n → R n × m are locally Lipschitz functions, x ∈ R n is the state, and u ∈ R m is the input.

Definition 1 (Control Barrier Function): Let h : R n → R be a continuously differentiable function, and define the set C = { x ∈ R n : h ( x ) ≥ 0 } . The function h is a CBF of C for the system (1) if there exists an extended class K ∞ function α such that, for all x ∈ C , there exists u ∈ R m satisfying ∇ h ( x ) ⊤ ( f ( x ) + g ( x u ) ) + α h ( ( x )) ≥ 0 . □

Suppose that a nominal controller u = k ( x ) is designed so that the system ˙ x = ˜ ( f x ) := f ( x ) + g ( x ) k ( x ) renders the origin globally asymptotically stable (this is without loss of generality). Consider the system

$$\dot { \mathbf x } = f ( \mathbf x ) + g ( \mathbf x ) v ( \mathbf x ), \quad \quad ( 2 ) \quad.$$

↦

where the map x → v ( x ) is defined as:

$$v ( \mathbf x ) = \arg \min _ { \theta \in \mathbb { R } ^ { m } } \| \theta \| _ { G ( \mathbf x ) } ^ { 2 } & & ( 3 ) \quad \text{sufficie} \\ r \, \dots.$$

$$\text{s.t.} \, \nabla h ( \mathbf x ) ^ { \top } ( f ( \mathbf x ) + g ( \mathbf x ) ( k ( \mathbf x ) + \theta ) ) + \alpha ( h ( \mathbf x ) ) \geq 0 \quad _ { \mathbf A \, \text{e} \, } ^ { \text{Le} \, }$$

with G : R n → R m m × continuously differentiable and positive definite for all x ∈ R n . We assume the following.

Assumption 1 (Origin in the interior of C ): The set { x ∈ R n : h ( x ) = 0 , ˜ ( f x ) = 0 n } is empty and h ( 0 n ) &gt; 0 . □

̸

Assumption 2 (Feasibility): There exists an extended class K ∞ function α such that g ( x ) ⊤ ∇ h ( x ) = 0 for all x in { x ∈ R n : h x ( ) ≥ 0 , ∇ h ( x ) ⊤ f ( x ) + α h ( ( x )) ≤ } 0 . □

̸

Assumption 2 ensures that (3) is feasible for all x ∈ R n and therefore v ( x ) is well-defined for all x ∈ R n . Moreover, under Assumption 2, and using arguments similar to [16, Lemma III.2], one can show that v ( x ) is locally Lipschitz. Assumption 2 also ensures that ∂h ∂ x ( x ) = 0 n for all x ∈ ∂ C . From [3, Thm. 2], it follows that the system (2) with the controller v ( x ) renders the set C forward invariant. Because of this feature, and because C is modeling a set of safe states, v ( x ) is typically referred to as safety filter .

## B. Problem Statement

We consider a control-affine dynamical system as in (1) and a safe set C ⊂ R n defined as the 0 -superlevel set of a differentiable function h : R n → R . Assume that h is a CBF of C for system (1), and that Assumptions 1 and 2 hold. Studying the dynamical behavior of (2) is challenging. Indeed, as noted in [7], it does not inherit the global asymptotic stability properties of the controller k , and can even have undesirable equilibria [8]-[11]. However, most of these works focus on studying conditions under which such undesirable equilibria exist or can be confined to specific regions of interest, but do not study dynamical properties of the closed-loop system. Hence the goal of this paper is as follows:

Problem 1: Given system (1) with a stabilizing nominal controller k ( x ) and the safety filter v ( x ) , characterize the dynamical properties of (2) (such as undesirable equilibria and their regions of attraction, limit cycles and region of attraction of the origin) and investigate how these properties are determined by the original closed-loop system ˙ x = f ( x ) + g ( x ) k ( x ) . □

In the following section, we consider the system (2) and characterize its undesirable equilibria. In Section IV, given the complexity of solving Problem 1, we then restrict our attention to linear planar systems.

## III. CHARACTERIZATION OF UNDESIRABLE EQUILIBRIA

We start by reformulating the expression for the unique optimal solution v ( x ) of the quadratic program (3). Let η ( x ) = ∇ h ( x ) T ( f ( x ) + g ( x ) k ( x )) + α h ( ( x )) . Then,

$$v ( { \mathbf x } ) = \begin{cases} { \mathbf 0 } _ { m }, & \text{if $\eta(\mathbf x) \geq 0$}, \\ \bar { u } ( { \mathbf x } ), & \text{if $\eta(\mathbf x) < 0$}, \end{cases}$$

where ¯( u x ) := -η ( x ) G ( x ) -1 g ( x ) ⊤ ∇ h ( x ) ∥ g ( x ) ⊤ ∇ h ( x ) ∥ 2 G ( x ) -1 . We use this expression in the following result, which provides a necessary and sufficient condition for undesirable equilibria of (2).

̸

Lemma 1: (Conditions for undesirable equilibria): Let Assumptions 1 and 2 be satisfied. Let p 0 ∈ R n be such that ˜ ( f p 0 ) = 0 n . Then, p 0 is an equilibrium of (2) if and only if there exists δ &lt; 0 such that

$$h ( \mathbf p _ { 0 } ) = 0 \text{ and }$$

$$\tilde { f } ( { \mathbf p } _ { 0 } ) = \delta g ( { \mathbf p } _ { 0 } ) G ( { \mathbf p } _ { 0 } ) ^ { - 1 } g ( { \mathbf p } _ { 0 } ) ^ { \top } \nabla h ( { \mathbf p } _ { 0 } ) \,. \quad \ \square$$

This result has the same flavor as [9, Theorem 2] and [10, Proposition 5.1], which characterize the undesirable equilibria for related, but different, safety filter designs.

By Lemma 1, we can define the set of potential undesirable equilibria of (2) as:

$$\mathcal { E } \coloneqq \{ \mathbf x \, \colon \, \exists \, \delta \in \mathbb { R } \ s. t. \, ( \mathbf x, \delta ) \text{ solves } ( 5 ) \} \,.$$

On the other hand, the set of undesirable equilibria is:

$$\hat { \mathcal { E } } \coloneqq \{ \boldsymbol x \, \colon \, \exists \, \delta < 0 \text{ s.t. } ( \boldsymbol x, \delta ) \text{ solves } ( 5 ) \} \subset \mathcal { E }.$$

The term undesirable stems from the fact that these equilibria are different from the origin, which is the equilibrium point where the system needs to be stabilized. By Lemma 1, it follows that determining the equilibrium points of system (2) is equivalent to solving (5) and checking the sign of δ . For a solution ( p 0 , δ p 0 ) to (5), we refer to δ p 0 as the indicator of p 0 , since the sign of δ p 0 determines whether p 0 is a new, undesirable , equilibrium of the system with the CBF filter. Additionally, we show that the value of the indicator is useful for determining the stability properties of the undesirable equilibrium. For a given undesirable equilibrium p 0 of (2), the indicator can be computed as δ p 0 = ∇ h ( p 0 ) ⊤ ˜ ( f p 0 ) ∥ g ( p 0 ) ⊤ ∇ h ( p 0 ) ∥ 2 G ( x ) -1 . In addition, Assumption 1 ensures that no solution of (5) has δ = 0 .

Under appropriate conditions, the next result shows that we can compute the Jacobian of ˜ ( f x ) + g ( x ) v ( x ) at x ∈ E ˆ and find one of its eigenvalues.

Lemma 2: (Jacobian at the undesirable equilibrium): Let Assumptions 1 and 2 be satisfied and assume that D = g ( x ) G ( x ) -1 g ( x ) ⊤ is a constant matrix, ˜ ( f x ) , α ( ) · are differentiable and h ( x ) is twice differentiable. For any x ∈ E ˆ , the Jacobian of ˜ ( f x ) + g ( x ) v ( x ) evaluated at x is

$$& J \, | _ { \mathbf x \in \hat { \mathcal { E } } } = J _ { \tilde { f } } - \frac { \ell \, \lambda \, \nu \, \cdots \nu \, \cdots \nu } { \nabla h ( \mathbf x ) ^ { \top } D \nabla h ( \mathbf x ) } [ J _ { \tilde { f } } + \alpha ^ { \prime } ( 0 ) \mathbf I _ { n } ] \\ & - \frac { D } { \nabla h ( \mathbf x ) ^ { \top } D \nabla h ( \mathbf x ) } [ H _ { h } \nabla h ( \mathbf x ) ^ { \top } \tilde { f } ( \mathbf x ) - \nabla h ( \mathbf x ) \tilde { f } ( \mathbf x ) ^ { \top } H _ { h } ], \quad \text{where} }$$

where J ˜ f is the Jacobian matrix of ˜ ( f x ) and H h is the Hessian of h ( x ) . Moreover, for any x ∈ E ˆ , it holds that

$$( J \, | _ { \mathbf x \in \mathcal { E } } ) ^ { \top } \nabla h ( \mathbf x ) = - \alpha ^ { \prime } ( 0 ) \nabla h ( \mathbf x ),$$

the algebraic multiplicity of -α ′ (0) is 1 , and all the other eigenvalues of J | x ∈E ˆ do not change when α ( ) · changes. □

The proof of Lemma 2 follows from a careful computation. Note that J always has an eigenvalue -α ′ (0) ; it follows that all the undesirable equilibria are degenerate if α ′ (0) = 0 , which complicates the stability analysis. If α ′ (0) &gt; 0 , the Jacobian evaluated at x ∈ E ˆ always has a negative eigenvalue. Lemmas 1 and 2 show that the extended K ∞ function α ( ) · does not play a role in the existence of undesirable equilibria. Additionally, changing α ( ) · will only affect one eigenvalue of the Jacobian evaluated at x ∈ E ˆ . The assumption that g ( x ) G ( x ) -1 g ( x ) ⊤ is constant is satisfied for several classes of systems, such as mechanical systems, like the ones considered in [7, Section III.B].

## IV. LTI PLANAR SYSTEMS WITH SAFETY FILTERS

Since Problem 1 is difficult to solve in general, here we provide a solution for it for planar LTI dynamics and ellipsoidal obstacles. Consider the LTI planar system

$$\dot { \mathbf x } = A \mathbf x + B \mathbf u, \quad \ \ \ ( 6 ) \quad \mathbf x ^ { + }$$

with x = [ x , x 1 2 ] ⊤ ∈ R 2 , u ∈ R m , with m ∈ { 1 2 , } , A ∈ R 2 × 2 , and with B ∈ R 2 × m full column rank. We make the following assumption on (6).

Assumption 3 (Stabilizability): The system (6) is stabilizable. Moreover, let u = -K x , K ∈ R 2 × m , be any stabilizing controller such that ˜ A = A -BK is Hurwitz. □

In this setup, the system (2) is then customized as follows:

$$\dot { \mathbf x } = F ( \mathbf x ) \coloneqq ( A - B K ) \mathbf x + B v ( \mathbf x ), \quad \ ( 7 ) \quad \text{follow}$$

where the safety filter is given by

$$v ( { \mathbf x } ) = \begin{cases} 0, & \text{if $\eta(\mathbf x) \geq 0$}, \\ - \frac { \eta ( { \mathbf x } ) G ( { \mathbf x } ) ^ { - 1 } B ^ { T } \nabla h ( { \mathbf x } ) } { \| B ^ { T } \nabla h ( { \mathbf x } ) \| _ { G ( { \mathbf x } ) ^ { - 1 } } ^ { 2 } }, & \text{if $\eta(\mathbf x) < 0$}. \end{cases} \quad ( 8 )$$

In the following, we show that the undesirable equilibria and their stability properties of (7) with ellipsoidal obstacles are equivalent to those of a system with circular obstacles.

Proposition 1: (Safety filters with ellipsoidal and circular obstacles have the same dynamical properties): Let x c ∈ R 2 , P ∈ R 2 × 2 positive definite, h ( x ) = ( x -x c ) T P ( x -x c ) -1 , C := { x ∈ R n : h ( x ) ≥ 0 } . Suppose that P = E E T , with E ∈ R 2 × 2 also positive definite, and define ˆ x c = E x c , ˆ (ˆ) h x = (ˆ x -ˆ ) x c T (ˆ x -ˆ ) x c -1 and ˆ C = { x ∈ R n : ˆ ( h x ) ≥ 0 } . Moreover, let ˆ = A EAE -1 , ˆ = B EB G , ˆ (ˆ) = x G E ( -1 ˆ) x and ˆ(ˆ) = η x ∇ ˆ (ˆ) h x T ( ˆ A -ˆ BKE -1 )ˆ + x α h ( ˆ (ˆ)) x . Consider the system

$$\dot { \hat { x } } = \hat { F } ( \hat { x } ) \coloneqq ( \hat { A } - \hat { B } K E ^ { - 1 } ) \hat { x } + \hat { B } \hat { v } ( \hat { x } ), \quad ( 9 )$$

where

$$& \tilde { \ } n _ { 1 }, \\ & \text{the} \quad & \hat { v } ( \hat { x } ) = \begin{cases} 0, & \text{if $\hat{\eta}(\hat{\boldsymbol x}) \geq 0$,} \\ - \frac { \hat { \eta } ( \hat { \boldsymbol x } ) \hat { G } ( \hat { \boldsymbol x } ) ^ { - 1 } ( \hat { \boldsymbol x } ) \hat { B } ^ { T } \nabla \hat { h } ( \hat { \boldsymbol x } ) } { \| \hat { B } ^ { T } \nabla \hat { h } ( \hat { \boldsymbol x } ) \| _ { \hat { G } ( \hat { \boldsymbol x } ) ^ { - 1 } } ^ { 2 } }, & \text{if $\hat{\eta}(\hat{\boldsymbol x}) < 0$} \end{cases} ( 1 0 )$$

Then,

- i) ˆ C is forward invariant under system (9) and C is forward invariant under system (7);
- ii) systems (9) and (7) are locally Lipschitz;
- iii) ( A,B ) is stabilizable if and only if ( ˆ A,B ˆ ) is stabilizable;
- iv) ˆ p ∈ R 2 is an undesirable equilibrium of (9) if and only if p := E -1 ˆ p is an undesirable equilibrium of (7);
- v) the Jacobian of ˆ F at ˆ p and the Jacobian of F at p are similar. □

Given that Proposition 1 ensures that undesirable equilibria for general ellipsoidal obstacles have the same stability properties as undesirable equilibria for circular obstacles, in the following we focus on studying the dynamical properties of safety filters for LTI systems and circular obstacles.

Accordingly, we consider the circular unsafe set:

$$\mathcal { C } = \{ \mathbf x \in \mathbb { R } ^ { 2 } \ \colon \, h ( \mathbf x ) = \| \mathbf x - \mathbf x _ { c } \| ^ { 2 } - r ^ { 2 } \geq 0 \},$$

̸

with x c ∈ R 2 the center. We take the extended class K ∞ function in Definition 1 to be linear and with slope α 0 &gt; 0 . We denote the eigenvalues of ˜ A as λ 1 , λ 2 ∈ C 2 . Let V ( x ) = x ⊤ Q x be the associated Lyapunov function, with a positive definite symmetric matrix Q , such that x ⊤ QA ˜ x &lt; 0 for all x = 0 2 . Additionally, we pick G ( x ) = B B ⊤ .

The first result rules out the existence of limit cycles.

Proposition 2: (Non-existence of limit cycles): Suppose that Assumptions 1-3 hold for the closed-loop system (7). Assume that for (7), ˆ E = { ˆ x ∗ } , with ˆ x ∗ a saddle point. Then, there exist α ∗ 1 &gt; 0 such that for any α s ( ) = α s 0 with α 0 ≥ α ∗ 1 , (7) does not have limit cycles in C . □

By combining the results in this section, we have the following.

Theorem 1 (Global behavior analysis): Suppose that the Assumptions 1-3 hold for the closed-loop system (7). Assume that ˆ = E { ˆ x ∗ } and ˆ x ∗ is a saddle point. Then, there exists α ∗ 2 &gt; 0 such that for any α s ( ) = α s 0 with α 0 ≥ α ∗ 2 , if W s (ˆ x ∗ ) denotes the global stable manifold of ˆ x ∗ it holds that:

- 1) if x 0 ∈ W s (ˆ x ∗ ) , then lim t →∞ x ( ; t x 0 ) = ˆ x ∗ ;

$$2 ) \text{ if } x _ { 0 } \not \in W _ { s } ( \hat { x } ^ { * } ), \text{ then } \overset { \iota \to \infty } { \underset { t \to \infty } { \lim } } x ( t ; x _ { 0 } ) = 0 _ { 2 }. \quad \quad \Box \quad \text{ stabiliz}$$

Remark 1 (Almost global asymptotic stability): The Stable Manifold Theorem [15, Ch. 2.7] ensures that if ˆ x ∗ is a saddle point in R 2 , the local stable manifold is 1 -dimensional. Therefore, it has measure of zero. Moreover, the global stable manifold must also have measure of zero. If this were not the case, solutions would have to intersect. However this is not possible due to the uniqueness of solutions. Hence { x 0 ∈ R n : lim t →∞ x ( ; t x 0 ) = 0 n } = S . It follows that the set of initial conditions whose associated trajectory converges to ˆ x ∗ has measure zero. □

## A. Under-actuated LTI Planar Systems

In the under-actuated case, we write

$$A = \begin{bmatrix} a _ { 1 1 } & a _ { 1 2 } \\ a _ { 2 1 } & a _ { 2 2 } \end{bmatrix}, \ B = \begin{bmatrix} b _ { 1 } \\ b _ { 2 } \end{bmatrix}, \ { \mathbf x } = \begin{bmatrix} x _ { 1 } \\ x _ { 2 } \end{bmatrix}, \quad ( 1 1 ) \quad \text{pictur} \\ \text{$\text{follow}$} \\ C o. }$$

Throughout this section, we denote x c = [ x c, 1 , x c, 2 ] ⊤ and let β = a 11 b 2 -b 1 a 21 , γ = a 22 b 1 -b 2 a 12 , and T 3 = -γx c, 2 + βx c, 1 and assume that k : R 2 → R is a linear stabilizing controller of the form k ( x ) = -K x = -k x 1 1 -k x 1 2 for some k , k 1 2 ∈ R . We note also that since in this case G is a scalar, (7) is independent of G .

The following results give conditions on h and system (11) that ensure that Assumptions 1 and 2 hold.

Lemma 3 (Conditions for Assumption 1): Assumption 1 holds if and only if ∥ x c ∥ 2 &gt; r 2 . □

The proof of Lemma 3 follows from the observation that ∥ x c ∥ 2 &gt; r 2 guarantees that the origin is safe.

Proposition 3 (Conditions for Assumption 2): Let α 0 &gt; 0 , T 1 := b 2 β + b 1 γ + 1 2 α 0 ( b 2 2 + b 2 1 ) , and T 2 := ( βx c, 1 -γx c, 2 ) 2 +2 α r T 0 2 1 . Suppose that r &gt; 0 , b 2 1 + b 2 2 &gt; 0 , T 1 &gt; 0 , and r √ b 2 2 + b 2 1 &gt; | T 3 | + √ T 2 2 T 1 . Then, Assumption 2 holds with the linear extended class K function α s ( ) = α s 0 . □

We next give a result that will be used later in the paper.

̸

Lemma 4 (Conditions for β and γ ): Let Assumption 3 hold, then γ 2 + β 2 &gt; 0 . Furthermore, suppose that the conditions in Proposition 3 hold. Then, r 2 ( γ 2 + β 2 ) -T 2 3 &gt; 0 . Moreover, if Assumption 1 holds, then γx c, 1 + βx c, 2 = 0 .

Next we characterize the undesirable equilibria of the closed-loop system (7) with (11).

Proposition 4: (Equilibria in Under-actuated Systems): Suppose that Assumptions 1, 3 and the conditions in Proposition 3 hold. Define p + := ( γz + , βz + ) , and p -:= ( γz -, βz -) , where

$$z _ { \pm } = \frac { \gamma x _ { c, 1 } + \beta x _ { c, 2 } \pm \sqrt { r ^ { 2 } ( \gamma ^ { 2 } + \beta ^ { 2 } ) - T _ { 3 } ^ { 2 } } } { \gamma ^ { 2 } + \beta ^ { 2 } }.$$

Then,

- i) if γx c, 1 + βx c, 2 &lt; 0 , p + is the only undesirable equilibrium of the closed-loop system (7) with (11);
- ii) if γx c, 1 + βx c, 2 &gt; 0 , p -is the only undesirable equilibrium of the closed-loop system (7) with (11).

̸

Note that by Lemma 4, γx c, 1 + βx c, 2 = 0 . Therefore Proposition 4 shows that for linear, planar, underactuated and stabilizable linear systems, (2) has exactly one undesirable equilibrium. Note also that the result in Proposition 4 is independent of the linear stabilizing controller k and the extended class K function α chosen.

The following result establishes that the undesirable equilibrium of the closed-loop system is always a saddle point.

Proposition 5: (Undesirable Equilibria are Saddle Points): Suppose that Assumptions 1, 3 and the conditions in Proposition 3 hold. Then there always exists one and only one undesirable equilibrium, which is a saddle point.

Note that the results in Propositions 4 and 5 are independent of the choice of weighting matrix G , nominal controller k or extended class K function α . The combination of Propositions 4 and 5 with Theorem 1 provide a complete picture of the under-actuated case, which we summarize as follows.

Corollary 1: (Characterization of trajectories for linear planar underactuated systems): Suppose that Assumptions 1, 3 and the conditions in Proposition 3 hold. Then, the closed-loop system (7) obtained from (11) has one and only one undesirable equilibrium ˆ x ∗ equal to either p + or p -. Additionally, there exists α ∗ 2 &gt; 0 such that for any α s ( ) = α s 0 with α 0 ≥ α ∗ 2 , if W s (ˆ x ∗ ) denotes the global stable manifold of ˆ x ∗ it holds that:

- if x 0 ∈ W s (ˆ x ∗ ) , then lim x ( ; t x 0 ) = ˆ x ∗ ;

$$\mathbf n _ { \sqcap } \, \left ( \frac { 1 } { 2 } \right ) \, \text{ if } x _ { 0 } \notin W _ { s } ( \hat { x } ^ { * } ), \, \text{ then } \sum _ { t = \infty } ^ { s + s \infty } x ( t ; x _ { 0 } ) = 0 _ { 2 }. \quad \right ) \, \right )$$

$$1 ) & \text{ if } x _ { 0 } \in W _ { s } ( \hat { x } ^ { * } ), \, \text{ then } \lim _ { t \to \infty } x ( t ; 2 \\ 2 ) & \text{ if } x _ { 0 } \not \in W _ { s } ( \hat { x } ^ { * } ), \, \text{ then } \lim _ { t \to \infty } x ( t ; 2 \\ n \ r. n. \ e \dots. \ e \ r \tau \ n \dots. \ e \dots.$$

B. Fully Actuated LTI Planar Systems

We now consider the case where B is invertible; in this case, Assumptions 2 and 3 are satisfied.

- 1) x c is an eigenvector of ˜ A : We start by considering two conditions for the case where x c is an eigenvector of ˜ A .

We have that there exists only one undesirable equilibrium and it is a degenerate equilibrium if and only if Condition 1 is true. If Condition 2 is true, there are two undesirable equilibria, one of which is a saddle point and the other one is a degenerate equilibrium.

Condition 1. λ 1 &lt; λ 2 &lt; 0 , ˜ A x c = λ 2 x c , ˜ A v 1 = λ 1 v 1 , v 2 = x c ∥ x c ∥ , 1 -( λ 1 -λ 2 ) 2 r 2 λ 2 2 ∥ x c ∥ 2 = 0 , ( v ⊤ 1 v 2 ) 2 = 1 -( λ 1 -λ 2 ) 2 r 2 λ 2 2 ∥ x c ∥ 2 . Condition 2. λ 1 &lt; λ 2 &lt; 0 , ˜ A x c = λ 2 x c , ˜ A v 1 = λ 1 v 1 , v 2 = x c ∥ x c ∥ , 1 -( λ 1 -λ 2 ) 2 r 2 λ 2 2 ∥ x c ∥ 2 = 0 , ( v ⊤ 1 v 2 ) 2 &gt; 1 -( λ 1 -λ 2 ) 2 r 2 λ 2 2 ∥ x c ∥ 2 .

If neither Condition 1 nor Condition 2 is true, we summarize the results about undesirable equilibria for the case that x c is an eigenvector of ˜ A in Tables I and II. We gather all the cases in the following result.

Proposition 6: (Characterization of undesirable equilibria): Let Assumptions 1 be satisfied and B be invertible. Given that ˜ A is stable and x c is an eigenvector of ˜ A , then one of the following is true:

|                                       |   SP |   DE |   ASE |
|---------------------------------------|------|------|-------|
| ( v ⊤ 1 v 2 ) 2 < 1 - r 2 λ 2 ∥ x 2 c |    1 |    0 |     0 |
| ( v ⊤ 1 v 2 ) 2 = 1 - r 2 λ 2 ∥ x 2 c |    1 |    1 |     0 |
| ( v ⊤ 1 v 2 ) 2 > 1 - r 2 λ 2 ∥ x 2   |    2 |    0 |     1 |

c

## TABLE I

˜ A STABLE, ˜ A v 2 = λ v 2 + v 1 , v 1 = x c ∥ x c ∥ , ˜ A x c = λ x c , ∥ v 2 ∥ = 1 . SP: SADDLE POINT, DE: DEGENERATE EQUILIBRIUM, ASE: UNDESIRABLE ASYMPTOTICALLY STABLE EQUILIBRIUM.

|                                                             |   SP |   DE |   ASE |
|-------------------------------------------------------------|------|------|-------|
| ( v ⊤ i v j ) 2 < 1 - ( λ i - λ j ) 2 r 2 λ 2 i ∥ x c ∥ 2 2 |    1 |    0 |     0 |
| ( v ⊤ i v j ) 2 = 1 - ( λ i - λ j ) 2 r λ 2 i ∥ x c ∥ 2     |    1 |    1 |     0 |
| ( v ⊤ i v j ) 2 > 1 - ( λ i - λ j ) 2 r 2 λ 2 i ∥ x c ∥ 2   |    2 |    0 |     1 |

## TABLE II

˜ A STABLE, ˜ A x c = λ i x c , v i = x c ∥ x c ∥ , ˜ A v j = λ j v j , ∥ v j ∥ = 1 , i, j = { 1 2 , } , { v i , v j } LINEARLY INDEPENDENT.

- (i) |E| = 2 , |E| ˆ = 1 , x ∈ E ˆ is a degenerate equilibrium.
- (ii) |E| = 2 , |E| ˆ = 1 , x ∈ E ˆ is a saddle point.
- (iii) |E| = 3 , |E| ˆ = 2 , one point in ˆ E is a saddle point and the other point in ˆ E is a degenerate equilibrium.
- (iv) |E| = 4 , |E| ˆ = 3 , two points in ˆ E are saddle points and the other point in ˆ E is asymptotically stable. □

Proposition 6 asserts that the number and the stability property of the undesirable equilibria are determined by the number of solutions of (5), if x c is an eigenvector of ˜ A .

Proposition 7: (Spectrum of ˜ A does not determine stability properties of undesirable equilibria): Let Assumption 1 be satisfied and B be invertible. Then for any given negative λ 1 and λ 2 , there exists K 1 and K 2 in the set { K : λ , λ 1 2 = spec ( A -BK ) } , such that there is an undesirable asymptotically stable equilibrium after applying the CBF filter with u = -K 1 x ; and there is only one undesirable equilibrium and it is a saddle point after applying the CBF filter with u = -K 2 x . □

Note that one can characterize the global stability properties of the origin based on the eigenvalues of A -BK . However, based on Proposition 7, the eigenvalues of A -BK do not fully determine the global stability property of the origin. On the other hand, Proposition 7 shows that there always exists a nominal controller u = -K x such that ˜ A has negative eigenvalues and the set of trajectories of (7) that do not converge to the origin has measure zero (cf. Theorem 1). Note that as shown in Lemma 2 and Tables I, II, the class K function only affects the rate of decay in the stable manifold of the undesirable equilibria and it does not affect the existence and stability of undesirable equilibria. Therefore, the choice of nominal controller u = -K x determines in which of the cases we fall into. Ideally, the controller should be designed so that there exists only one undesirable equilibrium and it is a saddle point.

2) x c is not an eigenvector of ˜ A : Next, we analyze the number of undesirable equilibria when x c is not an eigenvector of ˜ A . In this case, the analysis is more involved and we only study the stability properties of undesirable equilibria under some sufficient conditions.

Proposition 8 (Number of undesirable equilibria): Let Assumption 1 be satisfied and B be invertible. Given that G ( x ) = B B ⊤ and ˜ A is stable and x c is not an eigenvector of ˜ A , then 1 ≤ |E| ≤ ˆ 3 and |E \ E| ˆ ≥ 1 . In addition, if λ 1 ≤ λ 2 , there exists x ∈ E ˆ with indicator δ &lt; λ 1 2 . □

Combining Propositions 1, 4, 8 and Table I, II, it follows that applying the CBF filter to a LTI planar system (either under or fully actuated) with a linear stabilizing controller always introduces at least one undesirable equilibrium when the obstacle is ellipsoidal. By [17, Thm. 9.5] and Lemma 2, there exists at least one trajectory converging to the undesirable equilibrium. This result is consistent with [18], which states that given a local Lipschitz dynamical system and a compact unsafe set, if the safe set is forward invariant then there exists at least one trajectory that does not converge to the origin. Theorem 1 ensures that if there is only one undesirable equilibrium and it is a saddle point, then there is only one such trajectory and it corresponds to the global stable manifold of the undesirable equilibrium.

To analyze the stability of undesirable equilibria in the case that x c is not an eigenvector of ˜ A , we need to determine the eigenvalues of J | x ∈E ˆ . By Lemma 2, -α ′ (0) = -α 0 is an eigenvalue of J | x ∈E ˆ . The result in [13, Lemma 5] provides an expression for the other eigenvalue of J | x ∈E ˆ , and by leveraging it, we get the following result.

Proposition 9: (Sufficient conditions for undesirable equilibria): Let Assumption 1 be satisfied and B be invertible. Given that G ( x ) = B B A ⊤ , ˜ is stable with two real eigenvalues λ 1 &lt; λ 2 and x c is not an eigenvector of ˜ A , then there is no undesirable equilibrium with indicator δ ∈ { λ 1 2 , λ 2 2 } . Besides, let v 1 and v 2 be the eigenvectors associated with λ 1 and λ 2 , respectively, and v ⊤ 1 v 2 ≥ 0 , ∥ v 1 ∥ = ∥ v 2 ∥ = 1 ; and then we can write x c = β 1 v 1 + β 2 v 2 . Then, the following holds.

- i) If β 2 1 + β β 1 2 v ⊤ 1 v 2 ≥ 0 , then for any undesirable equilibrium x with indicator δ such that δ &lt; λ 1 2 , x is a saddle point.
- ii) If β β 1 2 v ⊤ 2 v 1 + β 2 2 ≥ 0 , then for any undesirable equilibrium x with indicator δ such that λ 2 2 &lt; δ &lt; 0 , x is asymptotically stable.
- iii) Define F 1 : R → R as:

$$\begin{array} { c c c } A & \cdot & \cdot \\ & F _ { 1 } ( \delta ) \coloneqq - | \lambda _ { 1 } - 2 \delta | ^ { 2 } | \lambda _ { 2 } - 2 \delta | ^ { 2 } r ^ { 2 } + | \lambda _ { 1 } \beta _ { 1 } | ^ { 2 } | \lambda _ { 2 } - 2 \delta | ^ { 2 } \\ & & + | \lambda _ { 2 } \beta _ { 2 } | ^ { 2 } | \lambda _ { 1 } - 2 \delta | ^ { 2 } \\ I, & & + 2 \text{Re} ( \lambda _ { 1 } ^ { * } \beta _ { 1 } ^ { * } \lambda _ { 2 } \beta _ { 2 } ( \lambda _ { 2 } - 2 \delta ) ^ { * } ( \lambda _ { 1 } - 2 \delta ) \nu _ { 1 } ^ { * } \nu _ { 2 } ). \\ & & \text{t} & \text{if the third order nolvnomial } \stackrel { d F _ { 1 } ( \delta ) } { \text{was only one real} } \end{array}$$

If the third order polynomial dF 1 ( δ ) dδ has only one real root 1 and β 2 1 + β β 1 2 v ⊤ 1 v 2 ≥ 0 , then there exists only one undesirable equilibrium and it is a saddle point. □

If | β 1 | ≫ | β 2 | and | v ⊤ 1 v 2 | is small, then the case β 2 1 + β β 1 2 v ⊤ 1 v 2 ≥ 0 is a generalized version of the case in the first row of Table II. If | β 2 | ≫| β 1 | (i.e., x c is 'essentially'

̸

1 For third-order polynomial ax 3 + bx 2 + cx + d , its discriminant is defined as 18 abcd -4 b 3 d + b 2 c 2 -4 ac 3 -27 a d 2 2 . If a = 0 and the discriminant is negative, the third-order polynomial only has one real root.

Fig. 1. Examples of trajectories of an LTI planar system with a safety filter for a circular obstacle; the figures show the vector fields, the undesirable equilibria, and the desired equilibrium (which is the origin). (a): Under-actuated system. (b)-(c)-(d): Fully actuated system, corresponding to the three rows of Table II respectively. In (a) and (b) the undesirable equilibrium is a saddle point. In (c) there is one degenerate equilibrium and one saddle point. In (d) there are three undesirable equilibria, one is asymptotically stable while the others are saddle points.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00958v2-figure-1769ca8b-4ee8-4569-91ed-6e68b66d0ca7.jpg)

eigenvector associated with λ 2 ) and λ 1 ≪ λ 2 , then the case β β 1 2 v ⊤ 2 v 1 + β 2 2 ≥ 0 is a generalized version of the case in the last row of Table II, as 1 -( λ 2 -λ 1 ) 2 r 2 λ 2 2 ∥ x c ∥ 2 &lt; 0 with λ 1 ≪ λ 2 .

## V. NUMERICAL EXPERIMENTS

As a first experiment, we consider the safety set C = { x : ∥ x -(3 , 2) ⊤ ∥ 2 -1 ≥ 0 } and the under-actuated system ˙ x = [ 4 2 1 1 ] x + [ ] 3 1 u with nominal controller u = -[ 3 -2 ] x . Once the CBF-based filter is applied, there is one undesirable equilibrium (2 2) , ⊤ , as guaranteed by Proposition 4. Examples of trajectories of the system with the safety filter are shown in Figure 1(a), along with the vector field, the spurious undesirable equilibrium, and the desirable equilibrium (which is the origin).

In Figures 1(b), (c) and (d), we consider a safety set C = { x : ∥ x -(2 , 0) ⊤ ∥ 2 -1 ≥ 0 } , and the integrator dynamics ˙ x = u as an example of (6).

In Figure 1(b), we show the results for the integrator dynamics with K = [ -5 0 0 -1 ] , G ( x ) = B B ⊤ , and the safety filter with α s ( ) = α s 0 , α 0 = 10 . There is one undesirable equilibrium (3 , 0) ⊤ . We note that for both the setups in Figures 1 (a) and (b), there is only one undesirable equilibrium and it is a saddle point. Only one trajectory converges to the undesirable equilibrium and all other trajectories converge to the origin.

In Figure 1(c), we show the results for (7) with K = [ -3 4 √ 2 0 -1 ] , G ( x ) = B B ⊤ and α s ( ) = α s 0 , α 0 = √

10 . There two undesirable equilibria, which are ( 5 3 , 2 2 3 ) ⊤ (degenerate equilibrium) and (3 , 0) ⊤ (saddle point). Only one trajectory converges to (3 , 0) ⊤ . The measure of the stable set of the degenerate equilibrium is positive (in fact, the measure is + ∞ ), although the degenerate equilibrium is unstable.

In Figure 1(d), we show that results for (7) with K = [ -1 0 0 -5 ] , G ( x ) = B B ⊤ and α s ( ) = α s 0 , α 0 = 10 . √ √

There are three undesirable equilibria: ( 5 2 , 3 2 ) ⊤ , ( 5 2 , -3 2 ) ⊤ and (3 0) , ⊤ ; the last one is asymptotically stable and the first two are saddle points. The two trajectories converging to ( 5 2 , √ 3 2 ) ⊤ , ( 5 2 , -√ 3 2 ) ⊤ and part of the obstacle constitute the boundary of the region of attraction of (3 , 0) ⊤ . Since the examples in Figure 1(b), (c) and (d) all satisfy that x c is an eigenvector of ˜ A , these results are consistent with Proposition 6 (ii), (iii), (iv), respectively.

## REFERENCES

- [1] N. Matni, A. D. Ames, and J. C. Doyle, 'Towards a theory of control architecture: A quantitative framework for layered multi-rate control,' arXiv preprint arXiv:2401.15185 , 2024.
- [2] P. Wieland and F. Allg¨ ower, 'Constructive safety using control barrier functions,' IFAC Proceedings Volumes , vol. 40, no. 12, pp. 462-467, 2007.
- [3] A. D. Ames, S. Coogan, M. Egerstedt, G. Notomista, K. Sreenath, and P. Tabuada, 'Control barrier functions: theory and applications,' in European Control Conference , Naples, Italy, 2019, pp. 3420-3431.
- [4] W. Xiao and C. Belta, 'High-order control barrier functions,' IEEE Transactions on Automatic Control , vol. 67, no. 7, pp. 3655-3662, 2022.
- [5] W. Xiao, C. G. Cassandras, and C. Belta, Safe Autonomy with Control Barrier Functions: Theory and Applications . Springer, 2023.
- [6] L. Wang, A. Ames, and M. Egerstedt, 'Safety barrier certificates for collisions-free multirobot systems,' IEEE Transactions on Robotics , vol. 33, no. 3, pp. 661-674, 2017.
- [7] W. S. Cortez and D. V. Dimarogonas, 'On compatibility and region of attraction for safe, stabilizing control laws,' IEEE Transactions on Automatic Control , vol. 67, no. 9, pp. 7706-7712, 2022.
- [8] M. F. Reis, A. P. Aguilar, and P. Tabuada, 'Control barrier functionbased quadratic programs introduce undesirable asymptotically stable equilibria,' IEEE Control Systems Letters , vol. 5, no. 2, pp. 731-736, 2021.
- [9] X. Tan and D. V. Dimarogonas, 'On the undesired equilibria induced by control barrier function based quadratic programs,' arXiv preprint arXiv:2104.14895 , 2021.
- [10] P. Mestres and J. Cort´ es, 'Optimization-based safe stabilizing feedback with guaranteed region of attraction,' IEEE Control Systems Letters , vol. 7, pp. 367-372, 2023.
- [11] Y. Yi, S. Koga, B. Gavrea, and N. Atanasov, 'Control synthesis for stability and safety by differential complementarity problem,' IEEE Control Systems Letters , vol. 7, pp. 895-900, 2023.
- [12] G. Notomista and M. Saveriano, 'Safety of dynamical systems with multiple non-convex unsafe sets using control barrier functions,' IEEE Control Systems Letters , vol. 6, pp. 1136-1141, 2022.
- [13] Y. Chen, P. Mestres, E. Dall'Anese, and J. Cort´ es, 'Characterization of the dynamical properties of safety filters for linear planar systems (extended manuscript),' https://drive.google.com/file/d/1YDd CdV8BVAr kAwEzOj01IeqpjhW7wS/view, 2024.
- [14] E. D. Sontag, Mathematical control theory: deterministic finite dimensional systems . Springer Science &amp; Business Media, 2013, vol. 6.
- [15] L. Perko, Differential Equations and Dynamical Systems , 3rd ed., ser. Texts in Applied Mathematics. New York: Springer, 2000, vol. 7.
- [16] M. Alyaseen, N. Atanasov, and J. Cort´ es, 'Continuity and boundedness of minimum-norm CBF-safe controllers,' arXiv preprint arXiv:2306.07398 , 2023.
- [17] G. Teschl, Ordinary Differential Equations And Dynamical Systems , ser. Graduate Studies in Mathematics. Providence, RI: American Mathematical Society, 2012, vol. 140.
- [18] P. Braun and C. M. Kellett, 'Comment on 'stabilization with guaranteed safety using control Lyapunov-barrier function',' Automatica , vol. 122, p. 109225, 2020.
Yiting Chen Pol Mestres Emiliano Dall'Anese Jorge Cortes
Characterization of the Dynamical Properties of Safety Filters for Linear Planar Systems
This paper studies the dynamical properties of closed-loop systems obtained
from control barrier function-based safety filters. We provide a sufficient and
necessary condition for the existence of undesirable equilibria and show that
the Jacobian matrix of the closed-loop system evaluated at an undesirable
equilibrium always has a nonpositive eigenvalue. In the special case of linear
planar systems and ellipsoidal obstacles, we give a complete characterization
of the dynamical properties of the corresponding closed-loop system. We show
that for underactuated systems, the safety filter always introduces a single
undesirable equilibrium, which is a saddle-point. We prove that all
trajectories outside the global stable manifold of such equilibrium converge to
the origin. In the fully actuated case, we discuss how the choice of nominal
controller affects the stability properties of the closed-loop system. Various
simulations illustrate our results.