Translations:Gerador Aleatório para Distribuição de Poisson/4/en
Ir para navegação
Ir para pesquisar
def pseudo_poisson(alpha, size=1):
"""
Poisson generator from uniform generator
"""
poisson = []
for _ in range(size):
# Sets seed based on the decimal portion of the current system clock
t = time.perf_counter()
seed = int(10**9*float(str(t-int(t))[0:]))
U = pseudo_uniform(seed=seed, size=5*alpha)
X,P,i = 0,1,0
while P >= np.exp(-alpha):
P = U[i]*P
X += 1
i += 1
poisson.append(X)
return np.array(poisson)