Program u C++

[cekaonica]

Pravi autor: cekaonica

Istina, ovisno o slucajevima koji se pojavljuju (najgori ili najbolji slucaj), ali prema nekom teoremu iz Teorije vjerojatnosti za jako veliki broj bacanja bi trebao dobiti priblizno jednaki broj pojavljivanja (neznatne razlike) za sve elementarne dogadjaje (polja na ruletu: 0, 00, 1, 2, ..., 36), odnosno aposteriorna vjerojatnost njihovog pojavljivanja trebala bi biti priblizno jednaka apriornoj vjerojatnosti pojavljivanja.

Evo koda koji to potvrdjuje (za 100000 iteracija max. devijacija je oko 4%, za milju je oko 1%, za 10 milja je oko 0,3%; dalje nisam isao jer treba malo vise vremena da se program izvrti, ali za beskonacno iteracija bi trebale biti jednake apriorne i aposteriorne vjerojatnosti)

#include <iostream>
#include <time.h>

using namespace std;

class Rulet
{
public:
Rulet()
{}

int GetGeneratedNumber()
{
int generated;
generated = (int) (rand()/(RAND_MAX+1.) * 38);
return generated;
}
};

int main(void)
{
srand(unsigned(time(NULL)));
int* roulette = new int[38];
int counter, secs;
int number;
int counterTemp;
int red = 0, black = 0, green = 0;
for (int i = 0; i < 38; i++)
{
roulette[i] = 0;
}
do
{
cout << "Upisite broj okretanja ruleta: ";
cin >> counter;
} while (counter <= 0);
counterTemp = counter;
cout << endl;

Rulet* r = new Rulet();
do
{
number = r->GetGeneratedNumber();
roulette[number]++;
} while (--counter != 0);
cout << endl;
for (int i = 0; i < 38; i++)
{
cout << ">> " << i << " --> " << ((double)roulette[i])/counterTemp << endl;
}
cout << endl;
cout << 1./38 << endl << endl;

double dev, devTemp;

dev = ((double)roulette[0])/counterTemp - 1./38;
if (dev < 0)
dev *= -1;
for (int i = 1; i < 38; i++)
{
devTemp = ((double)roulette[i])/counterTemp - 1./38;
if (devTemp < 0)
devTemp *= -1;
if (devTemp > dev)
dev = devTemp;
}

cout << "Najvece odstupanje aposteriorne od apriorne vjerojatnosti: " << dev*38*100 << "%" << endl << endl;
free(r);
return 0;
}