Monte-Carlo-Methode in Python

Habe ich versucht, die Python verwenden um ein Skript zu erstellen, dass mir erlaubt die Generierung von großen Nummern der Punkte, die für den Einsatz in der Monte-Carlo-Methode zu berechnen, eine Schätzung Pi. Das Skript habe ich bisher diese:

import math
import random
random.seed()

n = 10000

for i in range(n):
    x = random.random()
    y = random.random()
    z = (x,y)

    if x**2+y**2 <= 1:
        print z
    else:
        del z

So weit, ich bin in der Lage, zum generieren von alle die Punkte, die ich brauche, aber das, was ich möchte zu bekommen, ist die Anzahl der Punkte, die erzeugt wird, wenn der Skript-Ausführung für die Verwendung in einem späteren Berechnung. Ich bin nicht auf der Suche für unglaublich genaue Ergebnisse, gerade gut genug zu schätzen. Irgendwelche Vorschläge würde sehr geschätzt werden.

  • Wollen Sie zählen, wie viele random-Paare sind innerhalb des Kreises? Wenn das der Fall ist, verwenden Sie einfach einen counter...
InformationsquelleAutor Matt Johnson | 2012-11-19
  1. 15

    Wenn Sie dabei jede Art von heavy-duty-numerische Berechnung, wenn man bedenkt das lernen numpy. Dein problem ist im wesentlichen ein one-linear, mit einer numpy setup:

    import numpy as np

N   = 10000
pts = np.random.random((N,2))

# Select the points according to your condition
idx = (pts**2).sum(axis=1)  < 1.0
print pts[idx], idx.sum()

    Geben:

    [[ 0.61255615  0.44319463]
 [ 0.48214768  0.69960483]
 [ 0.04735956  0.18509277]
 ..., 
 [ 0.37543094  0.2858077 ]
 [ 0.43304577  0.45903071]
 [ 0.30838206  0.45977162]], 7854

    Die Letzte Zahl ist die Anzahl der Ereignisse gezählt, d.h. die Zählung der Punkte, deren radius kleiner als eins.

    InformationsquelleAutor Hooked
  2. 2

    Nicht sicher, ob dies ist, was Sie suchen, aber Sie können ausführen enumerate auf range und bekommen die position in Ihrem iteration:

    In [1]: for index, i in enumerate(xrange(10, 15)):
   ...:     print index + 1, i
   ...:
   ...:
1 10
2 11
3 12
4 13
5 14

    In diesem Fall index + 1 darstellen würde der aktuellen Punkt erstellt (index selbst wäre die Gesamtzahl der Punkte erstellt, die zu Beginn einer beliebigen iteration). Auch, wenn Sie mit Python 2.x, xrange ist in der Regel besser für diese Art von Iterationen, wie es lädt nicht die gesamte Liste in den Speicher, sondern greift auf einer as-needed basis.

    • +1 für die Erwähnung aufzählen, wusste das nicht vorher.
    • Kein problem - ich erinnere mich, es für die erste Zeit und denken, 'Ah, dachte, es war ein Weg, das zu tun.' Spaß haben es zu benutzen!
    InformationsquelleAutor RocketDonkey
  3. 1

    Fügen Sie einfach trifft variable vor der Schleife eine Initialisierung auf 0 und in deinem if-Anweisung Inkrement trifft von einem.

    Schließlich können Sie berechnen, PI-Wert mit trifft und n.

    import math
import random
random.seed()

n = 10000
hits = 0  # initialize hits with 0

for i in range(n):
    x = random.random()
    y = random.random()
    z = (x,y)

    if x**2+y**2 <= 1:
        hits += 1
    else:
        del z

# use hits and n to compute PI
    InformationsquelleAutor pm007
Schreibe einen Kommentar