hh模型仿真matlab程序,一个关于HH模型的基础问题

本帖最后由 summerworm 于 2014-4-19 20:43 编辑

用龙格-库塔法模拟神经元动作电位的发放写的代码……但是算不出来，值是无穷大，不知道哪里出了问题，求大神帮忙看看~

谢谢各位啦~如下为程序代码：

function y=HH_RungeKutta(dt,tfinal,Istim)

% input constant

vna=50; vk=-77; vl=-54.4 ;   %unit=mV

m=0.05; h=0.6; n=0.3; v=-70;

t=0; tin=50; i=0;

while t<=tfinal

t=t+dt;

i=i+1;

if t

signal=0;

else

signal=Istim;

end

[am,bm,ah,bh,an,bn,gna,gk,gl]=alphabeta(v,m,n,h);

k1n=an*(1-n)+bn*n;

k2n=an*(1-n-0.5*dt*k1n)+bn*(n+0.5*dt*k1n);

k3n=an*(1-n-0.5*dt*k2n)+bn*(n+0.5*dt*k2n);

k4n=an*(1-n-dt*k3n)+bn*(n+dt*k3n);

n=n+dt*(k1n+2*k2n+2*k3n+k4n)/6;

k1m=am*(1-m)+bm*m;

k2m=am*