Data & Program files for
Anticipatory Postural Control Emerges from a Predictive and Optimized Strategy for Movement Preparation
Anticipatory Postural Control Emerges from a Predictive and Optimized Strategy for Movement Preparation
Tetsuro Funato, Miho Ogawa, Akira Konosu, Dai Yanagihara
Motion & EMG data files (187 MB)
Program files (5 KB)
Note: Create "data" folder in the same folder with program files and place all data in the "measuredDataFiles.zip" to the "data" folder.
MotionData_subject[subject number].xlsx
comX, angAnkle, angHip, angForcePlate, posAnkleX: Time series of motion data. Rows correspond to time, and columns correspond to trials (60 trials in total). Each column shows the time series of motion measured by a motion capture system for each trial. The time series range in the data is 5 seconds before and after the start of the floor tilting. The measurement frequency is 300 Hz. "comX": horizontal COM displacement (mm), "angAnkle": ankle angle (rad), "angHip": hip angle (rad), "angForcePlate": force plate angle (rad), and "posAnkleX": horizontal ankle displacement (mm).
tiltEndIndx, Soundstart, Soundend: Time of end of floor tilting, start of cue, and stop of cue. Each column corresponds to trial (60 trials in total). Each value is represented as the index number of the measured data. "tiltEndIndx" is measured with 300 Hz, "Soundstart" and "Soundend" are measured with 1000 Hz.The actual time (seconds) is calculated by dividing the data index by 300 or 1000.
EMGData_subject[subject number].xlsx
RTA, LTA, RGAS, LGAS: Time series of Electromyography (EMG) data. Rows correspond to time, and columns correspond to trials (60 trials in total). Each column shows the time series of EMG for each trial. The time series range in the data is 5 seconds before and after the start of the floor tilting. The measurement frequency is 300 Hz. Each value represents %MVC, which is the measured EMG divided by maximum muscle activity. "RTA": right tibialis anterior (TA), "LTA": left TA, "RGAS": right gastrocnemius (GC), and "LGAS": left GC.
analysisMotion.m
The program loads measured motion data and analyzes the time series of each subject's movements and differences in movements between cue and no cue tasks.
Input data: MotionData_subject[subject number].xlsx
Output data: tResult{COM/Hip/Ank}.csv: Results of t-tests to determine significant differences in COM/hip/ankle changes in the CS-FS interval between cue and no cue tasks.
Output figures: {COM/Hip/Ankle}_subject{subject_number}.emf: Time series of COM/Hip/Ankle.
COMAnkleHipdiff.emf: The amount of change in COM/hip/ankle values for all subjects between CS and FS.
analysisEMG.m
The program loads measured EMG data and analyzes the time series of each subject's EMGs and differences in EMGs between cue and no cue tasks.
Input data: EMGData_subject[subject number].xlsx
Output data: tResult{TA/GC}.csv: Results of t-tests to determine significant differences in TA/GC activities in the CS-FS interval between cue and no cue tasks.
Output figures: {TA/GC}_subject{subject_number}.emf: Time series of TA/GCe.
cTAGC_ALL.emf: Coefficient for linear regression of TA/GC activities in the CS-FS interval for all subjects.
Note: Create "simResults" folder in the same folder with program files to store calculation results.
stand_ms.m
Main program.
Output files: simResults/simResultData_{simulation number}.csv: Simulation results (Time series of COM, Joint angles, Muscle activation, Floor angles, etc.)
simResults/simInfo_{simulation number}.csv: Properties used for simulation (Seed of the random number, COM Reference, Weight for the evaluation function)
diffEqLink2_ms.m
Differential equation of the system
odeEulerModel.m
Internal model used for MPC
standOutFunc_ms.m
Output function used for MPC
showResultFromData.m
Load simulation results in "simResults" folder and display results.