abclinuxu.cz AbcLinuxu.cz itbiz.cz ITBiz.cz HDmag.cz HDmag.cz abcprace.cz AbcPráce.cz
AbcLinuxu hledá autory!
Inzerujte na AbcPráce.cz od 950 Kč
Rozšířené hledání
×
dnes 16:00 | Nová verze

Byl vydán Mozilla Firefox 51.0. Z novinek lze upozornit například na upozorňování na přihlašování přes nešifrované spojení (HTTP), podporu pro přehrávání bezeztrátového formátu FLAC nebo podporu WebGL 2. Podrobné informace v poznámkách k vydání a na stránce věnované vývojářům. Řešeny jsou také bezpečnostní chyby.

Ladislav Hagara | Komentářů: 0
včera 17:25 | IT novinky

Do prodeje (Farnell) se dostal jednodeskový počítač Tinker Board (unboxing). Jedná se o konkurenci Raspberry Pi 3 od společnosti Asus. Porovnání (jpg) těchto počítačů například na CNXSoft. Cena Tinker Boardu je 55 £.

Ladislav Hagara | Komentářů: 14
včera 14:44 | Zajímavý projekt

Byla zveřejněna pravidla hackerské soutěže Pwn2Own 2017, jež proběhne od 15. do 17. března v rámci bezpečnostní konference CanSecWes ve Vancouveru. Soutěžit se bude o více než milion dolarů v pěti kategoriích. Letos se bude útočit i na Ubuntu. Jedná se již o 10. ročník této soutěže.

Ladislav Hagara | Komentářů: 2
včera 13:33 | Nová verze

Po sedmi měsících vývoje od vydání verze 5.7 byla vydána verze 5.8 (YouTube) toolkitu Qt. Z novinek lze zmínit například Qt Lite pro vestavěná zařízení. Nově jsou plně podporovány moduly Qt Wayland Compositor (YouTube) a Qt SCXML (YouTube). Současně byla vydána verze 4.2.1 integrovaného vývojového prostředí (IDE) Qt Creator.

Ladislav Hagara | Komentářů: 1
včera 11:52 | Pozvánky

Lednový Prague Containers Meetup se koná ve čtvrtek 26. ledna 2017 od 18:00 v Apiary, Pernerova 49, Praha 8. Přijďte se podívat na přednášky o Enterprise Kubernetes a Jenkins as a code.

little-drunk-jesus | Komentářů: 0
včera 11:40 | Pozvánky

Program letošního ročníku konference Prague PostgreSQL Developer Days, která se koná již 15. a 16. února 2017 na ČVUT FIT, Thákurova 9, Praha 6, byl dnes zveřejněn. Najdete ho na stránkách konference včetně anotací přednášek a školení. Registrace na konferenci bude otevřena zítra (24. ledna) v brzkých odpoledních hodinách.

TomasVondra | Komentářů: 0
22.1. 02:20 | Zajímavý článek

David Revoy, autor open source webového komiksu Pepper&Carrot nebo portrétu GNU/Linuxu, upozorňuje na svém blogu, že nový Inkscape 0.92 rozbíjí dokumenty vytvořené v předchozích verzích Inkscape. Problém by měl být vyřešen v Inkscape 0.92.2 [reddit].

Ladislav Hagara | Komentářů: 0
22.1. 02:02 | Komunita

Øyvind Kolås, hlavní vývojář grafických knihoven GEGL a babl, které využívá grafický program GIMP, žádá o podporu na Patreonu. Díky ní bude moci pracovat na vývoji na plný úvazek. Milník 1000 $, který by stačil na holé přežití, se již téměř podařilo vybrat, dalším cílem je dosažení 2500 $, které mu umožní běžně fungovat ve společnosti.

xkomczax | Komentářů: 12
21.1. 23:54 | Pozvánky

DevConf.cz 2017, již devátý ročník jedné z největších akcí zaměřených na Linux a open source ve střední Evropě, proběhne od pátku 27. ledna do neděle 29. ledna v prostorách Fakulty informačních technologií Vysokého učení technického v Brně. Na programu je celá řada zajímavých přednášek a workshopů. Letos je povinná registrace.

Ladislav Hagara | Komentářů: 0
21.1. 22:11 | Nová verze

Byla vydána verze 1.0.0 emulátoru terminálu Terminology postaveného nad EFL (Enlightenment Foundation Libraries). Přehled novinek v poznámkách k vydání.

Ladislav Hagara | Komentářů: 0
Jak se stavíte k trendu ztenčování přenosných zařízení (smartphony, notebooky)?
 (12%)
 (2%)
 (72%)
 (3%)
 (11%)
Celkem 395 hlasů
 Komentářů: 39, poslední dnes 19:30
Rozcestník
Reklama

Dotaz: Algoritmus na náhodné rozdělení dle gaussovy křivky

Jakub Lucký avatar 26.7.2009 23:06 Jakub Lucký | skóre: 40 | Praha
Algoritmus na náhodné rozdělení dle gaussovy křivky
Přečteno: 2562×
Programuji si takovou netriviální diskrétní simulaci a potřeboval bych poradit s jedním algoritmem.

Součástí mé simulace je odchod lidí z práce na oběd. Ten probíhá v přesně ohraničeném časovému úseku (dejme tomu 11:00 - 14:00) se špičkou ve 12:30... Tím nám vzniká jakási (dejme tomu symetrická) křivka, podle které bych potřeboval generovat náhodné časy odchodu na oběd.

Stručně řečeno: Potřebuji vygenerovat X náhodných čísel v nějakém intervalu, které při poskládání vytvoří gaussovu křivku... Jak na to?

Doufám, že je vysvětlení jasné, už mi z toho jde trochu hlava kolem...

Díky za pomoc
If you understand, things are just as they are; if you do not understand, things are just as they are. (Zen P.) Blogísek

Odpovědi

26.7.2009 23:34 Let_Me_Be | skóre: 20 | blog: cat /proc/idea/current | Brno
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Hledej implementaci generatoru nahodnych cisel, ktery dokaze generovat podle nastaveneho rozlozeni. Pro C++ to umi napriklad Boost.
Linked in profil - Můj web - Nemůžete vyhrát hádku s blbcem. Nejdřív vás stáhne na svoji úroveň a pak ubije zkušenostmi.
Jakub Lucký avatar 27.7.2009 00:34 Jakub Lucký | skóre: 40 | Praha
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
To mě ani nenapadlo, že bych dělal "reinventing the wheel"...

Má někdo tip na nějakou Python knihovnu která tohle umí?
If you understand, things are just as they are; if you do not understand, things are just as they are. (Zen P.) Blogísek
27.7.2009 09:40 Jan Martinek | skóre: 43 | blog: johny | Brno
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky

Řekněme, že chci vygenerovat deset hodnot s normálním rozdělením, střední hodnotou pět a "sigmou" (odmocninou z rozptylu) jedna:

 

>>> from scipy import *
>>> random.normal(5,1,10)
array([ 4.17654928,  5.42693759,  2.65380726,  6.30024096,  5.91105815,
        6.26047915,  3.8240613 ,  5.48209872,  5.05101277,  3.30617717])

A je to :-)

27.7.2009 06:48 pht | skóre: 48 | blog: pht
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Stačí když použijete normální generátor náhodných čísel od 0 do 1 a výsledek pošlete jako parametr do kvantilové funkce. U gausse je tato funkce poměrně složitá, popis algoritmu a implementace lze najít např http://home.online.no/~pjacklam/notes/invnorm/
In Ada the typical infinite loop would normally be terminated by detonation.
27.7.2009 09:30 ivan
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky

Hmm, a nestacilo by vygenerovat n cisel, secist je a pak vydelit n?.

 

27.7.2009 09:49 Filip Jirsák | skóre: 66 | blog: Fa & Bi
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Průměr n náhodných čísel bude mít zase stejný tvar rozdělení, jako původních n čísel, ne?
27.7.2009 10:12 x22
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Nie.
27.7.2009 10:39 Filip Jirsák | skóre: 66 | blog: Fa & Bi
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
No jo, ona náhodná čísla vygenerovaná v počítači asi budou patřit do nějakého omezeného intervalu…
27.7.2009 23:56 Martin | skóre: 10 | blog: Nádraží Perdido
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Důvod pro normální rozdělení průměrů posloupností stejně rozdělených náhodných veličin je v něčem jiném.
27.7.2009 10:28 pht | skóre: 48 | blog: pht
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Pokud si pamatuju dobře tak střední hodnota (="průměr") má normální rozdělení, takže by to asi tak šlo: vezmete např 1000 čísel s uniformním rozdělením od N do M (=rand()), spočtete střední hodnotu, máte jedno náhodné číslo s normálním rozdělením, vezmete dalších 1000, spočtete stř. h., máte další náhodné č., ... postupně vám vyjde gauss s mí=(M+N)/2 a sigma užnevímkolik.

Ale spotřebuje to nepoměrně víc náhodných čísel než kvantilová funkce, přes kterou se to obvykle dělá.
In Ada the typical infinite loop would normally be terminated by detonation.
27.7.2009 10:48 Martin | skóre: 10 | blog: Nádraží Perdido
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
To záleží na tom, co by ta čísla generovalo. Třeba takové rand() % n by bylo pro výše uvedenou simulaci pravděpodobně nepoužitelné. To by asi chtělo použít nějaký propracovanější generátor.

Pokud by nezáleželo až tak na jazyku, udělal bych to v R. Tam je aspoň trochu jistota, že to bude dávat rozumné hodnoty pro různá rozdělení, má to i vlastní programovací jazyk, dělá to grafy... Na simulace ideální věc. Na Pravděpodobnosti a statistice jsme v něm dělali zápočtový projekt, taky diskrétní simulaci.
27.7.2009 10:52 pht | skóre: 48 | blog: pht
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Třeba takové rand() % n by bylo pro výše uvedenou simulaci pravděpodobně nepoužitelné.
Proč?
In Ada the typical infinite loop would normally be terminated by detonation.
27.7.2009 11:10 Martin | skóre: 10 | blog: Nádraží Perdido
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Když jsem si s tím kdysi hrál, tak mi to nedávalo dostatečně rovnoměrně rozdělené hodnoty. Ale možná to tehdy bylo jenom nějaké zkriplené.

Určitě ale existují třeba pro takové C++ mnohem lepší a spolehlivější generátory.
27.7.2009 11:17 Martin | skóre: 10 | blog: Nádraží Perdido
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
(Samozřejmě záleží hlavně na tom, jak moc seriózní tu simulaci člověk potřebuje.)
27.7.2009 12:59 pht | skóre: 48 | blog: pht
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
IMHO stačí číst z /dev/urandom (na Linuxu).
In Ada the typical infinite loop would normally be terminated by detonation.
27.7.2009 23:57 Martin | skóre: 10 | blog: Nádraží Perdido
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Je to možné, takhle jsem to zatím nikdy nezkoušel. Ani vlastně nevím, jak /dev/urandom přesně funguje. Otestuju to a zkusím porovnat.
27.7.2009 12:33 podlesh | skóre: 38 | Praha
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
V dokumentaci výslovně varují před použitím modulo. Když už, tak rand() * N / RAND_MAX

27.7.2009 12:59 pht | skóre: 48 | blog: pht
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
To je ale jen v nějaké referenční implementaci.
In Ada the typical infinite loop would normally be terminated by detonation.
27.7.2009 13:56 podlesh | skóre: 38 | Praha
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Pravda, záleží na algoritmu a "Linux C Library use the same random number generator as random() and srandom(), so the lower-order bits should be as random as the higher-order bits"

27.7.2009 10:08 Vojtěch Horký | skóre: 39 | blog: Vojtův zápisník | Praha
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Tohle by mělo jít bez problémů simulovat v R (a přihodím ještě odkaz s obrázkem).
I am always ready to learn although I do not always like to be taught. (W. Churchill)
27.7.2009 13:09 Radovan
Rozbalit Rozbalit vše Dotaz nematematika

Jak je velký rozdíl v průběhu mezi Gaussovou křivkou a obyčejnou sinusoidou?

hodnota=int(maxvyska*(1-cos(index*pi/limit)))

X
X
XXXX
XXXXXXXXX
XXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXX
XXXXXXXXX
XXXX
X
X

27.7.2009 13:44 petr_p | skóre: 59 | blog: pb
Rozbalit Rozbalit vše Re: Dotaz nematematika
Rozdíl je přesně 1 :D (Za předpokladu, že normální rozložení bereme od −∞ do ∞ a sinus od 0 do π.)
27.7.2009 23:45 Radovan
Rozbalit Rozbalit vše Re: Dotaz nematematika

Mě šlo o tvar té křivky, ne o přesná čísla, interval toho mého grafu je <0;2π> Jen tak od oka bych totiž řekl, že ty průběhy jsou stejné...

28.7.2009 08:30 pht | skóre: 48 | blog: pht
Rozbalit Rozbalit vše Re: Dotaz nematematika
Tak to je opravdu jen od oka ;)
In Ada the typical infinite loop would normally be terminated by detonation.
28.7.2009 09:24 Radovan
Rozbalit Rozbalit vše Re: Dotaz nematematika

Právě proto se ptám :-D

28.7.2009 09:58 Jan Martinek | skóre: 43 | blog: johny | Brno
Rozbalit Rozbalit vše Re: Dotaz nematematika

Tato debata je naprosto mimo původní téma. Ale když máme tu okurkovou sezónu ...

Místo "sinusovky" jsem použil kosinus na druhou (neptej se proč). Nakreslím-li to do stejného grafu spolu s Gaussovkou, dostanu tohle:

kf.fyz.fce.vutbr.cz/pub/gauss_cos2.png

A na co že ses vlastně ptal?

 

28.7.2009 11:37 Radovan
Rozbalit Rozbalit vše Re: Dotaz nematematika

Jo, je to trochu off-topic, takže díky. Ta moje "kosinusovka" je nejen otočená a posunutá, ale i roztažená, takže ve vrcholech se s gaussovkou potkává. Ptal jsem se na tu spodní část, na tom tvém obrázku je ten rozdíl už vidět.

Jakub Lucký avatar 28.7.2009 12:16 Jakub Lucký | skóre: 40 | Praha
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Všem děkuji za návrhy...
If you understand, things are just as they are; if you do not understand, things are just as they are. (Zen P.) Blogísek
28.9.2012 15:10 Wrunx
Rozbalit Rozbalit vše Re: Algoritmus na náhodné rozdělení dle gaussovy křivky
Možná už to je OT, ale: #Python 2.7 gaussian_values = [random.gauss(1, .5) for i in range(1000)]

Založit nové vláknoNahoru

Tiskni Sdílej: Linkuj Jaggni to Vybrali.sme.sk Google Del.icio.us Facebook

ISSN 1214-1267   www.czech-server.cz
© 1999-2015 Nitemedia s. r. o. Všechna práva vyhrazena.