Metoda Monte Carlo

Monte Carlo je třída algoritmů pro simulaci systémů. Jde o stochastické metody používající pseudonáhodná čísla. Typicky využívány pro výpočet integrálů, zejména vícerozměrných, kde běžné metody nejsou efektivní. Metoda Monte Carlo má široké využití od simulací experimentů přes počítání určitých integrálů až třeba po řešení diferenciálních rovnic. Základní myšlenka této metody je velice jednoduchá, chceme určit střední hodnotu veličiny, která je výsledkem náhodného děje. Vytvoří se počítačový model toho děje a po proběhnutí dostatečného množství simulací se mohou data zpracovat klasickými statistickými metodami, třeba určit průměr a směrodatnou odchylku.

Historie

Metoda Monte Carlo byla formulována již ve 40. letech 20. +more století a svého využití se dočkala ještě v průběhu druhé světové války. Jejím zakladatelem byl Stanisław Marcin Ulam a John von Neumann, kteří v té době pracovali v americké Národní laboratoři Los Alamos, kde zkoumali chování neutronů, především je zajímalo, jaké množství neutronů projde různými materiály (např. nádrží vody).

Přes velké množství informací nebylo možné tento problém vyřešit teoreticky ani prakticky. K výsledku dopomohla až metoda Monte Carlo, kdy se autoři nechali inspirovat kolem rulety (odtud také název Monte Carlo). +more Bylo jim známo, že k pohlcení neutronu jiným atomem dojde v přibližně jednom případu ze sta. Každé roztočení rulety by simulovalo pohyb neutronu, pokud by se zastavila na dílku, který znázorňuje pohlcení neutronu, neutron by cestu neprošel. To by se opakovalo vždy tak dlouho, dokud by neutron nebyl pohlcen nebo dokud by úspěšně prošel celou cestu.

Ač je tento příklad velmi zjednodušený, podstatu metody Monte Carlo vystihuje. Simulace by byla pochopitelně velmi časově náročná, pokud by se skutečně pokaždé točilo ruletou. +more V té době však již vznikaly první počítače, které tyto simulace výrazně zkrátily.

Buffonova jehla a výpočet hodnoty π

Vzdáleným předchůdcem metody Monte Carlo je tzv. Buffonova jehla. +more Je to úloha z roku 1777, jejímž autorem je francouzský matematik Georges Louis Leclerc de Buffon. Jde o pokus, kdy se na arch papíru, který je rozdělen rovnoběžnými linkami, hází jehlou, která má stejnou délku, jako je vzdálenost mezi čárami. Lze spočítat, že pravděpodobnost, že jehla některou čáru protne, je 2/π. Z toho se dá na základě několika opakovaných pokusů odhadnout hodnota π.

Také dnes se princip metody Monte Carlo často demonstruje na následujícím odhadu hodnoty Ludolfova čísla π. (V praxi se ovšem metoda Monte Carlo k tomuto účelu nepoužívá, protože výpočet Ludolfova čísla jinými metodami je mnohem efektivnější. +more) Základem je čtverec, kterému je vepsaná čtvrtkružnice. Analogicky jako u Buffonovy jehly je možné jeho hodnotu zjistit náhodným házením drobnými předměty do prostoru čtverce, a výsledný poměr hodů do čtvrtkruhu ku počtu všech hodů dá odhad čtvrtiny čísla π. Rychlejší je ovšem napsat počítačový program, který házení předmětů simuluje pomocí generátoru pseudonáhodných čísel.

výpočet pí

obsah čtvrtkružnice: S_1 = \frac{\pi r^2}{4}

obsah čtverce: S_2 = r^2

jejich poměr je pak :\frac{S_1}{S_2} = \frac{\pi r^2}{4 r^2} = \frac{\pi}{4}

pak \pi = 4 \frac{S_1}{S_2}

Program v jazyce Python 3, který tento výpočet realizuje a je obohacený o výpočet standardní chyby, skutečné chyby a intervalu spolehlivosti odhadu, může vypadat takto:

import random, math def mc_pi(n=100): """ Odhad čísla pí na základě n pokusů metodou Monte Carlo. """ print ("Odhad Ludolfova čísla na základě", n, "pokusů:\n") v_ctvrtkruhu = 0 # čítač bodů, které padly do čtvrtkruhu for i in range(n): # opakování experimentu n-krát x, y = random. +morerandom, random. random # generuje náhodný bod if x * x + y * y.

Aplikace

Přesnost a efektivnost celého výpočtu metodou Monte Carlo pomocí výpočetní techniky je dána těmito faktory: * kvalitou generátoru náhodných čísel, resp. pseudonáhodných čísel * výběrem racionálního algoritmu výpočtu * kontrolou přesnosti získaného výsledku.

Metoda MC zahrnuje: * vytvoření modelu skutečného systému, se stejnými pravděpodobnostními charakteristikami jako má reálný systém (vliv náhody - náhodná čísla) * model musí zahrnovat veškeré relevantní skutečnosti, podstatně ovlivňující reálný systém * experimentování s modelem, mnohanásobné zkoumání chování modelu ** s pevným časovým krokem - sledujeme chování systému po určitých konstantních časových intervalech a zjišťujeme, zda došlo ke změnám) ** s proměnným časovým krokem - generujeme interval, po který v systému nedojde k žádným změnám

schéma monte carlo Schéma postupu metody Monte Carlo

Řešení problému metodou Monte Carlo můžeme rozdělit do tří kroků: * Rozbor problému a návrh modelu - z hlediska řešení problému se jedná o nejdůležitější krok. I když je MMC použitelná prakticky u všech problémů a její formulace není složitá, nalezení vhodného postupu může nezkušenému řešiteli dělat problémy. +more * Generování náhodných veličin, jejich transformace na veličiny s daným pravděpodobnostním rozdělením. Rychlost konvergence chyby výsledku k nulové hodnotě je u MMC rovna přibližně převrácené hodnotě odmocniny z počtu realizovaných pokusů N, z čehož plyne, že nepatří mezi metody nejefektivnější. * Statistické zpracování výsledků - hledaná hodnota je zpravidla dána některým z momentů statistických veličin, nejčastěji střední hodnotou.

Oblasti použití

Metoda Monte Carlo má širokou možnost využití. Obecně se dá říci, že je možné ji použít všude tam, kde je řešení možné nalézt pomocí mnohokrát opakovaných náhodných pokusů. +more Jelikož se jedná o metodu stochastickou, je nutné znát pravděpodobnostní rozdělení sledovaných veličin. Tyto problémy lze nalézt ve všech oborech, nejen v matematice, ale také v oblasti financí a obchodu, fyzice a fyzikální chemii, ve výpočetní technice a hrách apod.

Matematika

Metoda Monte Carlo je použitelná nejen k řešení jednoduchých určitých i vícerozměrných integrálů, parciálních diferenciálních rovnic, ale třeba také řešení systémů lineárních rovnic, či hledání kořenů rovnic.

Numerická matematika

Příkladem použití metody Monte Carlo je Rabinův-Millerův test prvočíselnosti. Používá se k rychlému rozhodnutí, zda je dané celé kladné číslo N prvočíslem nebo zda je složené. +more Počítá několik mocnin čísla b modulo N, první exponent je číslo N-1, každý další exponent je polovina předcházejícího, pokud je předcházející exponent sudý. Test končí ve chvíli, kdy je odpovídající mocnina různá od jedné i od minus jedné modulo N a číslo je správně klasifikováno jako složené, nebo v okamžiku, kdy buď už exponent nelze dělit dvěma, nebo se mocnina rovná minus jedné modulo N a tehdy se číslo označuje jako silné pseudoprvočíslo při bázi b. Pravděpodobnost, že silné pseudoprvočíslo není prvočíslem, je menší, než jedna osmina. Opakujeme-li tedy test s jinými volbami báze, a vždy dostáváme odpověď, že se jedná o silné pseudoprvočíslo, pravděpodobnost, že se skutečně jedná o prvočíslo, prudce vzrůstá.

Rozlišují se dvě varianty metody Monte Carlo, analogový a neanalogový model.

Analogový model

Musíme umět modelovat celou situaci na počítači, to znamená například znát všechna pravděpodobnostní rozložení zkoumaných jevů a fyzikální zákonitosti, kterými se řídí. Provedením této simulace získáme výsledek, realizaci jakési náhodné veličiny \xi. +more Tuto simulaci spustíme n-krát a získáme soubor historií x1 . xn . Odhad střední hodnoty \xi se určí :.

:\xi \approx \frac{1}{N}\sum_{i=1}^N x_i

a směrodatná odchylka se určí jednoduše jako

:\sigma = \sqrt{ \frac{\sum_{i=1}^n (\bar{x} - x_i)^2}{n-1}}

Neanalogový model

Tak se nazývá případ, kdy při výpočtu nepoužíváme model reálného děje. Například výpočet určitého integrálu, případně obsahu ohraničeného útvaru.

Matematická statistika

Metoda Monte Carlo se v matematické statistice používá například k počítání přesnosti odhadů (resampling, bootstrapping) nebo pro výpočty v bayesovské statistice (MCMC a jiné).

Přesnost metody Monte Carlo

K odhadu chyby výsledku získaného metodou Monte Carlo se většinou používá střední kvadratická chyba aritmetického průměru. Chyba výsledku získaného pomocí n historií je úměrná 1/\sqrt{n}. +more Takže aby se zlepšil výsledek o jeden řád musí se počet historií zvýšit alespoň o dva řády. Abychom získali výsledek s přesností na 6 desetinných míst, což odpovídá přesnosti jiných metod musíme získat 1012 historií.

Další použití metody

* Buffonova úloha * projekt Manhattan * Millerův-Rabinův test prvočíselnosti

Fyzika

Kromě již zmíněné fyzikální chemii se metoda Monte Carlo používá ke složitým výpočtům v kvantové chromodynamice, aerodynamice, dále ve fyzice polymerů, statistické fyzice, fyzice částic. Monte Carlo najde uplatnění i při předpovědi počasí.

Počítačová grafika

Zejména metoda Quasi Monte Carlo se využívá při renderování ve 3D modelech. Grafické programy mají funkci, která vytváří odrazy světla od různých ploch, Monte Carlo je generuje náhodně v rámci polokulové distribuce, Quasi Monte Carlo je vytváří pravidelně. +more Využití - počítačové hry, grafika, animace, filmové efekty, architektura apod. Jedna z metod renderování je radiozita, vychází ze zákona zachování energie a definuje se jako:.

B_i = E_i + R_i \int_j B_j F_ij

kde: Bi je radiozita plošky i. Ei je vyzařovaná energie této plošky. +more Ri je odrazivost plošky. integrál reprezentuje součet energií přicházejících na plošku i ze všech ostatních plošek. Fij je konfigurační faktor mezi ploškami i a j (vliv plošky j na plošku i).

Hazardní hry

Hazardní hry jsou založeny na náhodě, tak jako u jiných her je popsána systémem pravidel a je posloupností tahů. Právě výsledky tahů jsou u hazardních her zcela náhodné - náhodné tahy, u jiných her, kde jsou tahy ovlivněny vůlí hráče, je označujeme jako osobní tahy. +more Monte Carlo je ukázkou procesů, kdy systém osobních tahů bude simulovat systém náhodných tahů. Právě hazardní hry byly východiskem pro vznik a rozvoj teorie pravděpodobnosti. Monte Carlo resp. generátory náhodných a pseudonáhodných čísel umožňují mechanizovat hazardní hry, protože umožňují nastavit všem alternativám stejnou pravděpodobnost.

Finance a pojišťovnictví

V oblasti ekonomiky je možné metodu Monte Carlo využít pro oceňování opcí, investic a jiných finančních derivátů, analýzu rizika, pro zjištění optimální hodnoty portfolia, pro rozhodování či zajišťování.

Analýza rizika

Metoda Monte Carlo dává nejpřesnější pravděpodobnosti ve srovnání s jinými metodami. I přes značné výhody naráží na řadu překážek: * vysoká citlivost výsledků metody Monte Carlo k zákonům pravděpodobnostního rozdělení a typu závislostí vstupních proměnných; * i když současné programy umožňují vzít v úvahu zákony rozdělení pravděpodobnosti, provést korelaci mezi vstupními proměnnými a zhodnotit jejich spolehlivost, v praxi to obvykle není možné, protože ve většině případů analytici stanoví změny základních proměnných makro a mikro prostředí, vybírají zákony rozdělení pravděpodobnosti a statistické vztahy mezi proměnnými subjektivně. +more * V důsledku výše uvedených důvodů je přesnost výsledných odhadů do značné míry závislá na kvalitě základních předpokladů a vzájemných souvislostech vstupních proměnných, což může vést k významným chybám ve výsledcích. Z praktického hlediska je použitelnost metody Monte Carlo kvůli velkému množství zjednodušených předpokladů modelů značně omezená.

Pravděpodobnostní posudek spolehlivosti

Metoda Monte Carlo může být použita také pro pravděpodobnostní posudek spolehlivosti technických struktur (strojů a jejich částí, stavební konstrukce atd. ) Pravděpodobnostní posudek spolehlivosti je pak alternativou ke klasickému staršímu deterministickému posudku spolehlivosti. +more Jedním z průkopníků pravděpodobnostního posudku spolehlivosti technických struktur byl český inženýr Pavel Marek.

Další možnosti využití

Využití metody Monte Carlo je opravdu velmi široké, musí pouze splňovat výše uvedené předpoklady. Kromě již zmíněného se Monte Carlo používá také: * ve výzkumu polovodičů k modelování přenosu nabitých částic * ve studiu životního prostředí, v oblasti znečištění * tepelné štíty * při pátracích a záchranných činnostech - modely používané k předpovědím pohybu záchranných člunů nebo rozšiřování ropné skvrny na moři * předpovědi struktury bílkovin

Optimalizace

Optimalizačních problémů je celá řada, ve všech oborech, např. ve fyzice hledání optimální tloušťky tlakové nádoby, či v umělé inteligenci hledání optimální trajektorie robota apod. +more Optimalizační problémy se dají řešit buď analyticky, nebo pokud to možné není, je možné je řešit optimalizačními algoritmy (jako např. uvedené příklady). Optimalizační algoritmy převedou daný problém na matematický, jehož optimalizace vede k nalezení argumentů tzv. účelové funkce, což je cílem optimalizace. Optimalizační algoritmy hledají minimum dané účelové funkce.

Optimalizační metody je možné rozdělit do tří skupin: * deterministické - je třeba znát předběžné předpoklady, dávají pouze jedno jediné řešení ** horolezecký algoritmus ** metoda větví a mezí - klidně dá všechna řešení ** simplexová metoda ** lineární programování ** matematická analýza - . * stochastické - založeny na náhodě ** simulované žíhání ** Monte Carlo ** slepé prohledávání ** stochastická optimalizace ** stochastické tunelování ** greedy neboli hladový algoritmus - Drsná ignorance determinismu. +more Náhoda a spol. se někdy používají pro výběr z několika nejlepších kroků * smíšené ** genetické algoritmy ** matematické programování ** evoluční strategie.

Optimalizovat se dá i pomocí metastrategií.

Generátory náhodných čísel

Generátory náhodných čísel s definovaným stochastickým rozdělením jsou základem simulačních programů Monte Carlo. Generátory náhodných čísel fungují tak, že nejprve vygenerují posloupnost náhodných čísel s rovnoměrným rozdělením (primární generátor) a poté je z nich transformací vytvořena posloupnost s požadovaným rozdělením. +more Primární generátory jsou fyzikální nebo pseudonáhodné.

Fyzikální generátory náhodných čísel

Jako fyzikální generátor náhodných čísel se dá chápat vše, co má charakter náhodnosti, jako např. házení kostkou nebo mincí. +more Častěji se ale používají generátory založené na kombinaci radioaktivního zářiče a detektoru (vyzařuje částice v náhodném množství). Fyzikální generátory se ale málo používají, jelikož mohou být v závislosti na vnějších vlivech nestabilní, posloupnost čísel se nedá zopakovat a chování se může lišit v důsledku výrobních tolerancí.

Generátory pseudonáhodných čísel

Nejpoužívanějším{{Doplňte zdroj|aktualizovat :-}, taky by se hodil nějaký úvod}} generátorem pseudonáhodných čísel je tzv. kongruentní generátor. +more Posloupnost lze vyjádřit vztahem:.

X_{n+1}=(aX_n+c) mod\ m

kde mod m je celočíselný zbytek po dělení a, c, m jsou zvolené konstanty

Generovaná čísla jsou diskrétní z intervalu

Odkazy

Reference

Literatura

Hušek R. Ekonometrická analýza. Ekopress, Praha 1999.

Související články

Genetický algoritmus - další skupina numerických metod používající sekvence náhodných čísel.

Externí odkazy

Kategorie:Numerická matematika Kategorie:Teorie pravděpodobnosti Kategorie:Matematická statistika Kategorie:Finanční rizika Kategorie:Hazardní hry Kategorie:Algoritmy Kategorie:Počítačová grafika Kategorie:Optimalizace (matematika) Kategorie:Matematické problémy