Displaying 2 results from an estimated 2 matches for "izhikevich".
2013 Feb 10
4
A Hodgkin Huxley Model
Hi All
It has been suggested to me that the folks in this list might be able to
help me out of my misery.
As part of my learning curve to building a rather detailed model of a small
neurone complex I am implementing some existing models in R.
For the moment I want to implement the Izhikevich model as described in his
2003 paper. The equations are as follows:
v' = 0.04v^2 + 5v + 140 - u - I
u' = a(bv - u)
if v=30 then v<-c else u<-u+d
I don't know how to implement the if condition in my R code. I'm using
deSolve. Here is my code so far.
library(deSolve);
Izhikev...
2013 Feb 21
1
using and event in deSolve
...n the ball hits the ground (height = 0) Then
velocity (y2) is reversed
and reduced by 10 percent. The root function, y[1] = 0, triggers the event:
> root <- function(t, y, parms) y[1]
Firstly I couldn't see where y[1] became 0, but I implemented Izchikevich
as follows:
library(deSolve);
Izhikevich <- function(time, init, parms) {
with(as.list(c(init, parms)),{
dv <- (0.04*v^2)+(5*v)+140-u+I;
du <- a*(b*v-u);
#if (v>=30) v<-c else v<-u+d;
list( c(dv, du))
})}
parms=c( a=0.02, b=0.2, c=-65, d=2, I=10);
times=seq(from=1, to=1000, by=0.1);
init=c(v=-65, u=0....