Heteroclinic channels
This article is an orphan, as no other articles link to it. Please introduce links to this page from related articles; try the Find link tool for suggestions. (August 2024) |
Heteroclinic channels are ensembles of trajectories that can connect saddle equilibrium points in phase space.[1] Dynamical systems and their associated phase spaces can be used to describe natural phenomena in mathematical terms; heteroclinic channels, and the cycles (or orbits) that they produce, are features in phase space that can be designed to occupy specific locations in that space. Heteroclinic channels move trajectories from one equilibrium point to another. More formally, a heteroclinic channel is a region in phase space in which nearby trajectories are drawn closer and closer to one unique limiting trajectory, the heteroclinic orbit. Equilibria connected by heteroclinic trajectories form heteroclinic cycles and cycles can be connected to form heteroclinic networks.[2] Heteroclinic cycles and networks naturally appear in a number of applications, such as fluid dynamics,[3][4] population dynamics,[5] and neural dynamics.[6][7] In addition, dynamical systems are often used as methods for robotic control. In particular, for robotic control, the equilibrium points can correspond to robotic states, and the heteroclinic channels can provide smooth methods for switching from state to state.[8][9]
![Four saddle equilibrium points (gray) connected by heteroclinic channels (red).](http://upload.wikimedia.org/wikipedia/commons/thumb/3/34/12.LinkingSaddlesIntoACycle.jpg/220px-12.LinkingSaddlesIntoACycle.jpg)
Overview
editHeteroclinic channels (or heteroclinic orbits) are building blocks for a subset of dynamical systems that are built around connected saddle equilibrium points. Homoclinic channels/orbits join a single equilibrium point to itself, whereas heteroclinic channels join two different saddle equilibrium points in phase space. The connection is formed from the unstable manifold of the first saddle (“pushing away” from that point) to the stable manifold of the next saddle point (“pulling towards” this point). Combining at least three saddle equilibria in this way produces a heteroclinic cycle,[1] and multiple heteroclinic cycles can be connected into heteroclinic networks.
Heteroclinic channels have both spatial and temporal features in phase space. Spatial because they affect trajectories within a certain region around themselves,[1] and temporal because the parameters of a heteroclinic channel affect how much time a trajectory spends along that channel (or more specifically, how much time it spends around one of the saddle points).[9] The transient nature of heteroclinic channels is important for describing their “switching” nature. That is, some neighborhood around each equilibrium point can be defined as a separate state, and the heteroclinic channel itself presents a method of switching sequentially between these states.[8][10]
Heteroclinic "switching" is an important descriptor for natural phenomena, especially in neural dynamics. It has also been used as an approach for designing robotic control methods which cycle between states, whether those states are pre-defined behaviors[11] or transient states that lead to larger behaviors.[12]
History
editThe mathematical image described above – a series of states with a functional mechanism for switching between them – also describes a phenomenon known as winnerless competition (WLC). Winnerless competition describes the switching phenomenon between two competitive states and was identified by Busse & Heikes in 1980 when they were investigating the change of phases in a convection cycle.[3] However, the transient dynamics of WLC are widely agreed to first have been presented by Alfred J. Lotka, who first developed the concept to describe autocatalytic chemical reactions in 1910[13] and then developed an extended version in 1925 to describe ecological predator-prey relationships. In 1926, Vito Volterra independently published the same set of equations with a focus on mathematical biology, especially multi-species interactions.[14] These equations, now known as the Lotka-Volterra equations, are widely used as a mathematical model to describe transient heteroclinic switching dynamics.
Heteroclinic cycles which describe the transition between at least three states were first described by May & Leonard in 1975. They identified a special case of the Lotka-Volterra equations for population dynamics [5]. The re-emergence of heteroclinic cycles and the increased ability to do numerical computations as compared to the period of Lotka and Volterra, prompted a resurgence of interest in heteroclinic channels, cycles, and networks as mathematical models for transient sequential dynamics.
Heteroclinic channels have become models for neural dynamics. An example is Laurent et al. (2001) who described the neural responses of fish and insects to olfactory stimuli as a WLC system, where each stimulus and its response could be identified as a separate state within the space.[15] The responses could be modeled in this way because of their spatial and temporal properties, which aligned with the spatiotemporal nature of WLC. Rabinovich et al. (2001) & Afraimovich et al. (2004) used WLC networks (via the Fitzhugh-Nagumo & Lotka-Volterra models, respectively) to connect the mathematical concept of stable heteroclinic channels (SHCs) to transient neural dynamics more generally,[6][16] particularly other sensory processes and more abstract neural connections. Rabinovich et al. (2008) expanded this idea to larger cognitive dynamic systems, and large-scale brain networks.[9][17][18] Stable heteroclinic channels have also been used to model neuromechanical systems. The feeding structures and associated feeding processes (stages of swallowing) of marine mollusks have been analyzed using heteroclinic channels.[19][20]
Biological models have always been a source of inspiration for roboticists, especially those interested in robotic control. Since robotic control requires defining and sequencing the physical actions of the robot, models of neural dynamics can be very useful. An example of this can be found in central pattern generators, which are widely used for rhythmic robotic motion.[21] Heteroclinic channels have been used to replicate central pattern generators for robot control.[12] Similarly, dynamic movement primitives, another common robotic motion control system, have been adapted and made more flexible by using heteroclinic channels.[22] In more practical applications, stable heteroclinic channels have been directly used in the control of several biologically-inspired robots [11][23][24]
Concepts
editDynamical systems
editA dynamical system is a rule or set of rules that describe the evolution of a state (or a system of states) in time. The set of all possible states is called the state space. The phase space is the state space of a continuous system. Dynamical systems describe the state over time with mathematical equations, often ordinary differential equations. The current state at a particular time can be plotted as a point in phase space. The set of points over time can be plotted as a trajectory.
Stability
editA heteroclinic channel itself can be asymptotically stable. That is, any point near the vicinity of the channel is attracted to the heteroclinic cycle at the core of the channel. Both heteroclinic channels and cycles can be robust (or structurally stable) if, within a given parameter range, they maintain a given behavior; however, this is not required.
Stochasticity
editNoise is one input into a heteroclinic system to move it from one equilibrium to the next. The reason is that noise (or some other stochasticity) disturbs the system enough to move it into the vicinity of the next saddle equilibrium point in the sequence. The amount of noise required is inversely proportional to the “attractiveness” of the saddle points; the more attractive the stable part of the saddle is to the system state, the longer the trajectory will linger in its vicinity, and the more noise will be required to move the system’s state off of that attractive equilibrium point. There are also other ways of moving between the equilibrium points including parametric changes, or using sensory feedback.[19][25]
Control theory
editControl theory, in robotics, deals with the use of dynamical systems to control robotic systems. The goal of robotic control is to perform precise, coordinated actions using physical actuators in response to sensor input. Dynamical systems can be used to drive the robot to a desired state (or set of states) using sensor input to minimize actuator errors.
Mathematical definition
editAn equilibrium point in a dynamical system is a solution to the system of differential equations describing a trajectory that does not change with time. Equilibrium points can be described by their stability, which are often determined by the eigenvalues of the system’s Jacobian matrix. In general, the eigenvalues of a saddle point have non-zero real parts, at least one of the real parts is positive and at least one of the real parts is negative. Any eigenvalue with a negative real value indicates a stable manifold of the saddle which attracts trajectories, whereas any eigenvalue with a positive real value indicates the unstable manifold of the saddle which repels trajectories.
Phase space definition (from Heteroclinic orbit)
editLet be the ordinary differential equation describing a continuous dynamical system. If there are equilibria at and , then a solution is a heteroclinic connection from to if
as
and
as
This implies that the connection is contained in the stable manifold of and the unstable manifold of .
Neural dynamics examples
editNeural dynamics are the non-linear dynamics that describe neural processes, from single neurons to cognitive processes and large-scale neural systems.
Lotka-Volterra model
editThis model was first presented independently by Alfred J. Lotka for autocatalytic chemical reactions[13] and then again for biological species in competition by Vito Volterra from a mathematical biology perspective.[14] Originally, this model was only considered for two species: the two chemical species in the reaction, or a predator-prey situation in a shared environment.
The original equations were based on the logistic population equation, which is popularly used in ecology.
where is the size or concentration of a species at a given time, is the growth rate and is the carrying capacity of that species.
Lotka incorporated a term for the interaction between species and, with some generalization, the series of equations can be written as follows:
In this definition, is the size or concentration of the -th species and is the total number of species. The interaction between each species is described by the matrix .
May-Leonard expansion
editMay and Leonard expanded the Lotka-Volterra equations by investigating the system in which three species interact with each other (i.e., ). They found that for a system in which each equilibrium point is a saddle with an dimensional stable manifold, and the unstable manifold connects the points sequentially, the equation above can be re-written as follows:
Explicitly for , this becomes
where the coupling matrix, , is given by
In this model, the stability of the saddle equilibria can be easily determined. The stability requirements for the formation of a stable heteroclinic cycle are with either or .[5]
It was noted in this work that the system never asymptotically reaches any of the equilibrium points, but the amount of time the trajectory spends at each equilibrium point increases with time. In ecological terms, this suggests that a single population would eventually “beat out” the other two. May & Leonard noted that this is not a practical result in biology (and also see).[25]
Winnerless Competition framework
editThe “Winnerless Competition” framework (suggested by Laurent et al.)[15] allowed a single neuron and/or a collection of synchronized neurons to be encoded between "on" and "off". Laurent et al. investigated olfaction in fish and insects, particularly olfactory reception, and some of the postsynaptic structures in the odor sensory system. They found that the processing (or encoding) of perceived odors occurred over at least three timescales: fast, intermediate, and slow. They posited that an odor encoding system should be reproducible, which requires it to be insensitive to (or rapidly forget) any initial state. This is only possible if the dynamical system is strongly dissipative, that is, it settles on a state quickly and is insensitive to internal noise. Conversely, a useful odor encoding system should be sensitive to small variations in input, which requires the system to be active. An active system uses external sources to allow small variations in initial states to grow with time. The winnerless competition framework allowed a single neuron (or node) to encode a stimulus (the “fast” timescale), or many stimuli could be encoded via stimulus-specific trajectories (the “slow” timescale).
The winnerless competition system was described by
where and characterize the activities of stimulus-specific groups and , respectively, is the number of neurons being simulated, characterizes the strength of inhibition by and (i.e., their interactions with each other), and is the current input by a stimulus to .[15]
Winnerless competition required that the inhibitory connections in the matrix were asymmetrical and cyclic. For example, for , if then , and .
Overall, this description produces a heteroclinic channel composed of several heteroclinic orbits (trajectories).
Transient neural dynamics
editSensory encoding via heteroclinic orbits (which are facilitated by heteroclinic channels) as described by Laurent et al. was extrapolated beyond the olfactory system. Rabinovich et al.[6] explored winnerless competition as a spatiotemporal dynamical system corresponding to the activity of specific neurons or groups of neurons. They identified the added stimulus as the factor that would drive a trajectory from one node along the channel to the next. Without it, the system would reduce to a steady state in which one neuron (or neuronal group) was active whereas the others were quiescent.
Afraimovich et al.[16] also developed winnerless competition using connected saddle points in phase space as a model for transient, sequential neural activity. They outlined how the saddle points should be defined, the conditions for heteroclinic connections between them and the conditions for heteroclinic sequence stability. They performed numerical simulations of the dynamics of a network with N = 50 neurons and used Gaussian noise as the external input. They found that the movement of a trajectory along each connection was initiated by the noise, and the speed of switching from one saddle to the next depended on the noise level.
Cognitive dynamics
editThe sequential switching property of stable heteroclinic channels has been expanded to describe higher-level transient cognitive dynamics, particularly sequential decision making. Rabinovich et al.[17] first introduced this idea by applying the sequential switching that characterizes stable heteroclinic channels to the sequential decision making process seen in a fixed time game. The player takes sequential actions in a changing environment to maximize some reward. For a fixed time game, in order to maximize the reward, the player must encounter as many decision states as possible. This means that within a fixed amount of time, the trajectory must pass in the vicinity of as many saddle points, or nodes, as possible. When the trajectory reached the vicinity of a saddle point, a decision-making function was applied.
The reward was maximized by choosing appropriate system parameters. One of these was a decision-making rule that corresponded to the fastest motion away from the saddle, which was the shortest time to reach the next saddle. Additionally, there was an optimal level of additive noise; the noise was high enough that the trajectory could move away from each saddle quickly, but not so high that the trajectory would be directed off the cycle entirely.
A major point of this work was that, without significant external stimulus, the player was likely to find one of two extremes: ending decision-making quickly or reaching a cycle that runs through the entire allotted time. Behaviorally, this cycle translates to habit formation (on a cognitive level) and is sensitive to external stimuli that can change the trajectory’s direction at any time.
Rabinovich & Varona[18] described sequential memory in a similar way. They also introduced “chunking”, which describes how the brain groups sequential information items into chunks at different hierarchical levels. They used stable heteroclinic channels as a framework for building these chunks into high level heteroclinic networks.
Neuromechanical models
editHeteroclinic channels have also been used as a model for neuromechanical systems in animals, particularly the feeding structures in marine mollusks.[19][20] Shaw et al. (2015) investigated potential models for the feeding behavior of Aplysia californica. They found that heteroclinic channels could more accurately match features of actual experimental data than other models such as limit cycles. Lyttle et al. (2017) showed that both the heteroclinic model and the limit cycle model of the Aplysia californica’s feeding system grant different advantages and disadvantages, such as robustness to perturbations and flexibility to inputs. They also showed that a reasonable model of the animal’s behavior could be made by switching between these modes, heteroclinic and limit cycle, using external sensory input, providing a dynamical basis for understanding both robustness and flexibility in motor systems.
Robotic control examples
editMathematical expansions of the framework are required for robotic control applications.
For higher dimensional systems, the connection/inhibition matrix can be generalized as:
or formulations similar to this.[26]
Appropriate saddle values must be assigned to make the system dissipative. The strength of the saddle can be characterized by its two largest eigenvalues: the single unstable eigenvalue, , and the weakest stable eigenvalue, . The saddle value of the -th node can be defined as
If , the -th node is dissipative and stable, and if the entire cycle will be stable.[27]
SHCs for biologically inspired robotic control
editSHCs have been used directly to control robots, particularly biologically-inspired robotic systems. SHCs have also been used to adapt existing robot control frameworks. In both instances, the special properties of SHCs were used to improve the associated control tasks. Some examples include integrated contact sensing to modulate the additive SHC noise,[11] a combined Gaussian Mixture Model to inform SHC "switching",[24] a central pattern generator which was adapted to be temporally sensitive,[12] and a modified control framework which has an intuitive visualization property.[22]
Contact sensing for additive noise modulation
editSHCs make it possible to use sensory feedback for rapid choices in a high degree-of-freedom robotic system. For example, Daltorio et al. used SHCs as a controller for the simulated locomotion of a worm-like robot in a pipe.[11] The robot's structure consisted of 12 actuated body segments, each with one degree of freedom: segment length. Each segment coupled its height to the length such that as the length decreased, the height increased. This structure was used to simulate peristaltic locomotion, as the segments’ actuation was coordinated to form a peristaltic wave down the robot, with each segment contracting one after the other down the robot body.
For this system, each body segment was associated with a saddle point in the SHC system. The multi-dimensional connection matrix was constructed so that each point inhibited its neighbors except the point immediately after it. This asymmetry caused the active SHC node to move “backwards” down the robot structure, while the body moved forward.
The controller was tested in multiple pipe-shaped paths where contact sensors on the robot could provide information on the environment. Contact sensing information was used to modulate the noise added to the system, which in turn allowed the activation sequence to be altered. This was key for highly coordinated movement across all segments.
Gaussian mixture model to inform "switching"
editSHCs can be used to inform the switching among complex configurations. Petrič et al. used a combined Gaussian Mixture Model (GMM) and SHC system to control a spinal exoskeleton.[24] The exoskeleton was designed as a quasi-passive system that physically supports the user to different degrees depending on the current pose or movement of the user. Different functional poses/movements were identified as the nodes within the SHC system. GMMs were used to indicate what the additive inputs for each SHC node should be, which would drive the system from one pose to the next.
Temporally sensitive central pattern generator
editSHCs have been used as an alternative to central pattern generators for robotic control.[12] Horchler et al. used SHCs to produce an oscillator whose behavior near each node could be manipulated using system parameters: additive noise and saddle values. This produced a cyclic controller that could spend more time at a particular node when needed. The controller's responsiveness to external input was demonstrated by pausing and resetting the cycle using additive noise.
Intuitive visualization property
editRouse & Daltorio replaced the underlying attractor points of dynamic movement primitives, another biologically-inspired robotic control method, with the saddle points of SHCs.[22] This adaptive framework maintained the stability of the system. Additionally, it provided a visualization property which allowed the user to intuitively place saddle points in phase space to match a desired trajectory in the task space.
See also
editReferences
edit- ^ a b c Melbourne, Ian (1989-10-01). "Intermittency as a codimension-three phenomenon". Journal of Dynamics and Differential Equations. 1 (4): 347–367. Bibcode:1989JDDE....1..347M. doi:10.1007/BF01048454. ISSN 1572-9222. S2CID 121088204.
- ^ Kirk, V; Silber, M (1994-11-01). "A competition between heteroclinic cycles". Nonlinearity. 7 (6): 1605–1621. Bibcode:1994Nonli...7.1605K. doi:10.1088/0951-7715/7/6/005. ISSN 0951-7715. S2CID 250772340.
- ^ a b Busse, F. H.; Heikes, K. E. (1980-04-11). "Convection in a Rotating Layer: A Simple Case of Turbulence". Science. 208 (4440): 173–175. Bibcode:1980Sci...208..173B. doi:10.1126/science.208.4440.173. ISSN 0036-8075. PMID 17745536. S2CID 12241661.
- ^ Halcrow, J.; Gibson, J. F.; Cvitanović, P.; Viswanath, D. (February 2009). "Heteroclinic connections in plane Couette flow". Journal of Fluid Mechanics. 621: 365–376. arXiv:0808.1865. Bibcode:2009JFM...621..365H. doi:10.1017/S0022112008005065. ISSN 1469-7645. S2CID 884871.
- ^ a b c May, Robert M.; Leonard, Warren J. (1975-09-01). "Nonlinear Aspects of Competition Between Three Species". SIAM Journal on Applied Mathematics. 29 (2): 243–253. doi:10.1137/0129022. ISSN 0036-1399.
- ^ a b c Rabinovich, M.; Volkovskii, A.; Lecanda, P.; Huerta, R.; Abarbanel, H. D. I.; Laurent, G. (2001-07-20). "Dynamical Encoding by Networks of Competing Neuron Groups: Winnerless Competition". Physical Review Letters. 87 (6): 068102. Bibcode:2001PhRvL..87f8102R. doi:10.1103/PhysRevLett.87.068102. PMID 11497865.
- ^ Bick, Christian; Rabinovich, Mikhail I. (2010-03-01). "On the occurrence of stable heteroclinic channels in Lotka–Volterra models". Dynamical Systems. 25 (1): 97–110. doi:10.1080/14689360903322227. ISSN 1468-9367. S2CID 122024447.
- ^ a b Hansel, D.; Mato, G.; Meunier, C. (1993-11-01). "Clustering and slow switching in globally coupled phase oscillators". Physical Review E. 48 (5): 3470–3477. Bibcode:1993PhRvE..48.3470H. doi:10.1103/PhysRevE.48.3470. PMID 9961005.
- ^ a b c Rabinovich, Misha; Huerta, Ramon; Laurent, Gilles (2008-07-04). "Transient Dynamics for Neural Processing". Science. 321 (5885): 48–50. doi:10.1126/science.1155564. ISSN 0036-8075. PMID 18599763. S2CID 34129599.
- ^ Bick, Christian (2018-05-01). "Heteroclinic switching between chimeras". Physical Review E. 97 (5): 050201. arXiv:1703.03274. Bibcode:2018PhRvE..97e0201B. doi:10.1103/PhysRevE.97.050201. hdl:10871/32679. PMID 29906933. S2CID 49234225.
- ^ a b c d Daltorio, Kathryn A.; Horchler, Andrew D.; Shaw, Kendrick M.; Chiel, Hillel J.; Quinn, Roger D. (2013). "Stable Heteroclinic Channels for Slip Control of a Peristaltic Crawling Robot". In Lepora, Nathan F.; Mura, Anna; Krapp, Holger G.; Verschure, Paul F. M. J.; Prescott, Tony J. (eds.). Biomimetic and Biohybrid Systems. Lecture Notes in Computer Science. Vol. 8064. Berlin, Heidelberg: Springer. pp. 59–70. doi:10.1007/978-3-642-39802-5_6. ISBN 978-3-642-39802-5.
- ^ a b c d Horchler, Andrew D; Daltorio, Kathryn A; Chiel, Hillel J; Quinn, Roger D (2015-02-25). "Designing responsive pattern generators: stable heteroclinic channel cycles for modeling and control". Bioinspiration & Biomimetics. 10 (2): 026001. Bibcode:2015BiBi...10b6001H. doi:10.1088/1748-3190/10/2/026001. ISSN 1748-3190. PMID 25712192. S2CID 33000958.
- ^ a b Lotka, Alfred J. (1910-03-01). "Contribution to the Theory of Periodic Reactions". The Journal of Physical Chemistry. 14 (3): 271–274. doi:10.1021/j150111a004. ISSN 0092-7325.
- ^ a b Volterra, Vito (1926-10-01). "Fluctuations in the Abundance of a Species considered Mathematically1". Nature. 118 (2972): 558–560. Bibcode:1926Natur.118..558V. doi:10.1038/118558a0. ISSN 1476-4687. S2CID 4007591.
- ^ a b c Laurent, Gilles; Stopfer, Mark; Friedrich, Rainer W; Rabinovich, Misha I; Volkovskii, Alexander; Abarbanel, Henry DI (2001-03-01). "Odor Encoding as an Active, Dynamical Process: Experiments, Computation, and Theory". Annual Review of Neuroscience. 24 (1): 263–297. doi:10.1146/annurev.neuro.24.1.263. ISSN 0147-006X. PMID 11283312.
- ^ a b Afraimovich, V. S.; Zhigulin, V. P.; Rabinovich, M. I. (2004). "On the origin of reproducible sequential activity in neural circuits". Chaos: An Interdisciplinary Journal of Nonlinear Science. 14 (4): 1123–1129. Bibcode:2004Chaos..14.1123A. doi:10.1063/1.1819625. PMID 15568926 – via AIP.
- ^ a b Rabinovich, Mikhail I.; Huerta, Ramón; Varona, Pablo; Afraimovich, Valentin S. (2008-05-02). "Transient Cognitive Dynamics, Metastability, and Decision Making". PLOS Computational Biology. 4 (5): e1000072. Bibcode:2008PLSCB...4E0072R. doi:10.1371/journal.pcbi.1000072. ISSN 1553-7358. PMC 2358972. PMID 18452000.
- ^ a b Rabinovich, Mikhail I.; Varona, Pablo (2018). "Discrete Sequential Information Coding: Heteroclinic Cognitive Dynamics". Frontiers in Computational Neuroscience. 12: 73. doi:10.3389/fncom.2018.00073. ISSN 1662-5188. PMC 6137616. PMID 30245621.
- ^ a b c Shaw, Kendrick M.; Lyttle, David N.; Gill, Jeffrey P.; Cullins, Miranda J.; McManus, Jeffrey M.; Lu, Hui; Thomas, Peter J.; Chiel, Hillel J. (2015-02-01). "The significance of dynamical architecture for adaptive responses to mechanical loads during rhythmic behavior". Journal of Computational Neuroscience. 38 (1): 25–51. doi:10.1007/s10827-014-0519-3. ISSN 1573-6873. PMC 4544651. PMID 25182251.
- ^ a b Lyttle, David N.; Gill, Jeffrey P.; Shaw, Kendrick M.; Thomas, Peter J.; Chiel, Hillel J. (2017-02-01). "Robustness, flexibility, and sensitivity in a multifunctional motor control model". Biological Cybernetics. 111 (1): 25–47. doi:10.1007/s00422-016-0704-8. ISSN 1432-0770. PMC 5326633. PMID 28004255.
- ^ Ijspeert, Auke Jan (2008-05-01). "Central pattern generators for locomotion control in animals and robots: A review". Neural Networks. Robotics and Neuroscience. 21 (4): 642–653. doi:10.1016/j.neunet.2008.03.014. ISSN 0893-6080. PMID 18555958.
- ^ a b c Rouse, Natasha A.; Daltorio, Kathryn A. (April 2021). "Visualization of Stable Heteroclinic Channel-Based Movement Primitives". IEEE Robotics and Automation Letters. 6 (2): 2343–2348. doi:10.1109/LRA.2021.3061382. ISSN 2377-3766. S2CID 232216646.
- ^ Daltorio, Kathryn A; Boxerbaum, Alexander S; Horchler, Andrew D; Shaw, Kendrick M; Chiel, Hillel J; Quinn, Roger D (2013-08-27). "Efficient worm-like locomotion: slip and control of soft-bodied peristaltic robots". Bioinspiration & Biomimetics. 8 (3): 035003. Bibcode:2013BiBi....8c5003D. doi:10.1088/1748-3182/8/3/035003. ISSN 1748-3182. PMID 23981561. S2CID 21625970.
- ^ a b c Petrič, Tadej; Jamšek, Marko; Babič, Jan (2021). Lenarčič, Jadran; Siciliano, Bruno (eds.). "Exoskeleton Control Based on Network of Stable Heteroclinic Channels (SHC) Combined with Gaussian Mixture Models (GMM)". Advances in Robot Kinematics 2020. Springer Proceedings in Advanced Robotics. 15. Cham: Springer International Publishing: 341–348. doi:10.1007/978-3-030-50975-0_42. ISBN 978-3-030-50975-0. S2CID 224915139.
- ^ a b Barendregt, Nicholas W.; Thomas, Peter J. (2021-11-10). "Heteroclinic cycling and extinction in May-Leonard models with demographic stochasticity". arXiv:2111.05902 [q-bio.PE].
- ^ Ghosh, Souraj (2020-06-09). "Robustness of recurrent dynamics in disordered Lotka-Volterra system". EPL (Europhysics Letters). 130 (3): 38004. Bibcode:2020EL....13038004G. doi:10.1209/0295-5075/130/38004. ISSN 1286-4854. S2CID 225759661.
- ^ Afraimovich, Valentin S.; Rabinovich, Mikhail I.; Varona, Pablo (2004-04-01). "Heteroclinic contours in neural ensembles and the winnerless competition principle". International Journal of Bifurcation and Chaos. 14 (4): 1195–1208. arXiv:nlin/0304016. Bibcode:2004IJBC...14.1195A. doi:10.1142/S0218127404009806. ISSN 0218-1274. S2CID 12694109.
Further reading
edit- J. Lotka, Elements of Physical Biology. Baltimore: Williams and Wilkins, 1925.
- M. Krupa, “Robust Heteroclinic Cycles,” J. Nonlinear Sci, vol. 7, pp. 129–176, 1997.
- M. I. Rabinovich and M. K. Muezzinoglu, “Mutual Emotion-Cognition Dynamics,” Association for the Advancement of Artificial Intelligence, 07-Sep-2009. [Online]. Available: arXiv:0909.1144. [Accessed: 02-Nov-2021].
- Antonio Palacios (2007) Heteroclinic cycles. Scholarpedia, 2(1):2352.
- Eugene M. Izhikevich (2007) Equilibrium. Scholarpedia, 2(10):2014.
- James Meiss (2007) Dynamical systems. Scholarpedia, 2(2):1629.
- Frank Hoppensteadt (2006) Predator-prey model. Scholarpedia, 1(10):1563.
- Randall D. Beer (2009) Biologically inspired robotics. Scholarpedia, 4(4):1531.