Fink2014-model

MATLAB and python code to simulate the biomechanical model described in Fink et al. (2014)

Author

Akira Nagamori

File

1. biomechanical_model.m (MATLAB code)
2. biomechanical_model.ipynb (Python code in Jupyter notebook)

Model Descriptions

The model is simulated using its state-space representation. Inputs to flexor and extensor, I_e and I_f, are converted into extensor and flexor torques, T_e and T_f. These toruqes drive joint angle kinematics described by theta and theta_dot. Note that the model variables are dimensionless except for time. All parameters are the same as those in Fink et al., (2014). Feedback gain (h) and delay (t_delay) can be adjusted to manipulate effects of sensory feedback on joint kinematics.

1. State vector, x = [theta;theta_dot;T_e;T_e_dot;T_f;T_f_dot]
2. Input vector, u = [I_e;I_f]
3. State matrix, A
4. Delay matrix, A1
5. Input matrix, B

Reference

Fink AJ, Croce KR, Huang ZJ, Abbott LF, Jessell TM, Azim E. Presynaptic inhibition of spinal sensory feedback ensures smooth movement. Nature. 2014 May 1;509(7498):43-8. doi: 10.1038/nature13276. PMID: 24784215; PMCID: PMC4292914.