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
Last revisionBoth sides next revision
technical:poisson_spike_train [2020/10/13 19:36] – chunchungtechnical:poisson_spike_train [2022/10/10 05:54] – 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): def gen_spikes(r,dt):
 +    '''Generate spike train from Poisson rate
 + 
 +    Parameters
 +    ----------
 +    r:  Array of spike rates
 +    dt: Time step
 +    '''
     i = 0     i = 0
     s = 0     s = 0
     spks = []     spks = []
     while True:     while True:
-        s += -np.log(1-np.random.uniform())+        s += np.random.exponential()
         while s>r[i]*dt:         while s>r[i]*dt:
             s -= r[i]*dt             s -= r[i]*dt
Line 16: Line 22:
             continue             continue
         break         break
-    return spks+    return np.array(spks)
 </code> </code>