ikp.m

来自「基于Matlab的心肌细胞动作电位模型程序」· M 代码 · 共 20 行

function [I,gKp] = IKp(V)
% IKp    Plateau potassium current
%   I = IKp(V)
%    I   potassium current
%  gKp   conductance
%    V   membrane potential (mV)
% Reference: QN:Luo94a

global R T F K_o K_i;
GKp_ = 0.0183;
EKp = (R*T/F)*log(K_o/K_i);

numTimeStep = length(V);

for i=1:numTimeStep
  Kp = 1/(1+exp((7.488-V(i))/5.98));
  gKp(i) = GKp_*Kp;
  I(i) = gKp(i)*(V(i)-EKp);
end