Index exceeds array elements timei
显示 更早的评论
function membPoten = simulateAdEx(params)
%membPoten = simulateAdEx(params)
% Simulate a neuron according to the adaptive exponential model.
%
% INPUT: params - a structure containing simulation parameters,
% with the following required fields:
%
% params.C : capacitance
% params.gl : leak conductance
% params.El : leak voltage
% params.vt : resting membrane potential
% params.delT : spike width
% params.a : resonantor/integrator constant
% params.tauW : adaptation decay time
% params.b : adaptation jump (for updating w)
% params.vreset : voltage reset value
% params.simDur : simulation time in ms
% params.dt : integration time step in ms
% params.input : input current in pA
% params.EC50 : the drug concentration that produces a 50% maximal response
% params.Max : maximal inhibition
% params.H : Hill coefficient
% params.x : Drug concentration
%
% OUTPUT: membranePotential - A vector of the neuron's membrane
% potential in mV.
%
%
% Equations based on:
% https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2798047/
% Read more about this model:
% http://www.scholarpedia.org/article/Adaptive_exponential_integrate-and-fire_model
%
%
% vector to store membrane potential
membPoten = [params.vt; zeros(params.simDur/params.dt,1)];
v = params.vt;
w = 0;
vpeak = 0; % spike threshold
% simulation in loop over time
for timei = 1:params.simDur/params.dt-1
% set input for this time point
Iin = params.input(timei);
% integration
v = v + params.dt*(-params.gl*(v-params.El) + params.gl*params.delT*exp((v-params.vt)/params.delT) + Iin - w)/params.C;
% note: update w by pre-updated v
w = w + params.dt*(params.a*(membPoten(timei) - params.El) - w)/params.tauW;
% update and store variables
if v>=vpeak
v = params.vreset;
w = w + params.b;
else
% no AP, just store v
membPoten(timei+1) = v;
end
end
end % end function
Hello, I am getting an error in line 47 for index exceeds the number of array elements (1). Please help with this problem.
1 个评论
Walter Roberson
2022-3-9
As outside observers, we have no reason to believe that params.input is any particular size or datatype or range. You did not attach the inputs for the function.
My guess is that params.input is a scalar but the code expects there to be one input for each time step.
回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Neural Simulation 的更多信息
2022-3-9
2022-3-9
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
