Differences
This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| technical:poisson_spike_train [2020/10/13 19:36] – 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> | ||
| - | # generate spike train | + | def gen_spikes(r, |
| - | def gen_spikes(r, | + | ''' |
| + | |||
| + | 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> | while s> | ||
| s -= r[i]*dt | s -= r[i]*dt | ||
| Line 16: | Line 23: | ||
| continue | continue | ||
| break | break | ||
| - | return spks | + | return |
| </ | </ | ||