Discussion :

[mic0909]

Bonjour,
je suis étudiant en 2ème année de sciences physiques.
Au dernier cours on a parlé de distribution gaussienne, pour moi ca était trop vite et je n'ai pas su comprendre le programme qu'on a écrit....

Voici le 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
#include <iostream>
#include <cstdlib>
#include <cmath>
 
using namespace std;
 
double uniform (float a, float b) {
	return a+(float) rand()/RAND_MAX*fabs(b-a);  
}
 
double expo (double a) {
	double x = uniform(0,1);
	double R;
	while (x==1) {
		x = uniform(0,1);
	}
 
	R = -a*log(1-x);
	return R;
}
 
int main () {
	int N = 50;
	srand (15);
	for (int i = 0; i<N; i++) {
	cout << rand () << endl;	
	}
 
	return 0;
}

En un premier temps, si je comprends bien, on définit uniform, qui prend des nombres aléatoires et qui les distribue de manière uniforme. Par contre pour le expo j'ai vraiment du mal à voir ce qui se passe....

Merci d'avance pour votre aide !

[Davidbrcz]

1) utilise la balise code (symbole# dans l’éditeur) au lieu du gras.
2) Tu ne bloques pas sur du C++ mais sur des maths. Ton sujet aurait plus sa place sur ce forum
3) Mais comme je suis gentil je vais répondre.

uniform devrait te donner un nombre aléatoire entre a et b selon une loi .. uniforme, ie avec une densité de probabilité constante. En pratique, il y quelques problèmes

Pour transformer cette distribution en une distribution normale, il y a plusieurs méthodes.

• Une méthode peut être la transformé de Box-Muller.
• Toi, tu utilises le lemme suivant : si U suit une loi uniforme sur [0,1], alors X= F^(-1)(U) suit la loi de fonction de distribution F. On veut F une exponentielle, donc l'inverse de F c'est du log. Voir ici page 6 pour une démonstration

4) N'utilise JAMAIS du code que tu as écrit toi même sur ce genre de chose. Utilise std::normal_distribuion (introduit en C++11)