fdvdt103.m

利用Matlab模拟细胞HH 模型中的所有参数变化和动态Action Potential 生成过程

function dvdt = fdvdt103(nv, vv, iext);
% FDVDT103 - computes value of dV/dt for reduced, 2-variable HH model:
% dV/dt = {iext - Gna*minf(vv)^3*nvtohv(nv)*(V-Ena) - Gk*ninf^4*(V-Ek) - Gleak*(V-Eleak)}/C
% NOTE sethhparms must be called first to set up HH parameters.

% Parameters (mS/cm^2, mV, uFd/cm^2)
% hps = [gnamx=120; Ena=55; gkmx=36; Ek=-72; gleak=0.3; El=-49.4; C=1]
hps = gethhparms;

dvdt = (iext - hps(1)*minfhh(vv).^3.*nvtohv(nv).*(vv-hps(2)) - ...
          hps(3)*nv.^4.*(vv-hps(4)) - hps(5)*(vv-hps(6)))/hps(7);
return