Gledam i cudim se, sta li je ovde problem :)
Pola njih prica o generisaju slucajne promenljive sa uniformnom raspodelom, a pola o nekim transformacijama :).
E sada koliko sam ja shvatio problem je u tome da se napravi uniformni generator tacaka unutar nekog kruga.
Da budemo jos malo precizni, jel tako da je doupsteno koristiti neku funkciju koja daje uniformnu raspodelu unutar segmenta [0,1]. Takooooooooo jeeeeeeeee.
E da je u pitanu kvadrat, odavno bi se to resilo lako, lepo x=rand(); y=rand(); radi posao, ali mi ne zelimo kvadrat, zelimo krug.
Neko bi rekao: pa ajde daj uniformno po kvadratu, pa ce i na nekovom njegovom podskupu biti isto uniformno, pa cemo da secemo ako je izvan kruga, e to nece ici, jer ljudi zele algoritam u konacno mnogo koraka, bez petlji, jer ako u petlji cekamo da tacka padne unutar kruga, to moze da se desi posle bas puno koraka (jelte, ne mozemo da garantujemo da ce to da se desi bas u prvih 10 koraka), neko kakve je srece moze da nikada ne doceka, nema veze sto je uniformna raspodela, ona je uniformna u beskonacnosti, a konacno je diskretna :)))).
Dobro, ajde ne mora tako, sto ne bi iskoristili polarni kordinatni sistem, neka jedna promenjljiva ide od 0 do 2pi, a druga od 0 do 1 dobili bi lepo definisane tacke unutar kruga, ali ovakva raspodela nije bas uniformna, ili jeste :). Naravno tu su sada pale rasprave oko diskretnosti i nediskretnosti. da li je matrica ovakva ili ovakva ...
Nekako deluje logicno da ovako dobijen skup tacaka nije u uniformnoj rapodeli, jer blizu (0,0) ima vise tacaka gustina :)))).
E sada ispada da problem nije bas tako trivijalan, treba naci neko genijalno mapiranje segmenta [0,1] ili pak prozivoda vise segmenata [0,1] na krug. Dakle ajde da vidimo neke peanove krive ili sta vec :).
Dakle nije problem dobiti jednu uniformnu promenljivu iz segmenta [0,1], pitanje je samo da li statisticki zadovoljava uniformnost i na malim uzorcima :), nego je problem bas to, unutar kruga ....
CHUPCKO
Re: Generisanje raspodela