🚀 go-pugleaf

On 6 Oct 97 21:22:10 GMT, "yacko" <yacko@priv.onet.pl> wrote:

>
>	Czy ktoœ mo¿e mi pomóc? Próbuje stworzyc algorytm generuj¹cy liczby
>pseudolosowe odbiegaj¹cy dosc znacznie od tego uzywanego w narzedziach
>Borlanda.
>	Zalezy mi na tym by dla duzej serii dawa³ rozklad plaski, ale zdarzalyby
>mu sie wygenerowac czesciej serie np zlozone z duzych wartosci, malych itd.
>Po prostu "mniej doskona³y".
>	zainteresowanych prosze o kontakt
>-- 
Sprawa jest dosc prosta:
- musisz miec zadany rozklad zmiennej losowej opisany funkcja (moze
byc interpolacja kilku punktow). 
- rozklad zmiennej losowej musi spelniac warunek, zeby calka od 0 do
nieskonczonosci = 1 
- jezeli nie jest =1 to funkcje dzielisz przez tyle ile jest rowna
(tak, zeby w koncu calka wyszla 1)
- wyznaczasz dystrybuante rozkladu, ze wzoru:
F(x)=calka od 0 do x z f(x) dx
gdzie f(x) - rozklad, F(x) - dystrybuanta
- jezeli wszystko dobrze zrobiles to F(0) powinno sie rownac 0, a
F(nieskonczonosc) = 1.
- znajdujesz funkcje odwrotna do dystrybuanty. Jej dziedzina powinna
nalezec do przedzialu <0,1>
- piszesz procedure losujaca:
...
float moja_liczba_losowa() 
{
  float a=random()  // procedura generujaca liczby o rozkladzie
rownomiernym z przedzialu <0,1) (moze byc standardowa z C/Pascala)
  float moja_losowa=funkcja_odwrotna_do_dystrybuanty(a);
  return moja_losowa;
}

Jezeli nic nie pomylilem, to powinno byc OK :-)

Marek Matula
marekm@zeus.polsl.gliwice.pl