Mersenne Twister z 2 granicami

[oneiro]

Witam,

Potrzebuję zrobić generator liczb pseudolosowych albo mający ograniczenia z dołu i góry. Maksymalny przedział liczy to:

0 <= RAND <= 1

natomiast w większości przypadków będzie to

A [<=,<] RAND [<=,<] B

gdzie A<B oraz należą do przedziału (0,1). RAND oczywiście jest liczbą dziesiętną.
Ściągnąłem kod Mersenne Twister i mam tam opcje generowanie liczb z przedziałów
[0,1]
(0,1]
[0,1)
(0,1)

Mogę ograniczać liczby poprzez IF, ale mija się to z celem (dla dostatecznie wąskich przedziałów więcej czasu zajmie mi generowanie liczby niż właściwe obliczenia)

[Mironas]

Na początku programu (np. w OnCreate formy) daj:

KOD cpp:   �  UKRYJ  �

KOD cpp:   �  POKAŻ  �

  randomize();  // Uruchomienie generatora liczb losowych
 

Pojedynczą liczbę generujesz:

KOD cpp:   �  UKRYJ  �

KOD cpp:   �  POKAŻ  �

  int x = random(max);  // generuje liczbę całkowitą z zakresu 0÷(max-1)
 

I parę wskazówek:
Jeśli potrzebujesz wygenerować liczbę z zakresu 10÷20 to najpierw generujesz z zakresu 0÷10 a potem dodajesz 10.
Jeśli potrzebujesz liczbę zmiennoprzecinkową (float) z zakresu 0.0÷1.0 to generujesz liczbę w zakresie 0÷1000, zapisujesz do zmiennej float a na koniec dzielisz przez 1000.0

Powyższe wskazówki powinny Ci wystarczyć do przygotowania takiego generatora jaki potrzebujesz.

[Cyfrowy Baron]

Jeśli potrzebujesz wygenerować liczbę z zakresu 10÷20 to najpierw generujesz z zakresu 0÷10 a potem dodajesz 10.

Można tak, ale można również skorzystać z funkcji RandomRange(int AFrom, int ATo) - gdzie zmienna AFrom określa początkowy zakres losowania (włączenie), a ATo końcowy zakres losowania, przy czym jeżeli podamy np. 20 to zakres losowania kończy się na 19. Generatorem liczb dla RandomRange jest Randomise() - pisany z wielkiej litery nie z małej. Generator uruchamiamy tylko raz np. na starcie aplikacji w zdarzeniu OnCreate lub w konstruktorze klasy okna głównego. Jeżeli losowanie odbywa się na innym formularzu niż główny to generator i tak uruchamiamy raz na formularzu głównym. Generatora nie da się zatrzymać.

[oneiro]

Dziękuję Wam, na podstawie waszych postów znalazłem to czego szukałem (widać, źle wcześniej szukałem)

KOD cpp:   �  UKRYJ  �

KOD cpp:   �  POKAŻ  �

//To mi generuje od 0.0 do 1.0 włącznie:
float r = (float)rand()/(float)RAND_MAX;

//To generuje od 0.0 to górnej granicy X:
float r2 = (float)rand()/((float)RAND_MAX/X);

//To generuje od dolnej granicy LO do górnej HI
float r3 = LO + (float)rand()/((float)RAND_MAX/(HI-LO));
 

Znalazłem też fajną bibliotekę Boost.Random

KOD cpp:   �  UKRYJ  �

KOD cpp:   �  POKAŻ  �

  boost::mt19937 rng;
  boost::uniform_real<float> u(min, max);
  boost::variate_generator<boost::mt19937&, boost::uniform_real<float> > gen(rng, u);
  double wynik=gen();
 

Zależmy mi na MT (lub WELL - Well Equidistributed Long-period Linear), dlatego też chyba skorzystam z Boost.Random

[Cyfrowy Baron]

Funkcja rand() wymaga uruchomienia generatora srand. Jeżeli nie uruchomisz generatora to program zawsze będzie losował te same liczby:

KOD cpp:   �  UKRYJ  �

KOD cpp:   �  POKAŻ  �

#include <stdlib.h>     /* srand, rand */
#include <time.h>       /* time */

 srand ( time(NULL) );

Makro RAND_MAX to po prostu wartość: 32767, dlatego ten sam efekt uzyskasz za pomocą funkcji Random:

KOD cpp:   �  UKRYJ  �

KOD cpp:   �  POKAŻ  �

Randomize(); /* uruchomienie generatora */

float x = 32767;
float r = Random(x) / x;

//To generuje od dolnej granicy LO do górnej HI

Podobnie jak to:

KOD cpp:   �  UKRYJ  �

KOD cpp:   �  POKAŻ  �

Randomize(); /* uruchomienie generatora */

float x = 32767;
float r = RandomRange(x*2, x*4) / x;

w powyższym przykładzie zakres to: dolny = 2 (x*2); górny = 4 (x*4). Po prostu pierwszy argument funkcji to x=32767 pomnożony przez wartość dolną, a drugi adekwatnie pomnożony przez wartość górną. Wartość x nie musi się równać dokładnie 32767, wystarczy, że będzie większa od losowanego zakresu. Wszystko zależy od tego jak dokładna ma być wartość dziesiętna.