• Nie Znaleziono Wyników

Metody Monte Carlo

N/A
N/A
Protected

Academic year: 2021

Share "Metody Monte Carlo"

Copied!
29
0
0

Pełen tekst

(1)

Rachunek Prawdopodobieństwa

i

Statystyka Matematyczna

(2)

Prawdopodobieństwo i Statystyka Matematyczna

dr Jan Malinowski

Generator liczb losowych

(3)

Testy statystyczne dla generatorów liczb losowych

Generator rozkl. rownom.

1. Testy oparte na momentach rozkładu: • test zgodności wartości średniej; • test zgodności wariancji;

• momenty wyższe (asymetria eksces). 2. Test zgodności rozkładu χ2 (test częstości). 3. Test Kołmogorowa (dystrybuant).

4. Test zgodności rozkładu par.

5. Test OPSO (overlapping pairs sparse occupancy)

– dot. analizy nakładających się par liczb uzyskiwanych z generatora. 6. Testy serii.

(4)

Prawdopodobieństwo i Statystyka Matematyczna

dr Jan Malinowski Testy statystyczne dla generatorów liczb losowych

7. Testy zgodności rozkładów statystyk: • wartości ekstremalne i rozstęp;

• testy oparte na statystykach pozycyjnych; • test sum (y = x1 + x2 + ... + xm);

• test d2

• test urodzin dla spacji;

• test najmniejszej odległości w parach. 8. Testy kombinatoryczne:

• test pokerowy; • test kolekcjonera;

• test kolizji i test liczby pustych cel; • test permutacji;

• test oparty na rzędzie losowych macierzy binarnych. 9. Testowanie z pomocą zadań kontrolnych.

(5)
(6)

Prawdopodobieństwo i Statystyka Matematyczna

dr Jan Malinowski

(7)
(8)

Prawdopodobieństwo i Statystyka Matematyczna

dr Jan Malinowski

Test zgodności wartości średniej i wariancji

(9)
(10)

Prawdopodobieństwo i Statystyka Matematyczna

dr Jan Malinowski

Test OPSO (słów)

(11)

Test pokerowy

dla n = 5000 czwórek;

(12)

Prawdopodobieństwo i Statystyka Matematyczna

dr Jan Malinowski

Test pokerowy

Test pokerowy 4 bitowy.

Hipoteza zerowa (H0 Ciąg ma charakter losowy). poziom istotności α

Generujemy długi binarny ciąg losowy. Dzielimy go na bloki 4-bitowe.

n – liczba bloków, n > 80,

k – liczba możliwych wartości w bloku, dla bloków 4-bitowych k = 16 liczb czterobitowych)

(13)

Test pokerowy

Obliczamy

Dalej – jak zwykle, sprawdzamy,

(14)

Prawdopodobieństwo i Statystyka Matematyczna

dr Jan Malinowski

(15)
(16)

Prawdopodobieństwo i Statystyka Matematyczna

dr Jan Malinowski

Metoda odwracania dystrybuanty

Niech F(y) oznacza dystrybuantę

charakteryzującą pewien rozkład zmiennej losowej Y. Ciąg elementów próby yi zdefiniowany przez:

posiada rozkład f(y)

(17)
(18)

Prawdopodobieństwo i Statystyka Matematyczna

dr Jan Malinowski

(19)

Metoda eliminacji Von Neumana

losowanie liczb pseudolosowych

(20)

Prawdopodobieństwo i Statystyka Matematyczna

dr Jan Malinowski 1. Wygeneruj dwie liczby losowe x i y z rozkładu równomiernego,

odpowiednio

x ∈∈∈∈ U( a, b )

oraz

y ∈∈∈∈ U( 0, c )

2. Sprawdź czy y ≤ f( x )

Jeśli tak przyjmij x jako kolejny element generowanej próby losowej; Jeśli nie – powtórz krok 1.

Metoda eliminacji Von Neumana losowanie liczb pseudolosowych

(21)

Metoda eliminacji Von Neumana losowanie liczb pseudolosowych

z zadanego rozkładu

f

(

x

) a b x y f(x) y y y y a a x y b a c c y1 y2 x2 x1 x2 x1 x2 x1 x2 x1 y2 x2 x1 y1 y2 x2 x1 y1 y2 x2 x1 R ∈ U(0, 1) xi = a + R ⋅ (b – a) yi = R ⋅ c x1 – nie x y b a

(22)

Prawdopodobieństwo i Statystyka Matematyczna

(23)

Metoda eliminacji Von Neumana

szacowanie wartości całki oznaczonej

(24)

Prawdopodobieństwo i Statystyka Matematyczna

dr Jan Malinowski

Metoda eliminacji Von Neumana szacowanie wartości całki oznaczonej

a b x y f(x) y y y y a a x y b a c c y1 y2 x2 x1 x2 x1 x2 x1 x2 x1 y2 x2 x1 y1 y2 x2 x1 y1 y2 x2 x1 R ∈ U(0, 1) xi = a + R ⋅ (b – a) yi = R ⋅ c x y b a

=

b

a

dx

x

f

I

(

)

Ntotal = Ntotal + 1 jeśli f( xi ) ≥ yi → Nc = Nc + 1

(

)

(

c

b

a

)

N

N

I

total c

=

(25)

Metoda eliminacji Von Neumana szacowanie wartości całki oznaczonej

a b x y f(x) y y y y a a x y b a c c y1 y2 x2 x1 x2 x1 x2 x1 x2 x1 y2 x2 x1 y1 y2 x2 x1 y1 y2 x2 x1 x y b a

=

b

a

dx

x

f

I

(

)

S

N

N

I

total c

=

S

I

=

1

(

b

a

)

c

S

=

(26)

Prawdopodobieństwo i Statystyka Matematyczna

dr Jan Malinowski

(27)

Metoda eliminacji Von Neumana

szacowanie wartości liczby

π

π

π

π

(28)

Prawdopodobieństwo i Statystyka Matematyczna dr Jan Malinowski 0 1 0 1 x y

[ ]

(

)

4

1

2

/

1

lim

2

π

π

=

=

=

∞ → total o N o

N

N

S

S

Metoda eliminacji Von Neumana szacowanie wartości liczby

π

π

π

π

R ∈ U(0, 1) xi = R1 yi = R2 Ntotal = Ntotal + 1 jeśli p( xi , yi ) ∈ w okręgu → No = No + 1

(29)

Cytaty

Powiązane dokumenty

W roz- dziale 4 pokażemy, że test Dzhaparidze–Nikulina należy do ogólnej klasy testów wynikowych (ang. score tests), co tłumaczy m.in.. dobre własności

The default one-sample t-test task includes sample statistics for the variable death and the hypothesis test results.. Compute a Confidence Interval for

Innym przykładem związanym z analizowaniem i odszumianiem obrazów cy- frowych jest wykorzystanie metod MCMC w obróbce obrazów otrzymanych w tomografii komputerowej SPECT i PET

Za pomocą testu chi-kwadrat niezależności, przyjmując poziom istotności 0,05, zbadaj czy płeć telewidza i rodzaj seriali przez niego oglądanych są niezależne..

Zdziwiło go, że w każdej jest inna liczba ziarenek (mimo że są od jednego producenta!). Zweryfikuj hipotezę, że przeciętna liczba ziarenek ryżu w torebce to 8980... •

(ii) biorąc pod uwagę doskonałe wyniki testów migracji globalnej, jak również ekstremalne warunki testowania (10 dni w temperaturze +40º C ), czas pozostawiania w

ZauwaŜyłem, ze znacznie praktyczniejszym sposobem oceniania prawdo- podobieństwa ułoŜenia pasjansa jest wykładanie kart, czyli eksperymentowanie z tym procesem i po prostu

Kable do przesyłu danych CAT6A: PN-EN 50288-10-1 oraz IEC 61156-5 z późniejszymi zmanami - Przewody elektryczne wielożyłowe stosowane w cyfrowej i analogowej technice przesyłu