<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Robot. AI</journal-id>
<journal-title>Frontiers in Robotics and AI</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Robot. AI</abbrev-journal-title>
<issn pub-type="epub">2296-9144</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">632804</article-id>
<article-id pub-id-type="doi">10.3389/frobt.2021.632804</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Robotics and AI</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Adaptive and Energy-Efficient Optimal Control in CPGs Through Tegotae-Based Feedback</article-title>
<alt-title alt-title-type="left-running-head">Zamboni et&#x20;al.</alt-title>
<alt-title alt-title-type="right-running-head">Adaptive and Efficient Tegotae Feedback</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Zamboni</surname>
<given-names>Riccardo</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Owaki</surname>
<given-names>Dai</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/234417/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Hayashibe</surname>
<given-names>Mitsuhiro</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/84197/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<label>
<sup>1</sup>
</label>Politecnico di Milano, <addr-line>Milan</addr-line>, <country>Italy</country>
</aff>
<aff id="aff2">
<label>
<sup>2</sup>
</label>Department of Robotics, Graduate School of Engineering, Tohoku University, <addr-line>Sendai</addr-line>, <country>Japan</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/63830/overview">Claudius Gros</ext-link>, Goethe University Frankfurt, Germany</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/149990/overview">Eduardo J.&#x20;Izquierdo</ext-link>, Indiana University Bloomington, United&#x20;States</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/539496/overview">Guoyuan Li</ext-link>, NTNU &#xc5;lesund, Norway</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Dai Owaki, <email>owaki@tohoku.ac.jp</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Computational Intelligence in Robotics, a section of the journal Frontiers in Robotics and&#x20;AI</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>26</day>
<month>05</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>8</volume>
<elocation-id>632804</elocation-id>
<history>
<date date-type="received">
<day>24</day>
<month>11</month>
<year>2020</year>
</date>
<date date-type="accepted">
<day>03</day>
<month>05</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2021 Zamboni, Owaki and Hayashibe.</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Zamboni, Owaki and Hayashibe</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these&#x20;terms.</p>
</license>
</permissions>
<abstract>
<p>To obtain biologically inspired robotic control, the architecture of central pattern generators (CPGs) has been extensively adopted to generate periodic patterns for locomotor control. This is attributed to the interesting properties of nonlinear oscillators. Although sensory feedback in CPGs is not necessary for the generation of patterns, it plays a central role in guaranteeing adaptivity to environmental conditions. Nonetheless, its inclusion significantly modifies the dynamics of the CPG architecture, which often leads to bifurcations. For instance, the force feedback can be exploited to derive information regarding the state of the system. In particular, the <italic>Tegotae</italic> approach can be adopted by coupling proprioceptive information with the state of the oscillation itself in the CPG model. This paper discusses this policy with respect to other types of feedback; it provides higher adaptivity and an optimal energy efficiency for reflex-like actuation. We believe this is the first attempt to analyse the optimal energy efficiency along with the adaptivity of the Tegotae approach.</p>
</abstract>
<kwd-group>
<kwd>central pattern generator</kwd>
<kwd>sensory feedback</kwd>
<kwd>tegotae approach</kwd>
<kwd>efficiency</kwd>
<kwd>optimal control</kwd>
<kwd>learning</kwd>
<kwd>embodiment</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>The ability to efficiently move in complex environments is a key property for animals and their survival. This implies that many aspects of their morphology and central nervous system are shaped by constraints related to their locomotor skills. Animal locomotion is not generated merely from neural systems; instead, it is generated from the close interaction between neural systems, musculoskeletal systems, and the real-world environment (<xref ref-type="bibr" rid="B51">Pfeifer and Bongard, 2006</xref>; <xref ref-type="bibr" rid="B52">Pfeifer et&#x20;al., 2007</xref>). Thus, it is essential to elucidate the locomotion generation mechanism by analysing the interaction dynamics among these three systems and by analysing the neural systems themselves. Understanding these mechanisms is expected to result in contributions to biology and robotics by facilitating the design of durable and resilient robots that are energy-efficient.</p>
<p>Central pattern generators (CPGs) are neural circuits that are found in invertebrate (<xref ref-type="bibr" rid="B50">Pearson and Iles, 1973</xref>; <xref ref-type="bibr" rid="B8">B&#xe4;ssler and Wegner, 1983</xref>; <xref ref-type="bibr" rid="B9">B&#xe4;ssler, 1986</xref>) and vertebrate animals (<xref ref-type="bibr" rid="B56">Shik et&#x20;al., 1966</xref>; <xref ref-type="bibr" rid="B23">Grillner, 1975</xref>; <xref ref-type="bibr" rid="B24">Grillner, 1985</xref>). CPGs can produce rhythmic patterns of neural activity without receiving any rhythmic inputs. The term <italic>central</italic> indicates that the sensory feedback from the <italic>peripheral</italic> nervous system is not needed for generating the rhythms (<xref ref-type="bibr" rid="B41">Marder and Bucher, 2001</xref>; <xref ref-type="bibr" rid="B28">Ijspeert, 2008</xref>). Biological CPGs underlie many fundamental rhythmic activities such as chewing, breathing, and digesting. In addition, they also serve as the fundamental building blocks for locomotor neural circuits. From the perspective of control, they have several interesting characteristics such as a distributed control, the ability to deal with redundancies, the presence of fast control loops, and the ability to modulate the locomotion by using simple control signals. Owing to these properties, CPGs are considered to be transferred mathematical models. In addition, CPGs serve as the building blocks of robotic locomotion controllers and are being increasingly used in the robotics community (<xref ref-type="bibr" rid="B28">Ijspeert, 2008</xref>). To enable biologically inspired robotic control, the architecture of CPGs has been extensively adopted to generate periodic patterns for locomotion control owing to the properties of nonlinear oscillators (<xref ref-type="bibr" rid="B31">Kimura et&#x20;al., 1999</xref>; <xref ref-type="bibr" rid="B20">Fukuoka et&#x20;al., 2003</xref>; <xref ref-type="bibr" rid="B61">Tsujita et&#x20;al., 2003</xref>; <xref ref-type="bibr" rid="B3">Aoi and Tsuchiya, 2005</xref>; <xref ref-type="bibr" rid="B13">Buchli et&#x20;al., 2006</xref>; <xref ref-type="bibr" rid="B32">Kimura et&#x20;al., 2007</xref>; <xref ref-type="bibr" rid="B53">Righetti and Ijspeert, 2008</xref>; <xref ref-type="bibr" rid="B62">Wang et&#x20;al., 2011</xref>).</p>
<p>Although sensory feedback in CPGs is not necessary for generating rhythmic patterns, it plays a central role in guaranteeing adaptivity to the environmental conditions (<xref ref-type="bibr" rid="B28">Ijspeert, 2008</xref>).</p>
<p>Sensory feedback in CPGs for animal locomotion was first studied in the pioneering work on bipedal walking conducted by <xref ref-type="bibr" rid="B58">Taga et&#x20;al. (1991)</xref>, <xref ref-type="bibr" rid="B59">Taga (1994)</xref>, <xref ref-type="bibr" rid="B60">Taga (1995).</xref> In these studies, sensory information from the environment was fed back into the nervous system model to generate a walking pattern from the interaction among the nervous system model, musculoskeletal model, and environment (&#x201c;Global Entrainment&#x201d;). <xref ref-type="bibr" rid="B31">Kimura et&#x20;al. (1999)</xref>; <xref ref-type="bibr" rid="B20">Fukuoka et&#x20;al. (2003)</xref> proposed a model by integrating CPG and reflex mechanisms to realise uneven terrain quadruped walking. <xref ref-type="bibr" rid="B3">Aoi and Tsuchiya (2005)</xref>, <xref ref-type="bibr" rid="B4">Aoi and Tsuchiya (2006)</xref> focused on &#x201c;phase resetting&#x201d; (<xref ref-type="bibr" rid="B55">Schomburg et&#x20;al., 1998</xref>), a feedback mechanism found in animals, to include gait stabilisation in CPG-based control models. Aoi&#x2019;s group also applied the phase resetting feedback in CPGs to human-like musculoskeletal models of bipedal walking (<xref ref-type="bibr" rid="B5">Aoi et&#x20;al., 2010</xref>), quadrupedal gait transitions (<xref ref-type="bibr" rid="B6">Aoi et&#x20;al., 2011</xref>; <xref ref-type="bibr" rid="B7">Aoi et&#x20;al., 2013</xref>), and a hexapod walking model (<xref ref-type="bibr" rid="B2">Ambe et&#x20;al., 2018</xref>). <xref ref-type="bibr" rid="B57">Steingrube et&#x20;al. (2010)</xref>; <xref ref-type="bibr" rid="B39">Manoonpong et&#x20;al. (2010)</xref> proposed a modular neural control with bio-inspired CPG-based network and sensory feedback, demonstrating environmental adaptability, such as walking on uneven terrain and avoiding unknown obstacles, and then extended the models by introducing forward models (<xref ref-type="bibr" rid="B40">Manoonpong et&#x20;al., 2013</xref>; <xref ref-type="bibr" rid="B14">Dasgupta et&#x20;al., 2015</xref>), visual feedback (<xref ref-type="bibr" rid="B22">Goldschmidt et&#x20;al., 2014</xref>; <xref ref-type="bibr" rid="B25">Grinke et&#x20;al., 2015</xref>), muscle models (<xref ref-type="bibr" rid="B63">Xiong et&#x20;al., 2014</xref>; <xref ref-type="bibr" rid="B64">Xiong et&#x20;al., 2015</xref>), and so on. <xref ref-type="bibr" rid="B13">Buchli et&#x20;al. (2006)</xref>; <xref ref-type="bibr" rid="B44">Nachstedt et&#x20;al. (2017)</xref> proposed an adaptive frequency oscillator that could learn motion frequency adaptively and verified the generation of gait according to body characteristics. Furthermore, an interlimb coordination model that employed load information as sensory information and generated adaptive and diverse quadruped walking patterns was proposed (<xref ref-type="bibr" rid="B42">Maufroy et&#x20;al., 2010</xref>; <xref ref-type="bibr" rid="B21">Fukuoka et&#x20;al., 2015</xref>; <xref ref-type="bibr" rid="B45">Owaki and Ishiguro, 2017a</xref>). Sensory feedback inclusion significantly modifies the dynamics of the CPG&#x2019;s architecture, which often leads to bifurcations and other dynamic phenomena (<xref ref-type="bibr" rid="B6">Aoi et&#x20;al., 2011</xref>; <xref ref-type="bibr" rid="B62">Wang et&#x20;al., 2011</xref>; <xref ref-type="bibr" rid="B7">Aoi et&#x20;al., 2013</xref>).</p>
<p>To establish a systematic design principle of the sensory feedback in the CPGs to achieve biologically inspired robotic locomotion, a novel concept called &#x201c;Tegotae&#x201d; is proposed. Tegotae is a Japanese concept that describes the extent to which a perceived reaction matches the intended motor command. The potential of the Tegotae approach in reproducing animals&#x2019; locomotion and understanding the underlying mechanism has been previously demonstrated based on synthetic approaches. The Tegotae approach was first used by <xref ref-type="bibr" rid="B48">Owaki et&#x20;al. (2017)</xref> to develop a minimal model for interlimb coordination on hexapod robot locomotion with CPG-based control. <xref ref-type="bibr" rid="B29">Kano et&#x20;al. (2018)</xref> demonstrated gait transition between the concertina and scaffold-based locomotion in a snake model simulation with reflex-based control. <xref ref-type="bibr" rid="B30">Kano et&#x20;al. (2019)</xref> proposed the detailed design of the Tegotae function, particularly for motor commands, using the genetic algorithm (GA) to simulate a simple 1-D earthworm model with CPG-based control. <xref ref-type="bibr" rid="B49">Owaki et&#x20;al. (2021)</xref> demonstrated adaptive walking control on a biped model with CPG and reflex-based controllers.</p>
<p>The main contribution of this study is the construction of a specific proprioceptive feedback law through the so-called Tegotae approach (<xref ref-type="bibr" rid="B48">Owaki et&#x20;al., 2017</xref>). Together with a specific control policy, i.e. reflex-like actuation, it exploits it fruitfully based on the concept of embodied intelligence (<xref ref-type="bibr" rid="B51">Pfeifer and Bongard, 2006</xref>; <xref ref-type="bibr" rid="B52">Pfeifer et&#x20;al., 2007</xref>). Then, the feedback is applied to certain mechanical systems, i.e. hopping systems; is first considered for the simplest case of one leg, and is then extended to two legs. In such circumstances, the sensory feedback plays an important role in shaping the rhythmic patterns and ensuring coordination between the CPGs and body movements. This study demonstrates the adaptation processes as well as the acquirement of the different gait. In addition, it compares the analytical solution for the single-leg case with an optimal controller solution that is based on direct methods such as the multiple shooting methods (<xref ref-type="bibr" rid="B12">Bock and Plitt, 1984</xref>; <xref ref-type="bibr" rid="B15">Diehl et&#x20;al., 2005</xref>; <xref ref-type="bibr" rid="B19">Fagiano, 2019</xref>). This confirms the intuitions for the energy efficiency of the control policy. Finally, we extensively analyse the approach in relation with the considerations for learning and energy efficiency (<xref ref-type="bibr" rid="B26">Hayashibe and Shimoda, 2014</xref>).</p>
<p>The following section presents the materials and methods used in this study. First, we briefly describe the Tegotae approach. Second, we present the mathematical model for the Tegotae-based control. Third, we discuss the Tegotae approach based on the learning framework by comparing it with <italic>tacit learning</italic> as described in <xref ref-type="bibr" rid="B26">Hayashibe and Shimoda (2014)</xref>. Then, we present the simulation results to validate the Tegotae controller and then evaluate the energy efficiency. Finally, in <xref ref-type="sec" rid="s5">Section 5</xref>, we discuss the results and future&#x20;work.</p>
</sec>
<sec sec-type="methods" id="s2">
<title>2 Methods</title>
<sec id="s2-1">
<title>2.1 Tegotae Control</title>
<sec id="s2-1-1">
<title>2.1.1 Theory</title>
<p>The inclusion of feedback in the architecture of the CPG is a natural extension of these structures. However, any modification to the canonical form leads to a modification in the main dynamics, which may affect the effectiveness. This is achieved by considering a particular family of feedback functions in terms of the local effect of this inclusion on the dynamics of a neural oscillator. The approach to define these feedback functions is called the Tegotae approach, as described in <xref ref-type="bibr" rid="B48">Owaki et&#x20;al. (2017)</xref>. Tegotae is a novel concept that describes the extent to which a perceived reaction matches an expectation, or intention, of a controller. Tegotae stems not only from the reaction that is received from the environment, but also from the consistency between the perceived reaction and the intention or expectation of the controller, i.e. what the controller intends to do. In the case of matching, it is said that either &#x201c;good&#x201d; or &#x201c;bad&#x201d; Tegotae is obtained. In this manner, a cognitive meaning is added to the control framework, in which it denotes some actions as &#x201c;positive&#x201d; and others as &#x201c;negative&#x201d;. The objective is to maximise the Tegotae function. In this section, the Tegotae formalism is introduced. For the initial step of the investigation, Tegotae is quantified in the simplest mathematical form, i.e. a function that is based on the separation of the variables as follows.<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>u</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>e</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>Hereafter, the function <italic>T</italic> is referred to as the Tegotae function (T-function), which is a function that quantitatively measures the Tegotae. In <xref ref-type="disp-formula" rid="e1">Eq. (1)</xref>, <italic>u</italic> represents a control variable and <italic>e</italic> represents the sensory information obtained from multiple sensors that are embedded in the body. The T-function is expressed as the product between <italic>C</italic>(<italic>u</italic>) and <italic>S</italic>(<italic>e</italic>). The former expresses the intention of the controller, while the latter denotes the reaction obtained from the environment. <italic>T</italic> is designed such that it becomes more positive when an enhanced Tegotae is detected. Therefore, for a given T-function, the local sensory feedback <italic>f</italic> is designed in such a way that the control system modulates <italic>u</italic> to increase the amount of Tegotae received. Thus, with regard to the continuous-time systems, <italic>f</italic> is expressed simply as a mono-dimensional gradient system of the T-function <italic>T</italic> with respect to the control variable <italic>u</italic>, as follows.<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>With this formulation, it is possible to systematically design the decentralised controllers by only designing the T-functions that are required. When considering the CPGs&#x2019; framework, the <italic>i</italic>-th controller can be first defined as a generic Kuramoto oscillator (<xref ref-type="bibr" rid="B34">Kuramoto (1984)</xref>) of phase <italic>&#x3d5;</italic>
<sub>
<italic>i</italic>
</sub> without the coupling terms but with a specific external field <italic>f</italic>
<sub>
<italic>i</italic>
</sub> that consists of the local sensory feedback.<disp-formula id="e3">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>As a result, this equation leads to the following expression.<disp-formula id="e4">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>In <xref ref-type="bibr" rid="B48">Owaki et&#x20;al. (2017)</xref>, the T-function was expected to reproduce the hexapedal inter-limb coordination that was observed in insect locomotion by using the Kuramoto oscillators. For this reason, it was generally defined in the first case as follows.<disp-formula id="e5">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>sin</mml:mtext>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msubsup>
<mml:mi>N</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>V</mml:mi>
</mml:msubsup>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>where the sensory information <italic>e</italic> consists of the vertical ground reaction forces <italic>N</italic>
<sub>
<italic>i</italic>
</sub>
<sup>
<italic>v</italic>
</sup> that are acting on each leg. In the basic control of the hexapod robot in <xref ref-type="bibr" rid="B48">Owaki et&#x20;al. (2017)</xref>, the leg was controlled to be in the swing phase for <italic>&#x3d5;</italic>
<sub>
<italic>i</italic>
</sub> &#x3d; 0 to <italic>&#x3c0;</italic> and the stance phase for <italic>&#x3d5;</italic>
<sub>
<italic>i</italic>
</sub> &#x3d; <italic>&#x3c0;</italic> to 2<italic>&#x3c0;</italic> based on the function <italic>C</italic>&#x2009;(<italic>&#x3d5;</italic>
<sub>
<italic>i</italic>
</sub>) &#x3d; &#x2212;sin<italic>&#x3d5;</italic>
<sub>
<italic>i</italic>
</sub> In this formulation, <italic>T</italic>
<sub>
<italic>i</italic>
</sub> quantifies the Tegotae on the basis of the information that is only locally available at the corresponding leg. When the local controller intends to be in the stance phase (&#x2212;sin<italic>&#x3d5;</italic>
<sub>
<italic>i</italic>
</sub> &#x3e; 0) and receives a ground reaction force (<italic>N</italic>
<sub>
<italic>i</italic>
</sub>
<sup>
<italic>v</italic>
</sup> &#x3e; 0), <italic>T</italic>
<sub>
<italic>i</italic>
</sub> evaluates the situation as &#x201c;good&#x201d; Tegotae, and vice versa. As stated above, the reaction in <xref ref-type="disp-formula" rid="e1">Eq. 1</xref> is generic, and other types of reactions may be taken into account. In our study, the force passing through the body was taken into account, i.e. an elastic force. This definition is inspired by the Golgi tendon organ (<xref ref-type="bibr" rid="B43">Moore (1984)</xref>), which is a proprioceptive sensory receptor organ that senses changes in the muscle tension. The T-function is then defined for a generic <italic>i</italic>-th phase oscillator and the feedback signal is expressed as follows.<disp-formula id="e6">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x225c;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>F</mml:mi>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
<disp-formula id="e7">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>F</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mi>F</mml:mi>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>where <italic>&#x3c3;</italic> denotes a proportionality factor and <italic>F</italic> represents the force passing through the body. By the nature of <xref ref-type="disp-formula" rid="e6">Eq. (6)</xref>, it follows that this sensory feedback will be absent when there is no contact with the ground.</p>
</sec>
<sec id="s2-1-2">
<title>2.1.2 Tegotae Control Policy: Preliminary Design and Extensions for Reflex-like Actuation</title>
<p>In majority of the CPGs&#x2019; controllers, the actuator is driven by a proportional-integral-derivative (PID) control scheme, which compares the actual state of the physical system with the reference signal that was originated by the CPGs&#x2019; network (<xref ref-type="bibr" rid="B28">Ijspeert, 2008</xref>). One of our main contributions is to attempt to maintain the model-free control approach while taking into account some of the most recent considerations for the above embodied intelligence (<xref ref-type="bibr" rid="B51">Pfeifer and Bongard, 2006</xref>; <xref ref-type="bibr" rid="B52">Pfeifer et&#x20;al., 2007</xref>) and control by using neural-like dynamic systems and reflex-like motor control. <xref ref-type="bibr" rid="B13">Buchli et&#x20;al. (2006)</xref> demonstrates the manner in which the neuro-mechanical coupling provided by the feedback forces the secondary dynamics in the phase oscillator; our goal is to analyse and possibly exploit this effect. This study aims to use a critical point for the feedback dynamics, which is a minimum, or a specific section of it, to control the system. This section briefly describes the evolution of the Tegotae control policy towards its current form. In the former control policy law established by <xref ref-type="bibr" rid="B48">Owaki et&#x20;al. (2017)</xref>, a constant actuation force with the value <italic>A</italic> was used, and actuation was observed when the phase of the oscillator <italic>&#x3d5;</italic> was within a certain interval containing the selected critical point of the dynamics <italic>&#x3d5;</italic>
<sub>0</sub>.<disp-formula id="e8">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x394;</mml:mtext>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mtext>&#x394;</mml:mtext>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x21d2;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#x00B7;</mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>A</mml:mi>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>This implies that the force <italic>F</italic>
<sub>
<italic>a</italic>
</sub> &#x3d; <italic>A</italic> is applied when the phase <italic>&#x3d5;</italic>
<sub>
<italic>i</italic>
</sub> ranges from <italic>&#x3d5;</italic>
<sub>0</sub> &#x2212; &#x394;/2 to <italic>&#x3d5;</italic>
<sub>0</sub> &#x2b; &#x394;/2. It is apparent that a critical factor of this preliminary policy is the on-line adaptation of the values of <italic>&#x3d5;</italic>
<sub>0</sub> and &#x394; according to the evolution of the dynamics from the transient to the steady state (assuming it is reached), which is non-trivial. In the first instance, these values are considered to be a posteriori once the specific dynamic of the oscillator has been studied and maintained constantly throughout the entire simulation. The results obtained with this simple control policy are analysed in the monoped case study, which demonstrates how even this simple policy can guarantee good performance. Clearly, this policy can be made smoother by substituting the square wave with other types of functions such as bell-shaped trends.<disp-formula id="e9">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#x00B7;</mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:mtext>&#x394;</mml:mtext>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msqrt>
<mml:mtext>&#x394;</mml:mtext>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>Although this leads to an easier actuation and solves the numerical issues that are introduced by the switching controller, this control policy does not simplify the method of selection of the specific values of <italic>&#x3d5;</italic>
<sub>0</sub> and &#x394;. In contrast, the entire negative section that is centered around the minimum of the Tegotae feedback can correspond to a critical phase of the entire dynamics. The following expression can be considered.<disp-formula id="e10">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>This specifically indicates that the Tegotae is decreasing. By definition, the aim is to maximise it. It is clear how this area is the designated area to inject a certain force. In particular, this force is required to lead to the maximisation of the Tegotae, which is dependent on the case study. In this study, a positive force leading to a jump satisfies the requirements. Thus, following <xref ref-type="disp-formula" rid="e6">Eq. (6)</xref>, the final mathematical form for the reflex-like actuation that is newly proposed in this study is defined as follows.<disp-formula id="e11">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>min</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>min</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>cos</mml:mtext>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
</p>
<p>The reflex-like actuation is designed to be opposite in sign to the Tegotae feedback and disappear once the feedback is positive, indicating an increasing Tegotae (<xref ref-type="fig" rid="F1">Figure&#x20;1A</xref>). Thus, the negative sign can be attributed to fact that the force actuated in the feedback should be in a direction opposite to that of the force used as the feedback itself. This clearly reintroduces the numerical issues of the switching controller. However, it directly links the actuation and Tegotae feedback in a more biologically inspired reflex-like manner. It also assures an online adaptation to the variation of the dynamics since the Tegotae feedback corresponds to this variation itself, as shown in <xref ref-type="fig" rid="F1">Figure&#x20;1A</xref>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Tegotae approach: <bold>(A)</bold> The reflex-like actuation is designed to be opposite in sign to that of the Tegotae feedback and disappear once the feedback is positive, indicating an increasing Tegotae in <xref ref-type="disp-formula" rid="e11">Eq. 11</xref>. <bold>(B,C)</bold> Neuro-mechanical structure of the mono-dimensional hoppers. <bold>(B)</bold> Monopod: a mass is connected to a mass-less spring and a damper system. A linear actuator is parallel to the spring and damper and it determines the vertical thrust. The Kuramoto model for the phase oscillators was used as a model for the CPGs&#x2019; oscillator. <bold>(C)</bold> Biped: Two vertical hoppers are connected with a mechanical spring. Each hopper is controlled by using a decoupled Kuramoto oscillator with Tegotae feedback.</p>
</caption>
<graphic xlink:href="frobt-08-632804-g001.tif"/>
</fig>
</sec>
</sec>
<sec id="s2-2">
<title>2.2 Mechanical Model</title>
<sec id="s2-2-1">
<title>2.2.1 Monopod Model</title>
<p>First, a one-dimensional (1-D) hopping system was considered, which is characterised by a mass connected to a mass-less spring and a damper system (<xref ref-type="fig" rid="F1">Figure&#x20;1B</xref>). A linear actuator is parallel to the spring and damper and determines the vertical thrust. The Kuramoto model (<xref ref-type="bibr" rid="B34">Kuramoto, 1984</xref>) for the phase oscillators was used as a model for the CPGs&#x2019; oscillator, simplifying the analysis of the effects of the feedback. The integration of the ordinary differential equations (ODEs) was performed using MATLAB, which automatically stopped the integration when switching was detected. The initial step of the integration was set to 1<italic>e</italic>
<sup>&#x2212;3</sup>, which is equal to the maximum step of the integration. The evolution of a single phase of the oscillator <italic>&#x3d5;</italic> and the vertical height of the mass <italic>y</italic> is described by an ODE as follows.<disp-formula id="e12">
<mml:math id="m12">
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
<disp-formula id="e13">
<mml:math id="m13">
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#xa8;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>m</mml:mi>
</mml:mfrac>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>g</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
<disp-formula id="e14">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
<disp-formula id="e15">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>l</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>where <italic>f</italic>(<italic>&#x3d5;</italic>,<italic>F</italic>) is the sensory feedback in the CPG oscillator, while <italic>F</italic>
<sub>
<italic>k</italic>
</sub>(<italic>y</italic>), <inline-formula id="inf1">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, and <italic>F</italic>
<sub>
<italic>a</italic>
</sub>(<italic>&#x3d5;</italic>, &#x22C5;) represent the spring, damper, and actuator force, respectively. These three components are absent during the flight phase, assuming that there no forces that act from the environment.</p>
<p>As previously described, according to <xref ref-type="bibr" rid="B48">Owaki et&#x20;al. (2017)</xref>, the Tegotae sensory feedback <italic>f</italic>(<italic>&#x3d5;</italic>,<italic>F</italic>) is defined directly by the Tegotae function <italic>T</italic>(<italic>&#x3d5;</italic>,<italic>F</italic>), where we selected <italic>F</italic>&#x20;&#x3d; <italic>F</italic>
<sub>
<italic>k</italic>
</sub>(<italic>y</italic>).<disp-formula id="e16">
<mml:math id="m17">
<mml:mrow>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x225c;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>
<disp-formula id="e17">
<mml:math id="m18">
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(17)</label>
</disp-formula>with <italic>&#x3c3;</italic> being a proportionality factor. From <xref ref-type="disp-formula" rid="e11">Eq. (11)</xref>, <italic>F</italic>
<sub>
<italic>a</italic>
</sub> is described as follows:<disp-formula id="e18">
<mml:math id="m19">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>min</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>min</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(18)</label>
</disp-formula>
</p>
<p>Here, as a first step in the evaluation, we used the force passing through the spring <italic>F</italic>
<sub>
<italic>k</italic>
</sub>. An advantage of the Tegotae-based approach is that it can use different forces as sensory feedback. Further extensions may be a combination of many different forces. The novelty of this study lies in the reflex-like actuation equation and the validation of energetic optimality.</p>
</sec>
<sec id="s2-2-2">
<title>2.2.2 Biped Model</title>
<p>The effects of the Tegotae approach on a more complex mechanical and oscillatory system were also studied to prove its effectiveness and ability to sustain different patterns, which were also described by <xref ref-type="bibr" rid="B48">Owaki et&#x20;al. (2017)</xref>. The mechanical system was extended to a 1-D bipedal hopping robot as illustrated in <xref ref-type="fig" rid="F1">Figure&#x20;1C</xref>. The system corresponds to a slight modification of the previous case.<disp-formula id="e19">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mtext>cos</mml:mtext>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mrow>
<mml:mn>12</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mtext>sin</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(19)</label>
</disp-formula>
<disp-formula id="e20">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mtext>cos</mml:mtext>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mrow>
<mml:mn>21</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mtext>sin</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(20)</label>
</disp-formula>
<disp-formula id="e21">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#xa8;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:mi>m</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>m</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mi>g</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>12</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(21)</label>
</disp-formula>
<disp-formula id="e22">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#xa8;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:mi>m</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>m</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mi>g</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>21</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(22)</label>
</disp-formula>
</p>
<p>In Eqs (<xref ref-type="disp-formula" rid="e19">19</xref>, <xref ref-type="disp-formula" rid="e20">20</xref>), the Tegotae feedback is already taken into account, while the last term on the right-hand-side represents the weak-coupling between the phase oscillators (<xref ref-type="bibr" rid="B35">Kuramoto, 2003</xref>). In Eqs (<xref ref-type="disp-formula" rid="e21">21</xref>, <xref ref-type="disp-formula" rid="e22">22</xref>) the components are the same as those that are defined in <xref ref-type="disp-formula" rid="e12">Eq. (12)</xref>, which is from a simple additional elastic force that is introduced by the connecting spring <italic>F</italic>
<sub>
<italic>kij</italic>
</sub> &#x3d; <italic>k</italic>
<sub>
<italic>c</italic>
</sub>(<italic>y</italic>
<sub>
<italic>j</italic>
</sub>&#x2212;<italic>y</italic>
<sub>
<italic>i</italic>
</sub>). In contrast, the control policy was left unchanged with respect to the monopod case <xref ref-type="disp-formula" rid="e18">Eq.&#x20;(18)</xref>.</p>
</sec>
</sec>
</sec>
<sec id="s3">
<title>3 Tegotae in the Learning Framework</title>
<p>The Tegotae approach has certain interesting similarities with other learning frameworks, which motivates some of the intuitions for its energy efficiency. The adaptivity in the learning processes is typically defined for the parameters/weights of the controller/learning agent. In the Tegotae framework, although a further adaptation of the feedback coefficients <italic>&#x3c3;</italic> may be included, the main adaptation is induced by modifying the dynamics of the oscillators. This factor is taken into account in the comparison, since the eventual adaptation of the parameters is straightforward.</p>
<p>First, it is interesting to note how the Tegotae approach shares some similarities with the <italic>tacit learning</italic>, which is a learning framework that was introduced in <xref ref-type="bibr" rid="B10">Berenz et&#x20;al. (2014)</xref>; <xref ref-type="bibr" rid="B11">Berenz et&#x20;al. (2015).</xref> In tacit learning, the control law consists of an extension for the PD controller with a tacit learner block with the time frame (Lt). By using the scalar case for simplicity, the following expression can be obtained.<disp-formula id="e23">
<mml:math id="m24">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mi>u</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>k</mml:mi>
<mml:msubsup>
<mml:mi>x</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>T</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>q</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>q</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x222b;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>e</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtext>Lt</mml:mtext>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<label>(23)</label>
</disp-formula>where <italic>u</italic>, <italic>x</italic>
<sub>
<italic>c</italic>
</sub>, <italic>k</italic>, and <italic>e</italic> are respectively the control, the state variable that is expressed in the control space, the proportional and derivative gain, and any type of quantity that needs to be minimised. The learning process is obtained in the (Lt) block by accumulating the integral over the time of the quantity that needs to be minimised. On this basis, we neglect the proportional and derivative terms in this study.<disp-formula id="e24">
<mml:math id="m25">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mi>u</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>q</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>q</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x222b;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>e</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtext>Lt</mml:mtext>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<label>(24)</label>
</disp-formula>
</p>
<p>The function <italic>f</italic>(<italic>e</italic>) is recommended to have the form <italic>f</italic>(<italic>e</italic>) &#x3d; <italic>p</italic>(<italic>&#x3be;</italic>)<italic>a</italic>(<italic>e</italic>)<sup>
<italic>T</italic>
</sup>. In the one-dimensional case, <italic>a</italic>(<italic>e</italic>) can be a simple linear transformation <italic>a</italic>(<italic>e</italic>) &#x3d; <italic>ae</italic> and <italic>p</italic>(<italic>&#x3be;</italic>) is a periodic function of <italic>&#x3be;</italic>. Both&#x20;of these additional terms are selected to guarantee the following.<disp-formula id="e25">
<mml:math id="m26">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mfrac>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3be;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3be;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mtext>if</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mfrac>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3be;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3be;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mtext>if</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3c0;</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<label>(25)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e25">Eq. (25)</xref>, <italic>x</italic>
<sub>
<italic>e</italic>
</sub> represents the state variable that is expressed in the task space, in which the error <italic>e</italic> is minimized. In contrast, <italic>&#x3b1;</italic> is generically defined as the angle between <inline-formula id="inf2">
<mml:math id="m27">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>e</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> and <italic>D</italic>(<italic>e</italic>); the latter is the direction toward which <italic>e</italic> is minimized. In the one-dimensional case, <inline-formula id="inf3">
<mml:math id="m28">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>&#x2227;</mml:mo>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>. This formulation guarantees that min (<italic>f</italic>(<italic>e</italic>)) &#x3d; min(<italic>e</italic>). Now, let us consider the Tegotae framework. The objective is to construct feedback and not a feedforward controller. To do this, let us consider the factor that needs to be minimized that corresponds to <italic>e</italic>&#x20;&#x3d; &#x2212;<italic>F</italic>
<sub>
<italic>k</italic>
</sub>, the virtual variable <italic>&#x3be;</italic> to the physical variable <italic>&#x3d5;</italic>, and the error function <italic>a</italic>(<italic>e</italic>) &#x3d; <italic>&#x3c3;</italic>e. By neglecting the constant terms due to the integration, the feedback over the oscillator results in the following expressions.<disp-formula id="e26">
<mml:math id="m29">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mrow>
<mml:mi>u</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>q</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>sin</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x222b;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>e</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x222b;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>sin</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>&#x3d5;</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>cos</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtext>L</mml:mtext>
<mml:mi>&#x3d5;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<label>(26)</label>
</disp-formula>
</p>
<p>In the Tegotae framework, <italic>x</italic>
<sub>
<italic>e</italic>
</sub> &#x3d; &#x394;<italic>l</italic> represents the elongation speed of the spring length. This variable points towards the direction of the minimisation of the value of <italic>e</italic>&#x20;&#x3d; &#x2212;<italic>F</italic>
<sub>
<italic>k</italic>
</sub>. Thus, the following expression is obtained.<disp-formula id="e27">
<mml:math id="m30">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>sin</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mtext>sin</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mtext>if</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>sin</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mtext>sin</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mtext>if</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<label>(27)</label>
</disp-formula>
</p>
<p>This shows how the Tegotae approach is <italic>de facto</italic> obtaining a tacit learning feedback (Lt) as previously described. Nevertheless, this is achieved by accumulating the quantity that needs to be minimised for the integral of the state space variable that is directly from (L<italic>&#x3d5;</italic>). The integration over the state space frame <italic>&#x3d5;</italic> is coherent with the CPGs&#x2019; framework. The role of the oscillators is to provide a different time frame to the dynamics, which is reproduced by the linear transformation <italic>&#x3d5;</italic> &#x3d; <italic>&#x3c9;t</italic>. Thus, in the CPGs&#x2019; framework, the integration/derivation over the state variable of the oscillator <italic>&#x3d5;</italic> is conceptually equivalent to the integration over the time. Interestingly, it has been demonstrated in <xref ref-type="bibr" rid="B26">Hayashibe and Shimoda (2014)</xref> that this controller can guarantee energy efficiency during the task realisation in case the quantity that needs to be minimised is the actuation torque.</p>
</sec>
<sec sec-type="results" id="s4">
<title>4 Results</title>
<sec id="s4-1">
<title>4.1 Case1: Monoped</title>
<sec id="s4-1-1">
<title>4.1.1 Adaptation Transient and Energy Efficiency</title>
<p>The goal of the simulations is to analyse the effects of the different feedback in terms of the stability, transient periods, and power injection that is required from the actuator. Four different instances were taken into account for the sensory feedback dynamics, as illustrated in <xref ref-type="fig" rid="F2">Figure&#x20;2</xref>. Although <italic>f</italic>
<sub>2</sub> corresponds to the height of the jump, <italic>f</italic>
<sub>4</sub> is the force that passes through the spring. Then, <italic>f</italic>
<sub>1</sub> and <italic>f</italic>
<sub>3</sub> respectively represent the Tegotae feedback and the feedback that is proposed in <xref ref-type="bibr" rid="B13">Buchli et&#x20;al. (2006)</xref>. Interestingly, both of these share a neuro-mechanical coupling. It is evident that all of them introduce a strong polarisation with the critical points, which is defined as <italic>&#x3d5;</italic>
<sub>0</sub>. The mechanical parameters and the natural length of the spring are <italic>m</italic>&#x20;&#x3d; 0.1&#xa0;kg, <italic>k</italic>&#x20;&#x3d; 5&#xa0;N/m, <italic>c</italic>&#x20;&#x3d; 0.2&#xa0;Ns/m, and <italic>l</italic>
<sub>
<italic>0</italic>
</sub> &#x3d; 1&#xa0;m, respectively. The parameters of the oscillator are &#x3c9; &#x3d; 8&#xa0;rad/s and <italic>&#x3c3;</italic> &#x3d; 2, whose dimensionality is determined on the basis of the feedback law. The initial conditions are respectively <italic>y</italic>
<sub>
<italic>1</italic>
</sub> &#x3d; 0.7&#xa0;m, the velocity is null, and the angle of the oscillator is randomly selected to guarantee a certain robustness with respect to the initial conditions. The actuation parameters and the results of the simulations were obtained from the oscillations in the steady state and are reported in <xref ref-type="table" rid="T1">Table&#x20;1</xref>. The transient period &#x394;<italic>t</italic> is defined at the point at which the limit cycle is reached. The case <italic>f</italic>
<sub>4</sub> is unable to provide a stable orbit. Finally, it is evident that the introduction of the Tegotae feedback is optimal in terms of the synchronisation transient period. In addition, the energy efficiency <italic>E</italic>
<sub>
<italic>e</italic>
</sub> is defined by the limit cycle of the period <inline-formula id="inf4">
<mml:math id="m31">
<mml:mrow>
<mml:msup>
<mml:mi>T</mml:mi>
<mml:mi mathvariant="normal">&#x22c6;</mml:mi>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> with the actuation force <italic>F</italic>
<sub>
<italic>act</italic>
</sub> as follows.<disp-formula id="e28">
<mml:math id="m32">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:msup>
<mml:mi>T</mml:mi>
<mml:mi mathvariant="normal">&#x22c6;</mml:mi>
</mml:msup>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mo>,</mml:mo>
<mml:msup>
<mml:mi>T</mml:mi>
<mml:mi mathvariant="normal">&#x22c6;</mml:mi>
</mml:msup>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mi>E</mml:mi>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(28)</label>
</disp-formula>
<disp-formula id="e29">
<mml:math id="m33">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:munder>
<mml:mstyle displaystyle="true">
<mml:mo>&#x222b;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:msup>
<mml:mi>T</mml:mi>
<mml:mi mathvariant="normal">&#x22c6;</mml:mi>
</mml:msup>
</mml:mrow>
</mml:munder>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>h</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(29)</label>
</disp-formula>
</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Feedback dynamics over the phase <italic>&#x3d5;</italic>. The different lines represent four different instances for the sensory feedback dynamics. <italic>f</italic>
<sub>1</sub>: Tegotae feedback, <italic>f</italic>
<sub>2</sub>: height feedback, <italic>f</italic>
<sub>3</sub>: feedback in <xref ref-type="bibr" rid="B13">Buchli et&#x20;al. (2006)</xref>, <italic>f</italic>
<sub>4</sub>: force feedback.</p>
</caption>
<graphic xlink:href="frobt-08-632804-g002.tif"/>
</fig>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Comparison of performance index, transient period &#x394;<italic>t</italic>, energy efficiency <italic>E</italic>
<sub>
<italic>e</italic>
</sub>, and power injection <italic>J</italic>, for the feedback types on 1D hopping.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Feedback</th>
<th align="center">
<italic>f</italic>
<sub>1</sub>
</th>
<th align="center">
<italic>f</italic>
<sub>2</sub>
</th>
<th align="center">
<italic>f</italic>
<sub>3</sub>
</th>
<th align="center">
<italic>f</italic>
<sub>4</sub>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">A&#xa0;[N]</td>
<td align="center">4</td>
<td align="center">12</td>
<td align="char" char=".">12</td>
<td align="char" char=".">4</td>
</tr>
<tr>
<td align="left">
<italic>&#x3d5;</italic>
<sub>0</sub>, &#x394;&#xa0;[rad]</td>
<td align="center">1.75&#x3c0;0.1&#x3c0;</td>
<td align="center">1.96&#x3c0;0.1&#x3c0;</td>
<td align="center">1.96&#x3c0;0.1&#x3c0;</td>
<td align="center">1.75&#x3c0;0.1&#x3c0;</td>
</tr>
<tr>
<td align="left">&#x394;t&#xa0;[s]</td>
<td align="center">3</td>
<td align="center">4</td>
<td align="char" char=".">5</td>
<td align="center">
<inline-formula id="inf5">
<mml:math id="m34">
<mml:mi mathvariant="normal">&#x2204;</mml:mi>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td align="left">
<italic>E</italic>
<sub>
<italic>e</italic>
</sub>&#xa0;[m/Ws]</td>
<td align="center">1.50</td>
<td align="center">1.16</td>
<td align="char" char=".">1.15</td>
<td align="char" char=".">1.25</td>
</tr>
<tr>
<td align="left">J&#xa0;[W]</td>
<td align="center">5.49</td>
<td align="center">17.69</td>
<td align="char" char=".">20.15</td>
<td align="char" char=".">10.56</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Interestingly, to obtain a similar hopping in terms of the height, the cases <italic>f</italic>
<sub>2</sub> and <italic>f</italic>
<sub>3</sub> are required for a higher amplitude of the actuation&#x20;force.</p>
</sec>
<sec id="s4-1-2">
<title>4.1.2 Robustness and Adaptivity</title>
<p>Second, the case of the Tegotae approach <italic>f</italic>
<sub>1</sub> and the <italic>f</italic>
<sub>3</sub> case that is presented in <xref ref-type="bibr" rid="B13">Buchli et&#x20;al. (2006)</xref> were taken into account. In addition, the adaptivity was evaluated based on the dynamical change in the environment. In particular, at t &#x3d; 5&#xa0;s, the ground level was lowered from 0 to &#x2212;0.6&#xa0;m. The results are depicted in <xref ref-type="fig" rid="F3">Figure&#x20;3</xref>. It is evident that our approach can cope with these variations by performing a proper re-polarisation of the oscillator, even without the adaptation of <italic>&#x3c3;</italic>, <italic>&#x3d5;</italic>
<sub>0</sub>, or &#x394;. It is possible to notice how the Tegotae approach can quickly react to these variations, by modifying the force injection as shown in <xref ref-type="disp-formula" rid="e11">Eq.&#x20;(11)</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Dynamic environment and adaptation process. The ground level was lowered from 0 to &#x2212;0.6&#xa0;m at <italic>t</italic>&#x20;&#x3d; 5&#xa0;s. The upper and lower graphs depict the cases of <italic>f</italic>
<sub>1</sub>: Tegotae feedback and <italic>f</italic>
<sub>3</sub>: feedback in <xref ref-type="bibr" rid="B13">Buchli et&#x20;al. (2006)</xref>. The black and red lines represent the trajectories and force injected, respectively. The Tegotae approach can quickly react to these variations, by modifying the force injection as shown in <xref ref-type="disp-formula" rid="e11">Eq. 11</xref>. The initial state of the monopod robot was the equilibrium point of the spring-mass-damper system. Thus, the height is unchanged while no force is applied.</p>
</caption>
<graphic xlink:href="frobt-08-632804-g003.tif"/>
</fig>
</sec>
</sec>
<sec id="s4-2">
<title>4.2 Case2: Biped</title>
<sec id="s4-2-1">
<title>4.2.1 Gait</title>
<p>The objective in the biped case is to first obtain two different gaits, namely in-phase and anti-phase bipedal hopping. As already stated in <xref ref-type="bibr" rid="B48">Owaki et&#x20;al. (2017)</xref>; <xref ref-type="bibr" rid="B46">Owaki and Ishiguro (2017b)</xref>, for the architecture of the CPGs, the frequency of the oscillation <italic>&#x3c9;</italic> is a useful control variable that can be exploited to introduce a gait transition in the pattern generation. This frequency can be observed as one of the few high-level control variables that are required by CPG architectures, as already presented in <xref ref-type="bibr" rid="B28">Ijspeert (2008)</xref>. Interestingly, our Tegotae control policy can maintain these properties, even without introducing any oscillator couplings, i.e. <inline-formula id="inf6">
<mml:math id="m35">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mrow>
<mml:mn>12</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mrow>
<mml:mn>21</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>Two distinct gaits, in-phase hopping and anti-phase hopping, are reported in <xref ref-type="fig" rid="F4">Figures 4 (Top and Bottom)</xref>. The case of <xref ref-type="fig" rid="F4">Figure&#x20;4 (Top)</xref> is obtained with a frequency <italic>&#x3c9;</italic>
<sub>
<italic>in</italic>
</sub> &#x3d; 6&#xa0;rad/s, while the second case of <xref ref-type="fig" rid="F4">Figure&#x20;4 (Bottom)</xref> is obtained with <italic>&#x3c9;</italic>
<sub>
<italic>anti</italic>
</sub> &#x3d; 7.5&#xa0;rad/s. At first, we determined these parameters by trial and error. Then, we performed a study on the attractors of the dynamics via Lyapunov Exponents; however, this analysis is out of the scope of this article. The values of the mechanical parameters are generally equal to those in the monoped case, with the addition of a spring constant <italic>k</italic>
<sub>
<italic>c</italic>
</sub> &#x3d; 1. The feedback strength was <italic>&#x3c3;</italic> &#x3d; 2.4 to guarantee a higher vertical excursion. We considered a few <italic>&#x3c3;</italic> values, and found that the motion was stable for certain values, while it was unstable for others, suggesting that the value of <italic>&#x3c3;</italic> has an effect on the stability. However, the effect of <italic>&#x3c3;</italic> is not considered in this paper because it out of the scope of this study. The initial conditions are <italic>y</italic>
<sub>1</sub> &#x3d; 0.8&#xa0;m, <italic>y</italic>
<sub>2</sub> &#x3d; 0.7&#xa0;m, the velocities are null, and the angles of the oscillators are selected randomly to guarantee a certain robustness with respect to the initial conditions. These figures represent the mechanical section of the system (heights and forces) and the control section (phases and feedbacks), with the actuation force and Tegotae feedback, respectively.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Hopping gait patterns (Top) in-phase hopping: <italic>&#x3c9;</italic>
<sub>
<italic>in</italic>
</sub> &#x3d; 6&#xa0;rad/s (Bottom) anti-phase hopping: <italic>&#x3c9;</italic>
<sub>
<italic>anti</italic>
</sub> &#x3d; 7.5&#xa0;rad/s. The upper and lower graphs show the mechanical section (heights and forces) and control section (phases and feedbacks), respectively. The blue and red colors represent the left (1) and right legs (2), respectively.</p>
</caption>
<graphic xlink:href="frobt-08-632804-g004.tif"/>
</fig>
<p>Finally, it was evident that by changing the control variable from <italic>&#x3c9;</italic>
<sub>
<italic>in</italic>
</sub> to <italic>&#x3c9;</italic>
<sub>
<italic>anti</italic>
</sub>, it is possible to reproduce a gait transition, as depicted in <xref ref-type="fig" rid="F5">Figure&#x20;5</xref>. As demonstrated, the value is changed at <italic>t</italic>&#x20;&#x3d; 8&#xa0;s and the trend of the actuation forces and feedback are hidden for clarity reasons due to the presence of several transient sections. The motivations for these specific gaits are shown for the different values of <italic>&#x3c9;</italic> that are still an open point thus far. This also considers the fact that due to the random initialisation of the phase angles, the other gates are seldomly shown. These cases can be avoided by constructing a more robust architecture that can integrate several types of sensory feedback.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Hopping gait transition. The frequency <italic>&#x3c9;</italic> is changed from <italic>&#x3c9;</italic>
<sub>
<italic>in</italic>
</sub> to <italic>&#x3c9;</italic>
<sub>
<italic>anti</italic>
</sub> at <italic>t</italic>&#x20;&#x3d; 8&#xa0;s. The upper and lower graphs depict the height of each leg and phase sin<italic>&#x3d5;</italic>
<sub>
<italic>i</italic>
</sub> of each leg, respectively.</p>
</caption>
<graphic xlink:href="frobt-08-632804-g005.tif"/>
</fig>
</sec>
<sec id="s4-2-2">
<title>4.2.2 Robustness and Adaptivity</title>
<p>Finally, in equivalence to the monoped case, the way in which the control policy expressed in <xref ref-type="disp-formula" rid="e11">Eq. 11</xref> can sustain a change in the environmental conditions was also examined for the biped case. As depicted in <xref ref-type="fig" rid="F6">Figure&#x20;6A</xref>, the ground was first lowered to &#x2212;0.6&#xa0;m for both the legs as demonstrated in the monoped case. Meanwhile, the angular frequency was maintained equal to <italic>&#x3c9;</italic>
<sub>
<italic>in</italic>
</sub>. Second, as depicted in <xref ref-type="fig" rid="F6">Figure&#x20;6B</xref>, the ground was lowered again to &#x2212;0.6&#xa0;m for both the legs. Meanwhile, the angular frequency was equal to <italic>&#x3c9;</italic>
<sub>
<italic>anti</italic>
</sub>. The results confirm a good robustness of the control policy to the environmental conditions, which in this case is the ground&#x20;level.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Adaptation to a lower step (Top) In-phase hopping (Bottom) Anti-phase hopping. The ground level was lowered from 0 to &#x2212;0.6&#xa0;m at t &#x3d; 10&#xa0;s.</p>
</caption>
<graphic xlink:href="frobt-08-632804-g006.tif"/>
</fig>
</sec>
</sec>
<sec id="s4-3">
<title>4.3 Optimal Control for the Monoped Case</title>
<p>The optimisation was run for several values of the mass to validate the results for the different feedback dynamics. Meanwhile, all the other parameters were the same as described in the monoped case study. In contrast, the Tegotae controller was applied in <xref ref-type="disp-formula" rid="e11">Eq. (11)</xref> to exploit the adaptivity of the Tegotae feedback.</p>
<p>The values of the weights for the cost functions are reported in <xref ref-type="table" rid="T2">Table&#x20;2</xref> with respect to each simulation to determine the effectiveness of the weights. It follows that the actual effect of the weights is restricted to the power injection by the controller. Meanwhile, the optimal controller does not have access to the energy stored in the spring and the damping system or to the vertical excursion, as shown in <xref ref-type="sec" rid="s10">Supplementary Figures S1&#x2013;S3</xref> in the Supplementary Material (SM). In contrast, the ability of dynamically adapting to the mass changes of the Tegotae controller is verified by the optimal controller as well, as shown in <xref ref-type="fig" rid="F7">Figures 7</xref> (Top) to (Bottom). It is evident that the effect of the first term <inline-formula id="inf7">
<mml:math id="m36">
<mml:mrow>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is sufficient to reproduce, for three different values of masses, to reproduce the effects of the Tegotae control. This term corresponds to the energy consumption of the controller. Therefore, the Tegotae control and an optimal control that attempts to maximise the energy efficiency provide similar results for different masses, thereby validating our hypothesis. Further increments of the mass may require a change in the value of <italic>&#x3c3;</italic> or the use of a non-linear spring to avoid negative values of vertical movements.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Weight values for the cost functions and RMSE <italic>y</italic>, <inline-formula id="inf8">
<mml:math id="m37">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula>, and, <italic>q</italic> for MS.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Simulation</th>
<th align="center">
<italic>m</italic>
</th>
<th align="center">
<inline-formula id="inf9">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th align="center">
<italic>R</italic>
<sub>1</sub>
</th>
<th align="center">
<italic>L</italic>
<sub>1</sub>
</th>
<th align="center">RMSE <italic>y</italic>
</th>
<th align="center">RMSE <inline-formula id="inf10">
<mml:math id="m39">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th align="center">RMSE <italic>q</italic>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">MS 1</td>
<td align="char" char=".">0.1</td>
<td align="center">1<italic>e</italic>1</td>
<td align="center">&#x2212;1<italic>e</italic>1</td>
<td align="center">&#x2212;1<italic>e</italic>1</td>
<td align="char" char=".">0.03</td>
<td align="char" char=".">0.34</td>
<td align="char" char=".">0.58</td>
</tr>
<tr>
<td align="left">MS 2</td>
<td align="char" char=".">0.1</td>
<td align="char" char=".">0</td>
<td align="center">&#x2212;1<italic>e</italic>1</td>
<td align="center">&#x2212;1<italic>e</italic>1</td>
<td align="char" char=".">0.03</td>
<td align="char" char=".">0.34</td>
<td align="char" char=".">0.58</td>
</tr>
<tr>
<td align="left">MS 3</td>
<td align="char" char=".">0.1</td>
<td align="center">1<italic>e</italic>1</td>
<td align="char" char=".">0</td>
<td align="center">&#x2212;1<italic>e</italic>1</td>
<td align="char" char=".">0.03</td>
<td align="char" char=".">0.34</td>
<td align="char" char=".">0.58</td>
</tr>
<tr>
<td align="left">MS 4</td>
<td align="char" char=".">0.1</td>
<td align="center">1<italic>e</italic>1</td>
<td align="center">&#x2212;1<italic>e</italic>1</td>
<td align="char" char=".">0</td>
<td align="char" char=".">0.03</td>
<td align="char" char=".">0.34</td>
<td align="char" char=".">0.58</td>
</tr>
<tr>
<td align="left">MS 5</td>
<td align="char" char=".">0.3</td>
<td align="center">1<italic>e</italic>1</td>
<td align="char" char=".">0</td>
<td align="char" char=".">0</td>
<td align="char" char=".">0.03</td>
<td align="char" char=".">0.34</td>
<td align="char" char=".">0.58</td>
</tr>
<tr>
<td align="left">MS 6</td>
<td align="char" char=".">0.6</td>
<td align="center">1<italic>e</italic>1</td>
<td align="char" char=".">0</td>
<td align="char" char=".">0</td>
<td align="char" char=".">0.03</td>
<td align="char" char=".">0.34</td>
<td align="char" char=".">0.58</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Results of multiple shooting methods. The blue and solid red dotted lines represent the designed optimal controller (MS method) and Tegotae controller, respectively Top case MS1 in <xref ref-type="table" rid="T2">Table&#x20;2</xref>: <italic>m</italic>&#x20;&#x3d; 0.1 (Middle) case MS5 in <xref ref-type="table" rid="T2">Table&#x20;2</xref>: <italic>m</italic>&#x20;&#x3d; 0.3 Bottom case MS6 in <xref ref-type="table" rid="T2">Table&#x20;2</xref>: <italic>m</italic>&#x20;&#x3d; 0.6. Not only was the Tegotae control action extremely similar to the MS optimal control in all the cases, but also the position and velocity profiles demonstrated certain similarities.</p>
</caption>
<graphic xlink:href="frobt-08-632804-g007.tif"/>
</fig>
<p>Not only was the Tegotae control action extremely similar to the MS optimal control (see the <xref ref-type="sec" rid="s10">Supplementary Material</xref>) in all the cases, but also the position and velocity profiles demonstrated certain similarities. In all the MS cases, the root mean squared errors (RMSE) were found to be similar, as reported in <xref ref-type="table" rid="T2">Table&#x20;2</xref>, as expected from previous considerations. Finally, for all the cases considered in the MS examples, the energy efficiency of the optimal controller as expressed in <xref ref-type="disp-formula" rid="e29">Eq. (29)</xref> converged to a value similar to that of the Tegotae controller, whose value was determined considering 1&#xa0;m as the maximum height reached, for comparison purposes. The convergence is reported in <xref ref-type="fig" rid="F8">Figure&#x20;8</xref> for MS1 and leads to a final RMSE of 0.22. This seems to limit to the efficiency given the physical constraints of the system. Moreover, increasing the weight <italic>Q</italic> slightly increases the efficiency.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Energy-efficiency convergence in the MS method through comparison with the Tegotae feedback&#x20;case.</p>
</caption>
<graphic xlink:href="frobt-08-632804-g008.tif"/>
</fig>
<p>These results represent the MS case alone. The SS (see the <xref ref-type="sec" rid="s10">Supplementary Material</xref>) has several practical drawbacks, which motivates this choice. First, it requires extremely high weights for the sensitivity function of the final conditions and the smoothness of the control policy. The conditions are automatically satisfied by the continuity constraints in the MS. Second, the convergence is more difficult to obtain. The FHOC for the SS method is formulated by using the norm notation and the additional weights to guarantee a sensitivity to the final conditions and control policy.<disp-formula id="e30">
<mml:math id="m40">
<mml:mrow>
<mml:munder>
<mml:mrow>
<mml:mtext>min</mml:mtext>
</mml:mrow>
<mml:mi>q</mml:mi>
</mml:munder>
<mml:msubsup>
<mml:mrow>
<mml:mo>&#x222b;</mml:mo>
</mml:mrow>
<mml:mn>0</mml:mn>
<mml:mi>T</mml:mi>
</mml:msubsup>
<mml:msubsup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b3;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>y</mml:mi>
<mml:msubsup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>T</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>T</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(30)</label>
</disp-formula>Subject to<disp-formula id="equ1">
<mml:math id="m41">
<mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtext>Initial&#xa0;Value&#xa0;Constraints</mml:mtext>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="equ2">
<mml:math id="m42">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#xa8;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>q</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>g</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>q</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mi>t</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>T</mml:mi>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtext>ODE&#xa0;Constraint</mml:mtext>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="equ3">
<mml:math id="m43">
<mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>y</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>l</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>y</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2026;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>N</mml:mi>
<mml:mo>,</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtext>Inequality&#xa0;Constraints</mml:mtext>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>In our case, <italic>&#x3b3;</italic>
<sub>1</sub> &#x3d; 1<italic>e</italic>4 and <italic>F</italic>
<sub>2</sub> &#x3d; 1<italic>e</italic>10. As previously mentioned, these values are extremely high in comparison with the remaining weights in the cost function as presented in <xref ref-type="table" rid="T2">Table&#x20;2</xref>. Meanwhile, for the MS case, the weights remain the same as MS five in <xref ref-type="table" rid="T3">Table&#x20;3</xref>. Interestingly, it has not been a trivial fact to obtain similar results between the two optimal controllers. It is possible to obtain similar control trends with respect to the MS case, as shown in <xref ref-type="fig" rid="F9">Figure&#x20;9</xref>. (Top) and (Bottom) However, there are also cases that are similar to the Tegotae controller, as shown in <xref ref-type="sec" rid="s10">Supplementary Figures S5,S6</xref> in the SM; this is achieved by varying the values of the weights. For the SS case, the cost function is sensitive to the terms that are proper to the monopod cost function in <xref ref-type="disp-formula" rid="e29">Eq. 29</xref> and the spring&#x20;force.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Weights values for the cost functions for the MS-SS.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Simulation</th>
<th align="center">
<italic>m</italic>
</th>
<th align="center">
<inline-formula id="inf11">
<mml:math id="m44">
<mml:mrow>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th align="center">
<italic>R</italic>
<sub>2</sub>
</th>
<th align="center">
<italic>L</italic>
<sub>2</sub>
</th>
<th align="center">
<italic>S</italic>
<sub>2</sub>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">MS-SS 1</td>
<td align="char" char=".">0.3</td>
<td align="center">1<italic>e</italic>1</td>
<td align="center">&#x2212;1<italic>e</italic>1</td>
<td align="center">&#x2212;1<italic>e</italic>1</td>
<td align="center">1<italic>e</italic>3</td>
</tr>
<tr>
<td align="left">MS-SS 2</td>
<td align="char" char=".">0.6</td>
<td align="center">1<italic>e</italic>1</td>
<td align="center">&#x2212;1<italic>e</italic>1</td>
<td align="center">&#x2212;1<italic>e</italic>1</td>
<td align="center">1<italic>e</italic>3</td>
</tr>
<tr>
<td align="left">MS-SS 3</td>
<td align="char" char=".">0.4</td>
<td align="center">1<italic>e</italic>1</td>
<td align="center">&#x2212;1<italic>e</italic>2</td>
<td align="center">&#x2212;1<italic>e</italic>2</td>
<td align="center">1<italic>e</italic>3</td>
</tr>
<tr>
<td align="left">MS-SS 4</td>
<td align="char" char=".">0.6</td>
<td align="center">1<italic>e</italic>1</td>
<td align="center">&#x2212;1<italic>e</italic>2</td>
<td align="center">&#x2212;1<italic>e</italic>2</td>
<td align="center">1<italic>e</italic>3</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Results of multiple shooting-single shooting (Top) case MS-SS1 in <xref ref-type="table" rid="T3">:Table&#x20;3</xref> <italic>m</italic>&#x20;&#x3d; 0.3 (Bottom) case MS-SS2 in <xref ref-type="table" rid="T3">:Table&#x20;3</xref> <italic>m</italic>&#x20;&#x3d; 0.6.</p>
</caption>
<graphic xlink:href="frobt-08-632804-g009.tif"/>
</fig>
<p>The MS routine is solved by using the interior-point method that is provided by the MATLAB built-in function FMINCON. In contrast, the SS routine is solved by using the BFGS method and the SQP that is designed on the material, provided by <xref ref-type="bibr" rid="B19">Fagiano (2019)</xref>. With regard to the integration of the dynamics, the time interval was split into 40&#xa0;nodes with 2&#xa0;points per sub-interval for the MS case. Meanwhile, a sampling time of 0.01&#xa0;s was used for the SS case. In both cases, the integration of the dynamics was conducted using an explicit Runge-Kutta method with an order of four since the restricted dynamics were non-stiff. The step size was 0.01&#xa0;s in both the methods.</p>
</sec>
</sec>
<sec sec-type="discussion" id="s5">
<title>5 Discussion</title>
<p>The main contribution of this study is to propose a control policy with a reflex-like actuation (<xref ref-type="disp-formula" rid="e11">Eq. (11)</xref>) for the Tegotae-based feedback law in the CPG in such a way that the controller fruitfully exploits the <italic>embodiment</italic> (<xref ref-type="bibr" rid="B51">Pfeifer and Bongard, 2006</xref>; <xref ref-type="bibr" rid="B52">Pfeifer et&#x20;al., 2007</xref>). For the validation of the proposed method, we first demonstrated the energy efficiency of the monopod model as well as its robustness and adaptability using the controller. Then, we demonstrated the gait transition for the bipedal model with its robustness and adaptability. Based on the optimal control theory, we designed an optimal controller and then compared it with the Tegotae-based control input. The results indicate the Tegotae-based feedback with reflex-like actuation results for optimal and energy-efficient motion. This suggests the first evidence concerning the optimal energy efficiency for the Tegotae approach.</p>
<p>This study is the first attempt to analyse the optimal energy efficiency along with the adaptivity of the Tegotae approach. Previous studies (<xref ref-type="bibr" rid="B48">Owaki et&#x20;al., 2017</xref>) have mainly focused on the temporal (timing/phase) modulation in the oscillators by the Tegotae feedback on GPG-based models. The proposed reflex-like actuation can modulate the &#x201c;amplitude&#x201d; of the actuation via <italic>F</italic>
<sub>
<italic>a</italic>
</sub> function (<xref ref-type="disp-formula" rid="e11">Eq. (11)</xref>), depending on sensory feedback <italic>F</italic>
<sub>
<italic>k</italic>
</sub>. As presented in <xref ref-type="table" rid="T1">Table&#x20;1</xref>, in comparison with the previous methods, the introduction of the Tegotae feedback <italic>f</italic>
<sub>1</sub> was optimal in terms of the transient period for synchronisation and energy efficiency. The reflex-like pathway (<xref ref-type="fig" rid="F1">Figure&#x20;1A</xref>) resulted in a rapid response (fast control loop) on motion generation, leading to the first convergent time in <xref ref-type="table" rid="T1">Table&#x20;1</xref>. Furthermore, the proposed reflex-like actuation (<xref ref-type="disp-formula" rid="e11">Eq. (11)</xref>) induced by the Tegotae feedback in the CPG could generate an input (<xref ref-type="fig" rid="F7">Figures 7</xref>, <xref ref-type="fig" rid="F9">9</xref>) identical to that of the optimally designed controller, resulting in energy-efficient motion, as presented in <xref ref-type="table" rid="T1">Table&#x20;1</xref>. As discussed in <xref ref-type="sec" rid="s3">Section. 3</xref>, the Tegotae approach has similarities (<xref ref-type="disp-formula" rid="e26">Eq. (26)</xref>) with the <italic>tacit learning</italic> frameworks in <xref ref-type="bibr" rid="B26">Hayashibe and Shimoda (2014)</xref>. Energy efficiency is also achieved by the accumulation of a quantity that needs to be minimised when directly integrating the state variable. These facts suggest that our control policy, i.e. reflex-like actuation with the Tegotae-based proprioceptive feedback in the CPG, accomplishes optimal energy-efficient motion through the dynamical learning process along with the interaction between the controller, body, and environments (<xref ref-type="bibr" rid="B51">Pfeifer and Bongard, 2006</xref>; <xref ref-type="bibr" rid="B52">Pfeifer et&#x20;al., 2007</xref>).</p>
<p>The reflex-based leg coordination models (<xref ref-type="bibr" rid="B18">Ekeberg and Pearson, 2005</xref>; <xref ref-type="bibr" rid="B38">Manoonpong et&#x20;al., 2007</xref>; <xref ref-type="bibr" rid="B36">Lewinger and Quinn, 2011</xref>; <xref ref-type="bibr" rid="B54">Schilling et&#x20;al., 2013</xref>; <xref ref-type="bibr" rid="B16">D&#xfc;rr et&#x20;al., 2019</xref>) and reflex-like feedback integration into CPG (<xref ref-type="bibr" rid="B1">Ajallooeian et&#x20;al., 2013</xref>; <xref ref-type="bibr" rid="B17">Dzeladini et&#x20;al., 2014</xref>; <xref ref-type="bibr" rid="B37">Li et&#x20;al., 2014</xref>) have been studied in the past two&#xa0;decades. Pioneering research on &#x201c;event-driven&#x201d; reflex models in cats (<xref ref-type="bibr" rid="B18">Ekeberg and Pearson, 2005</xref>) and insects (<xref ref-type="bibr" rid="B36">Lewinger and Quinn, 2011</xref>; <xref ref-type="bibr" rid="B54">Schilling et&#x20;al., 2013</xref>; <xref ref-type="bibr" rid="B16">D&#xfc;rr et&#x20;al., 2019</xref>) has been conducted, successfully reproducing various aspects of animal inter- and intra-leg coordination during locomotion. <xref ref-type="bibr" rid="B38">Manoonpong et&#x20;al. (2007)</xref> demonstrated that a reflex-based neural controller could achieve stable and fast bipedal walking. Following the pioneering work integrating a CPG with reflex models (<xref ref-type="bibr" rid="B31">Kimura et&#x20;al., 1999</xref>), similar approaches have been proposed. <xref ref-type="bibr" rid="B1">Ajallooeian et&#x20;al. (2013)</xref>; <xref ref-type="bibr" rid="B37">Li et&#x20;al. (2014)</xref> also proposed to integrate a CPG with &#x201c;event-driven&#x201d; reflex models for adaptability against perturbations and environmental changes; One of characteristic approaches in this line, <xref ref-type="bibr" rid="B17">Dzeladini et&#x20;al. (2014)</xref> introduced CPG as feed-forward components in reflex-based neuromuscular models for human walking, confirming the idea of using CPGs as feedback predictors (<xref ref-type="bibr" rid="B33">Kuo (2002)</xref>) from the viewpoint of gait modulation. In our work, the CPG oscillator is not a feedback predictor, but can be considered as a representation of the movement (phase <italic>&#x3d5;</italic>
<sub>
<italic>i</italic>
</sub>), that is, an internal model. In the Tegotae approcah, the Tegotae function <italic>T</italic>
<sub>
<italic>i</italic>
</sub>(<italic>&#x3d5;</italic>
<sub>
<italic>i</italic>
</sub>,<italic>F</italic>
<sub>
<italic>k</italic>
</sub>) is defined as the product of the function of intended motor command <italic>C</italic>(<italic>&#x3d5;</italic>
<sub>
<italic>i</italic>
</sub>) and sensory information <italic>S</italic>(<italic>F</italic>
<sub>
<italic>k</italic>
</sub>); hence, our reflex-like actuation always modulates the motion based on the Tegotae feedback <italic>f</italic>
<sub>
<italic>i</italic>
</sub>, which increases the value of the Tegotae function <italic>T</italic>
<sub>
<italic>i</italic>
</sub>(<italic>&#x3d5;</italic>
<sub>
<italic>i</italic>
</sub>,<italic>F</italic>
<sub>
<italic>k</italic>
</sub>), leading to its adaptability and optimal energy efficiency, as mentioned in previous paragraph.</p>
<p>Past studies that have used the Tegotae approach (<xref ref-type="bibr" rid="B47">Owaki et&#x20;al., 2012</xref>; <xref ref-type="bibr" rid="B46">Owaki and Ishiguro, 2017b</xref>; <xref ref-type="bibr" rid="B48">Owaki et&#x20;al., 2017</xref>) have demonstrated adaptability and behavioural diversity for reproducing animal-like legged locomotion. For quadruped locomotion, the simple and local sensory feedback law in the CPG reproduced the adaptability against the change in mass distribution, which resulted in horse-like or primate-like walking patterns, and a spontaneous gait transition, from walking to trotting and galloping, in response to the locomotion speed. These studies for quadruped robots provide a basis for establishing a design scheme based on the Tegotae approach. For hexapod locomotion, <xref ref-type="bibr" rid="B48">Owaki et&#x20;al. (2017)</xref> designed a minimal model for the inter-limb coordination in a systematic manner based on the Tegotae concept, successfully reproducing the various aspects of the insect locomotion patterns, which includes adaptability to changes in the body properties, e.g. leg amputation. In line with these studies, this investigation also successfully reproduces the adaptability (<xref ref-type="fig" rid="F3">Figures 3</xref>, <xref ref-type="fig" rid="F6">6</xref>), and behavioural diversity (<xref ref-type="fig" rid="F4">Figures 4</xref>, <xref ref-type="fig" rid="F5">5</xref>) as well as the energy efficiency. As discussed in previous studies, in the Tegotae approach, the main aim of designing the Tegotae function is to consider the physical consistency of the action and reaction for the desired motion, and to design the Tegotae function such that its value increases in such cases. Once such a Tegotae function is designed, it is possible to modify the control variables in a situation-dependent manner by increasing the value of the Tegotae function as a feedback term <inline-formula id="inf12">
<mml:math id="m45">
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>S</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>. Therefore, the Tegotae approach enables the design of an autonomous decentralised controller in a systematic manner, by designing the Tegotae function in line with the desired motions.</p>
<p>This study proposes a reflex-like actuation for the Tegotae-based feedback law in the CPG. This is a significant contribution for the actuation and sensory feedback on the adaptation process to the environment and the optimisation process for energy efficiency. However, one of the limitations of this study is that we did not test the applicability of the Tegotae approach to the real-world environment with a physical robot. In addition, it is extremely difficult to perfectly model the dynamics in the real-world environment. One of the key aspects based on the Tegotae approach is the verification in the real world as shown in <xref ref-type="bibr" rid="B47">Owaki et&#x20;al. (2012)</xref>; <xref ref-type="bibr" rid="B46">Owaki and Ishiguro (2017b)</xref>; <xref ref-type="bibr" rid="B48">Owaki et&#x20;al. (2017)</xref>. Instead, we analysed the Tegotae control by using the optimal control theory and provided evidence concerning the optimal control input. Regarding the energy efficiency of tacit learning in the real-world environment, it has been verified by achieving a task with a redundant arm in <xref ref-type="bibr" rid="B27">Hayashibe and Shimoda (2018)</xref>. One potential future direction is to apply our control policy to a robot with more degrees of freedom that performs more complicated tasks. Our control policy is compatible with the force/torque-based control of a physical robot, which is a promising direction of study for future research.</p>
</sec>
</body>
<back>
<sec id="s6">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="sec" rid="s10">Supplementary Material</xref>, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s7">
<title>Author Contributions</title>
<p>DO and MH conceived the research and managed the data collection. RZ designed the model and controllers and conducted the simulations. RZ and DO conducted the analyses. All of the authors wrote the manuscript together.</p>
</sec>
<sec id="s8">
<title>Funding</title>
<p>We gratefully acknowledge the support from the JSPS KAKENHI (grant number JP17KK0109, 18H01399, and 20H04260).</p>
</sec>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of Interest</title>
<p>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.</p>
</sec>
<sec id="s10">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/frobt.2021.632804/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/frobt.2021.632804/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="DataSheet1.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ajallooeian</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Gay</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Tuleu</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Spr&#xf6;witz</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Ijspeert</surname>
<given-names>A. J.</given-names>
</name>
</person-group> (<year>2013</year>). &#x201c;<article-title>Modular Control of Limit Cycle Locomotion over Unperceived Rough Terrain</article-title>,&#x201d; in <conf-name>2013 IEEE/RSJ International Conference on Intelligent Robots and Systems</conf-name>, <conf-loc>Tokyo, Japan</conf-loc>, <conf-date>November 3&#x2013;7, 2013</conf-date>. <fpage>3390</fpage>&#x2013;<lpage>3397</lpage>. <pub-id pub-id-type="doi">10.1109/iros.2013.6696839</pub-id> </citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ambe</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Aoi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Timo</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Manoonpong</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Worgotter</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Matsuno</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Simple Analytical Model Reveals the Functional Role of Embodied Sensorimotor Interaction in Hexapod Gaits</article-title>. <source>PLOS One</source> <volume>13</volume>, <fpage>e0192469</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0192469</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aoi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Tsuchiya</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Locomotion Control of a Biped Robot Using Nonlinear Oscillators</article-title>. <source>Auton. Robot</source> <volume>19</volume>, <fpage>219</fpage>&#x2013;<lpage>232</lpage>. <pub-id pub-id-type="doi">10.1007/s10514-005-4051-1</pub-id> </citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aoi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Tsuchiya</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Stability Analysis of a Simple Walking Model Driven by an Oscillator with a Phase Reset Using Sensory Feedback</article-title>. <source>IEEE Trans. Robotics</source> <volume>22</volume>, <fpage>391</fpage>&#x2013;<lpage>397</lpage>. <pub-id pub-id-type="doi">10.1109/TRO.2006.870671</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aoi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ogihara</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Funato</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Sugimoto</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Tsuchiya</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Evaluating Functional Roles of Phase Resetting in Generation of Adaptive Human Bipedal Walking with a Physiologically Based Model of the Spinal Pattern Generator</article-title>. <source>Biol. Cybern.</source> <volume>102</volume>, <fpage>373</fpage>&#x2013;<lpage>387</lpage>. <pub-id pub-id-type="doi">10.1007/s00422-010-0373-y</pub-id> </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aoi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Yamashita</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Tsuchiya</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Hysteresis in the Gait Transition of a Quadruped Investigated Using Simple Body Mechanical and Oscillator Network Models</article-title>. <source>Phys. Rev. E.</source> <volume>83</volume>, <fpage>061909</fpage>. <pub-id pub-id-type="doi">10.1103/physreve.83.061909</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aoi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Katayama</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Fujiki</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Tomita</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Funato</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Yamashita</surname>
<given-names>T.</given-names>
</name>
<etal/>
</person-group> (<year>2013</year>). <article-title>A Stability-Based Mechanism for Hysteresis in the Walk-Trot Transition in Quadruped Locomotion</article-title>. <source>J.&#x20;R. Soc. Interf.</source> <volume>10</volume>, <fpage>20120908</fpage>. <pub-id pub-id-type="doi">10.1098/rsif.2012.0908</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>B&#xe4;ssler</surname>
<given-names>U.</given-names>
</name>
<name>
<surname>Wegner</surname>
<given-names>U. T. A.</given-names>
</name>
</person-group> (<year>1983</year>). <article-title>Motor Output of the Denervated Thoracic Ventral Nerve Cord in the Stick Insect Carausius Morosus</article-title>. <source>J.&#x20;Exp. Biol.</source> <volume>105</volume>, <fpage>127</fpage>&#x2013;<lpage>145</lpage>. <pub-id pub-id-type="doi">10.1242/jeb.105.1.127</pub-id> </citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>B&#xe4;ssler</surname>
<given-names>U.</given-names>
</name>
</person-group> (<year>1986</year>). <article-title>On the Definition of Central Pattern Generator and its Sensory Control</article-title>. <source>Biol. Cybern.</source> <volume>54</volume>, <fpage>65</fpage>&#x2013;<lpage>69</lpage>. <pub-id pub-id-type="doi">10.1007/BF00337116</pub-id> </citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Berenz</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Alnajjar</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Hayashibe</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Shimoda</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2014</year>). &#x201c;<article-title>Generalization of the Tacit Learning Controller Based on Periodic Tuning Functions</article-title>,&#x201d; in <conf-name>Proc. of 5th IEEE RAS/EMBS International Conference on Biomedical Robotics and Biomechatronics</conf-name>, <conf-loc>Sao Paulo, Brazil</conf-loc>, <conf-date>August 12&#x2013;15, 2014</conf-date>. <fpage>15144072</fpage>. <pub-id pub-id-type="doi">10.1109/biorob.2014.6913894</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Berenz</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Alnajjar</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Hayashibe</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Shimoda</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2015</year>). &#x201c;<article-title>Tacit Learning for Emergence of Task-Related Behaviour through Signal Accumulation</article-title>,&#x201d; in <source>Emergent Trends in Robotics and Intelligent Systems</source>, <publisher-loc>New York</publisher-loc>: <publisher-name>Springer</publisher-name>, <fpage>31</fpage>&#x2013;<lpage>38</lpage>. <pub-id pub-id-type="doi">10.1007/978-3-319-10783-7_3</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bock</surname>
<given-names>H. G.</given-names>
</name>
<name>
<surname>Plitt</surname>
<given-names>K. J.</given-names>
</name>
</person-group> (<year>1984</year>). <article-title>A Multiple Shooting Algorithm for Direct Solution of Optimal Control Problems &#x2a;</article-title>. <source>IFAC Proc. Volumes</source> <volume>17</volume>, <fpage>1603</fpage>&#x2013;<lpage>1608</lpage>. <pub-id pub-id-type="doi">10.1016/s1474-6670(17)61205-9</pub-id> </citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Buchli</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Iida</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Ijspeert</surname>
<given-names>A. J.</given-names>
</name>
</person-group> (<year>2006</year>). &#x201c;<article-title>Finding Resonance: Adaptive Frequency Oscillators for Dynamic Legged Locomotion</article-title>,&#x201d; in <conf-name>IEEE/RSJ International Conference on Intelligent Robots and Systems</conf-name>, <conf-loc>Beijing, China</conf-loc>, <conf-date>October 9&#x2013;15, 2006</conf-date>. <fpage>3903</fpage>&#x2013;<lpage>3909</lpage>. </citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dasgupta</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Goldschmidt</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Worgotter</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Manoonpong</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Distributed Recurrent Neural Forward Models with Synaptic Adaptation and Cpg-Based Control for Complex Behaviors of Walking Robots</article-title>. <source>Front. Neurorobot.</source> <volume>9</volume>, <fpage>10</fpage>. <pub-id pub-id-type="doi">10.3389/fnbot.2015.00010</pub-id> </citation>
</ref>
<ref id="B15">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Diehl</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Bock</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Diedam</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wieber</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2005</year>). &#x201C;<article-title>Fast Direct Multiple Shooting Algorithms for Optimal Robot Control</article-title>,&#x201D; in <source>Fast Motions in Biomechanics and Robotics Lecture Notes in Control and Information Sciences</source>. Editors <person-group person-group-type="editor">
<name>
<surname>Diehl</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Mombaur</surname>
<given-names>K.</given-names>
</name>
</person-group>, (<publisher-loc>Berlin, Heidelberg</publisher-loc>: <publisher-name>Springer</publisher-name>), <volume>340</volume>. <pub-id pub-id-type="doi">10.1007/978-3-540-36119-0_4</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>D&#xfc;rr</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Arena</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Cruse</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Dallmann</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Drimus</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Hoinville</surname>
<given-names>T.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Integrative Biomimetics of Autonomous Hexapedal Locomotion</article-title>. <source>Front. Neurorobot.</source> <volume>13</volume>, <fpage>88</fpage>. <pub-id pub-id-type="doi">10.3389/fnbot.2019.00088</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dzeladini</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>van den Kieboom</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ijspeert</surname>
<given-names>A. J.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>The Contribution of a Central Pattern Generator in a Reflex-Based Neuromuscular Model</article-title>. <source>Front. Hum. Neurosci.</source> <volume>8</volume>, <fpage>371</fpage>. <pub-id pub-id-type="doi">10.3389/fnhum.2014.00371</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ekeberg</surname>
<given-names>&#xd6;.</given-names>
</name>
<name>
<surname>Pearson</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Computer Simulation of Stepping in the Hind Legs of the Cat: An Examination of Mechanisms Regulating the Stance-To-Swing Transition</article-title>. <source>J.&#x20;Neurophysiol.</source> <volume>94</volume>, <fpage>4256</fpage>&#x2013;<lpage>4268</lpage>. <pub-id pub-id-type="doi">10.1152/jn.00065.2005</pub-id> </citation>
</ref>
<ref id="B19">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Fagiano</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2019</year>). <source>Constrained Numerical Optimization for Estimation and Control (Lecture Notes)</source>. <publisher-loc>Milan, Italy</publisher-loc>: <publisher-name>Politecnico di Milano</publisher-name>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fukuoka</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Kimura</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Cohen</surname>
<given-names>A. H.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Adaptive Dynamic Walking of a Quadruped Robot on Irregular Terrain Based on Biological Concepts</article-title>. <source>Int. J.&#x20;Robotics Res.</source> <volume>22</volume>, <fpage>187</fpage>&#x2013;<lpage>202</lpage>. <pub-id pub-id-type="doi">10.1177/0278364903022003004</pub-id> </citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fukuoka</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Habu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Fukui</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>A Simple Rule for Quadrupedal Gait Generation Determined by Leg Loading Feedback: a Modeling Study</article-title>. <source>Sci. Rep.</source> <volume>5</volume>, <fpage>8169</fpage>. <pub-id pub-id-type="doi">10.1038/srep08169</pub-id> </citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Goldschmidt</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Worgotter</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Manoonpong</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Biologically-inspired Adaptive Obstacle Negotiation Behavior of Hexapod Robots</article-title>. <source>Front. Neurorobot.</source> <volume>8</volume>, <fpage>3</fpage>. <pub-id pub-id-type="doi">10.3389/fnbot.2014.00003</pub-id> </citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Grillner</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>1975</year>). <article-title>Locomotion in Vertebrates: Central Mechanisms and Reflex Interaction</article-title>. <source>Physiol. Rev.</source> <volume>55</volume>, <fpage>247</fpage>&#x2013;<lpage>304</lpage>. <pub-id pub-id-type="doi">10.1152/physrev.1975.55.2.247</pub-id> </citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Grillner</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>1985</year>). <article-title>Neurobiological Bases of Rhythmic Motor Acts in Vertebrates</article-title>. <source>Science</source> <volume>228</volume>, <fpage>143</fpage>&#x2013;<lpage>149</lpage>. <pub-id pub-id-type="doi">10.1126/science.3975635</pub-id> </citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Grinke</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Tetzlaff</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Worgotter</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Manoonpong</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Synaptic Plasticity in a Recurrent Neural Network for Versatile and Adaptive Behaviors of a Walking Robot</article-title>. <source>Front. Neurorobot.</source> <volume>9</volume>, <fpage>11</fpage>. <pub-id pub-id-type="doi">10.3389/fnbot.2015.00011</pub-id> </citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hayashibe</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Shimoda</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Synergetic Motor Control Paradigm for Optimizing Energy Efficiency of Multijoint Reaching via Tacit Learning</article-title>. <source>Front. Comput. Neurosci.</source> <volume>8</volume>, <fpage>21</fpage>. <pub-id pub-id-type="doi">10.3389/fncom.2014.00021</pub-id> </citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hayashibe</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Shimoda</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Synergetic Learning Control Paradigm for Redundant Robot to Enhance Error-Energy Index</article-title>. <source>IEEE Trans. Cogn. Dev. Syst.</source> <volume>10</volume>, <fpage>573</fpage>&#x2013;<lpage>584</lpage>. <pub-id pub-id-type="doi">10.1109/tcds.2017.2697904</pub-id> </citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ijspeert</surname>
<given-names>A. J.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Central Pattern Generators for Locomotion Control in Animals and Robots: A Review</article-title>. <source>Neural networks</source> <volume>21</volume>, <fpage>642</fpage>&#x2013;<lpage>653</lpage>. <pub-id pub-id-type="doi">10.1016/j.neunet.2008.03.014</pub-id> </citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kano</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Yoshizawa</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Ishiguro</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Tegotae -based Decentralised Control Scheme for Autonomous Gait Transition of Snake-like Robots</article-title>. <source>Bioinspir. Biomim.</source> <volume>12</volume>, <fpage>046009</fpage>. <pub-id pub-id-type="doi">10.1088/1748-3190/aa7725</pub-id> </citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kano</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Wakimoto</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Sato</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Shinohara</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Ishiguro</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Designing Higher Fourier Harmonics of Tegotae Function Using Genetic Algorithm&#x2014;A Case Study with an Earthworm Locomotion</article-title>. <source>Bioinspir. Biomim.</source> <volume>14</volume>, <fpage>05400</fpage>. <pub-id pub-id-type="doi">10.1088/1748-3190/ab2fab</pub-id> </citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kimura</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Akiyama</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Sakurama</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Realization of Dynamic Walking and Running of the Quadruped Using Neural Oscillator</article-title>. <source>Autonomous Robots</source> <volume>7</volume>, <fpage>247</fpage>&#x2013;<lpage>258</lpage>. <pub-id pub-id-type="doi">10.1023/a:1008924521542</pub-id> </citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kimura</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Fukuoka</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Cohen</surname>
<given-names>A. H.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Adaptive Dynamic Walking of a Quadruped Robot on Natural Ground Based on Biological Concepts</article-title>. <source>Int. J.&#x20;Robotics Res.</source> <volume>26</volume>, <fpage>457</fpage>&#x2013;<lpage>490</lpage>. <pub-id pub-id-type="doi">10.1177/0278364907078089</pub-id> </citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kuo</surname>
<given-names>A. D.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>The Relative Roles of Feedforward and Feedback in the Control of Rhythmic Movements</article-title>. <source>Motor Control</source> <volume>6</volume>, <fpage>129</fpage>&#x2013;<lpage>145</lpage>. <pub-id pub-id-type="doi">10.1152/jn.00065.200510.1123/mcj.6.2.129</pub-id> </citation>
</ref>
<ref id="B34">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kuramoto</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>1984</year>). <source>Chemical Oscillations, Waves, and Turbulence</source>. <publisher-loc>New York</publisher-loc>: <publisher-name>Springer-Verlag</publisher-name>. <pub-id pub-id-type="doi">10.1007/978-3-642-69689-3</pub-id> </citation>
</ref>
<ref id="B35">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kuramoto</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2003</year>). &#x201C;<article-title>Collective behavior of coupled phase oscillators</article-title>,&#x201D; in <source>The Handbook of Brain Theory and Neural Networks</source>. Editors <person-group person-group-type="editor">
<name>
<surname>Arbib</surname>
<given-names>MA.</given-names>
</name>
</person-group> (<publisher-name>MIT Press</publisher-name>), <fpage>223</fpage>&#x2013;<lpage>226</lpage>.</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lewinger</surname>
<given-names>W. A.</given-names>
</name>
<name>
<surname>Quinn</surname>
<given-names>R. D.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Neurobiologically&#x2010;based Control System for an Adaptively Walking Hexapod</article-title>. <source>Ind. Robot</source> <volume>38</volume>, <fpage>258</fpage>&#x2013;<lpage>263</lpage>. <pub-id pub-id-type="doi">10.1108/01439911111122752</pub-id> </citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Hildre</surname>
<given-names>H. P.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>An Approach for Adaptive Limbless Locomotion Using a Cpg-Based Reflex Mechanism</article-title>. <source>J.&#x20;Bionic Eng.</source> <volume>11</volume>, <fpage>389</fpage>&#x2013;<lpage>399</lpage>. <pub-id pub-id-type="doi">10.1016/s1672-6529(14)60052-4</pub-id> </citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Manoonpong</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Geng</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Kulvicius</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Porr</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>W&#xf6;rg&#xf6;tter</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Adaptive, Fast Walking in a Biped Robot under Neuronal Control and Learning</article-title>. <source>PLoS Comput. Biol.</source> <volume>3</volume>, <fpage>e134</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pcbi.0030134</pub-id> </citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Manoonpong</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Pasemann</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Worgotter</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Sensor-driven Neural Control for Omnidirectional Locomotion and Versatile Reactive Behaviors of Walking Machines</article-title>. <source>Robot. Auton. Syst.</source> <volume>56</volume>, <fpage>256</fpage>&#x2013;<lpage>288</lpage>. <pub-id pub-id-type="doi">10.1016/j.robot.2007.07.004</pub-id> </citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Manoonpong</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Parlitz</surname>
<given-names>U.</given-names>
</name>
<name>
<surname>Worgotter</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Neural Control and Adaptive Neural Forward Models for Insect-like, Energy-Efficient, and Adaptable Locomotion of Walking Machines</article-title>. <source>Front. Neural Circuits</source> <volume>7</volume>, <fpage>12</fpage>. <pub-id pub-id-type="doi">10.3389/fncir.2013.00012</pub-id> </citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Marder</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Bucher</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Central Pattern Generators and the Control of Rhythmic Movements</article-title>. <source>Curr. Biol.</source> <volume>11</volume>, <fpage>R986</fpage>&#x2013;<lpage>R996</lpage>. <pub-id pub-id-type="doi">10.1016/s0960-9822(01)00581-4</pub-id> </citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Maufroy</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Kimura</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Takase</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Integration of Posture and Rhythmic Motion Controls in Quadrupedal Dynamic Walking Using Phase Modulations Based on Leg Loading/unloading</article-title>. <source>Auton. Robot</source> <volume>28</volume>, <fpage>331</fpage>&#x2013;<lpage>353</lpage>. <pub-id pub-id-type="doi">10.1007/s10514-009-9172-5</pub-id> </citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Moore</surname>
<given-names>J.&#x20;C.</given-names>
</name>
</person-group> (<year>1984</year>). <article-title>The Golgi Tendon Organ: a Review and Update</article-title>. <source>Am. J.&#x20;Occup. Ther.</source> <volume>38</volume>, <fpage>227</fpage>&#x2013;<lpage>236</lpage>. <pub-id pub-id-type="doi">10.5014/ajot.38.4.227</pub-id> </citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nachstedt</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Tetzlaff</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Manoonpong</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Fast Dynamical Coupling Enhances Frequency Adaptation of Oscillators for Robotic Locomotion Control</article-title>. <source>Front. Neurorobot.</source> <volume>11</volume>, <fpage>14</fpage>. <pub-id pub-id-type="doi">10.3389/fnbot.2017.00014</pub-id> </citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Owaki</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Ishiguro</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2017a</year>). <article-title>A Quadruped Robot Exhibiting Spontaneous Gait Transitions from Walking to Trotting to Galloping</article-title>. <source>Sci. Rep.</source> <volume>7</volume>, <fpage>291</fpage>&#x2013;<lpage>303</lpage>. <pub-id pub-id-type="doi">10.1038/s41598-017-00348-9</pub-id> </citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Owaki</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Ishiguro</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2017b</year>). <article-title>A Quadruped Robot Exhibiting Spontaneous Gait Transitions from Walking to Trotting to Galloping</article-title>. <source>Sci. Rep.</source> <volume>7</volume>, <fpage>277</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-017-00348-9</pub-id> </citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Owaki</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Kano</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Nagasawa</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Tero</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Simple Robot Suggests Physical Interlimb Communication Is Essential for Quadruped Walking</article-title>. <source>J.&#x20;R. Soc. Interf.</source> <volume>10</volume>, <fpage>20120669</fpage>. <pub-id pub-id-type="doi">10.1098/rsif.2012.0669</pub-id> </citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Owaki</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Goda</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Miyazawa</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ishiguro</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>A Minimal Model Describing Hexapedal Interlimb Coordination: The Tegotae-Based Approach</article-title>. <source>Front. Neurorobotics</source> <volume>11</volume>, <fpage>29</fpage>. <pub-id pub-id-type="doi">10.3389/fnbot.2017.00029</pub-id> </citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Owaki</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Horikiri</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Nishii</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ishiguro</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Tegotae-based Control Produces Adaptive Inter- and Intra-limb Coordination in Bipedal Walking</article-title>. <source>Front. Neurorobot.</source> <volume>15</volume>, <fpage>629595</fpage>. <pub-id pub-id-type="doi">10.3389/fnbot.2021.629595</pub-id> </citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pearson</surname>
<given-names>K. G.</given-names>
</name>
<name>
<surname>Iles</surname>
<given-names>J.&#x20;F.</given-names>
</name>
</person-group> (<year>1973</year>). <article-title>Nervous Mechanisms Underlying Intersegmental Co-ordination of Leg Movements during Walking in the Cockroach</article-title>. <source>J.&#x20;Exp. Biol.</source> <volume>58</volume>, <fpage>725</fpage>&#x2013;<lpage>744</lpage>. <pub-id pub-id-type="doi">10.1242/jeb.58.3.725</pub-id> </citation>
</ref>
<ref id="B51">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Pfeifer</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Bongard</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2006</year>). <source>How the Body Shapes the Way We Think: A New View of Intelligence</source>. <publisher-loc>Cambridge, MA</publisher-loc>: <publisher-name>The MIT Press</publisher-name>. <pub-id pub-id-type="doi">10.7551/mitpress/3585.001.0001</pub-id> </citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pfeifer</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Lungarella</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Iida</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Self-Organization, Embodiment, and Biologically Inspired Robotics</article-title>. <source>Science</source> <volume>318</volume>, <fpage>1088</fpage>&#x2013;<lpage>1093</lpage>. <pub-id pub-id-type="doi">10.1126/science.1145803</pub-id> </citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Righetti</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Ijspeert</surname>
<given-names>A. J.</given-names>
</name>
</person-group> (<year>2008</year>). &#x201c;<article-title>Pattern Generators with Sensory Feedback for the Control of Quadruped Locomotion</article-title>,&#x201d; in <conf-name>2008 IEEE International Conference on Robotics and Automation</conf-name>, <conf-loc>Pasadena, CA</conf-loc>, <conf-date>May 19&#x2013;23, 2008</conf-date>, <fpage>819</fpage>&#x2013;<lpage>824</lpage>. </citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schilling</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Hoinville</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Schmitz</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Cruse</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Walknet, a Bio-Inspired Controller for Hexapod Walking</article-title>. <source>Biol. Cybern.</source> <volume>107</volume>, <fpage>397</fpage>&#x2013;<lpage>419</lpage>. <pub-id pub-id-type="doi">10.1007/s00422-013-0563-5</pub-id> </citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schomburg</surname>
<given-names>E. D.</given-names>
</name>
<name>
<surname>Petersen</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Barajon</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Hultborn</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>Flexor Reflex Afferents Reset the Step Cycle during Fictive Locomotion in the Cat</article-title>. <source>Exp. Brain Res.</source> <volume>122</volume>, <fpage>339</fpage>&#x2013;<lpage>350</lpage>. <pub-id pub-id-type="doi">10.1007/s002210050522</pub-id> </citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shik</surname>
<given-names>M. L.</given-names>
</name>
<name>
<surname>Severin</surname>
<given-names>F. V.</given-names>
</name>
<name>
<surname>Orlovskii</surname>
<given-names>G. N.</given-names>
</name>
</person-group> (<year>1966</year>). <article-title>Control of Walking and Running by Means of Electrical Stimulation of the Mesencephalon</article-title>. <source>Electroencephalography Clin. Neurophysiol.</source> <volume>11</volume>, <fpage>659</fpage>&#x2013;<lpage>666</lpage>. </citation>
</ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Steingrube</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Timme</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>W&#xf6;rg&#xf6;tter</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Manoonpong</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Self-organized Adaptation of a Simple Neural Circuit Enables Complex Robot Behaviour</article-title>. <source>Nat. Phys.</source> <volume>6</volume>, <fpage>224</fpage>&#x2013;<lpage>230</lpage>. <pub-id pub-id-type="doi">10.1038/nphys1508</pub-id> </citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Taga</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Yamaguchi</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Shimizu</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>1991</year>). <article-title>Self-organized Control of Bipedal Locomotion by Neural Oscillators in Unpredictable Environment</article-title>. <source>Biol. Cybern.</source> <volume>65</volume>, <fpage>147</fpage>&#x2013;<lpage>159</lpage>. <pub-id pub-id-type="doi">10.1007/bf00198086</pub-id> </citation>
</ref>
<ref id="B59">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Taga</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>1994</year>). <article-title>Emergence of Bipedal Locomotion through Entrainment Among the Neuro-Musculo-Skeletal System and the Environment</article-title>. <source>Physica D: Nonlinear Phenomena</source> <volume>75</volume>, <fpage>190</fpage>&#x2013;<lpage>208</lpage>. <pub-id pub-id-type="doi">10.1016/0167-2789(94)90283-6</pub-id> </citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Taga</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>1995</year>). <article-title>A Model of the Neuro-Musculo-Skeletal System for Human Locomotion</article-title>. <source>Biol. Cybernetics</source> <volume>73</volume>, <fpage>97</fpage>&#x2013;<lpage>111</lpage>. <pub-id pub-id-type="doi">10.1007/s004220050166</pub-id> </citation>
</ref>
<ref id="B61">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tsujita</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Tsuchiya</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Onat</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Decentralized Autonomous Control of a Quadrupedal Locomotion Robot Using Oscillators</article-title>. <source>Artif. Life Robotics</source> <volume>5</volume>, <fpage>152</fpage>&#x2013;<lpage>158</lpage>. <pub-id pub-id-type="doi">10.1007/BF02481462</pub-id> </citation>
</ref>
<ref id="B62">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Tan</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2011</year>). &#x201c;<article-title>A Cpg-Based Sensory Feedback Control Method for Robotic Fish Locomotion</article-title>,&#x201d; in <conf-name>Proceedings of the 30th Chinese Control Conference (CCC) 2011</conf-name>, <conf-loc>Yantai, China</conf-loc>, <conf-date>July 22&#x2013;24, 2011</conf-date> </citation>
</ref>
<ref id="B63">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xiong</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>W&#xf6;rg&#xf6;tter</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Manoonpong</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Neuromechanical Control for Hexapedal Robot Walking on Challenging Surfaces and Surface Classification</article-title>. <source>Robotics Autonomous Syst.</source> <volume>62</volume>, <fpage>1777</fpage>&#x2013;<lpage>1789</lpage>. <pub-id pub-id-type="doi">10.1016/j.robot.2014.07.008</pub-id> </citation>
</ref>
<ref id="B64">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xiong</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Worgotter</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Manoonpong</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Adaptive and Energy Efficient Walking in a Hexapod Robot under Neuromechanical Control and Sensorimotor Learning</article-title>. <source>IEEE Trans. Cybern.</source> <volume>46</volume>, <fpage>2521</fpage>&#x2013;<lpage>2534</lpage>. <pub-id pub-id-type="doi">10.1109/TCYB.2015.2479237</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>