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 [2020/10/13 19:36] – 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>
-# generate spike train +def gen_spikes(r,dt,rng): 
-def gen_spikes(r,dt):+    '''Generate spike train from Poisson rate 
 +  
 +    Parameters 
 +    ---------- 
 +    r:   Array of spike rates 
 +    dt:  Time step 
 +    rng: Random number generator 
 +    '''
     i = 0     i = 0
     s = 0     s = 0
     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 16: Line 23:
             continue             continue
         break         break
-    return spks+    return np.array(spks)
 </code> </code>