Identification of the contribution of the ankle and hip joints to multi-segmental balance control
Identifieur interne : 001260 ( Pmc/Curation ); précédent : 001259; suivant : 001261Identification of the contribution of the ankle and hip joints to multi-segmental balance control
Auteurs : Tjitske Anke Boonstra [Pays-Bas] ; Alfred C. Schouten [Pays-Bas] ; Herman Van Der Kooij [Pays-Bas]Source :
- Journal of NeuroEngineering and Rehabilitation [ 1743-0003 ] ; 2013.
Abstract
Human stance involves multiple segments, including the legs and trunk, and requires coordinated actions of both. A novel method was developed that reliably estimates the contribution of the left and right leg (i.e., the ankle and hip joints) to the balance control of individual subjects.
The method was evaluated using simulations of a double-inverted pendulum model and the applicability was demonstrated with an experiment with seven healthy and one Parkinsonian participant. Model simulations indicated that two perturbations are required to reliably estimate the dynamics of a double-inverted pendulum balance control system. In the experiment, two multisine perturbation signals were applied simultaneously. The balance control system dynamic behaviour of the participants was estimated by Frequency Response Functions (FRFs), which relate ankle and hip joint angles to joint torques, using a multivariate closed-loop system identification technique.
In the model simulations, the FRFs were reliably estimated, also in the presence of realistic levels of noise. In the experiment, the participants responded consistently to the perturbations, indicated by low noise-to-signal ratios of the ankle angle (0.24), hip angle (0.28), ankle torque (0.07), and hip torque (0.33). The developed method could detect that the Parkinson patient controlled his balance asymmetrically, that is, the right ankle and hip joints produced more corrective torque.
The method allows for a reliable estimate of the multisegmental feedback mechanism that stabilizes stance, of individual participants and of separate legs.
Url:
DOI: 10.1186/1743-0003-10-23
PubMed: 23433148
PubMed Central: 3662596
Links toward previous steps (curation, corpus...)
- to stream Pmc, to step Corpus: Pour aller vers cette notice dans l'étape Curation :001260
Links to Exploration step
PMC:3662596Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Identification of the contribution of the ankle and hip joints to multi-segmental balance control</title>
<author><name sortKey="Boonstra, Tjitske Anke" sort="Boonstra, Tjitske Anke" uniqKey="Boonstra T" first="Tjitske Anke" last="Boonstra">Tjitske Anke Boonstra</name>
<affiliation wicri:level="1"><nlm:aff id="I1">Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500, The Netherlands</nlm:aff>
<country xml:lang="fr">Pays-Bas</country>
<wicri:regionArea>Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Schouten, Alfred C" sort="Schouten, Alfred C" uniqKey="Schouten A" first="Alfred C" last="Schouten">Alfred C. Schouten</name>
<affiliation wicri:level="1"><nlm:aff id="I1">Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500, The Netherlands</nlm:aff>
<country xml:lang="fr">Pays-Bas</country>
<wicri:regionArea>Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500</wicri:regionArea>
</affiliation>
<affiliation wicri:level="1"><nlm:aff id="I2">Department of Biomechanical Engineering, Delft University of Technology, Mekelweg 2, Delft, CD, 2628, The Netherlands</nlm:aff>
<country xml:lang="fr">Pays-Bas</country>
<wicri:regionArea>Department of Biomechanical Engineering, Delft University of Technology, Mekelweg 2, Delft, CD, 2628</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Van Der Kooij, Herman" sort="Van Der Kooij, Herman" uniqKey="Van Der Kooij H" first="Herman" last="Van Der Kooij">Herman Van Der Kooij</name>
<affiliation wicri:level="1"><nlm:aff id="I1">Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500, The Netherlands</nlm:aff>
<country xml:lang="fr">Pays-Bas</country>
<wicri:regionArea>Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500</wicri:regionArea>
</affiliation>
<affiliation wicri:level="1"><nlm:aff id="I2">Department of Biomechanical Engineering, Delft University of Technology, Mekelweg 2, Delft, CD, 2628, The Netherlands</nlm:aff>
<country xml:lang="fr">Pays-Bas</country>
<wicri:regionArea>Department of Biomechanical Engineering, Delft University of Technology, Mekelweg 2, Delft, CD, 2628</wicri:regionArea>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">23433148</idno>
<idno type="pmc">3662596</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3662596</idno>
<idno type="RBID">PMC:3662596</idno>
<idno type="doi">10.1186/1743-0003-10-23</idno>
<date when="2013">2013</date>
<idno type="wicri:Area/Pmc/Corpus">001260</idno>
<idno type="wicri:Area/Pmc/Curation">001260</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Identification of the contribution of the ankle and hip joints to multi-segmental balance control</title>
<author><name sortKey="Boonstra, Tjitske Anke" sort="Boonstra, Tjitske Anke" uniqKey="Boonstra T" first="Tjitske Anke" last="Boonstra">Tjitske Anke Boonstra</name>
<affiliation wicri:level="1"><nlm:aff id="I1">Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500, The Netherlands</nlm:aff>
<country xml:lang="fr">Pays-Bas</country>
<wicri:regionArea>Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Schouten, Alfred C" sort="Schouten, Alfred C" uniqKey="Schouten A" first="Alfred C" last="Schouten">Alfred C. Schouten</name>
<affiliation wicri:level="1"><nlm:aff id="I1">Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500, The Netherlands</nlm:aff>
<country xml:lang="fr">Pays-Bas</country>
<wicri:regionArea>Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500</wicri:regionArea>
</affiliation>
<affiliation wicri:level="1"><nlm:aff id="I2">Department of Biomechanical Engineering, Delft University of Technology, Mekelweg 2, Delft, CD, 2628, The Netherlands</nlm:aff>
<country xml:lang="fr">Pays-Bas</country>
<wicri:regionArea>Department of Biomechanical Engineering, Delft University of Technology, Mekelweg 2, Delft, CD, 2628</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Van Der Kooij, Herman" sort="Van Der Kooij, Herman" uniqKey="Van Der Kooij H" first="Herman" last="Van Der Kooij">Herman Van Der Kooij</name>
<affiliation wicri:level="1"><nlm:aff id="I1">Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500, The Netherlands</nlm:aff>
<country xml:lang="fr">Pays-Bas</country>
<wicri:regionArea>Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500</wicri:regionArea>
</affiliation>
<affiliation wicri:level="1"><nlm:aff id="I2">Department of Biomechanical Engineering, Delft University of Technology, Mekelweg 2, Delft, CD, 2628, The Netherlands</nlm:aff>
<country xml:lang="fr">Pays-Bas</country>
<wicri:regionArea>Department of Biomechanical Engineering, Delft University of Technology, Mekelweg 2, Delft, CD, 2628</wicri:regionArea>
</affiliation>
</author>
</analytic>
<series><title level="j">Journal of NeuroEngineering and Rehabilitation</title>
<idno type="eISSN">1743-0003</idno>
<imprint><date when="2013">2013</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><sec><title>Background</title>
<p>Human stance involves multiple segments, including the legs and trunk, and requires coordinated actions of both. A novel method was developed that reliably estimates the contribution of the left and right leg (i.e., the ankle and hip joints) to the balance control of individual subjects.</p>
</sec>
<sec><title>Methods</title>
<p>The method was evaluated using simulations of a double-inverted pendulum model and the applicability was demonstrated with an experiment with seven healthy and one Parkinsonian participant. Model simulations indicated that two perturbations are required to reliably estimate the dynamics of a double-inverted pendulum balance control system. In the experiment, two multisine perturbation signals were applied simultaneously. The balance control system dynamic behaviour of the participants was estimated by Frequency Response Functions (FRFs), which relate ankle and hip joint angles to joint torques, using a multivariate closed-loop system identification technique.</p>
</sec>
<sec><title>Results</title>
<p>In the model simulations, the FRFs were reliably estimated, also in the presence of realistic levels of noise. In the experiment, the participants responded consistently to the perturbations, indicated by low noise-to-signal ratios of the ankle angle (0.24), hip angle (0.28), ankle torque (0.07), and hip torque (0.33). The developed method could detect that the Parkinson patient controlled his balance asymmetrically, that is, the right ankle and hip joints produced more corrective torque.</p>
</sec>
<sec><title>Conclusion</title>
<p>The method allows for a reliable estimate of the multisegmental feedback mechanism that stabilizes stance, of individual participants and of separate legs.</p>
</sec>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Peterka, Rj" uniqKey="Peterka R">RJ Peterka</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kiemel, T" uniqKey="Kiemel T">T Kiemel</name>
</author>
<author><name sortKey="Zhang, Y" uniqKey="Zhang Y">Y Zhang</name>
</author>
<author><name sortKey="Jeka, Jj" uniqKey="Jeka J">JJ Jeka</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Peterka, Rj" uniqKey="Peterka R">RJ Peterka</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Van Asseldonk, Eh" uniqKey="Van Asseldonk E">EH van Asseldonk</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Van Soest, Aj" uniqKey="Van Soest A">AJ van Soest</name>
</author>
<author><name sortKey="Rozendaal, La" uniqKey="Rozendaal L">LA Rozendaal</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Horlings, Cg" uniqKey="Horlings C">CG Horlings</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Creath, R" uniqKey="Creath R">R Creath</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pinter, Ij" uniqKey="Pinter I">IJ Pinter</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Horak, Fb" uniqKey="Horak F">FB Horak</name>
</author>
<author><name sortKey="Nashner, Lm" uniqKey="Nashner L">LM Nashner</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Horak, Fb" uniqKey="Horak F">FB Horak</name>
</author>
<author><name sortKey="Macpherson, Jm" uniqKey="Macpherson J">JM Macpherson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fitzpatrick, R" uniqKey="Fitzpatrick R">R Fitzpatrick</name>
</author>
<author><name sortKey="Burke, D" uniqKey="Burke D">D Burke</name>
</author>
<author><name sortKey="Gandevia, Sc" uniqKey="Gandevia S">SC Gandevia</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Van Der Kooij, H" uniqKey="Van Der Kooij H">H van der Kooij</name>
</author>
<author><name sortKey="Van Asseldonk, E" uniqKey="Van Asseldonk E">E van Asseldonk</name>
</author>
<author><name sortKey="Van Der Helm, Fc" uniqKey="Van Der Helm F">FC van der Helm</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pintelon, R" uniqKey="Pintelon R">R Pintelon</name>
</author>
<author><name sortKey="Schoukens, J" uniqKey="Schoukens J">J Schoukens</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="De Vlugt, E" uniqKey="De Vlugt E">E de Vlugt</name>
</author>
<author><name sortKey="Schouten, Ac" uniqKey="Schouten A">AC Schouten</name>
</author>
<author><name sortKey="Van Der Helm, Fc" uniqKey="Van Der Helm F">FC van der Helm</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kuo, Ad" uniqKey="Kuo A">AD Kuo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Speers, Ra" uniqKey="Speers R">RA Speers</name>
</author>
<author><name sortKey="Kuo, Ad" uniqKey="Kuo A">AD Kuo</name>
</author>
<author><name sortKey="Horak, Fb" uniqKey="Horak F">FB Horak</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Alexandrov, Av" uniqKey="Alexandrov A">AV Alexandrov</name>
</author>
<author><name sortKey="Frolov, Aa" uniqKey="Frolov A">AA Frolov</name>
</author>
<author><name sortKey="Massion, J" uniqKey="Massion J">J Massion</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Alexandrov, Av" uniqKey="Alexandrov A">AV Alexandrov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Park, S" uniqKey="Park S">S Park</name>
</author>
<author><name sortKey="Horak, Fb" uniqKey="Horak F">FB Horak</name>
</author>
<author><name sortKey="Kuo, Ad" uniqKey="Kuo A">AD Kuo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kim, S" uniqKey="Kim S">S Kim</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fujisawa, N" uniqKey="Fujisawa N">N Fujisawa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stolze, H" uniqKey="Stolze H">H Stolze</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pickering, Rm" uniqKey="Pickering R">RM Pickering</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hely, Ma" uniqKey="Hely M">MA Hely</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Grimbergen, Ya" uniqKey="Grimbergen Y">YA Grimbergen</name>
</author>
<author><name sortKey="Munneke, M" uniqKey="Munneke M">M Munneke</name>
</author>
<author><name sortKey="Bloem, Br" uniqKey="Bloem B">BR Bloem</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Boonstra, Ta" uniqKey="Boonstra T">TA Boonstra</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Colnat Coulbois, S" uniqKey="Colnat Coulbois S">S Colnat-Coulbois</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Carpenter, Mg" uniqKey="Carpenter M">MG Carpenter</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Termoz, N" uniqKey="Termoz N">N Termoz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Maurer, C" uniqKey="Maurer C">C Maurer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Djaldetti, R" uniqKey="Djaldetti R">R Djaldetti</name>
</author>
<author><name sortKey="Ziv, I" uniqKey="Ziv I">I Ziv</name>
</author>
<author><name sortKey="Melamed, E" uniqKey="Melamed E">E Melamed</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Geurts, Ac" uniqKey="Geurts A">AC Geurts</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Van Der Kooij, H" uniqKey="Van Der Kooij H">H van der Kooij</name>
</author>
<author><name sortKey="Peterka, Rj" uniqKey="Peterka R">RJ Peterka</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Robinson, D" uniqKey="Robinson D">D Robinson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Van Der Kooij, H" uniqKey="Van Der Kooij H">H van der Kooij</name>
</author>
<author><name sortKey="De Vlugt, E" uniqKey="De Vlugt E">E de Vlugt</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zajac, Fe" uniqKey="Zajac F">FE Zajac</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nott, Cr" uniqKey="Nott C">CR Nott</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pasma, Jh" uniqKey="Pasma J">JH Pasma</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Koopman, Hfjm" uniqKey="Koopman H">HFJM Koopman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Koopman, B" uniqKey="Koopman B">B Koopman</name>
</author>
<author><name sortKey="Grootenboer, Hj" uniqKey="Grootenboer H">HJ Grootenboer</name>
</author>
<author><name sortKey="De Jongh, Hj" uniqKey="De Jongh H">HJ de Jongh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Brand, Ra" uniqKey="Brand R">RA Brand</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chandler, Rf" uniqKey="Chandler R">RF Chandler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Preuss, R" uniqKey="Preuss R">R Preuss</name>
</author>
<author><name sortKey="Fung, J" uniqKey="Fung J">J Fung</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jenkins, G" uniqKey="Jenkins G">G Jenkins</name>
</author>
<author><name sortKey="Watts, D" uniqKey="Watts D">D Watts</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kim, S" uniqKey="Kim S">S Kim</name>
</author>
<author><name sortKey="Atkeson, Cg" uniqKey="Atkeson C">CG Atkeson</name>
</author>
<author><name sortKey="Park, S" uniqKey="Park S">S Park</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="De Vlugt, E" uniqKey="De Vlugt E">E de Vlugt</name>
</author>
<author><name sortKey="Schouten, Ac" uniqKey="Schouten A">AC Schouten</name>
</author>
<author><name sortKey="Van Der Helm, Fc" uniqKey="Van Der Helm F">FC van der Helm</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kiemel, T" uniqKey="Kiemel T">T Kiemel</name>
</author>
<author><name sortKey="Elahi, Aj" uniqKey="Elahi A">AJ Elahi</name>
</author>
<author><name sortKey="Jeka, Jj" uniqKey="Jeka J">JJ Jeka</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kiemel, T" uniqKey="Kiemel T">T Kiemel</name>
</author>
<author><name sortKey="Zhang, Y" uniqKey="Zhang Y">Y Zhang</name>
</author>
<author><name sortKey="Jeka, Jj" uniqKey="Jeka J">JJ Jeka</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rocchi, L" uniqKey="Rocchi L">L Rocchi</name>
</author>
<author><name sortKey="Chiari, L" uniqKey="Chiari L">L Chiari</name>
</author>
<author><name sortKey="Horak, Fb" uniqKey="Horak F">FB Horak</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Van Der Kooij, H" uniqKey="Van Der Kooij H">H van der Kooij</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Roerdink, M" uniqKey="Roerdink M">M Roerdink</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Winter, Da" uniqKey="Winter D">DA Winter</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article" xml:lang="en"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">J Neuroeng Rehabil</journal-id>
<journal-id journal-id-type="iso-abbrev">J Neuroeng Rehabil</journal-id>
<journal-title-group><journal-title>Journal of NeuroEngineering and Rehabilitation</journal-title>
</journal-title-group>
<issn pub-type="epub">1743-0003</issn>
<publisher><publisher-name>BioMed Central</publisher-name>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">23433148</article-id>
<article-id pub-id-type="pmc">3662596</article-id>
<article-id pub-id-type="publisher-id">1743-0003-10-23</article-id>
<article-id pub-id-type="doi">10.1186/1743-0003-10-23</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Methodology</subject>
</subj-group>
</article-categories>
<title-group><article-title>Identification of the contribution of the ankle and hip joints to multi-segmental balance control</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" corresp="yes" id="A1"><name><surname>Boonstra</surname>
<given-names>Tjitske Anke</given-names>
</name>
<xref ref-type="aff" rid="I1">1</xref>
<email>t.a.boonstra@gmail.com</email>
</contrib>
<contrib contrib-type="author" id="A2"><name><surname>Schouten</surname>
<given-names>Alfred C</given-names>
</name>
<xref ref-type="aff" rid="I1">1</xref>
<xref ref-type="aff" rid="I2">2</xref>
<email>a.c.schouten@tudelft.nl</email>
</contrib>
<contrib contrib-type="author" id="A3"><name><surname>van der Kooij</surname>
<given-names>Herman</given-names>
</name>
<xref ref-type="aff" rid="I1">1</xref>
<xref ref-type="aff" rid="I2">2</xref>
<email>h.vanderkooij@utwente.nl</email>
</contrib>
</contrib-group>
<aff id="I1"><label>1</label>
Laboratory for Biomechanical Engineering, MIRA institute for biomechanical technology and technical medicine, University of Twente, Faculty of Engineering Technology, PO Box 217, Enschede, AE 7500, The Netherlands</aff>
<aff id="I2"><label>2</label>
Department of Biomechanical Engineering, Delft University of Technology, Mekelweg 2, Delft, CD, 2628, The Netherlands</aff>
<pub-date pub-type="collection"><year>2013</year>
</pub-date>
<pub-date pub-type="epub"><day>22</day>
<month>2</month>
<year>2013</year>
</pub-date>
<volume>10</volume>
<fpage>23</fpage>
<lpage>23</lpage>
<history><date date-type="received"><day>23</day>
<month>1</month>
<year>2012</year>
</date>
<date date-type="accepted"><day>24</day>
<month>1</month>
<year>2013</year>
</date>
</history>
<permissions><copyright-statement>Copyright © 2013 Boonstra et al.; licensee BioMed Central Ltd.</copyright-statement>
<copyright-year>2013</copyright-year>
<copyright-holder>Boonstra et al.; licensee BioMed Central Ltd.</copyright-holder>
<license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/2.0"><license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0">http://creativecommons.org/licenses/by/2.0</ext-link>
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<self-uri xlink:href="http://www.jneuroengrehab.com/content/10/1/23"></self-uri>
<abstract><sec><title>Background</title>
<p>Human stance involves multiple segments, including the legs and trunk, and requires coordinated actions of both. A novel method was developed that reliably estimates the contribution of the left and right leg (i.e., the ankle and hip joints) to the balance control of individual subjects.</p>
</sec>
<sec><title>Methods</title>
<p>The method was evaluated using simulations of a double-inverted pendulum model and the applicability was demonstrated with an experiment with seven healthy and one Parkinsonian participant. Model simulations indicated that two perturbations are required to reliably estimate the dynamics of a double-inverted pendulum balance control system. In the experiment, two multisine perturbation signals were applied simultaneously. The balance control system dynamic behaviour of the participants was estimated by Frequency Response Functions (FRFs), which relate ankle and hip joint angles to joint torques, using a multivariate closed-loop system identification technique.</p>
</sec>
<sec><title>Results</title>
<p>In the model simulations, the FRFs were reliably estimated, also in the presence of realistic levels of noise. In the experiment, the participants responded consistently to the perturbations, indicated by low noise-to-signal ratios of the ankle angle (0.24), hip angle (0.28), ankle torque (0.07), and hip torque (0.33). The developed method could detect that the Parkinson patient controlled his balance asymmetrically, that is, the right ankle and hip joints produced more corrective torque.</p>
</sec>
<sec><title>Conclusion</title>
<p>The method allows for a reliable estimate of the multisegmental feedback mechanism that stabilizes stance, of individual participants and of separate legs.</p>
</sec>
</abstract>
<kwd-group><kwd>Balance control</kwd>
<kwd>Closed-loop system identification</kwd>
<kwd>Multivariate systems</kwd>
<kwd>Asymmetry</kwd>
</kwd-group>
</article-meta>
</front>
<body><sec><title>Background</title>
<p>Maintaining an upright posture is a relatively easy task for healthy humans [<xref ref-type="bibr" rid="B1">1</xref>
,<xref ref-type="bibr" rid="B2">2</xref>
]. However, bipedal upright stance is inherently unstable, as small deviations from the upright posture result in disturbing torques due to gravity, which drives the system further away from upright posture [<xref ref-type="bibr" rid="B3">3</xref>
]. To stay upright, the body generates corrective torques to counteract the effects of internal (e.g., motor and sensory noise) and external (e.g., uneven surfaces) perturbations.</p>
<p>When postural deviations are small, the body is often simplified as an inverted pendulum pivoting at the ankles, which describes the so-called ankle strategy [<xref ref-type="bibr" rid="B3">3</xref>
-<xref ref-type="bibr" rid="B5">5</xref>
]. However, several studies demonstrated that human movement during stance is multi-segmental [<xref ref-type="bibr" rid="B6">6</xref>
-<xref ref-type="bibr" rid="B8">8</xref>
] and for example, the hips substantially contribute to upright stance (i.e., the hip strategy [<xref ref-type="bibr" rid="B9">9</xref>
]). Human balance control is a closed-loop multi-segmental process, i.e., sensory signals about the movement of the body are fed back to the central nervous system (CNS), and the CNS controls the muscles to generate adequate responses [<xref ref-type="bibr" rid="B10">10</xref>
,<xref ref-type="bibr" rid="B11">11</xref>
]. In a noisy closed-loop system, like human balance control, causality is difficult to determine and the dynamics of the different components (i.e., the body and the stabilizing mechanisms located in the CNS) affect both the input (joint angles) and output signals (joint torques). To “open” the loop and to separate the dynamics of the different components, the balance system need to be perturbed [<xref ref-type="bibr" rid="B11">11</xref>
,<xref ref-type="bibr" rid="B12">12</xref>
]. Furthermore, when estimating the dynamics in a <bold>multivariate</bold>
system, <bold>multiple</bold>
perturbations need to be applied [<xref ref-type="bibr" rid="B13">13</xref>
,<xref ref-type="bibr" rid="B14">14</xref>
].</p>
<p>Most studies investigating the multivariate nature of balance control do not take the multivariate noisy closed-loop nature into account, by either not using perturbations [<xref ref-type="bibr" rid="B15">15</xref>
-<xref ref-type="bibr" rid="B17">17</xref>
], or by using only one perturbation [<xref ref-type="bibr" rid="B18">18</xref>
-<xref ref-type="bibr" rid="B20">20</xref>
]. Only two studies investigated the multivariate nature of balance control by applying two perturbations [<xref ref-type="bibr" rid="B2">2</xref>
,<xref ref-type="bibr" rid="B21">21</xref>
].</p>
<p>Fujisawa and colleagues [<xref ref-type="bibr" rid="B21">21</xref>
] investigated the role of the hip joint to upright stance by applying pseudorandom perturbations (bandwidth 0 - 0.83 Hz) while manipulating the support surface width. Subsequently, an ARMAX model (with joint angles as input and joint torques as outputs) was fitted to the data to obtain the Frequency Response Functions (FRFs) of a two-segment model of balance control. Results showed an increase of balance contribution of the hip joint, when the support surface became narrower.</p>
<p>Jeka and colleagues [<xref ref-type="bibr" rid="B2">2</xref>
] identified neural feedback during upright stance in 18 subjects, while applying two mechanical perturbations (springs attached to a linear motor) and one sensory perturbation (visual scene rotations). By comparing the identified neural feedback (from joint angles to weighted electromyograms (EMGs) of the leg and trunk segments) with a large range of cost functions, it was concluded that the CNS stabilizes stance with near minimum muscle activation.</p>
<p>Ageing and many neurological diseases are associated with balance impairments and falls [<xref ref-type="bibr" rid="B22">22</xref>
]. Understanding the (patho)physiology of upright stance could aid to detect individuals with an increased risk of falls, help to design and evaluate intervention programs or monitor disease progression. Therefore, for clinical applications, it is very important to obtain a reliable individual estimate of balance control.</p>
<p>Of all neurological diseases, PD patients are at the highest risk of falling [<xref ref-type="bibr" rid="B22">22</xref>
-<xref ref-type="bibr" rid="B24">24</xref>
], but the pathophysiology of balance impairments in PD remains unclear [<xref ref-type="bibr" rid="B25">25</xref>
,<xref ref-type="bibr" rid="B26">26</xref>
]. Recently, it was suggested that one of the factors contributing to decreased balance control in PD patients, is impaired trunk control [<xref ref-type="bibr" rid="B27">27</xref>
,<xref ref-type="bibr" rid="B28">28</xref>
] or a decreased intersegmental coordination [<xref ref-type="bibr" rid="B29">29</xref>
,<xref ref-type="bibr" rid="B30">30</xref>
]. Another factor could be asymmetrical balance control, that is when one leg produces more force than the other leg to maintain an upright posture. Asymmetries in balance control have been rarely studied in PD, although it is an asymmetrical disease [<xref ref-type="bibr" rid="B31">31</xref>
]. One study [<xref ref-type="bibr" rid="B32">32</xref>
] found balance control asymmetries in 24% of the PD participants, indicating that balance asymmetries are important in PD.</p>
<p>Currently, there is no method available that can identify a multisegmental stabilizing mechanism of balance control on an individual level, separating the contribution of the joints of the left and right body side. We developed and evaluated a non-parametric MIMO (Multiple-Input-Multiple-Output) identification method based on the previously used non-parametric system SISO (Single-Input-Single-Output) identification method [<xref ref-type="bibr" rid="B4">4</xref>
]. To obtain a reliable individual response, periodic perturbations were applied, which have the advantage of having power at specific frequencies, decreasing the measurement time, and increasing the participants response. In addition, the stabilizing mechanisms were estimated based on left and right joint torques (contrary to weighted EMGs; [<xref ref-type="bibr" rid="B2">2</xref>
]), to be able to investigate balance control asymmetries.</p>
<p>In sum, our goal was to develop a method that can reliably estimate the stabilizing mechanisms of the closed-loop multivariate balance control system of <italic>individual</italic>
participants, making a distinction between the contribution of the left and right leg. The (clinical) applicability is demonstrated in an experiment that perturbed the balance of seven healthy participants and a PD patient with a novel device that can apply two independent mechanical perturbations.</p>
</sec>
<sec sec-type="methods"><title>Methods</title>
<sec><title>Model simulations</title>
<p>A two degree-of-freedom (DoF) balance control model (described extensively in Appendix A) was implemented in Matlab (The Mathworks, Natick, USA) and simulated with Simulink (equations were solved with a fifth order Dormand-Prince algorithm). The human balance model consisted of a two-segment human body with two actuators (ankle and hip), which were controlled using feedback of the joint angles (ankle and hip). In the model, no distinction between the left and right leg was made.</p>
<p>We perturbed the model with one and two perturbations. Two possible perturbation configurations of the two perturbations were evaluated: 1) external forces at the hip and shoulder (see Figure <xref ref-type="fig" rid="F1">1</xref>
), similar to push-pull rods and 2) a combination of platform forward-backward platform translations and a perturbation torque around the ankle (see Figure <xref ref-type="fig" rid="F2">2</xref>
). Also, simulations without and with pink sensor noise [<xref ref-type="bibr" rid="B33">33</xref>
] and white (measurement) noise were evaluated:</p>
<fig id="F1" position="float"><label>Figure 1</label>
<caption><p><bold>Multiple-Input-Multiple-Output closed-loop balance control system. </bold>
The body mechanics represent the dynamics of a double-inverted pendulum with the corrective ankle and hip torques as inputs and the joint angles as outputs. The stabilizing mechanisms represent the dynamics of the combination of active and passive feedback pathways of the concerned body(part) and generates a torque to correct for the deviation of upright stance. The balance control model can be perturbed with support surface movements (S<sub>x</sub>
), perturbation forces at the hip (Fp<sub>1</sub>
) or at the shoulder (Fp<sub>2</sub>
). Positive torques and positive angles are defined as counterclockwise.</p>
</caption>
<graphic xlink:href="1743-0003-10-23-1"></graphic>
</fig>
<fig id="F2" position="float"><label>Figure 2</label>
<caption><p><bold>Experimental set-up. </bold>
The participant stands on the dual forceplate (<bold>A</bold>
) embedded in the movement platform (<bold>B</bold>
). Two independent perturbations are applied with the movement platform (<bold>B</bold>
) and the pusher (<bold>C</bold>
) in the forward-backward direction. Interaction forces between the pusher (<bold>C</bold>
) and the participant are measured with a force sensor (<bold>D</bold>
). Actual falls are prevented by the safety harness (<bold>E</bold>
). Reflective spherical markers (<bold>F</bold>
) measure the movements of the participant.</p>
</caption>
<graphic xlink:href="1743-0003-10-23-2"></graphic>
</fig>
<p>i. two perturbation forces with noise (2F-N),</p>
<p>ii. two perturbation forces without noise (2F-Nn),</p>
<p>iii. platform perturbation and ankle torque with noise (PLT-N),</p>
<p>iv. platform perturbation and ankle torque without noise (PLT-Nn),</p>
<p>v. one perturbation and one perturbation round with noise (PL-N),</p>
<p>vi. one perturbation and one perturbation round without noise.</p>
<p>The characteristics of the perturbation signals are described in detail in the ‘Disturbance signals’ section and the amplitudes and power spectra are reported in Table <xref ref-type="table" rid="T1">1</xref>
. The input and output signals of the model were sampled at 120 Hz.</p>
<table-wrap position="float" id="T1"><label>Table 1</label>
<caption><p>Perturbation amplitudes, scaling of the perturbation signal, Goodness of Fit and Noise-to-Signal Ratios of the model simulations and the balance control experiment</p>
</caption>
<table frame="hsides" rules="groups" border="1"><colgroup><col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
</colgroup>
<thead valign="top"><tr><th align="center" valign="bottom"> <hr></hr>
</th>
<th colspan="6" align="center" valign="bottom"><bold>Simulations</bold>
<hr></hr>
</th>
<th colspan="3" align="center" valign="bottom"><bold>Experiment</bold>
<hr></hr>
</th>
</tr>
<tr><th align="center"> </th>
<th align="center"><bold>PLT-Nn</bold>
</th>
<th align="center"><bold>PLT-N</bold>
</th>
<th align="center"><bold>2F-Nn</bold>
</th>
<th align="center"><bold>2F-N</bold>
</th>
<th align="center"><bold>PL-Nn</bold>
</th>
<th align="center"><bold>PL-N</bold>
</th>
<th align="center"><bold>HC</bold>
</th>
<th align="center"><bold>HC mean</bold>
</th>
<th align="center"><bold>PD patient</bold>
</th>
</tr>
</thead>
<tbody valign="top"><tr><td align="center" valign="bottom">Signal Amplitude<hr></hr>
</td>
<td align="center" valign="bottom">PL: 0.06; T: 20Nm<hr></hr>
</td>
<td align="center" valign="bottom">PL: 0.06m; T: 16Nm<hr></hr>
</td>
<td align="center" valign="bottom">20N<hr></hr>
</td>
<td align="center" valign="bottom">20N<hr></hr>
</td>
<td align="center" valign="bottom">0.12m<hr></hr>
</td>
<td align="center" valign="bottom">0.12m<hr></hr>
</td>
<td align="center" valign="bottom">PL: 0.07; T: 18.4Nm<hr></hr>
</td>
<td align="center" valign="bottom">PL: 0.06m; T: 16Nm<hr></hr>
</td>
<td align="center" valign="bottom">PL: 0.06; T: 20Nm<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">Scaling<hr></hr>
</td>
<td align="center" valign="bottom">PL: 1/<italic>f</italic>
<hr></hr>
</td>
<td align="center" valign="bottom">PL: 1/<italic>f</italic>
<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
<td align="center" valign="bottom">1/<italic>f</italic>
<hr></hr>
</td>
<td align="center" valign="bottom">1/<italic>f</italic>
<hr></hr>
</td>
<td align="center" valign="bottom">PL: 1/<italic>f</italic>
<hr></hr>
</td>
<td align="center" valign="bottom">PL: 1/<italic>f</italic>
<hr></hr>
</td>
<td align="center" valign="bottom">PL: 1/<italic>f</italic>
<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom"><bold>GOF (1/Hz</bold>
<sup><bold>2</bold>
</sup>
<bold>)</bold>
<hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom"><inline-formula><mml:math id="M1" name="1743-0003-10-23-i2" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
<hr></hr>
</td>
<td align="center" valign="bottom">1.11<hr></hr>
</td>
<td align="center" valign="bottom">11.64<hr></hr>
</td>
<td align="center" valign="bottom">1.11<hr></hr>
</td>
<td align="center" valign="bottom">23.64<hr></hr>
</td>
<td align="center" valign="bottom">26.44<hr></hr>
</td>
<td align="center" valign="bottom">29.83<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom"><inline-formula><mml:math id="M2" name="1743-0003-10-23-i3" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
<hr></hr>
</td>
<td align="center" valign="bottom">0.25<hr></hr>
</td>
<td align="center" valign="bottom">4.99<hr></hr>
</td>
<td align="center" valign="bottom">0.25<hr></hr>
</td>
<td align="center" valign="bottom">7.16<hr></hr>
</td>
<td align="center" valign="bottom">21.56<hr></hr>
</td>
<td align="center" valign="bottom">22.59<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom"><inline-formula><mml:math id="M3" name="1743-0003-10-23-i4" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
<hr></hr>
</td>
<td align="center" valign="bottom">0.79<hr></hr>
</td>
<td align="center" valign="bottom">7.30<hr></hr>
</td>
<td align="center" valign="bottom">0.79<hr></hr>
</td>
<td align="center" valign="bottom">14.42<hr></hr>
</td>
<td align="center" valign="bottom">28.18<hr></hr>
</td>
<td align="center" valign="bottom">36.09<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom"><inline-formula><mml:math id="M4" name="1743-0003-10-23-i5" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
<hr></hr>
</td>
<td align="center" valign="bottom">0.25<hr></hr>
</td>
<td align="center" valign="bottom">2.94<hr></hr>
</td>
<td align="center" valign="bottom">0.25<hr></hr>
</td>
<td align="center" valign="bottom">4.38<hr></hr>
</td>
<td align="center" valign="bottom">12.51<hr></hr>
</td>
<td align="center" valign="bottom">14.55<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
<td align="center" valign="bottom">-<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom"><bold>NSR</bold>
<hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">T<sub>A</sub>
<hr></hr>
</td>
<td align="center" valign="bottom">0<hr></hr>
</td>
<td align="center" valign="bottom">0.54<hr></hr>
</td>
<td align="center" valign="bottom">0<hr></hr>
</td>
<td align="center" valign="bottom">0.43<hr></hr>
</td>
<td align="center" valign="bottom">0<hr></hr>
</td>
<td align="center" valign="bottom">1.17<hr></hr>
</td>
<td align="center" valign="bottom">0.55<hr></hr>
</td>
<td align="center" valign="bottom">0.24<hr></hr>
</td>
<td align="center" valign="bottom">0.11<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">T<sub>H</sub>
<hr></hr>
</td>
<td align="center" valign="bottom">0<hr></hr>
</td>
<td align="center" valign="bottom">0.39<hr></hr>
</td>
<td align="center" valign="bottom">0<hr></hr>
</td>
<td align="center" valign="bottom">0.10<hr></hr>
</td>
<td align="center" valign="bottom">0<hr></hr>
</td>
<td align="center" valign="bottom">0.59<hr></hr>
</td>
<td align="center" valign="bottom">0.08<hr></hr>
</td>
<td align="center" valign="bottom">0.28<hr></hr>
</td>
<td align="center" valign="bottom">0.06<hr></hr>
</td>
</tr>
<tr><td align="center"><sub>Өsway</sub>
</td>
<td align="center">-</td>
<td align="center">-</td>
<td align="center">-</td>
<td align="center">-</td>
<td align="center">-</td>
<td align="center">-</td>
<td align="center">0.10</td>
<td align="center">0.15</td>
<td align="center">0.08</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p><italic>GOF</italic>
goodness of fit, <italic>NSR</italic>
noise-to-signal ratio averaged over the two perturbation rounds, <italic>PL</italic>
platform, <italic>T</italic>
torque. Amplitude values are given in peak-to-peak values. <italic>f</italic>
are the frequencies of the perturbation signals.</p>
<p>The GOF is reported for each estimate of the model transfer function <inline-formula><mml:math id="M5" name="1743-0003-10-23-i1" overflow="scroll"><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mspace width="0.12em"></mml:mspace>
<mml:mtext mathvariant="normal">and</mml:mtext>
<mml:mspace width="0.24em"></mml:mspace>
<mml:msub><mml:mi>C</mml:mi>
<mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:msub>
<mml:msub><mml:mrow></mml:mrow>
<mml:mrow><mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:math>
</inline-formula>
. The median NSR is reported for the input (joint torques) and output (joint angles) for the model simulations and the balance control experiment, averaged over the two perturbation rounds. Six perturbation conditions were simulated: a platform translation+perturbation torque without and with noise (PLT-Nn and PLT-N), two perturbation forces with and without noise (2F-Nn and 2F-N) and a platform translation with and without noise (PL-N and PL-Nn). For the balance control experiment, the median NSRs of one representative healthy control (HC), the average of the healthy controls and of the PD patient are shown.</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec><title>Experiment</title>
<sec><title>Participants</title>
<p>Seven healthy participants (two female, mean age 65 yrs. std 5.7) and a PD patient (male, 57 yrs) participated in the study. The participants gave written informed consent prior to participation. The protocol was approved by the local medical ethics committee and was in accordance with the Declaration of Helsinki.</p>
</sec>
<sec><title>Apparatus and recording</title>
<p>Two independent perturbations were administered with a computer-controlled six DoF motion platform (Caren, Motek, Amsterdam, The Netherlands) and a custom-built actuated device able to apply perturbing forces in the anterior and posterior direction at the sacrum, called the pusher (see Figure <xref ref-type="fig" rid="F2">2</xref>
). The pusher was attached to the platform and actuated using a series elastic actuator [<xref ref-type="bibr" rid="B34">34</xref>
] controlled with an electro motor (Maxon motor ag, Sachseln, Switzerland). The pusher was force controlled using a custom-built controller in xPC (The Mathworks, Natick, USA) and had a bandwidth of 10 Hz and a maximum torque of 50 Nm. The gravitational pull due to the weight of the pusher was compensated for, such that the participants did not experience additional forces other than the perturbation force and a small force due to the pusher’s inertia. The interaction force in between the subject and the pusher was measured with a six DoF force transducer (ATI-Mini45-SI-580-20).</p>
<p>Body kinematics and the platform movements were measured using motion capture (Vicon Oxford Metrics, Oxford, UK) at a sample frequency of 120 Hz. Reflective spherical markers were attached to the following anatomical landmarks: the first metatarsal, calcaneus, medial malleolus, the sacrum, the manubrium and the last vertebrae of the cervical spine (C7). In addition, a cluster of three markers was attached to both anterior superior iliac spines on the pelvis. Furthermore, one additional marker was attached to the foot and two markers were attached to the lower leg (one on the tibia) to improve the estimation of the ankle joint rotational axis. Also, markers were attached to the knee (just below the lateral epicondyle) and shoulder joints (just in front of the acromion). Three markers were attached to the platform. Reactive forces from both feet were measured with a dual forceplate (AMTI, Watertown, USA), embedded in the motion platform. The signals from the dual forceplate, the 6 DoF force transducer and the perturbation of the pusher were sampled at 600 Hz and stored for further processing.</p>
</sec>
<sec><title>Procedure</title>
<p>During the experiment, participants stood with their arms folded in front of their chest on the dual forceplate and strapped to the pusher, with a strap band that opened with a click buckle, with their eyes open. They were instructed to maintain their balance without moving their feet, while multisine platform movements and multisine force perturbations were applied simultaneously in the forward-backward direction; see ‘Disturbance signals’. Participants wore a safety harness to prevent falling, but it did not constrain movements, provide support or orientation information in any way.</p>
<p>Before any data was recorded, the participants got acquainted to the perturbations. The experimenter determined the maximal amplitude the participant could withstand while keeping the feet flat on the floor and assessed whether the participant could withstand this amplitude for the total of four trials. Four double perturbation trials of 180s were recorded: in the first two trials, the perturbations had the same sign. In the other two trials, the perturbations had opposite signs. If needed, the participants were allowed rest in between trials.</p>
</sec>
</sec>
<sec><title>Disturbance signals</title>
<p>For both the model simulations and the experiment we used the same perturbation signal. During the model simulations the perturbation signal was used to either produce two perturbing forces, or a combination of a platform translation and a torque around the ankle (see ‘Model simulations’). In the experiment the perturbations were applied with a movement platform and an actuated backboard (see ‘Apparatus and recording’).</p>
<p>The perturbation signal was a multisine with a period of 34.13 s (equal to 2<sup>12</sup>
= 4096 samples at a sample rate of 120 Hz) [<xref ref-type="bibr" rid="B4">4</xref>
,<xref ref-type="bibr" rid="B35">35</xref>
]. This signal contained power at 112 frequencies in the range of 0.06–4.25 Hz. To increase the power at the excited frequencies the signal was divided into five frequency bands: 0.06-2.37 Hz (80 frequencies), 2.63-2.84 Hz (8 frequencies), 3.11-3.31 Hz (8 frequencies), 3.57-3.78 Hz (8 frequencies), 4.04-4.25 Hz (8 frequencies). The frequency points outside these frequency bands were not excited. The signal is unpredictable for participants, because the signal consists of many sinusoids. The power of the signal was optimized by crest optimization [<xref ref-type="bibr" rid="B13">13</xref>
].</p>
</sec>
<sec><title>Data analysis</title>
<p>The human body, i.e. the plant, is considered as a double-inverted pendulum, consisting of a leg and a Head-Arms-Trunk (HAT) segment. Stabilizing mechanisms generate ankle and hip torques based on sensory information of the joint angles, (see Figure <xref ref-type="fig" rid="F1">1</xref>
).</p>
<p>The stabilizing mechanisms have passive components such as muscle stiffness, generated by passive muscle properties and tonic activation. The active part incorporates the controller within the CNS (e.g. reflexive muscle activation), muscle activation dynamics, and time-delays, representing the neural signal conduction times.</p>
<p>Movements from the upper body segment will influence the movements of the lower body segment and vice versa due to mechanical coupling [<xref ref-type="bibr" rid="B36">36</xref>
,<xref ref-type="bibr" rid="B37">37</xref>
]. The stabilizing mechanisms have to deal with this mechanical coupling, which is especially expressed in coupling terms between ankles and hips (i.e., <inline-formula><mml:math id="M6" name="1743-0003-10-23-i6" overflow="scroll"><mml:mrow><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mspace width="0.12em"></mml:mspace>
<mml:mtext mathvariant="normal">and</mml:mtext>
<mml:mspace width="0.24em"></mml:mspace>
<mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
). The direct terms <inline-formula><mml:math id="M7" name="1743-0003-10-23-i7" overflow="scroll"><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mspace width="0.12em"></mml:mspace>
<mml:mtext mathvariant="normal">and</mml:mtext>
<mml:mspace width="0.24em"></mml:mspace>
<mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:math>
</inline-formula>
represent the corrective actions of the ankle and hip joint, based on the ankle and hip joint angle. In other words, this system is a multiple input (two joint angles) multiple output (two joint torques) system. When considering the corrective actions of both legs separately, two stabilizing mechanisms are defined; one for each leg [<xref ref-type="bibr" rid="B4">4</xref>
,<xref ref-type="bibr" rid="B38">38</xref>
].</p>
<p>For the model simulations, the perturbations, inputs (joint angles) and outputs (joint torques) of the model were determined for further processing. For the experiment these signals were calculated as described below.</p>
<p>From the recorded movement trajectories of the markers, the position of the center-of-mass (CoM) of the predefined segments and of the whole body and the position of the joints were estimated by custom written software [<xref ref-type="bibr" rid="B39">39</xref>
,<xref ref-type="bibr" rid="B40">40</xref>
]. In short, in each segment a local coordinate frame was determined on the basis of the position of anatomical landmarks, according to the method described by Koopman et al., 1995 [<xref ref-type="bibr" rid="B41">41</xref>
]. The mass, CoM position and the inertia tensor moment of the predefined segments (i.e., feet, legs and HAT) and the joint positions were determined with regression equations [<xref ref-type="bibr" rid="B41">41</xref>
,<xref ref-type="bibr" rid="B42">42</xref>
]. Subsequently, the CoMs were determined as the weighted sum of the separate segment CoM positions [<xref ref-type="bibr" rid="B40">40</xref>
]. From the static trial, the average distance in the sagittal plane from the ankle to the total body CoM (i.e., the length of the pendulum; l<sub>CoM</sub>
) was determined. The sway angle was calculated from l<sub>CoM</sub>
and the horizontal distance from the CoM to the mean position of the ankles. Forces and torques of the force plate and force sensor were filtered with a fourth-order low-pass Butterworth filter with a cut-off frequency of 8 Hz and subsequently resampled to 120 Hz. Forces and torques of the force plate were corrected for the inertia and mass of the top cover [<xref ref-type="bibr" rid="B43">43</xref>
]. On the basis of the corrected forces and torques and recorded body kinematics, ankle and hip joint torques were calculated with inverse dynamics [<xref ref-type="bibr" rid="B40">40</xref>
]. In addition, the applied platform perturbation was reconstructed from the platform markers.</p>
</sec>
<sec><title>Multiple-Input-Multiple-Output closed-loop system identification</title>
<p>To obtain a non-parametric spectral estimate of a two DoF multivariate closed-loop system we adopted a method described by Pintelon and colleagues stating that <bold>two</bold>
different combinations of a periodic excitation signal, <italic>D(k)</italic>
, in <bold>two</bold>
separate experiments should be applied [<xref ref-type="bibr" rid="B13">13</xref>
]. An optimal choice of <italic>D(f)</italic>
(maximizing det (<bold>P</bold>
(f)) using periodic excitation is given by:</p>
<p><disp-formula id="bmcM1"><label>(1)</label>
<mml:math id="M8" name="1743-0003-10-23-i8" overflow="scroll"><mml:mrow><mml:mi mathvariant="bold">P</mml:mi>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
<mml:mo>=</mml:mo>
<mml:mfenced open="[" close="]"><mml:mtable columnalign="center"><mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>1</mml:mn>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mn>1</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>1</mml:mn>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mrow><mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mfenced>
<mml:mi>D</mml:mi>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>With D(<italic>f</italic>
) the two perturbation signals. All calculations were performed in the frequency domain with <italic>f</italic>
the frequency in Hz. This means that in the first round, both inputs were excited with the same periodic excitation, while in the second round the sign of the second perturbation was changed. These perturbations excited the system and the system responded with movements (joint angles) and torques (corrective joint torques) at the frequencies of the perturbation signal. Corrective torques are the torques that restore the body’s equilibrium in response to motor and sensor noise and the perturbation signal. The estimate from the perturbation to inputs (joint angles) and outputs (corrective joint torques) was first obtained from:</p>
<p><disp-formula id="bmcM2"><label>(2)</label>
<mml:math id="M9" name="1743-0003-10-23-i9" overflow="scroll"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>G</mml:mi>
<mml:mo stretchy="true">^</mml:mo>
</mml:mover>
<mml:mi mathvariant="italic">py</mml:mi>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
<mml:mo>=</mml:mo>
<mml:mi>Y</mml:mi>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
<mml:msup><mml:mi>P</mml:mi>
<mml:mrow><mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>With <italic>Y(f),</italic>
a two-by-two matrix with on the first column the responses of the first perturbation round (i.e., ankle and hip joint angles or torques) and on the second column the responses of the second perturbation round. <inline-formula><mml:math id="M10" name="1743-0003-10-23-i10" overflow="scroll"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>G</mml:mi>
<mml:mo stretchy="true">^</mml:mo>
</mml:mover>
<mml:mi mathvariant="italic">py</mml:mi>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>
denotes the estimate of the cross spectral density (CSD) of the disturbance and the outputs (joint angles and joint torques).</p>
<p>Subsequently, the stabilizing mechanisms were estimated using the joint input-output approach [<xref ref-type="bibr" rid="B12">12</xref>
]:</p>
<p><disp-formula id="bmcM3"><label>(3)</label>
<mml:math id="M11" name="1743-0003-10-23-i11" overflow="scroll"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>C</mml:mi>
<mml:mo stretchy="true">^</mml:mo>
</mml:mover>
<mml:mi mathvariant="italic">θTc</mml:mi>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi mathvariant="italic">f</mml:mi>
</mml:mfenced>
<mml:mo>=</mml:mo>
<mml:mo>-</mml:mo>
<mml:msub><mml:mover accent="true"><mml:mi>G</mml:mi>
<mml:mo stretchy="true">^</mml:mo>
</mml:mover>
<mml:mtext mathvariant="italic">pTc</mml:mtext>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
<mml:msubsup><mml:mover accent="true"><mml:mi>G</mml:mi>
<mml:mo stretchy="true">^</mml:mo>
</mml:mover>
<mml:mi mathvariant="italic">pθ</mml:mi>
<mml:mrow><mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mspace width="0.48em"></mml:mspace>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>With <inline-formula><mml:math id="M12" name="1743-0003-10-23-i12" overflow="scroll"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>G</mml:mi>
<mml:mo stretchy="true">^</mml:mo>
</mml:mover>
<mml:mtext mathvariant="italic">pTc</mml:mtext>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>
and <inline-formula><mml:math id="M13" name="1743-0003-10-23-i13" overflow="scroll"><mml:mrow><mml:msubsup><mml:mover accent="true"><mml:mi>G</mml:mi>
<mml:mo stretchy="true">^</mml:mo>
</mml:mover>
<mml:mi mathvariant="italic">pθ</mml:mi>
<mml:mrow><mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mspace width="0.24em"></mml:mspace>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula>
the estimated CSD from the perturbations to the corrective torques and from the perturbations to the joint angles<italic>.</italic>
Note that C is a two-by-two matrix, see also Figure <xref ref-type="fig" rid="F1">1</xref>
. p is a vector with the two disturbances, θ(<italic>f</italic>
) is a vector with ankle and hip joint angles, and T<sub>c</sub>
<italic>(f)</italic>
is a vector with ankle and hip joint torques for each frequency <italic>f</italic>
; all expressed as Fourier coefficients. The method assumes that the system does not change between the two separate perturbation rounds, with the two different combinations of the periodic excitation signal (Eq. 1).</p>
<p>For both the model simulations and the experiment, data was obtained for eight response cycles of the perturbation signal for each perturbation round. Subsequently, the data were Fourier transformed and only the Fourier coefficients at the excited frequencies were used for further processing. These were averaged over the eight cycles, and the average Fourier coefficients were used to calculate the power- and cross spectral density (PSD and CSD, respectively). The PSDs and the CSDs were smoothed by averaging over four adjacent frequency points [<xref ref-type="bibr" rid="B44">44</xref>
]. The FRFs were calculated according to Eq. 2 and 3 to obtain a non-parametric spectral estimate of the total stabilizing mechanisms.</p>
<p>As the corrective torque which has to be delivered by the participants is dependent on gravity, all FRFs were normalized for the participants’ mass and length, i.e. the gravitational stiffness (<italic>mgl</italic>
), with <italic>m</italic>
the total body mass, <italic>l</italic>
the length of the pendulum (from the ankles to the Center of Mass; CoM), and <italic>g</italic>
the gravitational constant. The average FRF over all participants was obtained by taking the mean over the individual normalized FRFs. Note that, as we used a dual forceplate in the experiment, the experimentally obtained Fourier coefficients of the left and right FRFs were added to obtain the total FRFs.</p>
<sec><title>Reliability of the estimated MIMO frequency response functions</title>
<p>To determine whether the above described MIMO closed-loop system identification technique gives reliable estimates of the stabilizing mechanisms, several indicators were calculated (described below).</p>
</sec>
<sec><title>Goodness of estimate</title>
<p>For the model simulations, the goodness of fit (GOF) was determined by the object function [<xref ref-type="bibr" rid="B33">33</xref>
]. This function compared the theoretical Transfer Function (TF) as incorporated in the model (Appendix A) with the estimated FRF as obtained in the model simulations:</p>
<p><disp-formula id="bmcM4"><label>(4)</label>
<mml:math id="M14" name="1743-0003-10-23-i14" overflow="scroll"><mml:mrow><mml:mtext mathvariant="italic">GOF</mml:mtext>
<mml:mo>=</mml:mo>
<mml:msup><mml:mrow><mml:mstyle displaystyle="true"><mml:munder><mml:mo mathsize="big">∑</mml:mo>
<mml:mi>f</mml:mi>
</mml:munder>
<mml:mrow><mml:mfenced open="(" close=")"><mml:mfrac><mml:mrow><mml:mfenced open="|" close="|"><mml:mrow><mml:msub><mml:mi>H</mml:mi>
<mml:mtext mathvariant="italic">theoretical</mml:mtext>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
<mml:mo>-</mml:mo>
<mml:msub><mml:mi>H</mml:mi>
<mml:mtext mathvariant="italic">estimated</mml:mtext>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow><mml:mfenced open="|" close="|"><mml:mrow><mml:msub><mml:mi>H</mml:mi>
<mml:mtext mathvariant="italic">theoretical</mml:mtext>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
</mml:mfenced>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>Where a perfect estimation of the transfer function will result in a GOF of zero, i.e., the lower the GOF the better the estimation.</p>
</sec>
<sec><title>Noise-to-signal ratio</title>
<p>As a result of a periodic perturbation to the balance control system, the system’s response was periodic, while time-variant behavior and/or noise resulted in a stochastic, i.e. non-periodic response [<xref ref-type="bibr" rid="B35">35</xref>
]. The ratio of the non-periodic (also called the remnant) and periodic response is expressed by the noise-to-signal ratio (NSR):</p>
<p><disp-formula id="bmcM5"><label>(5)</label>
<mml:math id="M15" name="1743-0003-10-23-i15" overflow="scroll"><mml:mrow><mml:mi mathvariant="italic">NS</mml:mi>
<mml:msub><mml:mi>R</mml:mi>
<mml:mi>U</mml:mi>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
<mml:mo>=</mml:mo>
<mml:mfrac><mml:mrow><mml:msubsup><mml:mi>σ</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
</mml:mrow>
<mml:mfenced open="|" close="|"><mml:mrow><mml:msub><mml:mi>U</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:msup><mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>Where <italic>U</italic>
<sub><italic>p</italic>
</sub>
(<italic>f</italic>
) represents the periodic response and <italic>σ</italic>
<sub><italic>u</italic>
</sub>
<sup>2</sup>
(<italic>f</italic>
) the variance of the remnant. The NSR was calculated in the frequency domain with <italic>f</italic>
frequency in Hz. When multiple realizations are simulated or recorded, the periodic response is obtained by calculating the average over the realizations; the remnant can be estimated by calculating the variance over the realizations.</p>
<p>A small NSR indicates low variability of the system’s response to the perturbation signal over the multiple periods. This indicates time invariant behavior and a low presence of noise [<xref ref-type="bibr" rid="B35">35</xref>
]. As such, it gives insight into whether linear, time-invariant, system-identification methods can be used and it gives an estimate of how well the system is perturbed. More importantly, it quantifies how reliable the estimate of the stabilizing mechanism is. For example, a NSR of 1 indicates that the recorded data contains 50 percent response and 50 percent remnant. This means that describing the system as a deterministic linear time invariant (LTI) system (expressed by the estimated FRFs in this study) explains 50 percent of the recorded data.</p>
</sec>
<sec><title>Single-Input-Single-Output frequency response functions</title>
<p>To assess whether the two perturbation rounds in the experiment lead to a change in strategy of the participants (i.e., time-variant behavior), we estimated Single-Input-Single-Output (SISO) FRFs from sway angle to ankle torque and from sway angle to hip torque with the joint-input-output-method (Eq. 3; [<xref ref-type="bibr" rid="B4">4</xref>
]).</p>
</sec>
<sec><title>Coherence</title>
<p>For the SISO FRFs the (magnitude-squared) coherence was calculated between the input signal (perturbations) and output signals (ankle torques, hip torques, and body sway).</p>
<p><disp-formula id="bmcM6"><label>(6)</label>
<mml:math id="M16" name="1743-0003-10-23-i16" overflow="scroll"><mml:mrow><mml:msubsup><mml:mi>γ</mml:mi>
<mml:mrow><mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:msup><mml:mfenced open="|" close="|"><mml:msub><mml:mi>ϕ</mml:mi>
<mml:mrow><mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mfenced>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>.</mml:mo>
<mml:msup><mml:mfenced open="[" close="]"><mml:mrow><mml:msub><mml:mi>ϕ</mml:mi>
<mml:mrow><mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
<mml:mo>.</mml:mo>
<mml:msub><mml:mi>ϕ</mml:mi>
<mml:mrow><mml:mi>y</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
</mml:mrow>
</mml:mfenced>
<mml:mrow><mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>In which p represents the platform disturbance and y an output signal. <italic>φ</italic>
<sub><italic>p</italic>
,<italic>y</italic>
</sub>
, <italic>φ</italic>
<sub><italic>p</italic>
,<italic>p</italic>
</sub>
and <italic>φ</italic>
<sub><italic>y</italic>
,<italic>y</italic>
</sub>
are the CSDs and PSDs of the perturbation and output signals. By definition, coherence varies between 0 and 1, where coherence close to one indicates a low noise level and time-invariant behaviour.</p>
</sec>
</sec>
<sec><title>Balance contribution of the left and right body side</title>
<p>For the experiment we determined the relative contribution of each ankle and hip joint to the total amount of generated corrective torque to resist the perturbations by calculating the contribution of the gain and phase of each leg to the gain and phase of the total body [<xref ref-type="bibr" rid="B4">4</xref>
]:</p>
<p><disp-formula id="bmcM7"><label>(7)</label>
<mml:math id="M17" name="1743-0003-10-23-i17" overflow="scroll"><mml:mrow><mml:mtext mathvariant="italic">Contributio</mml:mtext>
<mml:msub><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
<mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="italic">FR</mml:mi>
<mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
<mml:mo>•</mml:mo>
<mml:mi mathvariant="italic">FR</mml:mi>
<mml:msub><mml:mi>F</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
</mml:mrow>
<mml:msup><mml:mfenced open="‖" close="‖"><mml:mrow><mml:mi mathvariant="italic">FR</mml:mi>
<mml:msub><mml:mi>F</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>f</mml:mi>
</mml:mfenced>
</mml:mrow>
</mml:mfenced>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mfrac>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>With FRF<sub>l,r</sub>
the left or right FRF and FRF<sub>t</sub>
the total FRF. The ● indicates the dot product of the FRFs. In this way the contribution of the left or right leg to the total balance control was expressed as a proportion. For example, a proportion of 0.8 for the left leg means that the left leg contributed for 80% to the total body stabilization. This was done for each separate MIMO FRF (see Eq. 2 and 3).</p>
</sec>
</sec>
<sec sec-type="results"><title>Results</title>
<sec><title>Model simulations</title>
<p>Table <xref ref-type="table" rid="T1">1</xref>
shows the GOF values and NSRs of the different simulations. In case of no sensor and measurement noise, a platform acceleration in combination with a perturbation torque around the ankle (PLT-Nn) gave the same results as two perturbation forces (2F-Nn). In these conditions, the small GOF values indicated that the stabilizing mechanisms were well estimated.</p>
<p>Adding sensor (pink) and measurement (white) noise to the simulations resulted in slightly worse estimations; the GOF values increased. However, the stabilizing mechanisms were still correctly estimated (see Figure <xref ref-type="fig" rid="F3">3</xref>
). Platform acceleration and a perturbation torque around the ankle (PLT-N) resulted in better estimations than two perturbation forces (2F-N).</p>
<fig id="F3" position="float"><label>Figure 3</label>
<caption><p><bold>Theoretical transfer function and estimated frequency response functions </bold>
<inline-formula><mml:math id="M18" name="1743-0003-10-23-i19" overflow="scroll"><mml:mstyle mathvariant="bold"><mml:mo stretchy="true">(</mml:mo>
<mml:mrow><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo mathvariant="bold">,</mml:mo>
<mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo mathvariant="bold">,</mml:mo>
<mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mspace width="0.12em"></mml:mspace>
<mml:mtext mathvariant="bold">and</mml:mtext>
<mml:mspace width="0.24em"></mml:mspace>
<mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">)</mml:mo>
</mml:mstyle>
</mml:math>
</inline-formula>
<bold>. </bold>
Results of the model simulations in the presence of pink and white noise during the condition with a platform acceleration and perturbation torque around the ankle are depicted (PLT-N; open dots) and during the condition with one perturbation (PL-N; solid grey dots). The bold solid line represents the model transfer function of the stabilizing mechanisms. Applying two independent perturbations in combination with a multivariate closed-loop system identification method resulted in a correct estimation of the stabilizing mechanisms (FRF 2 perturbations, indicated with open dots), whereas one perturbation resulted in an erroneous estimate (FRF 1 perturbation, indicated with solid grey dots).</p>
</caption>
<graphic xlink:href="1743-0003-10-23-3"></graphic>
</fig>
<p>Figure <xref ref-type="fig" rid="F3">3</xref>
shows the model transfer function and the estimated frequency response functions of the PLT-N and PL-N conditions. Clearly, applying two perturbations resulted in very good estimations (see also Table <xref ref-type="table" rid="T1">1</xref>
). However, applying only <bold>one</bold>
perturbation to estimate the stabilizing mechanisms of a MIMO system resulted in incorrect estimates, although the responses were time-invariant as shown by the NSRs values. In other words, applying one perturbation resulted in biased estimates of the stabilizing mechanisms. Note that the GOF values of the PL-Nn and PL-N simulations were of the same magnitude.</p>
</sec>
<sec><title>Experiment</title>
<sec><title>Time series</title>
<p>Figure <xref ref-type="fig" rid="F4">4</xref>
shows the perturbations and the response of one representative healthy participant to the applied perturbations. For the healthy controls, the average peak-to-peak amplitudes were 0.068 m (<italic>std</italic>
: 0.005) for the platform and 18.4Nm (<italic>std</italic>
: 1.05) for the pusher. For the PD patient the peak-to-peak amplitudes were 0.06 m and 20 Nm, see also Table <xref ref-type="table" rid="T1">1</xref>
. In general, the participant responded in a consistent fashion as indicated by the low standard deviation over the adjacent segments and corresponding low NSRs (Table <xref ref-type="table" rid="T1">1</xref>
). The average median NSRs (of both perturbation rounds) of all healthy participants were: 0.24 (ankle angle); 0.28 (hip angle); 0.15 (sway angle); 0.56 (ankle torque) and 0.35 (hip torque). Hence, on average, the healthy participants responded in a consistent fashion (see Figure <xref ref-type="fig" rid="F4">4</xref>
). Note, however, that the response of the sway angle was even more consistent than the joint angles. The PD patient was able to complete the balance control experiment without any problems. He could withstand the perturbations and the duration of the trials. The patient also responded in a time-invariant fashion. He had even slightly lower NSRs than the healthy controls.</p>
<fig id="F4" position="float"><label>Figure 4</label>
<caption><p><bold>Timeseries (left panels) and NSRs (right panels) of the first perturbation round of one representative participant. </bold>
From top to bottom: platform perturbation, pusher perturbation, ankle angle, hip angle, sway angle, ankle torque, and hip torque. The mean is depicted by the solid line and the standard deviation over the eight cycles by the grey area. The black line in the right panels depicts NSR=1. Ideally, the average NSR of the responses remains below one. The responses of the participant were consistent, as evidenced by small standard deviations over the adjacent segments and low NSRs. This means that a large part of the data is captured by the time-invariant MIMO system identification technique.</p>
</caption>
<graphic xlink:href="1743-0003-10-23-4"></graphic>
</fig>
</sec>
<sec><title>Single-Input-Single-Ouput frequency response functions</title>
<p>Figure <xref ref-type="fig" rid="F5">5</xref>
shows the SISO FRFs from sway angle to ankle torque and hip torque of each perturbation round of the healthy controls. Gain and phase of the FRFs and coherence of the joint torques were similar for both perturbation rounds. There was a small discrepancy between the gain of the FRFs at the lower and higher frequencies, but this can be attributed to a less periodic response at these frequencies (indicated by a decreased coherence), and hence the reliability of the FRFs decreased. In addition, the coherence of the sway angle was lower in the second perturbation round, especially at frequencies below 0.7 Hz. This resulted in a slightly worse estimate of the FRFs, compared to the first perturbation round. Similar results were found for the PD patient (data not shown).</p>
<fig id="F5" position="float"><label>Figure 5</label>
<caption><p><bold>Single-Input-Single-Output frequency response functions and coherences of first and second perturbation round. </bold>
The left panel depicts the FRF from the sway angle to the ankle joint torque; the right panel the FRF from the sway angle to the hip joint torque. The lower panels show the coherence between the perturbation, and the ankle joint torque, hip joint torque, and the sway angle. Similar gains and phases of the frequency response functions of the first and second perturbation round indicate that participants did not change their balance control strategy.</p>
</caption>
<graphic xlink:href="1743-0003-10-23-5"></graphic>
</fig>
</sec>
<sec><title>Multiple-Input-Multiple-Output frequency response functions</title>
<p>The MIMO FRFs of the stabilizing mechanisms of healthy controls and the PD patient are shown in Figure <xref ref-type="fig" rid="F6">6</xref>
. The total of four FRFs are shown: from ankle angle to ankle torque <inline-formula><mml:math id="M19" name="1743-0003-10-23-i20" overflow="scroll"><mml:mfenced open="(" close=")"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mfenced>
</mml:math>
</inline-formula>
, from ankle angle to hip torque <inline-formula><mml:math id="M20" name="1743-0003-10-23-i21" overflow="scroll"><mml:mfenced open="(" close=")"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mrow><mml:mi>A</mml:mi>
<mml:mo>→</mml:mo>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mfenced>
</mml:math>
</inline-formula>
,from hip angle to ankle torque <inline-formula><mml:math id="M21" name="1743-0003-10-23-i22" overflow="scroll"><mml:mfenced open="(" close=")"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mrow><mml:mi>H</mml:mi>
<mml:mo>→</mml:mo>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mfenced>
</mml:math>
</inline-formula>
and from hip angle to hip torque <inline-formula><mml:math id="M22" name="1743-0003-10-23-i23" overflow="scroll"><mml:mfenced open="(" close=")"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mrow><mml:mi>H</mml:mi>
<mml:mo>→</mml:mo>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mfenced>
</mml:math>
</inline-formula>
representing the multivariate stabilizing mechanisms of the participants.</p>
<fig id="F6" position="float"><label>Figure 6</label>
<caption><p><bold>Multiple-Input-Multiple-Output frequency response functions of the stabilizing mechanisms. </bold>
The solid line represents the average of the healthy participants, with the shaded area indicating the 95% confidence interval and the dotted line the Parkinson patient.</p>
</caption>
<graphic xlink:href="1743-0003-10-23-6"></graphic>
</fig>
<p>In general, for the healthy controls, the relationship between the ankle joint angle and the ankle joint torque remained roughly constant over the whole frequency range, whereas the hip joint gain was high at the low frequencies (<0.7 Hz), low at the mid frequencies (0.7-2 Hz) and increased at higher frequencies (>2 Hz). The coupling between the ankle joint angle to hip torque increased with frequency, whereas the <inline-formula><mml:math id="M23" name="1743-0003-10-23-i24" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
FRF remained roughly constant over the frequency range. The PD patient showed similar gain patterns for the <inline-formula><mml:math id="M24" name="1743-0003-10-23-i25" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
and <inline-formula><mml:math id="M25" name="1743-0003-10-23-i26" overflow="scroll"><mml:mrow><mml:mspace width="0.12em"></mml:mspace>
<mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
FRFs, but for the <inline-formula><mml:math id="M26" name="1743-0003-10-23-i27" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
and <inline-formula><mml:math id="M27" name="1743-0003-10-23-i28" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
FRFs the gains at the lower frequencies were much higher than those of the healthy controls, indicating an increased postural stiffness in the PD patient.</p>
<p>The phase of the <inline-formula><mml:math id="M28" name="1743-0003-10-23-i29" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
and <inline-formula><mml:math id="M29" name="1743-0003-10-23-i30" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
FRFs of both the HCs and the PD patient decreased with increasing frequency (i.e. a phase lag), indicating a neural time delay. The phase of the <inline-formula><mml:math id="M30" name="1743-0003-10-23-i31" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
showed a phase shift of about +360° around 1.5 Hz. For the phase the largest difference between the HCs and PD patient is found in the <inline-formula><mml:math id="M31" name="1743-0003-10-23-i32" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
FRF. Both groups showed a negative phase shift around 2 Hz, but this shift was larger in the PD patient.</p>
</sec>
<sec><title>Balance contribution of the left and right body side</title>
<p>In healthy controls both legs contributed equally to the body stabilization; the proportion for each FRF was 0.5. However, for the patient, the right leg contributed more to the balance control than the left leg, and this was the case for both the ankle and the hip joint. Note that the patient is clearly outside the 95% confidence interval (CI) of HCs (see Figure <xref ref-type="fig" rid="F7">7</xref>
). It can be seen that for the FRFs <inline-formula><mml:math id="M32" name="1743-0003-10-23-i33" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
and <inline-formula><mml:math id="M33" name="1743-0003-10-23-i34" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
the proportion between both legs was similar over the whole frequency range, whereas for the FRFs <inline-formula><mml:math id="M34" name="1743-0003-10-23-i35" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
and <inline-formula><mml:math id="M35" name="1743-0003-10-23-i36" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
the asymmetry decreased with increasing frequency and eventually disappeared above 1 Hz.</p>
<fig id="F7" position="float"><label>Figure 7</label>
<caption><p><bold>The average contribution of the right leg of the healthy controls and of the PD patient to each Multiple-Input-Multiple-Output frequency response function. </bold>
The average contribution of the right leg for the healthy controls (HC) is shown by the solid line (mean) and the grey area (95% confidence interval). The patient (dashed line) clearly controlled his balance asymmetrically, with the right leg producing more corrective torque than the left leg to resist the perturbations.</p>
</caption>
<graphic xlink:href="1743-0003-10-23-7"></graphic>
</fig>
</sec>
</sec>
</sec>
<sec sec-type="discussion"><title>Discussion</title>
<p>Bipedal upright stance is multisegmental and requires the coordinated activity of multiple joints, including the ankles and hips. We developed a system identification method to investigate the balance control contribution of the ankles and hips of the left and right leg separately of individual subjects. To investigate both DoF and their interactions in a feedback loop, two perturbations are required. Therefore, in our lab we developed an actuated pusher placed on a motion platform to provide these independent perturbations. Using model simulations, it was demonstrated that these perturbations can identify the system reliably. The (clinical) applicability of the method was demonstrated in seven healthy controls and a PD patient.</p>
<sec><title>Evaluation of the Multiple-Input-Multiple-Output identification method</title>
<sec><title>Model</title>
<p>The model simulations indicated that two independent perturbations are necessary to identify the stabilizing mechanisms of a two DoF MIMO system. Applying only one perturbation to the model resulted in biased estimates of the stabilizing mechanisms. This bias was not influenced by noise level. Two configurations with two independent perturbations were evaluated, that is, a combination of a platform acceleration and a disturbance torque around the ankle and two perturbation forces in the forward-backward direction. The stabilizing mechanism were estimated very well with both configurations, and the differences between the two approaches were small. After adding sensor and measurement noise to the simulation, the stabilizing mechanisms were still estimated well with both configurations, with a slightly worse result when perturbing with two forces. In short, the model simulations showed that the implemented MIMO identification method correctly identified the stabilizing mechanisms indicated by small differences between the theoretical transfer function of the model and the estimated frequency response functions, even in the presence of realistic levels of noise.</p>
</sec>
<sec><title>Experiment</title>
<p>In the experiment with human participants, the perturbations were applied with a motion platform and a custom-made actuated backboard (i.e., the first configuration of the model simulations, which gave slightly better results). The amplitudes of the applied perturbations were easy to withstand, both for the healthy participant as for the PD patient, making the method suitable for use in a large range of participants.</p>
<p>To be clinically relevant, the MIMO identification method should be able to reliably estimate the different segmental contributions to the total balance control of a single participant. First, the quality of our estimation was expressed in the NSR. A low NSR indicates that a large percentage of the data is captured by the estimated stabilizing mechanisms. The NSR gives the ratio between the (periodic) response to the perturbations and the remnant. As the method assumes a LTI system, remnants can be due to a) nonlinearities of the perturbed system, b) time-variant system behavior, c) unmeasured system noise, and d) measurement noise [<xref ref-type="bibr" rid="B35">35</xref>
]. Note that a low NSR does not necessarily mean that the system is linear.</p>
<p>The average NSRs of the healthy controls of the ankle and hip joint angle and of the ankle and hip torque were 0.24 and 0.28, 0.56 and 0.35, respectively. This results in an average NSR of 0.36, meaning that about 74% (i.e., 1/(1+0.36)) of the response is captured by the estimated stabilizing feedback mechanisms. The PD patient had lower NSRs than the average healthy subject, indicating less variability over the repetitions of the perturbation signal. Unfortunately, it is impossible to compare our results with respect to reliability of the estimates with other studies using two perturbations as these studies did not report any parameters quantifying reliability and only report averages over participants.</p>
<p>Secondly, the used MIMO identification scheme consists of two perturbation rounds: in one round, the perturbations have the same sign, while in the other round, the perturbations have opposite signs. To determine whether participants did not change their balance control behaviour in both perturbation rounds, SISO FRFs for the ankle and hip joints were determined. Gain, phase, and joint torque coherence were comparable for both perturbation rounds; hence the participants did not change their balance control strategies. Note that the sway angle and ankle and hip coherence was lower in the second perturbation round at the lower frequencies (<0.7 Hz), indicating less time invariant behavior (possibly due to fatigue) and/or a higher noise level. This could be due to not randomizing the perturbation rounds. However, this did not result into quantitatively different behavior, as the SISO FRFs were similar for both perturbation rounds.</p>
<p>The last requirement of the method is that it is able to distinguish between the balance contributions of the separate legs to investigate balance asymmetry. By using a dual forceplate and calculating the FRFs for each leg separately, balance control asymmetries were detected in a PD patient.</p>
<p>Perturbing with a different set-up (for example two push-pull rods), at different locations or by using different signals could have elicited different responses. This is not an artifact of the method, but reflects the adaptability of the nervous system [<xref ref-type="bibr" rid="B45">45</xref>
].</p>
<p>In sum, two unique features of the presented method are the applicability on the individual level and separation of balance contribution of each body side. The first was accomplished by applying multisine perturbation signals, which have the advantage of improving the estimation of FRFs because they concentrate signal power in a limited set of frequencies and are periodic. This results in reliable <italic>individual</italic>
results and shorter measurement times [<xref ref-type="bibr" rid="B35">35</xref>
,<xref ref-type="bibr" rid="B46">46</xref>
]. Secondly, by measuring the reaction forces of each foot with a dual forceplate balance control, asymmetries can be detected.</p>
</sec>
</sec>
<sec><title>Comparison with other multivariate methods</title>
<p>Most published studies investigating the multivariate nature of balance control do not use perturbations, or use only one perturbation [<xref ref-type="bibr" rid="B15">15</xref>
,<xref ref-type="bibr" rid="B16">16</xref>
,<xref ref-type="bibr" rid="B18">18</xref>
-<xref ref-type="bibr" rid="B20">20</xref>
]. However, the model simulations indicated that two independent perturbations are required to estimate the stabilizing mechanisms of a multivariate balance control system. Using only one perturbation in the model simulations gave biased and erroneous results.</p>
<p>The other two available methods [<xref ref-type="bibr" rid="B2">2</xref>
,<xref ref-type="bibr" rid="B21">21</xref>
] differ from the presented method in this paper. Our method is non-parametric, where Fujisawa [21] used a parametric method (ARMAX model structure). A parametric method has the advantage that it, theoretically, can better separate the measurement noise from the actual signals. However, knowledge about the structure of the system is required. A non-parametric method has the advantage that no prior knowledge of the system is needed [<xref ref-type="bibr" rid="B14">14</xref>
]. Despite the differences between the applied methods in [<xref ref-type="bibr" rid="B21">21</xref>
] and in this manuscript the obtained FRFs are in the same range for the low frequencies. That is the <inline-formula><mml:math id="M36" name="1743-0003-10-23-i37" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
FRF starts at unity gain which slightly increases with frequency, whereas the <inline-formula><mml:math id="M37" name="1743-0003-10-23-i38" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
<sub>,</sub>
<inline-formula><mml:math id="M38" name="1743-0003-10-23-i39" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
and <inline-formula><mml:math id="M39" name="1743-0003-10-23-i40" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
have lower gains at the lower frequencies (~0.05) but also increase with frequency. Note that we can compare the FRFs only in the low frequency range, as Fujisawa and colleagues [<xref ref-type="bibr" rid="B21">21</xref>
] used a perturbation signal with frequencies up to 0.83 Hz.</p>
<p>Recently, Jeka and colleagues [<xref ref-type="bibr" rid="B2">2</xref>
,<xref ref-type="bibr" rid="B47">47</xref>
] used an approach to investigate task goals for upright stance, similar to the one presented here. There are, however, a few fundamental differences. First of all, Kiemel et al. [<xref ref-type="bibr" rid="B2">2</xref>
] used filtered white noise as a perturbation signal, whereas we used multisines, which have the advantage of improving the individual estimate of the stabilizing mechanisms. Also, they define the input of the stabilizing mechanisms as weighted EMG signals of the anterior and posterior body sides, whereas we have used joint torques. Therefore, the method presented by Jeka and colleagues [<xref ref-type="bibr" rid="B2">2</xref>
,<xref ref-type="bibr" rid="B47">47</xref>
] focused on separating the reflexive from the total contribution of the stabilizing mechanisms by measuring EMG signals. This measurement and analysis of EMG signals can easily be added to our method. In order to detect balance control asymmetries, joint torques (which we have used in this study) are more suitable than EMG signals. Differences in EMG amplitude can for example, be due to different electrode placement on contralateral legs, different skin conductivity or due to different background activity, making these signals more prone to measurement artifacts. In addition, joint torques are more suitable for investigating joint stiffness.</p>
</sec>
<sec><title>Applications</title>
<sec><title>Multisegmental balance control</title>
<p>With a MIMO method and by applying multiple perturbations, multisegmental balance control strategies and the interplay between the joints can be investigated. This approach can be used to test hypotheses about the role of the different joints, and also of sensory information [<xref ref-type="bibr" rid="B48">48</xref>
].</p>
</sec>
<sec><title>Clinical applications</title>
<p>It has been suggested that PD patients have a decreased intersegmental coordination [<xref ref-type="bibr" rid="B29">29</xref>
,<xref ref-type="bibr" rid="B30">30</xref>
] or an increased hip stiffness [<xref ref-type="bibr" rid="B25">25</xref>
,<xref ref-type="bibr" rid="B28">28</xref>
]. In this study, the most pronounced difference between the healthy controls and the PD patients was found in the <inline-formula><mml:math id="M40" name="1743-0003-10-23-i41" overflow="scroll"><mml:msub><mml:mi>C</mml:mi>
<mml:mrow><mml:msub><mml:mi>θ</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>→</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
FRF, that is, the PD patient had a higher gain at the lower frequencies, indicating an increased hip stiffness. With our method we can now distinguish between coordination between the upper and lower body and of coordination of the upper/lower body separately. Further investigation of intersegmental coordination could lead to a better understanding of the pathophysiology of balance impairments in PD and possibly improve intervention programs.</p>
<p>The PD patient asymmetrically controlled his balance with both the ankle and the hip joint. This means that one leg contributed more to body stability than the other leg. Balance control asymmetries have been shown before in PD patients, both during quiet stance [<xref ref-type="bibr" rid="B32">32</xref>
,<xref ref-type="bibr" rid="B49">49</xref>
] and with a single DoF approach [<xref ref-type="bibr" rid="B50">50</xref>
], but not taking into account the contribution of the hip joint. The balance asymmetry was also present in the intersegmental coordination, i.e. in the FRFs from ankle angle to hip torque and from hip angle to ankle torque. Hence, our new method has the advantage of assessing balance control asymmetries during perturbations, considering the role of the hip joint and the interplay between joints. This creates the possibility of assessing differences in balance control contribution between distal and proximal joints, and it can be investigated whether balance asymmetries influence intersegmental coordination. Further necessary research in a large group of PD patients and healthy matched controls should demonstrate the (potential) clinical value of the new method.</p>
<p>Furthermore, balance control can also be asymmetrical in stroke patients [<xref ref-type="bibr" rid="B4">4</xref>
,<xref ref-type="bibr" rid="B51">51</xref>
]. During the recovering process, restoration of the paretic body side and/or compensation in the non-paretic body side may contribute to improved balance maintenance. Investigating balance control asymmetries provides the possibility of investigating different recovery and compensation strategies during the rehabilitation process.</p>
</sec>
</sec>
</sec>
<sec sec-type="conclusions"><title>Conclusions</title>
<p>We presented a new method to identify the multisegmental stabilizing mechanisms in human stance control using non-parametric system identification techniques and evaluated its performance. Model simulations showed that the newly presented method correctly and reliably estimated the balance control contribution of the ankle and hip joints and interactions between the segments. A balance control experiment showed the application in both healthy and pathological participants. Furthermore, the method can distinguish between the balance control contribution of each ankle and hip joint separately. Taken together, this can be used to create insights into the pathophysiology of postural instability and asymmetry in patients and possibly aid to develop and evaluate treatments.</p>
</sec>
<sec><title>Appendix A</title>
<sec><title>Human balance control model</title>
<p>To test the feasibility of the MIMO identification method, model simulations were performed. A two-DoF mechanical model with a stabilizing mechanism was derived, Figure 8 shows the general model; the derivation of the subsystems is described below.</p>
</sec>
<sec><title>Derivation of the equations of motion</title>
<p>Kane’s method (TMT method) was used to derive the equations of motions. In this method, the principles of virtual power are used to rewrite the Newton-Euler equations.</p>
<p>The first step is to define the degrees of freedom (generalized coordinates, q<sub>t</sub>
) of the model, which are in this case the support surface (S<sub>x</sub>
) and segment angles: θ<sub>1</sub>
and θ<sub>2.</sub>
</p>
<p>Secondly, a transformation matrix T<sub>j</sub>
, which describes the Cartesian coordinates of the center of mass of the segments (x,y coordinates) in the degrees of freedom of the system (i.e., the generalized coordinates S<sub>x</sub>
, θ<sub>1</sub>
and θ<sub>2</sub>
) is defined:</p>
<p><disp-formula id="bmcM8"><label>(8)</label>
<mml:math id="M41" name="1743-0003-10-23-i42" overflow="scroll"><mml:mrow><mml:msub><mml:mi>T</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfenced open="[" close="]"><mml:mtable columnalign="left"><mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:msub><mml:mi>S</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:mrow><mml:msub><mml:mi>S</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>d</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>*</mml:mo>
<mml:mo>cos</mml:mo>
<mml:mfenced open="(" close=")"><mml:msub><mml:mi>θ</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mfenced>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:mrow><mml:msub><mml:mi>d</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>*</mml:mo>
<mml:mo>sin</mml:mo>
<mml:mfenced open="(" close=")"><mml:msub><mml:mi>θ</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mfenced>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:msub><mml:mi>θ</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:mrow><mml:msub><mml:mi>S</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>l</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>*</mml:mo>
<mml:mo>cos</mml:mo>
<mml:mfenced open="(" close=")"><mml:msub><mml:mi>θ</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mfenced>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>*</mml:mo>
<mml:mo>cos</mml:mo>
<mml:mfenced open="(" close=")"><mml:msub><mml:mi>θ</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mfenced>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:mrow><mml:msub><mml:mi>l</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>*</mml:mo>
<mml:mo>sin</mml:mo>
<mml:mfenced open="(" close=")"><mml:msub><mml:mi>θ</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mfenced>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>*</mml:mo>
<mml:mo>sin</mml:mo>
<mml:mfenced open="(" close=")"><mml:msub><mml:mi>θ</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mfenced>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:msub><mml:mi>θ</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mfenced>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>With d<sub>1</sub>
and d<sub>2</sub>
the distance to the segment’s center of gravity from the distal point, l<sub>1</sub>
denotes the length of the leg segment. This matrix can be differentiated to the generalized coordinates to obtain the first (T<sub>j,t</sub>
) and second derivative (T<sub>j,tl</sub>
) of the transformation matrix.</p>
<p>The mass matrix is also defined:</p>
<p><disp-formula id="bmcM9"><label>(9)</label>
<mml:math id="M42" name="1743-0003-10-23-i43" overflow="scroll"><mml:mrow><mml:mi>M</mml:mi>
<mml:mo>=</mml:mo>
<mml:mtext mathvariant="italic">diag</mml:mtext>
<mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>m</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>m</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>m</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>I</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>m</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>m</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>I</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>With m<sub>f</sub>
, m<sub>1</sub>
, m<sub>2</sub>
, I<sub>1</sub>
and I<sub>2</sub>
, the mass of the foot, the segments masses and inertia’s respectively.</p>
<p>And we also define the external forces and moments (<italic>f</italic>
): this is the sum of the gravitational forces F<sub>grav, i</sub>
, the external forces F<sub>ext, i ,</sub>
and the joint torques (τ<sub>joint,i</sub>
):</p>
<p><disp-formula id="bmcM10"><label>(10)</label>
<mml:math id="M43" name="1743-0003-10-23-i44" overflow="scroll"><mml:mtable columnalign="center"><mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mi>f</mml:mi>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mo>=</mml:mo>
</mml:mtd>
<mml:mtd columnalign="center"><mml:msub><mml:mi mathvariant="normal">F</mml:mi>
<mml:mrow><mml:mtext mathvariant="normal">grav</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mo>+</mml:mo>
</mml:mtd>
<mml:mtd columnalign="center"><mml:msub><mml:mi mathvariant="normal">F</mml:mi>
<mml:mrow><mml:mtext mathvariant="normal">ext</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mo>+</mml:mo>
</mml:mtd>
<mml:mtd columnalign="center"><mml:msub><mml:mi mathvariant="normal">τ</mml:mi>
<mml:mrow><mml:mtext mathvariant="normal">joint</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mfenced open="[" close="]"><mml:mtable columnalign="center"><mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:msub><mml:mi>F</mml:mi>
<mml:mi mathvariant="italic">sx</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mrow><mml:mo>-</mml:mo>
<mml:msub><mml:mi>m</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>*</mml:mo>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mrow><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>τ</mml:mi>
<mml:mtext mathvariant="italic">ank</mml:mtext>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub><mml:mi>τ</mml:mi>
<mml:mtext mathvariant="italic">hip</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>τ</mml:mi>
<mml:mrow><mml:mi>p</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mrow><mml:mo>-</mml:mo>
<mml:msub><mml:mi>m</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>*</mml:mo>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:msub><mml:mi>τ</mml:mi>
<mml:mtext mathvariant="italic">hip</mml:mtext>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mfenced>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mo>=</mml:mo>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mfenced open="[" close="]"><mml:mtable columnalign="center"><mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mrow><mml:mo>-</mml:mo>
<mml:msub><mml:mi>m</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>*</mml:mo>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mrow><mml:mo>-</mml:mo>
<mml:msub><mml:mi>m</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>*</mml:mo>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mfenced>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mo>+</mml:mo>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mfenced open="[" close="]"><mml:mtable columnalign="center"><mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:msub><mml:mi>F</mml:mi>
<mml:mi mathvariant="italic">sx</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mrow><mml:mo>-</mml:mo>
<mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mi>p</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>*</mml:mo>
<mml:msub><mml:mi>p</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>*</mml:mo>
<mml:mtext>sin</mml:mtext>
<mml:mfenced open="(" close=")"><mml:msub><mml:mi>θ</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mfenced>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mfenced>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mo>+</mml:mo>
</mml:mtd>
<mml:mtd columnalign="center"><mml:mfenced open="[" close="]"><mml:mtable columnalign="center"><mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>τ</mml:mi>
<mml:mtext mathvariant="italic">ank</mml:mtext>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub><mml:mi>τ</mml:mi>
<mml:mtext mathvariant="italic">hip</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:msub><mml:mi>τ</mml:mi>
<mml:mtext mathvariant="italic">hip</mml:mtext>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mfenced>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
</p>
<fig id="F8" position="float"><label>Figure 8</label>
<caption><p><bold>Model of human balance control. </bold>
The model consists of body mechanics and a controller with intrinsic stiffness and damping (C<sub>pas</sub>
), an active proportional derivative controller (C<sub>act</sub>
), timedelays (δ), muscle activation dynamics (H<sub>act</sub>
), sensor and measurement noise. T<sub>ank</sub>
and T<sub>hip </sub>
denote the respective joint torques, θ<sub>1 </sub>
and θ<sub>2 </sub>
the joint angles. S<sub>x </sub>
and F<sub>pert </sub>
are the force and platform perturbations, respectively. g is the gravitational acceleration.</p>
</caption>
<graphic xlink:href="1743-0003-10-23-8"></graphic>
</fig>
<p>With g the gravitational acceleration (9.81m/s<sup>2</sup>
), <italic>τ</italic>
<sub><italic>p</italic>
1</sub>
= - <italic>Fp</italic>
<sub>1</sub>
* <italic>p</italic>
<sub>1</sub>
* sin(<italic>θ</italic>
<sub>1</sub>
), the perturbation torques due to the external perturbation force. Note that joint moments (τ<sub>ank</sub>
and τ<sub>hip</sub>
) are the inputs of the model (forward simulation). Note, that the internal forces are not incorporated in the force matrix; these are implicitly incorporated in the T<sub>j,t</sub>
and T<sub>j,tl</sub>
matrices.</p>
<p>The movement equations with the TMT method lead to the following expression:</p>
<p><disp-formula id="bmcM11"><label>(11)</label>
<mml:math id="M44" name="1743-0003-10-23-i45" overflow="scroll"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>M</mml:mi>
<mml:mo stretchy="true">¯</mml:mo>
</mml:mover>
<mml:mrow><mml:mtext mathvariant="italic">red</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">¨</mml:mo>
</mml:mover>
<mml:mo>=</mml:mo>
<mml:msub><mml:mover accent="true"><mml:mi>f</mml:mi>
<mml:mo stretchy="true">¯</mml:mo>
</mml:mover>
<mml:mrow><mml:mtext mathvariant="italic">red</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>With <inline-formula><mml:math id="M45" name="1743-0003-10-23-i46" overflow="scroll"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>M</mml:mi>
<mml:mo stretchy="true">¯</mml:mo>
</mml:mover>
<mml:mrow><mml:mtext mathvariant="italic">red</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msup><mml:mi>T</mml:mi>
<mml:mi>T</mml:mi>
</mml:msup>
<mml:mi mathvariant="italic">MT</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>
and T = (T<sub>i,t</sub>
), is the first derivative of the transformation matrix.</p>
<p><disp-formula id="bmcM12"><label>(12)</label>
<mml:math id="M46" name="1743-0003-10-23-i47" overflow="scroll"><mml:mrow><mml:mtext mathvariant="normal">And</mml:mtext>
<mml:mspace width="0.24em"></mml:mspace>
<mml:msub><mml:mover accent="true"><mml:mi>f</mml:mi>
<mml:mo stretchy="true">¯</mml:mo>
</mml:mover>
<mml:mrow><mml:mtext mathvariant="italic">red</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msup><mml:mi>T</mml:mi>
<mml:mi>T</mml:mi>
</mml:msup>
<mml:mfenced open="[" close="]"><mml:mstyle displaystyle="true"><mml:mo mathsize="big">∑</mml:mo>
<mml:mrow><mml:mi>f</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi mathvariant="italic">Mg</mml:mi>
</mml:mrow>
</mml:mstyle>
</mml:mfenced>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>With f the external forces, M the mass matrix and <inline-formula><mml:math id="M47" name="1743-0003-10-23-i48" overflow="scroll"><mml:mrow><mml:mi>g</mml:mi>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">tm</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
; this term corresponds to the centripetal term of the movement equations. In the case of forward simulation, we want to know the accelerations and the torques that are the inputs of the model:</p>
<p><disp-formula id="bmcM13"><label>(13)</label>
<mml:math id="M48" name="1743-0003-10-23-i49" overflow="scroll"><mml:mrow><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">¨</mml:mo>
</mml:mover>
<mml:mo>=</mml:mo>
<mml:msup><mml:mover accent="true"><mml:mi>M</mml:mi>
<mml:mo stretchy="true">¯</mml:mo>
</mml:mover>
<mml:mrow><mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mover accent="true"><mml:mi>f</mml:mi>
<mml:mo stretchy="true">¯</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>It can be deduced from Equation 11 and 12 that the equations of motion can be rewritten as:</p>
<p><disp-formula id="bmcM14"><label>(14)</label>
<mml:math id="M49" name="1743-0003-10-23-i50" overflow="scroll"><mml:mrow><mml:msub><mml:mi>M</mml:mi>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">¨</mml:mo>
</mml:mover>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">tm</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">grav</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">ext</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>τ</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">joint</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>For which we have used the index notation.</p>
</sec>
<sec><title>Drivers</title>
<p>In the case of a platform disturbance and a perturbation force, one of the generalized coordinates is specified, namely the platform movement (note that platform accelerations generate the perturbation). Therefore, the platform movement (<italic>S</italic>
<sub><italic>x</italic>
</sub>
) can be considered as a known degree of freedom, while the remaining degrees of freedom are unknown:</p>
<p><disp-formula id="bmcM15"><label>(15)</label>
<mml:math id="M50" name="1743-0003-10-23-i51" overflow="scroll"><mml:mrow><mml:msub><mml:mi>q</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfenced open="[" close="]"><mml:mtable columnalign="center"><mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:msub><mml:mi>q</mml:mi>
<mml:mi>u</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:msub><mml:mi>q</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mfenced>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>This results in a split up of the T<sub>j,t</sub>
matrix in a part for the known (T<sub>j,k</sub>
) and unknown degrees of freedom (T<sub>j,u</sub>
) and differentiation of these matrices results in <inline-formula><mml:math id="M51" name="1743-0003-10-23-i52" overflow="scroll"><mml:mrow><mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">kp</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
and <inline-formula><mml:math id="M52" name="1743-0003-10-23-i53" overflow="scroll"><mml:mrow><mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">un</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>u</mml:mi>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>n</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
, respectively. Note that the subscripts u and k denote the unknown and known coordinates, respectively.</p>
<p>By substitution of these matrices in Equation 14, we get:</p>
<p><disp-formula id="bmcM16"><label>(16)</label>
<mml:math id="M53" name="1743-0003-10-23-i54" overflow="scroll"><mml:mrow><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mi>M</mml:mi>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">¨</mml:mo>
</mml:mover>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">¨</mml:mo>
</mml:mover>
<mml:mi>u</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">un</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>u</mml:mi>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">kp</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mspace width="1em"></mml:mspace>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">grav</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">exter</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>τ</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">joint</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p><disp-formula><mml:math id="M54" name="1743-0003-10-23-i55" overflow="scroll"><mml:mrow><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msubsup><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:msubsup>
<mml:msub><mml:mi>M</mml:mi>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:msub>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">¨</mml:mo>
</mml:mover>
<mml:mi>u</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
</mml:mtd>
<mml:mtd><mml:msubsup><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:msubsup>
<mml:mo stretchy="true">(</mml:mo>
<mml:mrow><mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">grav</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">ext</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>τ</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">joint</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub><mml:mi>M</mml:mi>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:msub>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">¨</mml:mo>
</mml:mover>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd></mml:mtd>
<mml:mtd><mml:mrow><mml:mo>-</mml:mo>
<mml:msub><mml:mi>M</mml:mi>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:msub>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">un</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>u</mml:mi>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub><mml:mi>M</mml:mi>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:msub>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">kp</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">)</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p><disp-formula><mml:math id="M55" name="1743-0003-10-23-i56" overflow="scroll"><mml:mtable columnalign="left"><mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">¨</mml:mo>
</mml:mover>
<mml:mi>u</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mrow><mml:msubsup><mml:mi>M</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">red</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow><mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">red</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:msub><mml:mi>M</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">red</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mrow><mml:msubsup><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:msubsup>
<mml:msub><mml:mi>M</mml:mi>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:msub>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="left"><mml:mtd columnalign="left"><mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">red</mml:mtext>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mrow><mml:msubsup><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:msubsup>
<mml:mo stretchy="true">(</mml:mo>
<mml:mrow><mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">grav</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>F</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">ext</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>τ</mml:mi>
<mml:mrow><mml:mtext mathvariant="italic">joint</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub><mml:mi>M</mml:mi>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:msub>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">¨</mml:mo>
</mml:mover>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mspace width="7em"></mml:mspace>
<mml:mrow><mml:mrow><mml:mo>-</mml:mo>
<mml:msub><mml:mi>M</mml:mi>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:msub>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">un</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>u</mml:mi>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub><mml:mi>M</mml:mi>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:msub>
<mml:msub><mml:mi>T</mml:mi>
<mml:mrow><mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">kp</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:msub><mml:mover accent="true"><mml:mi>q</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
</p>
</sec>
<sec><title>Linearization and state space notation</title>
<p>Subsequently, the equations of motion are linearized by differentiating the equations of motion to the states of the system <inline-formula><mml:math id="M56" name="1743-0003-10-23-i57" overflow="scroll"><mml:mrow><mml:mo stretchy="true">(</mml:mo>
<mml:msubsup><mml:mi>θ</mml:mi>
<mml:mn>1</mml:mn>
<mml:mn>0</mml:mn>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:msubsup><mml:mover accent="true"><mml:mi>θ</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mn>1</mml:mn>
<mml:mn>0</mml:mn>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:msubsup><mml:mi>θ</mml:mi>
<mml:mn>2</mml:mn>
<mml:mn>0</mml:mn>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:msubsup><mml:mover accent="true"><mml:mi>θ</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mn>2</mml:mn>
<mml:mn>0</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>
; the joint angles and joint angle velocity) and to the system inputs and external disturbances <inline-formula><mml:math id="M57" name="1743-0003-10-23-i58" overflow="scroll"><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>τ</mml:mi>
<mml:mtext mathvariant="italic">ank</mml:mtext>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>τ</mml:mi>
<mml:mtext mathvariant="italic">hip</mml:mtext>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mover accent="true"><mml:mi>S</mml:mi>
<mml:mo stretchy="true">¨</mml:mo>
</mml:mover>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mfenced>
</mml:math>
</inline-formula>
with a first order Taylor approximation:</p>
<p><disp-formula id="bmcM17"><label>(17)</label>
<mml:math id="M58" name="1743-0003-10-23-i59" overflow="scroll"><mml:mrow><mml:msub><mml:mi>F</mml:mi>
<mml:mtext mathvariant="italic">lin</mml:mtext>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mi>f</mml:mi>
<mml:mfenced open="(" close=")"><mml:mi>a</mml:mi>
</mml:mfenced>
<mml:mo>+</mml:mo>
<mml:msup><mml:mi>f</mml:mi>
<mml:mo>′</mml:mo>
</mml:msup>
<mml:mfenced open="(" close=")"><mml:mi>a</mml:mi>
</mml:mfenced>
<mml:mo>*</mml:mo>
<mml:mfenced open="(" close=")"><mml:mrow><mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>Where a is the equilibrium point and x the deviation from the operating point. In our case, the equilibrium position is straight stance (segment angles are 90°, angular velocities, external perturbations, and joint torques are zero). Note that the centripetal term of the movement equations disappears with linearization.</p>
<p>Finally, the equations are rewritten in a state-space notation:</p>
<p><disp-formula id="bmcM18"><label>(18)</label>
<mml:math id="M59" name="1743-0003-10-23-i60" overflow="scroll"><mml:mtable columnalign="center"><mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mrow><mml:mover accent="true"><mml:mi>x</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>A</mml:mi>
<mml:mi mathvariant="italic">ss</mml:mi>
</mml:msub>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>B</mml:mi>
<mml:mi mathvariant="italic">ss</mml:mi>
</mml:msub>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr columnalign="center"><mml:mtd columnalign="center"><mml:mrow><mml:mi>y</mml:mi>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>C</mml:mi>
<mml:mi mathvariant="italic">ss</mml:mi>
</mml:msub>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>D</mml:mi>
<mml:mi mathvariant="italic">ss</mml:mi>
</mml:msub>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
</p>
<p>With <inline-formula><mml:math id="M60" name="1743-0003-10-23-i61" overflow="scroll"><mml:mrow><mml:mover accent="true"><mml:mi>x</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mo>=</mml:mo>
<mml:msub><mml:mover accent="true"><mml:mi>θ</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mover accent="true"><mml:mi>θ</mml:mi>
<mml:mo stretchy="true">¨</mml:mo>
</mml:mover>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mover accent="true"><mml:mi>θ</mml:mi>
<mml:mo stretchy="true">˙</mml:mo>
</mml:mover>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mover accent="true"><mml:mi>θ</mml:mi>
<mml:mo stretchy="true">¨</mml:mo>
</mml:mover>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
(thus the derivative of the systems states) and <italic>x</italic>
the states. U are the inputs of the system and represent the joint torques, perturbation force, and the platform perturbation.</p>
</sec>
<sec><title>Stabilizing mechanism</title>
<p>Passive muscle stiffness and damping were added to each joint; hence, only mono-arcticular muscles are added. The built-up of muscle force is described by the muscle activation dynamics and these were modeled as a second-order dynamical system:</p>
<p><disp-formula id="bmcM19"><label>(19)</label>
<mml:math id="M61" name="1743-0003-10-23-i62" overflow="scroll"><mml:mrow><mml:msub><mml:mi>H</mml:mi>
<mml:mtext mathvariant="italic">act</mml:mtext>
</mml:msub>
<mml:mfenced open="(" close=")"><mml:mi>ω</mml:mi>
</mml:mfenced>
<mml:mo>=</mml:mo>
<mml:mfrac><mml:mn>1</mml:mn>
<mml:mrow><mml:mfrac><mml:mn>1</mml:mn>
<mml:msubsup><mml:mi>ω</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mfrac>
<mml:msup><mml:mi>s</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mfrac><mml:mrow><mml:mn>2</mml:mn>
<mml:mi>β</mml:mi>
</mml:mrow>
<mml:msub><mml:mi>ω</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:mi>s</mml:mi>
<mml:mo>+</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>With <italic>s</italic>
= jω, the natural frequency (<italic>ω</italic>
<sub>n</sub>
) was set at 13.8 rad/s (≈2.2 Hz) and the relative damping (<italic>β</italic>
) was 0.7. As the sensory signals do not reach the CNS instantaneously (because of neural conduction times) a time delay is also modeled as a pure transport delay, see also Figure 8. Parameter values were taken from the literature [<xref ref-type="bibr" rid="B47">47</xref>
,<xref ref-type="bibr" rid="B52">52</xref>
], see Table <xref ref-type="table" rid="T2">2</xref>
.</p>
<table-wrap position="float" id="T2"><label>Table 2</label>
<caption><p>Human balance control model parameters</p>
</caption>
<table frame="hsides" rules="groups" border="1"><colgroup><col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
</colgroup>
<thead valign="top"><tr><th colspan="2" align="center"><bold>Anthropomorphic</bold>
</th>
<th colspan="2" align="center"><bold>Controller properties</bold>
</th>
</tr>
</thead>
<tbody valign="top"><tr><td align="center" valign="bottom"><bold>Parameter</bold>
<hr></hr>
</td>
<td align="center" valign="bottom"><bold>Value</bold>
<hr></hr>
</td>
<td align="center" valign="bottom"><bold>Parameter</bold>
<hr></hr>
</td>
<td align="center" valign="bottom"><bold>Value</bold>
<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">Mass<hr></hr>
</td>
<td align="center" valign="bottom">70 (kg)<hr></hr>
</td>
<td align="center" valign="bottom">Ankle stiffness<hr></hr>
</td>
<td align="center" valign="bottom">293 (Nm/rad)<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">m<sub>1</sub>
<hr></hr>
</td>
<td align="center" valign="bottom">22.54 (kg)<hr></hr>
</td>
<td align="center" valign="bottom">Ankle damping<hr></hr>
</td>
<td align="center" valign="bottom">2.2 (Nm s/rad)<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">m<sub>2</sub>
<hr></hr>
</td>
<td align="center" valign="bottom">47.46 (kg)<hr></hr>
</td>
<td align="center" valign="bottom">Hip stiffness<hr></hr>
</td>
<td align="center" valign="bottom">95 (Nm/rad)<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">Length<hr></hr>
</td>
<td align="center" valign="bottom">1.70 (m)<hr></hr>
</td>
<td align="center" valign="bottom">Hip damping<hr></hr>
</td>
<td align="center" valign="bottom">27.4 (Nm s/rad)<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">l<sub>1</sub>
<hr></hr>
</td>
<td align="center" valign="bottom">1 (m)<hr></hr>
</td>
<td align="center" valign="bottom">Time delay ankle<hr></hr>
</td>
<td align="center" valign="bottom">60 (ms)<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">l<sub>2</sub>
<hr></hr>
</td>
<td align="center" valign="bottom">0.7 (m)<hr></hr>
</td>
<td align="center" valign="bottom">Time delay hip<hr></hr>
</td>
<td align="center" valign="bottom">40 (ms)<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">d<sub>1</sub>
<hr></hr>
</td>
<td align="center" valign="bottom">0.57 (m)<hr></hr>
</td>
<td colspan="2" align="center" valign="bottom"> <hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">d<sub>2</sub>
<hr></hr>
</td>
<td align="center" valign="bottom">0.26 (m)<hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
<td align="center" valign="bottom"> <hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">I<sub>1</sub>
<hr></hr>
</td>
<td align="center" valign="bottom">0.65 (m)<hr></hr>
</td>
<td colspan="2" align="center" valign="bottom"><italic>Muscle dynamics</italic>
<hr></hr>
</td>
</tr>
<tr><td align="center" valign="bottom">I<sub>2</sub>
<hr></hr>
</td>
<td align="center" valign="bottom">0.44 (m)<hr></hr>
</td>
<td align="center" valign="bottom">Eigen frequency (ω<sub>n</sub>
)<hr></hr>
</td>
<td align="center" valign="bottom">4.4π (rad/s)<hr></hr>
</td>
</tr>
<tr><td align="center">g</td>
<td align="center">9.81 (m/s<sup>2</sup>
)</td>
<td align="center">Damping (β)</td>
<td align="center">0.7 (-)</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p>Anthropomorphic data is taken from [<xref ref-type="bibr" rid="B52">52</xref>
]. The controller properties values are taken from [<xref ref-type="bibr" rid="B47">47</xref>
].</p>
</table-wrap-foot>
</table-wrap>
<p>Based on the A and B matrices of the state space equations and intrinsic joint stiffness and damping, the steady-state linear quadratic regulator (LQR), was used to uniquely determine the components of the optimal state feedback matrix C<sub>act</sub>
(see Figure 8). For the optimization, activation dynamics and timedelays are not included in the feedback pathway. The feedback parameters are obtained by minimizing the cost function <italic>J</italic>
of the form:</p>
<p><disp-formula id="bmcM20"><label>(20)</label>
<mml:math id="M62" name="1743-0003-10-23-i63" overflow="scroll"><mml:mrow><mml:mi>J</mml:mi>
<mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true"><mml:mrow><mml:mo mathsize="big">∫</mml:mo>
<mml:mrow><mml:mfenced open="(" close=")"><mml:mrow><mml:msup><mml:mi>x</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi mathvariant="italic">Qx</mml:mi>
<mml:mo>+</mml:mo>
<mml:msup><mml:mi>u</mml:mi>
<mml:mo>'</mml:mo>
</mml:msup>
<mml:mi mathvariant="italic">Ru</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>With <italic>x</italic>
the system’s states and <italic>u</italic>
the system’s input. Q and R are diagonal matrices; the elements in Q are set to 1 and in R set to 10<sup>6</sup>
.</p>
<p>As the sensory signals do not reach the CNS instantaneously (because of neural conduction times) a time delay is included as a pure transport delay, see Figure 8. Parameter values were taken from the literature [<xref ref-type="bibr" rid="B6">6</xref>
], see Table <xref ref-type="table" rid="T2">2</xref>
.</p>
<p>Pink sensor noise was created off line, by scaling the power spectrum of a random timeseries by 1.2 such that at 1Hz, the power of the signal was 1.5*10<sup>-7</sup>
. Subsequently, this signal was added to the joint angles in the model simulations; this lead to a spontaneous sway of around 0.6° (peak-to-peak amplitude; comparable with human quiet stance data). Measurement noise was modeled as white noise (zero mean and 0.0001 variance) and added to the joint angles and joint torques.</p>
</sec>
</sec>
<sec><title>Abbreviations</title>
<p>CNS: Central nervous system; ARMAX: Auto regressive moving average eXogenous; FRF: Frequency response function; EMG: ElectroMyoGram; PD: Parkinson’s disease; SISO: Single input single output; MIMO: Multiple input multiple output; HAT: Head arms trunk; DoF: Degrees of freedom; CoM: Center of mass; GOF: Goodness of fit; NSR: Noise-to-signal ratio; LTI: Linear time invariant; USA: United States of America; m: Meter; Nm: Newton/m; Hz: Hertz; std: Standard deviation; CLSIT: Closed-loop-system-identification-technique; yrs: Years.</p>
</sec>
<sec><title>Competing interests</title>
<p>The authors declare no competing interests.</p>
</sec>
<sec><title>Authors’ contributions</title>
<p>TB performed the model simulations, participated in the conception and design of the experiment, performed the data acquisition, analyzed the data, drafted and revised the manuscript. AS has made a substantial contribution drafting and critically revising the manuscript. HvdK conceived of the study and participated in the design and the coordination of the study and critically revised the manuscript. All authors have read and approved the final manuscript.</p>
</sec>
<sec><title>Authors’ information</title>
<p>TB received the M.Sc. degree in human movement sciences of the Free University (Amsterdam, The Netherlands) in 2006. She is currently working towards her PhD degree at the laboratory of biomechanical engineering at the University of Twente, in close collaboration with the Radboud University medical Centre in Nijmegen. Her research focuses on developing new techniques and methods based on system identification to investigate and quantify human balance control. Concurrently, she is applying these new and adopted methods for patients with Parkinson’s disease in order to create more insight into the pathophysiology of this disease.</p>
<p>AC received the M.Sc. and Ph.D. degrees in mechanical engineering from the Delft University of Technology, The Netherlands, in 1999 and 2004, respectively. He holds a position as an assistant professor at the Delft University of Technology. He is co-founder of the Delft Laboratory for Neuromuscular Control. His research interest is the field of neuromuscular control and includes techniques to quantify the functional contribution of afferent feedback, neuromuscular modeling, haptic manipulators, and system identification. His research focuses on both able-bodied individuals and individuals suffering from movement disorders.</p>
<p>HvdK (1970, Rotterdam, the Netherlands) received his Phd with honors (cum laude) in 2000 and is professor in Biomechatronics and Rehabilitation Technology at the Department of Biomechanical Engineering at the University of Twente (0.8 fte), and Delft University of Technology (0.2 fte), the Netherlands. His expertise and interests are in the field of human motor control, adaptation, and learning, rehabilitation robots, diagnostic, and assistive robotics, virtual reality, rehabilitation medicine, and neuro computational modeling. Dr.ir. H. van der Kooij has published over 60 papers in the area of biomechatronics and human motor control, and has directed approximately € 6 million in research funding over the past 10 years. He is associate editor of IEEE TBME, member of IEEE EMBS technical committee of Biorobots and was member of several scientific program committees in the field of rehabilitation robotics, bio robotics, and assistive devices. At the UT he is founder and head of Rehabilitation robotics laboratory that developed powered exoskeletons for the rehabilitation of upper and lower extremities. He is founder and head of the Virtual Reality Human performance lab that combines robotic devices, motion capturing and virtual environments to asses and train human balance, walking, and hand-eye coordination.</p>
</sec>
</body>
<back><sec><title>Acknowledgements</title>
<p>The authors gratefully acknowledge the support of the BrainGain Smart Mix Programme of the Netherlands Ministry of Economic Affairs and the Netherlands Ministry of Education, Culture and Science. We thank Jantsje Pasma and Denise Engelhart for their assistance with the measurements. We also want to acknowledge the constructive suggestions of the anonymous reviewers.</p>
</sec>
<ref-list><ref id="B1"><mixed-citation publication-type="journal"><name><surname>Peterka</surname>
<given-names>RJ</given-names>
</name>
<article-title>Sensorimotor integration in human postural control</article-title>
<source>J Neurophysiol</source>
<year>2002</year>
<volume>88</volume>
<issue>3</issue>
<fpage>1097</fpage>
<lpage>118</lpage>
<pub-id pub-id-type="pmid">12205132</pub-id>
</mixed-citation>
</ref>
<ref id="B2"><mixed-citation publication-type="journal"><name><surname>Kiemel</surname>
<given-names>T</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>Y</given-names>
</name>
<name><surname>Jeka</surname>
<given-names>JJ</given-names>
</name>
<article-title>Identification of neural feedback for upright stance in humans: stabilization rather than sway minimization</article-title>
<source>J Neurosci</source>
<year>2011</year>
<volume>31</volume>
<issue>42</issue>
<fpage>15144</fpage>
<lpage>53</lpage>
<pub-id pub-id-type="doi">10.1523/JNEUROSCI.1013-11.2011</pub-id>
<pub-id pub-id-type="pmid">22016548</pub-id>
</mixed-citation>
</ref>
<ref id="B3"><mixed-citation publication-type="journal"><name><surname>Peterka</surname>
<given-names>RJ</given-names>
</name>
<article-title>Simplifying the complexities of maintaining balance</article-title>
<source>IEEE Eng Med Biol Mag</source>
<year>2003</year>
<volume>22</volume>
<issue>2</issue>
<fpage>63</fpage>
<lpage>8</lpage>
<pub-id pub-id-type="doi">10.1109/MEMB.2003.1195698</pub-id>
<pub-id pub-id-type="pmid">12733461</pub-id>
</mixed-citation>
</ref>
<ref id="B4"><mixed-citation publication-type="journal"><name><surname>van Asseldonk</surname>
<given-names>EH</given-names>
</name>
<article-title>Disentangling the contribution of the paretic and non-paretic ankle to balance control in stroke patients</article-title>
<source>Exp Neurol</source>
<year>2006</year>
<volume>201</volume>
<issue>2</issue>
<fpage>441</fpage>
<lpage>51</lpage>
<pub-id pub-id-type="doi">10.1016/j.expneurol.2006.04.036</pub-id>
<pub-id pub-id-type="pmid">16814283</pub-id>
</mixed-citation>
</ref>
<ref id="B5"><mixed-citation publication-type="journal"><name><surname>van Soest</surname>
<given-names>AJ</given-names>
</name>
<name><surname>Rozendaal</surname>
<given-names>LA</given-names>
</name>
<article-title>The inverted pendulum model of bipedal standing cannot be stabilized through direct feedback of force and contractile element length and velocity at realistic series elastic element stiffness</article-title>
<source>Biol Cybern</source>
<year>2008</year>
<volume>99</volume>
<issue>1</issue>
<fpage>29</fpage>
<lpage>41</lpage>
<pub-id pub-id-type="doi">10.1007/s00422-008-0240-2</pub-id>
<pub-id pub-id-type="pmid">18584202</pub-id>
</mixed-citation>
</ref>
<ref id="B6"><mixed-citation publication-type="journal"><name><surname>Horlings</surname>
<given-names>CG</given-names>
</name>
<article-title>Vestibular and proprioceptive influences on trunk movements during quiet standing</article-title>
<source>Neuroscience</source>
<year>2009</year>
<volume>161</volume>
<issue>3</issue>
<fpage>904</fpage>
<lpage>14</lpage>
<pub-id pub-id-type="doi">10.1016/j.neuroscience.2009.04.005</pub-id>
<pub-id pub-id-type="pmid">19362125</pub-id>
</mixed-citation>
</ref>
<ref id="B7"><mixed-citation publication-type="journal"><name><surname>Creath</surname>
<given-names>R</given-names>
</name>
<article-title>A unified view of quiet and perturbed stance: simultaneous co-existing excitable modes</article-title>
<source>Neurosci Lett</source>
<year>2005</year>
<volume>377</volume>
<issue>2</issue>
<fpage>75</fpage>
<lpage>80</lpage>
<pub-id pub-id-type="doi">10.1016/j.neulet.2004.11.071</pub-id>
<pub-id pub-id-type="pmid">15740840</pub-id>
</mixed-citation>
</ref>
<ref id="B8"><mixed-citation publication-type="journal"><name><surname>Pinter</surname>
<given-names>IJ</given-names>
</name>
<article-title>The dynamics of postural sway cannot be captured using a one-segment inverted pendulum model: a PCA on segment rotations during unperturbed stance</article-title>
<source>J Neurophysiol</source>
<year>2008</year>
<volume>100</volume>
<issue>6</issue>
<fpage>3197</fpage>
<lpage>208</lpage>
<pub-id pub-id-type="doi">10.1152/jn.01312.2007</pub-id>
<pub-id pub-id-type="pmid">18829852</pub-id>
</mixed-citation>
</ref>
<ref id="B9"><mixed-citation publication-type="journal"><name><surname>Horak</surname>
<given-names>FB</given-names>
</name>
<name><surname>Nashner</surname>
<given-names>LM</given-names>
</name>
<article-title>Central programming of postural movements: adaptation to altered support-surface configurations</article-title>
<source>J Neurophysiol</source>
<year>1986</year>
<volume>55</volume>
<issue>6</issue>
<fpage>1369</fpage>
<lpage>81</lpage>
<pub-id pub-id-type="pmid">3734861</pub-id>
</mixed-citation>
</ref>
<ref id="B10"><mixed-citation publication-type="book"><name><surname>Horak</surname>
<given-names>FB</given-names>
</name>
<name><surname>Macpherson</surname>
<given-names>JM</given-names>
</name>
<person-group person-group-type="editor">Shepard JT, Rowell LB</person-group>
<article-title><italic>Postural orientation and equilibrium</italic>
</article-title>
<source>Handbook of Physiology: Section 12, Exercise: Regulation and Integration of Multiple Systems</source>
<year>1996</year>
<publisher-name>New York: Oxford University Press</publisher-name>
<fpage>255</fpage>
<lpage>292</lpage>
</mixed-citation>
</ref>
<ref id="B11"><mixed-citation publication-type="journal"><name><surname>Fitzpatrick</surname>
<given-names>R</given-names>
</name>
<name><surname>Burke</surname>
<given-names>D</given-names>
</name>
<name><surname>Gandevia</surname>
<given-names>SC</given-names>
</name>
<article-title>Loop gain of reflexes controlling human standing measured with the use of postural and vestibular disturbances</article-title>
<source>J Neurophysiol</source>
<year>1996</year>
<volume>76</volume>
<issue>6</issue>
<fpage>3994</fpage>
<lpage>4008</lpage>
<pub-id pub-id-type="pmid">8985895</pub-id>
</mixed-citation>
</ref>
<ref id="B12"><mixed-citation publication-type="journal"><name><surname>van der Kooij</surname>
<given-names>H</given-names>
</name>
<name><surname>van Asseldonk</surname>
<given-names>E</given-names>
</name>
<name><surname>van der Helm</surname>
<given-names>FC</given-names>
</name>
<article-title>Comparison of different methods to identify and quantify balance control</article-title>
<source>J Neurosci Methods</source>
<year>2005</year>
<volume>145</volume>
<issue>1–2</issue>
<fpage>175</fpage>
<lpage>203</lpage>
<pub-id pub-id-type="pmid">15922036</pub-id>
</mixed-citation>
</ref>
<ref id="B13"><mixed-citation publication-type="book"><name><surname>Pintelon</surname>
<given-names>R</given-names>
</name>
<name><surname>Schoukens</surname>
<given-names>J</given-names>
</name>
<source>System identification, a frequency domain approach. Vol. 605</source>
<year>2001</year>
<publisher-name>New York: IEEE Press</publisher-name>
</mixed-citation>
</ref>
<ref id="B14"><mixed-citation publication-type="journal"><name><surname>de Vlugt</surname>
<given-names>E</given-names>
</name>
<name><surname>Schouten</surname>
<given-names>AC</given-names>
</name>
<name><surname>van der Helm</surname>
<given-names>FC</given-names>
</name>
<article-title>Closed-loop multivariable system identification for the characterization of the dynamic arm compliance using continuous force disturbances: a model study</article-title>
<source>J Neurosci Methods</source>
<year>2003</year>
<volume>122</volume>
<issue>2</issue>
<fpage>123</fpage>
<lpage>40</lpage>
<pub-id pub-id-type="doi">10.1016/S0165-0270(02)00303-5</pub-id>
<pub-id pub-id-type="pmid">12573472</pub-id>
</mixed-citation>
</ref>
<ref id="B15"><mixed-citation publication-type="journal"><name><surname>Kuo</surname>
<given-names>AD</given-names>
</name>
<article-title>Effect of altered sensory conditions on multivariate descriptors of human postural sway</article-title>
<source>Exp Brain Res</source>
<year>1998</year>
<volume>122</volume>
<issue>2</issue>
<fpage>185</fpage>
<lpage>95</lpage>
<pub-id pub-id-type="doi">10.1007/s002210050506</pub-id>
<pub-id pub-id-type="pmid">9776517</pub-id>
</mixed-citation>
</ref>
<ref id="B16"><mixed-citation publication-type="journal"><name><surname>Speers</surname>
<given-names>RA</given-names>
</name>
<name><surname>Kuo</surname>
<given-names>AD</given-names>
</name>
<name><surname>Horak</surname>
<given-names>FB</given-names>
</name>
<article-title>Contributions of altered sensation and feedback responses to changes in coordination of postural control due to aging</article-title>
<source>Gait Posture</source>
<year>2002</year>
<volume>16</volume>
<issue>1</issue>
<fpage>20</fpage>
<lpage>30</lpage>
<pub-id pub-id-type="doi">10.1016/S0966-6362(02)00003-6</pub-id>
<pub-id pub-id-type="pmid">12127183</pub-id>
</mixed-citation>
</ref>
<ref id="B17"><mixed-citation publication-type="journal"><name><surname>Alexandrov</surname>
<given-names>AV</given-names>
</name>
<name><surname>Frolov</surname>
<given-names>AA</given-names>
</name>
<name><surname>Massion</surname>
<given-names>J</given-names>
</name>
<article-title><italic>Biomechanical analysis of movement strategies in human forward trunk bending.</italic>
II. Experimental study</article-title>
<source>Biol Cybern</source>
<year>2001</year>
<volume>84</volume>
<issue>6</issue>
<fpage>435</fpage>
<lpage>43</lpage>
<pub-id pub-id-type="doi">10.1007/PL00007987</pub-id>
<pub-id pub-id-type="pmid">11417055</pub-id>
</mixed-citation>
</ref>
<ref id="B18"><mixed-citation publication-type="journal"><name><surname>Alexandrov</surname>
<given-names>AV</given-names>
</name>
<article-title>Feedback equilibrium control during human standing</article-title>
<source>Biol Cybern</source>
<year>2005</year>
<volume>93</volume>
<issue>5</issue>
<fpage>309</fpage>
<lpage>22</lpage>
<pub-id pub-id-type="doi">10.1007/s00422-005-0004-1</pub-id>
<pub-id pub-id-type="pmid">16228222</pub-id>
</mixed-citation>
</ref>
<ref id="B19"><mixed-citation publication-type="journal"><name><surname>Park</surname>
<given-names>S</given-names>
</name>
<name><surname>Horak</surname>
<given-names>FB</given-names>
</name>
<name><surname>Kuo</surname>
<given-names>AD</given-names>
</name>
<article-title>Postural feedback responses scale with biomechanical constraints in human standing</article-title>
<source>Exp Brain Res</source>
<year>2004</year>
<volume>154</volume>
<issue>4</issue>
<fpage>417</fpage>
<lpage>27</lpage>
<pub-id pub-id-type="doi">10.1007/s00221-003-1674-3</pub-id>
<pub-id pub-id-type="pmid">14618285</pub-id>
</mixed-citation>
</ref>
<ref id="B20"><mixed-citation publication-type="other"><name><surname>Kim</surname>
<given-names>S</given-names>
</name>
<article-title>Postural Feedback Scaling Deficits in Parkinson's disease</article-title>
<source>J Neurophysiol</source>
<year>2009</year>
</mixed-citation>
</ref>
<ref id="B21"><mixed-citation publication-type="journal"><name><surname>Fujisawa</surname>
<given-names>N</given-names>
</name>
<article-title>Human standing posture control system depending on adopted strategies</article-title>
<source>Med Biol Eng Comput</source>
<year>2005</year>
<volume>43</volume>
<issue>1</issue>
<fpage>107</fpage>
<lpage>14</lpage>
<pub-id pub-id-type="doi">10.1007/BF02345130</pub-id>
<pub-id pub-id-type="pmid">15742727</pub-id>
</mixed-citation>
</ref>
<ref id="B22"><mixed-citation publication-type="journal"><name><surname>Stolze</surname>
<given-names>H</given-names>
</name>
<article-title>Falls in frequent neurological diseases–prevalence, risk factors and aetiology</article-title>
<source>J Neurol</source>
<year>2004</year>
<volume>251</volume>
<issue>1</issue>
<fpage>79</fpage>
<lpage>84</lpage>
<pub-id pub-id-type="doi">10.1007/s00415-004-0276-8</pub-id>
<pub-id pub-id-type="pmid">14999493</pub-id>
</mixed-citation>
</ref>
<ref id="B23"><mixed-citation publication-type="journal"><name><surname>Pickering</surname>
<given-names>RM</given-names>
</name>
<article-title>A meta-analysis of six prospective studies of falling in Parkinson's disease</article-title>
<source>Mov Disord</source>
<year>2007</year>
<volume>22</volume>
<issue>13</issue>
<fpage>1892</fpage>
<lpage>900</lpage>
<pub-id pub-id-type="doi">10.1002/mds.21598</pub-id>
<pub-id pub-id-type="pmid">17588236</pub-id>
</mixed-citation>
</ref>
<ref id="B24"><mixed-citation publication-type="journal"><name><surname>Hely</surname>
<given-names>MA</given-names>
</name>
<article-title>The Sydney multicenter study of Parkinson's disease: the inevitability of dementia at 20 years</article-title>
<source>Mov Disord</source>
<year>2008</year>
<volume>23</volume>
<issue>6</issue>
<fpage>837</fpage>
<lpage>44</lpage>
<pub-id pub-id-type="doi">10.1002/mds.21956</pub-id>
<pub-id pub-id-type="pmid">18307261</pub-id>
</mixed-citation>
</ref>
<ref id="B25"><mixed-citation publication-type="journal"><name><surname>Grimbergen</surname>
<given-names>YA</given-names>
</name>
<name><surname>Munneke</surname>
<given-names>M</given-names>
</name>
<name><surname>Bloem</surname>
<given-names>BR</given-names>
</name>
<article-title>Falls in Parkinson's disease</article-title>
<source>Curr Opin Neurol</source>
<year>2004</year>
<volume>17</volume>
<issue>4</issue>
<fpage>405</fpage>
<lpage>15</lpage>
<pub-id pub-id-type="doi">10.1097/01.wco.0000137530.68867.93</pub-id>
<pub-id pub-id-type="pmid">15247535</pub-id>
</mixed-citation>
</ref>
<ref id="B26"><mixed-citation publication-type="journal"><name><surname>Boonstra</surname>
<given-names>TA</given-names>
</name>
<article-title>Gait disorders and balance disturbances in Parkinson's disease: clinical update and pathophysiology</article-title>
<source>Curr Opin Neurol</source>
<year>2008</year>
<volume>24</volume>
<issue>4</issue>
<fpage>461</fpage>
<lpage>471</lpage>
<pub-id pub-id-type="doi">10.1097/WCO.0b013e328305bdaf</pub-id>
<pub-id pub-id-type="pmid">18607208</pub-id>
</mixed-citation>
</ref>
<ref id="B27"><mixed-citation publication-type="journal"><name><surname>Colnat-Coulbois</surname>
<given-names>S</given-names>
</name>
<article-title>Management of postural sensory conflict and dynamic balance control in late-stage Parkinson's disease</article-title>
<source>Neuroscience</source>
<year>2011</year>
<volume>193</volume>
<fpage>363</fpage>
<lpage>9</lpage>
<pub-id pub-id-type="pmid">21627979</pub-id>
</mixed-citation>
</ref>
<ref id="B28"><mixed-citation publication-type="journal"><name><surname>Carpenter</surname>
<given-names>MG</given-names>
</name>
<article-title>Postural abnormalities to multidirectional stance perturbations in Parkinson's disease</article-title>
<source>J Neurol Neurosurg Psychiatry</source>
<year>2004</year>
<volume>75</volume>
<issue>9</issue>
<fpage>1245</fpage>
<lpage>54</lpage>
<pub-id pub-id-type="doi">10.1136/jnnp.2003.021147</pub-id>
<pub-id pub-id-type="pmid">15314109</pub-id>
</mixed-citation>
</ref>
<ref id="B29"><mixed-citation publication-type="journal"><name><surname>Termoz</surname>
<given-names>N</given-names>
</name>
<article-title>The control of upright stance in young, elderly and persons with Parkinson's disease</article-title>
<source>Gait Posture</source>
<year>2008</year>
<volume>27</volume>
<issue>3</issue>
<fpage>463</fpage>
<lpage>70</lpage>
<pub-id pub-id-type="doi">10.1016/j.gaitpost.2007.05.015</pub-id>
<pub-id pub-id-type="pmid">17644337</pub-id>
</mixed-citation>
</ref>
<ref id="B30"><mixed-citation publication-type="journal"><name><surname>Maurer</surname>
<given-names>C</given-names>
</name>
<article-title>Effect of chronic bilateral subthalamic nucleus (STN) stimulation on postural control in Parkinson's disease</article-title>
<source>Brain</source>
<year>2003</year>
<volume>126</volume>
<issue>Pt 5</issue>
<fpage>1146</fpage>
<lpage>63</lpage>
<pub-id pub-id-type="pmid">12690054</pub-id>
</mixed-citation>
</ref>
<ref id="B31"><mixed-citation publication-type="journal"><name><surname>Djaldetti</surname>
<given-names>R</given-names>
</name>
<name><surname>Ziv</surname>
<given-names>I</given-names>
</name>
<name><surname>Melamed</surname>
<given-names>E</given-names>
</name>
<article-title>The mystery of motor asymmetry in Parkinson's disease</article-title>
<source>Lancet Neurol</source>
<year>2006</year>
<volume>5</volume>
<issue>9</issue>
<fpage>796</fpage>
<lpage>802</lpage>
<pub-id pub-id-type="doi">10.1016/S1474-4422(06)70549-X</pub-id>
<pub-id pub-id-type="pmid">16914408</pub-id>
</mixed-citation>
</ref>
<ref id="B32"><mixed-citation publication-type="journal"><name><surname>Geurts</surname>
<given-names>AC</given-names>
</name>
<article-title>Assessment of postural asymmetry in mild to moderate Parkinson's disease</article-title>
<source>Gait Posture</source>
<year>2011</year>
<volume>33</volume>
<issue>1</issue>
<fpage>143</fpage>
<lpage>5</lpage>
<pub-id pub-id-type="doi">10.1016/j.gaitpost.2010.09.018</pub-id>
<pub-id pub-id-type="pmid">20951590</pub-id>
</mixed-citation>
</ref>
<ref id="B33"><mixed-citation publication-type="journal"><name><surname>van der Kooij</surname>
<given-names>H</given-names>
</name>
<name><surname>Peterka</surname>
<given-names>RJ</given-names>
</name>
<article-title>Non-linear stimulus-response behavior of the human stance control system is predicted by optimization of a system with sensory and motor noise</article-title>
<source>J Comput Neurosci</source>
<year>2011</year>
<volume>30</volume>
<issue>3</issue>
<fpage>759</fpage>
<lpage>778</lpage>
<pub-id pub-id-type="doi">10.1007/s10827-010-0291-y</pub-id>
<pub-id pub-id-type="pmid">21161357</pub-id>
</mixed-citation>
</ref>
<ref id="B34"><mixed-citation publication-type="other"><name><surname>Robinson</surname>
<given-names>D</given-names>
</name>
<source>Series Elastic Actuator Development for a Biomimetic Robot, in IEEE/ASME International Conference on Advance Intelligent Mechantronics</source>
<year>1999</year>
</mixed-citation>
</ref>
<ref id="B35"><mixed-citation publication-type="journal"><name><surname>van der Kooij</surname>
<given-names>H</given-names>
</name>
<name><surname>de Vlugt</surname>
<given-names>E</given-names>
</name>
<article-title>Postural responses evoked by platform pertubations are dominated by continuous feedback</article-title>
<source>J Neurophysiol</source>
<year>2007</year>
<volume>98</volume>
<issue>2</issue>
<fpage>730</fpage>
<lpage>43</lpage>
<pub-id pub-id-type="doi">10.1152/jn.00457.2006</pub-id>
<pub-id pub-id-type="pmid">17460106</pub-id>
</mixed-citation>
</ref>
<ref id="B36"><mixed-citation publication-type="journal"><name><surname>Zajac</surname>
<given-names>FE</given-names>
</name>
<article-title>Understanding muscle coordination of the human leg with dynamical simulations</article-title>
<source>J Biomech</source>
<year>2002</year>
<volume>35</volume>
<issue>8</issue>
<fpage>1011</fpage>
<lpage>8</lpage>
<pub-id pub-id-type="doi">10.1016/S0021-9290(02)00046-5</pub-id>
<pub-id pub-id-type="pmid">12126660</pub-id>
</mixed-citation>
</ref>
<ref id="B37"><mixed-citation publication-type="journal"><name><surname>Nott</surname>
<given-names>CR</given-names>
</name>
<article-title>All joint moments significantly contribute to trunk angular acceleration</article-title>
<source>J Biomech</source>
<year>2010</year>
<volume>43</volume>
<issue>13</issue>
<fpage>2648</fpage>
<lpage>52</lpage>
<pub-id pub-id-type="doi">10.1016/j.jbiomech.2010.04.044</pub-id>
<pub-id pub-id-type="pmid">20646711</pub-id>
</mixed-citation>
</ref>
<ref id="B38"><mixed-citation publication-type="journal"><name><surname>Pasma</surname>
<given-names>JH</given-names>
</name>
<article-title>Sensory reweighting of proprioceptive information of the left and right leg during human balance control</article-title>
<source>J Neurophysiol</source>
<year>2012</year>
<volume>108</volume>
<issue>4</issue>
<fpage>1138</fpage>
<lpage>48</lpage>
<pub-id pub-id-type="doi">10.1152/jn.01008.2011</pub-id>
<pub-id pub-id-type="pmid">22623486</pub-id>
</mixed-citation>
</ref>
<ref id="B39"><mixed-citation publication-type="book"><name><surname>Koopman</surname>
<given-names>HFJM</given-names>
</name>
<source>The three-dimensional analysis and prediction of human walking</source>
<year>1989</year>
<publisher-name>The Netherlands: University of Twente</publisher-name>
</mixed-citation>
</ref>
<ref id="B40"><mixed-citation publication-type="journal"><name><surname>Koopman</surname>
<given-names>B</given-names>
</name>
<name><surname>Grootenboer</surname>
<given-names>HJ</given-names>
</name>
<name><surname>de Jongh</surname>
<given-names>HJ</given-names>
</name>
<article-title>An inverse dynamics model for the analysis, reconstruction and prediction of bipedal walking</article-title>
<source>J Biomech</source>
<year>1995</year>
<volume>28</volume>
<issue>11</issue>
<fpage>1369</fpage>
<lpage>76</lpage>
<pub-id pub-id-type="doi">10.1016/0021-9290(94)00185-7</pub-id>
<pub-id pub-id-type="pmid">8522549</pub-id>
</mixed-citation>
</ref>
<ref id="B41"><mixed-citation publication-type="journal"><name><surname>Brand</surname>
<given-names>RA</given-names>
</name>
<article-title>A model of lower extremity muscular anatomy</article-title>
<source>J Biomech Eng</source>
<year>1982</year>
<volume>104</volume>
<issue>4</issue>
<fpage>304</fpage>
<lpage>10</lpage>
<pub-id pub-id-type="doi">10.1115/1.3138363</pub-id>
<pub-id pub-id-type="pmid">7154650</pub-id>
</mixed-citation>
</ref>
<ref id="B42"><mixed-citation publication-type="other"><name><surname>Chandler</surname>
<given-names>RF</given-names>
</name>
<source>Investigation of the intertial properties of the human body</source>
<year>1975</year>
</mixed-citation>
</ref>
<ref id="B43"><mixed-citation publication-type="journal"><name><surname>Preuss</surname>
<given-names>R</given-names>
</name>
<name><surname>Fung</surname>
<given-names>J</given-names>
</name>
<article-title>A simple method to estimate force plate inertial components in a moving surface</article-title>
<source>J Biomech</source>
<year>2004</year>
<volume>37</volume>
<issue>8</issue>
<fpage>1177</fpage>
<lpage>80</lpage>
<pub-id pub-id-type="doi">10.1016/j.jbiomech.2003.12.007</pub-id>
<pub-id pub-id-type="pmid">15212922</pub-id>
</mixed-citation>
</ref>
<ref id="B44"><mixed-citation publication-type="book"><name><surname>Jenkins</surname>
<given-names>G</given-names>
</name>
<name><surname>Watts</surname>
<given-names>D</given-names>
</name>
<source>Spectral Analysis and Its Applications</source>
<year>1969</year>
<publisher-name>San Francisco: Holden Day</publisher-name>
</mixed-citation>
</ref>
<ref id="B45"><mixed-citation publication-type="journal"><name><surname>Kim</surname>
<given-names>S</given-names>
</name>
<name><surname>Atkeson</surname>
<given-names>CG</given-names>
</name>
<name><surname>Park</surname>
<given-names>S</given-names>
</name>
<article-title>Perturbation-dependent selection of postural feedback gain and its scaling</article-title>
<source>J Biomech</source>
<year>2012</year>
<volume>45</volume>
<issue>8</issue>
<fpage>1379</fpage>
<lpage>86</lpage>
<pub-id pub-id-type="doi">10.1016/j.jbiomech.2012.03.001</pub-id>
<pub-id pub-id-type="pmid">22444347</pub-id>
</mixed-citation>
</ref>
<ref id="B46"><mixed-citation publication-type="journal"><name><surname>de Vlugt</surname>
<given-names>E</given-names>
</name>
<name><surname>Schouten</surname>
<given-names>AC</given-names>
</name>
<name><surname>van der Helm</surname>
<given-names>FC</given-names>
</name>
<article-title>Quantification of intrinsic and reflexive properties during multijoint arm posture</article-title>
<source>J Neurosci Methods</source>
<year>2006</year>
<volume>155</volume>
<issue>2</issue>
<fpage>328</fpage>
<lpage>49</lpage>
<pub-id pub-id-type="doi">10.1016/j.jneumeth.2006.01.022</pub-id>
<pub-id pub-id-type="pmid">16504304</pub-id>
</mixed-citation>
</ref>
<ref id="B47"><mixed-citation publication-type="journal"><name><surname>Kiemel</surname>
<given-names>T</given-names>
</name>
<name><surname>Elahi</surname>
<given-names>AJ</given-names>
</name>
<name><surname>Jeka</surname>
<given-names>JJ</given-names>
</name>
<article-title>Identification of the plant for upright stance in humans: multiple movement patterns from a single neural strategy</article-title>
<source>J Neurophysiol</source>
<year>2008</year>
<volume>100</volume>
<issue>6</issue>
<fpage>3394</fpage>
<lpage>406</lpage>
<pub-id pub-id-type="doi">10.1152/jn.01272.2007</pub-id>
<pub-id pub-id-type="pmid">18829854</pub-id>
</mixed-citation>
</ref>
<ref id="B48"><mixed-citation publication-type="journal"><name><surname>Kiemel</surname>
<given-names>T</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>Y</given-names>
</name>
<name><surname>Jeka</surname>
<given-names>JJ</given-names>
</name>
<article-title>Visual flow is interpreted relative to multisegment postural control</article-title>
<source>J Mot Behav</source>
<year>2011</year>
<volume>43</volume>
<issue>3</issue>
<fpage>237</fpage>
<lpage>46</lpage>
<pub-id pub-id-type="doi">10.1080/00222895.2011.568991</pub-id>
<pub-id pub-id-type="pmid">21534025</pub-id>
</mixed-citation>
</ref>
<ref id="B49"><mixed-citation publication-type="journal"><name><surname>Rocchi</surname>
<given-names>L</given-names>
</name>
<name><surname>Chiari</surname>
<given-names>L</given-names>
</name>
<name><surname>Horak</surname>
<given-names>FB</given-names>
</name>
<article-title>Effects of deep brain stimulation and levodopa on postural sway in Parkinson's disease</article-title>
<source>J Neurol Neurosurg Psychiatry</source>
<year>2002</year>
<volume>73</volume>
<issue>3</issue>
<fpage>267</fpage>
<lpage>74</lpage>
<pub-id pub-id-type="doi">10.1136/jnnp.73.3.267</pub-id>
<pub-id pub-id-type="pmid">12185157</pub-id>
</mixed-citation>
</ref>
<ref id="B50"><mixed-citation publication-type="journal"><name><surname>van der Kooij</surname>
<given-names>H</given-names>
</name>
<article-title>Detecting asymmetries in balance control with system identification: first experimental results from Parkinson patients</article-title>
<source>J Neural Transm</source>
<year>2007</year>
<volume>114</volume>
<issue>10</issue>
<fpage>1333</fpage>
<lpage>7</lpage>
<pub-id pub-id-type="doi">10.1007/s00702-007-0801-x</pub-id>
<pub-id pub-id-type="pmid">17703275</pub-id>
</mixed-citation>
</ref>
<ref id="B51"><mixed-citation publication-type="journal"><name><surname>Roerdink</surname>
<given-names>M</given-names>
</name>
<article-title>On the relative contribution of the paretic leg to the control of posture after stroke</article-title>
<source>Neurorehabil Neural Repair</source>
<year>2009</year>
<volume>23</volume>
<issue>3</issue>
<fpage>267</fpage>
<lpage>74</lpage>
<pub-id pub-id-type="pmid">19074685</pub-id>
</mixed-citation>
</ref>
<ref id="B52"><mixed-citation publication-type="book"><name><surname>Winter</surname>
<given-names>DA</given-names>
</name>
<source>Biomechanics and motor control of human movement</source>
<year>1990</year>
<edition>2</edition>
<publisher-name>New York: John Wiley & Sons Inc</publisher-name>
</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Ticri/CIDE/explor/HapticV1/Data/Pmc/Curation
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 001260 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Curation/biblio.hfd -nk 001260 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Ticri/CIDE |area= HapticV1 |flux= Pmc |étape= Curation |type= RBID |clé= PMC:3662596 |texte= Identification of the contribution of the ankle and hip joints to multi-segmental balance control }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Curation/RBID.i -Sk "pubmed:23433148" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Curation/biblio.hfd \ | NlmPubMed2Wicri -a HapticV1
![]() | This area was generated with Dilib version V0.6.23. | ![]() |