ip model

%23%20%2F%2F%2F%20script%0A%23%20requires-python%20%3D%20%22%3E%3D3.11%22%0A%23%20dependencies%20%3D%20%5B%22marimo%3E%3D0.23.9%22%2C%20%22numpy%3D%3D2.4.3%22%2C%20%22matplotlib%3D%3D3.10.8%22%2C%20%22pandas%3D%3D3.0.2%22%2C%20%22scipy%3D%3D1.17.1%22%2C%20%22urllib3%3D%3D2.6.3%22%2C%20%22cogve%20%40%20..%2Fpublic%2Fwheels%2Fcogve-0.0.0%252Bmarimo.bc1946d696db-py3-none-any.whl%22%2C%20%22notebook_resources%20%40%20..%2Fpublic%2Fwheels%2Fnotebook_resources-0.0.0%252Bmarimo.9f9afe92b910-py3-none-any.whl%22%5D%0A%23%20%0A%23%20%5Btool.marimo.export%5D%0A%23%20lock_kind%20%3D%20%22observed%22%0A%23%20%2F%2F%2F%0A%0Aimport%20marimo%0A%0A__generated_with%20%3D%20%220.25.2%22%0Aapp%20%3D%20marimo.App()%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%20The%20inverted%20pendulum%20model%20of%20the%20human%20standing%0A%0A%20%20%20%20%3E%20Marcos%20Duarte%2C%20Renato%20Naville%20Watanabe%3Cbr%3E%0A%20%20%20%20%3E%20%5BLaboratory%20of%20Biomechanics%20and%20Motor%20Control%5D(https%3A%2F%2Fbmclab.pesquisa.ufabc.edu.br%2F)%3Cbr%3E%0A%20%20%20%20%3E%20Federal%20University%20of%20ABC%2C%20Brazil%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20How%20to%20use%20this%20guide%0A%0A%20%20%20%20We%20will%20derive%20the%20small-angle%20inverted-pendulum%20relation%20and%20use%20it%20to%20estimate%20the%20horizontal%20projection%20of%20center%20of%20mass%20from%20center%20of%20pressure.%20Keep%20the%20foot%20assumptions%2C%20height%20estimate%20and%20filtering%20interpretation%20explicit.%20The%20recorded%20COP%20is%20an%20experimental%20input%2C%20not%20a%20direct%20COM%20measurement.%0A%0A%20%20%20%20Try%20the%20distributed%20**Exercises**%20before%20opening%20their%20collapsed%20worked%20solutions.%20Then%20use%20the%20**Checkpoint%20questions**%20and%20attempt%20the%20**Problems**%20independently%3B%20their%20panels%20contain%20final%20answers%20only.%20Predict%20the%20effect%20of%20a%20parameter%20change%20before%20running%20it%2C%20and%20check%20units%20and%20limiting%20cases.%0A%0A%20%20%20%20**Opening%20challenge.**%20A%20force%20plate%20shows%20center%20of%20pressure%20moving%20rapidly%20while%20the%20body%20sways%20slowly.%20Are%20these%20two%20traces%20measuring%20the%20same%20physical%20point%3F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Python%20setup%0A%0A%20%20%20%20The%20dependency%20declaration%20at%20the%20top%20of%20this%20file%20lists%20the%20packages%20used%20by%20this%20guide.%20Marimo%20runs%20cells%20in%20dependency%20order.%20Run%20the%20notebook%20in%20an%20editor%20to%20use%20its%20displays%20and%20controls%3B%20the%20computational%20examples%20below%20explain%20their%20inputs%20before%20using%20them.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20Despite%20the%20enormous%20complexity%20of%20the%20human%20body%2C%20part%20of%20the%20mechanical%20behavior%20of%20the%20human%20body%20during%20the%20standing%20still%20posture%2C%20namely%20the%20displacements%20of%20the%20center%20of%20gravity%20(%24COG%24)%20and%20center%20of%20pressure%20(%24COP%24)%20in%20the%20anterior-posterior%20direction%2C%20can%20be%20elegantly%20portraied%20by%20a%20physical-mathematical%20model%20of%20an%20inverted%20pendulum%20with%20rigid%20segments%20articulated%20by%20joints.%0A%0A%20%20%20%20Using%20such%20a%20model%2C%20it's%20possible%20to%20estimate%20the%20COG%20vertical%20projection%20(COGv)%20from%20the%20COP%20displacement.%20The%20Python%20function%20%60cogve.py%60%20(code%20at%20the%20end%20of%20this%20text)%20performs%20this%20estimation.%20The%20function%20signature%20is%3A%0A%20%20%20%20%60%60%60python%0A%20%20%20%20cogv%20%3D%20cogve(COP%2C%20freq%2C%20mass%2C%20height%2C%20show%3DFalse%2C%20ax%3DNone)%0A%20%20%20%20%60%60%60%0A%20%20%20%20Let's%20now%20derive%20the%20inverted%20pendulum%20model%20of%20the%20human%20standing%20posture%20implemented%20in%20this%20function.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Derivation%20of%20the%20inverted%20pendulum%20model%0A%0A%20%20%20%20In%20the%20most%20simple%20version%20of%20the%20model%2C%20the%20human%20body%20in%20the%20sagittal%20plane%20is%20reduced%20to%20a%20two-link%20body%20with%20a%20single%20inverted%20pendulum%20articulated%20by%20only%20one%20joint%20(representing%20the%20feet%2C%20with%20the%20rest%20of%20the%20body%20articulated%20by%20the%20ankle%20joint).%20Let's%20deduce%20the%20equations%20for%20such%20inverted%20pendulum%20model%20as%20the%20representation%20at%20the%20sagital%20plane%20of%20the%20human%20standing%20still%20posture.%20The%20inverted%20pendulum%20model%20and%20the%20correspondent%20free-body%20diagrams%20(FBDs)%20are%20shown%20in%20Figure%201.%0A%0A%20%20%20%20%3Cdiv%3E%3Cfigure%20aria-describedby%3D%22figure-caption-1%22%3E%3Cimg%20src%3D%22https%3A%2F%2Fgithub.com%2FBMClab%2FBMC%2Fblob%2Fmaster%2Fimages%2Finvpendulum.png%3Fraw%3D1%22%20width%3D400%20alt%3D%22onelink%22%2F%3E%3C%2Ffigure%3E%3C%2Fdiv%3E%0A%0A%20%20%20%20%3Cspan%20id%3D%22figure-caption-1%22%20class%3D%22figure-caption%22%20style%3D%22display%3Ablock%3Btext-align%3Acenter%3Bfont-style%3Aitalic%22%3E%3Cb%3EFigure%201.%3C%2Fb%3E%20Model%20of%20a%20two-link%20inverted%20pendulum%20and%20the%20external%20forces%20acting%20on%20it%20for%20the%20representation%20at%20the%20sagital%20plane%20of%20the%20human%20standing%20still%20posture%20and%20the%20corresponding%20free-body%20diagrams.%20%24COG%24%3A%20body%20center%20of%20gravity%3B%20%24COG_v%24%3A%20%24COG%24%20vertical%20projection%20(at%20the%20horizontal%20plane)%20in%20relation%20to%20the%20ankle%20joint%3B%20%24COP%24%3A%20body%20center%20of%20pressure%20in%20relation%20to%20the%20ankle%20joint%3B%20%24GRF%24%3A%20ground%20reaction%20force%20(typically%20measured%20by%20a%20force%20plate)%3B%20%24%5Calpha%24%3A%20angle%20of%20the%20body%20in%20relation%20to%20the%20vertical%20direction%3B%20%24m%24%3A%20mass%20of%20the%20body%20minus%20feet%3B%20%24g%24%3A%20acceleration%20of%20gravity%3B%20%24F_a%24%20and%20%24T_a%24%3A%20resultant%20force%20and%20torque%20at%20the%20ankle%20joint%3B%20%24h%24%3A%20height%20of%20the%20%24COG%24%20in%20relation%20to%20the%20ankle%20joint%3B%20%24m_f%24%20and%20%24h_f%24%3A%20mass%20and%20height%20of%20the%20feet.%3C%2Fspan%3E%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20The%20equations%20of%20motion%20for%20each%20FBD%20of%20the%20feet%20and%20rest-of-body%20segments%20at%20the%20sagittal%20plane%20(%24xy%24%20plane)%20can%20be%20expressed%20in%20the%20form%20of%20the%20Newton-Euler%20equations.%3Cbr%3E%0A%20%20%20%20%3Cbr%3E%0A%20%20%20%20%2F%2F%2F%20admonition%20%7C%20The%20Newton-Euler%20equations%0A%0A%20%20%20%20The%20%3Ca%20href%3D%22https%3A%2F%2Fen.wikipedia.org%2Fwiki%2FNewton%25E2%2580%2593Euler_equations%22%3ENewton-Euler%20equations%3C%2Fa%3E%20are%20a%20formalism%20to%20describe%20the%20combined%20translational%20and%20rotational%20dynamics%20of%20a%20rigid%20body.%3Cbr%3E%0A%20%20%20%20For%20a%20two-dimensional%20(at%20the%20%24xy%24%20plane)%20movement%2C%20their%20general%20form%20are%20given%20by%3A%3Cbr%3E%0A%20%20%20%20%3Cbr%20%2F%3E%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Csum%20%5Cmathbf%7BF%7D%20%3D%20m%20%5Cmathbf%7B%5Cddot%7Br%7D%7D_%7Bcm%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Csum%20%5Cmathbf%7BT%7D_z%20%3D%20I_%7Bcm%7D%20%5Cmathbf%7B%5Cddot%7B%5Calpha%7D%7D_z%0A%20%20%20%20%24%24%0A%0A%20%20%20%20Where%20the%20movement%20is%20considered%20around%20the%20center%20of%20mass%20(%24cm%24)%20of%20the%20body%2C%20%24%5Cmathbf%7BF%7D%24%20and%20%24%5Cmathbf%7BT%7D%24%20are%2C%20respectively%2C%20the%20forces%20and%20torques%20acting%20on%20the%20body%2C%20%24%5Cmathbf%7B%5Cddot%7Br%7D%7D%24%20and%20%24%5Cmathbf%7B%5Cddot%7B%5Calpha%7D%7D%24%20are%2C%20respectively%2C%20the%20linear%20and%20angular%20accelerations%2C%20and%20%24I%24%20is%20the%20body%20moment%20of%20inertia%20around%20the%20%24z%24%20axis%20passing%20through%20the%20body%20center%20of%20mass.%3Cbr%3E%0A%20%20%20%20It%20can%20be%20convenient%20to%20describe%20the%20rotation%20of%20the%20body%20around%20other%20point%20than%20the%20center%20of%20mass.%20In%20such%20cases%2C%20we%20express%20the%20moment%20of%20inertia%20around%20a%20reference%20point%20%24o%24%20instead%20of%20around%20the%20body%20center%20of%20mass%20and%20we%20will%20have%20an%20additional%20term%20to%20the%20equation%20for%20the%20torque%3A%3Cbr%3E%0A%20%20%20%20%3Cbr%20%2F%3E%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Csum%20%5Cmathbf%7BT%7D_%7Bz%2CO%7D%20%3D%20I_%7Bo%7D%20%5Cmathbf%7B%5Cddot%7B%5Calpha%7D%7D_z%20%2B%20%5Cmathbf%7Br%7D_%7Bcm%2Co%7D%5Ctimes%20m%20%5Cmathbf%7B%5Cddot%7Br%7D%7D_o%0A%20%20%20%20%24%24%0A%0A%20%20%20%20Where%20%24%5Cmathbf%7Br%7D_%7Bcm%2Co%7D%24%20is%20the%20position%20vector%20of%20the%20center%20of%20mass%20in%20relation%20to%20the%20reference%20point%20%24o%24%20and%20%24%5Cmathbf%7B%5Cddot%7Br%7D%7D_o%24%20is%20the%20linear%20acceleration%20of%20this%20reference%20point.%3Cbr%3E%0A%20%20%20%20%3Ca%20href%3D%22free_body_diagram.py%22%3ESee%20this%20notebook%20about%20free-body%20diagram%3C%2Fa%3E.%3Cbr%3E%0A%20%20%20%20%2F%2F%2F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20For%20the%20case%20of%20an%20inverted%20pendulum%20representing%20at%20the%20sagittal%20plane%20the%20body%20at%20the%20standing%20still%20posture%2C%20let's%20solve%20the%20Newton-Euler%20equations%20considering%20the%20rotation%20around%20the%20ankle%20joint%20(because%20it%20will%20simplify%20the%20problem)%20and%20let's%20adopt%20the%20simplifications%20that%20the%20feet%20don't%20move%20(including%20the%20ankle%2C%20so%20%24%5Cmathbf%7B%5Cddot%7Br%7D%7D_o%3D0%24%20in%20the%20Newton-Euler%20equation)%20and%20their%20mass%20are%20neglegible%20in%20relation%20to%20the%20mass%20of%20the%20rest%20of%20the%20body.%0A%0A%20%20%20%20For%20the%20feet%2C%20we%20have%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cbegin%7Barray%7D%7Bl%20l%7D%0A%20%20%20%20-F_%7Bax%7D%20%2B%20GRF_x%20%3D%200%20%5C%5C%0A%20%20%20%20%5C%5C%0A%20%20%20%20-F_%7Bay%7D%20%2B%20GRF_y%20%3D%200%20%5C%5C%0A%20%20%20%20%5C%5C%0A%20%20%20%20-T_a%20%2B%20COP%20%5Ccdot%20GRF_y%20%2B%20h_f%20%5Ccdot%20GRF_x%20%3D%200%0A%20%20%20%20%5Cend%7Barray%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20And%20for%20the%20rest%20of%20the%20body%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cbegin%7Barray%7D%7Bl%20l%7D%0A%20%20%20%20F_%7Bax%7D%20%3D%20m%5Cddot%7Bx%7D_%7Bcm%7D%20%5C%5C%0A%20%20%20%20%5C%5C%0A%20%20%20%20F_%7Bay%7D%20-%20mg%20%3D%20m%5Cddot%7By%7D_%7Bcm%7D%20%5C%5C%0A%20%20%20%20%5C%5C%0A%20%20%20%20T_a%20-%20COG_v%20%5Ccdot%20mg%20%3D%20I_a%20%5Cddot%7B%5Calpha%7D%0A%20%20%20%20%5Cend%7Barray%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20Where%20%24I_a%24%20is%20the%20moment%20of%20inertia%20of%20the%20whole%20body%20around%20the%20ankle%20joint.%3Cbr%3E%0A%20%20%20%20During%20the%20standing%20still%20posture%2C%20the%20%24GRF%24%20horizontal%20component%20is%20typically%20much%20smaller%20than%20the%20%24GRF%24%20vertical%20component%20and%20the%20torque%20of%20the%20former%20can%20be%20neglected.%20In%20addition%2C%20the%20magnitude%20of%20the%20%24GRF%24%20vertical%20component%20is%20approximately%20constant%20and%20equal%20to%20the%20body%20weight.%3Cbr%3E%0A%20%20%20%20Considering%20these%20approximations%2C%20the%20ankle%20joint%20torque%20using%20the%20equation%20for%20the%20feet%20is%20given%20by%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20T_a%20%5Capprox%20COP%20%5Ccdot%20mg%0A%20%20%20%20%24%24%0A%0A%20%20%20%20If%20now%20we%20substitute%20the%20ankle%20joint%20torque%20term%20in%20the%20equation%20for%20the%20torques%20calculated%20for%20the%20rest-of-body%20segment%2C%20we%20have%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20COP%20-%20COG_v%20%5Capprox%20%5Cfrac%7BI_a%7D%7Bmg%7D%20%5Cddot%7B%5Calpha%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20That%20is%2C%20the%20angular%20acceleration%20of%20the%20body%20is%20proportional%20to%20the%20difference%20between%20%24COP%24%20and%20%24COG_v%24%20displacements%20(with%20respect%20to%20the%20ankle%20joint%20position).%0A%0A%20%20%20%20We%20can%20continue%20with%20the%20deduction%20and%20now%20substitute%20the%20angular%20displacement%20by%20a%20term%20proportional%20to%20%24COG_v%24%20if%20we%20use%20the%20following%20trignometric%20relation%20(see%20figure%20above)%3A%20%24sin%20%5Calpha%3DCOG_v%2Fh%24.%20But%2C%20during%20the%20standing%20still%20posture%20%24%5Calpha%24%20is%20typicall%20very%20small%20and%20we%20can%20approximate%20%24sin%5Calpha%20%5Capprox%20%5Calpha%24%20and%20%24%5Calpha%20%5Capprox%20COG_v%2Fh%24.%20However%2C%20bear%20in%20mind%20that%20%24%5Calpha%24%20is%20defined%20as%20counterclockwise%20positive%20while%20%24COG_v%24%20is%20positive%20when%20pointing%20to%20the%20right%20direction.%20This%20means%20that%20in%20fact%20%24%5Calpha%20%5Capprox%20-COG_v%2Fh%24.%20As%20%24h%24%20is%20constant%2C%20the%20second%20derivative%20of%20%24%5Calpha%24%20with%20respect%20to%20time%20is%20simply%20the%20second%20derivative%20of%20%24COG_v%24%20divided%20by%20%24h%24.%0A%0A%20%20%20%20Finally%2C%20the%20last%20equation%20can%20be%20expressed%20in%20the%20following%20form%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20COG_v%20-%20COP%20%5Capprox%20%5Cfrac%7BI_a%7D%7Bmgh%7D%20%5Cddot%7BCOG%7D_v%0A%20%20%20%20%24%24%0A%0A%20%20%20%20Or%20simply%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20COG_v%20-%20COP%20%5Capprox%20k%20%5C%2C%20%5Cddot%7BCOG%7D_v%0A%20%20%20%20%24%24%0A%0A%20%20%20%20Where%20%24k%20%3D%20I_a%2F(mgh)%24.%3Cbr%3E%0A%20%20%20%20If%20the%20human%20body%20is%20represented%20as%20a%20rigid%20bar%2C%20its%20moment%20of%20inertia%20will%20be%20approximately%20equal%20to%20%241.33mh%5E2%24%2C%20so%20%24k%20%5Capprox%201.33h%2Fg%24.%3Cbr%3E%0A%20%20%20%20In%20turn%2C%20from%20the%20Newton-Euler%20equations%2C%20the%20horizontal%20acceleration%20of%20%24COG%24%20is%20equal%20to%20the%20horizontal%20component%20of%20%24GRF%24%20divided%20by%20the%20body%20mass%20and%20the%20equation%20above%20can%20be%20expressed%20as%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20COG_v%20-%20COP%20%5Capprox%20%5Cfrac%7Bk%7D%7Bm%7D%20GRF_x%0A%20%20%20%20%24%24%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20**Exercise%201%20%E2%80%94%20Read%20the%20acceleration%20sign.**%20For%20%24x_%7BCOM%7D-x_%7BCOP%7D%3DK%5Cddot%20x_%7BCOM%7D%24%20with%20%24K%3D0.1%24%20s%C2%B2%2C%20let%20%24x_%7BCOM%7D%3D0.01%24%20m%20and%20%24x_%7BCOP%7D%3D0.03%24%20m.%20Find%20COM%20acceleration.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%2F%2F%2F%20details%20%7C%20Solution%20to%20Exercise%201%0A%0A%20%20%20%20%24%5Cddot%20x_%7BCOM%7D%3D(0.01-0.03)%2F0.1%3D-0.2%24%20m%2Fs%C2%B2.%20Moving%20COP%20ahead%20of%20COM%20produces%20backward%20acceleration%20under%20these%20sign%20conventions.%20COP%20is%20a%20force-resultant%20location%2C%20not%20the%20body%20mass%20center.%0A%20%20%20%20%2F%2F%2F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Implications%20of%20the%20inverted%20pendulum%20model%0A%0A%20%20%20%20These%20two%20last%20equations%20express%20a%20very%20simple%20relation%20between%20body%20segment%20parameters%2C%20%24COG_v%24%2C%20%24COP%24%2C%20body%20acceleration%2C%20and%20horizontal%20force.%3Cbr%3E%0A%20%20%20%20Solely%20based%20on%20these%20equations%20it%20is%20possible%20to%20predict%20some%20interesting%20relations%20among%20the%20variables%20in%20these%20equations%2C%20which%20have%20been%20experimentally%20observed%3A%3Cbr%3E%0A%20%20%20%20-%20%24COG_v-COP%24%20is%20positively%20correlated%20with%20the%20horizontal%20ground%20reaction%20force%20in%20the%20anterior-posterior%20direction%20(Winter%20et%20al.%201998%3B%20Morasso%20and%20Schieppati%201999%3B%20Zatsiorsky%20and%20Duarte%202000)%3B%3Cbr%3E%0A%20%20%20%20-%20%24COG_v%24%20behaves%20as%20a%20low-pass%20filtered%20version%20of%20the%20%24COP%24%20signal%20and%20this%20fact%20has%20been%20used%20in%20a%20method%20to%20derive%20%24COG_v%24%20from%20the%20%24COP%24%20signal%20(Winter%201995%3B%20Caron%20et%20al.%201997%3B%20Morasso%20and%20Schieppati%201999).%20This%20method%20produces%20similar%20results%20as%20other%20methods%20(Lafond%20et%20al.%202004)%3B%3Cbr%3E%0A%20%20%20%20-%20For%20a%20continuously%20regulated%20inverted%20pendulum%20(like%20the%20standing%20still%20posture)%2C%20the%20common%20frequencies%20of%20%24COG_v%24%20and%20%24COP%24%20signals%20are%20in%20phase%20(Morasso%20and%20Schieppati%201999)%3B%3Cbr%3E%0A%20%20%20%20-%20When%20the%20horizontal%20force%20is%20zero%2C%20%24COG_v%24%20and%20%24COP%24%20coincide%20and%20this%20fact%20has%20been%20used%20as%20a%20method%20to%20derive%20%24COG_v%24%20from%20the%20%24COP%24%20displacement%20and%20the%20horizontal%20%24GRF%24%20(King%20and%20Zatsiorsky%201997%3B%20Zatsiorsky%20and%20Duarte%202000).%0A%0A%20%20%20%20Note%20that%20the%20four%20predictions%20made%20above%20are%20based%20entirely%20on%20the%20mechanical%20derivation%20of%20the%20inverted%20pendulum%20model.%20Nothing%20has%20been%20said%20about%20what%20type%20of%20neural%20control%20is%20being%20used%20for%20the%20regulation%20of%20the%20standing%20posture.%20This%20means%20that%20the%20statements%20above%20are%20consequence%20of%20the%20mechanical%20nature%20of%20the%20modeled%20phenomenon.%0A%0A%20%20%20%20Obviously%2C%20the%20most%20straightforward%20prediction%20of%20the%20single%20inverted%20pendulum%20model%20concerning%20the%20**kinematics%20of%20the%20segments**%20of%20the%20human%20body%20would%20be%20that%20we%20should%20observe%20merely%20the%20motion%20at%20the%20ankle%20joint%20and%20nothing%20in%20the%20other%20joints.%20This%20prediction%20has%20not%20been%20observed%20(see%20for%20example%2C%20Pinter%20et%20al.%202008%20and%20Gunther%20et%20al.%202009).%20During%20standing%20still%2C%20we%20seem%20to%20use%20all%20our%20joints%20and%20this%20is%20task%20dependent.%3Cbr%3E%0A%20%20%20%20However%2C%20one%20important%20point%20to%20consider%20is%20that%20even%20if%20the%20inverted%20pendulum%20fails%20as%20a%20suitable%20model%20of%20the%20**kinematics%20of%20the%20segments**%2C%20the%20inverted%20pendulum%20succeds%20as%20a%20model%20of%20the%20**kinematics%20of%20global%20body%20variables**%2C%20such%20as%20%24COG_v%24%20and%20%24COP%24%2C%20and%20their%20relation%20to%20kinetic%20variables%2C%20the%20external%20forces%20acting%20on%20the%20body.%0A%0A%20%20%20%20Certainly%20everyone%20agrees%20that%20the%20inverted%20pendulum%20model%20is%20insufficient%20to%20capture%20all%20the%20essential%20characteristics%20of%20the%20posture%20during%20standing.%20Nevertheless%2C%20the%20greatest%20power%20of%20the%20single%20inverted%20pendulum%20model%20is%20its%20simplicity%2C%20and%20it%20is%20somewhat%20surprising%20to%20note%20how%20much%20this%20simple%20model%20can%20capture%20the%20of%20the%20investigated%20phenomenon.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20**Exercise%202%20%E2%80%94%20Compare%20frequency%20components.**%20For%20a%20sinusoidal%20COP%20at%20angular%20frequency%20%24%5Comega%24%2C%20show%20the%20COM%2FCOP%20amplitude%20ratio%20implied%20by%20the%20linear%20model.%20Evaluate%20it%20at%20%24%5Comega%3D2%24%20rad%2Fs%20and%20%24K%3D0.1%24%20s%C2%B2.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%2F%2F%2F%20details%20%7C%20Solution%20to%20Exercise%202%0A%0A%20%20%20%20For%20a%20sinusoid%2C%20%24%5Cddot%20x%3D-%5Comega%5E2x%24.%20Thus%20%24(1%2BK%5Comega%5E2)x_%7BCOM%7D%3Dx_%7BCOP%7D%24%20and%20the%20ratio%20is%20%241%2F(1%2BK%5Comega%5E2)%3D1%2F1.4%5Capprox0.714%24.%20The%20estimate%20attenuates%20higher-frequency%20COP%20components%20more%20strongly.%0A%20%20%20%20%2F%2F%2F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Estimation%20of%20COGv%20from%20the%20COP%20signal%0A%0A%20%20%20%20Based%20on%20the%20inverted%20pendulum%20model%2C%20it's%20possible%20to%20estimate%20the%20%24COG_v%24%20displacement%20from%20the%20%24COP%24%20displacement%20after%20some%20mathematical%20manipulation%20we%20show%20next.%3Cbr%3E%0A%20%20%20%20Back%20to%20the%20relation%20between%20%24COG_v%24%20and%20%24COP%24%20displacements%2C%20it%20has%20the%20form%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20y(t)%20-%20x(t)%20%3D%20k%5C%2C%5Cddot%7By%7D(t)%0A%20%20%20%20%24%24%0A%0A%20%20%20%20Where%20%24y(t)%24%20stands%20for%20the%20%24COG_v%24%20signal%20and%20%24x(t)%24%20for%20the%20%24COP%24%20signal%2C%20which%20are%20functions%20of%20time%2C%20and%20%24k%20%3D%20I_a%2F(mgh)%24.%3Cbr%3E%0A%20%20%20%20The%20equation%20above%20is%20a%20linear%20ordinary%20differential%20equation%20of%20second%20order.%20This%20equation%20is%20solvable%20in%20the%20time%20domain%2C%20but%20if%20we%20transform%20it%20to%20the%20frequency%20domain%20using%20the%20Fourier%20transform%2C%20we%20will%20find%20a%20simpler%20relation%20between%20%24COG_v%24%20and%20%24COP%24.%0A%0A%20%20%20%20%3Cbr%3E%0A%20%20%20%20%2F%2F%2F%20admonition%20%7C%20The%20Fourier%20transform%0A%0A%20%20%20%20The%20%3Ca%20href%3D%22https%3A%2F%2Fen.wikipedia.org%2Fwiki%2FFourier_transform%22%3EFourier%20transform%3C%2Fa%3E%20is%20a%20mathematical%20operation%20to%20transform%20a%20signal%20which%20is%20function%20of%20time%2C%20%24g(t)%24%2C%20into%20a%20signal%20which%20is%20function%20of%20frequency%2C%20%24G(f)%24%2C%20and%20it%20is%20defined%20by%3A%3Cbr%3E%0A%20%20%20%20%3Cbr%20%2F%3E%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathcal%7BF%7D%5Bg(t)%5D%20%3D%20G(f)%20%3D%20%5Cint_%7B-%5Cinfty%7D%5E%7B%5Cinfty%7D%20g(t)%20e%5E%7B-i2%5Cpi%20ft%7D%20dt%0A%20%20%20%20%24%24%0A%0A%20%20%20%20Its%20inverse%20operation%20is%3A%3Cbr%3E%0A%20%20%20%20%3Cbr%20%2F%3E%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathcal%7BF%7D%5E%7B-1%7D%5BG(f)%5D%20%3D%20g(t)%20%3D%20%5Cint_%7B-%5Cinfty%7D%5E%7B%5Cinfty%7D%20G(f)%20e%5E%7Bi2%5Cpi%20ft%7D%20df%0A%20%20%20%20%24%24%0A%0A%20%20%20%20The%20function%20%24G(f)%24%20is%20the%20representation%20in%20the%20frequency%20domain%20of%20the%20time-domain%20signal%2C%20%24g(t)%24%2C%20and%20vice-versa.%20The%20functions%20%24g(t)%24%20and%20%24G(f)%24%20are%20referred%20to%20as%20a%20Fourier%20integral%20pair%2C%20or%20Fourier%20transform%20pair%2C%20or%20simply%20the%20Fourier%20pair.%3Cbr%3E%0A%20%20%20%20%3Ca%20href%3D%22http%3A%2F%2Fwww.thefouriertransform.com%2Ftransform%2Ffourier.php%22%3ESee%20here%20for%20an%20introduction%20to%20Fourier%20transform%3C%2Fa%3E%20and%20%3Ca%20href%3D%22http%3A%2F%2Fwww.thefouriertransform.com%2Fapplications%2Fdifferentialequations.php%22%3Esee%20here%20for%20the%20use%20of%20Fourier%20transform%20to%20solve%20differential%20equations%3C%2Fa%3E.%0A%20%20%20%20%2F%2F%2F%0A%0A%20%20%20%20%3Cbr%3E%0A%20%20%20%20Let's%20apply%20the%20Fourier%20transform%20to%20the%20differential%20equation%20with%20%24COG_v%24%20and%20%24COP%24%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20Y(j%5Comega)%20-%20X(j%5Comega)%20%3D%20-k%5C%2C%5Comega%5E2Y(j%5Comega)%0A%20%20%20%20%24%24%0A%0A%20%20%20%20Where%20we%20defined%20%24y(t)%20%5CLeftrightarrow%20Y(j%5Comega)%24%20and%20%24x(t)%20%5CLeftrightarrow%20X(j%5Comega)%24%20as%20the%20Fourier%20pairs%2C%20%24j%24%20is%20the%20imaginary%20unit%2C%20and%20%24%5Comega%24%20is%20the%20angular%20frequency%2C%20%242%5Cpi%20f%24.%3Cbr%3E%0A%20%20%20%20The%20reason%20why%20we%20use%20the%20Fourier%20transform%20is%20because%20we%20started%20with%20a%20second%20order%20differential%20equation%20and%20ended%20with%20the%20simple%20algebraic%20equation%20above.%20Rearranging%20the%20equation%20above%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cfrac%7BY(j%5Comega)%7D%7BX(j%5Comega)%7D%20%3D%20%5Cfrac%7B%5Comega_0%5E2%7D%7B%5Comega_0%5E2%20%2B%20%5Comega%5E2%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20Where%20%24%5Comega_0%20%3D%201%2F%5Csqrt%7Bk%7D%24.%3Cbr%3E%0A%20%20%20%20If%20we%20imagine%20a%20system%20where%20the%20%24COP%24%20is%20the%20input%20and%20the%20%24COG_v%24%20the%20output%2C%20the%20right%20side%20of%20the%20equation%20above%20is%20known%20as%20the%20%3Ca%20href%3D%22https%3A%2F%2Fen.wikipedia.org%2Fwiki%2FTransfer_function%22%3Etransfer%20function%3C%2Fa%3E%20of%20such%20system%2C%20the%20ratio%20between%20the%20output%20and%20the%20input.%3Cbr%3E%0A%20%20%20%20Analysing%20the%20transfer%20function%20given%20in%20the%20equation%20above%2C%20we%20see%20that%20it%20is%20of%20the%20type%20of%20a%20low-pass%20filter%20(and%20%24%5Comega_0%24%20is%20the%20cutoff%20frequency)%3B%20because%20of%20that%2C%20we%20can%20say%20that%20the%20%24COGv%24%20signal%20is%20a%20low-pass%20filtered%20version%20of%20the%20%24COP%24%20signal.%0A%0A%20%20%20%20We%20can%20implement%20such%20low-pass%20filter%20in%20order%20to%20determine%20the%20%24COG_v%24%20using%20the%20%24COP%24%20signal.%20For%20that%2C%20we%20simply%20have%20to%20estimate%20the%20Fourier%20transform%20of%20the%20%24COP%24%20signal%2C%20multiply%20by%20the%20transfer%20function%20(the%20right%20side%20of%20the%20equation%20above)%2C%20and%20calculate%20the%20inverse%20Fourier%20transform%20of%20this%20result.%20The%20Python%20function%20%60cogve.py%60%20(code%20at%20the%20end%20of%20this%20text)%20estimates%20the%20%24COG_v%24%20using%20the%20%24COP%24%20data%20based%20on%20this%20algorithm.%20Let's%20test%20this%20function%2C%20first%20we%20have%20to%20import%20the%20necessary%20Python%20libraries%20and%20configure%20the%20emvironment%3A%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_()%3A%0A%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20from%20notebook_resources%20import%20repo_file%0A%20%20%20%20except%20ModuleNotFoundError%20as%20_error%3A%0A%20%20%20%20%20%20%20%20if%20_error.name%20!%3D%20%22notebook_resources%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20raise%0A%20%20%20%20%20%20%20%20raise%20ModuleNotFoundError(%0A%20%20%20%20%20%20%20%20%20%20%20%20%22This%20notebook%20needs%20notebook_resources.py.%20Download%20it%20from%3A%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%22https%3A%2F%2Fraw.githubusercontent.com%2FBMClab%2FBMC%2Fmaster%2Ffunctions%2Fnotebook_resources.py%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Save%20it%20in%20the%20same%20folder%20as%20this%20notebook%2C%20then%20rerun%20this%20cell.%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20name%3D%22notebook_resources%22%2C%0A%20%20%20%20%20%20%20%20)%20from%20None%0A%20%20%20%20return%20(repo_file%2C)%0A%0A%0A%40app.cell%0Adef%20_(repo_file)%3A%0A%20%20%20%20%23%20Import%20the%20necessary%20libraries%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20%23%20'%25matplotlib%20inline'%20command%20supported%20automatically%20in%20marimo%0A%20%20%20%20import%20sys%0A%20%20%20%20sys.path.insert(0%2C%20str(repo_file('functions%2Fcogve.py').parent))%20%20%23%20add%20to%20pythonpath%0A%20%20%20%20from%20cogve%20import%20cogve%0A%0A%20%20%20%20return%20cogve%2C%20plt%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20Let's%20use%20stabilographic%20data%20from%20the%20public%20data%20set%20of%20human%20balance%20evaluations%20of%20Santos%20and%20Duarte%20(2016)%3A%20the%20center%20of%20pressure%20of%20a%20person%20standing%20still%20for%2060%20s%20with%20eyes%20open%20on%20a%20rigid%20surface%20(trial%20BDS00013)%2C%20measured%20at%20100%20Hz%20(the%20full%20data%20set%20is%20available%20at%20%5BGitHub%5D(https%3A%2F%2Fgithub.com%2Fdemotu%2Fdatasets%2Ftree%2Fmaster%2FBDS))%3A%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(repo_file)%3A%0A%20%20%20%20import%20pandas%20as%20pd%20%20%23%20use%20Pandas%20to%20read%20data%20from%20a%20website%0A%20%20%20%20fileUrl%20%3D%20repo_file('BDS%2Fdata%2FBDS00013.txt'%2C%20base_url%3D%22https%3A%2F%2Fraw.githubusercontent.com%2Fdemotu%2Fdatasets%2Fmaster%2F%22)%0A%20%20%20%20data%20%3D%20pd.read_csv(fileUrl%2C%20sep%3D'%5Ct')%20%20%23%20Pandas%20dataframe%0A%20%20%20%20COP%20%3D%20data%5B%5B'COPx%5Bcm%5D'%2C%20'COPy%5Bcm%5D'%5D%5D.to_numpy()%20%20%23%20ap%20and%20ml%20directions%2C%20in%20cm%0A%20%20%20%20freq%20%3D%20100%0A%20%20%20%20print('COP%20shape%3A%20'%2C%20COP.shape)%0A%20%20%20%20return%20(COP%2C)%0A%0A%0A%40app.cell%0Adef%20_(COP%2C%20cogve%2C%20plt)%3A%0A%20%20%20%20fig%2C%20ax%20%3D%20plt.subplots(1%2C%201%2C%20figsize%3D(8%2C%205))%0A%20%20%20%20cogv%20%3D%20cogve(COP%5B%3A%2C%200%5D%2C%20freq%3D100%2C%20mass%3D70%2C%20height%3D175%2C%20ax%3Dax%2C%20show%3DTrue)%20%20%23%20guess%20mass%2C%20height%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20**Exercise%203%20%E2%80%94%20Identify%20model%20limits.**%20Would%20the%20same%20single-link%20relation%20describe%20a%20deep%20squat%20with%20substantial%20hip%20and%20knee%20motion%3F%20Explain%20before%20trying%20the%20filter.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%2F%2F%2F%20details%20%7C%20Solution%20to%20Exercise%203%0A%0A%20%20%20%20No.%20The%20model%20assumes%20approximately%20fixed%20body%20geometry%2C%20a%20fixed%20ankle%20pivot%20and%20small%20sway.%20A%20deep%20squat%20changes%20COM%20height%2C%20segment%20orientations%20and%20angular%20momentum.%20A%20smooth%20filtered%20trace%20is%20not%20evidence%20that%20those%20assumptions%20hold.%0A%20%20%20%20%2F%2F%2F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Function%20cogve.py%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20repo_file)%3A%0A%20%20%20%20mo.ui.code_editor(repo_file(%22functions%2Fcogve.py%22).read_text()%2C%20language%3D%22python%22%2C%20disabled%3DTrue)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%0A%20%20%20%20return%20(mo%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Checkpoint%20questions%0A%0A%20%20%20%201.%20Is%20center%20of%20pressure%20defined%20in%20flight%3F%0A%20%20%20%202.%20Why%20does%20the%20estimated%20COM%20trace%20look%20smoother%3F%0A%20%20%20%203.%20What%20does%20the%20anthropometric%20height%20parameter%20influence%3F%0A%20%20%20%204.%20Can%20COP%20alone%20validate%20the%20inferred%20COM%20trajectory%3F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%2F%2F%2F%20details%20%7C%20Checkpoint%20answers%0A%0A%20%20%20%201.%20No%20meaningful%20contact%20COP%20exists%20without%20a%20supporting%20normal%20force.%0A%20%20%20%202.%20The%20model%20maps%20COP%20to%20COM%20through%20a%20low-pass%20frequency%20relation.%0A%20%20%20%203.%20The%20pendulum%20geometry%2Finertia%20estimate%20and%20therefore%20the%20filter%20frequency%20scale.%0A%20%20%20%204.%20No%3B%20independent%20kinematics%20or%20another%20COM%20estimate%20is%20needed.%0A%20%20%20%20%2F%2F%2F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Problems%0A%0A%20%20%20%201.%20At%20static%20equilibrium%20in%20this%20model%2C%20how%20are%20COP%20and%20projected%20COM%20related%3F%0A%0A%20%20%20%202.%20For%20%24K%3D0.2%24%20s%C2%B2%20and%20COP%E2%80%93COM%20difference%20%2B0.01%20m%2C%20find%20COM%20acceleration.%0A%0A%20%20%20%203.%20Give%20the%20amplitude%20ratio%20for%20%24K%3D0.1%24%20s%C2%B2%20at%201%20Hz.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%2F%2F%2F%20details%20%7C%20Final%20answers%20to%20Problems%0A%0A%20%20%20%201.%20They%20coincide%20horizontally.%0A%0A%20%20%20%202.%20%E2%88%920.05%20m%2Fs%C2%B2.%0A%0A%20%20%20%203.%20About%200.2021%20(dimensionless).%0A%20%20%20%20%2F%2F%2F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Go%20deeper%0A%0A%20%20%20%20Use%20%5BKineticsFundamentalConcepts%5D(https%3A%2F%2Fmolab.marimo.io%2Fgithub%2FBMClab%2FBMC%2Fblob%2Fmaster%2Fnotebooks%2Fkinetics_fundamental_concepts.py)%20to%20connect%20this%20guide%20with%20the%20surrounding%20mechanics%20and%20modeling%20material.%20Return%20to%20the%20assumptions%20behind%20each%20equation%20when%20applying%20it%20to%20a%20new%20experiment.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%23%20Video%20lectures%20on%20the%20Internet%0A%0A%20%20%20%20%5BMIT%20OpenCourseWare%5D(https%3A%2F%2Focw.mit.edu%2Fcourses%2F2-003sc-engineering-dynamics-fall-2011%2Fpages%2Ffree-body-diagrams-and-fictitious-forces%2F)%20%E2%80%94%20Free-body%20diagrams%3A%20practice%20selecting%20a%20body%20and%20frame%20before%20writing%20force%20and%20moment%20equations.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20References%0A%0A%20%20%20%20-%20Caron%20O%2C%20Faure%20B%2C%20et%20al.%20(1997)%20%5BEstimating%20the%20centre%20of%20gravity%20of%20the%20body%20on%20the%20basis%20of%20the%20centre%20of%20pressure%20in%20standing%20posture%5D(https%3A%2F%2Fpubmed.ncbi.nlm.nih.gov%2F9456386).%20J.%20Biomech.%2030%2C%201169-1171.%0A%0A%20%20%20%20-%20Lafond%20D%2C%20Duarte%20M%2C%20et%20al.%20(2004)%20%5BComparison%20of%20three%20methods%20to%20estimate%20the%20center%20of%20mass%20during%20balance%20assessment%5D(http%3A%2F%2Febm.ufabc.edu.br%2Fpublications%2Fmd%2FJB03.pdf).%20J.%20Biomech.%2037%2C%201421-1426.%0A%0A%20%20%20%20-%20King%20D%2C%20Zatsiorsky%20VM%20(1997)%20%5BExtracting%20gravity%20line%20displacement%20from%20stabilographic%20recordings%5D(http%3A%2F%2Fwww.sciencedirect.com%2Fscience%2Farticle%2Fpii%2FS0966636296011010).%20Gait%20%26%20Posture%206%2C%2027-38.%0A%0A%20%20%20%20-%20Morasso%20PG%2C%20Spada%20G%2C%20et%20al.%20(1999)%20%5BComputing%20the%20COM%20from%20the%20COP%20in%20postural%20sway%20movements%5D(http%3A%2F%2Fwww.sciencedirect.com%2Fscience%2Farticle%2Fpii%2FS0167945799000391).%20Human%20Movement%20Science%2018%2C%20759-767.%0A%0A%20%20%20%20-%20Santos%20DA%2C%20Duarte%20M%20(2016)%20%5BA%20public%20data%20set%20of%20human%20balance%20evaluations%5D(https%3A%2F%2Fdoi.org%2F10.7717%2Fpeerj.2648).%20PeerJ%2C%204%2C%20e2648.%0A%0A%20%20%20%20-%20Winter%20DA%20(1995)%20%5BA.B.C.%20(Anatomy%2C%20Biomechanics%20and%20Control)%20of%20Balance%20during%20Standing%20and%20Walking%5D(https%3A%2F%2Fbooks.google.com.br%2Fbooks%3Fid%3D0lSqQgAACAAJ%26).%20Waterloo%2C%20Waterloo%20Biomechanics.%0A%0A%20%20%20%20-%20Winter%20DA%2C%20Patla%20AE%2C%20et%20al.%20(1998)%20%5BStiffness%20control%20of%20balance%20in%20quiet%20standing%5D(https%3A%2F%2Fpubmed.ncbi.nlm.nih.gov%2F9744933).%20J.%20Neurophysiol.%2080%2C%201211-1221.%0A%0A%20%20%20%20-%20Zatsiorsky%20VM%2C%20Duarte%20M%20(2000)%20%5BRambling%20and%20trembling%20in%20quiet%20standing%5D(http%3A%2F%2Febm.ufabc.edu.br%2Fpublications%2Fmd%2FMC00.pdf).%20Motor%20Control%204%2C%20185-200.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A