210 likes | 345 Views
Ilustracja metody Monte Carlo obliczania całek podwójnych. Autorzy: Bernadetta Brzęczek Lech Groblewicz Tomasz Ożański Piotr Ratajczak Wojciech Hanus. Zasada działania Metody Monte Carlo. (-1,1). (1,1).
E N D
Ilustracja metody Monte Carlo obliczania całek podwójnych Autorzy: Bernadetta Brzęczek Lech Groblewicz Tomasz Ożański Piotr Ratajczak Wojciech Hanus
Zasada działania Metody Monte Carlo (-1,1) (1,1) Spróbujmy za pomocą powyższej metody obliczyć pole koła o R=1 i środku w punkcie P=(0,0) x2+y2=1 1. Na kole opisujemy kwadrat o wierzchołkach w punktach: 2. Losujemy n punktów z powierzchni opisanego kwadratu 3. Sprawdzamy czy wylosowane punkty należą do pola koła. (0,0) (-1,-1) (1,-1)
2 1 4 3 … Wynikiem losowania jest informacja, że z n wszystkich prób k było trafionych, zatem pole koła wynosi: gdzie P jest polem kwadratu opisanego na kole Przykład ten ilustruje ogólny sposób działania metody Monte Carlo k-1 n k
Przedstawienie doświadczenia Do prezentacji metody Monte Carlo obliczającej całki podwójne wykorzystaliśmy następujące całki: 1. 2. 3.
Dla każdej z nich wygenerowaliśmy metodą Monte Carlo po 100 wartości dla n z przedziału od 1 do 50 000 000. … … Przykładowa skrócona tabela wartości losowań dla trzech pierwszych i dwóch ostatnich prób w zależności od wielkości n dla całki nr 2.
3. Następnie przy użyciu kalkulatora obliczyliśmy prawidłowe wartości powyższych całek. Dla lepszego zrozumienia z jakimi całkami mamy do czynienia zamieszczone są poniżej także wykresy funkcji podcałkowych. 1. = 2669,851 2. 3. = 1,145521 1. 2.
Z każdej serii danych obliczyliśmy średnią, odchylenie standardowe oraz wartość błędu, czyli różnicę pomiędzy wartością prawidłową a średnią. Powyższe dane przedstawiliśmy za pomocą wykresów.
Poglądowe przedstawienie wzrostu dokładności obliczeń metodą Monte Carlo wraz ze wzrostem n dla całki nr 1
Poglądowe przedstawienie wzrostu dokładności obliczeń metodą Monte Carlo wraz ze wzrostem n dla całki nr 2
Poglądowe przedstawienie wzrostu dokładności obliczeń metodą Monte Carlo wraz ze wzrostem n dla całki nr 3
Dla wszystkich całek można zauważyć, że wraz ze wzrostem n (liczby losowanych punktów), maleje wartość błędu co potwierdzają także wykresy przedstawiające odchylenia standardowego. Na następnych wykresach przedstawione jest jak w zależności od n, na różnie wyskalowanych wykresach, dla poszczególnych całek zmienia się wartość obliczonej całki.
Wnioski: Obliczenia metodą Monte Carlo dają bardzo dokładne wyniki tylko wtedy, gdy rząd parametru n jest większy od 50 000. Dla takich prób dokładność pomiaru jest wyższa niż 99,93%