Abstract
We present and discuss a new numerical algorithm for solving the two-dimensional Neural Field Equation with space-dependent delays: c∂∂tV(x¯,t)=I(x¯,t)−V(x¯,t)+∫ΩK(|x¯−y¯|)S(V(y,t−τ(x¯,y¯))dy¯t∈[0,T], x¯∈Ω⊂ℝ2 used in Neuroscience to describe the evolution of a population of neurons and the interactions between them. A similar equation (without delays) was first introduced by Wilson and Cowan [1], and then by Amari [2]. Here V(x,t) represents the post-synaptic membrane potential at instant t and position x. The function I represents external sources of excitation and S describes the dependence between the firing rate of the neurons and their membrane potential (typically it is a function of sigmoidal type). The kernel function K(|x-y|) describes the connectivity between the neurons at positions × and y. The delay τ(x,y) takes into the consideration the time spent by an electrical signal to travel between these two positions. The new numerical method presented in this work is based on an implicit second order method for discretisation in time and uses Gaussian quadratures for space integration. We use low-rank methods to reduce the computational effort, which enables to reduce very significantly the dimensions of the involved matrices, without affecting the final accuracy of the method. This algorithm is targeted directly to the application in Neuroscience and Robotics. As an illustration of such applications, we consider a neural field described in [3], where the firing rate function has the form S(x)= 2/ (1+exp(- μ x)) (μ is a certain positive constant), and the connectivity function is given by K(x)= 1/(2 πξ12)1/2 exp(- ||x||2/ 2 ξ12) - A/(2 πξ22)1/2 exp(- ||x||2/ 2 ξ22), with A, ξ1, ξ2 given positive constants. It is known [3] that the stability of the trivial solution of this equation depends on the value of μ and on the parameters of the connectivity function. In particular, for each set of values of these parameters, it is possible to compute a bifurcation value μbif, such that if μ < μbif the zero solution is stable, and otherwise it is unstable. With the help of our algorithm, these bifurcation values can be efficiently computed. Some examples are displayed in Table Table11 . The presented values are for the case τ(x,y)=0. Table 1
Highlights
We present and discuss a new numerical algorithm for solving the two-dimensional Neural Field Equation with space-dependent delays: c
The function I represents external sources of excitation and S describes the dependence between the firing rate of the neurons and their membrane potential
The new numerical method presented in this work is based on an implicit second order method for discretisation in time and uses Gaussian quadratures for space integration
Summary
We present and discuss a new numerical algorithm for solving the two-dimensional Neural Field Equation with space-dependent delays: c. From 24th Annual Computational Neuroscience Meeting: CNS*2015 Prague, Czech Republic. K(|x − y|)S(V(y, t − τ (x, y))dyt ∈ [0, T], x ∈ ⊂ R2 used in Neuroscience to describe the evolution of a population of neurons and the interactions between them. A similar equation (without delays) was first introduced by Wilson and Cowan [1], and by Amari [2].
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.