Spatially-invariant opinion dynamics on the circle

Abstract

We propose a nonlinear opinion dynamics model for an agent making decisions about a continuous distribution of options in the presence of distributed input. Inspired by perceptual decision-making, we develop the theory for options distributed on the circle, representing, e.g., the set of possible heading directions in planar robotic navigation problems. Interactions among options are encoded through a spatially invariant kernel. We design the kernel to ensure that decision-making is robust, in the sense that only a finite subset of options can be favored over the continuum. We prove the spatial invariance of the model linearization and use this result to prove an opinion-forming bifurcation in the model with zero input. We then use space and time frequency domain analysis of the model linearization to infer the ultra-sensitive input-output behavior of the nonlinear dynamics with input. We illustrate our model’s versatility with a robotic navigation example in crowded spaces.

I INTRODUCTION

In perceptual decision-making, animals use sensory information to respond to their environment. Both the sensory information and the perceptual response are spatially distributed. For example, visual information, like object positions, is distributed both in the physical world and in its neural representation in the visual cortex [1]. Similarly, in robotics, an agent can use a distributed representation of its visual field and the objects in it to drive decisions, e.g., to navigate crowded spaces or track targets as in  [2, 3, 4, 5].

Motivated by perceptual decision-making, we propose a distributed nonlinear opinion dynamics model for an agent making decisions about a continuous distribution of options in the presence of continuously distributed input. We study option and input spaces defined by the unit circle, representing both a one-dimensional visual field and the set of possible heading directions. Our model generalizes the nonlinear opinion dynamics (NOD) of [6, 7] by extending it from a finite set to a continuous distribution of options. NOD has been used in a variety of applications, such as robotics and automated driving [8, 9, 10, 11, 12, 13]. These works leverage the fast-and-flexible decision-making of NOD resulting from a bifurcation in the dynamics [7].

The distributed model we present here inherits the same features as a result of a bifurcation in its infinite-dimensional dynamics. The proposed dynamics also allow for increased scalability and use in more complex applications. For instance, in an obstacle avoidance scenario as the one in [8], the representation of objects in the distributed NOD does not require prior knowledge of the number of objects to be avoided. In the finite-dimensional NOD this prior information is, in principle, fundamental to the formulation.

Our model and the bifurcation behavior we prove is closely related to the spatiotemporal neural field models of mathematical neuroscience [14, 15, 1, 16], which consider spatially-invariant integro-differential equations with a convolution kernel that captures excitatory and inhibitory interactions between different regions of the neural field. Response to input is considered in [16, 17, 18] but only for specific classes of inputs. Our approach differs from those works by using a spatiotemporal transfer function approach, which allows us to predict the response of the nonlinear model to arbitrary inputs from the response of its linearization.

For illustration, we consider a robot navigation scenario where the visual input corresponds to the locations of (possibly many) obstacle-free “gaps”. This scenario is inspired by classical empirical approaches exploring the use of neural fields for robotic navigation [2, 3, 4, 5]. Our model provides a theoretical foundation to such approaches that can be used for analysis, design, and control in more general scenarios.

Our contributions are as follows. First, we propose a nonlinear opinion dynamics model for an agent making decisions about a continuous distribution of options and in the presence of continuously distributed inputs. Second, we prove the spatial invariance of the model linearization. Third, we use spatial-invariance of the linearized dynamics to prove the existence of an opinion forming bifurcation for the model with zero input. Fourth, we use space and time frequency domain analysis of the model linearization to infer the input-output behavior of the nonlinear dynamics to arbitrary inputs. Fifth, we propose a framework for designing kernels for an application of our model to robotic navigation.

Notation and mathematical background are in Section II. We introduce the nonlinear opinion dynamics for a continuous distribution of options in Section III. We prove the spatial invariance of the model linearization in Section IV. In Section V we leverage spatial-invariance of the linearized dynamics to prove an opinion forming bifurcation in the model with zero input. In Section VI, we use spatial and time frequency analysis of the model linearization to infer the model’s input-output behavior. We propose a frequency-domain kernel design approach and illustrate the benefits of our approach for robotic navigation in Section VII. A discussion and final remarks are provided in Section VIII.

II MATHEMATICAL PRELIMINARIES

We denote the set of integer values as ${\mathbb{Z}}$, the set of non-negative integer values as ${\mathbb{N}}$, the set of real numbers as ${\mathbb{R}}$, and the set of complex numbers as ${\mathbb{C}}$. The unit circle is denoted by $\mathbb{S}^{1}$, i.e., $\mathbb{S}^{1}={\mathbb{R}}/{\mathbb{Z}}$. For $a,b\in{\mathbb{R}}$, the notation $a\nearrow b$ indicates the limit $a\to b$ with $a<b$. For a complex number $s=\sigma+i\omega$, the real and imaginary parts are denoted as $\Re(s)$ and $\Im(s)$, respectively. We represent the complex conjugate as $\bar{s}=\sigma-i\omega$, the modulus as $|s|=\sqrt{s\bar{s}}$ and the argument as $\arg(s)=\lim_{n\to\infty}n\Im(\sqrt[n]{s/|s|})$ for $n\in{\mathbb{N}}-\{0\}$.

The Hilbert space of square-integrable real functions on $\mathbb{S}^{1}$ is denoted by $L_{2}(\mathbb{S}^{1})$. Given $v,w\in L_{2}(\mathbb{S}^{1})$, their inner product is $\langle v,w\rangle=\int_{\mathbb{S}^{1}}v(\theta)w(\theta)d\theta$. The induced norm is $||v||=\langle v,v\rangle^{1/2}$.

We denote operators with capital roman letters. Let $A:L_{2}(\mathbb{S}^{1})\to L_{2}(\mathbb{S}^{1})$ be a linear operator. We let the set $\text{Sp}(A)=\{\lambda_{k}\}$ denote the point spectrum of $A$, if it is not empty. Each eigenvalue $\lambda_{k}\in\text{Sp}(A)$ satisfies $Av_{k}(\theta)=\lambda_{k}v_{k}(\theta)$, where $v_{k}\in L_{2}(\mathbb{S}^{1})$ denotes the eigenfunction corresponding to $\lambda_{k}$. We denote $\lambda_{\max}=\arg\!\max\{\Re(\lambda_{k})\}$ as the leading eigenvalue of $A$, and its corresponding eigenfunction, as $v_{\max}\in L_{2}(\mathbb{S}^{1})$, the leading mode.

Definition 1 (Differential operator):

Let $F:L_{2}(\mathbb{S}^{1})\to L_{2}(\mathbb{S}^{1})$ be a nonlinear operator. The differential of $F$ in the direction of $z$ at a point $z^{*}$, is

$$\textstyle A_{F}=D_{z}F(z^{*}):=\lim_{\epsilon\to 0}\frac{1}{\epsilon}\big{(}F(\epsilon z+z^{*})-F(z^{*})\big{)},$$

provided that the limit exists. If this limit exists, then $F$ is said to be differentiable in the direction $z$ at point ${z}^{*}$.

Definition 2 (Multiplication Operator):

A multiplication operator $M$ is defined by $[Mh](x):=M(x)h(x)$, where $h$ is in the domain of $M$. Multiplication operators are the infinite-dimensional equivalent of diagonal matrices.

Definition 3 (Spatial shift operator [19], [20]):

The spatial shift operator denoted by $T_{\psi}:L^{2}(\mathbb{S}^{1})\to L^{2}(\mathbb{S}^{1})$ is defined as $h(\theta)\mapsto[T_{\psi}h](\theta):=h(\theta-\psi)$ for $\psi\in\mathbb{S}^{1}$ and $h\in L^{2}(\mathbb{S}^{1})$.

Definition 4 (Spatially invariant operator [19], [20]):

An operator $F$ is spatially invariant if $T_{\psi}F=FT_{\psi}$.

We mainly work with a special class of spatially invariant linear operators, namely, spatial convolution operators

$$\textstyle[Az](\theta):=\int_{\mathbb{S}^{1}}W(\theta-\phi)z(\phi)d\phi,$$ (1)

where the convolution kernel $W:\mathbb{S}^{1}\to{\mathbb{R}}$.

Definition 5 (Spatially Invariant Linear System [19], [20]):

Consider a spatiotemporal input-output linear system. Let $u(\cdot,t),z(\cdot,t)\in L_{2}(\mathbb{S}^{1})$ be the scalar-valued input and output functions at time $t\in\mathbb{R}_{\geq 0}$, respectively. Let $\theta\in\mathbb{S}^{1}$ be the spatial coordinate. A linear system of the form

$$\textstyle\frac{\partial z}{\partial t}(\theta,t)=[Az](\theta,t)+[Bu](\theta,t),$$ (2)

is spatially invariant if the linear operators $A$, $B$ are spatially invariant.

Definition 6 (Spatial Fourier transform [19], [20]):

Let $f,g:\mathbb{S}^{1}\times{\mathbb{R}}_{\geq 0}$ be spatiotemporal fields with spatial and time coordinates $\theta\in\mathbb{S}^{1}$ and $t\in{\mathbb{R}}_{\geq 0}$. Suppose $f(\cdot,t),g(\cdot,t)\in L^{2}(\mathbb{S}^{1})$ for all $t\in{\mathbb{R}}_{\geq 0}$. The spatial Fourier transform maps $f(\theta,t)$ into its spatial Fourier coefficients

$$\textstyle\hat{f}(k,t):=\int_{\mathbb{S}^{1}}f(\theta,t)e^{-i2\pi k\theta}d\theta,$$ (3)

where $k\in{\mathbb{Z}}$ is the spatial frequency.

The spatial Fourier transform is a coordinate transformation that expresses $f(\theta,t)$ in terms of the spatial Fourier modes $\eta_{k}(\theta)=e^{i2\pi k\theta}$, i.e., the Fourier basis on $\mathbb{S}^{1}$, and the Fourier coefficients $\hat{f}(k,t)$. The inverse spatial Fourier transform can be used to recover $f(\theta,t)$ from its Fourier coefficients $\hat{f}(k,t)$ :

$$\textstyle f(\theta,t)=\sum_{k\in\mathbb{Z}}\hat{f}(k,t)e^{i2\pi k\theta}.$$ (4)

Parseval’s Identity [21] ensures that $\langle\hat{f},\hat{g}\rangle=\langle f,g\rangle$. The spatial Fourier transform operator is denoted by $\mathcal{F}(\cdot)$, and the inverse spatial Fourier transform operator by $\mathcal{F}^{-1}(\cdot)$.

The spatial Fourier transform (3) diagonalizes convolution operators [19], i.e., if $A$ is a convolution operator (1), then $\widehat{[Ah]}(k):=\hat{W}(k)\hat{h}(k)$, where $\hat{W}$ is the Fourier transform of the kernal of $A$. Thus, $A$ is mapped by $\mathcal{F}$ into a multiplication operator over the spatial frequency $k\in{\mathbb{Z}}$. For linear systems of the form (2), if $A$ and $B$ are convolution operators (1),

$$\frac{\partial\hat{z}(k,t)}{\partial t}=\hat{W}_{A}(k)\hat{z}(k,t)+\hat{W}_{B}(k)\hat{u}(k,t),$$ (5)

where $\hat{W}_{A}$ and $\hat{W}_{B}$ are the Fourier transforms of the kernels of $A$ and $B$, respectively. Following [19], we refer to (5) as the diagonalization of (2).

Definition 7 (Temporal Laplace Transform):

Let $z:\mathbb{S}^{1}\times{\mathbb{R}}_{\geq 0}$ be a spatiotemporal field with spatial and time coordinates $\theta\in\mathbb{S}^{1}$ and $t\in{\mathbb{R}}_{\geq 0}$, respectively. Then, the temporal Laplace transform maps $z(\theta,t)$ into

$$[\mathcal{L}z](\theta,s)=\int_{0}^{\infty}z(\theta,t)e^{-st}dt,$$ (6)

where $s\in{\mathbb{C}}$, whenever the integral exists.

III OPINION DYNAMICS ON THE CIRCLE

We propose a nonlinear opinion dynamics model for an agent making decisions about a continuous distribution of options on the circle. For every option $\theta\in\mathbb{S}^{1}$, $z(\theta,t)\in{\mathbb{R}}$ is the opinion of the agent for option $\theta$ at time $t$, where the more positive (negative) $z(\theta,t)$ is the more the agent favors (disfavors) option $\theta$. When $z(\theta,t)=0$, the agent is neutral about option $\theta$. The relationship between each option is encoded by the Lipschitz continuous kernel $W:\mathbb{S}^{1}\to{\mathbb{R}}$. A positive (negative) value of $W(\theta-\phi)$ corresponds to excitatory (inhibitory) interactions between the options $\theta$ and $\phi$. The opinion $z(\theta,t)$ evolves according to

$$\textstyle\tau\frac{{\partial}z}{\partial t}\!(\theta,t)\! = \textstyle\!-z(\theta,t)\!+\!\alpha\!\!\int\limits_{\mathbb{S}^{1}}\!W\!(\theta\!-\!\phi)S(z(\phi,t))d\phi\!+\!u(\theta,t)$$ (7)

$$= [Gz](\theta,t)+u(\theta,t),$$

where $u(\theta,t)\in{\mathbb{R}}$ is the input and $\tau\in{\mathbb{R}}_{>0}$ is the characteristic timescale. The parameter $\alpha\in{\mathbb{R}}_{>0}$ is the attention paid to option interactions, i.e., $\alpha$ models the agent commitment to forming strong opinions. The nonlinear nature of (7) comes from $S:{\mathbb{R}}\to{\mathbb{R}}$, a saturating function with $S(0)=0$, $S^{{}^{\prime}}(0)=1$.

IV SPECTRAL ANALYSIS OF LINEARIZATION

We study the spectrum of the linearization of (7) at the neutral equilibrium $z(\theta,t)=0$, $\forall\theta\in\mathbb{S}^{1}$. We first prove spatial invariance, a property that enables the linearized system to be diagonalized. Using the diagonalization, we then compute the eigenvalues and eigenfunctions of the linearized system and prove their relationship with the Fourier coefficients of the kernel and the spatial Fourier modes.

Lemma 1 (Spatial invariance of the model linearization):

Define the nonlinear operator in (7) as

$$\textstyle[Fz](\theta,t)=\int_{\mathbb{S}^{1}}W(\theta-\phi)S(z(\phi,t)).$$ (8)

The differential of $F$ in the direction $z$ at $z(\theta,t)=0$ is

$$\textstyle[A_{F}z](\theta,t)\!=\\!=\!\int_{\mathbb{S}^{1}}W(\theta\!-\!\phi)z(\phi,t)d\phi.$$ (9)

The linearization of (7) at the neutral equilibrium $z(\theta,t)=0$,

$$\begin{split}\textstyle\frac{{\partial}z}{\partial t}(\theta,t)\!&=\!\textstyle\frac{1}{\tau}\Big{(}\!-\!z(\theta,t)+\!\alpha[A_{F}z](\theta,t)\!+\!u(\theta,t)\!\Big{)}\!\\ &\textstyle=\!\frac{1}{\tau}([A_{G}z](\theta,t)+u(\theta,t)),\end{split}$$ (10)

is a spatially invariant system in the sense of Definition 2.

• Consider the expansion $S(\epsilon z)=\sum_{n=0}^{\infty}\frac{1}{n!}(\epsilon z)^{n}S^{(n)}(0)$. Then, we can express $D_{z}F(0)=\lim_{\epsilon\to 0}\frac{\epsilon\int_{\mathbb{S}^{1}}W(\theta-\phi)S^{\prime}(0)z(\phi)d\phi+\mathcal{O}(\epsilon^{2})}{\epsilon}$, where $\mathcal{O}(\epsilon^{2})$ denotes higher order terms in $\epsilon$. As $\epsilon\to 0$, the higher order terms vanish and we are left with (9). Note that $A_{F}$ is a spatial convolution operator, which is spatially invariant. Then, by linearity so is $A_{G}$. Therefore, by Definition 2, (10) is a spatially invariant system. ∎

As a consequence of Lemma 1, we can diagonalize the model linearization (10). Since (9) is a convolution operator, we use (5) to get

$$\textstyle\frac{{\partial}\hat{z}}{\partial t}(k,t)=\frac{1}{\tau}\big{(}-1+\alpha\hat{W}(k)\big{)}\hat{z}(k,t)+\frac{1}{\tau}\hat{u}(k,t).$$ (11)

Lemma 2 (Eigenvalues and eigenfunctions of the linearized system):

The eigenvalues $\lambda_{k}\in\text{Sp}(A_{G})$ of the linearized system (10) can be computed as

$$\textstyle\lambda_{k}=\frac{1}{\tau}(-1+\alpha\hat{W}(k)).$$ (12)

for $k\in{\mathbb{Z}}$. The corresponding eigenfunctions are the spatial Fourier modes $\eta_{k}(\theta)=e^{i2\pi k\theta}$.

• The form of the eigenvalues follows directly from the diagonalization (11) of the linearized dynamics (10). The eigenfunctions are the spatial Fourier modes because they form the basis of the Fourier transformation that is used to diagonalize the system. ∎

Lemma 12 reaffirms that, because of spatial invariance, the spatial Fourier modes are the eigenfunctions of the model linearization for any kernel design, provided it is spatially-invariant. Since the Fourier coefficients of the kernel determine the eigenvalues associated with each mode, they dictate which modes dominate. More precisely, if all modes are stable, i.e., $\Re(\lambda_{k})<0$ for all $k\in{\mathbb{Z}}$, spatiotemporal inputs $u(\theta,t)$ will be predominantly amplified along the Fourier modes with largest $\Re(\lambda_{k})$ as detailed in Section VI.

When the leading modes becomes unstable, that is, when the real part of their eigenvalues change from negative to positive through, e.g., an increase of the attention parameter $\alpha$, the nonlinear model (7) undergoes a bifurcation that enables robust opinion formation even in the absence of inputs. The leading Fourier modes determine the number of maxima of the stable steady-state opinion patterns emerging at the bifurcation, as detailed in the next section.

V OPINION-FORMING BIFURCATIONS

We revisit the results presented in [15], [16] for (7) with zero input. We prove that (7) undergoes a bifurcation and compute the bifurcation point. A local bifurcation occurs when the number and/or stability of the equilibrium solutions changes due to one or more eigenvalues of the model linearization crossing the imaginary axis as a parameter is varied. The state and parameter value at which this occurs is the bifurcation point. We study how opinion patterns emerge and the role of kernel $W$ and show a bistability that enables rapid formation of strong opinions. We make the following assumption to ensure the eigenvalues of (10) are real.

Assumption 1 (Symmetric kernels):

The kernel $W$ in (7)) is symmetric, i.e. $W(\psi)=W(-\psi)$. In particular, its Fourier coefficients $\hat{W}$ are real and $\hat{W}(k)=\hat{W}(-k)$.

Lemma 3 (Existence of a bifurcation point at the neutral equilibrium):

Consider (7). Let Assumption 1 hold. Let $\pm k_{\max}=\arg\max_{k\in{\mathbb{Z}}}\hat{W}(k)$, $k_{max}>0$, denote the spatial frequency corresponding to the largest $\hat{W}(k)$. Then, system (7) undergoes a bifurcation at the neutral equilibrium $z(\theta,t)=0$ and $\alpha^{*}=\frac{1}{\hat{W}(k_{\max})}$. In particular, for $0<\alpha<\alpha^{*}$ the neutral equilibrium is locally asymptotically stable and for $\alpha>\alpha^{*}$ the neutral equilibrium is unstable. The bifurcation branches emerging at bifurcation for $\alpha=\alpha^{*}$ are tangent to the subspace of $L^{2}(\mathbb{S}^{1})$ generated by $\cos(2\pi k_{\max}\theta)$ and $\sin(2\pi k_{\max}\theta)$. That is, steady-state opinion patterns along the bifurcation branches have the form $z(\theta,t)=A\cos(2\pi k_{\max}\theta)+B\sin(2\pi k_{\max}\theta)$ for $A,B\in{\mathbb{R}}$. In particular, the number of maxima exhibited in the opinion patterns forming at bifurcation is fixed by $k_{\max}$.

• From Lemma 12, $\lambda_{k}\in\text{Sp}(A_{G})$ are given by (12). We solve for $\lambda_{k}=\frac{1}{\tau}(-1+\alpha\hat{W}(k))=0$. Then, the first eigenvalue crossing occurs at $\alpha^{*}=\frac{1}{\hat{W}(k_{\max})}=\frac{1}{\hat{W}(-k_{\max})}$ with two eigenvalues crossing at $0$. For $\alpha<\alpha^{*}$ we have $\lambda_{k}<0$ $\forall k$ so the origin is stable. However, once $\alpha>\alpha^{*}$, there exist at least two positive eigenvalues so the origin will be unstable. For the linearization, at the bifurcation, $\lambda_{k}<0$, $\forall k\neq\pm k_{\max}$, so the corresponding spatial Fourier modes belong to the stable manifold. The spatial Fourier modes $\eta_{k_{\max}}(\theta)$ and $\eta_{-k_{\max}}(\theta)$ form a basis for the center manifold; hence, they determine the dominating bifurcation direction and emerging pattern of the model linearization. ∎

Figure 1: Bifurcation diagrams illustrating the effect of the shift value $\xi$ on the dynamics (7) with shifted sigmoid (13). Stable (unstable) branches of equilibria are shown as solid (dashed) lines.

The bifurcation of system (7) with zero input is studied in [15], [16]. There a perturbation analysis is used to show that the spatial pattern that appears is determined by the leading modes. There are infinitely many branches of non-zero equilibria which exhibit the same pattern with $k_{\max}$ maxima up to spatial translation. As discussed in [15][Theorem 4.1], [16][Section 4.2.1], the stability of the bifurcating branches can be computed as functions of $\hat{W}(k_{\max})$, $S^{\prime\prime}(0)^{2}$ and $S^{\prime\prime\prime}(0)$. Generally, when $S^{\prime\prime}(0)=0$ all of the non-zero branches of equilibria that bifurcate from the origin are stable. When $S^{\prime\prime}(0)\neq 0$, all non-zero branches bifurcating at the origin are unstable; however, due to higher order terms, stable branches of non-zero equilibria exist farther away from the origin. Fig. 1 illustrates $||z||$ at the equilibria as a function of the bifurcation parameter $\alpha$ with

$$\textstyle S(z)=\frac{\tanh(z-\xi)-\tanh(-\xi)}{\operatorname{sech}^{2}(\xi)},$$ (13)

a shifted hyperbolic tangent with $\xi\geq 0$ shift. Note that $S^{\prime\prime}(0)=0$ for $\xi=0$ and $S^{\prime\prime}(0)\neq 0$ $\forall\xi\neq 0$.

We see that for $\xi\neq 0$, there are regions below the bifurcation point for which the neutral and a non-neutral solution are both stable. This bistable region enables rapid formation of strong opinions in response to spatially distributed input, as discussed in Section VI. The patterns of opinion formation depend on the kernel, which can be designed. Fig. 2 shows how $k_{\max}$, the spatial frequency corresponding to the largest Fourier coefficient of the kernel, determines the number of maxima exhibited in the steady-state opinion pattern of the agent for (7) with zero-input. Spatial invariance ensures that for all initial conditions the solution converges to the same opinion pattern modulo a spatial translation (Fig. 2c).

VI DISTRIBUTED DECISION-MAKING ON THE CIRCLE WITH TUNABLE SENSITIVITY

We reintroduce distributed inputs to the model, and use its linearization, together with spatial and temporal frequency analysis, to infer the nonlinear input-output behavior.

We make the following assumptions.

Assumption 2 (Shifted sigmoid):

We assume $S^{\prime\prime}(0)\neq 0$ to ensure that a bistable region exists (see Fig. 1 for $\xi\neq 0$).

Assumption 3 (Input assumptions):

Inputs $u(\theta,t)\in L^{2}(\mathbb{S}^{1})$ for all $t\geq 0$. Furthermore, for all $\theta\in\mathbb{S}^{1}$, $u(\theta,\cdot):{\mathbb{R}}_{\geq 0}\to{\mathbb{R}}$ is slowly varying, that is, it is Lipschitz continuous with Lispschitz constant $0<l\ll\tau^{-1}$, for $\tau$ in (7).

The condition $l\ll\tau^{-1}$ implies that inputs vary much more slowly than the characteristic time constant of model (7). Hence, under Assumption 3, we can use the quasi-static input approximation and let $u(\theta,t)\equiv u_{h}(\theta)$.

For any diagonalized spatially distributed system of the form (5), if $[\mathcal{L}\hat{u}](k,s)$, $[\mathcal{L}\hat{z}](k,s)$ exist, then the transfer function $H(k,s)$ characterizes the input-output response in terms of the Laplace transforms of $\hat{z}(k,t)$ and $\hat{u}(k,t)$, i.e.,

$$[\mathcal{L}\hat{z}](k,s)=H(k,s)[\mathcal{L}\hat{u}](k,s).$$

By the Final Value Theorem, if the input is constant in time and  (5) is stable, then $\lim_{t\to\infty}\hat{z}(k,t)\!=\!sH(k,0)\Big{(}\frac{\hat{u}_{h}(k)}{s}\Big{)}\!=\!H(k,0)\hat{u}_{h}(k)$. This leads us to the following definition.

Figure 2: Influence of the kernel design on the steady-state opinion patterns of (7) with zero-input. (a) Two kernel designs. (b) Fourier coefficients of the kernel. Top: $\pm k_{\max}=\pm 1$. Bottom: $\pm k_{\max}=\pm 3$. (c) Steady-state opinion pattern $z(\theta,\infty)$, of dynamics (7) for initial conditions $z(\theta,0)$. The number of maxima of $z(\theta,\infty)$ equals $k_{\max}$ of the corresponding kernel. Parameters: $d=1$, $\alpha=0.98$, $p=3$, $\xi=0.7$.

Definition 8 (Spatial transfer function):

The spatial transfer function of (5) is

$$\textstyle\tilde{H}(k)=H(k,0).$$ (14)

Spatial transfer function $\tilde{H}(k)$ determines the steady-state output of (5) in response to input that is constant in time.

Theorem 4 (Spatial transfer function of (10)):

Let Assumptions 1–3 hold. Let $\pm k_{\max}=\arg\max_{k\in{\mathbb{Z}}}\hat{W}(k)$ denote the spatial frequencies corresponding to the largest $\hat{W}(k)$. Spatial transfer function (14) of the linearized model (10) is

$$\textstyle\tilde{H}(k)=\frac{\tau}{1-{\alpha}\hat{W}(k)},\quad k\in{\mathbb{Z}}.$$ (15)

In particular, for ${\alpha}\nearrow\alpha^{*}$, $\tilde{H}(\pm k_{\max})\to\infty$.

• From Lemma 12 we know the eigenvalues of (11), the diagonalization of (10). So we can compute ${H}(k,s)=\frac{1}{s-\lambda_{k}}$. Then $\tilde{H}(k)=\frac{\tau}{1-{\alpha}\hat{W}(k)}.$ As ${\alpha}\nearrow\alpha^{*}$, $H(\pm k_{\max},0)\to\infty$. ∎

Theorem 15 implies that close to bifurcation, i.e., for ${\alpha}\nearrow\alpha^{*}$, the input-output response of the linearized system is dominated by the leading spatial Fourier modes $\eta_{k_{\pm\max}}(\theta)$. This means that the alignment of $u(\theta)$ with $\eta_{\pm k_{\max}}(\theta)$ is the main determinant of the model response to inputs.

For the nonlinear system with ${\alpha}\nearrow\alpha^{*}$, the input-output response takes place in the bistable region. If $\langle\eta_{\pm k_{\max}}(\theta),{u_{h}(\theta)}\rangle=\hat{u}_{h}(k_{\max})\neq 0$, the input is aligned with the leading modes $\eta_{\pm k_{\max}}(\theta)$. These modes filter input nonlinearity and amplify the input because, by Theorem 15, the direction of these modes are ultrasensitive to input. The result is a steady-state opinion pattern with $k_{\max}$ maxima and $||z(\theta)||\gg||u_{h}(\theta)||$, as illustrated in the top row of Fig. 3. As shown in Fig. 3a (top), the Fourier coefficients of the input for $\pm k_{\max}=\pm 1$ are nonzero meaning the input is aligned with $\eta_{\pm k_{\max}}(\theta)=\eta_{\pm 1}(\theta)$. The input distribution in Fig. 3b (top) is small in magnitude (less than 0.01), while the resulting steady-state opinion pattern in Fig. 3c (top) has $k_{\max}=1$ maximum that is greater than 1.5 in magnitude.

If $\langle\eta_{\pm k_{\max}(\theta)},{u}_{h}(\theta)\rangle=0$, the input is unaligned with the leading modes $\eta_{\pm k_{\max}}(\theta)$ and the steady-state opinion pattern will not exhibit large maxima. I.e., because the input does not have a component along the ultrasensitive direction, by Theorem 15 it does not get amplified, as illustrated in the bottom row of Fig. 3. As shown in Fig. 3a (bottom), the Fourier coefficients of the input for $\pm k_{\max}=\pm 1$ are zero meaning the input is unaligned with $\eta_{\pm k_{\max}}(\theta)=\eta_{\pm 1}(\theta)$. The input distribution in Fig. 3b (bottom) is small in magnitude (less than 0.01) and the resulting steady-state opinion pattern in Fig. 3c (bottom) has no large maximum.

Our results show how even very small distributed inputs can trigger the formation of a strong opinion and whether or not this happens depends on the design of kernel $W$. Thus, $W$ can be designed to tune response to be ultrasensitive to inputs that matter for function and robust to inputs that don’t.

Figure 3: Input-output behavior of the dynamics (7) with input distributions aligned or unaligned with the Fourier mode corresponding to $\pm k_{\max}\!=\!\pm 1$. Top row: Aligned. Bottom row: Unaligned. (a) Magnitude of the Fourier coefficients of the input. (b) Input distribution. (c) Steady-state opinion pattern $z(0,\infty)$. Parameters: $d\!=\!1$, $\alpha\!=\!0.98$, $p\!=\!3$, $\xi\!=\!0.7$.

VII APPLICATION TO ROBOT NAVIGATION

We illustrate the benefits of our approach to perceptual decision-making with an application of the dynamics (7) to robot spatial navigation. We consider the case of a robot moving in a crowded space, such as an airport, where it must pass through gaps of different sizes (e.g., between people in a line) that may change over time. We assume the robot has a (visual) sensor so that it can perceive these gaps.

We specialize to a scenario where a robot finds itself trapped inside a circle of people and needs to choose and cross through a large enough gap between people. Choosing a gap is challenging as people may be distributed unevenly around the circle, resulting in multiple gap options, only a select few of which may be suitable for the robot to cross. Also, the size of the gaps may change over time due to people moving for their own purposes or in response to the robot, e.g., people may move to make space for the robot to cross.

In Section VII-A we present a framework for designing $W$ from its Fourier coefficients to allow the robot to select a single gap. We discuss four scenarios that demonstrate how our model can be used for fast-and-flexible decision-making in this robotic navigation problem. In Section VII-B two scenarios demonstrate the robot’s ability to choose a single gap, while in Section VII-C the two other scenarios show how the robot can quickly adapt to changes in gap sizes.

We take $\mathbb{S}^{1}$ to represent the circular visual field for the robot. Then an option $\theta\in\mathbb{S}^{1}$ represents the angle associated to a point in the visual field. The input $u(\theta,t)$ is the visual observation (e.g., pixel) at angle $\theta$ at time $t$. We let a point in the input distribution that reflects a gap be represented by $u(\theta,t)>0$, in blue in Fig. 4 (bottom). We assume changes in gaps occur slowly enough that Assumption 3 holds. The opinion $z(\theta,t)$ as shown in Fig. 4 (top), captures the robot’s preference over time for one gap, where the preference corresponds to the strongest opinion (in yellow).

VII-A Fourier-Based Kernel Design

We leverage the results of Theorem 15 to design kernel $W$ on $\mathbb{S}^{1}$ that imposes the desired opinion formation behavior in response to distributed input on $\mathbb{S}^{1}$. Options (angles) that are close (far) to each other should have an excitatory (inhibitory) interaction. And the opinion pattern should have a single maximum, so that the robot selects a single gap. From the results summarized in Section VI, to achieve this we design $W$ as follows and ensure that $\pm k_{\max}=\pm 1$:

$$\hat{W}(k)=\begin{cases}0\qquad&k=0,\\ \hat{f}(k)\qquad&\forall k\neq 0\end{cases}$$ (16)

where $\hat{f}(k):{\mathbb{Z}}\to{\mathbb{R}}$ is strictly decreasing, square-summable and symmetric. The strictly decreasing property ensures that $k_{\max}=\pm 1$ while the square-summability property ensures that the inverse Fourier transform of (16) exists and that $W(\theta-\phi)\in L_{2}(\mathbb{S}^{1})$, by Parseval’s identity [21]. Symmetry is required to satisfy Assumption 1. We take $\hat{W}(0)=0$ so that the homogeneous state is never a solution. Examples of functions that meet this criteria include $\hat{f_{k}}=1/k^{p}$ with $p>0.5$ and a Gaussian function of the form $\hat{f}(k)=e^{-(k-1)^{2}/p^{2}}$ where $p$ adjusts the width of the Gaussian. For the following simulations, we take $\hat{f}(k)=e^{-(k-1)^{2}/p^{2}}$.

VII-B Choosing the Best Gap and Avoiding Deadlock

We illustrate the model’s ability to pick the best among multiple gaps and to rapidly avoid deadlock when faced with two equally suitable gaps. We assume the people in the circle are not moving. In Fig. 4a (bottom) there are several gaps but one that is clearly wider than the others. In Fig. 4a (top), the robot forms a strong preference for the widest gap. In Fig. 4b (bottom), there are two equally wide gaps. In Fig. 4b (top), the robot forms strong opinions for one of the two widest gaps and avoids deadlock.

Figure 4: Decision-making of a robot selecting a gap through which to cross through a circle of non-moving people. Bottom row: gap distribution over time where gaps are indicated by $u(\theta,t)>0$ in blue. Top row: opinion pattern over time (strongest opinion in yellow). (a) One widest gap. (b) Two wide gaps of same size. Parameters: $d=1$, $u=0.96$, $\xi=0.35$, $p=3$.

VII-C Robustness and Responsiveness to Change

We illustrate the model’s robustness to unimportant change and responsiveness to important change in input. We assume the people in the circle are moving. In Fig. 5 (bottom), there is initially one very wide gap and one narrow gap. However, over time, the wide gap becomes narrower, and the narrower gap becomes wider. In Fig. 5a, the decrease in size of the initially wide gap is small enough that the robot can still fit through it and thus it does not change its choice. In Fig. 5b, the decrease in size of the initially wide gap is large enough that the robot changes its choice to the other emerging gap.

VIII DISCUSSION AND FINAL REMARKS

We presented a nonlinear opinion dynamics model for an agent making decisions about a continuous distribution of options in response to distributed input on the circle. We proved spatial invariance of the model linearization and a bifurcation of the model with zero input, which yields fast and flexible decision-making. We studied the input-output behavior of the model and design of the kernel. We demonstrated key advantages of the model in robot perceptual decision-making problem. In future work we aim to generalize the analysis to a continuous distribution of agents and extend the analysis to manifolds other than $\mathbb{S}^{1}$, e.g. ${\mathbb{R}}^{N}$. We also aim to implement this model for perceptual decision-making in robotics where the dynamics are in a closed loop with the physical dynamics of the agent.