Skip to Content

simulation de Monte Carlo avec le module Numpy

  1. import numpy as np
  2.  
  3. x1 = np.random.normal(0, stv, N)
  4. y1 = np.random.normal(0, stv, N)
  5. x2 = np.random.normal(100,stv, N)
  6. y2 = np.random.normal(0, stv, N)
  7. x3 = np.random.normal(0, stv, N)
  8. y3 = np.random.normal(100, stv, N)
  9. x4 = np.random.normal(100, stv, N)
  10. y4 = np.random.normal(100, stv, N)
  11.  
  12. # calcul des surfaces
  13. surfaces = np.sqrt((x2-x1)**2+(y2-y1)**2)*np.sqrt((x3-x1)**2+(y3-y1)**2)/2+np.sqrt((x4-x2)**2+(y4-y2)**2)*np.sqrt((x4-x3)**2+(y4-y3)**2)/2
  14.  
  15. # résultats d'une des simulations
  16. print mu, sigma
  17. 9999.1573234 199.631782423
  18. # résultats en tenant compte de la rêgle des chiffres significatifs
  19. print "Resultats :"
  20. print "Moyenne = ","%.f" % mu
  21. print "déviation standard = ", "%.f" % sigma
  22.  
  23. Resultats :
  24. Moyenne = 10000
  25. déviation standard = 200
  26.