Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision		 Previous revision
technical:poisson_spike_train [2021/09/21 09:05] – 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.log(1-np.random.uniform())+        s += rng.exponential()
         while s>r[i]*dt:         while s>r[i]*dt:
             s -= r[i]*dt             s -= r[i]*dt
Line 22: Line 23:
             continue             continue
         break         break
-    return spks+    return np.array(spks)
 </code> </code>