Discussion :

[Isabella83]

Bonjour à tous,
Je chercher à générer des nombres aléatoires dans les bornes d'un intervalle [0 .... mrna_length].
Après quelques recherches sur internet, j'ai fait comme ceci :

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
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>
 
void simulation(int mrna_length, int repeat);
 
int main()
{	
 
	int n_cut = 500 ;
	int repeat = 10;
	simulation(n_cut,repeat);
	return 0;
}
 
void simulation(int mrna_length, int repeat){
	int i,j,k;
	for(i=0;i<=repeat;i++){
                /* deux façons de générer des nombres aléatoires dans les bornes d'un intervalle */
		j = rand()%(mrna_length+1); 
		k= 1 + (int) ((mrna_length) * (rand() / (RAND_MAX + 1.0)));
		printf("valeur de j : %d\n",j);
		printf("valeur de k : %d\n",k);
	}
}

J'ai compilé avec la commande :

Code :

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

gcc -Wall -Wextra test.c -O2 -o programme

Ensuite je fais

Code :

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

./programme

.
Cela me sort bien des nombres aléatoires. Jusque là tout va bien.
Mais je suis très étonnée, car lorsque je fais à nouveau gcc puis ./mon programme, j'ai toujours la même liste de valeurs aléatoires, dans le même ordre etc.
Je ne comprends pas pourquoi ?
c'est surement très bête, mais je débute.
Merci pour votre aide !

[ternel]

C'est la magie de l'informatique.
La fonction rand() ne génère pas des nombres aléatoires.
Elle génère les valeurs successives d'une suite difficilement prédictible. C'est ce qu'on appelle des nombres pseudo-aléatoires.
En anglais, on parle de PRNG: pseudo random number generator/ion.

Il existe une fonction soeur à rand(): srand(seed)Cette fonction initialise la "seed" de la suite générée par rand().

A chaque fois que tu appelles srand(), la suite est réinitialisée.
Pour avoir une suite différente à chaque execution, il suffit d'utiliser une seed distincte à chaque fois.

La technique généralement utilisé, c'est de prendre l'heure actuelle.

Code :

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

srand(time(NULL));

Il faut exécuter cette instruction une seule fois, avant la première fois qu'est appelée rand().

[Isabella83]

D'accord je comprends mieux !
Effectivement en modifiant :

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
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>
 
void simulation(int mrna_length, int repeat);
 
int main()
{	
 
	int n_cut = 500 ;
	int repeat = 10;
	simulation(n_cut,repeat);
	return 0;
}
 
void simulation(int mrna_length, int repeat){
	int i,j,k;
	srand(time(NULL));
	for(i=0;i<=repeat;i++){
		j = rand()%(mrna_length+1);
		k= 1 + (int) ((mrna_length-1) * (rand() / (RAND_MAX + 1.0)));
		printf("valeur de j : %d\n",j);
		printf("valeur de k : %d\n",k);
	}
}

j'obtiens des nombres différents !

[picodev]

Bonjour,


attention, avec ce code tu auras k dans l'intervalle [1, mrna_length] et non [0, mnra_length]. A priori la seconde méthode te donnera des nombres mieux distribués dans cette intervalle que la méthode avec le modulo.

[Bktero]

Ne pas oublier qu'il y a une sur Developpez.com

http://c.developpez.com/faq/?page=Le...res-aleatoires

[Isabella83]

Ne pas oublier qu'il y a une sur Developpez.com

http://c.developpez.com/faq/?page=Le...res-aleatoires

Oui effectivement, merci !

[Sve@r]

La technique généralement utilisé, c'est de prendre l'heure actuelle.

Code :

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

srand(time(NULL));

Bonjour

Même avec l'heure utilisée comme graine, on peut obtenir une suite identique si le même programme est appelé plusieurs fois dans la même seconde. Peu probable certes mais néanmoins possible surtout dans les systèmes Unix/Linux multi utilisateurs (deux users appelant le programme en même temps...)

C'est pourquoi je prends souvent l'habitude de rajouter le pid (n° de processus) dans la graine. Ce n° étant unique à un instant "t", ça réduit les risques...

Code :

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

srand(time(NULL) ^ getpid());

Autre possibilité: prendre l'heure codée en millièmes de secondes...