Minimalizace kvadratické funkce s kvadratickými ... - FEI VŠB

Minimalizace kvadratické funkce s kvadratickými ... - FEI VŠB Minimalizace kvadratické funkce s kvadratickými ... - FEI VŠB

22.11.2014 Views

28 1 >> c = 1 c = 1 >> e = 0.001 e = 1.0000e−003 >> minimum(A,b,c,[−1;−2],e) ans = 0.7071 0.7071 −0.9142 Výpis 7: MNS s projekcí: Rěšení 1. příkladu Přidáním zvyšujícího se iterátoru do algoritmu lze získat počet iterací potřebných k výpočtu function [reseni] = minimum(A,b,c,bod0,e) ... it =0; while (norm(gradient)>e) it = it +1; ... end reseni=[bod(1) bod(2) vrat z (bod,A,b) it ]; Výpis 8: MNS s projekcí: Rěšení 1. příkladu s počítadlem iterací V našem případě bylo potřebných 24 iterací. Pokud bychom zkoušeli volit různé počáteční aproximace, zjistily bychom, že pro různé počáteční aproximace je rychlost konvergence (počet iterací potřebných pro provedení algoritmu) různá. Na obrázku níže jsou naznačeny některé body, u kterých je vyznačen potřebný počet iterací. Příklad 4.6 Minimalizujte kvadratickou funkci (12) vzhledem k množině (13), pokud A = V matlabu zavoláme funkci (6) [ 6 −1 −1 3 ] [ 1 , b = −2 >> minimum([6 −1; −1 3],[1;−2],2,[−2;−1],0.0001) ] , c = 2

29 Obrázek 13: MNS s projekcí: Rěšení 1. příkladu ans = −0.0917 −1.4112 0.1524 89.0000 Výpis 9: MNS s projekcí: Rěšení 2. příkladu Tedy minimum je v bodě x 0 = [−0.09; −1.41], f(x 0 ) = 0.15 a bylo potřebných 89 iterací. .

29<br />

Obrázek 13: MNS s projekcí: Rěšení 1. příkladu<br />

ans =<br />

−0.0917 −1.4112 0.1524 89.0000<br />

Výpis 9: MNS s projekcí: Rěšení 2. příkladu<br />

Tedy minimum je v bodě x 0 = [−0.09; −1.41], f(x 0 ) = 0.15 a bylo potřebných 89 iterací.<br />

.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!