salut à tous voila je voudrais savoir si c'est possible de calculer pi avec la methode de monte carlo avec un cercle entier
je m'explique , sur le lien dans les cours , on voit que l'algo utlise q'une portion de disque mais je voudrais determiner sans multiplier les tirDansLaPortion/Tir par 4
est ce possible ?
http://exo7.emath.fr/cours/livre-algorithmes.pdf => page 10
Bonjour,
Et bien oui c'est même très simple, tu fais le même dessin mais avec un cercle entier pour des coordonnées x et y allant de -1 à 1, si x2 + y2 lors du tirage, le point est sur le disque. Comme l'aire du disque est
tu as une approximation de
avec "nb de pts sur le disque" / "nb de tirages".
Pas exactement nb de pts sur le disque" / "nb de tirages < 1 donc comment tu peux obtenir pi?
La mesure du grand carré étant 4, il faut encore multiplier par 4 le résultat obtenu
oui donc obligé de multiplier part 4 , donc impossible a faire avec le ccercle entier ,c'est ce que je veux
algo de base pour UNE PORTION DE DISQUE:
import random
def hasard():
Tir=0
NbTirDansleDsique=0
while(Tir<1000000):
x=random.random()
y=random.random()
if((x)*(x)+(y)*(y)<=1):
NbTirDansleDsique=NbTirDansleDsique+1
Tir=Tir+1
Monpi=NbTirDansleDsique/Tir*4
print("Valeur experimentale e Pi :%0.3f"%Monpi)
if __name__ == '__main__':
hasard()
En fait je ne vois pas bien le problème puisque de toute façon faire avec un cercle ou un n-ième de cercle revient exactement au même, ça ne change ni le principe de l'algorithme ni la précision.
En plus comme l'a très justement fait remarquer lionel52 : nb de pts sur le disque" / "nb de tirages < 1. Donc le "sans multiplier les tirDansLaPortion/Tir par 4 " ça va pas trop être possible... et je ne crois pas que les machines de nos jours aient trop de difficultés à multiplier par 4.
ma question était juste de savoir si je peux faire sans multiplier par 4 , mon ordi peut le faire bien sûr le pb n'était pas là , donc => impossible c'est bien ça ?avec cercle de centre x= 0 et y= 0 pas possible de faire avec plus de pi /4
Vous devez être membre accéder à ce service...
Pas encore inscrit ?
1 compte par personne, multi-compte interdit !
Ou identifiez-vous :