Simulation Methods in Automotive Industry/cs

Přehled metod simulací
Metoda	Klíčová otázka	Zkratka
Metoda Monte Carlo	Jaká jsou rizika a pravděpodobnosti?	MC
Systémová dynamika	Jak se systém chová v čase?	SD
Agentní modelování	Jak vzniká makrochování z mikropravidel?	ABM
Diskrétní událostní simulace	Kolik, jak rychle, kde je úzké místo?	DES
Každá metoda pokrývá jinou vrstvu složitosti systému.

Úvod

V roce 2021 automobilový průmysl zažil krizi, jejíž příčina nebyla ani zdaleka zřejmá. Volkswagen, Toyota, Ford a desítky dalších výrobců byli nuceni zastavit výrobu a celkově nedodali na trh více než 7 milionů vozidel. Příčinou nebyla válka, pandemie ani přírodní katastrofy. Příčinou byl neregulovaný systém zpětných dodavatelských řetězců, který zjevně nebyl seřízen jako jeden mechanismus. Některé společnosti však disponovaly systémově-dynamickými modely svých dodavatelských řetězců, což jim umožnilo reagovat rychleji než ostatním.

Podobnými příklady jsou plné dějiny prakticky každé výrobní oblasti a nejen jí. Uvedený příklad ale zřetelně odráží problematiku automobilového průmyslu a vede k jádru tohoto článku.

V dnešní době automobilový průmysl není ani tak továrnou razící dopravní prostředky, jako spíše složitým systémem, jak technickým tak sociálním, který zahrnuje nepředstavitelné množství faktorů na něj působících. Strategie automobilových společností, dynamika trhů, logistika, chování spotřebitelů v souladu s regionálními specifiky a samozřejmě výrobní proces. To je jen malá část složek automobilového průmyslu, a každá z nich má pod sebou také různé skupiny problematik a detailů. Je zřejmé, že chaotické uspořádání takového systému může, a přesněji řečeno naprosto jistě povede ke kolapsu, a právě proto si ho lze jen těžko představit bez simulací.

V rámci tohoto článku se podíváme na hlavní metody simulací, které se v automobilovém průmyslu dnes používají, určíme jejich oblast působení, uspořádání, praktickou složku a výhody s nevýhodami. Zvláštní pozornost věnujeme metodám pracujícím s nejistotou a složitostí: Metodě Monte Carlo, Systémové dynamice, Agentnímu modelování a Diskrétní událostní simulaci.

Proč automobilový průmysl potřebuje simulace?

Dynamická složitost

Proč vůbec simulace potřebujeme? Odpovědí na tuto otázku je koncept který má název „dynamická složitost". Setkáváme se s ní v systémech kde chování není intuitivní a je těžko předvídatelné, přičemž problém není ve velkém počtu prvků ale v tom že mezi prvky existují nelineární zpětné vazby.

Automobilový průmysl je klasický příklad. Akce v každé fázi systému vede k ovlivnění ostatních. Například si představte že automobilová značka vydá povedený model za příjemnou cenu. To přímo povede k nárůstu prodejů, což zase přímo vede ke zvýšení vytížení továrny. Vytížení továrny vede ke snížení zásob ve skladech dealerů. To zase nutí továrnu navyšovat výrobu, a to vyžaduje zvýšení dodávek a potřebného množství komponentů. Výsledkem je že dodavatelé zdraží, což snižuje marži výrobce. A to je jen jeden scénář z mnoha.

Ve velkých a složitých systémech s propojenými fázemi mají různá manažerská rozhodnutí vždy důsledky, špatné i dobré. Často se stává že se důsledky projeví se zpožděním, a to na těch nejnečekanějších místech.

Simulace jsou určeny k tomu aby ukázaly chování rozhodnutí ve virtuálním světě, místo jejich reálného uskutečnění. To pomáhá zkoumat důsledky bez nutnosti nést riziko a ztráty. Samozřejmě pokud je simulace kvalitně zpracovaná. Ale i dobře detailní systém se všemi procesy a proměnnými se setkává s nejistotou, protože reálný svět je nepředvídatelný.

V metodologii simulací existují metody které nám umožňují zohledňovat i nejistotu, dovolují nám náhodnost nepřehlížet ale pracovat s ní systematicky.

Metoda Monte Carlo: řízení nejistoty

Podstata metody

💡 Zajímavost: Název metody pochází ze slavného kasina v Monaku, symbolu náhody a pravděpodobnosti. Samotná metoda ale nevznikla v nic tak zábavném: Stanislaw Ulam a John von Neumann ji vyvinuli ve 40. letech při práci na jaderných zbraních v rámci projektu Manhattan. Odtud se pak dostala do inženýrství, financí a řady dalších oborů, což u vojenských technologií není nic neobvyklého.

Metoda Monte Carlo patří k těm základnějším věcem ve výpočetním modelování. Princip je jednoduchý: místo jedné konkrétní hodnoty parametru, kterou stejně přesně neznáme, se vezme celý možný rozsah. Simulace pak proběhne mnohokrát, pokaždé s jinak náhodně vybranými hodnotami, a tím pádem dostaneme statistický obraz toho jak se systém chová. Ne jedno číslo, ale rozložení výsledků.

Podrobněji je metoda popsána v článku Metoda Monte Carlo.

Použití v automobilovém průmyslu

Prognózování poptávky a plánování výroby

Výroba se plánuje na základě očekávané poptávky, a poptávka po jednotlivých modelech se zase skládá z mnoha faktorů: nabídka konkurentů, trendy, úrokové sazby v bankách, regionální ceny paliv, infrastruktura (zejména pro elektromobily), terén a tak dále. Výrobce modelu, který plánuje uvést na trh, spouští tisíce iterací simulace, v každé budou klíčové parametry nabývat náhodných hodnot na základě historických rozdělení. Jako výsledek výrobce obdrží rozdělení pravděpodobností. Například s pravděpodobností 0,7 budou prodeje činit od 100 000 do 150 000 jednotek ročně.

Analýza rizik při uvedení nového modelu

Uvedení každého modelu na trh je vždy rizikem miliard investovaných prostředků, rizikem času stráveného vývojem, rizikem image značky a tak dále. I průměrně úspěšný model má horizont návratnosti 5 až 7 let. Metoda Monte Carlo může být využita pro prognózu rozdělení čisté současné hodnoty (NPV) projektu. Metoda zohledňuje nejistotu faktorů jako jsou objemy prodejů, výrobní náklady a tak dále. Výsledkem výrobce získá představu nejen o výnosnosti projektu, ale také o tom s jakou pravděpodobností projekt čeká neúspěch.

Analýza spolehlivosti dodavatelských řetězců

Zásobování komponenty od dodavatelů je nedílnou součástí průmyslu. Výrobní poruchy, zpoždění na celnicích, přetíženost tras, to vše je normální realitou automobilového průmyslu. Monte Carlo může pomoci modelovat možný deficit komponentů při různých kontextech a scénářích. Tímto způsobem automobiloví výrobci určují optimální pojistnou zásobu.

Výhody a omezení

Výhody	Omezení
Jednoduchá implementace, nevyžaduje řešení rovnic	„Odpadky na vstupu, odpadky na výstupu" — přesnost závisí na vstupních rozděleních
Snadno paralelizovatelná na moderních výpočetních systémech	Pro velmi vzácné události (p ≤ 10⁻⁶) nevhodná bez speciálních metod vzorkování
Univerzálně použitelná napříč odvětvími	Nevysvětluje proč se systém chová daným způsobem, pouze jak

Nástroje: Oracle Crystal Ball, @RISK od Palisade

Systémová dynamika: zpětné vazby a dlouhodobá dynamika trhu

Podstata metody

Na počátku 60. let vyvinul Jay Forrester na Massachusettském technologickém institutu metodu modelování systémové dynamiky. Pro tuto metodu platí tři neodmyslitelné koncepty: zásoby, toky a zpětné vazby.

• Zásoby jsou veličiny které se v průběhu času akumulují (například vozidla na skladech u dealerů).
• Toky jsou rychlosti změn množství zásob (například tempo prodeje vozidel).
• Zpětné vazby jsou příčinně-následkové vztahy mezi proměnnými.

Existují dva typy zpětných vazeb. Pozitivní (posilující, R) vede k exponenciálnímu růstu nebo poklesu. Negativní (vyvažující, B) směřuje systém k rovnováze.

Podrobněji je metoda popsána v článku Systémová dynamika.

Použití v automobilovém průmyslu

Modelování trhu elektromobilů

Elektromobily jsou v rámci automobilového průmyslu poměrně novou záležitostí. Jejich rozšíření není zdaleka přímočaré, hodně záleží třeba na tom, jestli má daná oblast vůbec dostatek dobíjecích stanic. Je to hezký příklad systému kde různé zpětné vazby působí současně a ovlivňují se navzájem.

Níže je znázorněn kauzální diagram zpětných vazeb trhu elektromobilů:

╔══════════════════════════════════════════════════════════════════╗
║         KAUZÁLNÍ DIAGRAM: TRH ELEKTROMOBILŮ                      ║
╠══════════════════════════════════════════════════════════════════╣
║                                                                  ║
║  POSILUJÍCÍ SMYČKA (R) — „Rozvoj infrastruktury"                 ║
║                                                                  ║
║  Prodeje EV ──(+)──► Příjmy z dobíjení ──(+)──► Investice        ║
║      ▲                                               │(+)        ║
║      │(+)                                            ▼           ║
║  Atraktivita EV ◄──────────────── Hustota dobíjecí sítě          ║
║                  (+)                                             ║
║                                                                  ║
║  VYVAŽUJÍCÍ SMYČKA (B) — „Zatížení sítě"                         ║
║                                                                  ║
║  Podíl EV ──(+)──► Zatížení el. sítě ──(+)──► Ceny elektřiny     ║
║      ▲                                               │(-)        ║
║      │                                               ▼           ║
║      └──────────────── Atraktivita EV ◄──────────────┘           ║
║                                  (-)                             ║
║                                                                  ║
║  (+) = stejný směr změny    (-) = opačný směr změny              ║
╚══════════════════════════════════════════════════════════════════╝

Efekt biče (Bullwhip effect)

Bullwhip effect patří k základním problémům dodavatelských řetězců. Funguje to takhle: z nějakého důvodu vzrostla poptávka po určitém modelu. Malý skok, ale každý článek řetězce na něj reaguje po svém. Dealer objedná víc vozidel než potřebuje, jen pro jistotu. Výrobce udělá totéž u dodavatele komponentů. Ten udělá totéž u svého dodavatele. Na konci řetězce původně malý výkyv vypadá jako obrovský skok, přestože reálná poptávka se změnila jen minimálně. SD model dokáže tento jev změřit a pomoct najít způsoby jak ho tlumit, třeba lepším řízením zásob nebo rychlejším předáváním informací mezi články řetězce.

Tento jev si lze zažít interaktivně v klasické simulační hře Beer Distribution Game, která byla původně vyvinuta právě na MIT.

Výhody a omezení

Výhody	Omezení
Nezbytná pro pochopení dlouhodobého chování složitých systémů	Pracuje pouze s průměrnými hodnotami — žádné individuální chování agentů
Vizualizace zpětných vazeb pomocí kauzálních diagramů	Méně vhodná pro přesnou optimalizaci výrobních procesů
Odhalí systémové příčiny problémů, nejen symptomy	Kalibrace modelu vyžaduje kvalitní historická data

Nástroje: Vensim, Stella Architect, Powersim

Agentní modelování: spotřebitelé, konkurence a doprava

Podstata metody

Metoda která popisuje systém jako autonomní agenty se nazývá Agent-based modeling. Každý agent interaguje s prostředím a ostatními agenty, má vlastní pravidla chování a stav. Mikro interakce agentů vedou k makrochování celého systému. Jedná se o emergentní výsledek, tedy takový který není předem naprogramován ale formuje se až v průběhu.

Pokud SD odpovídá na otázky „Kolik?" a „Jak rychle?", pak ABM odpovídá na otázky „Proč?" a „Jak to funguje?"

Podrobněji je metoda popsána v článku Agentní modelování.

Modelování automobilového trhu

Automobilový trh je především konkurenčním prostředím kde každý výrobce má vlastní strategii, politiku a intenzitu reakce na kroky ostatních. K tomu přistupují spotřebitelé s různými požadavky a preferencemi. ABM umožňuje toto prostředí modelovat.

Spotřebitelé-agenti jsou charakterizováni preferencemi jako je poměr ceny a kvality, příjmem, sklonem k vyhledávání nových alternativ (switching propensity) a dalšími faktory. Výrobci-agenti jsou naopak charakterizováni reakcemi na tržní podmínky, jako je přizpůsobování vlastností produktu nebo cenová politika.

Simulace silniční dopravy

Pokud uděláme krok stranou od prodeje automobilů, zvláštní pozornosti si zaslouží modelování dopravy — ABM je pro tuto oblast ideální. Každý automobil vystupuje jako agent s vlastními charakteristikami: rychlost, reakce na okolní faktory, vzdálenost k ostatním agentům a tak dále. Z drobných charakteristik agentů se skládá chování systému a vznikají takové fenomény jako dopravní zácpy, efekt harmoniky nebo nehodové úseky. Městské úřady aktivně využívají podobné simulace pro projektování infrastruktury, značení silnic, rozmístění semaforů a přechodů pro chodce.

Funkční ukázky ABM simulací dopravy i dalších systémů jsou dostupné v prostředí NetLogo, kde lze modely spustit přímo v prohlížeči.

Výhody a omezení

Výhody	Omezení
Jediná metoda schopná modelovat emergentní fenomény	Složitější na validaci — těžko prokazatelné že se agent chová správně
Zvládne simulaci i při neúplné znalosti systému	Méně vhodná pro přesnou kvantitativní optimalizaci
Přirozené modelování heterogenního chování (různí agenti, různá pravidla)	Výpočetně náročná při velkém počtu agentů

Nástroje: NetLogo, AnyLogic, Mesa (Python)

Diskrétní událostní simulace: výroba a logistika

Podstata metody

Poslední metoda simulací o které budeme v tomto článku hovořit je Discrete Event Simulation — modeluje systém jako posloupnost diskrétních událostí v čase. DES jako metoda je zvláště efektivní pro modelování výrobních linek, skladů nebo dopravních toků. Diskrétní událostní modely odpovídají na otázky o množství, například „kolik automobilů vyrobí linka za jednu směnu?".

Schéma typické DES výrobní linky:

╔══════════════════════════════════════════════════════════════════════════╗
║           SCHÉMA DES MODELU VÝROBNÍ LINKY                                ║
╠══════════════════════════════════════════════════════════════════════════╣
║                                                                          ║
║  ┌─────────┐   ┌──────────┐   ┌──────────┐   ┌──────────┐   ┌───────┐    ║
║  │ Vstup   │──►│Stanice 1 │──►│ Buffer 1 │──►│Stanice 2 │──►│Výstup │    ║
║  │  dílů   │   │(sváření) │   │          │   │(lakování)│   │       │    ║
║  └─────────┘   └──────────┘   └──────────┘   └──────────┘   └───────┘    ║
║                     │                              │                     ║
║               ┌─────┴──────┐                ┌─────┴──────┐               ║
║               │ BOTTLENECK?│                │ BOTTLENECK?│               ║
║               │Fronta čeká │                │Fronta čeká │               ║
║               └────────────┘                └────────────┘               ║
║                                                                          ║
║  Goldrattovo pravidlo: zvyšování efektivity mimo bottleneck              ║
║  celkovou výkonnost nezvýší.                                             ║
╚══════════════════════════════════════════════════════════════════════════╝

Podrobněji je metoda popsána v článku Diskrétní událostní simulace.

Použití v automobilovém průmyslu

Optimalizace výrobních linek

Nejočividnějším a nejklasičtějším použitím DES v automobilovém průmyslu je optimalizace výrobních linek. Úplný model bude zahrnovat všechny pracovní stanice, nárazníkové zóny, logiku priorit a také časy cyklů a variabilitu. Simulace pomůže odhalit bottlenecks — fáze ve kterých výkonnost naráží na strop a určují celkové maximum. To je kriticky důležité pro společnosti, podle jednoho z nejdůležitějších manažerských pravidel, Goldrattova pravidla: zvyšování efektivity ne v místě bottleneck celkovou výkonnost nezvýší. DES nám tak umožňuje najít páky působení na konkrétní lince.

Plánování nového závodu

Tato metoda je velmi užitečná při plánování otevření nové pobočky — navrhnout a spustit virtuálně nový závod a odhalit problémová místa ještě před tím než budou vynaloženy prostředky. Tak například BMW při plánování závodu ve městě Lipsko v Německu odhalilo kritická problémová místa a odstranilo je ještě ve fázi projektování.

Výhody a omezení

Výhody	Omezení
Nejpřesnější nástroj pro kvantitativní optimalizaci výrobních a logistických systémů	Vyžaduje přesná data, detailní popis struktury a parametry všech prvků
Přímočaré modelování procesů s jasně definovanými kroky	V systémech se silnými zpětnými vazbami méně vhodná než SD
Výborná vizualizace průtoku a čekacích dob	Nemodeluje emergentní nebo individuální chování agentů

Nástroje: SIMUL8, AnyLogic, FlexSim, Arena

Srovnání metod: co a kdy použít

Každá z metod má svou specifiku, ale v tak složitém a mnohotvárném odvětví jako je automobilový průmysl se tyto metody používají kombinovaně. Například výrobce elektromobilu za účelem analýzy strategie uvedení na trh může použít systémovou dynamiku pro pochopení dlouhodobé dynamiky trhu, agentní modelování pro pochopení chování spotřebitelů, Monte Carlo pro představení rizik z finanční stránky a DES pro optimalizaci řetězce výrobní linky.

Metoda	Klíčové otázky	Příklad použití v auto	Typický software
Monte Carlo	Jaká jsou rizika? Jaká je pravděpodobnost?	Hodnocení NPV nového modelu	Crystal Ball, @RISK
Systémová dynamika	Proč se systém chová tak jak se chová? Jaké jsou dlouhodobé trendy?	Dynamika trhu EV	Vensim, Stella
ABM	Jak vzniká makrochování z mikropravidel?	Konkurence na trhu, simulace dopravy	NetLogo, AnyLogic
DES	Kolik? Jak rychle? Kde je úzké místo?	Optimalizace výrobní linky	SIMUL8, Arena

Úlohy pro procvičení a pochopení

Úloha 1 (Monte Carlo). Výrobce automobilu plánuje uvést na trh nový model a potřebuje odhadnout očekávaný objem prodejů. Objem závisí na dvou nejistých faktorech: tempu růstu HDP, které se řídí normálním rozdělením se střední hodnotou 2 % a směrodatnou odchylkou 1 %, a ceně ropy, která se pohybuje rovnoměrně mezi 60 a 100 USD za barel. Pokud růst HDP překročí 2,5 % a cena ropy klesne pod 80 USD, očekává se prodej 120 000 vozů; v ostatních případech 80 000 vozů. Popište postup, jak pomocí simulace Monte Carlo odhadnout střední očekávaný objem prodejů a pravděpodobnost, že prodeje překročí 100 000 kusů.

Řešení: Postup: (1) Vygenerujte N = 10 000 náhodných dvojic (HDP, cena ropy) z daných rozdělení. (2) Pro každou dvojici určete objem prodejů podle podmínky. (3) Průměr všech N výsledků = střední očekávaný objem. (4) Podíl případů s prodejem > 100 000 = hledaná pravděpodobnost. Při daných parametrech (podmínka splněna cca ve 30 % případů) vychází střední objem přibližně 90 000 vozů a P(prodeje > 100 000) ≈ 0,30.

Úloha 2 (Systémová dynamika). Sestavte kauzální diagram pro následující situaci: vysoká poptávka po automobilech tlačí ceny nahoru, což poptávku postupně tlumí; zároveň vysoká poptávka motivuje výrobce k navýšení produkce, avšak s časovým zpožděním přibližně 18 měsíců, což nakonec vede k převisu nabídky. Určete typ každé zpětné vazby ve vašem diagramu.

Řešení: Zpětná vazba 1: Poptávka → (+) → Ceny → (-) → Poptávka. Typ: vyvažující (B). Zpětná vazba 2: Poptávka → (+) → Výroba (se zpožděním 18 měsíců) → (+) → Nabídka → (-) → Ceny → (-) → Poptávka. Typ: vyvažující (B) s časovým zpožděním, které způsobuje oscilace.

Úloha 3 (ABM). V modelu automobilového trhu vybírá agent-zákazník mezi dvěma značkami. Značka A: cena 25 000 €, kvalita 0,8. Značka B: cena 20 000 €, kvalita 0,6. Preference agenta je p = 0,3 (agent preferuje spíše kvalitu než nízkou cenu). Rozhodovací funkce má tvar: f = p × (cena / max_cena) + (1 − p) × (1 − kvalita). Agent vybere značku s nižší hodnotou f. Vypočítejte hodnotu f pro obě značky a určete, kterou si agent vybere.

Řešení: max_cena = 25 000 €. Značka A: f = 0,3 × (25000/25000) + 0,7 × (1 − 0,8) = 0,3 × 1 + 0,7 × 0,2 = 0,30 + 0,14 = 0,44. Značka B: f = 0,3 × (20000/25000) + 0,7 × (1 − 0,6) = 0,3 × 0,8 + 0,7 × 0,4 = 0,24 + 0,28 = 0,52. Agent si vybere Značku A (nižší f = 0,44), přestože je dražší — protože kvalita má pro tohoto agenta větší váhu.

Aktuální trendy a budoucnost

Integrace se strojovým učením

Žijeme v době kdy máme velmi výkonné nástroje pro provádění simulací. S ohledem na rozvoj strojového učení se novou iterací konceptu simulace staly trénované simulace (simulation intelligence). Je to nový směr kde jsou parametry systému kalibrovány na základě dat pomocí metod strojového učení, zejména je to aktuální pro metodu ABM.

Digitální dvojníci

Digital twin je přístup ve kterém paralelně, například s výrobní linkou, vzniká DES model který se v reálném čase aktualizuje pomocí senzorů na zařízení. To umožňuje rychle a přesně diagnostikovat problémy aktuálního stavu systému. Přehled konceptu digitálního dvojníka nabízí například IBM: What is a digital twin.

Autonomní doprava

Jelikož se v dnešní době modely strojového učení zavádějí do všech možných oblastí, dostávají se i do dopravy. V oblasti kde velmi mnoho prvků je regulováno různými mechanismy, jako jsou pravidla silničního provozu, nesení správní nebo trestní odpovědnosti, transparentní politiky priority a tak dále, podobná integrace vyvolává mnoho otázek. Nejzákladnější se kterou se setkávají regulátoři celého světa je, jak bezpečný je autonomní automobil. Zde se ABM simulace dopravy stávají nepostradatelným nástrojem pro výzkum a tvorbu důkazní základny rozhodnutí regulátorů.

Závěr

Automobilový průmysl stejně jako mnoho jiných je velmi úzce spjat se simulacemi. Tento druh spojení nevznikl kvůli trendům, vznikl z nutnosti předvídat chování složitých systémů — automobilový průmysl 21. století k takovým bez jakýchkoli pochyb patří. Je to obrovský, architektonicky složitý, konkurenčně náročný byznys který v sobě zahrnuje procesy absolutně různé povahy a uspořádání. Právě to je důvodem proč mnoho metod simulací systémů zaujalo své zasloužené místo v základech tohoto odvětví. Krize roku 2021, od které se tento článek odvíjel, je připomínkou toho, že složitost systému nevyhnutelně přináší nečekané a často nepříjemné překvapení. V takové realitě jsou simulace nejlepším nám známým způsobem jak je předvídat a být na ně připraveni. Pravděpodobně takovými zůstanou i v budoucnu, stávajíce se stále pokročilejšími a technologičtějšími.

Viz také

• Metoda Monte Carlo
• Systémová dynamika
• Agentní modelování
• Diskrétní událostní simulace
• Řízení dodavatelského řetězce
• Digitální dvojník
• Simulace (přehled)

Užitečné nástroje a zdroje

• AnyLogic — multiparadigmatická simulační platforma (MC, SD, ABM, DES v jednom)
• NetLogo — volně dostupné prostředí pro ABM s rozsáhlou knihovnou modelů
• Vensim — profesionální nástroj pro systémovou dynamiku
• Beer Distribution Game — interaktivní simulace efektu biče
• AnyLogic: Free Simulation Book — volně dostupná učebnice simulací
• Úvod do systémové dynamiky — MIT Sloan School of Management

Reference

• BANKS, Jerry, John S. CARSON, Barry L. NELSON a David M. NICOL. Discrete-Event System Simulation. 5. vyd. Upper Saddle River: Prentice Hall, 2010.
• ELMARAGHY, Hoda A., ed. Changeable and Reconfigurable Manufacturing Systems. London: Springer, 2010.
• FORRESTER, Jay W. Industrial Dynamics. Cambridge: MIT Press, 1961.
• RUBINSTEIN, Reuven Y. a Dirk P. KROESE. Simulation and the Monte Carlo Method. 3. vyd. Hoboken: Wiley, 2016.
• STERMAN, John D. Business Dynamics: Systems Thinking and Modeling for a Complex World. Boston: McGraw-Hill, 2000.
• WILENSKY, Uri a William RAND. An Introduction to Agent-Based Modeling. Cambridge: MIT Press, 2015.

---

Autor: Ulan Kurbanov