Advertisement
Guest User

Untitled

a guest
Mar 19th, 2019
55
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.57 KB | None | 0 0
  1. from random import random
  2. from math import log
  3.  
  4. def intensity(t):
  5. if 0<=t and 20>t:
  6. return 8
  7. elif 20<=t and t<40:
  8. return 4
  9. else:
  10. return 8
  11.  
  12. lamd = (8/20.+4/20.+8/20.)/3 # мат ожидание
  13. end_time = 60
  14.  
  15. current_time = 0
  16. events = []
  17.  
  18. while(True):
  19. r = random()
  20. current_time = current_time - 1 / lamd * log(r,2)
  21. if current_time > end_time:
  22. break
  23. r = random()
  24. if r <= intensity(current_time)/lamd:
  25. events.append(current_time)
  26.  
  27. print('Время прихода: ')
  28. for event in events:
  29. print("%.2f" % event)
  30.  
  31. print("Количество событий: ", len(events))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement