Discussion :

[integrale]

Bonjour

J' essai de faire un truc très simple avec python: je simule un lancer de pièces un certain nombre de fois, et je m' intéresse au résultat de l’expérience. Mon propos est d' afficher le nombre de fois (et ultérieurement la fréquence) que l' issue de l' expérience correspond au paramètre défini (pile, face).

Bien evidemment je me ramene a un tirage aléatoire parmi les deux premiers entiers (0 et 1) ...et cela me donen approximativement ce code

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
 
def lancerIssue(e):
		import random
		from random import randint
		R=[]
		X=0
		for i in range(1,4):
			X=random.randint(0,1)
			R=R+[X]
		print(R)
		s=0
		for j in range(len(R)):
			if e==R[j]:
				s=s+R.count(e)
		print(s)

e prend ses valeurs dans {0,1} pour {pile, face}

mais voila mon programme ne fait pas exactement ce que je voulais...si j' obtiens pile deux fois, le résultat qui s' affiche est 4 (cela semble logique , car je parcours le tableau du résultat et j' affiche le nombre de fois que je rencontre un élément du tableau).....voici une sortie d' ecran de l' execution de mon programme:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
 
 
>>> lancerIssue(0)
[0, 1, 0]
4
>>> lancerIssue(1)
[0, 0, 0]
0
>>> lancerIssue(1)
[0, 0, 1]
1
>>> lancerIssue(1)
[1, 0, 1]
4
>>> lancerIssue(1)
[0, 1, 1]
4
>>> lancerIssue(1)
[1, 0, 0]
1
>>> lancerIssue(1)
[0, 0, 1]
1
>>> lancerIssue(1)
[0, 0, 1]
1
>>> lancerIssue(1)
[0, 1, 1]
4
>>> lancerIssue(1)
[0, 1, 0]
1
>>> lancerIssue(1)
[0, 0, 1]
1
>>> lancerIssue(1)
[1, 0, 1]
4
>>> lancerIssue(1)
[1, 1, 1]
9

Je comprends pourquoi cela ne fonctionne pas....mais j' ai du mal a programmer la sortie du résultat attendu...Merci de me suggérer une indication.


Merci

[integrale]

Manifestement, selon la structure de mon code, je dois demander au programme de retourner la moyenne....

[integrale]

Je crois que j' ai trouvé mieux...il suffit que j' évite la boucle:
Et cela me donne cela:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
 
 
def lancerIssue(e):
		import random
		from random import randrange
		from random import randint
		R=[]
		X=0
		for i in range(1,4):
			X=random.randint(0,1)
			R=R+[X]
		print(R)
		s=0
		s=s+R.count(e)
		print(s)

[Julien N]

Salut,

Tu n'es pas obligé de "déclarer" ta variable avant de l'utiliser. Ansi ta liste de résultat peut être créée directement. Tu peux également faire une "liste compréhension", ce qui allège le code (attention, le mieux reste le plus simple, toujours):

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
def lancerIssue(e):
    import random
    from random import randrange
    from random import randint
 
    R = [random.randint(0,1) for i in range(3)]
    s = R.count(e)
 
    print 'Tirage: ',R
    print 'Occurence %i:' % (e),s
 
lancerIssue(1)

Je ne sais pas ce que tu comptes faire par la suite avec cette fonction, mais méfie toi des tirages "pseudo" aléatoires.

Les fonctions comme random (il me semble tout du moins) sont basées sur l'horloge du PC. Tu cherches à tirer rapidement des nombres, si rapidement que d'un tirage à l'autre le temps CPU ne varie presque pas. Tu tire alors des séries de 0 et 1 et le résultat n'est plus vraiment aléatoire.

Tu peux remarquer ce problème en faisant un large tirage (une vingtaine suffit).

Il existe de nombreuses infos sur le sujet sur le net. Renseigne toi au cas où si tu veux mon avis.


Ju

[dividee]

Je ne sais pas ce que tu comptes faire par la suite avec cette fonction, mais méfie toi des tirages "pseudo" aléatoires.

Les fonctions comme random (il me semble tout du moins) sont basées sur l'horloge du PC. Tu cherches à tirer rapidement des nombres, si rapidement que d'un tirage à l'autre le temps CPU ne varie presque pas. Tu tire alors des séries de 0 et 1 et le résultat n'est plus vraiment aléatoire.

Garder à l'esprit qu'il s'agit de tirages pseudo-aléatoires est un bon conseil, mais tu te trompes si tu penses qu'ils sont basés sur l'horloge du PC. Il s'agit simplement d'une séquence périodique d'états (çàd de nombres) où chaque état est entièrement déterminé par l'état précédent. Quand un nombre est demandé au générateur (avec randint par exemple), il avance la séquence et génère une nombre à partir de l'état résultant. Toute la complexité de l'affaire est que la séquence de nombres résultante ait les caractéristiques statistiques d'une distribution aléatoire uniforme, et que la période du générateur soit très longue (2**19937-1 dans le cas de Python).

Donc, peu importe s'il se passe très peu de temps entre deux tirages, la séquence de nombres générée est toujours la même. La seule chose qui change d'une exécution à l'autre, c'est la position initiale dans la séquence, et c'est là que peut intervenir l'horloge du PC.

[tyrtamos]

Bonjour,

On peut utiliser random.seed(x) pour repartir dans une nouvelle séquence aléatoire (sans argument, il reprend le temps de l'ordinateur).

Sinon, il existe d'autres algorithmes pour générer de bonnes séquences aléatoires. En cas de besoin, j'utilise Blum Blum Shub (https://fr.wikipedia.org/wiki/Blum_Blum_Shub) qui est plus lent mais qui est très bon, au point d'être reconnu en cryptographie. Il n'est d'ailleurs pas très compliqué à programmer en Python (mes essais ici: http://python.jpvweb.com/mesrecettes...id=genalea_bbs).

[Julien N]

Merci dividee pour la précision


J'avais déjà entendu parlé de la méthode par division de nombres premiers très grands, sans jamais avoir vu d'exemple concret.