Kapitola 1 Vícekriteriální rozhodování za jistoty Při řešení rozhodovacích problémů se často setkáváme s případy, kdy optimální rozhodnutí musí vyhovovat více než jednomu kritériu. Zadaná kritéria mohou mít kvantitativní i kvalitativní charakter (při koupi automobilu je rozhodující jak jeho cena, tak i vzhled), mohou být maximalizační i minimalizační (požadujeme, aby zakoupený automobil dosahoval co největší rychlosti a byl co nejlacinější) a mohou být i navzájem konfliktní (nízká cena výrobku je zpravidla spojena s jeho horší kvalitou). Úlohy vícekriteriálního rozhodování můžeme klasifikovat podle způsobu zadání množiny variant, které pro optimální rozhodnutí připadají v úvahu (jde o tzv. přípustné varianty). Je-li tato množina určena konečným seznamem variant, hovoříme o vícekriteriálním hodnocení variant. Je-li množina přípustných variant zadána podmínkami, které musí být při výběru optimální varianty splněny, jde o úlohy vícekriteriálního programování (též vícekriteriální nebo vektorové optimalizace). V těchto úlohách varianty rozhodnutí představují n-tice nezáporných čísel, které vyhovují daným omezujícím podmínkám a kterých může být nekonečně mnoho. Kritéria pro výběr nejvýhodnější varianty jsou vyjádřena účelovými funkcemi a musí být tedy pouze kvantitativní. Jednoduchým příkladem vícekriteriálního hodnocení variant je následující úloha, na které budou ilustrovány všechny používané pojmy a metody. Příklad 1. Uchazeč o zaměstnání se rozhoduje mezi firmami A, B, C, přičemž tato pracoviště posuzuje podle výše měsíčního platu (tis.Kč), doby strávené na cestě do zaměstnání (minuty), možnosti dalšího odborného růstu (hodnocení 1, 2, 3 pro malou, střední a velkou možnost) a začátku pracovní doby (hodiny:minuty). Potřebné údaje jsou uvedeny v tabulce
firma A firma B firma C
1.1
K1 30 22 26
K2 60 30 45
K3 2 1 3
K4 9:00 7:30 8:00
Základní pojmy
Rozhodnutí - výběr jedné nebo více variant z množiny všech přípustných variant. V našem případu je to výběr firmy, u které se uchazeč nechá zaměstnat. Rozhodovatel – subjekt, který má za úkol učinit rozhodnutí. Člověk, který vybírá vhodné zaměstnání. V úlohách vícekriteriální analýzy variant je dána konečná (diskrétní) množina m variant, které jsou hodnoceny podle n kritérií. Cílem je učinit rozhodnutí, která varianta je podle daných kritérií hodnocena nejlépe. Jedná se o tzv. optimální variantu. Varianty lze seřadit od nejlepší po nejhorší nebo je rozdělit na efektivní a neefektivní varianty. 1
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
2
Varianty (alternativy) – konkrétní rozhodovací možnosti, které jsou realizovatelné. V následujícím textu je budeme značit Ai (pro i = 1, 2, . . . , m ). Varianty jsou firmy A, B, C. Kritéria – hlediska, ze kterých jsou varianty posuzovány. Dále je budeme značit Kj (pro j = 1, 2, ..., n ). Kritéria jsou výše měsíčního platu, doba cesty do zaměstnání, možnost dalšího odborného růstu, začátek pracovní směny. Kriteriální matice - je-li hodnocení variant podle kritérií kvantifikováno, údaje uspořádáváme do kriteriální matice Y = (yi,j ). Prvky této matice vyjadřují hodnocení i-té varianty podle j-tého kritéria. Řádky odpovídají variantám, sloupce kritériím. Klasifikace kritérií dle povahy • maximalizační – nejlepší hodnoty mají nejvyšší hodnoty (výše měsíčního platu, možnost odborného růstu, začátek pracovní doby) • minimalizační – nejlepší hodnoty mají nejmenší hodnoty (doba cesty do zaměstnání) Vhodné je před hodnocením převést všechna kritéria na jeden typ. Pokud chceme například převést minimalizační kritérium na maximalizační, vybereme ve sloupci příslušného kritéria největší číslo a od tohoto čísla odečítáme ostatní kriteriální hodnoty v daném sloupci. Výsledkem je potom vzdálenost skutečné hodnoty od hodnoty nejhorší, čím je tato vzdálenost větší, tím lépe, kritérium je tedy maximalizační. V našem příkladu je minimalizační K2 , pokud hodnoty budeme chtít převést na maximalizační, vybereme největší číslo ve druhém sloupci (60), pro všechny prvky v druhém sloupci vypočteme rozdíl mezi nejhorší hodnotou (60) a ostatními hodnotami a získáme čísla 0, 30, 15. Víme tedy například, že při volbě varianty B budeme dojíždět do práce o 15 minut kratší dobu než u varianty A. Většina programů pro vícekriteriální hodnocení variant vyžaduje pouze zadání typu kritéria a program provede standardizaci všech kritérií na výnosový typ sám. dle kvantifikovatelnosti • kvantitativní – objektivně měřitelné údaje (výše měsíčního platu, doba cesty do zaměstnání, začátek pracovní doby) • kvalitativní – nelze objektivně měřit, varianty jsou hodnoceny slovně, proto je nutné užít k převedení slovního hodnocení různé bodovací stupnice či relativní hodnocení variant (možnost dalšího odborného růstu) Preference kritéria – důležitost kritéria v porovnání s ostatními kritérii. Vyjádření preference • aspirační úroveň – hodnota kritéria, které má být dosaženo (například uchazeč požaduje měsíční plat alespoň 25 tis. Kč, aby se nabídkou práce začal zabývat) • pořadí kritérií (ordinální informace o kritériích) - posloupnost kritérií od nejdůležitějšího po nejméně důležité • váhy kritérií – kardinální informace o kritériích; váha je hodnota z intervalu a vyjadřuje relativní důležitost kritéria v porovnání s ostatními • kompenzace kriteriálních hodnot – jsou vyjádřeny mírou substituce mezi kriteriálními hodnotami (možno vyrovnat špatné kriteriální hodnoty podle jednoho kritéria lepšími hodnotami podle jiného kritéria)
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
3
Varianty se speciálními vlastnostmi Dominovaná varianta – pokud jsou všechna kritéria maximalizační, varianta ai dominuje variantu aj pokud existuje alespoň jedno kritérium kl , že yil > yjl , přičemž pro ostatní kritéria platí (yi1 , yi2 , . . . , yin ) ≥ (yj1 , yj2 , . . . , yjn ). V našem případě taková varianta neexistuje, ale byla by to například taková, která by měla výši měsíčního platu 26 tisíc Kč, doba dojíždění by byla 60 minut, možnost odborného růstu by byla malá a začátek pracovní doby by byl v 7:30 Paretovská varianta, nedominovaná varianta – varianta, která není dominovaná žádnou jinou variantou. Ideální varianta – hypotetická či reálná varianta, která dosahuje ve všech kritériích nejlepší možné hodnoty. Taková varianta by dominovala všechny ostatní varianty. Ideální by bylo, kdyby měsíční plat byl 30 tis. Kč, doba dojíždění by byla 30 minut, u firmy by byla velká možnost odborného růstu a začátek pracovní doby by byl v 9:00. Taková skutečná firma (varianta) v seznamu možných variant není, proto je tato varianta pouze hypotetická. Kdyby taková skutečně existovala, uchazeč by si ji vybral a nemusel by hledat kompromisní řešení. Bazální varianta – hypotetická či reálná varianta, jejíž ohodnocení je nejhorší podle všech kritérií. Taková varianta by byla dominovaná ostatními variantami. V našem případě by taková varianta zahrnovala výši měsíčního platu 22 tis. Kč, dojíždění by trvalo 60 minut, možnost odborného růstu by byla malá a pracovní doba by začínala v 7:30. Taková varianta mezi skutečnými variantami není, uchazeč by ji jinak mohl rovnou vyřadit, byla by dominovaná. Kompromisní varianta – jediná nedominovaná varianta doporučená k řešení, vybraná podle různých pravidel - viz dále. Vlastnosti, které by měla mít kompromisní varianta: • nedominovanost – varianta nesmí být dominovaná jinou variantou • invariance vzhledem k pořadí kritérií – pořadí kritérií neovlivňuje výběr kompromisní varianty, • invariance vzhledem k měřítku kriteriálních hodnot – pokud ke všem prvkům přičteme stejné číslo (vynásobíme stejným číslem), množina vybraných variant nebo vybraná varianta se nesmí změnit, • nezávislost na identických hodnotách téhož kritéria – vyskytne-li se kritérium, jehož hodnoty jsou pro všechny varianty zhruba stejné, nesmí se změnit množina vybraných variant • invariance vzhledem k přidaným dominovaným variantám – přidáme-li do množiny variant dominovanou variantu, vybraná kompromisní varianta se nesmí změnit. • determinovanost – podle každého přístupu nejméně jedna varianta musí být vybrána jako kompromisní • jednoznačnost – zvolený postup dává jednoznačný výsledek, jednu variantu označí jako kompromisní
1.2
Metody stanovení vah kritérií
Většina metod vícekriteriálního rozhodování vyžaduje odlišení jednotlivých kritérií z hlediska jejich významnosti. Jednou z možností je číselné vyjádření této významnosti pomocí tzv. vah (čím je kritérium významnější, tím je jeho váha větší). Váhu kritéria Kj budeme značit vj , j = 1, 2, . . . , n, kde n je počet všech uvažovaných kritérií. Aby váhy kritérií, stanovené různými metodami, popř. různými experty, byly srovnatelné, vyjadřujeme je v normovaných hodnotách wj , které počítáme podle vztahu 1.2 vj wj = P , j = 1, 2, . . . , n (1.1) n vk k=1
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
4
Normované váhy představují nezáporná čísla, jejichž součet se rovná jedné. Rozdělení metod pro stanovení vah kritérií Metody na stanovení vah kritérií lze rozdělit podle informace, která je nutná ke stanovení vah. • rozhodovatel nemůže určit preference V případě, že rozhodovatel není schopen rozlišit důležitost jednotlivých kritérií, všem kritériím je přiřazena stejná váha. Máme-li tedy například pět kritérií (n = 5), každému z nich je přiřazena váha 0,2 (wj = n1 ). • rozhodovatel má ordinální informaci o kritériích V takovém případě je rozhodovatel schopen určit pořadí důležitosti kritérií. Mezi metody vyžadující ordinální informaci o kritériích patří metoda pořadí a Fullerova metoda • rozhodovatel má kardinální informace o kritériích rozhodovatel zná nejen pořadí, ale i rozestupy v pořadí preferencí mezi jednotlivými kritérii. Mezi metody založené na tomto principu patří bodovací metoda a Saatyho metoda Pro vlastní hodnocení variant stačí, aby si rozhodovatel vybral jednu metodu, tou spočítal váhy a s těmito váhami počítal dále. Metoda pořadí Rozhodovatel seřadí kritéria K1 , K2 , . . . , Kn od nejvýznamnějšího k nejméně významnému a takto uspořádaným kritériím přiřadí váhy n, n − 1, . . . , 2, 1. Pro normovanou váhu kritéria Kj s vahou vj pak platí vj , j = 1, 2, . . . , n. (1.2) wj = P n vk k=1
Řešený příklad 1. Spočítejte váhy kritérií z příkladu 1 metodou pořadí, pokud víte, že uchazeč preferuje kritéria v tomto pořadí: K1 , K3 , K2 , K4 Řešení. Například kritérium K1 je první v pořadí ze čtyř kritérií, proto přiřadíme tomuto kritérii 4 body. Celkový počet bodů přidělený všem kritériím je 10. Váha se pak vypočte jako podíl bodů přiřazených tomuto kritériu na celkovém počtu přiřazených bodů, tedy w1 = 4/10 = 0, 4. Stejně tak přiřazujeme body podle pořadí dalším kritériím a váhy jsou pak následující: K3 = 0, 3, K2 = 0, 2 a K4 = 0, 1. 2 Fullerova metoda Při větším počtu kritérií je výhodné srovnávat navzájem vždy pouze dvě kritéria, o kterých snáze rozhodneme, které je důležitější. Jednu z možností pro vyhodnocení těchto srovnání poskytuje tzv. Fullerův trojúhelník. Za předpokladu, že jednotlivá kritéria jsou pevně očíslována pořadovými čísly 1, 2, , n, Fullerův trojúhelník je tvořen dvojřádky, v nichž každá dvojice kritérií se vyskytne právě jednou (viz schéma). U každé dvojice hodnotitel zakroužkuje nebo jinak vyznačí číslo toho kritéria, které považuje za důležitější, takže pro kritérium Kj představuje počet zakroužkovaných čísel j počet jeho preferencí, který označíme fj . Protože při počtu kritérií n je počet párových srovnání roven kombinačnímu číslu , tj. pro normovanou váhu kritéria Kj platí wj = Schéma Fullerova trojúhelníku
fj n(n−1) 2
, j = 1, 2, . . . , n
(1.3)
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY 1 2
1 3 2 3
1 4 2 4
... ... ... ... ... n-2 n-1
5
1 n 2 n n-2 n n-1 n
Řešený příklad 2. Spočítejte váhy kritérií z příkladu 1 Fullerovou metodou, pokud víte, že uchazečovy preference jsou následující: K1 K3 K2 K4 . Řešení. Zapíšeme každou dvojici kritérií do Fullerova trojúhelníku a tučně označíme to kritérium, které je ve dvojici preferováno. 1 2
1 1 3 4 2 2 3 4 3 4
Pro každé kritérium spočítáme kolikrát je označené jako preferované před jiným kritériem Kritérium Počet preferencí Váha K1 3 1/2 K2 1 1/6 K3 2 1/3 K4 0 0 Celkem 6 1 2 Nevýhodou metody párového srovnávání je skutečnost, že nejméně důležité kritérium má nulovou váhu, i když nemusí jít o zcela bezvýznamné kritérium. Tento nedostatek lze odstranit tak, že četnost preferencí každého kritéria zvýšíme o 1 a jmenovatele zlomku ve vzorci (1.3) zvýšíme o n. Řešený příklad 3. Spočítejte váhy Fullerovou metodou tak, aby žádné z kritérií nemělo nulovou váhu. Preference jsou shodné se zadáním řešeného příkladu 2. Řešení. Vezmeme výsledky řešeného příkladu 2, navýšíme počet preferencí o jednotku, tím nám vzroste celkový počet porovnávání na 10 a váhy se pak počítají jako podíl počtu preferencí daného kritéria a celkového počtu porovnávání. Kritérium Počet preferencí Váha Navýšený počet preferencí K1 3 1/2 4 K2 1 1/6 2 K3 2 1/3 3 K4 0 0 1 Celkem 6 1 10
Upravená váha 0,4 0,2 0,3 0,1 1 2
Existují modifikace metody párového srovnávání kritérií, které připouštějí stejnou důležitost nebo nesrovnatelnost některých kritérií, ale těmi se v tomto kurzu neudeme zabývat. Pokud tedy budeme
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
6
chtít použít Fullerovu metodu, musíme být schopni pomocí relace preference kritéria uspořádat. Bodovací metoda Na rozdíl od metody pořadí, která vychází pouze z porovnání významnosti jednotlivých kritérií, při bodovací metodě se důležitost kritérií ohodnotí počtem bodů (čím je kritérium důležitější, tím má větší počet bodů). Bodovací stupnice může mít větší či menší rozsah – např. 1 až 5, 1 až 10 apod. Přidělený počet bodů se převádí na normovanou váhu dle vzorce (1.2). Zvláštním případem bodovací metody je alokace 100 bodů (zvaná též Metfesselova alokace), kdy mezi jednotlivá kritéria se v souladu s jejich důležitostí rozděluje 100 bodů. Normované váhy jsou potom stokrát menší než příslušný počet bodů. Řešený příklad 4. Předpokládejte, že vy sami jste uchazečem o zaměstnání z příkladu 1. Spočítejte váhy kritérií z příkladu 1 bodovací metodou. Řešení. Nejprve si musíme nejen srovnat kritéria podle pořadí preferencí, ale i určit sílu těchto preferencí. Jedno z možných bodových ohodnocení spolu s výslednými váhami pro takové bodové ohodnocení je v následující tabulce. Kritérium Počet bodů K1 50 K2 20 K3 25 K4 5 Celkem 100
Váha 0,5 0,2 0,25 0,05 1 2
Metoda kvantitativního párového srovnávání (Saatyho metoda) Kromě výběru preferovaného kritéria se určuje pro každou dvojici kritérií také velikost této preference (SAATY,1990). K vyjádření velikosti preferencí SAATY doporučuje bodovou stupnici: číselné 1 3 5 7 9
Vyjádření preferencí slovní kritéria jsou stejně významná první kritérium je slabě významnější než druhé první kritérium je silně významnější než druhé první kritérium je velmi silně významnější než druhé první kritérium je absolutně významnější než druhé
Pro citlivější vyjádření preferencí je možné použít i mezistupně (2, 4, 6, 8). Velikost preferencí i-tého kritéria proti j-tému můžeme uspořádat do Saatyho matice S, jejíž prvky sij představují odhady podílů vah kritérií (kolikrát je jedno kritérium významnější než druhé): sij ≈
vi , i, j = 1, 2, . . . , n vj
(1.4)
Matice S je čtvercová řádu r × n a pro prvky matice S platí sij =
1 , i, j = 1, 2, . . . , n sji
(1.5)
tedy matice S je reciproční. Na diagonále matice S jsou vždy hodnoty jedna (každé kritérium je samo sobě rovnocenné). Dříve než se počítají váhy jednotlivých kritérií, je nutné ověřit, zda zadaná matice párových porovnávání je konzistentní. Uvažujme ideální matici V = (vij ), pro jejíž prvky by platilo vhj = vhi vij pro i, j, h = 1, 2, . . . , n. Taková matice by byla dokonale konzistentní.
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
7
Příklad 2. Mějme matici párových srovnávání pro tři kritéria:
K1 K2 K3
K1 K 2 K 3 1 2 6 1/2 1 3 1/6 1/3 1
Prvky této matice jsou plně konzistentní, neboť platí s13 = s12 s23 nebo-li 6 = 2 · 3. Pokud je kritérií více, je téměř nemožné zadat odhady vah kritérií tak, aby matice byla dokonale konzistentní. V takovém případě se počítá míra konzistence. V rámci našeho kurzu nebudeme míru konzistence počítat, konzistenci budeme sledovat pouze ve výsledcích, které budou k dispozici díky užití softwaru. V případě zájmu o výpočet indexu konzistence doporučujeme literaturu. . . . . . . . Při stanovování vah můžeme vycházet z podmínky, že matice S by se měla od matice V lišit co nejméně. Potom minimalizujeme součet odchylek stejnolehlých prvků obou matic: n X n X vi 2 F = sij − ] → min (1.6) vj i=1 j=1 za podmínky
r P
vj = 1 a vj ≥ 0 pro i, j = 1, 2, . . . , n.
j=1
Další metoda, jak stanovit váhy je logaritmická metoda nejmenších čtverců. Řešíme F =
n X n X
[ln sij − (ln vi − ln vj )]2 → min
(1.7)
i=1 j=1
za podmínky
r P
vj = 1 a vj ≥ 0 pro i, j = 1, 2, . . . , n.
j=1
Saaty navrhl početně jednoduchý způsob, jak spočítat váhy. Řešením je normalizovaný geometrický průměr řádků matice S: " # n1 r Q sij j=1
wi = n P i=1
"
n Q
(1.8)
# n1 sij
j=1
Pokud je matice S plně konzistentní, popřípadě dostatečně konzistentní, váhy kritérií vypočítané podle vztahu 1.8 odpovídají požadavkům na jejich preferenci. Pokud matice S není konzistentní, je nutné upravit odhady důležitosti jednotlivých kritérií v původní matici S a tím zlepšit jejich konzistenci. Řešený příklad 5. Nejprve provedeme u našeho příkladu 1 porovnání důležitosti mezi všemi dvojicemi kritérií a tato porovnání uspořádáme do Saatyho matice párových srovnání. K1 K2 K3 K4
K1 K2 K 3 K4 1 5 3 7 1/5 1 1/2 4 1/3 2 1 5 1/7 1/4 1/5 1
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
8
Řešení. Vzhledem k tomu, že matice je dostatečně konzistentní (zjistili jsme v softwaru Sanna), nemusíme upravovat sílu preferencí mezi kritérii a přistoupíme k výpočtu vah podle vztahu 1.8. geometrický průměr vážený geometrický průměr (wi ) K1 3.20109 0.56774 K2 0.79527 0.14105 K3 1.35120 0.23965 K4 0.29072 0.05156 suma 5.638272125 1 Metoda postupného rozvrhu vah Při velkém počtu kritérií je vhodné seskupit kritéria do dílčích skupin podle příbuznosti jejich věcné náplně. Váhy jednotlivých kritérií pak určíme tak, že • stanovíme normované váhy jednotlivých skupin kritérií (pomocí některé z dříve uvedených metod) • stanovíme normované váhy každého kritéria v příslušné skupině • vynásobením vah skupin kritérií a vah jednotlivých kritérií v rámci každé skupiny zjistíme výsledné normované váhy kritérií Uvedený postup je ilustrován na následujícím příkladu: Výrobní firma uvažující o vytvoření společného podniku hledá partnera, který by co nejlépe vyhovoval těmto kritériím: K1 K2 K3 K4 K5 K6 K7 K8 K9
... ... ... ... ... ... ... ... ...
možnost kooperace v klíčových výrobních programech firmy kompatibilnost technologie technická úroveň partnera možnost dodávat zboží na trhy partnera možnost samostatně vystupovat na zahraničních trzích finanční výsledky hospodaření partnera v posledním roce velikost finančního vkladu partnera do společného podniku vzdálenost sídla partnera od sídla firmy image firmy partnera
Uvedený soubor kritérií lze podle jejich obsahové příbuznosti rozdělit do čtyř skupin: • kritéria výrobně-technologická (skupina S1 s kritérii K1 , K2 , K3 ) • kritéria obchodní (skupina S2 s kritérii K4 , K5 ) • kritéria finanční (skupina S3 s kritérii K6 , K7 ) • ostatní (skupina S4 s kritérii K8 , K9 ) Za předpokladu, že byly stanoveny váhy těchto skupin kritérií a váhy jednotlivých kritérií v rámci uvažovaných skupin, výše popsaným postupem byly spočítány a do posledního řádku následující tabulky zapsány výsledné váhy jednotlivých kritérií. Skupina kritérií Váhy skupin kritérií Kritéria Váhy kritérií Výsledn váhy kritérií
S1 S2 S3 S4 0,3 0,2 0,4 0,1 K1 K2 K3 K 4 K5 K 6 K7 K8 K9 0,5 0,2 0,3 0,5 0,5 0,3 0,7 0,4 0,6 0,15 0,06 0,09 0,1 0,1 0,12 0,28 0,04 0,06
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
9
Váhy kritérií patří k údajům subjektivního charakteru, závisejícím jednak na použití metody, jednak na hodnotiteli. Doporučuje se proto aplikovat více metod, zapojit více hodnotitelů a získané hodnoty průměrovat. Stejně jako lze kritéria rozdělit do skupin a každé skupině přidělit váhu, může hodnocení provádět více hodnotitelů, každému hodnotiteli je pak přidělena váha, kterou hodnotitel dělí mezi kritéria.
1.3
Metody stanovení pořadí variant
Cílem metod vícekriteriálního hodnocení variant je stanovení pořadí výhodnosti jednotlivých variant z hlediska zvolených kritérií, přičemž varianta s nejlepším umístěním představuje nejlepší kompromisní variantu. Metody pro výběr kompromisní varianty mezi nedominovanými variantami se liší přístupem k pojmu ”kompromisní varianta”, náročností a použitelností pro různé typy vícekriteriálních úloh. Výsledky získané různými metodami mají tedy subjektivní charakter a mohou se navzájem lišit. Metody je možné rozdělit podle toho, jaký typ informace vyžadují. • Metody vyžadující aspirační úrovně kriteriálních hodnot Do této skupina metod patří například konjunktivní metoda, disjunktivní metoda a metoda PRIAM. Informace o důležitosti kritérií je vyjádřena aspirační úrovní kritérií. Porovnávají se kriteriální hodnoty všech variant s aspiračními úrovněmi všech kritérií. Obvykle se rozdělí skupina variant na dvě skupiny. Varianty, které mají horší kriteriální hodnoty, než je nastavená aspirační úroveň (neakceptovatelné, neefektivní) a varianty, které mají lepší nebo stejné kriteriální hodnoty, než je aspirační úroveň (akceptovatelné, efektivní). Při dostatečném zpřísnění aspiračních úrovní může v množině akceptovatelných variant zůstat varianta jediná, kterou označíme jako kompromisní. • Metody vyžadující ordinální informace o variantách podle každého kritéria Jsou to například metoda pořadí, lexikografická metoda, permutační metoda, metoda ORESTE • Metody vyžadující kardinální informace o variantách podle každého kritéria Tato skupina metod se dále rozděluje na dílčí podskupiny podle principu, na kterém jsou hodnocení založena. Existují tyto základní přístupy: – maximalizace užitku (metoda váženého součtu, metoda bázické varianty, metoda AHP, metoda bodovací) – minimalizace vzdálenosti od ideální varianty, popř. maximalizace vzdálenosti od bazální varianty (TOPSIS) – preferenční relace (ELECTRE, PROMETHEE) – metody založené na mezní míře substituce (metoda postupné substituce) Některé zmíněné metody budou vysvětleny dále, jiné popisují např. Fiala (1994) nebo Brožová (2003). Konjunktivní a disjunktivní metoda Při aplikaci těchto metod je nutné, aby byly známé aspirační úrovně všech kritérií a kardinální ohodnocení variant podle jednotlivých kritérií. Podle aspirační úrovně rozdělíme varianty na akceptovatelné a neakceptovatelné. V případě konjunktivní metody připustíme pouze varianty, které splňují všechny aspirační úrovně. Nejprve předpokládejme, že všechna kritéria jsou maximalizační. Ze všech možných variant - alternativ Ai budou pak v množině akceptovatelných variant M ty varianty, které podle všech posuzovaných hledisek mají alespoň hodnotu předem stanovené aspirační úrovně zj nebo hodnotu větší. Matematicky můžeme zapsat
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
M = { Ai |yij ≥ zj , ∀j ∈ 1, 2, . . . , n}.
10
(1.9)
Nyní předpokládejme, že jsou všechna kritéria minimalizační. Ze všech možných variant - alternativ Ai budou pak v množině akceptovatelných variant M ty varianty, které podle všech posuzovaných hledisek mají nejvýše hodnotu předem stanovené aspirační úrovně zj nebo hodnotu menší. Matematicky můžeme zapsat M = { Ai |yij ≤ zj , ∀j ∈ 1, 2, . . . , n}.
(1.10)
Běžně se v úlohách vyskytují jak maximalizační, tak minimalizační kritéria. Potom můžeme říci, že do množiny M (množiny akceptovatelných variant) patří pouze varianty, které podle všech kritérií dosahují předem stanovenou aspirační úroveň nebo jsou ještě lepší. V případě disjunktivní metody připustíme varianty, které splňují alespoň jeden požadavek. V případě maximalizačních kritérií platí, že ze všech možných variant - alternativ Ai budou pak v množině akceptovatelných variant M ty varianty, které alespoň podle jednoho posuzovaného hlediska mají alespoň hodnotu předem stanovené aspirační úrovně zj nebo hodnotu větší. Matematicky můžeme zapsat M = { Ai |yij ≥ zj , ∃j ∈ 1, 2, . . . , n}, (1.11) V případě minimalizačních kritérií platí, že ze všech možných variant - alternativ Ai budou pak v množině akceptovatelných variant M ty varianty, které alespoň podle jednoho posuzovaného hlediska mají alespoň hodnotu předem stanovené aspirační úrovně zj nebo hodnotu menší. Matematicky můžeme zapsat M = { Ai |yij ≤ zj , ∃j ∈ 1, 2, . . . , n}, (1.12) Pokud jsou požadavky vyjádřeny aspiračními úrovněmi příliš přísné, je množina akceptovatelných variant prázdná. V takovém případě je nutno zadat nové, mírnější aspirační úrovně. A naopak, budou-li požadavky mírné, množina variant bude příliš velká. Pak je nutné aspirační úrovně zpřísnit. Řešený příklad 6. Vezměme u našeho příkladu 1 tyto aspirační úrovně: z1 = 25, z2 = 45, z3 = 2, z4 = 8. Určete, které varianty jsou akceptovatelné a které nikoliv. Řešení. Pokud problém budeme řešit konjunktivní metodou, jediná varianta, která vyhovuje všem aspiračním úrovním, je varianta C. Pokud k řešení přistoupíme metodou disjunktivní, pak alespoň z jednoho hlediska vyhovují všechny varianty. 2 Metoda PRIAM (Programme utilisatnt lIntelligence Artificiele en Multicritere) Tato metoda je založena na postupném prohledávání množiny variant v s krocích, aby bylo nalezeno jediné nedominované řešení. Požadované informace jsou ohodnocení variant podle jednotlivých kritérií. Každá varianta Ai je zobrazena vektorem kriteriálních hodnot yi ∈ Y (jeden z řádků kriteriální matice). Aspirační úroveň j-tého kritéria v s-tém kroku je označena zjs a změny aspirační úrovně j-tého v s-tém kroku ∆zjs . Rozhodovatel navrhne první aspirační úroveň kritérií (0)
(0)
z(0) = (z1 , z2 , . . . , zn(0) ).
(1.13)
Zpravidla se aspirační úroveň v nultém kroku stanoví jako nejhorší hodnota podle každého kritéria. Tímto ”sítem” projdou všechny varianty. Nebo-li pro všechny varianty (jejich kriteriální hodnoty) platí yi ≥ z(0) . (1.14)
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
11
Obecně do dalšího kola projdou varianty, které splňují požadované aspirační úrovně. Počet variant splňující tyto aspirační úrovně udává číslo d. Vzhledem k hodnotě d rozhodovatel mění aspirační úroveň kritérií pro krok s + 1 z(s+1) = z(s) + ∆z(s) . (1.15) Podle hodnoty čísla d (počet variant, které splňují požadavky) nastávají tři případy: • d > 1 - rozhodovatel mění aspirační úroveň tak, aby snížil počet akceptovatelných variant • d = 1 - je nalezena kompromisní varianta • d = 0 - neexistuje žádná přijatelná varianta, hledá se nejbližší varianta k zadaným aspiračním úrovním; v takovém případě se pro každou variantu vypočte odchylka od aspirační úrovně podle vztahu 1.16 (s) n X |zj − yij | , (1.16) yj∗ j=1 kde yj∗ pro j = 1, 2, . . . , n jsou ideální kriteriální hodnoty. Tento vztah nelze použít v případě, že ideální kriteriální hodnota je nula. Jako přijatelnou variantu vybereme variantu s nejmenším podílem odchylky od aspiračních úrovní kritérií na ideální kriteriální hodnotě. Řešený příklad 7. Vyberte kompromisní variantu v příkladu 1 metodou PRIAM. Řešení. Výchozí aspirační úroveň u našeho příkladu je z(0) = (22; 60; 1; 7.5). Této aspirační úrovni vyhovují všechny varianty. Zpřísníme tedy aspirační úroveň nejprve u nejdůležitějšího kritéria o 3. Nebo-li δz(1) = (3; 0; 0; 0) a aspirační úrovně jsou pak z(1) = (25; 60; 1; 7.5). Tím je vyřazena varianta B, zbývají stále dvě varianty. Zpřísníme aspirační úroveň například u kritéria odborný růst o 1, tedy δz(2) = (0; 0; 1; 0). Aspirační úroveň bude po této úpravě z(2) = (25; 60; 2; 7.5). Stále zbývají dvě varianty, které stanovené aspirační úrovně splňují. Proto zpřísníme aspirační úroveň u druhého kritéria δz(3) = (0; 10; 0; 0). Po této úpravě je aspirační úroveň z(3) = (25; 50; 2; 7.5) a tomu vyhovuje už jen varianta C. 2 Poznámka. Metoda PRIAM a metody konjunktivní a disjunktivní jsou doporučovány zejména k ”předvýběru” variant, které se pak dále hodnotí jinými metodami. Dalo by se říci, že se těmtito metodami odstranily vyloženě špatné a nevyhovující varianty a ostatní se pak dále vyhodnocují. Metoda pořadí Metoda pořadí je založena na převedení kriteriální matice na matici pořadí. To znamená, že postupně se podle všech kritérií přiřadí variantám jejich pořadí. Pokud nejsou známé preference kritérií, pouze se sečtou pro každou variantu všechna pořadí. Nejlepší varianta má tento součet nejnižší. Pokud jsou známé preference kritérií (váhy), lze vypočítat vážené pořadí variant, opět nejlepší varianta má tento součet nejnižší. Pro náš příklad je matice pořadí bez vah následující:
A B C S váhami pak:
K1 1 3 2
K2 3 1 2
K3 2 3 1
K4 1 3 2
Součet pořadí 7 10 7
Pořadí 1. 2. 1.
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
12
K1 K2 K 3 K4 Vážený součet pořadí Pořadí A 0.5 0.6 0.5 0.05 1.65 1. B 1.5 0.2 0.75 0.15 2.6 3. C 1 0.4 0.25 0.1 1.75 2. váhy 0.5 0.2 0.25 0.05 Metoda lexikografická Tato metoda vychází z předpokladu, že největší vliv na výběr kompromisní varianty má nejdůležitější kritérium. V případě, že existuje více variant se stejným hodnocením podle nejdůležitějšího kritéria, přichází v úvahu druhé nejdůležitější kritérium. Algoritmus se zastaví ve chvíli, kdy je vybraná jediná varianta, nebo když jsou vyčerpána všechna uvažovaná kritéria. Kompromisní varianty jsou potom všechny ty, které zůstaly stejně hodnoceny po zařazení posledního kritéria. Tato metoda by se dala přirovnat k vyhledávání ve slovníku. Příklad 3. V našem příkladu 1 si zvolíme jako nejdůležitější např. kritérium mzda. Z tohoto hlediska je nejlepší varianta A, proto bychom doporučili vybrat tuto variantu. V případě, že by mzda byla stejná u více variant, pak bychom přistoupili k druhému nejdůležitějšímu kritériu a podle něj bychom rozhodli o výběru vhodné varianty. 4 Metoda bodovací Při této metodě rozhodovatel přiřadí každému prvku rozhodovací matice určitý počet bodů ze zvolené stupnice, a to tak, že lepší hodnotě kritéria přiřadí větší počet bodů. Maximálně (minimálně) možný počet bodů přiřazený nejlepší (nejhorší) hodnotě kritéria musí být pro všechna kritéria stejný, přičemž může jít o hypoteticky stanovená čísla, která se v žádné variantě nevyskytují. Výhodné je opatřit bodovou stupnici pro každé kritérium slovním popisem. Řešený příklad 8. Vyhodnoťte varianty z příkladu 1 bodovací metodou. Řešení. Bodovací stupnice, kterou jsme navrhli pro náš příklad je následující: Body 1 2 3
K1 K2 K3 K4 méně než 25 tis. Kč nad 50 minut malá možnost do 8 hodin včetně h25; 28) (40; 50i střední možnost (8; 9i 28 tis. Kč a více 40 minut a méně velká možnost po 9. hodině
Na základě bodovací stupnice jsme každé kriteriální hodnotě přiřadili příslušný počet bodů a potom jsme opět (po pronásobení vahami) sečetli pro každou variantu body. Nejlepší je varianta s nejvyšším součtem. K1 K2 K3 K4 Body Pořadí A 3 1 2 2 2.3 1. B 1 3 1 1 1.4 3. C 2 2 3 1 2.2 2. Váhy 0.5 0.2 0.25 0.05 Poznámka. Čím širší si zvolíme bodovací stupnici, tím přesnější můžeme získat výsledky. Metoda váženého součtu Při vícekriteriálním hodnocení variant můžeme každé hodnotě kritéria Kj přiřadit její užitek, tedy můžeme vytvořit dílčí užitkovou funkci uj , která pro variantu Ai nabývá hodnoty uj Ai = uij ; i = 1, 2, . . . , m; j = 1, 2, . . . , n.
(1.17)
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
13
Definičním oborem této funkce je interval mezi nejlepší a nejhorší hodnotou příslušného kritéria. Oborem funkčních hodnot je interval h0, 1i. Tato metoda je vhodná především pro kvantitativní kritéria. Předpokládá lineární závislost užitku na hodnotách kritéria, přičemž nejhorší hodnotě j -tého kritéria (budeme značit dj ) přiřadíme hodnotu 0 a nejlepší hodnotě (budeme značit hj ) užitek 1. Pro dílčí užitek uij hodnoty yij platí uij =
yij − dj ; i = 1, 2, . . . , m; j = 1, 2, . . . , n. hj − dj
(1.18)
Pro jednotlivé varianty vypočteme agregovanou funkci užitku podle vztahu u(Ai ) =
n X
wj uij
(1.19)
j=1
Varianty pak seřadíme podle hodnot u(Ai ). Nejlepší varianta má tuto hodnotu největší. Řešený příklad 9. Vyhodnoťte metodou váženého součtu varianty z příkladu 1. Řešení. Nejprve tedy vybereme pro každé kritérium nejlepší (hj ) a nejhorší (dj ) hodnotu. Varianty K1 A 30 B 22 C 26 hj 30 dj 22
K2 60 30 45 30 60
K3 2 1 3 3 1
K4 9 7.5 8 9 7.5
Nyní přepočteme kriteriální hodnoty podle vztahu (1.18), potom vypočteme agregovanou funkci užitku pro každou variantu podle vztahu (1.19) a varianty seřadíme podle této hodnoty od nejlepší po nejhorší. Použité váhy jsou váhy získané bodovací metodou. K1 K2 K3 K4 u(Ai ) pořadí A 1 0 0.5 1 0.675 1. B 0 1 0 0 0.2 3. C 0.5 0.5 1 0.333 0.616666667 2. váhy 0.5 0.2 0.25 0.05 2 Metoda bázické varianty Za bazickou variantu je považována varianta, která dosahuje nejlepších či předem stanovených hodnot z hlediska všech kritérií. Vytvoření užitkové funkce s využitím bazické varianty spočívá v porovnávání hodnot důsledků jednotlivých variant s odpovídajícími hodnotami v bazické variantě. Označíme-li (b) yj hodnotu j-tého kritéria v bazické variantě, pro užitek kritéria výnosového typu při volbě i-té varianty platí uij =
yij (b)
; i = 1, 2, . . . , m; j = 1, 2, . . . , n.
(1.20)
yj
a u kritéria nákladového typu je dílčí užitek dán vztahem (b)
yj uij = ; i = 1, 2, . . . , m; j = 1, 2, . . . , n. yij
(1.21)
Pro jednotlivé varianty opět spočítáme agregované funkce užitku a podle jejich hodnot varianty seřadíme.
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
14
Poznámka. Pozor na nulové kriteriální hodnoty u bázické varianty! V případě, že by bázická varianta měla nulovou hodnotu, tato metoda nelze použít. Řešený příklad 10. Vyhodnoťte varianty z příkladu ?? metodou bázické varianty. Předpokládejte, že bázická varianta je varianta ideální, to znamená, že podle každého kritéria dosahuje nejlepšího možné výsledku. Řešení. Nejprve opět přepočteme podle vztahů (1.20) a (1.21) kriteriální matici, a stejně jako u metody WSA spočítáme hodnoty agregovaného užitku a varianty seřadíme podle těchto hodnot sestupně. K1 K2 K3 K4 u(Ai ) pořadí A 1 0.5 0.0.667 1 0.817 2. B 0.733 1 0.333 0.833 0.692 3. C 0.867 0.667 1 0.889 0.861 1. váhy 0.5 0.2 0.25 0.05 AHP Metoda AHP (Analytic Hierarchy Process) byla navržena prof. Saatym v roce 1980. Při řešení rozhodovacích problémů je třeba brát v úvahu všechny prvky, které ovlivňují výsledek analýzy, vazby mezi nimi a intenzitu, s jakou na sebe vzájemně působí. Rozhodovací problém lze znázornit jako hierarchickou strukturu. Je to lineární struktura obsahující s-úrovní, přičemž každá z těchto úrovní zahrnuje několik prvků. Uspořádání jednotlivých úrovní je vždy od obecného ke konkrétnímu. Pro obecnou úlohu vícekriteriálního hodnocení variant může být hierarchie následující: 1. úroveň - cíl vyjednávání 2. úroveň - experti, kteří se na hodnocení podílí 3. úroveň - kritéria vyhodnocování 4. úroveň - posuzované varianty 1.1
Obrázek 1.1: Hierarchická struktura úlohy vícekriteriálního rozhodování
Obdobným způsobem, jako mezi kritérii při určování vah kritérií Saatyho metodou, lze určit vztahy mezi všemi komponentami na každé úrovni hierarchie. Pokud máme čtyřúrovňovou hierarchii, tzn. jeden cíl, h expertů, n kritérií a m variant, bude na druhé úrovni hierarchie jedna matice párového srovnávání o rozměrech h × h. Na třetí úrovni bude h matic o rozměrech n × n a na čtvrté úrovni n matic o rozměrech m × m . Pomocí propočtů (viz Saatyho metoda pro výpočet vah kritérií) v těchto maticích si varianty ”rozdělují” hodnotu váhy příslušného kritéria (kritéria si pak ”rozdělují” váhy příslušného experta). Hodnoty, které získáme se nazývají preferenční indexy variant z hlediska všech kritérií. Pokud tedy sečteme tyto preferenční indexy z hlediska všech kritérií, získáme hodnocení varianty z pohledu všech expertů a z hlediska všech kritérií.
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
15
Pro ilustraci této metody použijeme váhy získané Saatyho metodou. Předpokládáme, že do hodnocení je zapojen pouze jeden expert, nebo-li jeden hodnotitel porovnává kritéria mezi sebou (tím získá váhy) a potom porovnává jednotlivé varianty mezi sebou podle jednotlivých kritérií. Všechna porovnání postupně pro plat, dojíždění, odborný růst a začátek pracovní doby jsou zaznamenána v následujících tabulkách: Kritéria K1 K2 K3 K4 K1 1 5 3 7 K2 1/5 1 1/2 4 K3 1/3 2 1 5 K4 1/7 1/4 1/5 1
Plat A B A 1 7 B 1/7 1 C 1/5 2
Dojíždění A B C
A 1 7 3
Růst A A 1 B 7 C 3
Začátek A A 1 B 7 C 3
Geometrický průměr Vážený geometrický průměr 3.20109 0.56774 0.79527 0.14105 1.35120 0.23965 0.29072 0.05156
C Geometrický průměr Vážený geometrický průměr 5 3.27107 0.73959 1/2 0.41491 0.09381 1 0.73681 0.16659
B C Geometrický průměr Vážený geometrický průměr 1/7 1/3 0.36246 0.08795 1 3 2.75892 0.66942 1/3 1 1 0.24264
B C Geometrický průměr Vážený geometrický průměr 1/7 1/3 1 0.24986 1 3 0.38157 0.09534 1/3 1 2.62074 0.65481
B C Geometrický průměr Vážený geometrický průměr 1/7 1/3 2.62074 0.65481 1 3 0.38157 0.09534 1/3 1 1 0.24986
Kritéria jsou celkem čtyři a každé z nich má svou váhu. Tato váha musí být dále rozdělena mezi jednotlivé varianty. My známe váhy jednotlivých kritérií a známe také váhy variant podle těchto kritérií. Výsledné váhy každé varianty podle každého kritéria jsou v další tabulce.
A B C Váhy kritérií
K1 0.73959 0.09381 0.16659 0.56774
K2 0.36246 0.66942 0.24264 0.14105
K3 K4 Součet hodnocení Pořadí 0.24986 0.6548 0.525943481 1. 0.09534 0.09534 0.175444971 3. 0.65481 0.24986 0.298611548 2. 0.23965 0.05156
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
16
Pro každou variantu sčítáme hodnocení podle kritéria vynásobené váhou tohoto kritéria. Součty dílčích hodnocení jsou seřazeny sestupně. Varianty jsou seřazeny A C B. TOPSIS(Technique for Order Preference by Similarity to Ideal Solution) Metoda TOPSIS je založena na výběru varianty, která je nejblíže k ideální variantě a nejdále od bazální varianty. Předpokládá se maximalizační charakter všech kritérií. Pokud nejsou všechna kritéria maximalizační, je nutné na maximalizační převést. Postup při metodě TOPSIS lze popsat v následujících krocích: • převedení všech kritérií na maximalizační • vytvoření normalizované kriteriální matice R = rij podle vztahu yij rij = r m P
; i = 1, 2, . . . , m; j = 1, 2, . . . , n.
(1.22)
yij2
i=1
Sloupce v matici R představují vektory jednotkové délky. • Dále se převede kriteriální matice R na normalizovanou kriteriální matici Z tak, že každý sloupec matice R vynásobíme vahou odpovídajícího kritéria podle vztahu zij = wj rij
(1.23)
• pomocí prvků matice Z se vytvoří ideální varianta (h1 , h2 , . . . , hn ) a bazální varianta (d1 , d2 , . . . , dn ), kde hj = max zij ; j = 1, 2, . . . , n (1.24) i
dj = min zij ; j = 1, 2, . . . , n i
• Vzdálenost od ideální varianty se počítá podle vztahu v uX u n + di = t (zij − hj )2 ; i = 1, 2, . . . , m
(1.25)
(1.26)
j=1
Vzdálenost od bazální varianty se počítá podle vztahu v uX u n − di = t (zij − hj )2 ; i = 1, 2, . . . , m
(1.27)
j=1
• Relativní ukazatel vzdálenosti variant od bazální varianty se vypočte podle vztahu ci =
d− i ; i = 1, 2, . . . , m d+ + d− i i
Varianty jsou uspořádány podle nerostoucích hodnot ci Řešený příklad 11. Vyhodnoťe metodou TOPSIS varianty z příkladu 1. Řešení. Všechna kritéria převedeme na maximalizační typ:
A B C
K1 30 22 26
K2 0 30 15
K3 2 1 3
K4 9 7.5 8
(1.28)
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
17
Vytvoříme normalizovanou kriteriální matici podle vztahu (1.22):
A B C
K1 K2 K3 K4 0.660978974 0 0.534522484 0.634416639 0.484717914 0.894427191 0.267261242 0.528680533 0.572848444 0.447213595 0.801783726 0.563925901
Kriteriální matici znormujeme váhami získanými bodovací metodou podle vztahu (1.23) a pak vybereme pro každé kritérium nejvyšší (ideální) a nejnižší (bazální) hodnotu:
A B C hj dj
K1 K2 K3 K4 0.330489487 0 0.133630621 0.031720832 0.242358957 0.178885438 0.06681531 0.026434027 0.286424222 0.089442719 0.200445931 0.028196295 0.330489487 0.178885438 0.200445931 0.031720832 0.242358957 0 0.06681531 0.026434027
Podle vztahu (1.26) vypočteme vzdálenost od ideální varianty, podle (1.27) vzdálenost od bazální varianty a podle (1.28) relativní vzdálenost od bazální varianty. Varianty se řadí sestupně podle relativní vzdálenosti od bazální varianty.
A B C
d+ d− ci Pořadí i i 0.19095624 0.110721391 0.367018894 3. 0.160162678 0.178885438 0.527610771 2. 0.099770587 0.166739306 0.625640213 1. 2
Metoda ELECTRE I Cílem této metody je rozdělit množinu variant na dvě indiferenční třídy, na efektivní a neefektivní varianty. Předpokladem pro využití této metody je znalost kriteriální matice, vektoru normalizovaných vah a stanovení dvou prahových hodnot, a to prahu preference a prahu dispreference. Ohodnocení varianty Ai podle kritéria j označíme symbolem yij ; i = 1, 2, . . . , m; j = 1, 2, . . . , n. Pro každou dvojici variant Ai , Ah ; i, h = 1, 2, . . . , m pak určíme množinu Cih = { j|yij ≥ yhj ; j = 1, 2, . . . , n}; i, h = 1, 2, . . . , m
(1.29)
obsahující indexy kritérií, z jejichž hlediska je varianta Ai ohodnocena alespoň tak dobře jako varianta Ah a množinu Dih = { j|yij < yhj ; j = 1, 2, . . . , n}; i, h = 1, 2, . . . , m, (1.30) která obsahuje indexy zbývajících kritérií, ve kterých je varianta Ai horší než varianta Ah . Na základě množiny Cih a normalizovaného vektoru vah pro každou dvojici variant Ai , Ah určíme číslo cih představující součet vah kritérií, z jejichž hlediska je varianta Ai stejně tak dobrá nebo lepší než varianta Ah . V případě, že Cih není prázdná množina, X cih = wj ; i, h = 1, 2, . . . , m, (1.31) j∈Cih
pokud je Cih prázdná množina, pak cih = 0.
(1.32)
Hodnota cih představuje stupeň preference varianty Ai před variantou Ah . Platí, že cih ∈ h0; 1i.
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
18
Dále je nutné vypočítat pro každou dvojici Ai , Ah stupeň dispreference dij . Pokud je Dih prázdná množina, pak dih = 0, v opačném případě max (yij − yhj )
dih =
j∈Dih
max(yij − yhj )
; i, h = 1, 2, . . . , m.
(1.33)
h
Rovněž dih ∈ h0; 1i. Pro určení celkové preference P mezi dvojicí variant musí rozhodovatel zadat práh preference c∗ a práh dispreference d∗ . Platí, že varianta Ai je preferována před variantou Ah tehdy, když platí cih ≥ c∗ ∧ dih ≤ d∗ .
(1.34)
Párové preference jsou pak zapsány v matici P = (pih ), kde pih = 1 pokud Ai je preferována před Ah , v opačném případě je pih = 0. Varianty jsou rozděleny na efektivní a neefektivní. Za efektivní jsou považovány varianty, které jsou preferovány alespoň před jednou variantou a neexistuje k nim žádná preferující varianta. Výsledek závisí na stanovených hodnotách prahů preference a dispreference. Doporučuje se vycházet z průměrných hodnot v maticích preference a dispreference a potom je postupně buď zpřísňovat (v případě mnoha efektivních variant) a nebo zmírňovat (v případě žádné efektivní varianty). Poznámka. Ve výše popsaném postupu vyhodnocování variant metodou ELECTRE porovnáváme mezi sebou kriteriální hodnoty v různých měřítkách (peníze, vzdálenost, body, atd.), což není úplně vpořádku, protože měřítko má vliv na výsledky. Lepší by bylo nejprve kriteriální matici znormovat a teprve potom dělat párové srovnávání. Postup výpočtu s nenormovanou maticí jsme použili proto, že software, který používáme při výuce rovněž počítá s nenormovanou maticí. Nejprve je nutné pro každou dvojici zjistit množinu C (viz (1.29))a množinu D (viz (1.30)), například pro dvojici A, B je množina CA,B =K1 , K3 , K4 , množina DA,B =K2 . Varianta A je lepší podle kritérií 1, 3 a 4, horší je podle kritéria 2. Naopak budeme-li porovnávat B, A, pak CB,A =K2 , množina DB,A =K1 , K3 , K4 . Varianta B je lepší než A podle kritéria 2 a horší podle kritérií 1, 3, 4. Dále podle vztahů (1.31) a (1.33) sestavíme matici preferencí C a dispreferencí D. − 0.8 0.55 C = 0.2 − 0.2 0.45 0.8 − Matice je v našem příkladu typu 3 × 3, protože mezi sebou porovnáváme tři varianty. Výsledek porovnání zapisujeme v matici na pozici: první z dvojice porovnávaných hledáme v řádku, druhou z dvojice porovnávaných ve sloupci. Předpokládáme pořadí variant A, B, C. Proto například pro dvojici A, B píšeme výsledek porovnávání na pozici 1, 2 (1. řádek, 2. sloupec). Pro dvojici A, B je pak cA,B = 0.5 + 0.25 + 0.05 = 0.8, pro dvojici B, A je cB,A = 0.2
− 1 1 D = 0.267 − 0.267 0.267 1 − 8 30 = 1, pro dvojici B, A je dB,A = 30 = 0.267 Pro dvojici A, B je pak dA,B = 30 ∗ Nyní můžeme zadat práh preference (např. c = 0.4), práh dispreference (např. d∗ = 0.6) a na základě vztahu (1.34) sestavit matici preferencí:
− 0 0 P= 0 − 0 1 0 −
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
19
Z matice můžeme vyčíst, že pouze varianta C je efektivní, je preferována alespoň před jednou variantou a neexistuje k nim žádná preferující varianta.V matici P je pouze u varianty C v řádku jenička a zároveň ve sloupci samé nuly.
1.4
Analýza citlivosti preferenčního pořadí variant
Pořadí výhodnosti daných variant rozhodování, stanovené některou z metod vícekriteriálního hodnocení variant, závisí především na vahách jednotlivých kritérií a na použité metodě. Zkoumání citlivosti preferenčního pořadí variant na stanovení důležitosti jednotlivých kritérií spadá do oblasti experimentování na modelech, kdy vícekriteriální hodnocení variant provádíme při měnících se vahách kritérií. Jedině tehdy, kdy při těchto změnách se preferenční pořadí variant (nebo aspoň varianta s nejvyšším oceněním) nemění, není toto pořadí citlivé na nepřesnost stanovených vah a lze říci, že význam jednotlivých kritérií byl rozhodovatelem správně oceněn. Jestliže preferenční pořadí variant je značně citlivé na změny vah kritérií, je nutné jejich spolehlivost zvýšit. Závislost preferenčního pořadí variant na použité metodě jejich vícekriteriálního hodnocení je dána tím, že jednotlivé metody vycházejí z různých, zpravidla zjednodušujících předpokladů a že v důsledku různého přístupu k pojmu „kompromisní variantaÿ využívají různé výpočetní postupy. Doporučuje se proto při vícekriteriálním hodnocení variant uplatnit více metod a ověřit citlivost preferenčního pořadí variant vzhledem k použitým metodám. Jedině tu variantu, která zůstává stále na 1. místě při použití libovolné metody, lze považovat za nejvýhodnější. Vliv použit metody vícekriteriálního hodnocení variant na jejich preferenční pořadí lze dokumentovat na příkladu. V následující tabulce je uvedeno pořadí jednotlivých pracovišť z hlediska daných kritérií, stanovené různými metodami: Metoda Pořadí Bodovací WSA A 1. 1. 1. B 3. 3. 3. C 2. 2. 2.
Bázická AHP 2. 1. 3. 3. 1. 2.
TOPSIS 3. 2. 1.
Z uvedené tabulky je patrná značná citlivost preferenčního pořadí variant na aplikované metodě. Nicméně varianta A je na prvním místě nejčastěji, na třetím místě je pouze při použití metody TOPSIS, proto bychom ji označili jako kompromisní variantu. Jiná možnost využití výsledků získaných různými metodami spočívá v tom, že se na pořadí pracovišť můžeme dívat jako na prvky kriteriální matice a znovu můžeme použít některou z metod vícekriteriálního hodnocení variant. Příklad 4. Vezměme pořadí variant získaná různými metodami jako prvky kriteriální matice a vyhodnoťme je metodou pořadí. Nebudeme uvažovat váhy, takže pouze sečteme pořadí v řádcích. Nejlépe vyjde varianta A, u které je součet pořadí 9, potom je varianta C se součtem 10 a nejhůře dopadla varianta B se součtem 17. 4 Před volbou varianty určené k realizaci je nutné vzít v úvahu také možný vliv rizika a nejistoty na výsledky vícekriteriálního hodnocení variant. Dosud jsme předpokládali, že důsledky variant vzhledem k jednotlivým kritériím, tj. prvky kriteriální matice, jsou jednoznačně určeny. Jestliže tyto důsledky mají náhodný charakter, buď opět pomocí analýzy citlivosti zkoumáme vliv možných změn prvků kriteriální matice na preferenční pořadí variant, nebo použijeme některou z metod vícekriteriálního hodnocení variant za rizika a nejistoty. Tyto metody popisuje např. Fotr 1992.
1.5
PŘÍPADOVÁ STUDIE - výběr osobního automobilu
Rodina stojí před rozhodnutím, jaké auto si má pořídit. Rozhoduje se podle sedmi kritérií: K1 cena (tis.Kč ), K2 výkon (kW ), K3 spotřeba (l ), K4 zavazadlový prostor (l ), K5 maximální rychlost
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
20
(km/hod ), K6 poruchovost (%) a K7 vzhled. Kritérium vzhled je kvalitativní kritérium, je ohodnoceno slovně, a to velmi spokojen, spokojen spokojen s výhradami a nespokojen. Toto slovní hodnocení je převedeno na kvantitativní (body). Poruchovost automobilů je čerpána z údajů TÜV pro automobily staré dva až tři roky. Aspirační úrovně kritérií, která si rodina stanovila, splňují tyto automobily: Škoda Octavia, Toyota Avensis, Peugeot, Ford Focus a Renault Mégane. Hodnoty jednotlivých kritérií jsou uvedeny v následující tabulce. Typ auta Škoda Octavia Renault Mégane Peugeot 407 Toyota Avensis Ford Focus
K1 K2 K3 K4 K5 K6 K7 658.4 103 5.5 1620 208 5.5 2 646.6 96 5.7 1600 200 5.9 1 784 100 5.6 1365 189 11.2 4 699.9 93 5.7 1500 200 4.2 3 717.9 100 5.7 1525 203 5 2
STANOVENÍ VAH KRITÉRIÍ Nejprve je nutné stanovit si, jak jsou která kritéria preferovaná a jakou má které kritérium váhu. Preference si ”hlava rodiny” stanovila takto: K6 K4 K3 K1 K2 K7 K5 . Na základě těchto preferencí lze spočítat váhy různými metodami. Metoda pořadí Kritérium K1 K2 K3 K4 K5 K6 K7 Pořadí 4 5 3 2 7 1 6 Body 4 3 5 6 1 7 2 Váhy 0.143 0.107 0.179 0.214 0.036 0.25 0.071 Metoda bodovací Kritérium K1 K2 K3 K4 K5 K6 K7 Body 15 10 20 23 2 25 5 Váhy 0.15 0.1 0.2 0.23 0.02 0.25 0.05 Fullerova metoda I u této metody dodržujeme stejné preferenční uspořádání, jako u předchozích metod, to je K6 K4 K3 K1 K 2 K7 K5 . Kritérium Váhy neupravené Váhy upravené
K1 K2 K3 K4 K5 K6 K7 0.143 0.095 0.190 0.238 0 0.286 0.048 0.143 0.107 0.179 0.214 0.036 0.25 0.071
Metoda bodovací - více hodnotitelů Vzhledem k tomu, že se jedná o rodinné auto, můžeme dát možnost všem členům rodiny vyjádřit svůj názor a převést ho na váhy. Rodina má čtyři členy, každý z nich má k dispozici například 100 bodů, které rozdělí podle svého uvážení mezi kritéria, viz následující tabulka:
Hodnotitel Hodnotitel Hodnotitel Hodnotitel
1 2 3 4
K1 20 10 10 15
K2 10 30 15 10
K3 20 10 15 20
K4 20 10 15 23
K5 5 10 20 2
K6 20 5 15 25
K7 5 15 20 5
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY
21
Váhy se pak vypočítají stejně, jako u prosté bodovací metody. Hodnotitel Hodnotitel Hodnotitel Hodnotitel
1 2 3 4
K1 0.20 0.10 0.10 0.15
K2 0.10 0.30 0.15 0.10
K3 0.20 0.10 0.15 0.20
K4 0.20 0.10 0.15 0.23
K5 0.05 0.10 0.20 0.02
K6 0.20 0.05 0.15 0.25
K7 součet 0.05 1 0.15 1 0.20 1 0.05 1
Pokud by například rozhodoval jen první hodnotitel, váhy kritérií bychom vyčetli v prvním řádku tabulky a použili bychom je k dalším výpočtům. Součet vah všech kritérií je rovný jedné. Vzhledem k tomu, že hodnotitelé jsou čtyři a každý z nich má stejnou váhu (w = 0.25), musíme váhy kritérií ještě přepočítat (pronásobit váhou rozhodovatele). Výsledné váhy jsou pak v násĺedující tabulce:
Hodnotitel Hodnotitel Hodnotitel Hodnotitel Váhy
1 2 3 4
K1 K2 K3 K4 K5 0.05 0.025 0.05 0.05 0.0125 0.025 0.075 0.025 0.025 0.05 0.025 0.0375 0.0375 0.0375 0.025 0.0375 0.025 0.05 0.0575 0.005 0.1375 0.1625 0.1625 0.17 0.0925
K6 0.05 0.0125 0.0375 0.0625 0.1625
K7 Součet 0.0125 0.25 0.0375 0.25 0.05 0.25 0.0125 0.25 0.1125 1
HODNOCENÍ VARIANT Známe-li váhy jednotlivých kritérií, můžeme přistoupit k vlastnímu hodnocení automobilů. Pro další hodnocení použijeme váhy vypočtené posledním způsobem, t.j. bodovací metodou s více hodnotiteli. Auta byla vybrána tak, aby splňovala aspirační úrovně všech kritérií, všechny varianty - automobily jsou nedominované. Metoda pořadí Varianta Škoda Octavia Renault Mégane Peugeot 407 Toyota Avensis Ford Focus Váhy
K1 K2 K3 K4 K5 K6 K7 Vážené pořadí 2 1 1 1 1 3 3.5 1.74375 1 4 4 2 3 4 5 3.2675 5 2.5 2 5 5 5 1 3.65625 3 5 4 4 4 1 2 3.3125 4 2.5 4 3 2 2 3.5 3.02 0.1375 0.1625 0.1625 0.17 0.0925 0.1625 0.1125
Poznámka. V případě, že varianty nabývají podle jednoho kritéria stejnou hodnotu, použijeme v tabulce pořadí tzv. průměrné pořadí. Metodou pořadí byly automobily uspořádány takto: Škoda Octavia Ford Focus Renault Mégane Toyota Avensis Peugeot 407. Metoda váženého součtu - WSA Před vlastním hodnocením variant touto metodou musíme vybrat vždy nejlepší a nejhorší hodnotu pro každé kritérium. Varianta K1 K2 K3 K4 K5 K6 K7 Nejlepší hodnota kritéria (H) 646.6 103 5.5 1620 208 4.2 4 Nejhorší hodnota kritéria (D) 784 93 5.7 1365 189 11.2 1 Hodnocení metodou WSA:
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY Varianta Škoda Octavia Renault Mégane Peugeot 407 Toyota Avensis Ford Focus Váhy
22
K1 K2 K3 K4 K5 K6 K7 Vážený užitek 0.914 1 1 1 1 0.814 0.333 0.883 1 0.3 0 0.922 0.579 0.757 0 0.52 0 0.7 0.5 0 0 0 1 0.3075 0.612 0 0 0.529 0.579 1 0.667 0.465 0.481 0.7 0 0.627 0.737 0.886 0.333 0.536 0.1375 0.1625 0.1625 0.17 0.0925 0.1625 0.1125
Z tabulky vyplývá, že varianty jsou metodou WSA seřazeny takto: Škoda Octavia Ford Focus Renault Mégane Toyota Avensis Peugeot 407. Metoda bázické varianty Varianta Škoda Octavia Renault Mégane Peugeot 407 Toyota Avensis Ford Focus Váhy
K1 K2 K3 0.982 1 1 1 0.932 0.965 0.825 0.97 0.982 0.924 0.903 0.965 0.9 0.97 0.965 0.1375 0.1625 0.1625
K4 K5 K6 K7 Vážený užitek 1 1 0.764 0.5 0.903 0.988 0.962 0.712 0.25 0.846 0.843 0.909 0.375 1 0.831 0.926 0.962 1 0.75 0.924 0.941 0.976 0.84 0.25 0.881 0.17 0.0925 0.1625 0.1125
Podle metody bázické varianty je pořadí automobilů následující: Toyota Avensis Škoda Octavia Ford Focus Renault Mégane Peugeot 407. Bodovací metoda Nejprve je nutné určit, kolik bodů náleží kriteriálním hodnotám: Body 5 4 3 2 1
K1 K2 K3 h630 − 650) h102 − 105) h5 − 5.5) h650 − 670) h99 − 102) h5.5 − 5.7) h670 − 690) h96 − 99) h5.7 − 6) h690 − 710) h93 − 96) h6 − 6.5) h710 − 730) h90 − 93) h6.5 − 7)
K4 h1700 − 1800) h1600 − 1700) h1450 − 1600) h1350 − 1450) h1200 − 1350)
K5 K6 K7 h205 − 210) h4 − 5) 5 h202 − 205) h5 − 6) 4 h199 − 202) h6 − 7) 3 h190 − 199) h7 − 10) 2 h180 − 190) h10 − 12) 1
Původním kriteriálním hodnotám přiřadíme podle bodovací stupnice body, které potom pronásobíme váhami a sečteme pro každou variantu. Varianty jsou seřazeny sestupně. Body Škoda Octavia Renault Mégane Peugeot 407 Toyota Avensis Ford Focus Váhy
K1
K2
K3
K4
K5
K6
K7
4 5 4 4 5 4 2 5 3 3 4 3 4 1 1 4 4 2 1 1 4 2 2 3 3 3 5 3 1 4 3 3 4 4 2 0.1375 0.1625 0.1625 0.17 0.0925 0.1625 0.1125
Vážený bodů 4.03 3.3825 2.4825 3.025 3.03
součet
Bodovací metodou jsou automobily seřazeny: Škoda Octavia Renault Mégane Ford Focus Toyota Avensis Peugeot 407. Metoda TOPSIS
KAPITOLA 1. VÍCEKRITERIÁLNÍ ROZHODOVÁNÍ ZA JISTOTY Varianta Škoda Octavia Renault Mégane Peugeot 407 Toyota Avensis Ford Focus Váhy
K1
K2
K3
0.08 0.076 0.145 0.088 0.071 0 0 0.074 0.077 0.054 0.069 0 0.042 0.074 0 0.1375 0.1625 0.1625
K4
K5
K6
23 K7
0.081 0.043 0.076 0.039 0.08 0.041 0.071 0.02 0.068 0.039 0 0.077 0.075 0.041 0.093 0.058 0.076 0.042 0.083 0.039 0.17 0.0925 0.1625 0.1125
Relativní vzdálenost 0.811 0.418 0.386 0.432 0.377
Metodou TOPSIS jsme určili pořadí variant Škoda Octavia Toyota Avensis Renault Mégane Ford Focus Peugeot 407. Různými metodami jsme dospěli k různým výsledkům, viz následující tabulka. Metoda Škoda Octavia Renault Mégane Peugeot 407 Toyota Avensis Ford Focus
Pořadí Bodovací WSA 1 1 1 3 2 3 5 5 5 4 4 4 2 3 2
TOPSIS 1 3 4 2 5
ORESTE Bázické varianty 1 2 3 3 5 5 4 1 2 3
Ze souhrnu je patrné, že nejlepších výsledků podle daných kritérií dosahovala varianta první - Škoda Octavia. Renault Mégane, Ford Focus a Toyota Avensis dosahují v průměru velmi podobných výsledků, Peugeot 407 v našem hodnocení zaostává za ostatními automobily.