Simplexový algoritmus

Šablona:Různé významy Simplexový algoritmus nebo také simplexová metoda je algoritmus pro řešení úlohy lineárního programování, který byl poprvé popsán George Dantzigem. Algoritmus efektivně prohledává tzv. základní řešení úloh lineárního programování, kterých je konečný počet a hledá mezi nimi řešení optimální. Optimální řešení je takové základní řešení, které poskytuje nejlepší hodnotu účelové funkce. Metoda souvisí s vlastnostmi polytopu v N dimenzionálním euklidovském prostoru. Řešená úloha je tak i graficky interpretovatelná – hledají se vlastně co nejvzdálenější vrcholy polytopu.

Popis řešených úloh zpracovávaných simplexovým algoritmem

Úloha lineárního programování je taková úloha, která hledá extrém zadané kriteriální neboli účelové funkce při zadaných omezujících podmínkách. V ekonomických úlohách jsou přidány podmínky nezápornosti proměnných modelů z důvodu interpretace proměnných.

Mějme úlohu lineárního programování ve standardním tvaru:

Maximalizujte účelovou funkci z$𝐳={𝐜}^T𝐱$

vůči omezením $𝐀𝐱\le 𝐛,𝐱\ge 0$

• A je matice strukturních koeficientů o velikosti m řádků x n sloupců (jinými slovy máme m omezení, každé max. n proměnných),
• vektor x je vektor proměnných modelu s n složkami,
• vektor b je vektor omezení nebo také vektor pravých stran a má m složek.
• $c^T$ je transponovaný vektor cenových koeficientů

Průběh algoritmu klasické simplexové metody

Algoritmus simplexové metody

Následující popis simplexového algoritmu neřeší speciální případy. Těmi jsou zejména:

• Neomezenost úlohy
• Degenerovanost úlohy
• Perturbace úlohy
• Nepřípustnost úlohy

1. Start – přepíšeme na ekvivalentní pomocnou úlohu

Zavedeme nové proměnné, tzv. přídatné proměnné, které vyrovnají omezující podmínky z nerovnic na rovnice. Přídatné proměnné definujeme následujícím způsobem:

${𝐱}_{n+i}={𝐛}_i-\sum _{j=1}^n{𝐚}_{ij}{𝐱}_j$

${𝐱}_{n+i}\ge 0,i=1...m$

Úlohu přepíšeme na tzv. ekvivalentní soustavu rovnic ve tvaru:

Maximalizujte ${𝐝}^{T}x$ vzhledem k omezením:

$𝐁𝐱\mathbf{=}𝐛$

$𝐱\ge 0$

kde:

$𝐁=(𝐀|𝐈),𝐈$ – je jednotková matice (identity)

$𝐝=(𝐜|\mathrm{𝟎})$

$𝐱=({𝐱}_{\mathrm{𝟏}}{𝐱}_{\mathrm{𝟐}}\mathbf{.}\mathbf{.}\mathbf{.}{𝐱}_{𝐧}\mathbf{|}{𝐱}_{𝐧\mathbf{+}\mathrm{𝟏}}\mathbf{.}\mathbf{.}\mathbf{.}{𝐱}_{𝐧\mathbf{+}𝐦})$

Při výpočtu je vhodné používat zápis v tzv. simplexové tabulce, která má následující výchozí tvar:

$A$	$I$	$b$
$-c^T$	$0$	$0$

kde:

A je matice strukturních koeficientů velikosti m x n

I je jednotková matice řádu m x m

b je vektor pravých stran o velikosti m x 1

c^T je transponovaný vektor cenových koeficientů

2. Získání výchozího přípustného řešení

Pokud položíme všechny přídatné proměnné v jednotlivých omezeních pravým stranám těchto omezení a všechny ostatní proměnné se budou rovnat nule, získáme výchozí přípustné řešení.

Nenulové proměnné nazveme tzv. bazickými proměnnými, nulové proměnné jsou tzv. nebazické. Dá se používat i termínů základní a nezákladní proměnné.

3. Test optima a výběr vstupující proměnné

Řešení je optimální právě tehdy, pokud v maximalizační úloze platí, že všechny redukované a stínové ceny jsou větší rovno nule nebo v minimalizační úloze platí, že všechny redukované a stínové ceny jsou menší rovno nule. Pokud je daná podmínka splněna, řešení je optimální a hodnota účelové funkce je nejvyšší možná.

Pokud je uvedená podmínka v daném typu úlohy porušena, řešení není optimální, tedy hodnotu účelové funkce lze zlepšit. Vybereme proměnnou, která nejvíce porušuje kritérium optimality, odstraníme ji z báze a nahradíme výhodnější proměnnou.

Nová vstupující proměnná se nachází v tzv. klíčovém řádku k, kde $x_k$ určíme následujícím způsobem:

v maximalizační úloze k je index řádku, kde $x_k=\min z_j,j=1..n$

v minimalizační úloze k je index řádku, kde $x_k=\max z_j,j=1..n$

4. Výběr vystupující proměnné

Vystupující proměnnou určíme z čísla t_i přiřazenému každému omezení s $a_{ik}>0$, kde k je index klíčového sloupce. Vybíráme

$\min t_i=\frac{b_i}{a_{ik}}$

index i je index klíčového řádku, který označíme q. Potom $a_{qk}$ je klíčový prvek.

5. Transformace tabulky

Tabulku upravíme pomocí Gaussovy eliminační metody, kdy pro klíčový řádek platí:

• $a_{qk}^{new}=1$
• $a_{qj}^{new}=\frac{a_{qj}}{a_{qk}}$
• $b_i^{new}=\frac{b_i}{a_{qk}}$

a ostatní řádky upravíme tak, aby v klíčovém řádku vznikl jednotkový vektor s jedničkou v klíčovém řádku. Gaussova eliminace se týká i vektoru pravé strany a řádku účelové funkce v simplexové tabulce.

Po transformaci pokračujeme na test optima.

Tabulka se po transformaci mění výpočtem na následující tvar:

$B_s^{-1}A$	$B_s^{-1}$	$B_s^{-1}b$
$c_B^T{B}_s^{-1}A-c^T$	$c_B^T{B}_s^{-1}$	$c_B^T{B}_s^{-1}b$

kde

$B_s^{-1}$ je inverzní matice báze

$c_B^T$ je transponovaný vektor cen bazických proměnných

Demonstrační úloha

Maximalizujte funkci:

$𝐳={𝐱}_1+{𝐱}_2$

vůči omezením

$-{𝐱}_1+{𝐱}_2\le 1$

${𝐱}_1\le 3$

${𝐱}_2\le 2$

${𝐱}_1,{𝐱}_2\ge 0$

Řešení

1. Zavedeme nové proměnné a zapíšeme do tabulky:

${𝐱}_3=1-(-1{𝐱}_1+1{𝐱}_2)$

${𝐱}_4=3-(1{𝐱}_1+0{𝐱}_2)$

${𝐱}_5=2-(0{𝐱}_1+1{𝐱}_2)$

$z^{\prime }={𝐱}_1+{𝐱}_2$

Po úpravě:

${𝐱}_3=1+1{𝐱}_1-1{𝐱}_2$

${𝐱}_4=3-1{𝐱}_1$

${𝐱}_5=2-1{𝐱}_2$

$z^{\prime }={𝐱}_1+{𝐱}_2$

2. Hledáme kandidáty do báze:

1. Vezměme si např. x₁ – mohli bychom i x₂ (obě proměnné mají v účelové funkci kladné a v tabulce záporné znaménko), ale je to pomalejší.

2. Jediná rovnice, která nám klade na x₁ omezení je druhá rovnice. Podle ní je x₁ maximálně 3. (Kdybychom zvolili x₂, nejpřísnější omezení by plynulo z první rovnice: x₂ je max. 1)

3. vyjádříme x₁ z druhé rovnice, dosadíme x₁ do zbývajících rovnic a do účelové funkce. Dopracujeme se k následující soustavě:

${𝐱}_1=3-{𝐱}_4$

${𝐱}_3=4-{𝐱}_4-{𝐱}_2$

${𝐱}_5=2-{𝐱}_2$

$z^{\prime }=3-{𝐱}_4+{𝐱}_2$

4. Nyní nám již zbývá jediná proměnná s kladným znaménkem v účelové funkci a současně se záporným znaménkem v rovnicích z tabulky – je to x₂

5. Postupujme obdobně jako před chvílí. Nejpřísnější omezení klade poslední rovnice – x₂ je podle ní max. dva. Vyjádřeme tedy x₂ z poslední rovnice a dosaďme do účelové funkce a ostatních rovnic:

${𝐱}_2=2-{𝐱}_5$

${𝐱}_3=2-{𝐱}_4+{𝐱}_5$

${𝐱}_1=3-{𝐱}_4$

$z^{\prime }=5-{𝐱}_4-{𝐱}_5$

6. Další analogická úprava není možná – účelová funkce nemůže být dále maximalizována. Víme (předpoklady úlohy), že x₄ ani x₅ není menší, než nula. Maximální hodnota účelové funkce z’ je menší, nejvýše rovna 5 pro vektor (x₁,…,x₅) = (3,2,2,0,0).

Literatura

• Lagová M., Jablonský J., Lineární modely, Praha, 2004, nakladatelství Oeconomica, Šablona:ISBN
• Jablonský J., Programy pro matematické modelování, Praha, 2007, nakladatelství Oeconomica, Šablona:ISBN
• Jablonský J., Operační výzkum – kvantitativní metody pro ekonomické rozhodování, Praha, 2002, Professional Publishing, Šablona:ISBN

Externí odkazy

• Šablona:Commonscat
• Aplet využívající simplexového algoritmu