Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
technical:poisson_spike_train [2022/10/10 05:54] – chunchungtechnical:poisson_spike_train [2022/10/10 06:03] (current) – chunchung
Line 1: Line 1:
 ======Inhomogeneous Poisson spike train generation====== ======Inhomogeneous Poisson spike train generation======
 <code python> <code python>
-def gen_spikes(r,dt):+def gen_spikes(r,dt,rng):
     '''Generate spike train from Poisson rate     '''Generate spike train from Poisson rate
    
     Parameters     Parameters
     ----------     ----------
-    r:  Array of spike rates +    r:   Array of spike rates 
-    dt: Time step+    dt:  Time step 
 +    rng: Random number generator
     '''     '''
     i = 0     i = 0
Line 13: Line 14:
     spks = []     spks = []
     while True:     while True:
-        s += np.random.exponential()+        s += rng.exponential()
         while s>r[i]*dt:         while s>r[i]*dt:
             s -= r[i]*dt             s -= r[i]*dt