A simulation method for the firing sequences of motor units (Jiang et al 2006)

 Download zip file 
Help downloading and running models
Accession:83320
" ... a novel model based on the Hodgkin–Huxley (HH) system is proposed, which has the ability to simulate the complex neurodynamics of the firing sequences of motor neurons. The model is presented at the cellular level and network level, and some simulation results from a simple 3-neuron network are presented to demonstrate its applications." See paper for more and details.
Reference:
1 . Jiang N, Englehart KB, Parker PA (2007) A simulation method for the firing sequences of motor units. J Electromyogr Kinesiol 17:527-34 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Spinal motoneuron;
Cell Type(s):
Channel(s): I Na,t; I K;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: MATLAB;
Model Concept(s): Activity Patterns; Simplified Models;
Implementer(s): Jiang, Ning [ning.jiang at unb.ca];
Search NeuronDB for information about:  I Na,t; I K;
function inp=gene_inp(iapp,stdin,duration,step_size,cdfun)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%This routine generate the inp current to a MN in the network.
%inputs:
%iapp: individual input mean (bias current, or base current)
%stdin: individual input standard deviation
%step_size: step_size of the ODE solver
%cdfun: common drive func
%
%Written by Ning Jiang, Institute of Biomedical Engineering, Univesity of New
%Brunswick, NB, Canada, E3B 5A3
%Email: ning.jiang@unb.ca
%
%Date: Nov 19, 2006
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

filter_len=16;
%put a limit on the possible value of inp current
SCALING=1;
%maximum inp current value
inp_inf=30;


b=hann(filter_len);
b=b/sum(b);
rstd=stdin/sqrt(sum(b.^2));
inp=randn(1,duration*1000/step_size+filter_len+1)*rstd+iapp;
inp=filter(b,1,inp);
inp=inp(filter_len+1:length(inp));
inp=inp+cdfun;
%limiting the input current
if SCALING
    inp=30*tanh(inp/inp_inf);
end
    

Loading data, please wait...