Toto je starší verze dokumentu!


U příkladů dávám vždy FIXME zkontrolovat aby ho zkontroloval alespoň ještě jeden další člověk. Prosím prvního člověka, kterému příklad vyjde stejně aby toto SMAZAL ať je vždy kontrola! Dííkes

Ze script

Kdo co máte ze sript, zkuste se poďelit aspoň vám to ostatní opraví.

Kapitola 1 Úvod

1.4 Cviceni (str. 9)

1.1 Najdete (uvahou) co nejjednoduzsi popis nasl. mnozin

Zadani a)

Graph

Reseni

Graph

Zadani b)

Graph

Reseni

Graph

Zadani c)

Graph

Reseni

Graph

Zadani d)

Graph

Reseni

Graph

Zadani e)

Graph

Reseni

Graph

Zadani f)

Graph

Reseni

Graph

Zadani g)

Graph

Reseni

Graph

Zadani h)

Graph

Reseni

Graph

1.2 Mejme nozinu bodu v rovine

Graph

Nacrtnete nasledujici mnoziny

Zadani a)

Graph

Reseni

Je to usecka od (-1,0) do (1,0)

Zadani b)

Graph

Reseni

Tohle je utvar ze dvou usecek (-1,2) az (1,2) a (-1,-2) az (1,-2) a dvou polokruznic, ktere je spojuji.

Kapitola 2 Vektory a Matice

Kapitola 3 Skalární součin a ortogonalita

Kapitola 4 Metoda nejmenších čtverců

Kapitola 5 Vlastní čísla a kvadratické formy

Kapitola 6 Množiny a zobrazení v euklidovských prostorech

Kapitola 7 Derivace

Kapitola 8 Analytické podmínky na lokální extrémy

Kapitola 9 Numerické algoritmy na hledání volných lokálních extrémů

Kapitola 10 Konvexní množiny

Kapitola 11 Lineární programování

Kapitola 12 Simplexová metoda

Kapitola 13 Dualita v lineárním programování

Ze zkoušek

Zkouška 23.1.2012

Minimální bodová hranice na Ečko byla snížena z 25 na 20/50. Úlohu 5 nikdo nedal, tak se do hodnocení nepočítala. Histogram úspěšnosti:

#A=2    #B=4    #C=11    #D=19    #E=32    #F=41

1) Řeka

Zadání

Stojíte na břehu řeky široké 1km. Chcete se dostat ke stanu, který je na protějším břehu 1km po proudu. Chůzí po břehu se pohybujete 3km/h, plavete 2km/h. Proud řeky zanedbáme. Je to absolutní volej. Za jaký nejkratší čas se dostanete ke stanu?

Řešení

Po kraji řeky ujdeme nějakou vzdálenost, pak vlezeme do řeky a budeme plavat šikmo, přímo ke stanu.

Úsek který jdeme pěšky podél řeky si označíme Graph zbytek vzdálenosti podél řeky bude tedy . Podle pythagorovy věty budeme plavat vzdálenost Graph.

FIXME Přidat obrázek

Sestavíme si rovnici pro čas:

Graph

Hledáme minimum této funkce, proto jí zderivujeme a prfní derivaci položíme rovnou nule:

Graph

Graph

Graph

Graph

Graph

Graph

Graph

Graph

Graph

Odlogaritmujeme:

Graph

Graph

Graph

Graph

Na kalkulačce:

Graph

2) Vážené nejmenší čtverce

Zadání

Hldedejte minimum x metodou vážených nejmenších čtverců:

Graph

  1. Napište funkci v maticové formě Graph
  2. Najděte řešení v maticové podobě Graph

Řešení

Matice Graph je ctvercova matice, ktera ma na diagonale prvky Graph, jinak vsude nuly.

Graph

FIXME zkontrolovat

Tady nevim. Ja to zderivoval, polozil rovno nule a vyjadril z toho vektor Graph.

DERIVACE JE BUĎ

Graph

FIXME NEBO zda se mi, ze derivujeme slozenou funkci, a ze tedy v predchozi verzi chybi „derivace vnitrni funkce“ (Ax-b), coz by melo byt Acko) tedy derivace by mela vypadat takto: Graph Covynato? Dal pocitam s prvni variantou derivace.

Matice Graph je diagonalni, proto:

Graph

Roztransponujeme a roznasobime:

Graph

Polozime rovno nule:

Graph

Graph

Graph

FIXME zkontrolovat reseni

Řešil jsem to pro druhou a podle mě správnější variantu derivace Graph

upravím protože Graph

Položím derivaci rovno nule:

Graph

Graph

Graph

Graph

Graph

Graph což je stejné jako vyšlo v předchozím případě..

Souhlas s druhym resenim, vyjde to tak, i kdyz si f(x) nejdrive roznasobim a az pak derivuji

3) Nejkratší spojnice mimoběžek

Zadání

Máme dvě mimoběžky v prostoru, každá dána dvěma body Graph. Hledejte příčku, tj. nejkratší spojnici mezi mimoběžkami metodou nejmenších čtverců.

Graph

Řešení

Kdybych si to byl byval namaloval, mohl to byt celkem v pohode priklad =(

Nasleduji myslenkove pochody a az nakonec vysledna forma.

Abychom mohli „pochodovat“ po prvni primce nadefinujeme si totok:

Graph

Kde Graph je smerovy vektor prvni primky a Graph je promenna, kterou tento vektor natahujeme. Analogicky pro druhou primku:

Graph

Tim jsme si nadefinovali po jednom bodu na kazde primce. Rozdil techto bodu je vektor, jehoz velikost chceme ze zadani minimalizovat. Body tedy od sebe odecteme:

Graph

A zbavime se nehezkych minusu:

Graph

Nyni se to pokusime dostat do tvaru Graph

Graph

Odtud uz by to melo byt videt.

Graph

Graph

Graph

Vypočítané koeficienty a se pak dosadí zpět do rovnic pro jednotlivé přímky (první dvě rovnice) abychom dostali body příčky, kterou hledáme.

FIXME zkontrolovat spravnost

4) GPS, nejmenší čtverce

Zadání

Máme souřadnice satelitů Graph. V reálném světě by Graph, my počítáme obecné Graph. Naše souřadnice je Graph. Známe vzdálenosti od jednotlivých satelitů Graph.

Minimalizujeme:

Graph

  1. Napište matlabskou funkci x = gps(A,y,x0) která řeší úlohu čistou Gauss-Newtonovou metodou, kde Graph a kde Graph je počáteční odhad.
  2. Je funkce f konvexní? Napište důkaz z definice.
  3. Kde není funkce f diferencovatelná?
  4. Bude nediferencovatelnost dělat problém v reálu?

Řešení

Navrh alespon na prvni cast - kod v matlabu (mozna nejaka indexace nemusi sedet - snad jsem to ale nepopletl)

5) Brambůrky, lupínky s pomocníkem

:!::!: Příklad nikdo nedal. nebyl započítán do bodování. :!::!:

Zadání

Děláme brambůrky a lupínky. Na kilo brambůrků spotřebujem 2kg brambor a 0.4 litry oleje. Na kilo lupínků spotřebujem 1.5kg brambor a 0.2 litrů oleje. Brambůrky prodáváme za 120,- za kilo, lupínky za 76,- za kilo.

Máme pomocníka. Pomocník si nechá zaplatit 10,- za vyrobený kilogram zboží. Nad 30 kilogramů si ale nechává platit 20,- za kilogram. Nad 40 kilogramů čehokoli vyrobeného za den dáme pomocníkovi 40,- za kilo. (Funkce kolik peněz pomocník dostane, závislá na množství vyrobeného zboží je taková dvakrát nahoru zlomená čára. )

Předpokládáme, že na začátku dne máme nakoupeno 100kg brambor a 10litrů oleje. Co se nespotřebuje to se vyhodí. (Cenu surovin jsem si nezapisoval, páč je nám k ničemu).

  1. vyrobte LP
  2. simplex tabulka
  3. jeden krok simplex metody

Řešení

Kdyby nebyla pomocníkova cena dvakrát zlomená, vypadala by LP snad takhle:

max  120 l  +  76 h  -  10 (l+h)
z.p.   2 l  + 1.5 h  <=  100
     0.4 l  + 0.2 h  <=   10

FIXME zkontrolovat

Myslim, ze se to da resit takto: (vcetne te „lomene“ fce, kterou je trik radsi nelamat)

6) Maxima, minima, nic-ima

Zadání

Máme funkci Graph na množině Graph. Určete lokální extrémy a zda jsou maxima, nebo mimnima.

Řešení

Parciální derivace položíme rovny nule:

Graph

Graph

Z první rovnice:

Graph

Dosadíme do druhé:

Graph

Odtud:

Graph

A tedy:

Graph

Graph

Toto je sice stacionární bod zadané funkce, ale ani nemusíme řešit, jestli je to extrém. Ono totiž Graph a tím nám to vypadne z množiny, pro kterou máme příklad řešit.

FIXME zkontrolovat

FIXME dodělat Teď mi ale došlo, že musíme teda najít, kterým směrem ta fce soupá a padá a hranice toho „definičního oboru“ označit za extrémy.

lagrange nerovnosti „hrubou silou“ zkusit to vyřešit pomocí Langrangeových multiplikátorů vyzkoušením všech možností (4)

nějaké rozumné extrémy mi vyšly jen při aktivaci omezení x-y=1, při x-y=-1 vyjdou komplexní kořeny nerovnice a při aktivaci obou nemá soustava řešení

mám tedy 2 řešení Graph

vypočítám si Hessian Graph dosadím a pro první bod mi vyjde indefinitní a stejně tak pro druhou, oba tedy budou nejspíše sedla..

7) Linprog, duál, podm.kompl.

Zadání

Máme zadaný vektor a číslo Graph a řešíme Graph

  1. řešte úvahou
  2. napište duál
  3. napište podmínky komplementarity

Řešení

Úvaha: Stačí si představit co dělají podmínky. Součet všech čísel v x je menší nebo rovno k, a my maximalizujeme, tedy bude rovno k. Navíc jednotlivá čísla v x jsou od nuly do jedné. Vezmeme tedy k jedniček a narveme je do x na stejné souřadnice, jako je ve vektoru c, k-největších čísel.

\\edit Tady je podle mě zrada, je třeba vzít v úvahu, že Graph může být třeba celé záporné a tedy výsledek bude Graph

Doporučuju se podívat na str. 94 do skript. Rovnice úplně dole vyjadřují téměř přesně tuhle úlohu.

Graph

Primár pak vypadá takto:

Graph

Graph

A od něj duál:

Graph

Graph

\\edit tohle je podle mě špatně a zneménka mají být opačně…

Graph

Podmínky komplementarity pak stačí opsat ze skript z 13.2

Zkouška 6.2.2012

1) Středoškolská úloha - žebřík

Zadání

Zeď od které se svažuje zem pod úhlem alfa (zadaným ve formě směrnice k = tg(alfa)). Kam máme postavit žebřík (vodorovná vzdálenost x) délky 1 (jedna) tak, aby dosáhl co nejvýš?

|
|\          /
| \        /
|  \      /
|   \    /
|    \  /
|     \/
|     /
|    /
|   / )
|  /   )
| / alfa)
+-----------
|--x--|

Řrešení

Rozdělíme si výšku do které dosáhne žebřík na dvě části tak, že vedeme rovnoběřku s osou x patou žebříku.

Tím nám vzniknou dva trojúhelníky.

Horní trojúhelník má šikmou stranu žebřík, spodní stranu o délce x a vlevo je svislá strana, tedy část maximalizované výšky. Z tohoto trojúhelníka vyjádříme výšku podle pythagorovy věty jako Graph

Spodní trojúhelník má šikmou stěnu na „svahu“ vodorovná je opět rovna x a svislá je druhá část výšky. druhou část dopočítáme podle tan(alfa) Graph

Graph

Toto se zderivuje polozi rovno nule a vyjadri se z toho x.

2) Extremy

Zadani

Najte extremy a reknete co to je

Graph

Reseni

Parcialni derivace (jakobian) polozime rovno nule … vyjdou body (0,0) a (1,1)

Druha derivace (hessian) vyjde jako ( 1 -1 ; -1 2y ) a po dosazeni vyjde pro (0,0) indefinitni a (1,1) pozitivne definitni, tedy (1,1) je lok. minimum a (0,0) sedlo.

3) Simplex

Zadani

Na lisu na plasty vyrabime bud hadicky nebo trubicky.

hadicky  = 25 kc/kg
trubicky = 30 kc/kg

vyrobime 200 kg/h hadicek   maximalne ale 6000kg
vyrobime 140 kg/h trubicek  maximalne ale 4000kg

mame maximalne 40 hodin provozu lisu

Vse ostatni taknejak ignorujem =))

Reseni

Pomerne hezke reseni (diky Petru Nobstovi) bylo prevest si vsechno na hodiny a koruny. LPcko pak vypada takhle:

max  5000*t1  +  4200*t2
z.p.      t1              + u1 = 30
                      t2  + u2 = 28,57    <<<--- 4000 / 140
           t1 +       t2  + u3 = 40

Tak jak to je se to namlati do simplexove tabulky, (bacha je to maximum, tak nejlip prevest na minimum).

Nakonec vyjde, ze vyrabime

hadicek   30hod  *  200kg/h  =  6000 kg
trubicek  10hod  *  140kg/h  =  1400 kg

4) lagaranze

Zadani

minimalizute a^Tx za podminky xTCx = 1 , C je symetricka pozitivne definitni matice.

Graph

b) jak se zmeni kdyz misto Graph bude Graph

c) graficka reprezentace pro n = 2

courses/a4b33opt/vypocitane_priklady.1391357075.txt.gz · Poslední úprava: 2025/01/03 18:24 (upraveno mimo DokuWiki)
Nahoru
chimeric.de = chi`s home Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0