-
Data: 2016-10-11 21:47:33
Temat: Re: kontynuacja generatory: mersen vs ranlux
Od: bartekltg <b...@g...com> szukaj wiadomości tego autora
[ pokaż wszystkie nagłówki ]On 11.10.2016 17:44, Borneq wrote:
> W dniu 09.10.2016 o 15:26, M.M. pisze:
>> Pierwszy wniosek: Do tej pory nigdy nie zaobserwowałem aby MT
>> oblał stabilny test. Odpaliłem kilka testów dla ranluxa,
>
> Test urodzinowy to przybliżenie:
> po pierwsze jest rozkład Poissona a ang. wiki mówi że ma być wykładniczy
> po drugie - wartość oczekiwana to m^3/4n a gdy m=n to wychodzi n^2/4
> podczas gdy maksymalnie może być n. Z tego wynika że ten wzór to
> przybliżenie.
Z tego wynika, zę mieszasz dwa testy o podobnej tej samej nazwie!
Jeden mówi*) o rozkładzie odstępów (i to w przybliżeniu ciagłym,
a w rzeczywistośći losujemy z dyskretnego zbioru!)
Drugi o liczbie _powtórzeń_ na liście odstępów.
m^3 /(4n) dotyczy tego drugiego (też jest tylko asymptotyczne).
*) do tego nieźle oszukuje;> różnice wspołrzednych będą
liczbami z wykładniczego o gęstości ? jesli wygenerujemy
na tym przedziale proces Poissona. Np losując liczbę
N z rozkładu Poissona (? * dlugość przedzialu), po czym
losując jednostajnie N liczb z tego przedziału.
Wtedy różnice pomiędzy kolejnymi (posortowanymi) liczbami
są niezależnymi zmiennymi z rozkłądu wykładniczego.
Jeśli mamy _zadane_ N i losujemy N liczb jednostajnie na odcinku,
ich różnice:
nie są niezależne!
różnice x1-0, x2-x1, x3-x2.. 1-x_n
mają rozkład Dirichleta Dir(1,1,1,1,1,1).
Zapewne dąży do tego, co trzeba, ale też jest
to rozkłąd przybliżony, i to dość nieźle**, ale nadal
tylko asymptotyczny, tak jak w przypadku testu na
powtórzenia.
**)
n=100000000;
r=rand(n,1); %jednorodny na [0,1]
k=1000;
r=sort(r);
[h,xx]=hist(diff(r),k);
r=[];
semilogy(xx,h/max(xx)/n*k,'.',xx, n*exp(-(n)*xx),'-')
https://www.dropbox.com/s/jlltp0p2chh18bo/b_test.eps
?dl=0
Układa się doćś ładnie na krzywej odpowiedniego exp,
ale oczywisćie nie idealnie.
Losowe rzeczy, które się otworzyły z googla
https://www.jstatsoft.org/article/view/v007i03/tufte
sts.pdf
https://www.cs.indiana.edu/~kapadia/project2/node21.
html
https://en.wikipedia.org/wiki/Dirichlet_distribution
pzdr
bartekltg
Następne wpisy z tego wątku
Najnowsze wątki z tej grupy
- Nowa ustawa o ochronie praw autorskich - opis problemu i szkic ustawy
- Alg. kompresji LZW
- Popr. 14. Nauka i Praca Programisty C++ w III Rzeczy (pospolitej)
- Arch. Prog. Nieuprzywilejowanych w pełnej wer. na nowej s. WWW energokod.pl
- 7. Raport Totaliztyczny: Sprawa Qt Group wer. 424
- TCL - problem z escape ostatniego \ w nawiasach {}
- Nauka i Praca Programisty C++ w III Rzeczy (pospolitej)
- testy-wyd-sort - Podsumowanie
- Tworzenie Programów Nieuprzywilejowanych Opartych Na Wtyczkach
- Do czego nadaje się QDockWidget z bibl. Qt?
- Bibl. Qt jest sztucznie ograniczona - jest nieprzydatna do celów komercyjnych
- Co sciaga kretynow
- AEiC 2024 - Ada-Europe conference - Deadlines Approaching
- Jakie są dobre zasady programowania programów opartych na wtyczkach?
- sprawdzanie słów kluczowych dot. zła
Najnowsze wątki
- 2025-03-19 Warszawa => JavaScript / Node / Fullstack Developer <=
- 2025-03-19 Gdańsk => PHP Developer <=
- 2025-03-19 Kraków => MS Dynamics 365BC/NAV Developer <=
- 2025-03-19 Aresztowany na rok "powinien podziękować za to, że miał możliwość przebywania w zakładzie karnym, bo tam jego stan zdrowia się poprawił"
- 2025-03-19 Chrzanów => Specjalista ds. public relations <=
- 2025-03-19 China-Kraków => Key Account Manager IT <=
- 2025-03-19 Wrocław => Konsultant wdrożeniowy Comarch XL (Logistyka, WMS, Produk
- 2025-03-19 Prezydent Duda śmie szkalować bodnaturę niepokalanie poczętą, dwóch pełnomocników Geralda B. i "standardy Tuskie"
- 2025-03-18 Tesla na złom
- 2025-03-18 Ziobrotura 3.0 będzie w prawie przesłuchać "świadka" Tuska bez adwokata w sprawach Sienkiewicza/Bodnara/...?
- 2025-03-18 Produkty ,,Made in Germany" wciąż na topie - art. na www.dw.com
- 2025-03-18 ulaskawienia
- 2025-03-18 Gdynia => Sales Executive / KAM <=
- 2025-03-18 42 MILIARDY ZŁOTYCH ZYSKU W ROK. DLACZEGO BANKI TYLE ZARABIAJĄ W POLSCE?
- 2025-03-17 Nie matura lecz chęć szczera ...