Il link esatto è :
http://manitou.gel.ulaval.ca/~parizeau/doc/a00092.html
Quel sito mostra un insieme di file c++ in cui sono definite varie classi, scritte da Marc Parizeau, un professore di un 'universita francese.
In particolare ha definito una classe NumberGenerator che offre una funzione getGaussianFloat() per avere i numeri con distribuzione gaussiana; la funzione getGaussianFloat genera numeri gaussiani con media zero e deviazione standard pari a 1; se vuoi puoi cambiare questi valori passandoli alla funzione: getGaussianFloat(media, devStand);
Per usare la classe fai cosi: creati un file NumberGenerator.cpp in cui metterai il codice che ti metto dopo, poi nella stessa cartella in cui c'è il NumberGenerator.cpp mettici un altro file NumberGenerator.hpp, in cui inserirai il codice che ti metto dopo. Quindi per usare il generatore di numeri casuali basta includerti il file header NumberGenerator.hpp, creare un oggetto NumberGenerator e invocarci la funzione getGaussianFloat():
****************************************Codice PHP:#include "NumberGenerator.hpp"
....
....
...
//Creo il NumberGenerator
NumberGenerator gen;
double d;
//Genero il numero casuale
d = gen.getGaussianFloat();
CODICE PER IL FILE NumberGenerator.cpp
****************************************
Codice PHP:#include "NumberGenerator.hpp"
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <iostream>
#include <fstream>
using namespace std;
NumberGenerator rand1;
NumberGenerator::NumberGenerator(void)
{
ifstream lStream("rand1.seed");
if(lStream.good()) lStream >> idum;
else idum = time(0);
for(int j=NTAB+7; j >= 0; --j)
{
long k = idum/IQ;
idum = IA*(idum-k*IQ)-IR*k;
if(idum < 0) idum += IM;
if(j < NTAB) iv[j] = idum;
}
iy = iv[0];
iset = 0;
}
NumberGenerator::NumberGenerator(unsigned long inSeed
)
{
idum = inSeed;
for(int j=NTAB+7; j >= 0; --j)
{
long k = idum/IQ;
idum = IA*(idum-k*IQ)-IR*k;
if(idum < 0) idum += IM;
if(j < NTAB) iv[j] = idum;
}
iy = iv[0];
iset = 0;
}
NumberGenerator::~NumberGenerator(void)
{
ofstream lStream("rand1.seed");
if(!lStream.good()) cerr << "Unable to create file \"rand1.seed\"!" << endl;
else lStream << idum;
}
int NumberGenerator::getUniformInteger(int inFirst,
int inLast
)
{
int lNumber = (int) (inFirst + (inLast - inFirst + 1) * getUniformFloat());
if(lNumber > inLast) lNumber = inLast;
return lNumber;
}
double NumberGenerator::getUniformFloat(double inFirst,
double inLast
)
{
double lTmp, lNumber;
long k = idum/IQ;
idum = IA*(idum-k*IQ)-IR*k;
if(idum < 0) idum += IM;
int j = (int) iy/NDIV;
iy = iv[j];
iv[j] = idum;
if((lTmp=AM*iy) > RNMX) lNumber = RNMX;
else lNumber = lTmp;
return inFirst + (inLast - inFirst) * lNumber;
}
double NumberGenerator::getGaussianFloat(double inMean,
double inStdDev
)
{
double fac, rsq, v1, v2;
if(iset == 0) {
do {
v1 = 2.0*getUniformFloat() - 1.0;
v2 = 2.0*getUniformFloat() - 1.0;
rsq = v1*v1+v2*v2;
} while(rsq >= 1.0 || rsq == 0);
fac = sqrt(-2.0*log(rsq)/rsq);
gset = v1*fac;
iset = 1;
return v2*fac*inStdDev + inMean;
} else {
iset = 0;
return gset*inStdDev + inMean;
}
}
****************************************
CODICE PER IL FILE NumberGenerator.hpp
****************************************
Codice PHP:#ifndef Data_NumberGenerator_hpp_
#define Data_NumberGenerator_hpp_
#define IA 16807
#define IM 2147483647
#define AM (1.0/IM)
#define IQ 127773
#define IR 2836
#define NTAB 32
#define NDIV (1+(IM-1)/NTAB)
#define EPS 1.2e-7
#define RNMX (1.0-EPS)
class NumberGenerator
{
public:
NumberGenerator(void);
NumberGenerator(unsigned long inSeed);
~NumberGenerator(void);
int operator()(unsigned long inValue) {return getUniformInteger(0, inValue-1);}
unsigned int getCurrentSeed(void) {return idum;}
bool getUniformBool(void) {return getUniformInteger(0,1);}
int getUniformInteger(int inFirst, int inLast);
double getUniformFloat(double inFirst=0., double inLast=1.0);
double getGaussianFloat(double inMean=0, double inStdDev=1);
private:
long iy;
long iv[NTAB];
long idum;
int iset;
double gset;
};
extern NumberGenerator rand1;
#endif
Rispondi quotando