Differences
This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| technical:poisson_spike_train [2021/09/21 09:05] – chunchung | technical: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, | + | def gen_spikes(r, |
| - | ''' | + | ''' |
| Parameters | Parameters | ||
| ---------- | ---------- | ||
| - | r: | + | r: Array of spike rates |
| - | dt: time step | + | dt: |
| + | 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> | while s> | ||
| s -= r[i]*dt | s -= r[i]*dt | ||
| Line 22: | Line 23: | ||
| continue | continue | ||
| break | break | ||
| - | return spks | + | return |
| </ | </ | ||