## Safety-Critical Control with Offline-Online Neural Network Inference

Junhui Zhang, Member, IEEE , Sze Zheng Yong, Member, IEEE and Dimitra Panagou, Senior Member, IEEE

Abstract -This paper presents a safety-critical control framework for an ego agent moving among other agents. The approach infers the dynamics of the other agents, and incorporates the inferred quantities into the design of control barrier function (CBF)-based controllers for the ego agent. The inference method combines offline and online learning with radial basis function neural networks (RBFNNs). The RBFNNs are initially trained offline using collected datasets. To enhance the generalization of the RBFNNs, the weights are then updated online with new observations, without requiring persistent excitation conditions in order to enhance the applicability of the method. Additionally, we employ adaptive conformal prediction to quantify the estimation error of the RBFNNs for the other agents' dynamics, generating prediction sets to cover the true value with high probability. Finally, we formulate a CBF-based controller for the ego agent to guarantee safety with the desired confidence level by accounting for the prediction sets of other agents' dynamics in the sampleddata CBF conditions. Simulation results are provided to demonstrate the effectiveness of the proposed method.

function can be approximated by a RBFNN arbitrarily closely when the number of neurons is sufficiently large. To train the RBFNN, one can use either offline or online approaches [10], [11]. For offline learning methods, generalizing to unseen data and behaviors during the online implementation process is a challenge. On the other hand, training RBFNNs online poses challenges in achieving convergence of parameters, especially without the persistent excitation (PE) condition [12].

Index Terms -Control barrier function, inference, radial basis function network, offline-online learning, conformal prediction.

## I. INTRODUCTION

Ensuring safe motion for autonomous agents is an ongoing and challenging task. Control Barrier Functions (CBFs) [1]-[3] enforce the satisfaction of safety constraints along the system trajectories based on conditions that rely on the system model. Earlier work assumes that the dynamics of other agents are known within certain bounds [4], [5], which may be limiting or conservative for real-world applications. In this work, we develop an algorithm to infer the other agents' dynamics and integrate the learned dynamics into sampled data CBF designs.

Inferring the dynamics of other agents can be seen as a system identification problem. Unknown dynamics can be learned via non-parametric, data-driven approaches [6], or via parametric methods [7]. When there is no prior knowledge about the structure of the unknown dynamics, a three-layer neural network known as radial basis function neural network (RBFNN) can be employed [8] with the advantage that the universal approximation theorem [9] ensures that any continuous

This work was partially sponsored by the Office of Naval Research (ONR), under grant number N00014-20-1-2395 and NSF grant CNS-2312007. The views and conclusions contained herein are those of the authors only and should not be interpreted as representing those of ONR, the U.S. Navy, NSF or the U.S. Government.

Junhui Zhang is with the Department of Robotics at University of Michigan, Ann Arbor, MI 48109, USA (e-mail: junhuiz@umich.edu)

Sze Zheng Yong is with the Department of Mechanical and Industrial Engineering, Northeastern University, Boston, MA 02115, USA (e-mail: s.yong@northeastern.edu).

Dimitra Panagou is with the Department of Robotics and the Department of Aerospace Engineering at University of Michigan, Ann Arbor, MI 48109, USA (e-mail: dpanagou@umich.edu).

Our first aim in this work is to develop an offline-online inference learning method to infer other agents' dynamics. We use collected data to train the RBFNNs offline, and then the weights are updated during the online process. The offline-trained RBFNNs offer a 'warm start' for the onlinelearning process, ensuring that the estimation error remains relatively small before sufficient data are collected for training the RBFNNs online. On the other hand, the online-learning process enhances the generalization of offline trained RBFNNs since it keeps tuning their weights using new observations.

Our next aim is to quantify the estimation error of the learned dynamics so that we can use them in CBF-based control. Conformal prediction is a promising method to quantify uncertainty for arbitrary prediction algorithms [13]. In [14], adaptive conformal prediction (ACP) is proposed that does not rely on exchangeability of data points, although the confidence is only guaranteed in an average sense. This motivates us to employ ACP in the estimation error quantification of our offline-online neural network inference. Specifically, prediction sets are produced online by ACP such that the true value is covered within the prediction sets with a desired probability. Note that in [15], prediction sets are embedded into a model predictive control framework, while in the current work, the prediction sets are utilized within a sampled-data CBF design.

The main contributions of this work are summarized here: (1) We propose an offline-online method that utilizes RBFNNs to infer the unknown dynamics of other agents. (2) We leverage adaptive conformal prediction to online quantify the estimation error of our offline-online trained RBFNNs. (3) We incorporate prediction sets produced by adaptive conformal prediction into a sampled-data CBF framework, ensuring that the ego agent maintains safety with a desired confidence level.

Notations: The set of real numbers and non-negative real numbers are denoted by R and R + , respectively. A continuous function β : ( -b, a ) → -∞ ( , + ∞ ) is an extended class K function for a, b ∈ R + , denoted by β ∈ K , if it is strictly increasing and β (0) = 0 . ∥ · ∥ and ∥ · ∥ ∞ represent 2and infinity-norm, respectively. For a ∈ R , ⌈ a ⌉ is the smallest integer greater than or equal to a . Prob [ A ] is the probability of an event A L h x . f ( ) is the Lie derivative of h along f at x .

## II. PRELIMINARIES AND PROBLEM FORMULATION

1) Preliminaries: Consider a control affine system,

$$\dot { x } = f ( x ) + g ( x ) u,$$

where x ∈ X ⊂ R n and u ∈ U ⊂ R m are the state and control vectors of (1). f : X → R n and g : X → R n × m are Lipschitz-continuous functions.

Definition 1. [1] Let a safety set S , i.e., the set of safe states, for the system (1) be defined as S = { x ∈ X : h x ( ) ≥ 0 } , where h : X → R is a continuously differentiable function. The function h is a control barrier function for system (1) on the set S if there exists β ∈ K , such that

$$\sup _ { u \in \mathbb { R } ^ { m } } \{ L _$$

2) Problem Formulation: In this work, we consider an ego agent with dynamics described by

$$\dot { x } _ { e } = f ( x _ { e } ) + g ( x$$

where x e ∈ X e ⊂ R n e and u e ∈ U e ⊂ R m e are the state and control input vectors of ego agent, respectively. f : X → e R n e and g : X e → R n e × m e are Lipschitz-continuous functions. The set of other agents is denoted by V = { 1 , · · · , N} . Let x i ∈ X i ⊂ R n i denote the state of the i th agent in V , and x = [ x T e , x T 1 , x T 2 , · · · , x T N ] ∈ R n , n = n e + ∑ N i =1 n i , which includes all the agents' states. Suppose that each agent i ∈ V observes x and designs its action based on the observations. Then, the dynamics of agent i are ˙ x i = F i ( x ) , where F i is unknown to ego agent.

The ego agent is assigned to track its reference trajectory x r e with reference control u r e , and avoid collisions with other agents. One can define functions h i : R n e × R n i → R to encode the distance of ego agent w.r.t agent i ∈ V , so that the set S i = { x e ∈ X e , x i ∈ X i : h i ( x , x e i ) ≥ 0 } encodes the safety set for the ego agent and agent i . The CBF condition of the ego agent w.r.t. an agent i ∈ V is

$$\begin{array} {$$

where β i ∈ K .

The objectives of this work are: Problem 1 : How can the ego agent infer other agents' dynamics? Problem 2 : How can the ego agent quantify the estimation error of the inference algorithm for other agents' dynamics? Problem 3 : How to design a safe controller for the ego agent incorporating inference of other agents, ensuring that the ego agent safely tracks its reference trajectory as closely as possible?

## III. METHODOLOGY

In our setup, the ego agent observes states of all agents and updates its control signal at time instants t k = k ∆ t, k ∈ { 0 1 2 , , , . . . } , where ∆ t &gt; 0 is a fixed and sufficiently small sampling time interval. A zero-order hold (ZOH) control law is utilized so that u e ( ) t = u e ( t k ) , for t ∈ [ t k , t k +1 ) . More specifically, at each sampling time instant t k , the ego agent observes states of all agents x ( t k ) , and needs to update its control signal u e ( t k ) based on a CBF-based controller. In

Fig. 1: The proposed framework of sampled data CBF-based controller with offline-online neural network inference: RBFNNs are utilized to learn other agents' dynamics and are trained in an offlineonline manner. Adaptive conformal prediction (ACP) quantifies the estimation error online, and the prediction sets based on ACP are incorporated into the sampled-data CBF-based controller.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00918v1-figure-33d5b287-fc5f-4995-b94a-b0724d47270f.jpg)

Section III-A, we study how to estimate ˙ x i ( t k ) at sampling time t k using online observations and offline collected data so that this estimate can be used in the CBF-based controller. Moreover, we characterize in Section III-B the estimation error with high average confidence using adaptive conformal prediction. Finally, since the satisfaction of (3) at sampling instants only is not sufficient for safety, we design CBF conditions for sampled-data systems [16] with the incorporation of estimates of other agents' dynamics and their estimation errors in Section III-C. This framework is summarized in Fig. 1 and to facilitate this, we make the following assumptions.

Assumption 1. For the ego agent, there exists δ e,k ≥ 0 such that ∥ x e ( ) t -x e ( t k ) ∥ ≤ δ e,k for ∀ t ∈ [ t k , t k +1 ) , i.e., the change of ego agent's state during the sampling time interval is bounded. 1

Assumption 2. For agent i ∈ V , there exist δ i,k ≥ 0 and η i,k ≥ 0 such that ∥ x i ( ) t -x i ( t k ) ∥ ≤ δ i,k for ∀ t ∈ [ t k , t k +1 ) . In addition, we assume that F i is Lipschitz-continuous, i.e., that ∥ ˙ x i ( ) t -˙ x i ( t k ) ∥ = ∥ F i ( x ( )) t -F i ( x ( t k )) ∥ ≤ c F i ∥ x ( ) t -x ( t k ) ∥ ≤ η i,k = c F i ( ∑ N j δ j,k + δ e,k ) , for t ∈ [ t k , t k +1 ) , where c F i is a known Lipschitz constant for F i .

Assumption 3. The historical dynamics of other agents (i.e., ˙ x i with a one-step time delay) can be computed from measured state trajectory, e.g., using finite difference via ˙ x i ( t k -1 ) = x i ( t k ) -x i ( t k -1 ) ∆ t . For simplicity, we assume the approximation error is zero. If desired, the approximation error can also be quantified and incorporated in a straightforward manner.

## A. Offline-Online Inference

First, we develop an offline-online framework to infer other agents' dynamics for addressing Problem 1 .

1) Offline Learning : To collect a dataset offline, the ego agent moves in a multi-agent environment with any safe policy.

1 The bound δ e,k can be quantified as in [16], that is, ∥ x e ( ) t -x e ( t k ) ∥ ≤ σ t ( -t k ) ≤ σ ∆ t , where σ = sup x e ∈X e ,u e ∈U e ∥ f ( x e ) + g x ( e ) u e ∥ .

The other agents act following their intentions in reaction to the ego agent. The states of all agents are observed and collected. This means we have x ( t k ) , where k ∈ { 0 1 , , · · · , N } . The corresponding dynamics of other agents are also obtained by Assumption 3, that is, ˙ x i ( t k ) , for all k ∈ { 0 1 , , · · · , N - } 1 .

Using the offline dataset, we first estimate the functional relationship F i ( x ) between the dynamics of other agents ˙ x i and the states x . From the well-known Universal Approximation Theorem [9], when the number of neurons is sufficiently large, the radial basis function neural network (RBFNN) can approximate any continuous function arbitrarily closely; thus, we use a RBFNN to learn the functional relationship F i ( x ) as

$$\hat { \dot { x } } _ { i } = \hat { W } _ { i$$

where ˆ ˙ x i is the estimate for ˙ x i , ˆ W i ∈ R ( M i +1) × n i is the weight. Φ ( i x ) = [ e -| x -c i, 1 | 2 2 ρ i, 2 1 , · · · , e -| x -c i,M | 2 2 ρ i,M i 2 , 1 ] T is the Gaussian basis function with centers C i = [ c i, 1 , . . . , c i,M i ] ∈ R n × M i and widths ρ = [ ρ , . . . , ρ ] T ∈ R M i .

i i, 1 i,M i

The goal of training the RBFNN is to obtain the optimal weights such that the average of the estimation errors ϵ i ( t k ) = ˆ ˙ x i ( t k ) -˙ x i ( t k ) across all samples is minimized, i.e.,

$$\min _ { W _ { i } } \frac { 1 } { N } \sum _ { k = 0 } ^ { N - 1 } \| \epsilon _ { i } ( t _ { k } ) \| ^ { 2 }. \quad \begin{array} { c c c c c } & \underset { ( 5 ) } { \text{pred} } \\ & \text{det} \\ & \text{desir} \end{array}.$$

The RBFNN can be trained offline with methods such as gradient-based methods [17] and the orthogonal leastsquares methods [18]. The centers in the basis functions can be optimized with weights together during the training process or determined by implementing clustering approaches on collected datasets, e.g., the K -means method. 2 The widths of radial basis functions are manually tuned, e.g., close to 1. 2) Online Learning : The trained RBFNNs may not accurately estimate behaviors of the other agents that do not belong to the collected dataset. To enhance the generalization of RBFNNs, the weights obtained from offline training are utilized as a warm initialization ˆ W t i ( 0 ) and updated with new observations during the online process. For the online learning algorithm, the convergence of weights usually relies on the conditions of Persistent Excitation (PE), which are challenging to check online. In [19], the concurrent learning method was proposed, guaranteeing convergence without PE conditions, provided that the recorded data (either offline or online) contains as many linearly independent elements as the dimensions of the basis function. The key idea is that the weights are updated using both instantaneous observed data and the recorded data D z,i = { x ( m ) | m = 1 , · · · , ξ i } that has as many linearly independent elements as the dimension of the basis function Φ i , when it becomes available, as follows:

$$& \text{$u$ as $function $v_{i}$, when $n$ occurs avaance, as $n$rows.} \quad \text{time-ser} \\ & \hat { W } _ { i } ( t _ { k } ) = & \text{$a$ recur} \\ & \begin{cases} \hat { W } _ { i } ( t _ { k - 1 } ) - \Gamma _ { 1 } \Phi _ { i } ( \mathbf x ( t _ { k - 1 } ) ) \epsilon _ { i } ( t _ { k - 1 } ) \Delta t, \text{ rank} ( Z _ { i } ) < M _ { i } + 1, & \text{consiste} \\ \hat { W } _ { i } ( t _ { k - 1 } ) - [ \Gamma _ { 1 } \Phi _ { i } ( \mathbf x ( t _ { k - 1 } ) ) \epsilon _ { i } ( t _ { k - 1 } ) \\ + \sum _ { m = 1 } ^ { \xi _ { i } } \Gamma _ { 2 } \Phi _ { i } ( \mathbf x ( m ) ) \epsilon _ { i } ( m ) ] \Delta t, & \text{otherwise}, & \text{from $[1$} \\ ( 6 ) \end{cases}$$

(6)

2 Note that offline training occurs before the system deployment, so it is acceptable for it to take a longer time, if required.

where Z i = [Φ ( i x (1)) , · · · , Φ ( i x ( ξ i ))] , the estimation error ϵ i ( t k -1 ) is defined in above (5), ϵ i ( m ) is the estimation error for the dynamics of agent i at the m th data point of the recorded data D z,i , and Γ 1 and Γ 2 are positive constants.

In the above, we adopt a standard adaptive law from the adaptive control literature before we collect enough data for D z,i , that is, when rank ( Z i ) &lt; M i + 1 , Z i = [Φ ( i x (1)) , · · · , Φ ( i x ( ξ i ))] and switch to concurrent learning when rank ( Z i ) = M i +1 . As discussed in [20, Theorem 1], when the recorded data satisfies the above rank condition, the weights converge to a small bounded set around its ideal value.

## B. Adaptive Conformal Prediction

Next, we aim to quantify the estimation error of the proposed offline-online inference method and solve Problem 2 .

At sampling time t k , the ego agent observes the state x ( t k ) , and the other agents' dynamics at t k -1 are calculated by Assumption 3. The weight estimation ˆ W i of RBFNN is updated online under (6) with data point ( x ( t k -1 ) , x ˙ ( t k -1 )) at t k -1 and recorded data D z . Subsequently, the estimate ˆ ˙ x i ( t k ) of other agents' dynamics at time t k are obtained using the RBFNN (4); however, no estimation error in the prediction of the RBFNN is available. Conformal prediction determines a prediction set that covers the true value with a desired probability. The original conformal prediction assumes exchangeability of data points. However, in our problem, the state x ( t k ) depends on the state x ( t k -1 ) and agents' actions at t k -1 , violating the exchangeability condition. In [14], an adaptive conformal prediction (ACP) method was further proposed that does not rely on the exchangability condition, with coverage confidence guarantees in an average sense. We employ the ACP to quantify the estimation error of dynamics at t k estimated by the RBFNN.

Let D c,i ( t k ) = { ( x ( t k -( L +1 -j ) ) , x ˙ i ( t k -( L +1 -j ) )) , j ∈ { 1 2 , , . . . , L }} denote the set consisting of the most recent L data points with observed states and calculated historical dynamics by Assumption 3 at the sampling time t k . At the subsequent sampling time, the oldest sample is removed, and a new sample is added to the set. In the online conformal prediction framework, D c,i ( t k ) acts as the calibration dataset. That is, the most recent L data points are utilized to predict the estimation error for time t k .

After updating the weight ˆ W i using (6) at time t k , the nonconformity scores, i.e., estimation errors of the RBFNN on the calibration dataset are calculated as:

$$\cdot \, \text{s.s.s.s.} \} \text{ that } \quad R _ { i, j } = \| \dot { x } _ { i } ( t _ { k - ( L + 1 - j ) } ) - \hat { W } _ { i } ^ { T } \Phi _ { i } ( \mathbf x ( t _ { k - ( L + 1 - j ) } ) ) \| _ { \infty }, \ ( 7 )$$

where j ∈ { 1 2 , , · · · , L } . Due to the dependence between time-series data points, inspired by ACP in [14], we use a recursively updated failure probability α i ( t k ) instead of consistently using a fixed target failure probability α i ∈ (0 1) , . From [15], the width of the prediction set for ˆ ˙ x i ( t k ) is the (1 -α i ( t k )) th quantile of the sequence of R i, 1 , · · · , R i,L , ∞ , and when α i ( t k ) ∈ [ 1 L +1 , 1) , the width can be calculated by:

$$\text{so it is} \quad q _ { i } ( t _ { k } ) = \lceil ( 1 - \alpha _ { i } ( t _ { k } ) ) ( L + 1 ) \rceil \text{th smallest of } \{ R _ { i, m } \} _ { m = 1 } ^ { L },$$

with some special cases discussed in Remark 1. Then, the prediction set C i ( t k ) is:

$$C _ { i } ( t _ { k } ) = \{ \dot { x } _ { i } \ | \ \| \dot { x } _ { i } - \hat { \dot { x } } _ { i } ( t _ { k } ) \| _ { \infty } \leq q _ { i } ( t _ { k } ) \}. \quad ( 9 ) \ \text{Theor} \.$$

For the adaptive law of miscoverage level α i ( t k ) , we define err i ( t k -1 ) = 1 if ˙ x i ( t k -1 ) / ∈ C i ( t k -1 ) , and err i ( t k -1 ) = 0 otherwise, where ˙ x i ( t k -1 ) is calculated by Assumption 3. Then, the α i ( t k ) increases if the prediction set successfully covers the true value at the previous time step and decreases otherwise. The update law of α i ( t k ) is given by

$$\alpha _ { i } ( t _ { k } ) = \alpha _ { i } ( t _ { k - 1 } ) + \gamma ( \alpha _ { i } - e r r _ { i } ( t _ { k - 1 } ) ), \quad ( 1 0 )$$

where γ is the learning rate. Then, from [14], [15], the prediction set in (9) is guaranteed to cover the true dynamics of other agents with high average probability:

Proposition 1. [14], [15] Consider a time horizon [ t 0 , T ] , with sampling time instants { t 0 , t 1 , t 2 , · · · , t K } . For each sampling time t k , where k ∈ { 1 2 , , · · · , K } , the adaptive failure probability α i ( t k ) is updated using recursion (10) with a learning rate γ &gt; 0 , a fixed target failure probability α i ∈ (0 1) , and an initial value α i ( t 0 ) , then, it holds that

$$\frac { 1 } { K } \sum _ { k = 0 } ^ { K - 1 } P r o b [ \dot { x } _ { i } ( t _ { k } ) \in C _ { i } ( t _ { k } ) ] \geq 1 - \alpha _ { i } - \mathcal { P } _ { i } ( K ), \ \ ( 1 1 ) \quad \text{with} \ \mathcal { D } _ { i }$$

$$\text{where } \mathcal { P } _ { i } ( K ) = \frac { \alpha _ { i } ( t _ { 0 } ) + \gamma } { K \gamma }, \, a n d \lim _ { K \to + \infty } \mathcal { P } _ { i } ( K ) = 0.$$

Remark 1. From the recursion (10) , it is possible for α i ( t k ) to fall outside the range [ 1 L +1 , 1) during the online process, causing the (8) to fail in finding the widths of prediction sets in these instants. To address this, we enforce q i ( t k ) = R i, max when α i ( t k ) &lt; 1 L +1 , ensuring a sufficiently large prediction set to definitely cover the true value. The value of R i, max can be determined by the potential maximum value of the dynamics of agent i in real-world applications. For example, if the dynamics represent the speed of an agent, R i, max can be set to the maximum speed limit of agent i . Similarly, we enforce q i ( t k ) = R i, min = 0 when α i ( t k ) ≥ 1 .

Remark 2. Note that the result in Proposition 1 is independent of the calibration dataset used, including the choice of a moving horizon of length L in Section III-B. This can be observed from scrutinizing its proof in [14, Lemma 4.1] that holds as long as err i ( t k -1 ) is binary, regardless of the choice of calibration set that determines it. However, the choice of calibration dataset does affect the width of the prediction sets and how to optimally select this with a constraint on the calibration dataset size is a subject of ongoing research.

## C. CBF-based Controller with Inference

To solve Problem 3 , given the prediction set at t k by ACP, we can further find a prediction set for all t ∈ [ t k , t k +1 ) . By Assumption 2, the change of agent i 's dynamics during the sampling time interval is bounded, that is, ∥ ˙ x i ( ) t -˙ x i ( t k ) ∥ ≤ η i,k , for t ∈ [ t k , t k +1 ) . Then, from (9), we can define

$$\mathcal { C } _ { i, k } = \{ \dot { x } _ { i } \ | \ \| \dot { x } _ { i } - \hat { \dot { x } } _ { i } ( t _ { k } ) \| _ { \infty } \leq q _ { i } ( t _ { k } ) + \eta _ { i, k } \}, \quad ( 1 2 ) \ \text{ apply } z$$

to ensure that ˙ x i ( ) t ∈ C i,k , ∀ t ∈ [ t k , t k +1 ) , if ˙ x i ( t k ) ∈ C i ( t k ) .

Using this, we derive sampled-data CBF conditions to guarantee safety with high average probability.

Theorem 1. Consider a time horizon [ t 0 , T ] with sampling time instants { t 0 , t 1 , t 2 , · · · , t K } and safety set S i for i ∈ V , let c f,i , c g,i and c β,i be Lipschitz constants for L h f i , L h g i and β i ( h i ) with respect to the 2-norm, respectively. Under Assumptions 1 to 3, if the ZOH control of ego agent at each t k ∈ { t 0 , t 1 , t 2 , · · · , t K } satisfies the following sampled-data CBF conditions:

$$\begin{smallmatrix} L _ { f } h _ { i } ( x _ { e } ( t _ { k } ), x _ { i } ( t _ { k } ) ) + L _ { g } h _ { i } ( x _ { e } ( t _ { k } ), x _ { i } ( t _ { k } ) ) u _ { e } ( t _ { k } ) \\ + M _ { i, k } \geq - \beta _ { i } ( h _ { i } ( x _ { e } ( t _ { k } ), x _ { i } ( t _ { k } ) ) ) + \phi _ { i, k }, \end{smallmatrix}$$

with ϕ i,k = ( c f,i + c g,i u e, max + c β,i )( δ e,k + δ i,k ) , u e, max = min u e ∈U e , ∥ u e ∥≤ δ δ and

$$M _ { i, k } = \min _ { \dot { x } _ { i } \in \mathcal { C } _ { i, k } } \{ \frac { \partial h _ { i } } { \partial x _ { i } } \dot { x } _ { i } \} = \sum _ { r = 1 } ^ { n _ { i } } b _ { i, k, r }, \quad ( 1 4 )$$

with b i,k,r = min { µ i,r ( ˆ ˙ x i,r ( t k ) -q i ( t k ) -η i,k ) , µ i,r ( ˆ ˙ x i,r ( t k )+ q i ( t k )+ η i,k ) } , where µ i,r is the r th column of ∂h i ∂x i and ˆ ˙ x i,r ( t k ) is the r th row of ˆ ˙ x i ( t k ) , then safety is guaranteed with a probability of at least 1 -α i -P i ( K ) on average, i.e.,

$$\ u u y \quad \frac { 1 } { K } \sum _ { k = 0 } ^ { K - 1 } P r o b [ ( x _ { e } ( t ), x _ { i } ( t ) ) \in S _ { i }, \forall t \in [ t _ { k }, t _ { k + 1 } ) ] \\ \geq 1 - \alpha _ { i } - \mathcal { P } _ { i } ( K ),$$

$$\text{with} \, \mathcal { P } _ { i } ( K ) = \frac { \alpha _ { i } ( t _ { 0 } ) + \gamma } { K \gamma } \ a n d \lim _ { K \to + \infty } \mathcal { P } _ { i } ( K ) = 0.$$

Proof: If ˙ x i ( t k ) ∈ C i ( t k ) , k &lt; K , then for all t ∈ [ t k , t k +1 ) , it holds that,

$$& \cdots \succ \nu, \, \cdots \, \widetilde { \nu } _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \nu } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rp } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rq } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rho } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rp } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rtype } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpr } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rPC } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpic } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \r pc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rphi } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rpc } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rclap } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rclip } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rphi } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \ralg } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rarity } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \racity } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \rangle } \nu _ { \r$$

The argument t in x e ( ) t , x i ( ) t , and ˙ x i ( ) t is omitted for brevity. Thus, under the CBF-ZOH control satisfying (13), if ˙ x i ( t k ) ∈ C i ( t k ) , we obtain ˙ h i ( x , x e i ) ≥ -β i ( h i ( x , x e i )) for all t ∈ [ t k , t k +1 ) and all k . Finally, combining with Proposition 1, the conclusion holds.

Based on Theorem 1, the sampled-data CBF-based optimization problem is constructed for each sampling time t k (this argument is omitted below for brevity) as:

$$\max \max \min \colon & \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \colon & \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$ } \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \  \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \text{$\infty$} \, \end{cases}$$

where i ∈ V , u r e is the reference control for ego agent. To summarize, the control input u e is generated by solving the CBF-based optimization (16) at each sampling time t k and we apply zero-order hold (ZOH) for t ∈ [ t k , t k +1 ) , which enables the ego agent to track its reference trajectory as closely as possible with guaranteed safety with high average confidence.

## IV. CASE STUDY

The proposed method is implemented in simulation to demonstrate its effectiveness. The scenario involves one ego agent and two other agents moving in a two-dimensional space. The objective of the ego agent is to track its reference trajectory while ensuring safety w.r.t. the other agents. Agent 1 is moving towards the ego agent. Agent 2 moves close to ego agent while going away from agent 1. The ego agent does not know the intents of the other agents, but it can measure the positions of other agents at each sampling time. Our inference algorithm is utilized to learn their dynamics. In the first step, offline datasets are collected by navigating the ego agent using a CBF-based controller that forces the agent to track various reference trajectories, such as circles, sine waves, and spirals, while remaining safe with respect to other agents. The two other agents move with single integrator model following their intentions in reaction to the ego agent.

The ego agent dynamics are given as ˙ p x = v cos( θ ) , ˙ p y = v sin( θ ) , ˙ θ = ω , where p x , p y are its position coordinates, and v and ω are its linear and angular velocities that serve as the control inputs. Let p x,i , p y,i , i = 1 2 , denote the positions of two other agents. We use RBFNNs with 8 neurons to capture the behaviors of agents 1 and 2, respectively. The widths in Gaussian basis function are set as 0.85. By employing newrb in MATLAB neural network toolbox, the RBFNNs are trained offline using the collected datasets, where the weights and centers are optimized to minimize the estimation errors.

We set the offline-obtained weights as initial values ˆ W t i ( 0 ) of the online weights, and fix the offline-obtained centers C i of RBFNNs. Then, the update law (6) is utilized to tune the weights ˆ W i online. We set the initial miscoverage/target failure probability as α i ( t 0 ) = α i = 0 01 . , and the learning rate as γ = 0 002 . . The prediction sets corresponding to the RBFNN at time t k are obtained by adaptive conformal prediction in Section III-B. For the ego agent, the CBF for the unicycle model with dynamic obstacles in [21] is chosen as safety sets and the control signal is obtained from the sampled-data CBFbased optimization problem in (16) at t k with ZOH. 3

Fig. 2: Estimation of the dynamics of agent 1 in the x direction.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00918v1-figure-f97a5233-3944-44f1-92ea-6cfe39d6f2e9.jpg)

The estimates of other agents' dynamics during online process are depicted in Figs. 2-5. We observe that the prediction sets cover the real values at all times, ensuring safety through the CBF-based control that incorporates prediction sets. We also compare the performance of the offline-online inference method via the estimation error against purely offline and online learning approaches, see Figs. 6-7.

Fig. 3: Estimation of the dynamics of agent 1 in the y direction.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00918v1-figure-18abcf7b-434b-4fd2-8be7-0679508a0d08.jpg)

Agent 2 in the X direction

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00918v1-figure-a0bf62b2-35e1-4c21-b68f-0b06e85bb99c.jpg)

Fig. 4: Estimation of the dynamics of agent 2 in the x direction.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00918v1-figure-f6d89f49-14b9-4dc8-a659-30a459f335b5.jpg)

Fig. 5: Estimation of the dynamics of agent 2 in the y direction.

Fig. 6: Comparison of estimation performance for the dynamics of agent 1 (proposed vs offline only vs online only).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00918v1-figure-22bc8f55-75bd-4d3c-b1e1-1964c4f0cd67.jpg)

Fig. 7: Comparison of estimation performance for dynamics of agent 2 (proposed vs offline only vs online only)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00918v1-figure-6faeffb5-9c66-4bc2-9d8b-588f92b7a382.jpg)

The widths of the prediction sets are shown in Fig. 8. Compared to earlier work in [5], which utilizes a conservative constant bound and does not attempt to infer other agents'

## Widths of prediction intervals

Fig. 8: Widths of prediction sets for agents 1 and 2 using adaptive conformal prediction.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00918v1-figure-1a72f846-c092-4061-97d9-41c5de32c361.jpg)

dynamics, our proposed method imposes less conservatism in the CBF-based conditions. Finally, the trajectories of the ego agent, agents 1 and 2 are depicted in Fig. 9, where the colors transitioning from dark to light represent the trajectories of the agents over time. The reference trajectory of the ego agent violates the safety constraint. Nevertheless, the ego agent minimizes the distance between its trajectory and the reference trajectory while ensuring safety, ultimately maintaining a distance of 1.1 from the agents 1 and 2. Collision is avoided during the entire process. In addition, we ran 20 simulation trials with different initial positions, and observed that collisions did not occur, which is well within the fixed target failure probability of α i = 0 01 . . This demonstrates the effectiveness of our method.

Fig. 9: Trajectories of ego agent, and agents 1 and 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00918v1-figure-5e30f1ec-61fb-431c-98c2-2c4df0bbf406.jpg)

Fig. 10: Distances between ego agent and agents 1 and 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00918v1-figure-68ab7c3e-4391-42bf-88e6-45da311e25a7.jpg)

## V. CONCLUSIONS

This paper presented an inference method for an ego agent to estimate the dynamics of other agents, and its integration within a safe control synthesis framework. Offline-trained RBFNNs are utilized to learn the dynamics of other agents, and their weights are updated online with new observations. Adaptive conformal prediction is utilized to quantify the estimation error of the RBFNNs for other agents' dynamics by generating prediction sets covering the true value with desired average confidence level. Subsequently, a sampled-data CBF-based optimization problem is formulated and solved to guarantee safety with the desired average confidence level by incorporating the prediction set. Finally, a case study demonstrated the effectiveness of the proposed method. Future work includes extending the methods to distributed scenarios involving cooperative/non-cooperative agents.

## VI. ACKNOWLEDGEMENTS

We would like to thank Hardik Parwana and Taekyung Kim for inspiring and insightful discussions on this work.

## REFERENCES

- [1] A. D. Ames, X. Xu, J. W. Grizzle, and P. Tabuada, 'Control barrier function based quadratic programs for safety critical systems,' IEEE Transactions on Automatic Control , vol. 62, no. 8, pp. 3861-3876, 2016.
- [2] K. Garg, J. Usevitch, J. Breeden, M. Black, D. Agrawal, H. Parwana, and D. Panagou, 'Advances in the theory of control barrier functions: Addressing practical challenges in safe control synthesis for autonomous and robotic systems,' Annual Reviews in Control , vol. 57, 2024.
- [3] J. Zhang, L. Ding, X. Lu, and W. Tang, 'A novel real-time control approach for sparse and safe frequency regulation in inverter intensive microgrids,' IEEE Transactions on Industry Applications , 2023.
- [4] A. Mustafa and D. Panagou, 'Adversary detection and resilient control for multiagent systems,' IEEE Transactions on Control of Network Systems , vol. 10, no. 1, pp. 355-367, 2022.
- [5] H. Parwana, A. Mustafa, and D. Panagou, 'Trust-based rate-tunable control barrier functions for non-cooperative multi-agent systems,' in IEEE 61st Conference on Decision and Control , 2022, pp. 2222-2229.
- [6] Z. Jin, M. Khajenejad, and S. Z. Yong, 'Robust data-driven control barrier functions for unknown continuous control affine systems,' IEEE Control Systems Letters , vol. 7, pp. 1309-1314, 2023.
- [7] B. T. Lopez, J.-J. E. Slotine, and J. P. How, 'Robust adaptive control barrier functions: An adaptive and data-driven approach to safety,' IEEE Control Systems Letters , vol. 5, no. 3, pp. 1031-1036, 2020.
- [8] M.-J. Lee and Y.-K. Choi, 'An adaptive neurocontroller using RBFN for robot manipulators,' IEEE Transactions on Industrial Electronics , vol. 51, no. 3, pp. 711-717, 2004.
- [9] J. Park and I. W. Sandberg, 'Universal approximation using radial-basisfunction networks,' Neural computation , vol. 3, no. 2, pp. 246-257, 1991.
- [10] M. Pazouki, Z. Wu, Z. Yang, and D. P. M¨ oller, 'An efficient learning method for rbf neural networks,' in 2015 International joint conference on neural networks (IJCNN) . IEEE, 2015, pp. 1-6.
- [11] A.-K. Seghouane and N. Shokouhi, 'Adaptive learning for robust radial basis function networks,' IEEE Transactions on Cybernetics , vol. 51, no. 5, pp. 2847-2856, 2019.
- [12] T. Zheng and C. Wang, 'Relationship between persistent excitation levels and rbf network structures, with application to performance analysis of deterministic learning,' IEEE Transactions on cybernetics , vol. 47, no. 10, pp. 3380-3392, 2017.
- [13] A. N. Angelopoulos and S. Bates, 'Conformal prediction: A gentle introduction,' Foundations and Trends® in Machine Learning , vol. 16, no. 4, pp. 494-591, 2023.
- [14] I. Gibbs and E. Candes, 'Adaptive conformal inference under distribution shift,' Advances in Neural Information Processing Systems , vol. 34, pp. 1660-1672, 2021.
- [15] A. Dixit, L. Lindemann, S. X. Wei, M. Cleaveland, G. J. Pappas, and J. W. Burdick, 'Adaptive conformal prediction for motion planning among dynamic agents,' in Learning for Dynamics and Control Conference . PMLR, 2023, pp. 300-314.
- [16] J. Breeden, K. Garg, and D. Panagou, 'Control barrier functions in sampled-data systems,' IEEE Control Systems Letters , vol. 6, pp. 367372, 2021.
- [17] H.-G. Han, M.-L. Ma, and J.-F. Qiao, 'Accelerated gradient algorithm for RBF neural network,' Neurocomputing , vol. 441, pp. 237-247, 2021.
- [18] S. Chen, C. Cowan, and P. Grant, 'Orthogonal least squares learning algorithm for radial basis function networks,' IEEE Trans. Neural Netw , vol. 2, no. 2, pp. 302-309, 1991.
- [19] G. Chowdhary and E. Johnson, 'Concurrent learning for convergence in adaptive control without persistency of excitation,' in 49th IEEE Conference on Decision and Control , 2010, pp. 3674-3679.
- [20] O. Djaneye-Boundjou and R. Ord´ o˜ nez, 'Gradient-based discrete-time concurrent learning for standalone function approximation,' IEEE Transactions on Automatic Control , vol. 65, no. 2, pp. 749-756, 2019.
- [21] G. Wu and K. Sreenath, 'Safety-critical control of a planar quadrotor,' in 2016 American control conference (ACC) . IEEE, 2016, pp. 2252-2258.
Junhui Zhang Sze Zheng Yong Dimitra Panagou
Safety-Critical Control with Offline-Online Neural Network Inference
This paper presents a safety-critical control framework for an ego agent
moving among other agents. The approach infers the dynamics of the other
agents, and incorporates the inferred quantities into the design of control
barrier function (CBF)-based controllers for the ego agent. The inference
method combines offline and online learning with radial basis function neural
networks (RBFNNs). The RBFNNs are initially trained offline using collected
datasets. To enhance the generalization of the RBFNNs, the weights are then
updated online with new observations, without requiring persistent excitation
conditions in order to enhance the applicability of the method. Additionally,
we employ adaptive conformal prediction to quantify the estimation error of the
RBFNNs for the other agents' dynamics, generating prediction sets to cover the
true value with high probability. Finally, we formulate a CBF-based controller
for the ego agent to guarantee safety with the desired confidence level by
accounting for the prediction sets of other agents' dynamics in the
sampled-data CBF conditions. Simulation results are provided to demonstrate the
effectiveness of the proposed method.