
Preprint
performance: the left panel shows the learned partitions and the right panel shows the ground truth
target function (solid black) and the prediction (dashed cian). In each partition, a set of monomials
with the maximal degree 2 is fitted optimally by solving local linear least-squares problems.
2.2 PARTITION-OF-UNITY-BASED NEURAL ORDINARY DIFFERENTIAL EQUATIONS
Now, we introduce the proposed partition-of-unity-based neural ordinary differential equations,
where the model parameters are represented as a POUNet: Θ(s)∈RnΘ:
Θ(s;α, π) =
npart
X
i=1
φi(s;π)pi(s) =
npart
X
i=1
φi(s;π)
npoly
X
j=1
αi,j ψj(s),(1)
where sdenotes a set of variables whose domains are expected to have a set of partitions (e.g.,
scan be a depth variable in depth-continuous neural network architectures), φi(s;π)∈Rde-
notes a partition of unity network, parameterized by π,ψj(s)∈Rdenotes a polynomial basis,
and α·,j ∈RnΘdenote the polynomial coefficients. Thus, collectively, there is a set of parameters
α= (α1, . . . , αnpart )with αi= [αi,1· · · , αi,npoly ]∈Rnθ×npoly . In the following, we present a couple
example cases of the types of the variables s.
Temporally varying dynamics / depth variance As in the typical settings of NODEs, when an
MLP is considered to parameterize the velocity function, f(·; Θ), the model parameters can be rep-
resented as a set of constant-valued variables, Θ = {(W`,b`)}L
`=1, where W`and b`denote weights
and biases of the `-th layer. As opposed to the depth-invariant NODE parameters Θ, POUNODEs
represent depth-variant NODEs (or non-autonomous dynamical systems) by setting the model pa-
rameters as
Θ(t) = {(W`(t),b`(t))}L
`=1,
where tdenotes the time variable or the depth of the neural network and represent, and by represent-
ing Θ(t)as a POUNet as in Eq. (1) with s=t.
Spatially varying dynamics Another example dynamical systems that can be represented by
POUNODEs is a class of dynamical systems whose dynamics modes are defined differently on
different spatial regions. In this case, the model parameters can be set as spatially-varying ones:
Θ(x) = {(W`(x),b`(x))}L
`=1.
and can be represented as a POUNet as in Eq. (1) with s=x.
Remark 2.1. Although not numerically tested in this study, the idea of representing the evolution
of model parameters via POUNets can be applied to different neural network architectures, e.g.,
POU-Recurrent Neural Networks (POU-RNNs).
3 USE CASES
This section exhibits example use cases where the benefits of using POUNODE can be pronounced.
All implementations are based on PYTORCH (Paszke et al., 2019) and the TORCHDIFFEQ library
(Chen et al., 2018) for the NODEs capability.
For all following experiments, we consider a POUNet, Φ = {φi}npart
i=1, based on a radial basis function
(RBF) network (Broomhead & Lowe, 1988; Billings & Zheng, 1995); for each partition, there is an
associated RBF layer, defined by its center and shape parameter, and then the output of the RBF
layers is normalized to satisfy the partition-of-unity property (refer to Appendix for more details).
3.1 SYSTEM IDENTIFICATION OF A HYBRID SYSTEM
As a first set of use cases, we apply POUNODEs for data-driven dynamics modeling. In particular,
we aim to learn a dynamics model for a hybrid system, where the different dynamics models are
mixed in a single system: a system consisting of multiple smooth dynamical flows (SDFs), each
of which is interrupted by sudden changes (e.g., jump discontinuities or distributional shifts) (Van
Der Schaft & Schumacher, 2000).
3