HiDeS: a higher-order-derivative-supervised neural ordinary differential equation for multi-robot systems and opinion dynamics

Abstract

Addressing the limitations of current neural ordinary differential equations (NODEs) in modeling and predicting complex dynamics, this paper introduces a higher-order-derivative-supervised (HiDeS) NODE. This approach extends traditional NODE frameworks by incorporating higherorder derivatives and their interactions into the modeling process, thereby enabling the capture of intricate system behaviors. Unlike conventional NODEs that utilize only the state vector as a supervised signal, HiDeS NODE employs both the state vector and its higher-order derivatives, significantly enhancing predictive accuracy. We demonstrate the superiority of HiDeS NODE through extensive experiments in multi-robot systems and opinion dynamics, areas where existing models struggle due to the complexity of interactions and dynamics. Our results indicate that HiDeS NODE offers improved modeling and prediction capabilities, paving the way for new applications in dynamic systems. This research not only proposes an expressive and predictive framework for dynamic systems but also marks the first application of NODEs to the fields of multi-robot systems and opinion dynamics, suggesting broad potential for future interdisciplinary work. The code is available at https://github.com/MengLi-Thea/HiDeS-A-Higher-Order-Derivative-Supervised-Neural-Ordinary-Differential-Equation.

1 Introduction

As a learnable model parameterized by θ ∈ ℝⁿ, a standard neural ordinary differential equation (NODE) x.=ϕθ(x,t) x.=ϕθ(x,t) is particularly adept at representing complex and nonlinear dynamics (Chen et al., 2018; Liufu et al., 2024), where x ∈ ℝ^d is the state at time t, x.=dx/dt x.=dx/dt denotes the time derivative of x, and ϕ(x, t) is a vector field with ϕ ∈ (ℝ^d×ℝ → ℝ^d) being a function of x and t. Its strength lies in processing time-variant data and adaptively learning from it. This modeling flexibility renders NODEs great potential for the intricate nature of dynamic systems (Hua et al., 2023; Wang et al., 2023; Jin et al., 2024), enabling a more nuanced understanding of complex dynamic systems.

Despite these strengths, the standard NODE encounters expressivity limitations, failing to model functions like NOT operations (Kidger, 2021; Xu et al., 2023). The NOT operation [i.e., (0, 1) → (1, 0)] involves trajectories that necessarily intersect, presenting a challenge for standard NODEs that cannot model intersecting trajectories due to their first-order nature. The NODE with momentum, which can be regarded as a second-order ODE, improves the expressive capability (Sander et al., 2021): x¨=c0x.+c1ϕθ(x,t) x¨=c0x.+c1ϕθ(x,t) , where c₀ and c₁ are constants. Nonetheless, it can only express limited dynamics due to the linear relationship of x. x. and ϕ_θ(x, t). Besides, it cannot model interactions between x. x. and x. The second-order NODE (SONODE) presented in Norcliffe et al. (2020) seeks to address this limitation by modeling the interactions between x. x. and x. However, SONODE cannot model interactions between higher-order derivatives and x, and the supervised signal used in training is only the ground-truth value of x, which confines its scope and limits its prediction capability.

To surmount these challenges, we propose a higher-order-derivative-supervised NODE (HiDeS NODE) that is able to model interactions between higher-order derivatives and x. This approach not only expands the expressive range of NODEs but also enhances predictive ability through employing the state vector and its higher-order derivatives as supervised signals, surpassing the modeling and predicting performance of existing NODEs.

This paper evaluates the effectiveness of the HiDeS NODE in the realms of multi-robot systems and opinion dynamics, key areas of dynamic systems, both domains that inherently involve complex interactions and communication (Granha et al., 2022). In multi-robot systems, conventional analytic solutions fall short in high-dimensional control tasks (Károly et al., 2021), such as multi-robot grasping and motion control. NODEs, in contrast, offer a promising avenue for modeling and controlling complex dynamic interactions in a continuous, efficient, and adaptable manner in multi-robot systems. Regarding opinion dynamics research, the primary objective is to decipher the underlying mechanisms and influences that catalyze shifts in opinions. Existing methodologies for learning opinion dynamics overlook the critical prior knowledge that opinion dynamics can be described as an ODE formulated as x.=ϕ(x,t) x.=ϕ(x,t) . ODEs are particularly well-suited for modeling the fluid nature of opinion dynamics due to their inherent capacity to capture the dynamics of evolving systems. However, contemporary models employed in learning opinion dynamics underutilize this foundational knowledge. This oversight hampers their ability to effectively capture the nuanced and intricate nature of opinion evolution. Furthermore, the complexities inherent in the evolution of opinions present considerable challenges to the application of existing NODEs in both modeling and forecasting the trajectories of opinion dynamics. The HiDeS NODE conquers these aspects, providing a more effective tool for understanding and predicting opinion evolution.

To bridge these gaps, we propose a new NODE, termed HiDeS NODEs, for modeling and predicting tasks in multi-robot control and opinion dynamics. Figure 1 illustrates the framework of the HiDeS NODE, and Table 1 qualitatively demonstrates the HiDeS NODE's superiority compared with existing NODEs.

Figure 1

Figure 1. Framework of this paper. (Top) Evolution from first to third-order NODEs, highlighting their progressively sophisticated ability to model complex trajectories where higher orders allow for more intricate behaviors. (Middle) Predicting precision with and without using higher-order derivatives as supervised signals, showing the latter's superior approximation of ground truth. (Bottom) Practical applications of HiDeS NODEs in multi-robot systems and opinion dynamics.

Table 1

Table 1. Comparisons among different NODEs.

The contributions of this paper are demonstrated as follows:

• We propose the HiDeS NODE, a novel approach for modeling the intricacies of dynamics. The HiDeS NODE excels in modeling and predicting interactions among higher-order derivatives within dynamic systems. This advancement provides a more accurate and nuanced representation of dynamic systems.

• The HiDeS NODE integrates higher-order derivatives as supervised signals, significantly enhancing the ability to predict dynamical behaviors.

• We examine the versatility and effectiveness of the proposed HiDeS NODE through its application in two distinct yet complex fields: Multi-robot systems and opinion dynamics. In these fields, the model's ability to capture and predict intricate system dynamics is evaluated.

• To our knowledge, this is the first time that the NODE is introduced for opinion dynamics and multi-robot-system control. Application of the proposed HiDeS NODE to these fields unveils new avenues for both the advancement of NODE methodologies and the nuanced modeling of opinion dynamics and multi-robot-system control.

2 Related work

In this section, we briefly review three lines of research that are close to our work: NODEs, multi-robot-system control methods, and opinion dynamics modeling.

2.1 NODEs

The intersection of neural networks and differential equations, especially interpreting residual networks (ResNets) as discretized ODEs, spurs the development of NODEs (Weinan, 2017; Cui et al., 2023; Ruiz-Balet and Zuazua, 2023). NODEs integrate black-box ODE solvers and neural networks to parameterize the hidden state's derivative. This integration substantially advances time-series modeling, offering robust function approximation and handling of irregular data (Chen et al., 2018; Kidger, 2021). However, standard NODEs encounter representational constraints without dimensionality augmentation, constraining their universal approximation capabilities for certain functions (Dupont et al., 2019).

Research pivots toward higher-dimensional NODEs to overcome these limitations. Momentum-enhanced ResNets, representing second-order NODE extensions, exhibit enhanced capability in modeling non-homeomorphic dynamics and demonstrate improved convergence properties (Sander et al., 2021). In parallel, augmented NODEs, by expanding the solution space, facilitate the learning of more complex functions through simpler dynamic flows, thereby sidestepping the limitations of the vector field's general-representation property (Kidger, 2021). Nonetheless, augmented NODEs introduce challenges in interpretability and alter the loss landscape's structure (Norcliffe et al., 2020). A specific iteration of augmented NODEs, termed second-order NODEs (SONODEs) (Norcliffe et al., 2020), captures more intricate behaviors by integrating second-order dynamics, effectively combining the principles of coupled augmented NODEs. Additionally, the advent of heavy ball NODEs (HBNODEs) (Xia et al., 2021) marks a significant advancement. HBNODEs incorporate the classical momentum accelerated gradient descent method and adeptly mitigate the vanishing gradient problem, thereby enhancing the model's capacity in learning long-term dependencies in sequential data (Xia et al., 2021).

2.2 Multi-robot-system control

Multi-robot systems provide significant benefits in tasks that demand the duplication of effort, risk reduction, or adaptability, offering distinct advantages over single-robot systems (Hichri et al., 2022; Kwa et al., 2022). Multi-robot-system control methods can be categorized into deterministic methods with fixed forms and learning-based methods (Pierpaoli et al., 2021). However, deterministic methods lack flexibility and adaptability in dynamic or unpredictable environments (Liu et al., 2023). In order to overcome these defects, learning methods are increasingly applied to multi-robot control problems. Adaptation methods, for instance, are proposed to enhance trajectory prediction efficiency in multi-agent systems (Aydemir et al., 2023). Furthermore, the parameter-adaptive learning methods are improved through iterative parametric learning controllers (Yu and Chen, 2023). Additionally, neural network-based adaptive learning methods are utilized to learn unknown fault functions, ensuring cooperative tracking in distributed multi-robot systems (Khalili et al., 2020). Despite these advancements, existing methods often fall short in naturally and efficiently modeling the dynamics that often can be described as an ODE, a gap that NODEs can potentially fill.

2.3 Opinion dynamics modeling

Opinion dynamics studies how opinions form and evolve over time through interactions with individuals and environments. Researchers propose various mathematical models to understand and predict the dynamics of opinions. These include continuous-time models such as the DeGroot model (Wu et al., 2023), Hegselmann-Krause model, and bounded confidence model (Kolarijani et al., 2021), as well as discrete-time models like the Ising model, Voter model, and Friedkin and Johnsen model (Baumann et al., 2020; Ao and Jia, 2023; Peng et al., 2023). However, these models, with their fixed forms, lack the flexibility to model the evolution of opinions independently.

In response to these limitations, researchers leverage advances in neural networks to utilize their nonlinear relation approximation ability for learning complex opinion dynamics. An early approach introduces a linear influence model that learns edge influence strength from real data (De et al., 2014). Unlike traditional models, this linear model represents a foundational step in opinion dynamics learning methods, but its simplicity fails to capture the complexity of societal opinion dynamics. Furthering this exploration, SLANT (De et al., 2016; Zhu et al., 2020) introduces a linear model of latent opinions driven by stochastic differential equations (SDEs) using historical, fine-grained event data. Subsequently, SLANT+ (Kulkarni et al., 2017) extends this model with a nonlinear generative model and a network-guided recurrent neural network (RNN) architecture. This model underscores the importance of nonlinearity in designing opinion dynamics models. However, the RNN architecture it relies on faces the challenge of the vanishing gradient problem, hindering long-term predictions of opinion flow. Learnable opinion dynamics model (LODM) (Monti et al., 2020) emerges as a learnable generalization of an opinion dynamics model, combining the causal interpretability of traditional agent-based models with data-driven approaches. Additionally,Okawa and Iwata (2022) introduces the sociologically-informed neural network (SINN), a novel hybrid approach that integrates sociological and social psychological theories with data-driven neural networks to model and predict opinion dynamics in social networks. Despite these advances, current models do not fully exploit the prior knowledge of differential equations in opinion evolution, nor do they effectively model higher-order derivatives.

3 Materials and methods

In this section, formal descriptions and analyses of the proposed HiDeS NODE are provided. Table 2 presents the main symbols and notations used throughout this paper to ensure clarity and ease of understanding.

Table 2

Table 2. Main symbols and notations.

3.1 Formulation of the HiDeS NODE

The HiDeS NODE has two unique features for modeling and predicting opinion evolution. The first is that the HiDeS NODE is a higher-order NODE that is able to model interactions of higher-order derivatives of the opinion variable x. The second is that the HiDeS NODE adopts higher-order derivatives as supervisory signals to predict opinion evolution better. The HiDeS NODE is described as Equation (1):

⎡⎣⎢⎢⎢⎢⎢⎢⎢x˙(t)x¨(t)x...(t)⋮x(c)(t)⎤⎦⎥⎥⎥⎥⎥⎥⎥=ϕθ⎛⎝⎜⎜⎜⎜⎜⎜⎜⎡⎣⎢⎢⎢⎢⎢⎢⎢x(t)x˙(t)x¨(t)⋮x(c−1)(t)⎤⎦⎥⎥⎥⎥⎥⎥⎥,t⎞⎠⎟⎟⎟⎟⎟⎟⎟, (1) [x˙(t)x¨(t)x...(t)⋮x(c)(t)]=ϕθ([x(t)x˙(t)x¨(t)⋮x(c−1)(t)],t), (1)

where x(t) ∈ ℝ^d is a time-varying vector representing the opinion of d individuals; t∈[t∨,tˆ] t∈[t∨,t^] is the time; Vectors x.(t),x¨(t),x...(t) x.(t),x¨(t),x...(t) , and x^(c)(t) correspond to the first, second, third, and c-th order time derivatives of x, respectively; The function ϕ:ℝ^cd×ℝ → ℝ^cd is parameterized by a neural network with the parameter θ ∈ ℝⁿ. Note that [x.(t),x¨(t),x...(t),…,x(c)(t)]⊤∈Rcd [x.(t),x¨(t),x...(t),…,x(c)(t)]⊤∈ℝcd is a concatenation of higher-order derivatives, where the superscript ^⊤ means a transpose of a vector, and we call a HiDeS NODE with up to c-th order time derivatives in this concatenation as the HiDeS-c NODE. To enhance readability and avoid redundancy, we may omit “(t)” in certain contexts where the time dependency is understood and does not affect the meaning or clarity of the mathematical expressions.

Remark 1. One advantage of a HiDeS NODE is that it is able to model nonlinear interactions between higher-order derivatives and x. In practice, multiple higher-order derivatives and x can interact with each other. For example, there can be terms like x.⊗x¨ x.⊗x¨ in the vector field, where ⊗ is the Hadamard product.

It can be seen that the standard NODE (Chen et al., 2018) is a HiDeS-1 NODE, and if we just focus on the formulation, SONODE (Norcliffe et al., 2020) can be regarded as a HiDeS-2 NODE. In fact, the HiDeS-2 NODE distinguishes itself from SONODE due to its unique training process.

3.2 Training of the HiDeS NODE

Existing variants of the NODE utilize the ground-truth value of x(t) as the label for training. Differently, the HiDeS NODE adopts the entire [x(t),x.(t),x¨(t),…,x(c−1)(t)]⊤ [x(t),x.(t),x¨(t),…,x(c-1)(t)]⊤ as the label (the model's prediction [x.(t),x¨(t),x...(t),…,x(c)(t)]⊤ [x.(t),x¨(t),x...(t),…,x(c)(t)]⊤ is integrated first to get [x(t),x.(t),x¨(t),…,x(c−1)(t)]⊤ [x(t),x.(t),x¨(t),…,x(c-1)(t)]⊤ ). This approach is beneficial for predicting the future evolution of x(t). Training the entire [x(t),x.(t),x¨(t),…,x(c−1)(t)]⊤ [x(t),x.(t),x¨(t),…,x(c-1)(t)]⊤ gives accurate approximations of all these variables. Since the prediction of the next time step for the HiDeS NODE relies on the entire [x(t),x.(t),x¨(t),…,x(c−1)(t)]⊤ [x(t),x.(t),x¨(t),…,x(c-1)(t)]⊤ , the training strategy of the HiDeS NODE leads to a better prediction performance compared to only training with the ground-truth value of x(t). When the model is predicting the next x(t), utilizing only the ground-truth value of x(t) as the label may lead to an inaccurate result because the basic information it relies on is inaccurate. The inclusion of the derivatives ensures that the model is sensitive to not just the position or condition at a given time but also to the trends and patterns of change, which are critical for forecasting. An explanation is illustrated in Figure 1.

3.3 Inexpressible trajectories of the HiDeS NODE

The superior expressive capability of the HiDeS NODE comes from two aspects.

The first is that lower-order NODEs have limitations in modeling trajectories that require the representation of higher-order dynamics. Consider a trajectory that requires an abrupt change in its acceleration (second derivative of x), which is not expressible in a first-order system but can be expressed in a second-order system [x.(t),x¨(t)]⊤=ϕθ([x(t),x.(t)]⊤,t) [x.(t),x¨(t)]⊤=ϕθ([x(t),x.(t)]⊤,t) . Similarly, trajectories requiring changes in the third derivative (jerk) are not expressible in a second-order system but can be captured in a third-order system, and so on. As a result, there exist trajectories that can not be expressed by [x.(t),x¨(t),x...(t),…,x(c)(t)]⊤=ϕθ([x(t),x.(t),x¨(t),…,x(c−1)(t)]⊤,t) [x.(t),x¨(t),x...(t),…,x(c)(t)]⊤=ϕθ([x(t),x.(t),x¨(t),…,x(c-1)(t)]⊤,t) but can be expressed by [x.(t),x¨(t),x...(t),…,x(c+1)(t)]⊤=ϕθ([x(t),x.(t),x¨(t),…,x(c)(t)]⊤,t) [x.(t),x¨(t),x...(t),…,x(c+1)(t)]⊤=ϕθ([x(t),x.(t),x¨(t),…,x(c)(t)]⊤,t) .

The second origin of the superior expressive capability of the HiDeS NODE is that it alleviates the restriction that trajectories cannot cross. One major limitation of the standard NODE is that trajectories under different initial conditions cannot intersect, which constrains its expressive capability. In the following, we show how this constraint is able to be eliminated by the HiDeS NODE.

Theorem 1 (Inexpressible trajectories of a HiDeS-c NODE). Assume that the function ϕθ(ω,t):Rcd×R→Rcd ϕθ(ω,t):ℝcd×ℝ→ℝcd with t∈[t∨,tˆ] t∈[t∨,t^] is Lipschitz continuous w.r.t. ω ∈ ℝ^cd. Consider a HiDeS-c NODE governed by Equation (2):

[x.(t),x¨(t),x...(t),…,x(c)(t)]⊤      =ϕθ([x(t),x.(t),x¨(t),…,x(c−1)(t)]⊤,t),  (2) [x.(t),x¨(t),x...(t),…,x(c)(t)]⊤     =ϕθ([x(t),x.(t),x¨(t),…,x(c-1)(t)]⊤,t), (2)

where x(t) ∈ ℝ^d is the state vector, ϕθ:Rcd×R→Rcd ϕθ:ℝcd×ℝ→ℝcd is a continuously differentiable function parameterized by neural network parameters θ. For any two initial conditions ω(t∨)=[x(t∨),x.(t∨),…,x(c−1)(t∨)]⊤ ω(t∨)=[x(t∨),x.(t∨),…,x(c-1)(t∨)]⊤ and ω˜(t∨)=[x˜(t∨),x˜.(t∨),…,x˜(c−1)(t∨)]⊤ ω~(t∨)=[x~(t∨),x~.(t∨),…,x~(c-1)(t∨)]⊤ , trajectories that require ω(t) and ω˜(t) ω~(t) to cross over the interval [t∨,tˆ] [t∨,t^] are inexpressible by a HiDeS-c NODE.

Proof. Define the extended state vector ω ∈ ℝ^cd as ω(t)=[x(t),x.(t),x¨(t),…,x(c−1)(t)]⊤ ω(t)=[x(t),x.(t),x¨(t),…,x(c-1)(t)]⊤ . The HiDeS-c NODE can be represented as Equation (3):

ω.=ϕθ(ω(t),t). (3) ω.=ϕθ(ω(t),t). (3)

Given that ϕ_θ is Lipschitz continuous w.r.t. ω, the Picard-Lindelöf theorem (Anil Kumar et al., 2022; Zhang et al., 2023) assures the existence of a unique solution ω(t) for a given initial condition ω(0). This uniqueness implies that for any two distinct initial conditions ( ω(t∨) ω(t∨) ) and ω˜(t∨) ω~(t∨) , the resulting trajectories ω(t) and ω˜(t) ω~(t) do not cross over [t∨,tˆ] [t∨,t^] . As a result, trajectories that require ω(t) and ω˜(t) ω~(t) to cross over the interval [t∨,tˆ] [t∨,t^] are inexpressible by a HiDeS-c NODE. The proof is thus completed. ⎕

Remark 2. From Theorem 1, it can be seen that as the order c increases, the degree of freedom for avoiding the crossing of x(t) and x˜(t) x~(t) increases. In practice, there are some trajectories cross in the x-t space or in the phase space [i.e., x- x. x. - x¨ x¨ -…-x^(c−1)-t space], so the HiDeS-c NODE provides a better capability to model these dynamics compared with the standard NODE (a HiDeS-1 NODE) (Chen et al., 2018) and other second-order NODEs (HiDeS-2 NODEs) (Norcliffe et al., 2020; Sander et al., 2021). An intuitive understanding is that additional dimensions provide new directions to make trajectories elude each other, which is illustrated in Figure 1. Consider two actual evolution trajectories under different initial conditions, x^*(t) and x˜∗(t) x~*(t) , which may intersect at some time points if there are no restrictions. However, for c = 1, the trajectories x(t) and x˜(t) x~(t) generated by the standard NODE cannot intersect due to its nature of first-order ODEs. This limitation means they cannot accurately approximate x^*(t) and x˜∗(t) x~*(t) in cases where the actual trajectories intersect. Theorem 1 from our manuscript implies similar limitations for higher-order NODEs, but with increasing order c, the trajectories have more freedom, reducing the limitations.

3.4 The HiDeS NODE's utilization of historical information

Due to the introduction of higher-order derivatives, the HiDeS NODE implicitly uses historical state information for predicting the next state. The reason is that higher-order derivatives can be approximated by historical states. For example, the first derivative x. x. at the k-th moment t_k can be approximated as x.(tk)≈(x(tk)−x(tk−1))/Δt x.(tk)≈(x(tk)-x(tk-1))/Δt . Similarly, the second derivative x¨ x¨ can be approximated as Equation (4):

x¨(tk)≈(x(tk)−x(tk−1)Δt)−(x(tk−1)−x(tk−2)Δt)Δt. (4) x¨(tk)≈(x(tk)-x(tk-1)Δt)-(x(tk-1)-x(tk-2)Δt)Δt. (4)

Iteratively, we have Equation (5):

x(c)(tk)≈∑i=0c(−1)i( ci)x(tk−i)/(Δtc), (5) x(c)(tk)≈∑i=0c(−1)i( ic)x(tk−i)/(Δtc), (5)

where the binomial coefficient ( ci)=c!/(i!(c−i)!) ( ic)=c!/(i!(c-i)!) represents the combinatorial number of ways to choose i elements from a set of c elements. Consequently, ϕθ([x(t),x.(t),x¨(t),…,x(c−1)(t)]⊤,t) ϕθ([x(t),x.(t),x¨(t),…,x(c-1)(t)]⊤,t) in the HiDeS NODE (Equation 1) can be approximated as a function of historical states as in Equation (6):

[x.(t),x..(t),x...(t),…,x(c)(t)]⊤      =ϕθ([x(t),x.(t),x..(t),…,x(c−1)(t)]⊤,t)   ≈ ϕv(x(tk), x(tk−1),…, x(tk−c),tk), (6) [x.(t),x..(t),x...(t),…,x(c)(t)]⊤=ϕθ([x(t),x.(t),x..(t),…,x(c−1)(t)]⊤,t) ≈ ϕv(x(tk), x(tk−1),…, x(tk−c),tk), (6)

where ϕϑ:Rd×Rd×⋯×R→Rcd ϕϑ:ℝd×ℝd×⋯×ℝ→ℝcd is a function parameterized by ϑ. The HiDeS NODE's utilization of historical information could enhance the prediction of the next state.

3.5 Implementation

We provide Algorithm 1 to show the process of constructing, training, and using a HiDeS NODE. Besides, structures of HiDeS-3 NODEs for multi-robot systems and for opinion dynamics are shown in Figures 2A, B, respectively. In Figure 2, the inputs to the system are a concatenation of the initial state x(0), the initial velocity x.(0) x.(0) , and the initial acceleration x¨(0) x¨(0) , which are fed into an ODE solver alongside the parametrized function ϕθ([x(t),x.(t),x¨(t)]⊤) ϕθ([x(t),x.(t),x¨(t)]⊤) to compute the state x(t), velocity x.(t) x.(t) , and acceleration x¨(t) x¨(t) at time t.

Algorithm 1

Algorithm 1. Algorithm of HiDeS NODE.

Figure 2

Figure 2. Structure of HiDeS-3 NODEs. (A) Structure of HiDeS-3 NODE for multi-robot systems. (B) Structure of HiDeS-3 NODE for opinion dynamics.

4 Results

In this section, we conduct experiments to evaluate the effectiveness of our models, HiDeS-2 NODE and HiDeS-3 NODE, by comparing them with baseline models [standard NODE (Chen et al., 2018) and SONODE (Norcliffe et al., 2020)] on two applications: multi-robot control and opinion dynamics. These baseline models have the same configurations in terms of network architecture, optimizer, epochs, and learning rate, ensuring a fair comparison. Notably, our models utilize higher-order derivatives as supervised signals, crucial for accurately capturing the intricate, nonlinear evolution of opinions over time. The implementation details of our experiments are as follows.

4.1 Experimental settings

4.1.1 Settings for multi-robot-system control

In multi-robot-system control, each NODE block is composed of three fully connected layers, each succeeded by a Tanh activation function, as shown in Figure 2A. The NODE block undergoes forward propagation 200 times, evolving from t = 0 to t = 20, in order to develop a deep model. A weighted loss ℓ that emphasizes the trajectory's later stages is applied: ℓW=∑kˆk=0(k/kˆ)pℓ(ω(tk),ω∗(tk)) ℓW=∑k=0k^(k/k^)pℓ(ω(tk),ω*(tk)) , where kˆ k^ is the total number of steps, p>0 is a scalar, ℓ(·, ·) is a loss function, and ω∗(tk) ω*(tk) is the ground truth of ω(t_k). In simulations, p is taken as 4. All models are trained for 1,000 epochs using the Adam optimizer and a cosine annealing scheduler with a base learning rate of 0.01.

4.1.2 Settings for opinion dynamics

In simulations of opinion dynamics, each block of NODEs consists of three fully connected layers, each followed by an exponential linear unit (ELU) activation function, as shown in Figure 2B. The block of NODEs loops in the forward propagation 50 times from t = 0 to t = 5 to form a deep model. To enhance the training stability, we incorporate residual connections. The optimizer employed is NAdam (Dozat, 2016; Li et al., 2022), with an initial learning rate of 0.01, modulated using a cosine annealing learning rate scheduler (Jin et al., 2022). We train the models over 1,000 and 2,000 epochs, respectively. The loss function is the mean square error between the predicted and actual x(t),x.(t),x¨(t),…,x(c−1)(t) x(t),x.(t),x¨(t),…,x(c-1)(t) . Given the higher dimensions of HiDeS NODEs and SONODEs compared to standard NODEs, we introduce an auxiliary loss to ensure a fair comparison. The auxiliary loss is evaluated based on predicted and actual x(t), and it is only used for comparisons rather than for training purposes.

4.2 Target chasing of multi-robot system

In this section, the NODE models are used to control a multi-robot system to chase a moving target. The target's trajectory is described by x₁ = 0.5t and x₂ = 2sin(0.5t+2). During the training phase, this trajectory serves as the ground truth to minimize the total distance between the robots and the target. In the test phase, the target location is unknown. The transition from training to test simulates a scenario in tracking processes where, despite initially having knowledge of the target location, the information regarding the target's position is lost from a certain moment onward.

As tracking problems in reality often occur under finite energy consumption, we introduce an inequality constraint on the total energy consumption of all agents: e ≤ e_max, where e represents the energy and e_max is the predetermined energy ceiling. This constraint is implemented during training through a regularization term as in Equation (7):

ℓ∨=ℓ(ω,ω∗)+max{emax,e}. (7) ℓ∨=ℓ(ω,ω*)+max{emax,e}. (7)

During training, if e>e_max, the term max{e_max, e} encourages a reduction in e; If e ≤ e_max, then this term does not affect e. Since the tracking occurs on a horizontal plane, potential energy is not considered; thus, e=∑rimiv2i/2 e=∑irmivi2/2 , where r is the total number of robots. In the simulations, we set each robot's mass as equal, with the total mass being 2 kg (therefore, e=∑riv2i=∑dj=0∑ri=0ẋ2ij e=∑irvi2=∑j=0d∑i=0rẋij2 ), and set e_max = 5,000 J.

4.2.1 Chasing trajectories with given target trajectory

Figures 3 and 4 respectively illustrate the trajectories of multiple robots and a chased target at different moments in time without energy constraints. By comparing Figures 3A and 4A, it can be observed that the chasing speed of the proposed HiDeS-3 NODE is significantly faster than the standard NODE. At the end of the tracking phase (e.g., at t = 20 s), both the standard NODE and HiDeS-3 NODE successfully reach the target. Examination of the various subfigures in Figure 4 reveals that the trajectories of the HiDeS-3 NODE exhibit typical characteristics of high-order dynamic systems similar to those seen with higher-order optimizers (Su et al., 2016; An et al., 2018) and proportional-integral-derivative (PID) controllers (Huba et al., 2023), such as rapid convergence and overshooting. This is attributed to HiDeS-3 NODE being a high-order dynamic system, as demonstrated by Equation 1. Figures 5 and 6 present a superficially similar performance between the standard NODE and the HiDeS NODE when subject to energy constraints. However, a distinct contrast emerges during the subsequent testing phase, which operates without a predefined target position.

Figure 3

Figure 3. Trajectories of multiple robots using standard NODE in target-chasing task without energy constraint. This figure depicts the paths of multiple robots (indicated by dots) as they follow a dynamically moving target (denoted by stars) after 1,000 epochs of training. (A) t = 2 s. (B) t = 4 s. (C) t = 5 s. (D) t = 10 s. (E) t = 15 s. (F) t = 20 s.

Figure 4

Figure 4. Trajectories of multiple robots using HiDeS-3 NODE (ours) in target-chasing task without energy constraint. This figure depicts the paths of multiple robots (indicated by dots) as they follow a dynamically moving target (denoted by stars) after 1,000 epochs of training. (A) t = 2 s. (B) t = 4 s. (C) t = 5 s. (D) t = 10 s. (E) t = 15 s. (F) t = 20 s.

Figure 5

Figure 5. Trajectories of multiple robots using standard NODE in target-chasing task with energy constraint. This figure depicts the paths of multiple robots (indicated by dots) as they follow a dynamically moving target (denoted by stars) after 1,000 epochs of training. (A) t = 2 s. (B) t = 4 s. (C) t = 5 s. (D) t = 10 s. (E) t = 15 s. (F) t = 20 s.

Figure 6

Figure 6. Trajectories of multiple robots using HiDeS-3 NODE (ours) in target-chasing task with energy constraint. This figure depicts the paths of multiple robots (indicated by dots) as they follow a dynamically moving target (denoted by stars) after 1,000 epochs of training. (A) t = 2 s. (B) t = 4 s. (C) t = 5 s. (D) t = 10 s. (E) t = 15 s. (F) t = 20 s.

4.2.2 Predicted trajectories without given target position

In the test phase, the target position is not given to the model. Figure 7 presents a comparison of the predicted trajectories of multi-robots in a target-chasing scenario where the target position is not provided. The comparison is between the trajectories generated by the standard NODE and the HiDeS NODE, with and without the imposition of energy constraints. The figure clearly demonstrates that the HiDeS NODE offers superior performance over the standard NODE. Specifically, the trajectories predicted by the standard NODE show significant deviations from the target (Figures 7A, C). In contrast, those predicted by the HiDeS NODE closely align with the target's trajectory. Although the HiDeS NODE with energy constraints shows slight deviations due to restricted velocity, it still significantly outperforms the standard NODE (Figures 7B, D). This suggests that the advantages of the HiDeS NODE may not solely be attributed to increased velocity but could also derive from additional information, such as curvature, which is inferred by the high-order supervised signals, as illustrated in the middle of Figure 1.

Figure 7

Figure 7. Trajectories of multiple robots in target-chasing task after losing information of target. This figure depicts the paths of multiple robots (indicated by dots) as they follow a dynamically moving target (denoted by stars). Solid curves denote training phase (trajectory of target is known), and dashed curves denote test phase (trajectory of target is unknown). (A) Standard NODE; without energy constraint. (B) HiDeS-3 NODE (ours); without energy constraint. (C) Standard NODE; with energy constraint. (D) HiDeS-3 NODE (ours); with energy constraint.

4.2.3 Energy cost

Figure 8 draws a parallel of the total energy expenditure of multiple robots engaged in a target-chasing task, contrasting the standard NODE with the HiDeS NODE. The figure shows that the implementation of energy constraints has a significant impact. Without energy constraints, both the standard NODE and the HiDeS NODE incur substantial energy costs after training, reaching up to 7,000 and 100,000 J, respectively (Figures 8A,B). However, with the application of energy constraints, both NODEs manage to keep the energy expenditure no more than 5,000 J after 1,000 epochs of training (Figures 8C, D).

Figure 8

Figure 8. Total energy cost of multiple robots in target-chasing task during 1,000 epochs of training. Energy constraint is set such that total energy cost ≤5,000 J. (A) standard NODE; without energy constraint. (B) HiDeS-3 NODE (ours); without energy constraint. (C) Standard NODE; with energy constraint. (D) HiDeS-3 NODE (ours); with energy constraint.

4.3 Modeling and predicting of opinion dynamics

In this section, simulations of NODEs on modeling and predicting opinion dynamics are conducted.

4.3.1 One dimension, multiple initial conditions

We respectively present experimental results conducted over 1,000 and 2,000 epochs in Figures 9, 10, and complemented by Table 3. The results indicate that our HiDeS NODEs surpass the standard NODE and SONODE in capturing the subtleties of opinion dynamics. Note that the SONODE and HiDeS-2 have the same hidden layer dimensions. This superiority is evident from the more accurate approximations of actual opinion dynamics during both training and testing phases (Figures 9 and 10) and lower auxiliary losses in consistent iterations (Table 3). For instance, Figure 9 shows that the prediction curves of standard NODE diverge from the ground truth traces in both training and testing since the trajectories from different initial conditions can not cross. While SONODE performs better than standard NODE in the learning stage, its predictive ability remains inferior in the testing phase. In contrast, our models' prediction curves closely align with the ground truth. Extending the epochs to 2,000 shows that both models perform better in the learning stage than under 1,000 epochs, with our models demonstrating remarkable superiority in the testing phase. This empirical observation aligns with our theoretical analyses that using higher-order derivatives as supervised signals enhances the predictive capacity. Furthermore, as Table 3 shows, the proposed HiDeS-3 NODE exhibits significantly lower training and testing losses compared to the baseline models. Although the HiDeS-2 NODE exhibits a slightly higher training loss than SONODE, its test loss is substantially lower, indicating superior generalization ability, a key goal of neural networks.

Figure 9

Figure 9. Visualization of learning results with standard NODE, SONODE, HiDeS-2 NODE, and HiDeS-3 NODE. All NODEs are trained for 1,000 epochs. (A) Standard NODE; opinion evolution; training. (B) Training loss w.r.t. epoch. (C) Standard NODE; opinion evolution; test. (D) SONODE; Opinion evolution; training. (E) Training loss w.r.t. epoch. (F) SONODE; opinion evolution; test. (G) HiDeS-2 NODE; opinion evolution; training. (H) Training loss w.r.t. epoch. (I) HiDeS-2 NODE; opinion evolution; test. (J) HiDeS-3 NODE; opinion evolution; training. (K) Training loss w.r.t. epoch. (L) HiDeS-3 NODE; opinion evolution; test.

Figure 10

Figure 10. Visualization of learning results with standard NODE, SONODE, HiDeS-2 NODE, and HiDeS-3 NODE. All NODEs are trained for 2,000 epochs. (A) Standard NODE; opinion evolution; training. (B) Training loss w.r.t. epoch. (C) Standard NODE; opinion evolution; test. (D) SONODE; Opinion evolution; training. (E) Training loss w.r.t. epoch. (F) SONODE; opinion evolution; test. (G) HiDeS-2 NODE; opinion evolution; training. (H) Training loss w.r.t. epoch. (I) HiDeS-2 NODE; opinion evolution; test. (J) HiDeS-3 NODE; opinion evolution; training. (K) Training loss w.r.t. epoch. (L) HiDeS-3 NODE; opinion evolution; test.

Table 3

Table 3. Comparisons of training and test losses among SOTA NODEs and HiDeS NODE.

4.3.2 Multiple dimension, one initial condition

Figure 11 presents modeling and predicting results on diverse types of individuals' intra- and inter-group interactions for opinion dynamics. The top row (Figures 11A–D) represents ground truth values, while the bottom row (Figures 11E, F) shows predictions generated by the HiDeS-3 NODE. Each subfigure illustrates different combinations of consensus and dissensus within and between groups, highlighting the model's performance in capturing extensive opinion dynamics.

Figure 11

Figure 11. Opinion dynamics modeling and predicting using HiDeS-3 NODE for different consensus and dissensus scenarios. (A) Intra-group consensus and inter-group consensus, ground truth. (B) Intra-group consensus and inter-group dissensus, ground truth. (C) Intra-group dissensus and inter-group consensus, ground truth. (D) Intra-group dissensus and inter-group dissensus, ground truth. (E) Intra-group consensus and inter-group consensus, predicted by HiDeS-3 NODE. (F) Intra-group consensus and inter-group dissensus, predicted by HiDeS-3 NODE. (G) Intra-group dissensus and inter-group consensus, predicted by HiDeS-3 NODE. (H) Intra-group dissensus and inter-group dissensus, predicted by HiDeS-3 NODE.

4.3.3 Training dynamics

The learning results of opinion evolution are shown in Figure 12 with the increase of epochs during a 2,000-epoch training. It is clear that the learning results get more accurate and fine-grained with the increasing epochs. Specifically, the predicted results show significant deviation from the ground truth under 1–200 epochs (Figures 12B, C), while the more granular learning results are presented with the increase of epochs (Figures 12D–H), achieving more accurate predictions for opinion evolution.

Figure 12

Figure 12. Opinion dynamics learning results during 2,000-epoch training. (A) Epoch = 1. (B) Epoch = 100. (C) Epoch = 200. (D) Epoch = 500. (E) Epoch = 1,000. (F) Epoch = 1,100. (G) Epoch = 2,000. (H) Ground truth.

5 Discussion

This paper has proposed the HiDeS NODE, a higher-order-derivative-supervised NODE, as a novel approach for modeling and predicting complex dynamics in multi-robot systems and opinion dynamics. This framework excels in capturing interactions among higher-order derivatives and the state vector, significantly enhancing modeling precision over existing NODE methodologies. The introduction of higher-order derivatives as supervised signals in the HiDeS NODE brings a superior predicting ability. Applications of the HiDeS NODE in multi-robot systems and opinion dynamics have demonstrated its effectiveness. To our knowledge, this is the first initiative that introduces NODEs into multi-robot systems and opinion dynamics. Applying the HiDeS NODE to these fields opens new avenues for broader applications in various intricate and dynamic systems.

The broader impact of the HiDeS NODE extends into numerous fields where dynamic systems play a crucial role. For instance, it has the potential to offer refined predictions of climate change effects or pollution dispersion. In healthcare, the HiDeS NODE could lead to breakthroughs in understanding the dynamics of disease spread or patient response to treatments, enabling personalized medicine. The adaptability and advanced modeling capabilities of the HiDeS NODE position it as a versatile tool capable of addressing complex problems across various domains.

Despite its potential, the HiDeS NODE faces limitations such as computational demands, particularly as the order of derivatives increases, making real-time applications challenging. The model's accuracy is heavily reliant on the quality and quantity of data, which can be a significant constraint in environments where data is sparse or noisy. Addressing these challenges will be essential for the HiDeS NODE's successful application across different fields.

A valuable future direction is to utilize real-world data to validate our model's performance in practical scenarios. Additionally, enhancing the robustness of the HiDeS NODE to noisy data presents a promising direction for future research.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

MLi: Conceptualization, Formal analysis, Methodology, Software, Writing – original draft. WB: Data curation, Funding acquisition, Resources, Supervision, Validation, Visualization, Writing – review & editing. LC: Funding acquisition, Project administration, Resources, Software, Validation, Visualization, Writing – review & editing. MLiu: Data curation, Formal analysis, Investigation, Resources, Supervision, Validation, Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The handling editor LJ declared a shared affiliation with the author MLi at the time of review.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

An, W., Wang, H., Sun, Q., Xu, J., Dai, Q., and Zhang, L. (2018). “A PID controller approach for stochastic optimization of deep networks,” in Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (Los Alamitos, CA: IEEE Computer Society), 8522–8531.

Anil Kumar, N., Patrick, S., Hong, W., and Hur, P. (2022). Control framework for sloped walking with a powered transfemoral prosthesis. Front. Neurorobot. 15:790060. doi: 10.3389/fnbot.2021.790060

Ao, Y., and Jia, Y. (2023). Agents attraction competition in an extended Friedkin-Johnsen social network. IEEE Transact. Cont. Netw. Syst. 10, 1100–1112. doi: 10.1109/TCNS.2022.3220709

Aydemir, G., Akan, A. K., and Güney, F. (2023). “Adapt: efficient multi-agent trajectory prediction with adaptation,” in Proceedings of the IEEE/CVF International Conference on Computer Vision (Los Alamitos, CA: IEEE Computer Society), 8295–8305.

Baumann, F., Lorenz-Spreen, P., Sokolov, I. M., and Starnini, M. (2020). Modeling echo chambers and polarization dynamics in social networks. Phys. Rev. Lett. 124:048301. doi: 10.1103/PhysRevLett.124.048301

Chen, R. T., Rubanova, Y., Bettencourt, J., and Duvenaud, D. K. (2018). Neural ordinary differential equations. Adv. Neural Inf. Process. Syst. 31, 6572–6583. doi: 10.5555/3327757.3327764

Cui, W., Zhang, H., Chu, H., Hu, P., and Li, Y. (2023). On robustness of neural ODEs image classifiers. Inf. Sci. 632, 576–593. doi: 10.1016/j.ins.2023.03.049

De, A., Bhattacharya, S., Bhattacharya, P., Ganguly, N., and Chakrabarti, S. (2014). “Learning a linear influence model from transient opinion dynamics,” in Proceedings of the 23rd ACM International Conference on Conference on Information and Knowledge Management (New York, NY: ACM), 401–410.

De, A., Valera, I., Ganguly, N., Bhattacharya, S., and Gomez Rodriguez, M. (2016). Learning and forecasting opinion dynamics in social networks. Adv. Neural Inf. Process. Syst. 29, 397–405. doi: 10.5555/3157096.3157141

Dozat, T. (2016). “Incorporating Nesterov momentum into Adam,” in International Conference on Learning Representations (OpenReview.net).

Dupont, E., Doucet, A., and Teh, Y. W. (2019). Augmented neural ODEs. Adv. Neural Inf. Process. Syst. 32, 3140–3150. doi: 10.5555/3454287.3454569

Granha, M. F., Vilela, A. L., Wang, C., Nelson, K. P., and Stanley, H. E. (2022). Opinion dynamics in financial markets via random networks. Proc. Nat. Acad. Sci. U. S. A. 119:e2201573119. doi: 10.1073/pnas.2201573119

Hichri, B., Gallala, A., Giovannini, F., and Kedziora, S. (2022). Mobile robots path planning and mobile multirobots control: a review. Robotica 40, 4257–4270. doi: 10.1017/S0263574722000893

Hua, C., Cao, X., Liao, B., and Li, S. (2023). Advances on intelligent algorithms for scientific computing: an overview. Front. Neurorobot. 17:1190977. doi: 10.3389/fnbot.2023.1190977

Huba, M., Vrancic, D., and Bistak, P. (2023). Series PID control with higher-order derivatives for processes approximated by IPDT models. IEEE Transact. Automat. Sci. Eng. doi: 10.1109/TASE.2023.3296201. [Epub ahead of print].

Jin, J., Zhao, L., Chen, L., and Chen, W. (2022). A robust zeroing neural network and its applications to dynamic complex matrix equation solving and robotic manipulator trajectory tracking. Front. Neurorobot. 16:1065256. doi: 10.3389/fnbot.2022.1065256

Jin, L., Liu, L., Wang, X., Shang, M., and Wang, F.-Y. (2024). Physical-informed neural network for MPC-based trajectory tracking of vehicles with noise considered. IEEE Transact. Intell. Vehicl. doi: 10.1109/TIV.2024.3358229. [Epub ahead of print].

Károly, A. I., Galambos, P., Kuti, J., and Rudas, I. J. (2021). Deep learning in robotics: survey on model structures and training strategies. IEEE Transact. Syst. Man Cybernet. Syst. 51, 266–279. doi: 10.1109/TSMC.2020.3018325

Khalili, M., Zhang, X., Cao, Y., Polycarpou, M. M., and Parisini, T. (2020). Distributed fault-tolerant control of multiagent systems: an adaptive learning approach. IEEE Transact. Neural Netw. Learn. Syst. 31, 420–432. doi: 10.1109/TNNLS.2019.2904277

Kidger, P. (2021). On Neural Differential Equations (PhD thesis). Oxford: University of Oxford.

Kolarijani, M. A. S., Proskurnikov, A. V., and Esfahani, P. M. (2021). Macroscopic noisy bounded confidence models with distributed radical opinions. IEEE Trans. Automat. Contr. 66, 1174–1189. doi: 10.1109/TAC.2020.2994284

Kulkarni, B., Agarwal, S., De, A., Bhattacharya, S., and Ganguly, N. (2017). “SLANT+: a nonlinear model for opinion dynamics in social networks,” in IEEE International Conference on Data Mining (Piscataway, NJ: IEEE), 931–936.

Kwa, H. L., Leong Kit, J., and Bouffanais, R. (2022). Balancing collective exploration and exploitation in multi-agent and multi-robot systems: a review. Front. Robot. AI 8:771520. doi: 10.3389/frobt.2021.771520

Li, Z., Liu, F., Yang, W., Peng, S., and Zhou, J. (2022). A survey of convolutional neural networks: analysis, applications, and prospects. IEEE Transact. Neural Netw. Learn. Syst. 33, 6999–7019. doi: 10.1109/TNNLS.2021.3084827

Liu, L., Wang, X., Yang, X., Liu, H., Li, J., and Wang, P. (2023). Path planning techniques for mobile robots: review and prospect. Expert Syst. Appl. 227:120254. doi: 10.1016/j.eswa.2023.120254

Liufu, Y., Jin, L., Shang, M., Wang, X., and Wang, F.-Y. (2024). ACP-incorporated perturbation-resistant neural dynamics controller for autonomous vehicles. IEEE Transact. Intell. Vehicl. doi: 10.1109/TIV.2023.3348632. [Epub ahead of print].

Monti, C., De Francisci Morales, G., and Bonchi, F. (2020). “Learning opinion dynamics from social traces,” in Proceedings of the ACM SIGKDD International Conference on Knowledge Discovery & Data Mining (New York, NY: ACM), 764–773.

Norcliffe, A., Bodnar, C., Day, B., Simidjievski, N., and Liò, P. (2020). On second order behaviour in augmented neural ODEs. Adv. Neural Inf. Process. Syst. 33, 5911–5921. doi: 10.5555/3495724.3496220

Okawa, M., and Iwata, T. (2022). “Predicting opinion dynamics via sociologically-informed neural networks,” in Proceedings of the ACM SIGKDD Conference on Knowledge Discovery and Data Mining (New York, NY: ACM), 1306–1316.

Peng, Y., Zhao, Y., and Hu, J. (2023). On the role of community structure in evolution of opinion formation: a new bounded confidence opinion dynamics. Inf. Sci. 621, 672–690. doi: 10.1016/j.ins.2022.11.101

Pierpaoli, P., Li, A., Srinivasan, M., Cai, X., Coogan, S., and Egerstedt, M. (2021). A sequential composition framework for coordinating multirobot behaviors. IEEE Transact. Robot. 37, 864–876. doi: 10.1109/TRO.2020.3036628

Ruiz-Balet, D., and Zuazua, E. (2023). Neural ODE control for classification, approximation, and transport. SIAM Rev. 65, 735–773. doi: 10.1137/21M1411433

Sander, M. E., Ablin, P., Blondel, M., and Peyré, G. (2021). “Momentum residual neural networks,” in International Conference on Machine Learning (London: PLMR), 9276–9287.

Su, W., Boyd, S., and Candès, E. J. (2016). A differential equation for modeling Nesterov's accelerated gradient method: theory and insights. J. Mach. Learn. Res. 17, 1–43. doi: 10.5555/2946645.3053435

Wang, T., Wang, X., Lee-Sarwar, K. A., Litonjua, A. A., Weiss, S. T., and Sun, Y. (2023). Predicting metabolomic profiles from microbial composition through neural ordinary differential equations. Nat. Mach. Intell. 5, 284–293. doi: 10.1038/s42256-023-00627-3

Weinan, E. (2017). A proposal on machine learning via dynamical systems. Commun. Math. Stat. 1, 1–11. doi: 10.1007/s40304-017-0103-z

Wu, Z., Zhou, Q., Dong, Y., Xu, J., Altalhi, A. H., and Herrera, F. (2023). Mixed opinion dynamics based on degroot model and hegselmann-krause model in social networks. IEEE Transact. Syst. Man Cybernet. Syst. 53, 296–308. doi: 10.1109/TSMC.2022.3178230

Xia, H., Suliafu, V., Ji, H., Nguyen, T., Bertozzi, A., Osher, S., et al. (2021). Heavy ball neural ordinary differential equations. Adv. Neural Inf. Process. Syst. 34, 18646–18659. doi: 10.48550/arXiv.2110.04840

Xu, S., Xu, T., Li, D., Yang, C., Huang, C., and Wu, X. (2023). A robot motion learning method using broad learning system verified by small-scale fish-like robot. IEEE Trans. Cybern. 53, 6053–6065. doi: 10.1109/TCYB.2023.3269773

Yu, X., and Chen, T. (2023). Distributed iterative learning control of nonlinear multiagent systems using controller-based dynamic linearization method. IEEE Transact. Cybernet. doi: 10.1109/TCYB.2023.3281479. [Epub ahead of print].

Zhang, Y., Li, P., Xu, C., Peng, X., and Qiao, R. (2023). Investigating the effects of a fractional operator on the evolution of the enso model: bifurcations, stability and numerical analysis. Fract. Fract. 7:602. doi: 10.3390/fractalfract7080602

Zhu, L., He, Y., and Zhou, D. (2020). Neural opinion dynamics model for the prediction of user-level stance dynamics. Inf. Process. Manag. 57:102031. doi: 10.1016/j.ipm.2019.03.010