Data and code for: Feedback control stabilizing the center of mass can be identified in unperturbed, upright standing
Authors/Creators
Description
Here we provide the code and data associated with our research article titled “Feedback control stabilizing the center of mass can be identified in unperturbed, upright standing.”
In this study, we propose a model relating corrective ground reaction forces to preceding center of mass position and velocity. We first checked the model’s validity by simulating an inverted pendulum model with known intrinsic and feedback parameters, to see whether and to what extend we could effectively identify the feedback parameters. Next, we tested the model in 15 young adult volunteers (21± 4 yrs, 63 ± 9 kg, 1.70 ± 0.10 m, 13 females, 2 males) in normal standing.
By running the provided code, you will be able to reproduce the results presented in the paper. Please refer to the instructions and table below for detailed information about file contents and usage.
Model simulation results
Run the following scripts in order:
1. P1_MS01_RunSimulation
2. P1_MS02_Analysis
3. P1_MS03_plotFigure
4. P1_MS04_TypicalSimFigure
To reproduce the results shown in the manuscript, run all scripts except P1_MS01_RunSimulation, since it includes random perturbations that will generate slightly different outcomes each time. If you want to generate new simulation data, make sure to run P1_MS01_RunSimulation first.
Human experimental results
Run the following scripts in order:
1. P1_S01_CalcModelInput
2. P1_S02_CalcFeedback_ParaFig
3. P1_S03_CalcSway
To reproduce the manuscript results, run these scripts in order.
|
Folders |
Content |
||
|
Model_Simulation |
code |
P1_MS01_RunSimulation |
Scripts to get all combinations of model inputs: ‘all_combinations.mat’, and simulation results: ‘sim_results_XX.mat’. |
|
P1_MS02_Analysis |
Script to get the linear regression results shown in Table 2, and the unrealistic sway simulations for each parameter value shown in Table 3. |
||
|
P1_MS03_plotFigure |
Script to get the figures illustrating the relationship between estimated and input parameters, as shown in Figure 4. |
||
|
P1_MS04_TypicalSimFigure |
Script to get the model simulation results with default input parameters, as well as estimated human control parameters, as shown in Figure 3 and 9. |
||
|
MS_function |
Functions which are used by the scripts |
||
|
all_combinations.mat |
All possible input combinations. |
||
|
sim_results_01.mat … sim_results_10.mat |
Simulation results, the number ‘01’, ‘02’, …, ‘10’ represents the number of the simulation. |
||
|
sim_default.mat |
Simulation results under default input parameters |
||
|
sim_human.mat |
Simulation results under the estimated human control parameters |
||
|
Human_Experiment |
code |
P1_S01_CalcModelInput |
Script to get the feedback model input, you will get the COM kinematics and Force data, save in the ‘result’ folder. |
|
P1_S02_CalcFeedback_ParaFig |
Script to get the feedback parameters and related figures; for each subject, you will get the feedback results of each trial and average value of all trials under normal standing conditions. |
||
|
P1_S03_CalcSway |
Script to get the sway results: mean, SD, RMS of COM,VCOM,COP, VCOP, and body sway. |
||
|
=VU 3D model new= |
Necessary functions and folders to run the main scripts. |
||
|
P1_functions |
|||
|
ndigital |
|||
|
AllSub_Info.xlsx |
Basic subjects’ information: age, height, mass, gender, null force number and force plates coordination records. |
||
|
Remove_Standing.xlsx |
Remove noise setup, will be used in ‘P1_S01_CalcModelInput’. |
||
|
StandingInfo_AllSub.xlsx |
Standing trial numbers. |
||
|
data |
Sub01, Sub02, …, Sub18 |
NormalStanding.xlsx: trial information XX.ndf: motion capture data XX.afp: force plate data |
|
Files
Feedback control stabilizing the center of mass can be identified in unperturbed, upright standing.zip
Files
(3.2 GB)
| Name | Size | Download all |
|---|---|---|
|
md5:6ce3b74192960903b506c3a27444a4cc
|
3.2 GB | Preview Download |
Additional details
Related works
- Documents
- Dataset: 10.1101/2024.12.20.629656 (DOI)