ˇ ´ vysoke ´ uc ˇen´ı technicke ´ v Praze Cesk e ´ Fakulta elektrotechnicka
´ PRACE ´ DIPLOMOVA Agentn´ı model pˇ repravn´ı kontroly ve veˇ rejn´ e dopravˇ e
Praha, 2012
Autor: Bc. Jan S´ alus
Podˇ ekov´ an´ı Na tomto m´ıstˇe bych chtˇel podˇekovat vˇsem, kteˇr´ı mi umoˇznili dokonˇcen´ı t´eto pr´ace. Zvl´aˇstˇe vedouc´ımu pr´ace Ing. Michalu Jakobovi Ph.D. za podnˇetn´e rady a v´ ybornou spolupr´aci.
i
ii
Abstrakt Tato diplomov´a pr´ace je zamˇeˇrena na pouˇzit´ı netradiˇcn´ı simulaˇcn´ı techniky, kterou je agentn´ı modelov´an´ı. Navazuje na v´ yvoj softwarov´ ych prostˇredk˚ u pro simulaci dopravn´ıch syst´em˚ u AgentPolis a d´ale rozˇsiˇruje modul FareEvasion, kter´ y simuluje syst´em veˇrejn´e dopravy. Je zde vyhodnocena sada ilustraˇcn´ıch simulac´ı prov´adˇen´ ych na veˇrejn´e dopravn´ı s´ıti Los Angeles. Testy jsou nejprve zamˇeˇreny na jednu linku a n´aslednˇe sleduj´ı dynamiku v cel´e s´ıti, ovˇeˇruj´ı adaptaci pasaˇz´era na prob´ıhaj´ıc´ı kontroly. V pr´aci je navrˇzen pl´anovac´ı algoritmus vyuˇz´ıvaj´ıc´ı cenovou funkci a zkuˇsenosti cestuj´ıc´ıho k hled´an´ı v´ yhodn´e cesty prostˇredky veˇrejn´e dopravy, vˇcetnˇe zohlednˇen´ı faktor˚ u ovlivˇ nuj´ıc´ıch jejich rozhodnut´ı. Implementuje model cestuj´ıc´ıho jako adaptivn´ıho agenta interaguj´ıc´ıho s prostˇred´ım, coˇz je vyuˇziteln´e pro eliminov´an´ı u ´niku zisk˚ u dopravn´ıch spoleˇcnost´ı.
iii
Abstract This diploma thesis is focused on using innovative simulation technique known as agent-based modeling. It follows the development of a platform and tools for agent-based transport systems simulation AgentPolis and extends FareEvasion module that simulates public transportation networks. There is a set of illustrative simulations evaluated out to the public transport network in the City of Los Angeles. Tests are first focused on one transportation line and then follow the dynamics across the entire network verify the passenger’s adaptation to the ongoing control. The work proposed planning algorithm using a cost function and the passenger experience convenient way to search for public transport, including consideration of factors influencing their decision. Implements the model as a passenger adaptive agent interacting with the environment, which is useful for eliminating the leakage of profits of transportation companies.
iv
v
vi
Obsah Seznam obr´ azk˚ u
xi
Seznam tabulek
xiii
´ 1 Uvod
1
2 Agentn´ı modelov´ an´ı
3
2.1
Co je agentn´ı modelov´an´ı . . . . . . . . . . . . . . . . . . . . . . . . . . .
3
2.2
Srovn´an´ı pˇr´ıstup˚ u . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4
2.3
V´ yhody agentn´ıho modelov´an´ı . . . . . . . . . . . . . . . . . . . . . . . .
4
2.4
Pouˇzit´ı pˇr´ıstupu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5
2.5
Shrnut´ı . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5
3 Souvisej´ıc´ı pr´ ace
7
3.1
Anal´ yza prostˇred´ı cestuj´ıc´ıch . . . . . . . . . . . . . . . . . . . . . . . . .
7
3.2
Modelov´an´ı rozhodov´an´ı cestuj´ıc´ıch . . . . . . . . . . . . . . . . . . . . .
8
3.3
Shrnut´ı . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
8
4 Specifikace a formalizace probl´ emu 4.1
Reprezentace prostˇred´ı pasaˇz´era . . . . . . . . . . . . . . . . . . . . . . .
9 9
4.1.1
Stavy v prostˇred´ı pasaˇz´era . . . . . . . . . . . . . . . . . . . . . .
11
4.1.2
Akce v prostˇred´ı pasaˇz´era . . . . . . . . . . . . . . . . . . . . . .
11
4.1.3
Mnoˇzina akc´ı j´ızdy metrem . . . . . . . . . . . . . . . . . . . . . .
12
4.1.4
Mnoˇzina akc´ı pˇrechodu mezi stanicemi . . . . . . . . . . . . . . .
12
4.1.5
Mnoˇzina akc´ı vstup do stanice . . . . . . . . . . . . . . . . . . . .
13
4.1.6
Mnoˇzina akc´ı v´ ystup ze stanice . . . . . . . . . . . . . . . . . . .
13
4.2
Vlastnosti pasaˇz´era . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
13
4.3
Pl´an . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
14
vii
4.4
4.5
Model ocenˇen´ı cestovn´ıho pl´anu . . . . . . . . . . . . . . . . . . . . . . . ˇ 4.4.1 Casov´ a sloˇzka ceny . . . . . . . . . . . . . . . . . . . . . . . . . .
15
4.4.2
Ekonomick´a sloˇzka ceny . . . . . . . . . . . . . . . . . . . . . . .
16
4.4.3
Psychologick´a sloˇzka ceny . . . . . . . . . . . . . . . . . . . . . .
18
Shrnut´ı . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
18
5 Uvaˇ zov´ an´ı pasaˇ z´ era 5.1
15
19
Pl´anov´an´ı trasy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
19
5.1.1
Metoda vˇetv´ı a mez´ı . . . . . . . . . . . . . . . . . . . . . . . . .
20
5.1.2
Popis pl´anov´an´ı . . . . . . . . . . . . . . . . . . . . . . . . . . . .
22
5.1.3
Korektnost algoritmu . . . . . . . . . . . . . . . . . . . . . . . . . ´ Uplnost algoritmu . . . . . . . . . . . . . . . . . . . . . . . . . . .
23
5.2
Uˇcen´ı se o pˇrepravn´ı kontrole . . . . . . . . . . . . . . . . . . . . . . . .
23
5.3
Shrnut´ı . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
24
5.1.4
6 Implementace 6.1
6.2
25
Tˇr´ıdy pro model pasaˇz´era . . . . . . . . . . . . . . . . . . . . . . . . . .
25
6.1.1
Vytvoˇren´ı pasaˇz´era . . . . . . . . . . . . . . . . . . . . . . . . . .
25
6.1.2
Vlastnosti pasaˇz´era . . . . . . . . . . . . . . . . . . . . . . . . . .
26
6.1.3
27
6.1.4
Pl´an pasaˇz´era . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ˇ Zivotn´ ı cyklus pasaˇz´era . . . . . . . . . . . . . . . . . . . . . . . .
6.1.5
Pl´anovaˇc . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
28
6.1.6
Stav . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
28
6.1.7
Ukl´ad´an´ı informac´ı . . . . . . . . . . . . . . . . . . . . . . . . . .
29
Tˇr´ıdy pro model inspektora . . . . . . . . . . . . . . . . . . . . . . . . .
29
6.2.1
Vytvoˇren´ı inspektora . . . . . . . . . . . . . . . . . . . . . . . . . ˇ Zivotn´ ı cyklus inspektora . . . . . . . . . . . . . . . . . . . . . . .
29
Shrnut´ı . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
31
6.2.2 6.3
7 Experimenty 7.1
7.2
23
27
29
33
Nastaven´ı pro experimenty . . . . . . . . . . . . . . . . . . . . . . . . . .
33
7.1.1
Konfigurace prostˇred´ı . . . . . . . . . . . . . . . . . . . . . . . . .
33
7.1.2
Nastaven´ı sc´en´aˇre . . . . . . . . . . . . . . . . . . . . . . . . . . .
35
7.1.3
Metriky k vyhodnocen´ı . . . . . . . . . . . . . . . . . . . . . . . .
37
7.1.4
Nastaven´ı simulace . . . . . . . . . . . . . . . . . . . . . . . . . .
38
Experimenty s jednotlivci
. . . . . . . . . . . . . . . . . . . . . . . . . . viii
39
7.2.1
Chybn´e poˇca´teˇcn´ı nastaven´ı v prostˇred´ı bez inspektora . . . . . .
39
7.2.2
Chybn´e poˇca´teˇcn´ı nastaven´ı v prostˇred´ı s inspektorem . . . . . . .
41
Experimenty s populac´ı . . . . . . . . . . . . . . . . . . . . . . . . . . . .
43
7.3.1
Experimenty na jedn´e lince . . . . . . . . . . . . . . . . . . . . .
43 46
7.4
Experimenty na cel´e s´ıti . . . . . . . . . . . . . . . . . . . . . . . V´ ypoˇcetn´ı a pamˇet’ov´e n´aroky . . . . . . . . . . . . . . . . . . . . . . . .
7.5
Shrnut´ı . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
49
7.3
7.3.2
49
8 Z´ avˇ er
51
Literatura
54
A N´ azev pˇ r´ılohy
I
B Obsah pˇ riloˇ zen´ eho CD
VII
ix
x
Seznam obr´ azk˚ u 4.1
Logaritmick´a z´avislost na ekonomick´em stavu pasaˇz´era . . . . . . . . . .
17
6.1
Z´avislosti mezi platformami . . . . . . . . . . . . . . . . . . . . . . . . .
26
6.2
Stavov´ y diagram ˇzivotn´ıho cyklu pasaˇz´era . . . . . . . . . . . . . . . . .
27
6.3
Stavov´ y diagram prov´adˇen´ı pl´anu . . . . . . . . . . . . . . . . . . . . . .
28
6.4
Diagram deterministick´eho ˇzivotn´ıho cyklu inspektora . . . . . . . . . . .
30
6.5
Diagram nedeterministick´eho ˇzivotn´ıho cyklu inspektora . . . . . . . . .
31
6.6
Zobrazen´ı grafick´eho prostˇred´ı simulaˇcn´ı platformy
. . . . . . . . . . . .
32
7.1
Sch´ema dopravn´ı s´ıtˇe . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
35
7.2
V´ ysledek experimentu pro ekonomicky zaloˇzen´eho pasaˇz´era . . . . . . . .
40
7.3
V´ ysledek experimentu pro ˇcasovˇe zaloˇzen´eho pasaˇz´era . . . . . . . . . . .
40
7.4
V´ ysledek experimentu pro psychologicky zaloˇzen´eho pasaˇz´era . . . . . . .
41
7.5
V´ ysledek v prostˇred´ı s inspektorem na vstupn´ı stanici . . . . . . . . . . .
42
7.6
V´ yvoj poˇctu zakoupen´ ych l´ıstk˚ u v prostˇred´ı jedn´e linky . . . . . . . . . .
44
7.7
Ust´alen´e stavy r˚ uzn´ ych nastaven´ıch situac´ı (jedna linka) . . . . . . . . .
45
7.8
Histogramy poˇctu pasaˇz´er˚ u a jejich zmˇen pl´an˚ u (jedna linka) . . . . . . .
45
Zn´azornˇen´ı situace na pˇrestupn´ı stanici . . . . . . . . . . . . . . . . . . . 7.10 Ust´alen´e stavy r˚ uzn´ ych nastaven´ıch situac´ı (cel´a s´ıt’) . . . . . . . . . . . . 7.11 Histogramy poˇctu pasaˇz´er˚ u a jejich zmˇen pl´an˚ u (cel´a s´ıt’) . . . . . . . . .
47
A.1 V´ yvoj poˇctu zakoupen´ ych l´ıstk˚ u v prostˇred´ı cel´e s´ıtˇe . . . . . . . . . . . .
IV
7.9
xi
48 48
xii
Seznam tabulek 7.1
Tabulka fixn´ıch parametr˚ u prostˇred´ı . . . . . . . . . . . . . . . . . . . . .
34
7.2
Tabulka fixn´ıch parametr˚ u normalizaˇcn´ıch konstant pasaˇz´er˚ u . . . . . . .
37
7.3
Tabulka nastaven´ı experiment˚ u v prostˇred´ı s jednotlivci . . . . . . . . . .
39
7.4
Tabulka nastaven´ı experiment˚ u jedn´e linky . . . . . . . . . . . . . . . . .
43
7.5
Tabulka nastaven´ı experiment˚ u cel´e s´ıtˇe . . . . . . . . . . . . . . . . . . .
46
7.6
Poˇcet proj´ıˇzdˇej´ıc´ıch pasaˇz´er˚ u stanicemi pˇred um´ıstˇen´ım inspektora . . . .
47
7.7
Poˇcet proj´ıˇzdˇej´ıc´ıch pasaˇz´er˚ u stanicemi po um´ıstˇen´ım inspektora . . . . .
47
7.8
V´ ypoˇcetn´ı n´aroky simulace . . . . . . . . . . . . . . . . . . . . . . . . . .
49
A.1 Pˇr´ıjmy pr˚ umˇern´e tˇr´ıˇclenn´e dom´acnosti v LA . . . . . . . . . . . . . . . .
I
A.2 Zmˇeny pl´an˚ u v pˇredv´ıdateln´em prostˇred´ı, 0%, jedna linka . . . . . . . . .
II
A.3 Zmˇeny pl´an˚ u v nepˇredv´ıdateln´em prostˇred´ı, 0%, jedna linka . . . . . . .
II
A.4 Zmˇeny pl´an˚ u v pˇredv´ıdateln´em prostˇred´ı, 1.1%, jedna linka . . . . . . . .
II
A.5 Zmˇeny pl´an˚ u v nepˇredv´ıdateln´em prostˇred´ı, 1.1%, jedna linka . . . . . A.6 Zmˇeny pl´an˚ u v pˇredv´ıdateln´em prostˇred´ı, 0%, cel´a s´ıt’ . . . . . . . . . . A.7 Zmˇeny pl´an˚ u v nepˇredv´ıdateln´em prostˇred´ı, 0%, cel´a s´ıt’ . . . . . . . . . A.8 Zmˇeny pl´an˚ u v pˇredv´ıdateln´em prostˇred´ı, 0.7%, cel´a s´ıt’ . . . . . . . . .
.
III
.
IV
.
V
.
V
A.9 Zmˇeny pl´an˚ u v nepˇredv´ıdateln´em prostˇred´ı, 0.7%, cel´a s´ıt’ . . . . . . . . .
V
xiii
xiv
Kapitola 1 ´ Uvod Tato diplomov´a pr´ace navazuje na v´ yzkum skupiny p˚ usob´ıc´ı na katedˇre poˇc´ıtaˇc˚ u elekˇ eho vysok´eho uˇcen´ı technick´eho v Praze pod veden´ım Ing. Mitrotechnick´e fakulty Cesk´ chala Jakoba, Ph.D.. V´ yzkum t´eto skupiny je zamˇeˇren na v´ yvoj softwarov´ ych prostˇredk˚ u pro simulaci dopravn´ıch syst´em˚ u nesouc´ı n´azev AgentPolis1 (Jakob, M., Moler, Z., Komenda, A., Yin, Z., Jiang, A., X., Johnson, M., P., Pechoucek, M. a Tambe, M., 2012). Tato pr´ace rozˇsiˇruje modul simuluj´ıc´ı syst´em veˇrejn´e dopravy FareEvasion, kter´ y je nadstavbou projektu AgentPolis. Syst´em dopravn´ı pˇrepravy je pro jejich spoleˇcnosti velmi n´akladn´ y. Nav´ıc jako kaˇzd´ y podobn´ y syst´em trp´ı i tento neˇza´douc´ım jevem, kter´ ym je u ´nik zisk˚ u. Tento jev nelze u ´plnˇe odstranit nicm´enˇe dopravn´ı spoleˇcnosti se ho snaˇz´ı eliminovat. Zp˚ usob eliminace u ´niku zisk˚ u vˇzdy z´avis´ı na typu prostˇred´ı, ale vˇzdy je potˇreba m´ıt o nˇem dostatek informac´ı, ze kter´ ych vypl´ yv´a jeho chov´an´ı. S re´aln´ ym probl´emem pˇrepravn´ı kontroly je velmi sloˇzit´e pracovat, vyhodnocovat a vylepˇsovat. Zjednoduˇsen´ı pˇrin´aˇs´ı simulaˇcn´ı technika - modelov´an´ı, kter´e nav´ıc pom´ah´a zlepˇsit pochopen´ı chov´an´ı syst´emu. Pomoc´ı techniky modelov´an´ı je moˇzn´e vytvoˇrit model syst´emu, na kter´em lze prov´adˇet simulace a n´azornˇe uk´azat i skryt´e chov´an´ı syst´emu, jeho reakce na vloˇzen´a data (napˇr. kontroluj´ıc´ı inspektoˇri) nebo m˚ uˇze poskytovat informace, kter´e pˇrispˇej´ı k eliminov´an´ı u ´niku zisku dopravn´ı spoleˇcnosti napˇr. slab´a m´ısta kontroly → rozmist’ov´an´ı inspektor˚ u .... Pouˇzit´ı tradiˇcn´ıch pˇr´ıstup˚ u (dynamick´e syst´emy, diskr´etn´ı ud´alosti) na probl´em pˇrepravn´ı kontroly je velmi obt´ıˇzn´e. S v´ yvojem technologie a n´ar˚ ustem v´ ypoˇcetn´ıho v´ ykonu lze vyuˇz´ıt netradiˇcn´ı simulaˇcn´ı techniku - agentn´ı modelov´an´ı. Pˇr´ıstup agentn´ıho 1
http://agents.fel.cvut.cz/projects/agentpolis/
1
´ KAPITOLA 1. UVOD
2
modelov´an´ı je odliˇsn´ y, nepopisuje jiˇz syst´em jako celek, ale pouze jeho z´akladn´ı jednotku tj. agenta, coˇz b´ yv´a ˇcasto jednoduch´e. Principem techniky je mnoho agent˚ u, kteˇr´ı vz´ajemnou interakc´ı vyv´ıjej´ı chov´an´ı odpov´ıdaj´ıc´ı modelovan´emu syst´emu. V syst´emu pˇrepravn´ı kontroly vystupuje mnoho jednotek, kter´ ymi jsou cestuj´ıc´ı a inspektoˇri. Jednotliv´ı cestuj´ıc´ı interaguj´ı s prostˇred´ım i s inspektory, coˇz pˇr´ımo souvis´ı s pˇr´ıstupem agentn´ıho modelov´an´ı a bylo tedy v´ yhodn´e vyuˇzit´ı pr´avˇe t´eto techniky. Pomoc´ı simulaˇcn´ı techniky agentn´ıho modelov´an´ı byl vytvoˇren model syst´emu pˇrepravn´ı kontroly. Kde popisovan´ y pasaˇz´er, jako autonomn´ı agent, samovolnˇe uvaˇzuje tedy vol´ı zp˚ usob cestov´an´ı v˚ uˇci aktu´aln´ımu prostˇred´ı, na kter´e je schopen se pr˚ ubˇeˇznˇe adaptovat tj. uˇcit se ho. C´ılem bylo pochopit chov´an´ı a adaptaci pasaˇz´er˚ u pˇri mˇen´ıc´ıch se podm´ınk´ach v prostˇred´ı. Pouˇzit´a simulaˇcn´ı technika agentn´ıho modelov´an´ı je popisovan´a v druh´e kapitole t´eto pr´ace, spoleˇcnˇe s jej´ımi v´ yhodami, pouˇzit´ım a porovn´an´ım s ostatn´ımi technikami modelov´an´ı. Pouˇzit´ı t´eto metody je velmi v´ yhodn´e aplikovat na syst´emy, ve kter´ ych jedn´a v´ıce nez´avisl´ ych jednotek jako napˇr. lid´e v dopravn´ım syst´emu. Tˇret´ı kapitola pr´ace struˇcnˇe popisuje souvisej´ıc´ı ˇcl´anky, kter´e byly pouˇzity pro z´ısk´an´ı z´akladn´ıho pohledu do veˇrejn´ ych dopravn´ıch syst´em˚ u. Dalˇs´ı kapitola pr´ace ”Specializace a formalizace probl´emu” popisuje prostˇred´ı dopravn´ıho syst´emu, kter´e pasaˇz´er vn´ım´a, a d´ale popisuje model pasaˇz´era. Model obsahuje vlastnosti pasaˇz´era, prostˇred´ı a cenovou funkci, kterou urˇcuje kolik ho dan´a cesta bude st´at. Vlastnosti agenta a jeho chov´an´ı je definovan´e pomoc´ı sady provediteln´ ych akc´ı a n´aslednˇe jsou pl´anem pops´ana vˇsechna jeho rozhodnut´ı. Kapitola ”Uvaˇzov´an´ı pasaˇz´era”se jiˇz zab´ yv´a, jak´ ym zp˚ usobem se pasaˇz´er rozhoduje pro nejv´ yhodnˇejˇs´ı pl´an a d´ale pokraˇcuje popisem zp˚ usobu jeho adaptace na prostˇred´ı. Pro v´ ybˇer nejv´ yhodnˇejˇs´ıho pl´anu se vyuˇz´ıv´a metody ”Vˇetv´ı a mez´ı”. Praktick´a ˇc´ast pr´ace zaˇc´ın´a popisem implementace vytvoˇren´eho modelu, rozhodov´an´ım a adaptac´ı z pˇredchoz´ıch dvou kapitol. Bylo ovˇeˇreno, ˇze zvolen´ y algoritmus je korektn´ı, a protoˇze se jedn´a o koneˇcn´e prostˇred´ı a algoritmus nekonˇc´ı s prvn´ım nalezen´ ym ˇreˇsen´ım, ale hled´a dalˇs´ı moˇznosti, je tak´e u ´pln´ y. V posledn´ı kapitole ”Exterimenty” se ovˇeˇruje implementace a jsou vyhodnoceny zaj´ımav´e metriky, kter´e ukazuj´ı chov´an´ı modelovan´eho syst´emu v r˚ uzn´ ych stavech prostˇred´ı. Experimenty ovˇeˇruj´ı, ˇze implementace syst´emu dopravn´ı spoleˇcnosti se chov´a intuitivn´ım zp˚ usobem. Dynamika cel´eho syst´emu vˇzdy konverguje do vyv´aˇzen´eho stavu.
Kapitola 2 Agentn´ı modelov´ an´ı Pˇrepravn´ı kontrola se ˇrad´ı mezi zn´am´e syst´emy re´aln´eho svˇeta. Obecnˇe je s tˇemito syst´emy velmi obt´ıˇzn´e nˇekdy dokonce nemoˇzn´e d´ale pracovat, vyhodnocovat je nebo s nimi experimentovat. Pro snaˇzˇs´ı pr´aci s probl´emy re´aln´eho svˇeta se vyuˇz´ıv´a simulaˇcn´ı technika zvan´a modelov´an´ı. Modelov´an´ı zahrnuje proces pˇrev´adˇen´ı syst´emu z re´aln´eho svˇeta do svˇeta simulac´ı v podobˇe modelu, proces rozboru a optimalizace, pˇrev´adˇen´ı ˇreˇsen´ı zpˇet do re´aln´eho syst´emu. Existuj´ı tˇri hlavn´ı pˇr´ıstupy v simulaˇcn´ım modelov´an´ı - dynamick´e syst´emy, diskr´etn´ı ud´alosti a agentn´ı modelov´an´ı1 . Narozd´ıl od prvn´ıch dvou tradiˇcn´ıch pˇr´ıstup˚ u je agentn´ı modelov´an´ı relativnˇe nov´a technika nal´ezaj´ıc´ı ˇsirok´e uplatnˇen´ı.
2.1
Co je agentn´ı modelov´ an´ı
Agentn´ı modelov´an´ı uvaˇzuje syst´em sloˇzen´ y z mnoha z´akladn´ıch jednotek a zab´ yv´a se, oproti jin´ ym pˇr´ıstup˚ um, pouze pops´an´ım t´eto jednotky, kter´a se naz´ yv´a agent. V prostˇred´ı y modelu se z´aroveˇ n pohybuje vˇetˇs´ı mnoˇzstv´ı agent˚ u, kteˇr´ı mezi sebou interaguj´ı2 . Kaˇzd´ agent um´ı samostatnˇe posoudit svoji situaci vzhledem k prostˇred´ı, ve kter´em se pohybuje, a na z´akladˇe mnoˇziny z´akladn´ıch pravidel sestav´ı odpov´ıdaj´ıc´ı rozhodnut´ı. Jednotka, nesouc´ı n´azev agent, by mˇela b´ yt autonomn´ı (jedn´a v prostˇred´ı bez pˇr´ım´eho vlivu ostatn´ıch, ˇr´ıd´ı sv´e akce a vnitˇrn´ı stav), proaktivn´ı (schopnost agenta uj´ımat se iniciativy, ovlivˇ novat prostˇred´ı k dosaˇzen´ı vlastn´ıch c´ıl˚ u), reaktivn´ı (reaguje na zmˇeny v prostˇred´ı), vn´ımaj´ıc´ı 1 2
Dynamic Systems, Discrete Events, Agen-based modeling Jednoduˇse ˇreˇceno, je agentn´ı modelovac´ı syst´em tvoˇren´ y agenty a vztahy mezi nimi.
3
´ ´I KAPITOLA 2. AGENTN´I MODELOVAN
4
okoln´ı prostˇred´ı, schopnost uˇcen´ı, soci´aln´ı vlastnosti (interaguje s ostatn´ımi agenty, kooperac´ı, komunikac´ı, koordinac´ı ˇci konkurenc´ı) atd. (Borshchev, A. a Filippov, A., 2004).
2.2
Srovn´ an´ı pˇ r´ıstup˚ u
Podstatn´ y rozd´ıl mezi pˇr´ıstupy je tvoˇren popisem glob´aln´ıho chov´an´ı neboli dynamiky syst´emu. V modelech pˇr´ıstupu ”dynamick´ ych syst´em˚ u” a ”diskr´etn´ıch ud´alost´ı” existuje objekt, ve kter´em je toto chov´an´ı definov´ano. V ”agentn´ım” pˇr´ıstupu nen´ı ˇza´dn´ y takov´ y objekt, kter´ y by ho s´am definoval, a definuje se pouze lok´aln´ı chov´an´ı tj. chov´an´ı jednotlivce. Glob´aln´ı chov´an´ı vyplyne z jedn´an´ı mnoha (des´ıtek, stovek, tis´ıc˚ u, mili´on˚ u) jednotlivc˚ u ˇzij´ıc´ıch ve stejn´em prostˇred´ı, komunikuj´ıc´ıch mezi sebou a interaguj´ıc´ıch s prostˇred´ım3 (Borshchev, A. a Filippov, A., 2004).
2.3
V´ yhody agentn´ıho modelov´ an´ı
V´ yhody, kter´e agentn´ı modelov´an´ı poskytuje, lze shrnout v n´asleduj´ıc´ım seznamu (Bonabeau, E., 2002): • Agentn´ı modelov´an´ı pokr´ yv´a ˇsirokou ˇsk´alu probl´em˚ u - je flexibiln´ı. • Dokonce i jednoduch´ y agentn´ı model m˚ uˇze projevovat sloˇzit´e (komplexn´ı) vzory chov´an´ı4 , a poskytovat hodnotnou informaci o dynamice re´aln´eho svˇeta, kterou syst´em napodobuje. • Pouˇzit´a technika 5 je jednoduˇse pouˇziteln´a a teoreticky pochopiteln´a a proto z´ıskala znaˇcnou popularitu. • I prost´a pravidla jednotlivce mohou v´est k ucelen´emu chov´an´ı skupiny. • S agentn´ım modelov´an´ım lze vytvoˇrit model, kter´ y se bude jevit jako realistiˇctˇejˇs´ı. 3
Proto se agentn´ı modelov´ an´ı tak´e naz´ yv´ a botton-up modelov´an´ı. Emergence phenomena 5 Technika modelov´ an´ı nejmenˇs´ı ˇc´ asti syst´emu, kter´a je ˇcasto jednoduch´a. 4
ˇ ´I PR ˇ ´ISTUPU 2.4. POUZIT
2.4
5
Pouˇ zit´ı pˇ r´ıstupu
Moˇznost aplikace agentn´ıho modelov´an´ı je ˇsirok´a. Obecnˇe se d´a pouˇz´ıt v probl´emech, kde vystupuje v´ıce samostatn´ ych jednotek, na kter´ ych lze uk´azat kolektivn´ı chov´an´ı tj. na soci´aln´ıch, politick´ ych a ekonomick´ ych oborech. Pouˇzit´ı (Bonabeau, E., 2002) napˇr. u tok˚ u (evakuace, doprava, . . . ), trhu (burza, strategick´e simulace, . . . ), organizace (provozn´ı rizika, organizaˇcn´ı design) a . . . .
2.5
Shrnut´ı
Z v´ yˇse zm´ınˇen´eho lze soudit, ˇze agentn´ı pˇr´ıstup se skuteˇcnˇe hod´ı na probl´emy, v nichˇz jsou modelovan´e jednotky lid´e. Syst´em pˇrepravn´ı kontroly je pˇr´ımo tento druh probl´emu. Autonomn´ı agenti jsou tvoˇreni podle pasaˇz´er˚ u/inspektor˚ u, kteˇr´ı dennˇe prov´adˇej´ı urˇcitou cestu/kontrolu. Pasaˇz´er popisovan´ y v kapitole 4 a jeho rozhodov´an´ı v kapitole 5 m´a povˇedom´ı o prostˇred´ı, reaguje na zmˇeny v nˇem prob´ıhaj´ıc´ı a konkurenˇcnˇe interaguje s inspektorem v podobˇe pˇrepravn´ıch kontrol.
6
´ ´I KAPITOLA 2. AGENTN´I MODELOVAN
Kapitola 3 Souvisej´ıc´ı pr´ ace Pˇrepravn´ı kontrolou se zab´ yvaj´ı jistˇe vˇsechny dopravn´ı spoleˇcnosti, kter´e maj´ı za c´ıl omezit u ´nik zisk˚ u v syst´emu. Pro vytvoˇren´ı odpov´ıdaj´ıc´ı simulace prostˇred´ı pˇrepravn´ı kontroly bylo potˇreba z´ıskat informace o probl´emu, ˇcehoˇz bylo dosaˇzeno prostudov´an´ım nˇekolika prac´ı zab´ yvaj´ıc´ıch se danou problematikou. Tato ˇc´ast obsahuje struˇcn´e reˇserˇse souvisej´ıc´ıch v´ yzkum˚ u. Souvisej´ıc´ı pr´ace, kter´e se zab´ yvaj´ı celou oblast´ı u ´nik˚ u zisk˚ u v pˇrepravn´ım syst´emu, se mohou rozdˇelit do nˇekolika skupin. Prvn´ı skupinu tvoˇr´ı pr´ace, kter´e se zab´ yvaj´ı anal´ yzou probl´emu a na z´akladˇe poznatk˚ u uˇcin´ı opatˇren´ı, druhou jsou ˇcl´anky zab´ yvaj´ıc´ı se jednoduch´ ym n´avrhem rozhodov´an´ı pasaˇz´era.
3.1
Anal´ yza prostˇ red´ı cestuj´ıc´ıch
V pr´aci (Leung, L., F., 2003), kter´a spad´a do prvn´ı skupiny, se autor zab´ yv´a pˇr´ıˇcinami a v´aˇznost´ı u ´niku zisku v Hong Kongsk´em dopravn´ım syst´emu. Autor popisuje probl´emy odbavovac´ıch mechanism˚ u a to jak v uzavˇren´em, tak otevˇren´e syst´emu, d´ale pomoc´ı sestaven´ ych dotazn´ık˚ u z´ısk´a anal´ yzu prostˇred´ı (kdo a proˇc poruˇsuje) a shrnuje r˚ uzn´e typy metod u ´nik˚ u zisku (jak´ ym zp˚ usobem). ˇ anek (Killias, M., Scheidegger, D. a Nordenson, P., 2009) analyzuje Cl´ prostˇred´ı skuteˇcn´ ym experimentem v pˇr´ımˇestsk´e oblasti a to pˇrid´an´ım k n´ahodn´ ym kontrol´am i kontrolu vˇsech cestuj´ıc´ıch v kritick´e dobˇe. V´ ysledky experimentu ukazuj´ı, ˇze jistota trestu funguje jako odstraˇsuj´ıc´ı prostˇredek. 7
´ KAPITOLA 3. SOUVISEJ´IC´I PRACE
8
3.2
Modelov´ an´ı rozhodov´ an´ı cestuj´ıc´ıch
Tato skupina sdruˇzuje ˇcl´anky, kter´e se zab´ yvaj´ı vytvoˇren´ım z´akladn´ıho rozhodov´an´ı ˇ anky (Boyd, C., Martini, Ch., Rickard, J. a Russell, A., 1989), cestuj´ıc´ıch. Cl´ (Kooreman, P., 1993) a (Piliavin, I., Thornton, C., Gartner, R. a Matsueda, R., L., 1986) popisuj´ı vytvoˇren´ı rozhodnut´ı mezi j´ızdou bez j´ızdenky a s n´ı. Ke kaˇzd´e moˇznosti pˇridˇeluj´ı ohodnocen´ı ve formˇe uˇzitku. To znamen´a, ˇze urˇcuj´ı jak´ y uˇzitek jim pˇrinese jej´ı vykon´an´ı. Shodnˇe ˇreˇs´ı volbu nezakoupen´ı j´ızdenky, kterou popisuj´ı jako ˇ anek (Piliavin, I., Thornmoˇznost obdrˇzen´ı pokuty (s urˇcitou pravdˇepodobnost´ı). Cl´ ton, C., Gartner, R. a Matsueda, R., L., 1986) poukazuje na racion´aln´ı volbu, ˇze se vˇzdy zvol´ı v´ yhodnˇejˇs´ı varianta. V ˇcl´anku (Kooreman, P., 1993) vytv´aˇrej´ı sloˇzitˇejˇs´ı model. Do uˇzitku v podobˇe neline´arn´ı funkce pˇridali poˇc´ateˇcn´ı kapit´al a jeho hodnotu poˇc´ıtali jako rozd´ıl uˇzitku pˇred a po proveden´ı dan´e volby. V ˇcl´anku (Boyd, C., Martini, Ch., Rickard, J. a Russell, A., 1989) se po sestaven´ı modelu pasaˇz´era, pokouˇseli zjistit, pro jakou velikost rizika - tedy poˇctu revizor˚ u v dopravn´ım syst´emu se zaruˇc´ı nejvyˇsˇs´ı zisk spoleˇcnosti. Udrˇzov´an´ım urˇcit´e velikosti rizika pˇri pˇrepravˇe m˚ uˇze zaruˇcit vˇetˇs´ı pˇr´ıjem (z pokut) spoleˇcnosti, neˇz kdyby bylo pˇr´ıliˇs inspektor˚ u nebo naopak m´alo.
3.3
Shrnut´ı
´ Unik zisk˚ u je nevyhnutelnˇe spojen se syst´emem veˇrejn´e dopravn´ı pˇrepravy a nem˚ uˇze b´ yt plnˇe eliminov´an. Existuje moˇznost jeho omezen´ı napˇr. n´ahodn´ ymi kontrolami a vysok´ ymi pokutami, coˇz m´a psychologick´ y zastraˇsovac´ı efekt. Uveden´e ˇcl´anky stav´ı svoje modely na znalosti pravdˇepodobnosti kontroly. Maj´ı urˇcenou tuto pravdˇepodobnost z pohledu dopravn´ı spoleˇcnosti, kter´a je d´ana poˇctem inspektor˚ u. Nezab´ yvaj´ı se t´ım, jak tuto pravdˇepodobnost vid´ı pasaˇz´eˇri. ˇ ank˚ Cl´ u zab´ yvaj´ıc´ıch se u ´nikem ze zisk˚ u v pˇrepravn´ı kontrole nen´ı mnoho. Informace z nich z´ıskan´e jsou ovˇsem dostaˇcuj´ıc´ı k vytvoˇren´ı sloˇzitˇejˇs´ıho1 modelu, jehoˇz n´avrhem tzn. popisem prostˇred´ı a rozhodov´an´ı pasaˇz´er˚ u se tato pr´ace zab´ yv´a v n´asleduj´ıc´ıch dvou kapitol´ach, kter´e popisuj´ı nejenom moˇznost volby j´ızdenky, ale i obch´azen´ı nebezpeˇcn´ ych stanic/´ usek˚ u a dokonce volby ˇcasu v´ yjezdu. 1
Obsahuje v´ıce moˇznost´ı voleb (j´ızdenka, trasa nebo ˇcas v´ yjezdu), kter´e racion´aln´ı pasaˇz´er m´a.
Kapitola 4 Specifikace a formalizace probl´ emu Pouˇz´ıvan´e agentn´ı modelov´an´ı vyˇzaduje pops´an´ı chov´an´ı agenta a prostˇred´ı, ve kter´em se pohybuje. V t´eto kapitole se popisuje prostˇred´ı, vlastnosti agenta a jeho chov´an´ı definovan´e pomoc´ı sady provediteln´ ych akc´ı. N´aslednˇe jsou pl´anem popsan´a vˇsechna jeho rozhodnut´ı a nakonec pasaˇz´erova funkce ceny, kter´a ohodnocuje zvolen´ y pl´an tj. sekvence jeho rozhodnut´ı.
4.1
Reprezentace prostˇ red´ı pasaˇ z´ era
Pasaˇz´er se pohybuje v prostˇred´ı veˇrejn´e dopravn´ı s´ıtˇe s jedn´ım typem dopravn´ıho prostˇredku, kter´ ym je metro. Jedna linka metra se skl´ad´a ze stanic a odliˇsnˇe dlouh´ ych cest (spoj˚ u) mezi nimi. Dopravn´ı s´ıt’ se vˇetˇsinou skl´ad´a z v´ıce linek uzp˚ usoben´ ych tak, aby se v nˇejak´e stanici prot´ınaly. Takov´ ym stanic´ım, kde se prot´ınaj´ı dvˇe a v´ıce linek, se ˇr´ık´a pˇrestupn´ı stanice. Reprezentace prostˇred´ı neboli sestaven´ı jeho modelu se inspirovalo sekvenˇcn´ımi rozhodovac´ımi probl´emy (Vidal, J., M., 2010), kde stejnˇe jako zde celkov´a cena cesty1 pasaˇz´era z´avis´ı na jeho sekvenˇcn´ıch rozhodnut´ıch (napˇr. volba trasy, j´ızdenky a kdy pojede). V pasaˇz´erovˇe prostˇred´ı v definici 4.1 je stanoven´e deterministick´e prov´adˇen´ı akc´ı a nedeterministick´e z´ısk´av´an´ı odmˇen.
1
Pasaˇz´era jeho cesta nˇeco stoj´ı. U sekvenˇcn´ıch rozhodovac´ıch probl´em˚ u se tato funkce naz´ yv´a uˇzitek.
9
´ KAPITOLA 4. SPECIFIKACE A FORMALIZACE PROBLEMU
10
Definice 4.1 (Prostˇ red´ı pasaˇ z´ era): Prostˇred´ı pasaˇz´era
ε
je
definov´ano
n-tic´ı
⟨S, A, T, pk , t, cj , cp ⟩: • S . . . mnoˇzina vˇsech stanic metra, • A . . . mnoˇzina vˇsech moˇznost´ı dopravy mezi stanicemi, • T (s, a, s′ ) . . . pˇrechodov´a funkce, kde s ∈ S a ∈ A a s′ ∈ S, • pk (a, t) . . . pravdˇepodobnost kontrolov´an´ı pˇri akci a ∈ A v ˇcase t, • t . . . ˇcas, • cj . . . cena j´ızdenky a • cp . . . velikost pokuty.
I
Pˇrechodov´a funkce s oznaˇcen´ım T (s, a, s′ )2 v tomto prostˇred´ı ud´av´a pouˇzitelnost akc´ı, kter´a z d˚ uvodu deterministick´ ych akc´ı nab´ yv´a pouze dvou hodnot. Stoprocentn´ı pravdˇepodobnosti, pokud je akce a dostupn´a ve stavu s, a nulov´e, pokud nen´ı dostupn´a, tzn. vˇzdy kdyˇz lze akci prov´est, tak u ´spˇeˇsnˇe provede napˇr. pasaˇz´er se pˇresune do jin´e stanice. Tedy obor hodnot nab´ yv´a pouze {0, 1} a definiˇcn´ı obor pˇrechodov´e funkce je tvoˇren mnoˇzinou vˇsech stav˚ u a akc´ı. Pravdˇepodobnost chycen´ı pk (a, t) zast´av´a v modelu funkci nedeterministick´eho z´ısk´av´an´ı odmˇeny za dosaˇzen´ı stavu, neboli po u ´spˇeˇsn´em dokonˇcen´ı akce. Pasaˇz´er m˚ uˇze b´ yt na sv´e cestˇe s urˇcitou pravdˇepodobnost´ı kontrolov´an inspektorem na vstupn´ı, v´ ystupn´ı a pˇrestupn´ı stanici nebo v soupravˇe mezi nimi. Protoˇze pasaˇz´er m˚ uˇze z´ıskat jen z´apornou odmˇenu pˇri pokutov´an´ı ˇci koupi j´ızdenky, pojem odmˇena jiˇz nen´ı na m´ıstˇe a byl nahrazen pojmem penalizace. Aby nebylo nutn´e rozliˇsovat kontrolu ve stanici a v soupravˇe, jsou vˇsechny moˇznosti sjednoceny pouze do akc´ı, kter´ ymi jsou vstup a v´ ystup ze stanice, pˇrejezdy mezi nimi a pˇrechody pˇeˇsky viz. d´ale Akce v prostˇred´ı pasaˇz´era sekce 4.1.2. Pasaˇz´er potom dost´av´a penalizaci po u ´spˇeˇsn´em dokonˇcen´ı kaˇzd´e akce. Velikost penalizace z´avis´ı na pasaˇz´erovˇe volbˇe. Pokutov´an´ı se totiˇz net´ yk´a pasaˇz´era, kter´ y m´a platnou j´ızdenku, ten obdrˇz´ı nˇejakou malou penalizaci (dle ceny j´ızdenky). Pokud pasaˇz´er nem´a platnou j´ızdenku, obdrˇz´ı s urˇcitou pravdˇepodobnost´ı, kter´a odpov´ıd´a 2
U sekvenˇcn´ıch rozhodovac´ıch probl´em˚ u oznaˇcuje pravdˇepodobnost dosaˇzen´ı stavu s′ , pokud se agent
nach´az´ı ve stavu s a provede akci a neboli pravdˇepodobnost u ´spˇeˇsn´eho dokonˇcen´ı akce (Vidal, J., M., 2010).
ˇ ´I PASAZ ˇERA ´ 4.1. REPREZENTACE PROSTRED
11
pokutov´an´ı, vyˇsˇs´ı penalizaci (velikost pokuty) nebo nez´ısk´a ˇza´dnou penalizaci, pokud nebyl kontrolov´an. Z pk (a, t) je vidˇet, ˇze penalizace jsou v pasaˇz´erovˇe prostˇred´ı dynamick´e tj. promˇenliv´e v ˇcase. Tedy je zˇrejm´ y pravdˇepodobnostn´ı charakter penalizaˇcn´ı funkce, pravdˇepodobnost, s jakou bude pasaˇz´er kontrolov´an v ˇcase t ve stavu s nebo pˇri prov´adˇen´ı akce a. D´ale je definov´ana mnoˇzina stav˚ u (stanic), do kter´ ych se m˚ uˇze dostat pomoc´ı jednotliv´ ych akc´ı (pˇrechod˚ u/pˇrejezd˚ u).
4.1.1
Stavy v prostˇ red´ı pasaˇ z´ era
Stavy v pasaˇz´erovˇe prostˇred´ı reprezentuj´ı skuteˇcn´e stanice veˇrejn´e dopravn´ı s´ıtˇe a stavy pˇred vstupy do jednotliv´ ych stanic. Aby jednotliv´e linky a stanice byly dobˇre identifikovateln´e a pozdˇeji tak´e pro lepˇs´ı identifikaci akc´ı pˇrechod˚ u mezi r˚ uzn´ ymi stanicemi/linkami, je zavedeno znaˇcen´ı stav˚ u INix ve stanici a OU Tix pˇred stanic´ı, kde i je poˇradov´e ˇc´ıslo stanice na lince x. Mnoˇzina vˇsech stav˚ u je vyj´adˇrena rovnic´ı 4.1. S ≡ {INix , OU Tix | x ∈ L ∧ n ∈ {1 . . . len(x)}},
(4.1)
kde L je mnoˇzina vˇsech linek metra (napˇr. L ≡ {a, b, c, . . . }) a funkce len(x) je d´elka linky metra x.
4.1.2
Akce v prostˇ red´ı pasaˇ z´ era
Z pasaˇz´erova pohledu existuj´ı ˇctyˇri akce, kter´e m˚ uˇze vykon´avat na cestˇe do c´ılov´eho stavu. Prvn´ı MOVE - pˇrejezd metrem z jedn´e stanice do n´asleduj´ıc´ı stanice, druh´a WALK pˇrechod pˇeˇsky mezi r˚ uzn´ ymi stanicemi, tˇret´ı ENTRY - vstup do stanice a posledn´ı EXIT - v´ ystup ze stanice. V´ ysledn´a mnoˇzina vˇsech akc´ı, kter´a je zn´azornˇena v rovnici 4.2, je sjednocen´ım tˇechto ˇctyˇr mnoˇzin akc´ı. A ≡ AMOVE ∪ AWALK ∪ AENTRY ∪ AEXIT ,
(4.2)
kde AMOVE je mnoˇzina akc´ı urˇcuj´ıc´ı j´ızdy metrem, AWALK je mnoˇzina akc´ı urˇcuj´ıc´ı pˇrechody mezi stanicemi, AENTRY je mnoˇzina akc´ı definuj´ıc´ı vstup do stanice a AEXIT je mnoˇzina akc´ı definuj´ıc´ı v´ ystupy ze stanic. Kaˇzd´a akce je pevnˇe sv´az´ana s konkr´etn´ım stavem. Pokud agent bude cht´ıt pouˇz´ıt akci, kter´a je spojena s jin´ ym stavem neˇz, ve kter´em se pasaˇz´er nach´az´ı, pak se akce
´ KAPITOLA 4. SPECIFIKACE A FORMALIZACE PROBLEMU
12
neprovede viz. pˇrechodov´ y model v´ yˇse.
4.1.3
Mnoˇ zina akc´ı j´ızdy metrem
Mnoˇzina akc´ı je definov´ana v rovnici 4.3 a urˇcuje pro vˇsechny linky metra pˇrejezdy mezi vˇsemi stanicemi na dan´e lince. Jednotliv´a akce MOVE je urˇcena z´apisem x x MOVE{INix → INi+1 } a definuje pˇrejezd metrem ze stavu INix do stavu INi+1 na
lince x.
x AMOVE ≡ {MOVE{INix → INi+1 }, x MOVE{INi+1 → INix } | x ∈ L ∧ i ∈ {1 . . . len(x)}}, (4.3)
kde L je mnoˇzina vˇsech linek metra a funkce len(x) je d´elka linky metra x.
4.1.4
Mnoˇ zina akc´ı pˇ rechodu mezi stanicemi
Pˇrechody mezi stanicemi tvoˇr´ı mnoˇzinu akc´ı viz. rovnice 4.4, kter´e jsou definov´any WALK{INix → INjy } s vyuˇzit´ım pˇrestupu a WALK{OU Tix → OU Tjy } pro venkovn´ı pˇrechod. Obˇe akce znamenaj´ı pˇrechod mezi stavem INix /OU Tix a stavem INjy /OU Tjy . Linky x a y pro venkovn´ı pˇrestup mohou b´ yt identick´e, ovˇsem potom mus´ı b´ yt r˚ uzn´e stanice i a j (i ̸= j). Moˇznost, kdy x = y a i = j je pro agenta neracion´aln´ı volba, kter´a by znamenala z˚ ustat ve stejn´em stavu. Akce vnitˇrn´ıho pˇrestupu je pot´e umoˇznˇena pouze na pˇrestupn´ı stanici.
AWALK ≡ {WALK{INix → INjy } | P (INix , INjy ) ̸= ∅} ∪ ∪ {WALK{OU Tix → OU Tjy } | x ∈ L ∧ y ∈ L∧ ∧ i ∈ {1 . . . len(x)} ∧ j ∈ {1 . . . len(y)}}, (4.4) kde L je mnoˇzina vˇsech linek metra, funkce len(x) je d´elka linky metra x, funkce P (INix , INjy ) vrac´ı pˇrestupy mezi stanicemi INix a INjy , pro x ∈ L ∧ y ∈ L ∧ i ∈ {1 . . . len(x)} ∧ j ∈ {1 . . . len(y)}.
ˇERA ´ 4.2. VLASTNOSTI PASAZ
4.1.5
13
Mnoˇ zina akc´ı vstup do stanice
Vˇsechny vstupy do stanic (akce ENTRY{OU Tix }) tvoˇr´ı dalˇs´ı mnoˇzinu akc´ı viz. rovnice 4.5, kter´e pasaˇz´er m˚ uˇze pouˇz´ıt. AENTRY ≡ {ENTRY{OU Tix } | x ∈ L ∧ i ∈ {1 . . . len(x)}},
(4.5)
kde L je mnoˇzina vˇsech linek metra a funkce len(x) je d´elka linky metra x.
4.1.6
Mnoˇ zina akc´ı v´ ystup ze stanice
Pasaˇz´er tak´e mus´ı m´ıt moˇznost nˇejak´ ym zp˚ usobem vystoupit z kaˇzd´e stanice metra. Tuto moˇznost pokr´ yv´a mnoˇzina akc´ı (akce EXIT{INix }) definuj´ıc´ı v´ ystup ze stanice v rovnici 4.6. AEXIT ≡ {EXIT{INix } | x ∈ L ∧ i ∈ {1 . . . len(x)}},
(4.6)
kde L je mnoˇzina vˇsech linek metra a funkce len(x) je d´elka linky metra x.
4.2
Vlastnosti pasaˇ z´ era
Vlastnosti pasaˇz´era, zapsan´e v definici 4.2, jsou tvoˇreny poˇc´ateˇcn´ı stanic´ı, tedy m´ıstem odkud bude pl´anovat svoji cestu, c´ılovou stanic´ı, moˇzn´ ym ˇcasem odjezdu, nutn´ ym ˇcasem dojezdu a jeho ekonomick´ ym stavem3 . Ekonomick´ y stav je pasaˇz´er˚ uv d˚ uleˇzit´ y parametr pro rozhodnut´ı, zda si koup´ı ˇci nekoup´ı j´ızdenku na svoji trasu. Definice 4.2 (Vlastnosti pasaˇ z´ era): Vlastnosti pasaˇz´era α jsou definov´any pˇetic´ı ⟨s0 , g, k, t0 , td ⟩ • s0 . . . poˇc´ateˇcn´ı stanice, • g . . . c´ılov´a stanice, • k . . . ekonomick´ y stav pasaˇz´era (kapit´al), • t0 . . . moˇzn´ y ˇcas v´ yjezdu, 3
Inspirov´ ano ˇcl´ ankem (Kooreman, P., 1993).
´ KAPITOLA 4. SPECIFIKACE A FORMALIZACE PROBLEMU
14
• td . . . ˇcas, kdy mus´ı pasaˇz´er bezpodm´ıneˇcnˇe dorazit do c´ılov´e stanice.
4.3
I
Pl´ an
Pl´an pasaˇz´era pˇredstavuje jiˇz konkr´etn´ı ”cestovn´ı pl´an”, kter´ y si pasaˇz´er zvol´ı. Obsahuje zvolenou trasu, kterou pojede, volbu ˇcasu v´ yjezdu a typ j´ızdenky. Typem j´ızdenky se mysl´ı, zda pasaˇz´er zvol´ı jednu z moˇznost´ı - koupit ˇci nekoupit j´ızdenku. Volba ˇcasu v´ yjezdu je pro pasaˇz´era tak´e d˚ uleˇzit´a, m˚ uˇze se tak vyhnout nˇejak´ ym kontrol´am, kter´e jsou ˇcasovˇe promˇenn´e (napˇr. nˇekteˇr´ı inspektoˇri se pohybuj´ı souˇcasnˇe se soupravou), a tedy pˇr´ıpadn´emu pokutov´an´ı. Jak jiˇz bylo zm´ınˇeno, prostˇred´ı je ˇcasovˇe promˇenn´e, v ˇcase se mˇen´ı pravdˇepodobnosti kontrol, protoˇze inspektoˇri nejsou st´ale na stejn´em m´ıstˇe. Definice 4.3 (Reprezentace trasy): Reprezentace trasy tr je uspoˇr´adan´a n-tice (a0 , a1 , . . . , an ), kde ∀i ∈ {0, . . . , n} : ai ∈ {MOVE, WALK, ENTRY, EXIT}
Definice 4.4 (Cestovn´ı pl´ an): Cestovn´ı pl´an π je definov´an trojic´ı ⟨tr, ts , j⟩, kde • tr . . . zvolen´a trasa (definice 4.3), • ts . . . zvolen´ y ˇcas v´ yjezdu a • j . . . volba j´ızdenky ANO/NE.
I
V pˇr´ıkladu reprezentace trasy 4.7 akce WALK{INa6 → INc4 } znamen´a pˇrechod pˇeˇsky na pˇrestupn´ı stanici. Obch´azen´ı nebo pˇrech´azen´ı mezi stanicemi se v reprezentaci projev´ı v´ ystupem ze stanice pˇred ch˚ uz´ı a potom vstup zpˇet do stanice: {. . . EXIT{INa6 }, WALK{OU Ta6 → OU Tc4 }, ENTRY{OU Ta6 }, . . . }.
{ENTRY{OU Ta3 }, MOVE{INa3 → INa4 }, . . . . . . , WALK{INa6 → INc4 }, MOVE{INc4 → INc5 }, . . . , EXIT{INc8 }} (4.7)
ˇ ´I CESTOVN´IHO PLANU ´ 4.4. MODEL OCENEN
4.4
15
Model ocenˇ en´ı cestovn´ıho pl´ anu
Model pasaˇz´era lze jednoduˇse pochopit jako ”agenta” vykon´avaj´ıc´ıho cestu v transportn´ım syst´emu veˇrejn´e dopravy. Jeho u ´kolem je zvolit trasu, kterou pojede, ˇcas a jak´ ym zp˚ usobem pojede ⇒ pl´an. Zaveden´ı cenov´e hodnoty pl´anu neboli pasaˇz´erov´ ych rozhodnut´ı v dan´em prostˇred´ı a vlastnostech vyjadˇruje rovnice 4.8. Pasaˇz´er podle cenov´e hodnoty m˚ uˇze rozliˇsovat kvalitu r˚ uzn´ ych pl´an˚ u a zvolit ten nejv´ yhodnˇejˇs´ı. U (π | α, ε) ≡ U (⟨tr, ts , j⟩ | ⟨s0 , g, k, t0 , td ⟩, ⟨S, A, T, pk , t, cj , cp ⟩)
(4.8)
Pokud pasaˇz´er m´a pl´an trasy v prostˇred´ı, pak m˚ uˇze spoˇc´ıtat pravdˇepodobnost kontroly pˇri vykon´av´an´ı tohoto pl´anu. Pravdˇepodobnost pokutov´an´ı je doplˇ nkem k pravdˇepodobnosti situace, kdy pasaˇz´er na cel´e cestˇe nebude kontrolov´an. Tato pravdˇepodobnost se poˇc´ıt´a jako souˇcin pravdˇepodobnost´ı4 nekonrolov´an´ı pˇri vykon´av´an´ı jednotliv´e akce. Jelikoˇz v prostˇred´ı je jen pravdˇepodobnost kontroly pˇri prov´adˇen´ı akce, jedn´a se zase o doplnˇek k t´eto pravdˇepodobnosti v ˇcase t. Cel´ y vztah je zn´azornˇen´ı v rovnici 4.9 p(kontroly | π, ε) = 1 −
n ∏
(1 − pk (ai , t)),
(4.9)
i=1
kde n je d´elka trasy v pl´anu a pk (a, t) ∈ ε je pravdˇepodobnost kontroly pˇri akci a v ˇcase t viz. definice prostˇred´ı 4.1 a ai ∈ π. V´ yˇse zm´ınˇen´a cenov´a funkce definov´ana v rovnici 4.10, pro rozliˇsov´an´ı kvality pl´an˚ u, je rozdˇelena do jednotliv´ ych sloˇzek urˇcuj´ıc´ı ˇcasovou, ekonomickou a psychologickou sloˇzku pl´anu.
U (π) = Ut (π) + Ue (π) + Up (π),
(4.10)
kde Ut tvoˇr´ı ˇcasovou sloˇzku ceny, Ue ekonomickou sloˇzku a Up psychologickou sloˇzku.
4.4.1
ˇ Casov´ a sloˇ zka ceny
Cestov´an´ı mezi stanicemi (prov´adˇen´ı akce) zabere pasaˇz´erovi urˇcit´ y ˇcas. S moˇznost´ı obch´azen´ı nebezpeˇcn´ ych stanic mus´ıme zaˇclenit i ˇcasovou penalizaci, doby str´aven´e na cestˇe pˇri pˇresunu. Kdyby v modelu nebyla zaˇclenˇena ˇcasov´a penalizace, bylo by pro 4
Jedn´ a se o nez´avisl´e jevy, proto souˇcin.
´ KAPITOLA 4. SPECIFIKACE A FORMALIZACE PROBLEMU
16
ˇ pasaˇz´era nejv´ yhodnˇejˇs´ı napˇr. doj´ıt celou cestu nebo obej´ıt podstatnˇe velk´ yu ´sek. Casov´ a penalizace je urˇcena rovnic´ı 4.11. { Ut (π) =
0 + τ (tr) − t0 , td − τ (tr) ≥ 0 ρ + τ (tr) − t0 , td − τ (tr) < 0
,
(4.11)
kde tr ∈ π je trasa, t0 je doba moˇzn´eho startu, τtr je doba doraˇzen´ı do c´ıle v rovnici 4.12, td je nejzazˇs´ı term´ın doraˇzen´ı do c´ıle a pokuta ρ za nedodrˇzen´ı term´ınu.
τ (tr) = t0 +
n ∑
t ai ,
(4.12)
i=1
kde tr je trasa tj. posloupnost akc´ı, n d´elka pl´anu a tai je doba prov´adˇen´ı akce ai ∈ tr.
4.4.2
Ekonomick´ a sloˇ zka ceny
Hypot´ezou, proˇc jezd´ı lid´e bez j´ızdenek, se stal jejich ekonomick´ y stav (neboli ”poˇca´teˇcn´ı”majetek dan´eho jedince), cena j´ızdenky a velikost moˇzn´e pokuty. Tento fakt je moˇzno zohlednit v modelu zaveden´ım denn´ıho rozpoˇctu pasaˇz´era v hodnotˇe k. Vytvoˇren´a funkce ceny je z´avisl´a na ekonomick´em stavu jedince a zn´azornˇen´a v rovnici 4.13. Ue′ (k) =
log(k) , NORM EKO
(4.13)
kde NORM EKO je normalizaˇcn´ı konstanta, kter´a je urˇcena z pr˚ uzkumu ekonomick´e situace simulovan´eho prostˇred´ı. Cena ekonomick´e sloˇzky je z´avisl´a na v´ ybˇeru j´ızdenky. Pokud si pasaˇz´er koup´ı j´ızdenku, tak se v uˇzitkov´e hodnotˇe nepoˇc´ıt´a s pravdˇepodobnost´ı pokutov´an´ı. Naopak se s n´ı poˇc´ıt´a, kdyˇz si pasaˇz´er j´ızdenku nekoup´ı. Moˇznosti ekonomick´e sloˇzky jsou zn´azornˇeny v rovnici 4.14. { UE (π | α, ε) =
Ue′ (k) − Ue′ (k − cj ), iff π ≡ ⟨tr, ts , true⟩
∑ a∈A
pk (a, t) · (Ue′ (k) − Ue′ (k − cp )), iff π ≡ ⟨tr, ts , true⟩
,
(4.14)
kde pk (a, t) ∈ ε je pravdˇepodobnost chycen´ı pˇri prov´adˇen´ı akce a ∈ tr, trasa tr ∈ π, k je ekonomick´ y stav pasaˇz´era, cj je cena j´ızdenky a Ue′ je funkce ekonomick´e z´avislosti popsan´e v´ yˇse.
ˇ ´I CESTOVN´IHO PLANU ´ 4.4. MODEL OCENEN
17
ˇ ast cenov´e funkce UE s cestovn´ım pl´anem π ≡ ⟨tr, ts , true⟩ (s j´ızdenkou) v rovnici C´ 4.14 vych´az´ı z kombinace uˇzitkov´e funkce U (N E) = −cj 5 (Boyd, C., Martini, Ch., Rickard, J. a Russell, A., 1989) a z´avislosti pasaˇz´era na jeho ekonomick´em stavu. Je vidˇet, ˇze hodnota cenov´e funkce n´akupu l´ıstku je rozd´ıl hodnot ekonomick´eho stavu pˇred n´akupem a po n´akupu j´ızdenky (Kooreman, P., 1993). Ze znalosti vztahu ekonomick´eho stavu viz. rovnice 4.13 vypl´ yv´a jej´ı logaritmick´a z´avislost. Tedy pokud je pasaˇz´er ekonomicky slabˇs´ı je pro nˇeho tento rozd´ıl vyˇsˇs´ı, neˇz kdyby pasaˇz´er mˇel vyˇsˇs´ı poˇc´ateˇcn´ı kapit´al (viz. obr´azek 4.16 ). Pokud tuto myˇslenku pˇrevedeme do skuteˇcnosti, kaˇzd´eho v´ıce mrz´ı, kdyˇz plat´ı pomˇerovˇe vˇetˇs´ı ˇca´stku z kapit´alu a naopak m´enˇe tr´ap´ı, kdyˇz mus´ı zaplatit ˇca´st menˇs´ı.
5 4.5 4
Uzitek [−]
3.5 3 2.5 2 1.5 1 0.5 0
0
20
40
60
80
100
Kapital [$]
Obr´azek 4.1: Logaritmick´a z´avislost na ekonomick´em stavu pasaˇz´era
Druh´a ˇca´st cenov´e funkce UE s cestovn´ım pl´anem π ≡ ⟨tr, ts , true⟩ (bez j´ızdenky neboli riskov´an´ı) tvoˇr´ı souˇcet cen prov´adˇen´ ych akc´ı z trasy cestovn´ıho pl´anu. Cena prov´adˇen´ı akce byla odvozena z uˇzitkov´e funkce U (E) = −pr · cp 7 (Boyd, C., Martini, Ch., 5
Uˇzitkov´ a hodnota j´ızdy s j´ızdenkou z ˇcl´anku (Boyd, C., Martini, Ch., Rickard, J. a Russell, A.,
1989). 6 Na obr´azku je logaritmick´ a z´avislost na ekonomick´em stavu pasaˇz´era a ukazuje, jak stejn´ y pokles kapit´alu znamen´a jinou cenu. 7 Uˇzitkov´ a hodnota riskov´ an´ı v ˇcl´ anku (Boyd, C., Martini, Ch., Rickard, J. a Russell, A., 1989). Kde pr je pravdˇepodobnost kontroly.
´ KAPITOLA 4. SPECIFIKACE A FORMALIZACE PROBLEMU
18
Rickard, J. a Russell, A., 1989) a v rovnici byla nahrazena velikost pokuty rozd´ılem ekonomick´eho stavu pˇred a po zaplacen´ı pokuty pˇri prov´adˇen´ı akce.
4.4.3
Psychologick´ a sloˇ zka ceny
D˚ uleˇzitou vlastnost´ı pˇri rozhodov´an´ı o koupi j´ızdenky je psychika pˇri pˇr´ıpadn´em pokutov´an´ı. Zpravidla ˇc´ım vyˇsˇs´ı postaven´ı (tj. ekonomick´ y stav) pasaˇz´er m´a, t´ım je pro nˇeho v´ıce nepˇr´ıjemn´e, ˇze by mohl b´ yt vidˇen pˇri pokutov´an´ı. Tento fakt je vyj´adˇren´ y v rovnici 4.15 a takt´eˇz z´avis´ı na pasaˇz´erovˇe volbˇe l´ıstku. { 0, pokud π ≡ ⟨tr, ts , true⟩ UPSY (π | α, ε) = , 2 k , pokud π ≡ ⟨tr, ts , f alse⟩ NORM PSY
(4.15)
kde k je ekonomick´ y stav a N ORM P SY je normalizaˇcn´ı konstanta urˇcena z pr˚ uzkumu prostˇred´ı.
4.5
Shrnut´ı
V popisu modelu pasaˇz´era se uvaˇzuj´ı pouze jednor´azov´e j´ızdenky. Dlouhodob´e j´ızdenky by znamenaly sloˇzitˇejˇs´ı model. Dalˇs´ım zjednoduˇsen´ım modelu je nez´avislost jednotliv´ ych cest tzn. ˇze pasaˇz´er˚ uv kapit´al pro dalˇs´ı cestu nen´ı ovlivnˇen´ y zakoupen´ım j´ızdenky nebo pokutov´an´ım. Samostatn´ y model popsan´ y v t´eto kapitole nestaˇc´ı pro rozhodov´an´ı pasaˇz´era. V n´asleduj´ıc´ı kapitole se popisuje model spoleˇcnˇe s uˇcen´ım se prostˇred´ı. Jeho dobr´a znalost je d˚ uleˇzit´a pro tvorbu nejv´ yhodnˇejˇs´ıho pl´anu.
Kapitola 5 Uvaˇ zov´ an´ı pasaˇ z´ era V pˇredch´azej´ıc´ı kapitole byl pops´an model pasaˇz´era a jeho prostˇred´ı. To samostatnˇe nestaˇc´ı pro pasaˇz´erovo samovoln´e jedn´an´ı a je potˇreba mechanismus k tvorbˇe posloupnosti rozhodnut´ı a dosaˇzen´ı c´ıl˚ u pasaˇz´era. Mechanismem, pˇri v´ ybˇeru sekvenc´ı rozhodnut´ı, je pl´anovac´ı algoritmus, kter´ y by pasaˇz´erovi staˇcil v pˇr´ıpadˇe, ˇze by znal vˇsechny pˇresn´e pravdˇepodobnosti kontroly. Bez pˇresn´ ych znalost´ı je pro spr´avn´e chov´an´ı poˇzadov´ana schopnost adaptace na prostˇred´ı, kter´a se popisuje na konci kapitoly.
5.1
Pl´ anov´ an´ı trasy
Pasaˇz´er na urˇcit´e ekonomick´e u ´rovni mus´ı kaˇzd´ y den vykonat urˇcit´e cesty. Na kaˇzdou cestu m´a vyhrazen´ y ˇcasov´ y interval ⟨t0 , td ⟩. Dobu, kdy nejdˇr´ıve m˚ uˇze zaˇc´ıt svoji cestu t0 a dobu td , kdy mus´ı b´ yt bezpodm´ıneˇcnˇe v c´ılov´em m´ıstˇe. Pasaˇz´er mus´ı napl´anovat trasu z m´ısta A (poˇc´ateˇcn´ıho stanice) do m´ısta B (c´ılov´e stanice) tak, aby cena jeho pl´anu byla minim´aln´ı tj. minimalizuje cenovou funkci. Z kapitoly 4 sekce Cestovn´ı pl´an jsou zn´amy moˇznosti rozhodov´an´ı pasaˇz´era: zda si koup´ı/nekoup´ı j´ızdenku, jakou zvol´ı trasu a ˇcas, kdy zaˇcne s realizac´ı pl´anu. Prostˇred´ı, ve kter´em se pasaˇz´er pohybuje definovan´e v sekci Reprezentace prostˇred´ı pasaˇz´era kapitoly 4, m˚ uˇze b´ yt tak´e n´azornˇe pops´ano stavov´ ym diagramem. Pasaˇz´er m´a mnoˇzinu dostupn´ ych stav˚ u, kde kaˇzd´ y z nich obsahuje mnoˇzinu akc´ı, kter´e m˚ uˇze v dan´em stavu aplikovat (vyuˇz´ıt). 19
ˇ AN ´ ´I PASAZ ˇERA ´ KAPITOLA 5. UVAZOV
20
Pasaˇz´er nehled´a jen nˇejak´ y pl´an, se kter´ ym se dostane z poˇca´teˇcn´ıho stavu do c´ılov´eho. Hled´a nejlepˇs´ı moˇzn´ y pl´an, kter´ y m´a nejmenˇs´ı hodnotu cenov´e funkce. Prostˇred´ı je pozorovateln´e, pasaˇz´er m´a tedy k dispozici nˇejakou informaci o prostˇred´ı, ve kter´em se pohybuje.
5.1.1
Metoda vˇ etv´ı a mez´ı
K pl´anov´an´ı v pozorovateln´em stavov´em diagramu se vyuˇz´ıvaj´ı informovan´e pl´anovac´ı algoritmy. Nejl´epe vyhovuj´ıc´ım pl´anovac´ım algoritmem pro tento probl´em je metoda Vˇetv´ı a Mez´ı (Clausen, J., 1999) s prohled´av´an´ım do hloubky. V´ yhoda metody vˇetv´ı a mez´ı s prohled´av´an´ım do hloubky je ta, ˇze rychle1 najde nˇejak´e ˇreˇsen´ı, podle kter´eho m˚ uˇze d´ale proˇrez´avat vˇetve prohled´avan´eho stromu, ve kter´em se tvoˇr´ı jiˇz jen horˇs´ı pl´any, neˇz je doposud nalezen´ y. Toto proˇrez´av´an´ı znaˇcnˇe urychluje v´ ypoˇcet a tedy nalezen´ı nejlepˇs´ı moˇzn´e trasy resp. pl´anu. Existuje nˇekolik pˇredpoklad˚ u, kter´e je tˇreba splnit pro nasazen´ı algoritmu na tento pl´anovac´ı probl´em. Jedn´ım je pˇredpoklad diskretizovan´eho ˇcasu a druh´ ym je fakt, ˇze z kaˇzd´eho stavu lze zrekonstruovat pl´an. Jeho vykon´an´ım se agent dostane do poˇzadovan´eho stavu, a proto se kaˇzd´ y stav m˚ uˇze prezentovat jako pl´an. Pseudok´od je zn´azornˇen v algoritmu 18. Aby proˇrez´av´an´ı algoritmu bylo co moˇzn´a nejefektivnˇejˇs´ı, jsou u prohled´av´an´ı preferov´any v´ ybˇery akc´ı j´ızdy metrem, ˇc´ımˇz se rychle nalezne (prohled´av´an´ı do hloubky s backtrakingem) dostateˇcnˇe kvalitn´ı ˇreˇsen´ı. Pot´e se jiˇz horˇs´ı ˇreˇsen´ı zahazuj´ı neboli proˇrez´avaj´ı dˇr´ıve a zbyteˇcnˇe se d´ale nerozv´ıjej´ı, coˇz urychluje v´ ypoˇcet. Metoda vˇetv´ı a mez´ı byla zvolena i z d˚ uvodu, ˇze nebyla nalezena pˇr´ıpustn´a heuristika pro optim´aln´ı algoritmus A*. Metoda s prohled´av´an´ım do hloubky a preferov´an´ım j´ızdy
1
Pˇri spr´avn´em poˇrad´ı akc´ı.
´ ´ ´I TRASY 5.1. PLANOV AN metrem je dostateˇcnˇe efektivn´ı. Algorithm 1: Metoda vˇetv´ı a mez´ı s prohled´av´an´ım do hloubky Data: Prostˇred´ı ε ≡ ⟨S, A, T, Pkontroly , t, cj , cp ⟩;Vlastnosti α ≡ ⟨s0 , g, k, t0 , d⟩ Result: Pl´an π ≡ ⟨trasa, ts , j⟩; 1
Doln´ı mez ←− ∞;
2
foreach ts ::=
3
t0 : d do
foreach j ∈ {true, f alse} do
4
Otevˇren´a mnoˇzina←− s0 ;
5
Uzavˇren´a mnoˇzina←− ∅;
6
while Open ̸= ∅ do
7
P ≡ ⟨trasa, ts , j⟩ ←−Prvn´ı z otevˇren´e mnoˇziny();
8
Otevˇren´a mnoˇzina←−Otevˇren´a mnoˇzina \ P ;
9 10
Uzavˇren´a mnoˇzina←−Uzavˇren´a mnoˇzina ∪ P ; if (g ∈ trasa) ∧ (Doln´ı mez > Ohodnot’(P )) then
11
π ←− P ;
12
Doln´ı mez←−Ohodnot’(P );
13
else
14
P ′ ←−Rozˇsiˇr pl´an(P );
15
P ′ ←− P ′ \ Uzavˇren´a mnoˇzina;
16
foreach p ∈ P ′ do
17 18
if Doln´ı Mez > Ohodnot’(p) then Pˇridej na zaˇc´atek otevˇren´e mnoˇziny(p);
21
ˇ AN ´ ´I PASAZ ˇERA ´ KAPITOLA 5. UVAZOV
22
5.1.2
Popis pl´ anov´ an´ı
Pasaˇz´eˇri si kaˇzd´ y den pˇred moˇzn´ ym zaˇca´tkem sv´e cesty zaˇc´ınaj´ı pl´anovat sv´e trasy. Informace o prostˇred´ı a jejich vlastnostech, kter´e maj´ı k dispozici, se st´avaj´ı vstupem pro pl´anovac´ı algoritmus. Vstupy jsou tedy definovan´e v sekci 4.1 a 4.2 kapitoly 4 tj. vlastn´ı znaj´ı kapit´al, odkud, kam a kdy jet, z prostˇred´ı znaj´ı cenu l´ıstku, velikost pokuty, dopravn´ı s´ıt’ a pravdˇepodobnost kontroly v n´ı. 1. Pl´anovat pasaˇz´er zaˇc´ın´a ze stavu pˇred poˇca´teˇcn´ı stanic´ı v ˇcase startu. 2. Ze zn´am´eho prostˇred´ı zjist´ı, jak´e akce m˚ uˇze z t´eto pozice prov´est. Tato ˇc´ ast je vˇetv´ıc´ım m´ıstem. Velikost vˇetv´ıc´ıho faktoru a poˇrad´ı vyb´ıran´ych akc´ı urˇcuje sloˇzitost v´ypoˇctu, kter´a m˚ uˇze b´yt aˇz O(bd ), kde b je vˇetv´ıc´ı faktor a d hloubka prohled´ avan´eho stromu. 3. Zvol´ı si (dalˇs´ı) jednu akci, kterou si pˇrid´a do prozat´ımn´ıho pl´anu. K ostatn´ım akc´ım se vr´at´ı pozdˇeji, protoˇze prohled´av´ a do hloubky. Pokud jiˇz proˇsel vˇsechny akce, odebere z pl´anu naposled pˇridanou akci a pokraˇcuje bodem 3. Pokud vyzkouˇsel vˇsechny moˇzn´e akce a prozat´ımn´ı pl´an neobsahuje ˇz´adnou akci → prohled´av´an´ı konˇc´ı a vrac´ı nejlepˇs´ı nalezen´ y pl´an. Tento pl´an m˚ uˇze b´yt pr´azdn´y pouze v pˇr´ıpadˇe, ˇze neexistuje. 4. Ohodnot´ı prozat´ımn´ı pl´an souˇctem vˇsech cen akc´ı v pl´anu viz. kapitola 4 sekce 4.4. 5. Zkontroluje, zda cena pl´anu nen´ı horˇs´ı neˇz cena doposud nejlepˇs´ıho nalezen´eho pl´anu. Pokud je horˇs´ı, akci z pl´anu odstran´ı a pokraˇcuje bodem 3. V tomto m´ıstˇe nast´av´ a moˇzn´e proˇrez´an´ı prohled´ avan´eho stromu. Proˇr´ızne celou vˇetev pod pl´anem, kter´y jiˇz m´a horˇs´ı cenu, neˇz doposud nalezen´y. 6. Pokud m´a nov´ y pl´an lepˇs´ı cenu, tak ji pasaˇz´er pomyslnˇe provede a posune se do dalˇs´ıho stavu (ve stanici nebo pˇred stanic´ı). 7. Zkontroluje, zda je dosaˇzen´ y stav c´ılov´ y. Pokud ano, uloˇz´ı ho jako doposud nejlepˇs´ı nalezen´ y, odstranˇen´ı posledn´ı akce z pl´anu a pokraˇcuje bodem 3. Pokud nen´ı c´ılov´ ym m´ıstem, pokraˇcuje dalˇs´ım rozv´ıjen´ım pl´anu tedy bodem 2. Tento postup pasaˇz´er opakuje pro r˚ uzn´e ˇcasy v´ yjezd˚ u poˇc´ınaj´ıc´ı moˇzn´ ym startem. Pro kaˇzd´ y ˇcas v´ yjezdu zkouˇs´ı nal´ezt pl´an s j´ızdenkou i bez j´ızdenky. Postupem v´ yˇse proch´az´ı cel´ y stavov´ y prostor tam, kde jsou pl´any s lepˇs´ı cenou neˇz je doposud nejlepˇs´ı nalezen´ y pl´an.
ˇ ´I SE O PREPRAVN ˇ ´I KONTROLE 5.2. UCEN
23
V´ ystup tvoˇr´ı sestaven´ y ”cestovn´ı pl´an” tj. posloupnost akc´ı, ˇcas v´ yjezdu, a volbu j´ızdenky viz. sekce 4.3 kapitoly 4.
5.1.3
Korektnost algoritmu
Algoritmus postupn´ ym rozˇsiˇrov´an´ım pl´an˚ u tvoˇr´ı nov´e pl´any. Vr´at´ı tedy pl´an vedouc´ı do c´ılov´e stanice nebo ˇz´adn´ y, pokud cesta neexistuje. Zvolen´ y algoritmus je proto korektn´ı.
5.1.4
´ Uplnost algoritmu
Protoˇze se jedn´a o koneˇcn´ y stavov´ y diagram s koneˇcnou mnoˇzinou stav˚ u a akc´ı, a protoˇze algoritmus nekonˇc´ı s prvn´ım nalezen´ ym ˇreˇsen´ım, ale hled´a dalˇs´ı moˇzn´a ˇreˇsen´ı, je algoritmus tak´e u ´pln´ y. Je u ´pln´ y i vzhledem k jeho proˇrez´av´an´ı. Algoritmus proˇrez´av´a pouze takov´a ˇreˇsen´ı, kter´e jiˇz nemohou dos´ahnout ˇreˇsen´ı s lepˇs´ı cenou, neˇz je doposud nalezen´e.
5.2
Uˇ cen´ı se o pˇ repravn´ı kontrole
Pl´anovac´ı algoritmus hled´a ˇreˇsen´ı s nejniˇzˇs´ı cenou v aktu´aln´ım prostˇred´ı, kter´e m´a k dispozici. Model vˇsak m˚ uˇze obsahovat zkreslen´e informace a v´ ysledky pl´anov´an´ı pak mohou b´ yt chybn´e. Pasaˇz´er proto potˇrebuje schopnost uˇcen´ı se spr´avn´emu modelu a to jak pˇri informac´ıch chybˇej´ıc´ıch, tak zastaral´ ych vlivem dynamiky syst´emu. Pro adaptaci na (ne)deterministicky vyv´ıjej´ıc´ı se prostˇred´ı nen´ı uvaˇzov´ano pln´e uˇcen´ı s posilov´an´ım, protoˇze v mˇen chyb´ı f´aze explorace. Pasaˇz´er vyuˇz´ıv´a exploitaci v podobˇe prost´eho statistick´eho uˇcen´ı, kter´e slouˇz´ı k pozn´av´an´ı a pˇredv´ıd´an´ı jev˚ u v re´aln´ ych podm´ınk´ach. Do uˇcen´ı jsou nav´ıc zaˇclenˇeny apriorn´ı pravdˇepodobnosti metodou fiktivn´ıch dat, kter´a je pˇrid´ana z d˚ uvodu pasaˇz´erovy poˇc´ateˇcn´ı informace o prostˇred´ı. pk (a, t) =
n1 , n1 + n0
(5.1)
kde n1 je poˇcet kolikr´at byl pasaˇz´er pˇri prov´adˇen´ı akce a v ˇcase t kontrolov´an a n0 poˇcet kdy ne. Vztah 5.1 ud´av´a v´ ypoˇcet pravdˇepodobnosti kontroly pˇri prov´adˇen´ı akce a v ˇcase t. Pˇrid´an´ı apriorn´ı pravdˇepodobnosti spoˇc´ıv´a pouze v poˇc´ateˇcn´ım nastaven´ı n1 a n0 podle
ˇ AN ´ ´I PASAZ ˇERA ´ KAPITOLA 5. UVAZOV
24
z´ıskan´ ych fiktivn´ıch dat viz. rovnice 5.2 a 5.3. n1 =
T ∑
ft
(5.2)
(1 − ft ),
(5.3)
Metoda fiktivn´ıch dat je kombinac´ı ”expertn´ı”
informace, u kter´e se apriorn´ı
t=1
n0 =
T ∑ t=1
kde f ∈ {0, 1} je vektor fiktivn´ıch dat. pravdˇepodobnost zad´av´a na z´akladˇe zkuˇsenost´ı experta, a apriorn´ıch dat namˇeˇren´ ych do zaˇca´tku odhadov´an´ı. Metoda spoˇc´ıv´a v sestavov´an´ı fiktivn´ıch pokus˚ u expertem tak, aby data vych´azej´ıc´ı z pokus˚ u odpov´ıdala pˇredstav´am experta (Nagy, I., Nedoma, P., Kratky, M., Pavelkova, L. a Ettler, P., 2002). Horizont, kter´ y ud´av´a poˇcet fiktivn´ıch dat, ud´av´a ve statistice d˚ uvˇeryhodnost experta tj. nakolik pasaˇz´er vˇeˇr´ı z´ıskan´ ym dat˚ um. N´asledn´a u ´prava pravdˇepodobnosti pk (a, t) se odhaduje jako aritmetick´ y pr˚ umˇer. Pasaˇz´er upravuje pravdˇepodobnosti kontroly, na z´akladˇe pozorov´an´ı, vˇzdy po dokonˇcen´ı cel´eho pl´anu. Pasaˇz´er m˚ uˇze upravovat pouze ty pravdˇepodobnosti, o kter´ ych z´ıskal novou informaci tj. pro vˇsechny akce z pl´anu a ∈ π v ˇcase t. Kr´atk´ y horizont odpov´ıd´a slabˇs´ı apriorn´ı pravdˇepodobnosti tzn. ˇze se uˇcen´ y model bude rychleji adaptovat na zmˇeny prostˇred´ı neboli rychleji uprav´ı chybn´e informace o prostˇred´ı. S kr´atk´ ym horizontem nemaj´ı pˇr´ıliˇs velkou v´ahu informace od experta. Pˇri nepravdiv´e expertn´ı informaci, bude zpoˇc´atku budouc´ı odhadovan´a hodnota pravdˇepodobnˇe chybn´a, ale pozorov´an´ı bude postupnˇe korigovat ˇspatnou apriorn´ı znalost.
5.3
Shrnut´ı
V t´eto kapitole bylo k modelu pasaˇz´era pˇrid´ano i jeho rozhodov´an´ı v podobˇe hled´an´ı pl´anu s nejniˇzˇs´ı cenou a zp˚ usob pasaˇz´erovy adaptace na prostˇred´ı. Lze tedy pˇrej´ıt od teoretick´e ˇc´asti k praktick´e.
Kapitola 6 Implementace Model pasaˇz´era popsan´ y v kapitole 4, pasaˇz´erovo pl´anov´an´ı a uˇcen´ı popsan´e v kapitole 5, bylo implementov´ano jako rozˇs´ıˇren´ı st´avaj´ıc´ıho frameworku FareEvasion. Implementace agentn´ıho modelu pˇrepravn´ı kontroly je tedy n´astavbou na existuj´ıci framework FareEvasion, kter´ y tvoˇr´ı n´astavbu agentn´ı simulaˇcn´ı platformy AgentPolis slouˇz´ıc´ı pro simulaci, verifikaci a experimentov´an´ı agentn´ıch model˚ u. FareEvasion, AgentPolis i platforma A-lite, na kter´e je AgentPolis postaven, jsou naps´any v programovac´ım jazyku Java. Implementace rozˇs´ıˇren´ı tˇechto simulaˇcn´ıch platforem zachov´av´a tento objektovˇeorientovan´ y jazyk. Uk´azka struktury rozˇs´ıˇren´ı a z´avislost´ı mezi platformami je zobrazena na obr´azku 6.1. V t´eto ˇc´asti zpr´avy je bl´ıˇze pops´ana implementace rozˇs´ıˇren´ı platformy agentn´ım modelem pˇrepravn´ı kontroly. Jsou zde popsan´e tˇr´ıdy pouˇzit´e v implementaci.
6.1
Tˇ r´ıdy pro model pasaˇ z´ era
Vytvoˇren´e tˇr´ıdy pro chov´an´ı dˇr´ıve popsan´eho pasaˇz´era byly inspirov´any z tˇr´ıd, kter´e umoˇzn ˇuj´ı pohyb z´akladn´ıho pasaˇz´era v projektu FareEvasion.
6.1.1
Vytvoˇ ren´ı pasaˇ z´ era
Vytvoˇren´ı pasaˇz´era je implementov´ano ve tˇr´ıdˇe PassengerInit a spoˇc´ıv´a v nastaven´ı poˇc´ateˇcn´ıch vlastnost´ı pasaˇz´era implementovan´e ve tˇr´ıdˇe PassengerStatus. D´ale se zde inicializuje ˇzivotn´ı cyklus pasaˇz´era a vytvoˇr´ı se pl´anovaˇc trasy. 25
26
KAPITOLA 6. IMPLEMENTACE
Obr´azek 6.1: Z´avislosti mezi platformami
6.1.2
Vlastnosti pasaˇ z´ era
Vlastnosti pasaˇz´era jsou um´ıstˇeny ve tˇr´ıdˇe PassengerStatus a obsahuj´ı informace shrnut´e v definici 4.2 v kapitole 4. D´ale jsou ve vlastnostech pasaˇz´era uloˇzeny informace o prostˇred´ı z definice prostˇred´ı 4.1. Tyto informace maj´ı sice zpoˇca´tku nastaven´e vˇsichni pasaˇz´eˇri stejnˇe, ale d´ale si je kaˇzd´ y upravuje vlastn´ım zp˚ usobem podle sv´ ych zkuˇsenost´ı z j´ızd. Z d˚ uvodu ˇcasov´e spojitosti, existuje nekoneˇcnˇe r˚ uzn´ ych stav˚ u prostˇred´ı. Pro jejich koneˇcnost byla diskretizov´ana ˇcasov´a osa mˇen´ıc´ıho se prostˇred´ı tj. vˇsechny moˇzn´e stavy z urˇcit´eho ˇcasov´eho rozmez´ı se nahrad´ı jedn´ım reprezentantem. T´ım je zaruˇcena koneˇcnost stav˚ u prostˇred´ı. Pasaˇz´er napˇr. ˇcasov´ yu ´sek od 12:30 do 12:40 vn´ım´a jako jeden stav, tedy prostˇred´ı se v tomto ˇcasov´em intervalu pro nˇeho nezmˇen´ı. Jelikoˇz si v t´eto tˇr´ıdˇe pasaˇz´er udrˇzuje vlastn´ı pˇredstavu o prostˇred´ı, bylo tak´e nutn´e implementovat jeho uˇcen´ı. S uˇcen´ım, kter´e bylo definov´ano v kapitole 5, je spojen´a konstanta horizontu fiktivn´ıch j´ızd, urˇcuj´ıc´ı vˇerohodnost poˇc´ateˇcn´ıch informac´ı. Uˇcen´ı neboli adaptace na aktu´aln´ı prostˇred´ı prob´ıh´a u kaˇzd´eho navˇst´ıven´eho stavu prostˇred´ı1 . Tˇr´ıda implementuje i model pasaˇz´era, kter´eho se pl´anovaˇc, volan´ y z ˇzivotn´ıho cyklu 1
Stav prostˇred´ı je ch´ ap´ an jako stanice ˇci u ´sek v dan´em ˇcasov´em rozmez´ı
ˇ ´IDY PRO MODEL PASAZ ˇERA ´ 6.1. TR
27
pasaˇz´era dotazuje, na uˇzitkov´e hodnoty aktu´aln´ıho prohled´avan´eho stavu.
6.1.3
Pl´ an pasaˇ z´ era
V implementaci tˇr´ıdy PassengerPlan je uchov´av´an jiˇz cel´ y pl´an pasaˇz´era.
6.1.4
ˇ Zivotn´ ı cyklus pasaˇ z´ era
ˇ Zivotn´ ı cyklus pasaˇz´era zaˇc´ın´a ”narozen´ım agenta” tedy metodou born, ve kter´e se pasaˇz´erovi nastav´ı kaˇzdodenn´ı cesta. Pokraˇcov´an´ı v denn´ım cyklu pasaˇz´era je pl´anov´an´ı t´eto trasy zobrazen´e na stavov´em diagramu v obr´azku 6.2.
Obr´azek 6.2: Stavov´ y diagram ˇzivotn´ıho cyklu pasaˇz´era
Po pl´anov´an´ı pasaˇz´er ˇcek´a na ˇcas, kdy m˚ uˇze zaˇc´ıt plnit sv˚ uj pl´an, kter´ y pˇri dosaˇzen´ı spr´avn´eho ˇcasu zaˇcne plnit. Na z´avˇer jednoho cyklu pasaˇz´er vyhodnot´ı svou cestu t´ım, ˇze si uprav´ı informace o prostˇred´ı (uˇc´ı se). Implementace prov´adˇen´ı pasaˇz´erova pl´anu je shodn´a se stavov´ ym diagramem zobrazen´em v obr´azku 6.3.
28
KAPITOLA 6. IMPLEMENTACE
Obr´azek 6.3: Stavov´ y diagram prov´adˇen´ı pl´anu
6.1.5
Pl´ anovaˇ c
Pl´anovaˇc cesty vyuˇz´ıv´a metody Vˇetv´ı a mez´ı popsan´e v kapitole 5. Pro pl´anov´an´ı pouˇz´ıv´a oznaˇcen´ı cel´e linky ROUTE ID.
6.1.6
Stav
Stav, ve kter´em se pasaˇz´er nach´az´ı, je implementov´an ve tˇr´ıdˇe PlannerNode. Tento stav uchov´av´a informace o vˇsech pˇredchoz´ıch stavech a akc´ıch, kter´e vedly k dosaˇzen´ı tohoto stavu. Na kaˇzd´ y stav se nahl´ıˇz´ı jako na ˇc´asteˇcn´e nebo u ´pln´e ˇreˇsen´ı, tud´ıˇz je potˇreba zn´at ohodnocen´ı nebo-li pasaˇz´erovu cenu dosaˇzen´ı tohoto stavu. Stavu vyuˇz´ıv´a pl´anovaˇc, kter´ y prostˇrednictv´ım tˇr´ıdy PassengerStatus, tento stav ohodnot´ı uˇzitkovou funkc´ı pl´anuj´ıc´ıho pasaˇz´era.
ˇ ´IDY PRO MODEL INSPEKTORA 6.2. TR
6.1.7
29
Ukl´ ad´ an´ı informac´ı
Pro vyhodnocen´ı v´ ysledk˚ u kaˇzd´e simulace se v jej´ım pr˚ ubˇehu ukl´adaj´ı informace o pasaˇz´erovi (typ pasaˇz´era), jeho pl´anu (trasa, j´ızdenka) a vyhodnocen´ı pl´anu po jeho dokonˇcen´ı (kontrola, pokutov´an´ı). Na konci kaˇzd´e simulace se vˇsechny informace uloˇz´ı do souboru s hodnotami oddˇelen´ ymi ˇca´rkou (csv soubor˚ u)2 .
6.2
Tˇ r´ıdy pro model inspektora
Chov´an´ı inspektora v projektu FareEvasion se v´ yraznˇe neliˇs´ı od tohoto inspektora, proto byly pˇrevzaty z nˇeho. Upraven byl zejm´ena ˇzivotn´ı cyklus inspektora, ve kter´em se rozliˇsuje v´ ybˇer strategie kontroly.
6.2.1
Vytvoˇ ren´ı inspektora
Vytvoˇren´ı a nastaven´ı vlastnost´ı inspektora prob´ıh´a ve tˇr´ıdˇe inspektorInit. Vlastnost´ı inspektora je poˇcet pasaˇz´er˚ u, kter´e inspektor zkontroluje bˇehem jedn´e minuty a pomˇer ud´avaj´ıc´ı kolik procent pasaˇz´er˚ u m˚ uˇze inspektor maxim´alnˇe zkontrolovat na jednom u ´seku tj. mezi sousedn´ımi stanicemi. Pˇri vytv´aˇren´ı inspektora se tak´e spust´ı jeho ˇzivotn´ı cyklus.
6.2.2
ˇ Zivotn´ ı cyklus inspektora
Do modelu pˇrepravn´ı kontroly byly implementov´any dva typy ˇzivotn´ıch cykl˚ u inspektora. Tyto ˇzivotn´ı cykly se liˇs´ı pouze ve v´ ybˇeru strategie kontroly pro n´asleduj´ıc´ı den. V jednom ˇzivotn´ım cyklu je v´ ybˇer strategie kontroly na dalˇs´ı den shodn´ y se strategi´ı v pˇredchoz´ım dni, jak je uk´az´ano ve stavov´em diagramu 6.43 . Takov´e chov´an´ı inspektor˚ u je pˇredv´ıdateln´e, proto prostˇred´ı pouze s t´ımto typem inspektor˚ u je deterministick´e. Kaˇzd´ y inspektor bude m´ıt kaˇzd´ y den stejnou strategii kontroly tzn. bude kontrolovat na shodn´ ych m´ıstech, ve stejn´ y ˇcas, kaˇzd´ y den. 2 3
CSV: Comma-separated values - hodnoty oddˇelen´e ˇc´arkami. Obr´azek byl pˇrekreslen ze zpr´avy (Moler, Z., n.d.).
30
KAPITOLA 6. IMPLEMENTACE
Obr´azek 6.4: Diagram deterministick´eho ˇzivotn´ıho cyklu inspektora
Druh´ ym ˇzivotn´ım cyklem se vyb´ır´a nov´a strategie kontroly kaˇzd´ y den. Stavov´ y diagram 6.54 ukazuje chov´an´ı inspektora s t´ımto druhem ˇzivotn´ıho cyklu. Prostˇred´ı, v nˇemˇz se vyskytuj´ı tito inspektoˇri, je nedeterministick´e. V´ ybˇer strategie kontroly prob´ıh´a na zaˇca´tku kaˇzd´eho dne v´ yˇse uveden´ ym zp˚ usobem. Kaˇzd´ y inspektor m´a k dispozici seznam moˇzn´ ych strategi´ı, kde kaˇzd´a strategie kontroly obsahuje pravdˇepodobnost, s kterou se bude vykon´avat. Pravdˇepodobnosti vˇsech strategi´ı tvoˇr´ı hustotu pravdˇepodobnosti, ze kter´e se vyb´ıraj´ı jednotliv´e strategie pro inspektora. Kaˇzd´a strategie inspektora se m˚ uˇze skl´adat z jedn´e nebo v´ıce akc´ı, kter´e inspektor mus´ı bˇehem dne prov´est napˇr. kontrola ve stanici, kontrola na nˇejak´em u ´seku a kontrola v jin´e stanici.
4
Obr´azek byl inspirov´ an zpr´avou (Moler, Z., n.d.).
6.3. SHRNUT´I
31
Obr´azek 6.5: Diagram nedeterministick´eho ˇzivotn´ıho cyklu inspektora
6.3
Shrnut´ı
K ovˇeˇren´ı implementace byly provedeny a vyhodnoceny experimenty, kter´e jsou rozeps´any v n´asleduj´ıc´ı kapitole. K jejich vyhodnocen´ı bylo potˇrebn´e pr˚ ubˇeˇznˇe ukl´adat informace uveden´e v´ yˇse. K vyhodnocen´ı v´ ysledk˚ u simulace (vykreslen´ı graf˚ u, . . . ), byl pouˇzit v´ ypoˇcetn´ı n´astroj MATLAB. Kv˚ uli prov´adˇen´ı v´ıce experiment˚ u, byl opˇet v programovac´ım jazyku java implementov´an spouˇstˇec´ı program, kter´ y postupnˇe prov´adˇel dan´e experimenty. Z d˚ uvodu rychl´eho vyhodnocov´an´ı v´ ysledk˚ u bylo do nˇeho implementov´ano pˇredzpracov´an´ı dat pro MATLAB5 . Projekt FareEvasion je implementovan´e grafick´e rozhran´ı, kter´e je uk´az´ano na obr´azku 6.6(a) pro jednu linku a na 6.6(b) pro celou s´ıt’, zobrazuj´ıc´ı dr´ahy linek dopravn´ı s´ıtˇe, moˇznosti pˇeˇs´ıho pˇrestupu, pohybuj´ıc´ı se dopravn´ı soupravy, rozm´ıstˇen´ı inspektor˚ u a pohyb dopravn´ı soupravy s pasaˇz´ery. 5
Pˇreveden´ı MATLAB k´odu na jar soubor, kter´ y lze spustit javou. Tento jar soubor vol´a MATLAB
knihovny, bez kter´ ych se neobejde.
32
KAPITOLA 6. IMPLEMENTACE
(a) Simulace jedn´e linky
(b) Simulace cel´e s´ıtˇe
Obr´azek 6.6: Zobrazen´ı grafick´eho prostˇred´ı simulaˇcn´ı platformy
Kapitola 7 Experimenty V pˇredchoz´ıch kapitol´ach bylo definov´ano prostˇred´ı, model pasaˇz´era a zp˚ usob pasaˇz´erovy adaptace na okoln´ı prostˇred´ı. Adaptac´ı je myˇsleno uˇcen´ı se nov´emu (aktu´aln´ıho) prostˇred´ı nebo nauˇcen´ı se na zmˇeny v prostˇred´ı. V t´eto kapitole je provedena a vyhodnocena sada experiment˚ u ovˇeˇruj´ıc´ı pˇredpoklady1 chov´an´ı vytvoˇren´eho modelu. Experimenty jsou rozdˇeleny do dvou ˇc´ast´ı. Prvn´ı ˇc´ast tvoˇr´ı z´akladn´ı experimenty pro simulaci chov´an´ı jednotliv´ ych typ˚ u pasaˇz´er˚ u a druhou ˇc´ast tvoˇr´ı experimenty s populac´ı pasaˇz´er˚ u a inspektor˚ u. V druh´e ˇc´asti jsou rozdˇeleny podle velikosti prostˇred´ı (jedna linka, cel´a s´ıt’).
7.1 7.1.1
Nastaven´ı pro experimenty Konfigurace prostˇ red´ı
Vˇetˇsina parametr˚ u, potˇrebn´ ych ke konfiguraci prostˇred´ı, je pops´ana v kapitole 4 v sekci ”reprezentace prostˇred´ı pasaˇz´era”. Jsou jimi vˇsechny stanice dopravn´ı s´ıtˇe reprezentovan´e stavy a vˇsechny zp˚ usoby pˇrepravy mezi nimi v podobˇe akc´ı. Jedn´ım z parametr˚ u, kter´e nejsou viditelnˇe zanesen´e v reprezentaci prostˇred´ı je poˇcet linek. Prostˇred´ı, kter´e je pouˇz´ıvan´e frameworkem FareEvasion (v kapitole 6) je podobn´e dopravn´ı s´ıti americk´eho mˇesta Los Angeles viz. obr´azek 7.12 . V prostˇred´ı, kter´e je tvoˇreno geograficky rozm´ıstˇen´ ymi stanicemi a trasami mezi nimi, je umoˇznˇeno cestov´an´ı pouze jedn´ım dopravn´ım prostˇredkem, kter´ ym je metro. Samozˇrejmˇe je umoˇznˇena i pˇreprava mezi stani1 2
Intuitivn´ı chov´ an´ı. http://www.metro.net/riding metro/maps/images/rail map.pdf
33
34
KAPITOLA 7. EXPERIMENTY
cemi ch˚ uz´ı. Z v´ ypoˇcetn´ı n´aroˇcnosti v simulaci byly uvaˇzov´any jen d˚ uleˇzit´e pˇeˇs´ı pˇrechody mezi stanicemi, kter´ ymi jsou kr´atk´e pˇrechody pˇrev´aˇznˇe v okol´ı pˇrestupn´ıch stanic a pˇrestupy mezi sousedn´ımi stanicemi3 . Ned˚ uleˇzit´e pˇrechody, kter´e zpomaluj´ı v´ ypoˇcet, jsou pˇr´ıliˇs dlouh´e a pro pasaˇz´ery nev´ yhodn´e. Konfiguraci prostˇred´ı tedy urˇcuj´ı tyto parametry: • Seznam vˇsech stanic, • seznam vˇsech pˇr´ım´ ych cest mezi stanicemi, • zp˚ usoby pˇrepravy, • kapacita dopravn´ıho prostˇredku, • rychlost cestov´an´ı v dopravn´ım prostˇredku, • rychlost ch˚ uze, • diskretizace ˇcasu, • cena j´ızdenky a • velikost pokuty. Kde pro experimenty jsou zafixov´any parametry uveden´e v tabulce 7.1. Tabulka 7.1: Tabulka fixn´ıch parametr˚ u prostˇred´ı
Parametry
Hodnota
Parametry kontroly Kontrol za minutu [-]
2
Max. kontrol na u ´seku [%]
0.1
Parametry pˇrepravy Cena j´ızdenky [$]
4
Velikost pokuty [$]
50
Kapacita prostˇredku [os.] Rychlost prostˇredku Rychlost ch˚ uze [ km ] h 3
[ km h
]
50 40 4.3
Pasaˇz´er je schopen kombinac´ı tˇechto pˇrechod˚ u uj´ıt v´ıce stanic. Tato varianta m˚ uˇze b´ yt m´enˇe v´ yhodn´a,
neˇz by byla ta skuteˇcn´ a.
7.1. NASTAVEN´I PRO EXPERIMENTY
35
Obr´azek 7.1: Sch´ema dopravn´ı s´ıtˇe
7.1.2
Nastaven´ı sc´ en´ aˇ re
Sc´en´aˇr obsahuje nˇekolik parametr˚ u, kter´e je potˇreba nastavit pro pˇribl´ıˇzen´ı k re´aln´ ym v´ ysledk˚ u tˇechto experiment˚ u. Volby hodnot tˇechto parametr˚ u tedy urˇcuj´ı, jak se bu-
36
KAPITOLA 7. EXPERIMENTY
dou vyv´ıjet mˇeˇren´e metriky, tedy poˇzadovan´e (vyhodnocovan´e) v´ ystupy z experimentu. Nˇekter´e z parametr˚ u sc´en´aˇre z˚ ustanou, at’ uˇz z v´ ypoˇcetn´ıch nebo jin´ ych d˚ uvod˚ u, nemˇenn´e. Seznam parametr˚ u sc´en´aˇre lze rozdˇelit do dvou skupin. Jedna skupina je tvoˇrena nastaven´ım inspektor˚ u tedy parametry: • Poˇcet inspektor˚ u, • rozvrh (strategie) inspektor˚ u, • rozdˇelen´ı pravdˇepodobnosti v´ ybˇeru strategie, • poˇcet kontrol inspektorem za minutu a • kolik procent pasaˇz´er˚ u maxim´alnˇe zkontroluje inspektor na jednom u ´seku, Prvn´ı dva parametry se liˇs´ı experiment od experimentu. Sada strategi´ı inspektor˚ u a jejich pravdˇepodobnostn´ı rozdˇelen´ı byla z´ısk´ana od dopravn´ıho podniku Los Angeles d´ale v pˇriloˇzen´em CD. Pokud nen´ı v experimentech uvedeno jinak, je vyuˇz´ıv´ano t´eto sady strategi´ı s v´ ybˇerem podle jejich rozdˇelen´ı pravdˇepodobnosti. Druh´a ˇca´st tvoˇr´ı vlastnosti druh´eho typu agenta umoˇzn ˇuj´ıc´ı simulaci, kter´ ym je pasaˇz´er. Parametry jsou popsan´e v kapitole 4 v definici 4.2. • poˇcet cestuj´ıc´ıch pasaˇz´er˚ u, • sada cest pasaˇz´er˚ u: – Poˇc´ateˇcn´ı stanice, – c´ılov´a stanice, – moˇzn´ y ˇcas v´ yjezdu, – nutn´ y ˇcas dojezdu do c´ılov´e stanice, • pravdˇepodobnostn´ı rozdˇelen´ı v´ ybˇeru cesty, • parametry uˇzitku pasaˇz´era: – Kapit´al pasaˇz´era, – penalizace za pˇrekroˇcen´ı ˇcasu nutn´eho dojezdu, – parametry z prostˇred´ı (cena j´ızdenky a velikost pokuty), • parametry urˇcuj´ıc´ı pasaˇz´er˚ uv typ:
7.1. NASTAVEN´I PRO EXPERIMENTY
37
– preferuj´ıc´ı ˇcasov´ y uˇzitek, – preferuj´ıc´ı ekonomick´ y uˇzitek, – preferuj´ıc´ı psychologick´ y uˇzitek, • apriorn´ı znalost prostˇred´ı a • horizont (pro fiktivn´ı data od experta). Cesty pasaˇz´er˚ u jsou stejnˇe jako v pˇredchoz´ı sadˇe (jako strategie inspektor˚ u) z´ısk´any od dopravn´ı spoleˇcnosti Los Angeles. Denn´ı rozpoˇ cet pasaˇ z´ era Z d˚ uvodu podobnosti dopravn´ı s´ıtˇe s mˇestem Los Angeles byly vygenerov´any pasaˇz´er˚ um denn´ı rozpoˇcty dle ekonomick´eho rozloˇzen´ı v LA v roce 2009, tabulka A.1. Nastaven´ı parametr˚ u urˇ cuj´ıc´ı typ pasaˇ z´ era Podle pozorov´an´ı jednotliv´ ych sloˇzek cenov´e funkce pasaˇz´era, kter´a p˚ usob´ı na chov´an´ı pasaˇz´era, byly nastaveny ˇr´adovˇe vlivy jednotliv´ ych sloˇzek. Pasaˇz´eˇri se tedy dˇel´ı do skupin podle toho jak´a sloˇzka u nich pˇrevaˇzuje. Fixn´ı nastaven´ı parametr˚ u, kter´a jsou zobrazena v tabulce 7.2, z˚ ust´avaj´ı ve vˇsech experimentech shodn´a. Tabulka 7.2: Tabulka fixn´ıch parametr˚ u normalizaˇcn´ıch konstant pasaˇz´er˚ u
Typ pasaˇ z´ era
ˇ Casov´ y uˇ zitek Ekonom. uˇ zitek Psycholog. uˇ zitek
ˇ Casov´ y [-]
5 · 106
100
106
Ekonomick´ y [-]
2 · 107
5 · 10−1
106
Psychologick´ y [-]
2 · 107
100
5 · 102
7.1.3
Metriky k vyhodnocen´ı
Pro vyhodnocen´ı experiment˚ u je potˇreba bˇehem experimentu ukl´adat r˚ uzn´e parametry simulace pˇrepravy, ze kter´ ych lze n´aslednˇe tvoˇrit z´avˇery o chov´an´ı syst´emu. • U modelu pˇrepravn´ı kontroly je oˇcek´avanou a sledovanou statistikou v´ yvoj poˇctu neplat´ıc´ıch pasaˇz´er˚ u. Dopravn´ı podnik se ˇcasto snaˇz´ı sledovat, jakou ztr´atu jim neplat´ıc´ı pasaˇz´eˇri zp˚ usobili nebo pˇr´ıpadnˇe zisk z udˇelen´ ych pokut v pˇr´ıpadˇe nasazen´ı inspektor˚ u do dopravn´ı s´ıtˇe.
38
KAPITOLA 7. EXPERIMENTY • Pozorovat dynamiku syst´emu tj. za jak dlouhou dobu se pasaˇz´eˇri nauˇc´ı rozm´ıstˇen´ı inspektor˚ u. • Sledov´an´ı velikosti u ´niku zisk˚ u (v podobˇe zakoupen´ ych l´ıstk˚ u), poˇctu kontrol a udˇelen´ ych pokut v ust´alen´em syt´emu v z´avislosti na poˇctu inspektor˚ u v nˇem se pohybuj´ıc´ıch. Dopravn´ı spoleˇcnosti toto pozorov´an´ı umoˇzn ˇuje maximalizovat zisk. • Sledovat, jak se pasaˇz´eˇri zachovaj´ı pˇri statick´e kontrole. Kolik pasaˇz´er˚ u proj´ıˇzd´ı pˇrestupn´ı stanic´ı pˇred a po pˇrid´an´ı inspektora.
7.1.4
Nastaven´ı simulace
Pro vˇetˇsinu experiment˚ u byla z v´ ypoˇcetn´ıch n´arok˚ u pouˇzita jen jedna linka. Bylo by zaj´ımav´e pozorovat, jak moc se pasaˇz´eˇri adaptuj´ı ve velk´e s´ıti, tedy s moˇznost´ı pˇrestup˚ u. Cel´a s´ıt’, kter´a je pro experimenty k dispozici nen´ı natolik komplexn´ı, aby se v n´ı uk´azaly moˇznosti obj´ıˇzdˇen´ı. Prostˇred´ı agentn´ı kontroly bylo simulov´ano ve vˇsech experimentech po dobu 200 dn´ı. V simulovan´em prostˇred´ı byla z v´ ypoˇcetn´ıch d˚ uvod˚ u zvolena velikost populace pro jednu linku rovna 300 a pro celou s´ıt’ 120. D´ale pro jednoduˇsˇs´ı vyhodnocov´an´ı m´a kaˇzd´ y pasaˇz´er v prostˇred´ı urˇcenou jedinou cestu, kterou bude dennˇe jezdit. Protoˇze stavov´ y prostor, kter´ y pokr´ yv´a simulovan´e prostˇred´ı pˇres cel´ y den, je velk´ y, bylo zavedeno zjednoduˇsen´ı v podobˇe simulov´an´ı cest zaˇc´ınaj´ıc´ıch mezi ˇsestou a osmou hodinou rann´ı. Pro simulaci bylo pouˇzito diskretizov´an´ı ˇcasu po 10 minutov´ ych intervalech.
7.2. EXPERIMENTY S JEDNOTLIVCI
7.2
39
Experimenty s jednotlivci
Tato ˇc´ast experiment˚ u se zab´ yv´a z´akladn´ım chov´an´ım jednotliv´ ych typ˚ u pasaˇz´er˚ u a jejich adaptac´ı na zmˇeny prostˇred´ı. Je evidentn´ı, ˇze pasaˇz´er zvol´ı sv˚ uj nejv´ yhodnˇejˇs´ı pl´an pˇri nastaven´ı spr´avn´ ych parametr˚ u prostˇred´ı tj. pravdˇepodobnosti chycen´ı z´ıskan´e od experta. Pasaˇz´er by v tomto pˇr´ıpadˇe zvolil jeden pl´an (nejv´ yhodnˇejˇs´ı), kter´eho by se drˇzel kaˇzd´ y den (po celou dobu simulace). Stejn´a situace, kdy pasaˇz´er sv˚ uj pl´an nemˇen´ı, je bez schopnosti pasaˇz´era uˇcen´ı se zmˇen´am. Pasaˇz´er v t´eto situaci v´ı, jak´ y je jeho nejv´ yhodnˇejˇs´ı pl´an a jiˇz nezkoum´a, jestli jsou jeho informace spr´avn´e. Protoˇze v modelu nen´ı f´aze explorace, pasaˇz´er nebude n´ahodnˇe zkouˇset, zda existuje nˇejak´a jin´a cesta, kter´a by pro nˇej byla v urˇcit´e situaci v´ yhodnˇejˇs´ı. Zmˇenu pl´anu mohou zp˚ usobit pouze expertem ˇspatnˇe nastaven´e poˇca´teˇcn´ı pravdˇepodobnosti nebo zmˇena v prostˇred´ı. Pasaˇz´er si zvol´ı aktu´alnˇe nejv´ yhodnˇejˇs´ı pl´an a vlivem uˇcen´ı muˇze doch´azet i k postupn´emu zhorˇsov´an´ı ceny pl´anu aˇz do doby, kdy se pro pasaˇz´era stane v´ yhodnˇejˇs´ı jin´ y. V t´eto ˇc´asti experiment˚ u jsou nastaveny dalˇs´ı fixn´ı parametry sc´en´aˇre a to pro poˇcet pasaˇz´er˚ u roven tˇrem, jedinou moˇznou cestu a to ze stanice oznaˇcen´e ˇc. 63 do stanice ˇc. 60, horizont expertn´ıch znalost´ı roven 100 a denn´ı kapit´al pasaˇz´er˚ u pro prvn´ı ze dvou experiment˚ u ˇcin´ı 25$ a druh´ y 75$. Kapit´al je takto stanoven z d˚ uvodu pozorov´an´ı chov´an´ı ´ pasaˇz´era s menˇs´ım/vˇetˇs´ım kapit´alem, neˇz je velikost pokuty. Udaje o experimentech, ve kter´ ych jsou rovnomˇernˇe rozdˇelen´e typy pasaˇz´er˚ u, shrnuje tabulka 7.3. Tabulka 7.3: Tabulka nastaven´ı experiment˚ u v prostˇred´ı s jednotlivci
7.2.1
Experiment
Cesta
Poˇ cet pasaˇ z´ er˚ u Horizont Kapit´ al
ˇc. 1
63 - 60
3
100
25
ˇc. 2
63 - 60
3
100
75
Chybn´ e poˇ c´ ateˇ cn´ı nastaven´ı v prostˇ red´ı bez inspektora
Nastaven´a situace v tomto experimentu spoˇc´ıv´a v zobrazen´ı chov´an´ı pasaˇz´er˚ u na ˇspatn´e poˇc´ateˇcn´ı nastaven´ı pravdˇepodobnost´ı kontroly. V situaci se u ´myslnˇe nastav´ı informace pasaˇz´er˚ um (expertn´ı informace), ˇze na u ´seku mezi stanicemi ˇc. 62 a ˇc. 61 s pravdˇepodobnost´ı 10% dojde ke kontrole inspektorem. Avˇsak ve sc´en´aˇri nen´ı nastaven´ y
40
KAPITOLA 7. EXPERIMENTY
ˇz´adn´ y inspektor. C´ıle tohoto experimentu maj´ı zjistit, jak se budou chovat jednotliv´e typy pasaˇz´er˚ u ve stanoven´em sc´en´aˇri.
(a) Pro kapit´ al 25$
(b) Pro kapit´al 75$
Obr´azek 7.2: V´ ysledek experimentu pro ekonomicky zaloˇzen´eho pasaˇz´era
V´ ysledky nejsou pˇr´ıliˇs pˇrekvapiv´e4 ekonomicky slabˇs´ı pasaˇz´er si na poˇca´tku zvol´ı jako nejv´ yhodnˇejˇs´ı pl´an obej´ıt nebezpeˇcn´ yu ´sek. Pasaˇz´er se t´ım p´adem nikdy nedozv´ı, ˇze na obch´azen´em u ´seku nehroz´ı kontrola. Pr˚ ubˇeh je vidˇet na obr´azku 7.2(a). Pasaˇz´er, kter´ y m´a naopak vyˇsˇs´ı poˇc´ateˇcn´ı kapit´al, zvol´ı koupi j´ızdenky a pˇr´ımou cestu metrem. Pr˚ ujezdem zjiˇst’uje, ˇze mˇel p˚ uvodnˇe k dispozici chybnou informaci a upravuje si pravdˇepodobnosti kontroly. T´ım se mu vylepˇsuje uˇzitkov´a hodnota j´ızdy bez j´ızdenky viz. obr´azek 7.2(b).
(a) Pro kapit´ al 25$
(b) Pro kapit´al 75$
Obr´azek 7.3: V´ ysledek experimentu pro ˇcasovˇe zaloˇzen´eho pasaˇz´era 4
Takov´ y byl u ´mysl.
7.2. EXPERIMENTY S JEDNOTLIVCI
41
ˇ Casovˇ e zaloˇzen´ y pasaˇz´er vol´ı pˇr´ımou cestu metrem s j´ızdenkou, protoˇze se mu nevyplat´ı u ´sek obch´azet. Tak´e se uˇc´ı nov´e informace o prostˇred´ı a ˇcasem pˇrejde na j´ızdu bez j´ızdenky viz. obr´azky 7.3. V porovn´an´ı tˇechto dvou pr˚ ubˇeh˚ u s rozd´ıln´ ym denn´ım kapit´alem je vidˇet, jak pasaˇz´erovi s niˇzˇs´ım kapit´alem trv´a d´ele, neˇz zaˇcne jezdit bez j´ızdenky. Jelikoˇz pasaˇz´er m´a niˇzˇs´ı denn´ı rozpoˇcet, neˇz je velikost pokuty, tak by pˇriˇsel o vˇse a nav´ıc by z˚ ustal dluˇzen. Mus´ı si b´ yt jistˇejˇs´ı, ˇze ho na dan´em u ´seku inspektor nechyt´ı.
(a) Pro kapit´ al 25$
Obr´azek 7.4: V´ ysledek
(b) Pro kapit´al 75$
experimentu
pro
psychologicky
zaloˇzen´eho
pasaˇz´era
Psychologicky zaloˇzen´ y pasaˇz´er s niˇzˇs´ım kapit´alem zareaguje na situaci stejn´ ym zp˚ usobem jako ekonomicky zaloˇzen´ y pasaˇz´er se stejn´ ym kapit´alem viz. 7.4 (a). Naopak u pasaˇz´era s vyˇsˇs´ım denn´ım rozpoˇctem pˇrevl´adne psychika pokutov´an´ı a zvol´ı j´ızdenku. Na obr´azku 7.4 (b) je vidˇet, ˇze za simulovanou dobu se mu nezmˇenila v´ yhodnost tohoto pl´anu.
7.2.2
Chybn´ e poˇ c´ ateˇ cn´ı nastaven´ı v prostˇ red´ı s inspektorem
Nastaven´a situace je opaˇcn´a, neˇz v pˇredchoz´ım pˇr´ıpadˇe. Pasaˇz´er opˇet dost´av´a mylnou informaci o prostˇred´ı a zpoˇca´tku se domn´ıv´a, ˇze z poˇca´teˇcn´ıho do c´ılov´eho m´ısta nedoch´az´ı k ˇza´dn´ ym kontrol´am. Avˇsak je zde um´ıstˇen jeden inspektor na vstupn´ı stanici tedy na stanici oznaˇcen´e ˇc. 63. Inspektor na stanici stoj´ı od zapoˇcet´ı pasaˇz´erovy j´ızdy 20 minut.
42
KAPITOLA 7. EXPERIMENTY
(a) Dvacet minut kontroly
(b) Patn´act minut kontroly
Obr´azek 7.5: V´ ysledek v prostˇred´ı s inspektorem na vstupn´ı stanici
V t´eto ˇc´asti se vˇsichni pasaˇz´eˇri zachovaj´ı obdobn´ ym zp˚ usobem. Pasaˇz´eˇri si zpoˇc´atku mysl´ı, ˇze na jeho cestˇe nehroz´ı ˇza´dn´a pokuta, zvol´ı tedy pl´an bez j´ızdenky. Vlivem uˇcen´ı se mu tento pl´an v d˚ usledku kontrol zhorˇs´ı a zkouˇsej´ı tedy nal´ezt jin´ y pl´an. Nakonec vˇsichni postupnˇe pˇrejdou na pl´an s j´ızdenkou viz. obr´azek 7.5(a), kter´ y je reprezentantem. Obch´azen´ı stanice je pro nˇe nev´ yhodn´e z d˚ uvodu vˇetˇs´ı vzd´alenosti a t´ım p´adem tento pl´an obsahuje penalizaci za poruˇsen´ı nutn´eho ˇcasu dojezdu. Ze stejn´eho d˚ uvodu se jim nevyplat´ı posunut´ı zaˇc´atku odjezdu aˇz na dobu, kdy inspektor zmiz´ı. Z d˚ uvodu otestov´an´ı posunut´ı ˇcasu odjezdu byla sn´ıˇzena doba kontroly inspektora na 15 minut. Jak ukazuje vzorov´ y pr˚ ubˇeh na obr´azku 7.5(b) vˇsichni pasaˇz´eˇri si postupnˇe posunuli ˇcasy v´ yjezdu aˇz na v´ yhodnˇejˇs´ı ˇcas a tedy mohli jet bez j´ızdenky.
7.3. EXPERIMENTY S POPULAC´I
7.3
43
Experimenty s populac´ı
Tato sada, kter´a je rozdˇelena na experimenty s jednou linkou a s celou s´ıt´ı, jiˇz pracuje s populac´ı pasaˇz´er˚ u a inspektor˚ u. V experimentech se kv˚ uli zjiˇstˇen´ı dynamiky syst´emu pouˇz´ıvaj´ı dvˇe r˚ uzn´a poˇca´teˇcn´ı nastaven´ı pasaˇz´erova vn´ım´an´ı prostˇred´ı. Jedno s chybnou apriorn´ı informac´ı, kter´a je rovna 0% kontrole a druh´a s pˇresnou, ale uniformˇe rozdˇelenou apriorn´ı informac´ı, kter´a ˇcin´ı na jedn´e lince 1.1% a na cel´e s´ıti 0.7% pravdˇepodobnost kontroly kdekoliv. Pravdˇepodobnosti kontroly byly vypoˇc´ıt´any ze strategi´ı inspektor˚ u. Pasaˇz´eˇri s chybnou apriorn´ı informac´ı se zpoˇc´atku domn´ıvaj´ı, ˇze nemohou b´ yt v syst´emu kontrolovan´ı. Varianta s pˇresnou ale uniformn´ı pravdˇepodobnost´ı je v´ıce realistick´a, pasaˇz´eˇri totiˇz tuˇs´ı, ˇze existuje nˇejak´a pravdˇepodobnost kontroly ale neznaj´ı jej´ı pˇresn´e rozloˇzen´ı.
7.3.1
Experimenty na jedn´ e lince
Parametry pro tuto sadu experiment˚ u zobrazuje tabulka 7.4. Experimenty byly prov´adˇeny jak pro prostˇred´ı s deterministick´ ym5 , tak i nedeterministick´ ym6 rozvrhem inspektor˚ u a pro kaˇzd´e prostˇred´ı s r˚ uznou apriorn´ı informac´ı (0%, 1.1%). Pro kaˇzd´e nastaven´ı prostˇred´ı bylo provedeno ˇsest experiment˚ u (pro r˚ uzn´e poˇcty inspektor˚ u) s pˇeti opakov´an´ımi. V´ ysledky jsou tvoˇreny pr˚ umˇerem tˇechto pˇeti bˇeh˚ u. Tabulka 7.4: Tabulka nastaven´ı experiment˚ u jedn´e linky
Parametry
Hodnoty experiment˚ u
Populace pasaˇz´er˚ u
300
Horizont znalost´ı
100
Populace inspektor˚ u 1 5
10
20
35
50
Dynamika syst´emu je vidˇet v grafech na obr´azc´ıch 7.6, na kter´ ych lze pozorovat v´ yvoj zakoupen´ ych l´ıstk˚ u v pr˚ ubˇehu cel´e simulace. Je moˇzn´e pozorovat, ˇze vˇsechny bˇehy, kter´e jsou urˇceny nastaven´ım situace, konverguj´ı k vyv´aˇzen´emu stavu syst´emu. V prostˇred´ı, s deterministick´ ym rozvrhem inspektor˚ u, je moˇzn´e pozorovat, ˇze m´alo statick´ ych inspektor˚ u nestaˇc´ı pokr´ yt velk´ y prostor. Pasaˇz´eˇri se dan´emu stavu prostˇred´ı pˇrizp˚ usob´ı a je zˇrejm´e, ˇze k pokryt´ı prostˇred´ı a zv´ yˇsen´ı pˇr´ıjm˚ u je nutn´e m´ıt hodnˇe kontrol. 5 6
Inspektor si kaˇzd´ y den prov´ ad´ı stejnou strategii kontroly. Je tedy pˇredv´ıdateln´e jeho chov´an´ı. Inspektor si kaˇzd´ y den vol´ı jinou strategii kontroly.
44
KAPITOLA 7. EXPERIMENTY
V porovn´an´ı s jejich nedeterministick´ ym rozvrhem jiˇz menˇs´ı poˇcet inspektor˚ u staˇc´ı k
300
300
250
250
200
200 Pocet listku [ks]
Pocet listku [ks]
pokryt´ı dan´eho prostˇred´ı.
150 1 inspektor 5 inspektoru 10 inspektoru 20 inspektoru 35 inspektoru 50 inspektoru
100
50
0
0
50
100 Cas [den]
150
150
100
1 inspektor 5 inspektoru 10 inspektoru 20 inspektoru 35 inspektoru 50 inspektoru
50
0
200
0
50
100 Cas [den]
150
200
(b) Chybn´a (0%) apriorn´ı informace s nedetermi-
tick´ ym rozvrhem inspektor˚ u
nistick´ ym rozvrhem inspektor˚ u
300
300
250
250
200
200 Pocet listku [ks]
Pocet listku [ks]
(a) Chybn´a (0%) apriorn´ı informace s determinis-
150
1 inspektor 5 inspektoru 10 inspektoru 20 inspektoru 35 inspektoru 50 inspektoru
100
50
0
0
50
100 Cas [den]
150
150
1 inspektor 5 inspektoru 10 inspektoru 20 inspektoru 35 inspektoru 50 inspektoru
100
50
200
0
0
50
100 Cas [den]
150
200
(c) Pˇresn´a (1.1%) uniformn´ı apriorn´ı informace s (d) Pˇresn´a (1.1%) uniformn´ı apriorn´ı informace s deterministick´ ym rozvrhem inspektor˚ u
deterministick´ ym rozvrhem inspektor˚ u
Obr´azek 7.6: V´ yvoj poˇctu zakoupen´ ych l´ıstk˚ u v prostˇred´ı jedn´e linky
Ust´alen´e stavy poˇctu koupen´ ych j´ızdenek, kontrol a udˇelen´ ych pokut jsou zn´azornˇeny v grafech 7.7. Hodnoty l´ıstk˚ u, pokut a kontrol v grafu jsou v´ ysledkem pr˚ umˇeru za 20 posledn´ıch dn´ı. Zejm´ena na grafu 7.7(b), kde je zn´azornˇen´e nepˇredv´ıdateln´e prostˇred´ı, je vidˇet saturov´an´ı. Za bodem zlomu, kter´ ym je urˇcit´ y poˇcet inspektor˚ u, nen´ı d´ale moˇzn´e zvyˇsovat zisk spoleˇcnosti z d˚ uvodu poctiv´eho chov´an´ı pasaˇz´er˚ u tj. n´ahodn´e kontroly funguj´ı podle oˇcek´av´an´ı. V grafech jsou zobrazeny pr˚ ubˇehy dvou prostˇred´ı, kter´e mˇely odliˇsnou apriorn´ı informaci.
7.3. EXPERIMENTY S POPULAC´I
45
300
350
300 Pocet listku,konrol a pokut [ks]
Pocet listku,konrol a pokut [ks]
250
200
150
100
Pocet koupenych listku (0%) Pocet prepravnich kontrol (0%) Pocet udelenych pokut (0%) Pocet koupenych listku (1.1%) Pocet prepravnich kontrol (1.1%) Pocet udelenych pokut (1.1%)
50
0
0
10
20 30 Pocet inspektoru [ks]
40
250
200
150 Pocet koupenych listku (0%) Pocet prepravnich kontrol (0%) Pocet udelenych pokut (0%) Pocet koupenych listku (1.1%) Pocet prepravnich kontrol (1.1%) Pocet udelenych pokut (1.1%)
100
50
0
50
0
10
20 30 Pocet inspektoru [ks]
40
50
(a) Prostˇred´ı s deterministick´ ym rozvrhem inspek- (b) Prostˇred´ı s nedeterministick´ ym rozvrhem intor˚ u
spektor˚ u
Obr´azek 7.7: Ust´alen´e stavy r˚ uzn´ ych nastaven´ıch situac´ı (jedna linka)
Chov´an´ı pasaˇz´er˚ u lze uk´azat na poˇctu zmˇen, kter´e jsou donuceni udˇelat vlivem adaptace na prostˇred´ı. Na grafech 7.8, kde jsou zobrazeny pr˚ umˇern´e poˇcty7 zmˇen pasaˇz´er˚ uv dan´em prostˇred´ı, je moˇzn´e pozorovat, ˇze pasaˇz´eˇri provedou jednu zmˇenu v deterministick´em prostˇred´ı naopak v´ıce zmˇen v nedeterministick´em, coˇz je logick´e. V histogramu 7.8(a) sice pˇrevl´ad´a poˇcet pasaˇz´er˚ u, kteˇr´ı neudˇelaj´ı ˇz´adnou zmˇenu. Toto je zp˚ usobeno pr˚ umˇerem pˇres vˇsechny pr˚ ubˇehy, protoˇze pˇri mal´em poˇctu inspektor˚ u nemus´ı hodnˇe pasaˇz´er˚ u udˇelat ˇza´dnou zmˇenu. Histogram 7.8(c) ukazuje chov´an´ı v nedeterministick´em prostˇred´ı a zde je vidˇet, ˇze v´ıce pasaˇz´er˚ u zmˇenilo v´ıcekr´at sv˚ uj pl´an. V histogramech 7.8(b) a (d) pak opˇet pˇrevl´ad´a chov´an´ı beze zmˇeny, coˇz je zp˚ usobeno apriorn´ı znalost´ı a
200
180
180
140 120 100 80 60 40 20 0
160 140 120 100 80 60 40 20
0
1
2
3
4
5
6
7
Pocet zmen [−]
(a)
8
9
10
0
250
Pocet pasazeru [ks]
200
160
Pocet pasazeru [ks]
180
Pocet pasazeru [ks]
Pocet pasazeru [ks]
vˇetˇsina pasaˇz´er˚ u se jiˇz rozhodla pro spr´avn´ y pl´an.
160 140 120 100 80 60 40
200
150
100
50
20 0
1
2
3
4
5
6
7
Pocet zmen [−]
(b)
8
9
10
0
0
1
2
3
4
5
6
7
Pocet zmen [−]
(c)
8
9
10
0
0
1
2
3
4
5
6
7
8
9
10
Pocet zmen [−]
(d)
Obr´azek 7.8: Histogramy poˇctu pasaˇz´er˚ u a jejich zmˇen pl´an˚ u (jedna linka) 7
Jedn´ a se o pr˚ umˇer vˇsech nastaven´ı poˇct˚ u inspektor˚ u v dan´em deterministick´em/nedeterministick´em
prostˇred´ı. Podrobn´e informace o zmˇen´ach pasaˇz´er˚ u jsou uvedeny v tabulk´ach A.2, A.3, A.4 a A.5 v pˇr´ıloze A.
46
KAPITOLA 7. EXPERIMENTY Histogram (a) zobrazuje kolik pasaˇz´er˚ u udˇelalo zmˇeny v pl´anu v prostˇred´ı s determis-
nistick´ ym rozvrhem inspektor˚ u a chybnou apriorn´ı informac´ı (0%), (b) v prostˇred´ı s determisnistick´ ym rozvrhem inspektor˚ u a pˇresnou apriorn´ı informac´ı (1.1%), (c) v prostˇred´ı s nedetermisnistick´ ym rozvrhem inspektor˚ u a chybnou apriorn´ı informac´ı (0%) a (d) v prostˇred´ı s nedetermisnistick´ ym rozvrhem inspektor˚ u a pˇresnou apriorn´ı informac´ı (1.1%).
7.3.2
Experimenty na cel´ e s´ıti
Tato sada experiment˚ u byla tak´e prov´adˇena jak pro prostˇred´ı s deterministick´ ym, tak i nedeterministick´ ym rozvrhem inspektor˚ u. Kaˇzd´e prostˇred´ı bylo spouˇstˇeno pro dvˇe r˚ uzn´e apriorn´ı informace (0%, 0.7%). Pro kaˇzd´e nastaven´ı bylo provedeno ˇsest experiment˚ u pro r˚ uzn´ y poˇcet inspektor˚ u v prostˇred´ı a protoˇze se zvˇetˇsilo simulovan´e prostˇred´ı, byl zv´ yˇsen i jejich poˇcet. Parametry experiment˚ u jsou uvedeny v tabulce 7.5.
Tabulka 7.5: Tabulka nastaven´ı experiment˚ u cel´e s´ıtˇe
Parametry
Hodnoty experiment˚ u
Populace pasaˇz´er˚ u
120
Horizont znalost´ı
100
Populace inspektor˚ u 5
10
25
50
75
100
V cel´e s´ıti se nab´ız´ı vyzkouˇset, jak se zachovaj´ı pasaˇz´eˇri na pˇrestupn´ı stanici v pˇr´ıpadˇe, ˇze se na n´ı pravidelnˇe objevuje inspektor. Bylo sledov´ano, jak se v ˇcase zmˇen´ı situace na pˇrestupn´ı a okoln´ıch stanic´ıch (poˇcet pasaˇz´er˚ u, kteˇr´ı pˇres n´ı proj´ıˇzdˇej´ı) pˇred a po um´ıstˇen´ı inspektora. Situace, kter´a je v okol´ı pˇrestupn´ı stanice, je zn´azornˇena na obr´azku 7.9. Vˇsem pasaˇz´er˚ um byla nastavena stejn´a trasa ze stanice ˇc. 21 do ˇc. 62. Na pˇrestupn´ı stanici ˇc. 12 byl v dobu jejich pr˚ ujezdu (pˇrestupu) um´ıstˇen inspektor na 10 minut (ˇcasov´e rozmez´ı 8:50 aˇz 9:00). Pr˚ ujezdnost uk´azan´a v tabulce 7.6 se pˇres pˇrestupn´ı stanici zmˇenila a v´ ysledn´a pr˚ ujezdnost je zn´azornˇena v tabulce 7.7. Z tabulek je vidˇet, ˇze vˇetˇsina pasaˇz´er˚ u vol´ı pˇrestup mezi stanicemi ˇc. 12 a ˇc. 57 a j´ızdu bez j´ızdenky. Postupnˇe se nauˇc´ı, ˇze v danou dobu na stanici je inspektor a tedy zvol´ı posunut´ı ˇcasu odjezdu, nˇekteˇr´ı obch´azej´ı stanici a ostatn´ı si koup´ı j´ızdenku.
7.3. EXPERIMENTY S POPULAC´I
47
Obr´azek 7.9: Zn´azornˇen´ı situace na pˇrestupn´ı stanici
Tabulka 7.6: Poˇcet proj´ıˇzdˇej´ıc´ıch pasaˇz´er˚ u stanicemi pˇred um´ıstˇen´ım inspektora
V´ ystupy ze stanice
Vstupy do stanice
ˇ Cas
ˇc. 11
ˇc. 12
ˇc. 57
ˇc. 56
8:50 - 9:00
4
116
0
0
9:00 - 9:10
0
0
116
0
9:40 - 9:50
0
0
0
4
Tabulka 7.7: Poˇcet proj´ıˇzdˇej´ıc´ıch pasaˇz´er˚ u stanicemi po um´ıstˇen´ım inspektora
V´ ystupy ze stanice
Vstupy do stanice
ˇ Cas
ˇc. 11
ˇc. 12
ˇc. 57
ˇc. 56
8:50 - 9:00
4
23
0
0
9:00 - 9:10
0
39
62
0
9:20 - 9:30
0
0
54
0
9:40 - 9:50
0
0
0
4
Dynamika syst´emu v cel´e dopravn´ı s´ıt’ je obdobn´a jako dynamika jedn´e linky (viz. pˇr´ıloha ...). Na grafech 7.10, ve kter´ ych jsou hodnoty pr˚ umˇerem za 20 posledn´ıch dn´ı, jsou
48
KAPITOLA 7. EXPERIMENTY
zobrazeny v´ ysledky v ust´alen´em stavu tj. poˇcet koupen´ ych j´ızdenek, kontrol a udˇelen´ ych pokut v z´avislosti na poˇctu inspektor˚ u v syst´emu. Na grafu 7.10 vlevo jsou pr˚ ubˇehy s deterministick´ ym rozvrhem inspektor˚ u a vpravo s nedeterministick´ ym. V porovn´av´an´ı v´ ysledk˚ u z obou prostˇred´ı je moˇzn´e si vˇsimnout, ˇze v nepˇredv´ıdateln´em prostˇred´ı je v´ıce pasaˇz´er˚ u donuceno si koupit j´ızdenku, neˇz v prostˇred´ı pˇredv´ıdateln´em.
100
100 Pocet listku,konrol a pokut [ks]
120
Pocet listku,konrol a pokut [ks]
120
80
60 Pocet koupenych listku (0%) Pocet prepravnich kontrol (0%) Pocet udelenych pokut (0%) Pocet koupenych listku (0.7%) Pocet prepravnich kontrol (0.7%) Pocet udelenych pokut (0.7%)
40
20
0
0
20
40 60 Pocet inspektoru [ks]
80
80
60 Pocet koupenych listku (0%) Pocet prepravnich kontrol (0%) Pocet udelenych pokut (0%) Pocet koupenych listku (0.7%) Pocet prepravnich kontrol (0.7%) Pocet udelenych pokut (0.7%)
40
20
0
100
0
20
40 60 Pocet inspektoru [ks]
80
100
(a) Prostˇred´ı s deterministick´ ym rozvrhem inspek- (b) Prostˇred´ı s nedeterministick´ ym rozvrhem intor˚ u
spektor˚ u
Obr´azek 7.10: Ust´alen´e stavy r˚ uzn´ ych nastaven´ıch situac´ı (cel´a s´ıt’)
Histogramy 7.11 ukazuj´ı zmˇeny pr˚ umˇern´eho poˇctu pasaˇz´er˚ u v pˇredv´ıdateln´em obr´azek
45
40 35 30 25 20 15 10 5 0
40 35 30 25 20 15 10 5
0
1
2
3
4
5
6
7
Pocet zmen [−]
(a)
8
9
10
0
0
1
2
3
4
5
6
7
Pocet zmen [−]
(b)
8
9
10
70
70
60
60
Pocet pasazeru [ks]
50
45
Pocet pasazeru [ks]
50
Pocet pasazeru [ks]
Pocet pasazeru [ks]
vlevo a nepˇredv´ıdateln´em prostˇred´ı. I je podobn´a situace jako v prostˇred´ı jedn´e linky.
50 40 30 20 10 0
50 40 30 20 10
0
1
2
3
4
5
6
7
Pocet zmen [−]
(c)
8
9
10
0
0
1
2
3
4
5
6
7
Pocet zmen [−]
8
9
10
(d)
Obr´azek 7.11: Histogramy poˇctu pasaˇz´er˚ u a jejich zmˇen pl´an˚ u (cel´a s´ıt’)
Histogram (a) zobrazuje, kolik pasaˇz´er˚ u udˇelalo zmˇeny v pl´anu v prostˇred´ı s determisnistick´ ym rozvrhem inspektor˚ u a chybnou apriorn´ı informac´ı (0%), (b) v prostˇred´ı s determisnistick´ ym rozvrhem inspektor˚ u a pˇresnou apriorn´ı informac´ı (0.7%), (c) v prostˇred´ı s nedetermisnistick´ ym rozvrhem inspektor˚ u a chybnou apriorn´ı informac´ı (0%) a (d) v prostˇred´ı s nedetermisnistick´ ym rozvrhem inspektor˚ u a pˇresnou apriorn´ı informac´ı (0.7%).
´ ˇ ´I A PAME ˇ TOV ˇ ´ NAROKY ´ 7.4. VYPO CETN E
7.4
49
V´ ypoˇ cetn´ı a pamˇ et’ov´ e n´ aroky
Experimenty byly prov´adˇeny na notebooku vybaven´ y procesorem Intel Core i5 2520M s frekvenc´ı 2,5 GHz (2 j´adra / 4 vl´akna) a pamˇet´ı 4GB DDR3. Dobu v´ ypoˇctu a pamˇet’ov´e n´aroky jsou uvedeny v tabulce 7.8, kter´e ovˇsem z´avis´ı na nastaven´ı simulace (poˇcet pasaˇz´er˚ u a inspektor˚ u). Tabulka 7.8: V´ ypoˇcetn´ı n´aroky simulace
Prostor
Pr˚ umˇern´a v´ ypoˇcetn´ı doba
Pr˚ umˇern´e n´aroky na pamˇet’
Jedna linka Cel´a s´ıt’
1 minuta
0.5 (max 0.5) GB
24 minut
1.2 (max 1.5) GB
V´ ypoˇcetn´ı i pamˇet’ov´ y n´ar˚ ust na cel´e s´ıti, oproti jedn´e lince, je zp˚ usoben vˇetˇs´ım prostorem k prohled´an´ı. V´ ypoˇcetn´ı doba stoupla pˇrev´aˇznˇe kv˚ uli vˇetˇs´ımu faktoru vˇetven´ı a vˇetˇs´ı hloubce prohled´avan´eho stromu. Pamˇet’ov´a n´aroˇcnost se zvˇetˇsila v d˚ usledku vˇetˇs´ı s´ıtˇe s vyˇsˇs´ım poˇctem dopravn´ıch prostˇredk˚ u a pouˇz´ıv´an´ım v´ıce objekt˚ u v pl´anov´an´ı.
7.5
Shrnut´ı
Experimenty ovˇeˇruj´ı, ˇze implementace syst´emu dopravn´ı spoleˇcnosti se chov´a intuitivn´ım zp˚ usobem tj. pasaˇz´eˇri se adaptuj´ı na aktu´aln´ı prostˇred´ı. Dynamika cel´eho syst´emu - pˇri pouˇzit´ı cel´e dopravn´ı s´ıtˇe - je obdobn´a jako dynamika jedn´e linky. Syst´em vˇzdy konverguje do vyv´aˇzen´eho stavu. Z ust´alen´ ych stav˚ u vypl´ yv´a, ˇze v pˇredv´ıdateln´em prostˇred´ı je k sn´ıˇzen´ı u ´niku zisk˚ u zapotˇreb´ı v´ıce inspektor˚ u, neˇz v prostˇred´ı nepˇredv´ıdateln´em, napˇr. v´ ysledn´a hladina u ´niku zisk˚ u jedn´e linky je pro 35 deterministick´ ych inspektor˚ u srovnateln´a s 5 - 10 nedeterministik´ ymi inspektory. N´ahodn´ ym pohybem inspektoˇri pokr´ yvaj´ı vˇetˇs´ı prostor a udrˇzuj´ı rozloˇzenou m´ıru pravdˇepodobnosti kontroly. Pravdˇepodobnost kontroly nebude pˇri statick´ ych kontrol´ach norm´alnˇe rozloˇzen´a. Statiˇct´ı inspektoˇri na sv´ ych m´ıstech zvyˇsuj´ı pravdˇepodobnost kontroly na u ´kor ostatn´ıch m´ıst, ve kter´ ych se pasaˇz´eˇri mohou nauˇcit cestovat, tedy bez j´ızdenky bez vˇetˇs´ıho rizika.
50
KAPITOLA 7. EXPERIMENTY
Kapitola 8 Z´ avˇ er Touto diplomovou prac´ı jsem z´ıskal pˇrehled o siln´e simulaˇcn´ı technice agetn´ıho modelov´an´ı. Experiment´alnˇe jsem vyhodnotil fungov´an´ı modelu s pˇrihl´ednut´ım vlivu jednotliv´ ych faktor˚ u, ovlivˇ nuj´ıc´ıch rozhodnut´ı pasaˇz´era - ekonomick´e, psychologick´e a ˇcasov´e faktory. Do navrˇzen´eho pl´anovac´ıho algoritmu, kter´ y hled´a v´ yhodn´e cesty prostˇredky veˇrejn´e dopravy, jsem implementoval model cestuj´ıc´ıho jako adaptivn´ıho agenta, kter´ y byl integrov´an do simulaˇcn´ı platformy. Pˇrestoˇze je uˇc´ıc´ı model pro pouˇzit´e dopravn´ı s´ıtˇe, kter´e nejsou pˇr´ıliˇs sloˇzit´e, dostaˇcuj´ıc´ı, existuje moˇznost vylepˇsen´ı. Pˇridat do uˇcen´ı f´azi explorace, tedy cel´e uˇcen´ı s posilov´an´ım. N´ahodn´e prozkoum´av´an´ı jin´ ych cest m˚ uˇze v´est k odhalen´ı lepˇs´ıho pl´anu. Dalˇs´ım vylepˇsen´ım pasaˇz´erovy adaptace na prostˇred´ı je pˇrid´an´ı interakce i mezi pasaˇz´ery napˇr. pomoc´ı soci´aln´ıch s´ıt´ı nebo vymˇen ˇov´an´ı informac´ı pˇr´ımo pˇri cestov´an´ı. Takov´a interakce by vedla k prozkoum´an´ı vˇetˇs´ıho tˇreba i nenavˇst´ıven´eho prostoru, coˇz by bylo v´ yhodn´e pˇri cestov´an´ı do m´ıst, kde pasaˇz´er ˇcasto nejezd´ı. Pro re´aln´e pouˇzit´ı a vyhodnocov´an´ı agentn´ıho modelu pˇrepravn´ı kontroly, je nutn´e v simulaci nastavit skuteˇcn´e hodnoty a parametry prostˇred´ı tj. poˇcty pasaˇz´er˚ u atd.. Simulace prov´adˇen´a v t´eto pr´aci je pouze ilustrativn´ı, ovˇeˇruje spr´avnost a intuitivn´ı chov´an´ı navrˇzen´eho modelu. V uveden´e simulaci je m´alo pasaˇz´er˚ u, k pˇretˇeˇzov´an´ı inspektor˚ u nedoch´az´ı, proto je i vˇetˇs´ı pravdˇepodobnosti kontroly, coˇz implikuje rychlejˇs´ı konvergenci syst´emu a ve v´ ysledku m´enˇe udˇelen´ ych pokut. Dopravn´ı spoleˇcnosti by mohly, tuto simulaci vyuˇz´ıvat k maximalizaci zisk˚ u tj. optimalizov´an´ı poˇctu inspektor˚ u v syst´emu, d´ale k ovˇeˇren´ı a vytv´aˇren´ı strategi´ı sv´ ych inspektor˚ u. Ze sledov´an´ı dynamiky syst´emu by mohly zjiˇst’ovat, po jak´e dobˇe je dobr´e zmˇenit rozloˇzen´ı nebo cel´e strategie inspektor˚ u, protoˇze se jiˇz pasaˇz´eˇri nauˇcili jejich chov´an´ı. Simul´ator poskytuje ovˇeˇren´ı, ˇze pˇri statick´ y kontrol´ach je potˇrebn´e pouˇz´ıt vˇetˇs´ı mnoˇzstv´ı 51
´ ER ˇ KAPITOLA 8. ZAV
52
inspektor˚ u, coˇz sniˇzuje zisky dopravn´ı spoleˇcnosti. Ze simulace dopravn´ıho syt´emu je moˇzn´e z´ıskat informace o stanic´ıch ˇci u ´sec´ıch, na kter´ ych se nejv´ıce pohybuj´ı neplat´ıc´ı pasaˇz´eˇri, tedy m˚ uˇze odhalit slab´a m´ısta v kontrole dopravn´ı s´ıtˇe. Problematika pˇrepravn´ı kontroly v s´ıt´ıch veˇrejn´e dopravy a u ´niky zisk˚ u dopravn´ıch spoleˇcnost´ı, kter´e zp˚ usobuj´ı neplat´ıc´ı cestuj´ıc´ı, je aktu´aln´ı t´ema. U simulovan´eho modelu jsou zohlednˇeny pouze hlavn´ı faktory ovlivˇ nuj´ıc´ı rozhodov´an´ı pasaˇz´er˚ u, kter´ ych je ve skuteˇcnosti v´ıce. S v´ yvojem nov´ ych technologi´ı jako jsou napˇr. chytr´e telefony a aplikace v nich, vznik´a vˇetˇs´ı a rychlejˇs´ı informovanost o prostˇred´ı veˇrejn´e dopravy a t´ım i rychlejˇs´ı adaptace cestuj´ıc´ıch. I na takovou situaci se vyvinut´ y model, d´ıky vyuˇzit´ı agentn´ıho modelov´an´ı, d´a snadno pˇrizp˚ usobit.
Literatura Bonabeau, E. (2002), ‘Agent-based modeling: Methods and techniques for simulating human systems’, 99, 7280–7287. Borshchev, A. a Filippov, A. (2004), From system dynamics and discrete event to practical agent based modeling: Reasons, techniques, tools, in ‘The 22nd International Conference of the System Dynamics Society’. Boyd, C., Martini, Ch., Rickard, J. a Russell, A. (1989), ‘Fare evasion and noncompliance’, Journal of transport economics and policy . Clausen, J. (1999), Branch and bound algoritms - principles and examples, Technical report, Deparment of Computer Science, University of Copenhagen. Jakob, M., Moler, Z., Komenda, A., Yin, Z., Jiang, A., X., Johnson, M., P., Pechoucek, M. a Tambe, M. (2012), Agentpolis: Towards a platform for fully agent-based modeling of multi-modal transportation (demonstration), in ‘12th International Conference on Autonomous Agents and Multiagent Systems (to appear)’. Killias, M., Scheidegger, D. a Nordenson, P. (2009), ‘The effects of increasing the certainty of punishment’, Criminology 6. Kooreman, P. (1993), ‘Fare evasion as a result of expected utility maximisation’, Journal of transport economics and policy . Leung, L., F. (2003), A study of fare evasion in railway systems in hong kong, Master’s thesis, The University of Hong Kong. Moler, Z. (n.d.), Fare evasion abm req spec and analysis example, Technical report, ˇ Katedra poˇc´ıtaˇc˚ u, FEL CVUT. Nagy, I., Nedoma, P., Kratky, M., Pavelkova, L. a Ettler, P. (2002), ‘O bayesovsk´em uˇcen´ı’, Automa 7, 56–61. 53
54
LITERATURA
Piliavin, I., Thornton, C., Gartner, R. a Matsueda, R., L. (1986), ‘Crime, deterrence, and rational choice’, American Sociological Review 51, 101–119. Vidal, J., M. (2010), Fundamentals of multiagent system with netlogo examples, in ‘Fundamentals of Multiagent System’.
Pˇ r´ıloha A N´ azev pˇ r´ılohy
Rozmez´ı pˇ r´ıjm˚ u
Poˇ cet dom´ acnost´ı
M´enˇe neˇz $10,000:
169738
$10,000 aˇz $14,999:
96131
$15,000 aˇz $19,999:
91434
$20,000 aˇz $24,999:
90634
$25,000 aˇz $29,999:
83895
$30,000 aˇz $34,999:
79625
$35,000 aˇz $39,999:
68487
$40,000 aˇz $44,999:
62754
$45,000 aˇz $49,999:
54614
$50,000 aˇz $59,999:
91959
$60,000 aˇz $74,999:
106186
$75,000 aˇz $99,999:
107198
$100,000 aˇz $124,999:
61105
$125,000 aˇz $149,999:
33453
$150,000 aˇz $199,999:
32418
$200,000 a v´ıce:
46978
Tabulka A.1: Pˇr´ıjmy pr˚ umˇern´e tˇr´ıˇclenn´e dom´acnosti v LA
I
ˇ ´ILOHA A. NAZEV ´ ˇ ´ILOHY PR PR
II
Tabulka A.2: Zmˇeny pl´an˚ u v pˇredv´ıdateln´em prostˇred´ı, 0%, jedna linka
Poˇcet pasaˇz´er˚ u Poˇ cet inspektor˚ u
0
1
2
3
4
5 6
7
8
9
10 a v´ıce
1 inspektor
288
10
0
2
0
0
0
0
0
0
0
5 inspektor˚ u
243
51
0
4
0
2
0
0
0
0
0
10 inspektor˚ u
200
88
1
8
0
1
0
0
0
0
2
25 inspektor˚ u
129 155
5
7
1
0
1
0
0
0
2
35 inspektor˚ u
65
208
6
9
5
3
1
0
0
0
3
50 inspektor˚ u
50
230
5
4
3
2
1
2
0
0
3
Tabulka A.3: Zmˇeny pl´an˚ u v nepˇredv´ıdateln´em prostˇred´ı, 0%, jedna linka
Poˇcet pasaˇz´er˚ u Poˇ cet inspektor˚ u
0
1
2
3
4
5
6
7
8
9
10 a v´ıce
1 inspektor
168
79
15
19
8
4
5
2
0
0
0
5 inspektor˚ u
68
156
10
38
4
12
5
2
1
2
2
10 inspektor˚ u
39
198
8
32
4
7
1
3
0
1
7
25 inspektor˚ u
20
234
10
16
7
3
0
4
1
0
5
35 inspektor˚ u
14
247
9
14
3
3
1
1
0
1
7
50 inspektor˚ u
11
253
11
12
2
4
0
0
1
0
6
Tabulka A.4: Zmˇeny pl´an˚ u v pˇredv´ıdateln´em prostˇred´ı, 1.1%, jedna linka
Poˇcet pasaˇz´er˚ u Poˇ cet inspektor˚ u
2
3
4
5 6
7
8
9
10 a v´ıce
193 107
0
0
0
0
0
0
0
0
0
5 inspektor˚ u
195
97
3
4
1
0
0
0
0
0
0
10 inspektor˚ u
197
96
0
5
0
0
0
0
1
0
1
25 inspektor˚ u
197
91
6
3
1
0
1
0
0
0
1
35 inspektor˚ u
203
84
9
2
1
0
0
0
0
0
1
50 inspektor˚ u
196
96
3
2
2
0
0
0
0
0
1
1 inspektor
0
1
III Tabulka A.5: Zmˇeny pl´an˚ u v nepˇredv´ıdateln´em prostˇred´ı, 1.1%, jedna linka
Poˇcet pasaˇz´er˚ u Poˇ cet inspektor˚ u
0
1
2
3
4
5
6
7
8
9
10 a v´ıce
1 inspektor
239 22
19
4
9
4
0
1
2
0
0
5 inspektor˚ u
237 27
8
14
3
4
2
0
1
1
3
10 inspektor˚ u
215 50
11
14
2
4
1
2
0
1
0
25 inspektor˚ u
203 79
2
11
0
1
1
1
0
0
2
35 inspektor˚ u
197 81
11
8
0
1
0
0
0
0
2
50 inspektor˚ u
196 93
2
5
1
0
0
0
1
0
2
ˇ ´ILOHA A. NAZEV ´ ˇ ´ILOHY PR PR
120
120
100
100
80
80 Pocet listku [ks]
Pocet listku [ks]
IV
60
5 inspektoru 10 inspektoru 25 inspektoru 50 inspektoru 75 inspektoru 100 inspektoru
40
20
0
0
50
100 Cas [den]
150
60
5 inspektoru 10 inspektoru 25 inspektoru 50 inspektoru 75 inspektoru 100 inspektoru
40
20
0
200
0
50
100 Cas [den]
150
200
(b) Chybn´a (0%) apriorn´ı informace s nedetermi-
tick´ ym rozvrhem inspektor˚ u
nistick´ ym rozvrhem inspektor˚ u
120
120
100
100
80
80 Pocet listku [ks]
Pocet listku [ks]
(a) Chybn´a (0%) apriorn´ı informace s determinis-
60
5 inspektoru 10 inspektoru 25 inspektoru 50 inspektoru 75 inspektoru 100 inspektoru
40
20
0
0
50
100 Cas [den]
150
60
5 inspektoru 10 inspektoru 25 inspektoru 50 inspektoru 75 inspektoru 100 inspektoru
40
20
0
200
0
50
100 Cas [den]
150
200
(c) Pˇresn´a (0.7%) uniformn´ı apriorn´ı informace s (d) Pˇresn´a (0.7%) uniformn´ı apriorn´ı informace s deterministick´ ym rozvrhem inspektor˚ u
deterministick´ ym rozvrhem inspektor˚ u
Obr´azek A.1: V´ yvoj poˇctu zakoupen´ ych l´ıstk˚ u v prostˇred´ı cel´e s´ıtˇe
Tabulka A.6: Zmˇeny pl´an˚ u v pˇredv´ıdateln´em prostˇred´ı, 0%, cel´a s´ıt’
Poˇcet pasaˇz´er˚ u Poˇ cet inspektor˚ u
0
1
2
3
4
5
6
7
8
9
10 a v´ıce
5 inspektor˚ u
98
15
1
0
0
1
0
0
1
0
4
10 inspektor˚ u
80
30
2
2
0
0
0
0
0
0
6
25 inspektor˚ u
42
53
5
4
2
1
2
1
0
1
9
50 inspektor˚ u
30
62
6
4
2
3
5
0
2
0
6
75 inspektor˚ u
22
65
6
8
5
2
3
1
1
0
7
100 inspektor˚ u
20
70
5
6
5
1
4
1
2
0
6
V Tabulka A.7: Zmˇeny pl´an˚ u v nepˇredv´ıdateln´em prostˇred´ı, 0%, cel´a s´ıt’
Poˇcet pasaˇz´er˚ u Poˇ cet inspektor˚ u
0
1
2
3
4
5
6
7
8
9
10 a v´ıce
5 inspektor˚ u
27 48
4
16
4
4
2
4
0
2
9
10 inspektor˚ u
21 59
2
10
6
4
4
2
1
1
10
25 inspektor˚ u
17 68
3
9
2
1
3
4
1
3
9
50 inspektor˚ u
15 71
5
7
3
4
0
4
2
0
9
75 inspektor˚ u
15 73
6
7
3
2
0
2
3
1
8
100 inspektor˚ u
15 73
7
5
5
0
2
2
3
0
8
Tabulka A.8: Zmˇeny pl´an˚ u v pˇredv´ıdateln´em prostˇred´ı, 0.7%, cel´a s´ıt’
Poˇcet pasaˇz´er˚ u Poˇ cet inspektor˚ u
0
1
2
3
4
5
6
7
8
9
10 a v´ıce
5 inspektor˚ u
97 15
1
2
0
0
0
0
1
0
4
10 inspektor˚ u
80 30
2
2
0
0
0
0
0
0
6
25 inspektor˚ u
40 57
4
2
2
3
0
1
1
0
10
50 inspektor˚ u
30 62
6
4
2
3
5
0
2
0
6
75 inspektor˚ u
22 65
6
8
5
2
3
1
1
0
7
100 inspektor˚ u
21 70
5
6
5
2
2
1
2
0
6
Tabulka A.9: Zmˇeny pl´an˚ u v nepˇredv´ıdateln´em prostˇred´ı, 0.7%, cel´a s´ıt’
Poˇcet pasaˇz´er˚ u 1
2
3
4
5
6
7
8
9
10 a v´ıce
5 inspektor˚ u
28 47
3
13
3
8
5
1
1
0
11
10 inspektor˚ u
21 57
4
12
5
6
0
3
2
2
8
25 inspektor˚ u
17 67
4
11
2
3
0
2
4
1
9
50 inspektor˚ u
16 72
7
5
3
2
1
3
1
1
9
75 inspektor˚ u
15 73
6
7
3
2
0
2
3
1
8
100 inspektor˚ u
15 73
8
5
7
0
2
0
1
1
8
Poˇ cet inspektor˚ u
0
VI
ˇ ´ILOHA A. NAZEV ´ ˇ ´ILOHY PR PR
Pˇ r´ıloha B Obsah pˇ riloˇ zen´ eho CD K t´eto pr´aci je pˇriloˇzeno CD, na kter´em jsou uloˇzeny zdrojov´e k´ody. ˇ anky: • Cl´ • Experimenty: • Obr´azky: – Grafy – Diagramy • Zdrojov´e k´ody: – Pomocn´e k´ody: ∗ K´ody pro pˇredzpracov´an´ı ∗ K´ody pro generov´an´ı dat ∗ Spouˇstˇec´ı k´od ∗ Vyhodnocovac´ı k´od – Pˇrepravn´ı kontrola – Spustiteln´ y soubor
VII