[c] génération uniforme de nombres aléatoires
Résolu
leti
-
A DISTRIBUTION DE DÉPART EST UNIFORME, DONC -
A DISTRIBUTION DE DÉPART EST UNIFORME, DONC -
Bonjour,
je ne m'y connais pas du tout en programmation de toute sorte, mais je dois passer quelques mois sur c et j'ai un souci avec les nombres aléatoires.
Je dois lancer un grand nombre de nombres aléatoires réels, et au final, j'ai un plus grand nombres de chiffres tirés à 0. Est-ce que zero est plus probable que les autres chiffres?
j'ai utilisé srand(time(NULL)); et (rand()/(double)RAND_MAX )*(20); et je dois tirer un million de nombres.
merci d'avance
je ne m'y connais pas du tout en programmation de toute sorte, mais je dois passer quelques mois sur c et j'ai un souci avec les nombres aléatoires.
Je dois lancer un grand nombre de nombres aléatoires réels, et au final, j'ai un plus grand nombres de chiffres tirés à 0. Est-ce que zero est plus probable que les autres chiffres?
j'ai utilisé srand(time(NULL)); et (rand()/(double)RAND_MAX )*(20); et je dois tirer un million de nombres.
merci d'avance
A voir également:
- [c] génération uniforme de nombres aléatoires
- Generateur mot de passe - Télécharger - Sécurité
- Generation ryzen - Guide
- Réinitialiser chromecast 1ere génération - Guide
- Code binaire des nombres - Guide
- Nombre de jours entre deux dates excel - Guide
7 réponses
Bonjour,
Non, il n'y a absolument aucune raison que 0 soit restiué plus de fois.
En utilisant : 20.0*rand()/(RAND_MAX+1.0);
Il faudrait voir le programme en entier, peut-être une petite erreur ailleurs...
Non, il n'y a absolument aucune raison que 0 soit restiué plus de fois.
En utilisant : 20.0*rand()/(RAND_MAX+1.0);
Il faudrait voir le programme en entier, peut-être une petite erreur ailleurs...
Salut.
à ma connaissance il n'y a pas de raison que 0 sorte plus souvent.
à lire : https://linux.die.net/man/3/rand
Sinon, comme tu travail sous Linux, il existe un générateur de nombre aléatoire intégré qu'il est possible d'utiliser en lisant le fichier /dev/urand
https://linux.die.net/man/4/random
à ma connaissance il n'y a pas de raison que 0 sorte plus souvent.
à lire : https://linux.die.net/man/3/rand
Sinon, comme tu travail sous Linux, il existe un générateur de nombre aléatoire intégré qu'il est possible d'utiliser en lisant le fichier /dev/urand
https://linux.die.net/man/4/random
rand() est un entier inférieur à RAND_MAX, donc quand tu fais la division ça devient 0...
Pour que ça marche tu dois caster rand()
Remarque : avoir fait le cast sur RAND_MAX aurait du suffire.
La confiance n'exclut pas le contrôle
Pour que ça marche tu dois caster rand()
double alea = (double) rand()/RAND_MAX*20;
Remarque : avoir fait le cast sur RAND_MAX aurait du suffire.
La confiance n'exclut pas le contrôle
En utilisant time(NULL), tu genere tes nombres aléatoire selon un nombre de seconde, donc il se peu que tu ai plusieur fois le meme nombre a la suite si tu en tire plusieur a la suite car, si je ne me trompe pas, le type time_t est un entier.
Cependant il nous faudrai plus d'info concernant les nombres qui se repete (s'il se repete a la suite, regulierement, si c'est toujour 0 etc)
Cependant il nous faudrai plus d'info concernant les nombres qui se repete (s'il se repete a la suite, regulierement, si c'est toujour 0 etc)
Vous n’avez pas trouvé la réponse que vous recherchez ?
Posez votre question
bonjour,
j'ai du me tromper ailleurs alors sans doute... :(
Je dois étudier une diffusion (N=1/(4 pi K t )^3/2*exp(-r²/4Kt) avec comme variables r et t rentré par l'utilisateur. le reste sont des constantes.
La distribution de départ est uniforme, donc je dois tirer au hasard un (r2,theta) par rapport au centre d'un cercle de rayon 20.
La diffusion est regardée à partir d'un point se trouvant à 8 du centre, donc la formule est celle d'al kashi :
r²=r2²+8²-2*8*r2*cos(theta).
enfin je regarde la contribution de chaque anneau autour du point d'observation. C'est à ce moment que j'obtient une irrégularité au niveau de 8, alors que la courbe devrait être lisse.
voilà, j'espère avoir été claire sur la modélisation, voilà ce que j'ai fait. (je rappelle aque c'est mon premier programme, donc soyez indulgents....)
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
//calcul de la diffusion dimension n;
double diff(double r, double t)
{
double D0,delta,N0;
double D,N,R;
int n;
D0=0.1 ;//kpc²/Myr
delta=0.5;
N0=1;
R=1; //GeV
n=2;
D=D0*pow(R,delta);
N=N0/pow(4*M_PI*D*t,(n/2)) *exp (-r*r/(4*D*t));
return N;
}
int main(int argc, char *argv[])
{
srand(time(NULL));
FILE* contribution = NULL;
contribution = fopen("contribution.txt", "w");
double deltat,tobs,dimension;
double rmax,tetamax,Nmax;
double N,r2,r,SN,teta;
double rsol;
int i,numero,j,Nb;
double a;
double tableau [280]; //contribution par rapport au rayon
deltat= atoi (argv[1]); //durée d'étude
tobs= atoi (argv[2]); //temps de l'observateur
//initialisation des paramètres
for (j=1;j<=280;j++)
{
tableau[j]=0;
}
Nmax=0;
tetamax=0;
rmax=0;
rsol=8.; //kpc
Nb=deltat/(pow(10,-4));
//calcul de toutes les SN
for (i=1;i<=Nb;i++)
{ a=0;
SN=(rand()/(double)RAND_MAX )*(deltat); //temps d'explosion de la ieme SN
if(SN<tobs) //observation à un temps SN< tems total d'étude
{
//position de la ieme SN
r2=(rand()/(double)RAND_MAX )*(20);
teta=(rand()/(double)RAND_MAX )*(2*M_PI);
r=sqrt(pow(r2,2)-(2*rsol*r2*cos(teta))+pow(rsol,2));
//calcul de N a la ieme SN
N=diff(r2,SN);
//contribution par rapport au rayon
for (j=1;j<=280;j++)
{
if (r>a && r<=(a+0.1))
{
tableau[j]=tableau[j]+N;
}
a=a+0.1;
}
// recherche de la SN qui contribue le plus
if (N>Nmax)
{
Nmax=N;
tetamax=teta;
rmax=r;
numero=i;
}
}
}
a=0;
for (j=1;j<=279;j++)
{
a=a+0.1;
fprintf(contribution,"%f %f\n",a,tableau[j]);
}
fclose(contribution);
return 0;
}
j'ai du me tromper ailleurs alors sans doute... :(
Je dois étudier une diffusion (N=1/(4 pi K t )^3/2*exp(-r²/4Kt) avec comme variables r et t rentré par l'utilisateur. le reste sont des constantes.
La distribution de départ est uniforme, donc je dois tirer au hasard un (r2,theta) par rapport au centre d'un cercle de rayon 20.
La diffusion est regardée à partir d'un point se trouvant à 8 du centre, donc la formule est celle d'al kashi :
r²=r2²+8²-2*8*r2*cos(theta).
enfin je regarde la contribution de chaque anneau autour du point d'observation. C'est à ce moment que j'obtient une irrégularité au niveau de 8, alors que la courbe devrait être lisse.
voilà, j'espère avoir été claire sur la modélisation, voilà ce que j'ai fait. (je rappelle aque c'est mon premier programme, donc soyez indulgents....)
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
//calcul de la diffusion dimension n;
double diff(double r, double t)
{
double D0,delta,N0;
double D,N,R;
int n;
D0=0.1 ;//kpc²/Myr
delta=0.5;
N0=1;
R=1; //GeV
n=2;
D=D0*pow(R,delta);
N=N0/pow(4*M_PI*D*t,(n/2)) *exp (-r*r/(4*D*t));
return N;
}
int main(int argc, char *argv[])
{
srand(time(NULL));
FILE* contribution = NULL;
contribution = fopen("contribution.txt", "w");
double deltat,tobs,dimension;
double rmax,tetamax,Nmax;
double N,r2,r,SN,teta;
double rsol;
int i,numero,j,Nb;
double a;
double tableau [280]; //contribution par rapport au rayon
deltat= atoi (argv[1]); //durée d'étude
tobs= atoi (argv[2]); //temps de l'observateur
//initialisation des paramètres
for (j=1;j<=280;j++)
{
tableau[j]=0;
}
Nmax=0;
tetamax=0;
rmax=0;
rsol=8.; //kpc
Nb=deltat/(pow(10,-4));
//calcul de toutes les SN
for (i=1;i<=Nb;i++)
{ a=0;
SN=(rand()/(double)RAND_MAX )*(deltat); //temps d'explosion de la ieme SN
if(SN<tobs) //observation à un temps SN< tems total d'étude
{
//position de la ieme SN
r2=(rand()/(double)RAND_MAX )*(20);
teta=(rand()/(double)RAND_MAX )*(2*M_PI);
r=sqrt(pow(r2,2)-(2*rsol*r2*cos(teta))+pow(rsol,2));
//calcul de N a la ieme SN
N=diff(r2,SN);
//contribution par rapport au rayon
for (j=1;j<=280;j++)
{
if (r>a && r<=(a+0.1))
{
tableau[j]=tableau[j]+N;
}
a=a+0.1;
}
// recherche de la SN qui contribue le plus
if (N>Nmax)
{
Nmax=N;
tetamax=teta;
rmax=r;
numero=i;
}
}
}
a=0;
for (j=1;j<=279;j++)
{
a=a+0.1;
fprintf(contribution,"%f %f\n",a,tableau[j]);
}
fclose(contribution);
return 0;
}
Bon, je ne vois rien qui pourrait mettre la merde.
juste quelques remarque : les tableaux en C vont de 0 à N-1 et non de 0 à N, du coup tu risques des erreurs (c'est une grosse erreur en C).
attention à la division d'entiers dans diff() (n/2). Lorsque tu veux utiliser des doubles, utilise un point à la fin : 3./2. te donnera 1.5, c'est une bonne méthode à prendre. Aussi, au lieu d'écrire pow(10,-4), écrit simplement 1e-4.
pour afficher les double j'utilise %e, il me semble que %f c'est pour les float.
ton programme n'est pas simple pour quelqu'un qui n'est pas du domaine. explique nous où ça merde dans ton programme. Il est possible que tu es une erreur d'algorithme (et pas de programmation).
Utilise les balises de code pour mettre ton source :
La vrai soumission c'est quand les esclaves s'inquiètent du cours du coton.
Char Snipeur
juste quelques remarque : les tableaux en C vont de 0 à N-1 et non de 0 à N, du coup tu risques des erreurs (c'est une grosse erreur en C).
attention à la division d'entiers dans diff() (n/2). Lorsque tu veux utiliser des doubles, utilise un point à la fin : 3./2. te donnera 1.5, c'est une bonne méthode à prendre. Aussi, au lieu d'écrire pow(10,-4), écrit simplement 1e-4.
pour afficher les double j'utilise %e, il me semble que %f c'est pour les float.
ton programme n'est pas simple pour quelqu'un qui n'est pas du domaine. explique nous où ça merde dans ton programme. Il est possible que tu es une erreur d'algorithme (et pas de programmation).
Utilise les balises de code pour mettre ton source :
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
//calcul de la diffusion dimension n;
double diff(double r, double t)
{
double D0,delta,N0;
double D,N,R;
int n;
D0=0.1 ;//kpc²/Myr
delta=0.5;
N0=1;
R=1; //GeV
n=2;
D=D0*pow(R,delta);
N=N0/pow(4*M_PI*D*t,(n/2)) *exp (-r*r/(4*D*t));
return N;
}
int main(int argc, char *argv[])
{
srand(time(NULL));
FILE* contribution = NULL;
contribution = fopen("contribution.txt", "w");
double deltat,tobs,dimension;
double rmax,tetamax,Nmax;
double N,r2,r,SN,teta;
double rsol;
int i,numero,j,Nb;
double a;
double tableau [280]; //contribution par rapport au rayon
deltat= atoi (argv[1]); //durée d'étude
tobs= atoi (argv[2]); //temps de l'observateur
//initialisation des paramètres
for (j=1;j<=280;j++)
{
tableau[j]=0;
}
Nmax=0;
tetamax=0;
rmax=0;
rsol=8.; //kpc
Nb=deltat/(pow(10,-4));
//calcul de toutes les SN
for (i=1;i<=Nb;i++)
{
a=0;
SN=(rand()/(double)RAND_MAX )*(deltat); //temps d'explosion de la ieme SN
if (SN<tobs) //observation à un temps SN< tems total d'étude
{
//position de la ieme SN
r2=(rand()/(double)RAND_MAX )*(20);
teta=(rand()/(double)RAND_MAX )*(2*M_PI);
r=sqrt(pow(r2,2)-(2*rsol*r2*cos(teta))+pow(rsol,2));
//calcul de N a la ieme SN
N=diff(r2,SN);
//contribution par rapport au rayon
for (j=1;j<=280;j++)
{
if (r>a && r<=(a+0.1))
{
tableau[j]=tableau[j]+N;
}
a=a+0.1;
}
// recherche de la SN qui contribue le plus
if (N>Nmax)
{
Nmax=N;
tetamax=teta;
rmax=r;
numero=i;
}
}
}
a=0;
for (j=1;j<=279;j++)
{
a=a+0.1;
fprintf(contribution,"%f %f\n",a,tableau[j]);
}
fclose(contribution);
return 0;
}
La vrai soumission c'est quand les esclaves s'inquiètent du cours du coton.
Char Snipeur
merci déjà pour toutes ces précisions.
Je ne comprends pas exactement ce que tu entend par erreur d'algorithme...:)
mon programme fonctionne, je peux tracer une courbe une fois qu'il a terminé. Seulement cette courbe n'est pas lisse, elle a un pic en 8. Lorsque dans la définition de r (r²=r2²+8²+2*8*r2*cos), j'enlève le 8² (juste par test, la fonction est exacte), l'irrégularité disparait, ce qui me poussait à croire à un problème avec rand.
Lorsque je tire deux nombres aléatoires et que je les fait agir ensemble dans une fonction, peut-être y a t il un effet qui fait cette augmentation en 8?
en tous cas, merci déjà, je suis contente de n'avoir pas fait de grosse erreur dans le code!!
Je ne comprends pas exactement ce que tu entend par erreur d'algorithme...:)
mon programme fonctionne, je peux tracer une courbe une fois qu'il a terminé. Seulement cette courbe n'est pas lisse, elle a un pic en 8. Lorsque dans la définition de r (r²=r2²+8²+2*8*r2*cos), j'enlève le 8² (juste par test, la fonction est exacte), l'irrégularité disparait, ce qui me poussait à croire à un problème avec rand.
Lorsque je tire deux nombres aléatoires et que je les fait agir ensemble dans une fonction, peut-être y a t il un effet qui fait cette augmentation en 8?
en tous cas, merci déjà, je suis contente de n'avoir pas fait de grosse erreur dans le code!!
