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)