• Nebyly nalezeny žádné výsledky

Matematika 5 FSV UK, ZS 2018-19 Miroslav Zelený

N/A
N/A
Protected

Academic year: 2022

Podíl "Matematika 5 FSV UK, ZS 2018-19 Miroslav Zelený"

Copied!
99
0
0

Načítání.... (zobrazit plný text nyní)

Fulltext

(1)

Matematika 5

FSV UK, ZS 2018-19

Miroslav Zelený

(2)

1. Stabilita ˇrešení soustav diferenciálních rovnic 2. Úvod do variaˇcního poˇctu

3. Globální extrémy

4. Teorie optimálního ˇrízení 5. R ˚uzné

(3)

1. Stabilita ˇrešení soustav diferenciálních rovnic

(4)

Budeme uvažovat rovnici

x0(t) = f(x(t)), (1.1) kdef je zobrazení definované na neprázdné otevˇrené množin ˇeG ⊂Rn s hodnotami vRn. Soustavy tvaru (1.1) se nazývajíautonomní.

V dalším budeme pˇredpokládat, žef ∈ C1(G,Rn), tj. složkyf1, . . . ,fn zobrazeníf jsou tˇrídyC1 naG.

Definice

Rekneme, žeˇ a∈G jestacionární bod rovnice(1.1), jestližef(a) = o.

(5)

Budeme uvažovat rovnici

x0(t) = f(x(t)), (1.1) kdef je zobrazení definované na neprázdné otevˇrené množin ˇeG ⊂Rn s hodnotami vRn. Soustavy tvaru (1.1) se nazývajíautonomní.

V dalším budeme pˇredpokládat, žef ∈ C1(G,Rn), tj.

složkyf1, . . . ,fn zobrazeníf jsou tˇrídyC1 naG.

Definice

Rekneme, žeˇ a∈G jestacionární bod rovnice(1.1), jestližef(a) = o.

(6)

Budeme uvažovat rovnici

x0(t) = f(x(t)), (1.1) kdef je zobrazení definované na neprázdné otevˇrené množin ˇeG ⊂Rn s hodnotami vRn. Soustavy tvaru (1.1) se nazývajíautonomní.

V dalším budeme pˇredpokládat, žef ∈ C1(G,Rn), tj.

složkyf1, . . . ,fn zobrazeníf jsou tˇrídyC1 naG.

Definice

Rekneme, žeˇ a∈G jestacionární bod rovnice(1.1), jestližef(a) =o.

(7)

Definice

Rekneme, že stacionární bodˇ a∈G rovnice (1.1) je

stabilní, jestliže pro každéε >0 existujeδ >0 takové, že pro každé maximální ˇrešeníx rovnice (1.1) spl ˇnujícíkx(0)−ak< δ platí:

(a) definiˇcní obor ˇrešeníx obsahuje intervalh0,+∞); (b) kx(t)−ak< εprot ∈ h0,+∞);

nestabilní, jestliže není stabilní,

asymptoticky stabilní, jestliže je stabilní a navíc existuje∆>0 takové, že pro každé maximální ˇrešení x rovnice (1.1) spl ˇnujícíkx(0)−ak<∆platí

limt→+∞x(t) = a.

(8)

Definice

Rekneme, že stacionární bodˇ a∈G rovnice (1.1) je stabilní, jestliže pro každéε >0 existujeδ >0 takové, že pro každé maximální ˇrešeníx rovnice (1.1) spl ˇnujícíkx(0)−ak< δ platí:

(a) definiˇcní obor ˇrešeníx obsahuje intervalh0,+∞); (b) kx(t)−ak< εprot ∈ h0,+∞);

nestabilní, jestliže není stabilní,

asymptoticky stabilní, jestliže je stabilní a navíc existuje∆>0 takové, že pro každé maximální ˇrešení x rovnice (1.1) spl ˇnujícíkx(0)−ak<∆platí

limt→+∞x(t) = a.

(9)

Definice

Rekneme, že stacionární bodˇ a∈G rovnice (1.1) je stabilní, jestliže pro každéε >0 existujeδ >0 takové, že pro každé maximální ˇrešeníx rovnice (1.1) spl ˇnujícíkx(0)−ak< δ platí:

(a) definiˇcní obor ˇrešeníx obsahuje intervalh0,+∞);

(b) kx(t)−ak< εprot ∈ h0,+∞);

nestabilní, jestliže není stabilní,

asymptoticky stabilní, jestliže je stabilní a navíc existuje∆>0 takové, že pro každé maximální ˇrešení x rovnice (1.1) spl ˇnujícíkx(0)−ak<∆platí

limt→+∞x(t) = a.

(10)

Definice

Rekneme, že stacionární bodˇ a∈G rovnice (1.1) je stabilní, jestliže pro každéε >0 existujeδ >0 takové, že pro každé maximální ˇrešeníx rovnice (1.1) spl ˇnujícíkx(0)−ak< δ platí:

(a) definiˇcní obor ˇrešeníx obsahuje intervalh0,+∞);

(b) kx(t)−ak< εprot ∈ h0,+∞);

nestabilní, jestliže není stabilní,

asymptoticky stabilní, jestliže je stabilní a navíc existuje∆>0 takové, že pro každé maximální ˇrešení x rovnice (1.1) spl ˇnujícíkx(0)−ak<∆platí

limt→+∞x(t) = a.

(11)

Definice

Rekneme, že stacionární bodˇ a∈G rovnice (1.1) je stabilní, jestliže pro každéε >0 existujeδ >0 takové, že pro každé maximální ˇrešeníx rovnice (1.1) spl ˇnujícíkx(0)−ak< δ platí:

(a) definiˇcní obor ˇrešeníx obsahuje intervalh0,+∞);

(b) kx(t)−ak< εprot ∈ h0,+∞);

nestabilní, jestliže není stabilní,

asymptoticky stabilní, jestliže je stabilní a navíc existuje∆>0 takové, že pro každé maximální ˇrešení x rovnice (1.1) spl ˇnujícíkx(0)−ak<∆platí

limt→+∞x(t) = a.

(12)

V ˇeta 1.1

Necht’A∈M(n×n).

Stacionární bod o rovnicex0 =Ax je asymptoticky stabilní, práv ˇe když<λ <0pro každé vlastní ˇcísloλ maticeA.

Stacionární bod o rovnicex0 =Ax je stabilní, práv ˇe když<λ ≤0pro každé vlastní ˇcísloλmaticeAa pokud<λ =0, pak násobnostλje rovna

n−h(λI−A).

(13)

V ˇeta 1.1

Necht’A∈M(n×n).

Stacionární bod o rovnicex0 =Ax je asymptoticky stabilní, práv ˇe když<λ <0pro každé vlastní ˇcísloλ maticeA.

Stacionární bod o rovnicex0 =Ax je stabilní, práv ˇe když<λ ≤0pro každé vlastní ˇcísloλmaticeAa pokud<λ =0, pak násobnostλje rovna

n−h(λI−A).

(14)

V ˇeta 1.2 (Ljapunov)

Necht’f ∈ C1(G,Rn)aaje stacionární bod rovnice (1.1).

Oznaˇcme

A= ∂fi

∂xj

(a)

i=1..n,j=1..n

. Pak platí:

Jestliže každé vlastní ˇcíslo maticeAmá zápornou reálnou ˇcást, paka je asymptoticky stabilní bod rovnice (1.1).

Jestliže alespo ˇn jedno vlastní ˇcíslo maticeAmá kladnou reálnou ˇcást, pak jeanestabilní bod rovnice (1.1).

(15)

V ˇeta 1.2 (Ljapunov)

Necht’f ∈ C1(G,Rn)aaje stacionární bod rovnice (1.1).

Oznaˇcme

A= ∂fi

∂xj

(a)

i=1..n,j=1..n

. Pak platí:

Jestliže každé vlastní ˇcíslo maticeAmá zápornou reálnou ˇcást, paka je asymptoticky stabilní bod rovnice (1.1).

Jestliže alespo ˇn jedno vlastní ˇcíslo maticeAmá kladnou reálnou ˇcást, pak jeanestabilní bod rovnice (1.1).

(16)

λ1>0, λ2>0, λ16=λ2

(17)

λ1<0, λ2<0, λ16=λ2

(18)

λ1<0< λ2

(19)

λ1C\R,1>0

(20)

λ1C\R,1<0

(21)

0 0.5 1 1.5 2 2.5 3

y

-1 -0.5 0.5 1 1.5 2

x

(22)

0 0.5 1 1.5 2 2.5 3

y

-1 -0.5 0.5 1 1.5 2

x

(23)

0 0.5 1 1.5 2 2.5 3

y

-1 -0.5 0.5 1 1.5 2

x

(24)

0 0.5 1 1.5 2 2.5 3

y

-1 -0.5 0.5 1 1.5 2

x

(25)

2. Úvod do varia ˇcního po ˇctu

(26)

2.1 Derivování funkcí na vektorových prostorech

Definice

Necht’X je vektorový prostorF :X →R,a∈X,h∈X. Derivací funkceF v bod ˇeave sm ˇeruh rozumíme

δF(a,h) = lim

t→0

F(a+th)−F(a)

t ,

pokud limita existuje vlastní.

(27)

Definice

Necht’X je reálný vektorový prostor,M ⊂X,a∈M aF je reálná funkce definovaná alespo ˇn naM. ˇRekneme, že a jebod minima(resp. bod maxima) funkceF na množin ˇe M, jestliže pro každéx ∈M platíF(x)≥F(a)(resp.

F(x)≤F(a)).

V ˇeta 2.1

Necht’ X je vektorový prostor, F :X →Ra a∈X . Jestliže má F v bod ˇe a extrém (tj. minimum nebo maximum), pak pro každé h∈X platí, žeδF(a,h)neexistuje nebo je rovna nule.

(28)

Definice

Necht’X je reálný vektorový prostor,M ⊂X,a∈M aF je reálná funkce definovaná alespo ˇn naM. ˇRekneme, že a jebod minima(resp. bod maxima) funkceF na množin ˇe M, jestliže pro každéx ∈M platíF(x)≥F(a)(resp.

F(x)≤F(a)).

V ˇeta 2.1

Necht’ X je vektorový prostor, F :X →Ra a∈X . Jestliže má F v bod ˇe a extrém (tj. minimum nebo maximum), pak pro každé h∈X platí, žeδF(a,h)neexistuje nebo je rovna nule.

(29)

2.2 Derivování integrálu

Definice

Necht’M ⊂Rn. ˇRekneme, že funkcef jestejnom ˇern ˇe spojitá naM, jestliže platí

∀ε >0∃δ >0∀x,y ∈M,||x −y||< δ : |f(x)−f(y)|< ε.

V ˇeta 2.2

Necht’ K ⊂Rn je kompaktní a f :K →Rje spojitá na K . Potom f je stejnom ˇern ˇe spojitá na K .

(30)

2.2 Derivování integrálu

Definice

Necht’M ⊂Rn. ˇRekneme, že funkcef jestejnom ˇern ˇe spojitá naM, jestliže platí

∀ε >0∃δ >0∀x,y ∈M,||x −y||< δ : |f(x)−f(y)|< ε.

V ˇeta 2.2

Necht’ K ⊂Rn je kompaktní a f :K →Rje spojitá na K . Potom f je stejnom ˇern ˇe spojitá na K .

(31)

V ˇeta 2.3

Necht’ f: (a,b)×(c,d)→Rje spojitá a∂1f (= parciálni derivace podle první prom ˇenné) je také spojitá na (a,b)×(c,d). Necht’ϕ: (a,b)→(c,d)je funkce, která má v každém bod ˇe vlastní derivaci. Necht’ x0∈(c,d).

Položme

K(y) = Z ϕ(y)

x0

f(y,x)dx pro y ∈(a,b). Potom má funkce K v každém bod ˇe intervalu(a,b) vlastní derivaci a platí

K0(y) = Z ϕ(y)

x0

1f(y,x)dx+f(y, ϕ(y))ϕ0(y), y ∈(a,b).

(32)

V ˇeta 2.3

Necht’ f: (a,b)×(c,d)→Rje spojitá a∂1f (= parciálni derivace podle první prom ˇenné) je také spojitá na (a,b)×(c,d). Necht’ϕ: (a,b)→(c,d)je funkce, která má v každém bod ˇe vlastní derivaci. Necht’ x0∈(c,d).

Položme

K(y) = Z ϕ(y)

x0

f(y,x)dx pro y ∈(a,b).

Potom má funkce K v každém bod ˇe intervalu(a,b) vlastní derivaci a platí

K0(y) = Z ϕ(y)

x0

1f(y,x)dx+f(y, ϕ(y))ϕ0(y), y ∈(a,b).

(33)

V ˇeta 2.3

Necht’ f: (a,b)×(c,d)→Rje spojitá a∂1f (= parciálni derivace podle první prom ˇenné) je také spojitá na (a,b)×(c,d). Necht’ϕ: (a,b)→(c,d)je funkce, která má v každém bod ˇe vlastní derivaci. Necht’ x0∈(c,d).

Položme

K(y) = Z ϕ(y)

x0

f(y,x)dx pro y ∈(a,b).

Potom má funkce K v každém bod ˇe intervalu(a,b) vlastní derivaci a platí

K0(y) = Z ϕ(y)

x0

1f(y,x)dx+f(y, ϕ(y))ϕ0(y), y ∈(a,b).

(34)

2.3 Základní úloha variaˇcního poˇctu

Formulace základní úlohy varia ˇcního po ˇctu (P1).

Dáno:T ∈R,T >0;A,Z ∈R;F ∈ C2(h0,Ti ×R×R)

Hledáme takovéy ∈ C1(h0,Ti),y(0) = A,y(T) =Z, že hodnota

V(y) = Z T

0

F(s,y(s),y0(s))ds je minimální (resp. maximální).

(35)

2.3 Základní úloha variaˇcního poˇctu

Formulace základní úlohy varia ˇcního po ˇctu (P1).

Dáno:T ∈R,T >0;A,Z ∈R;F ∈ C2(h0,Ti ×R×R) Hledáme takovéy ∈ C1(h0,Ti),y(0) = A,y(T) =Z, že hodnota

V(y) = Z T

0

F(s,y(s),y0(s))ds je minimální (resp. maximální).

(36)

V ˇeta 2.4 (nutná podmínka pro extrém)

Necht’ y je bodem extrému pro (P1). Pak y je ˇrešením rovnice

(ER1) Fy(t,y,y0) = d

dtFy0(t,y,y0).

(37)

Lemma 2.5

Necht’ funkceϕ∈ C(h0,Ti)je nezáporná a RT

0 ϕ(s)ds=0. Pakϕ=0nah0,Ti.

Lemma 2.6 (základní lemma variaˇcního poˇctu)

Necht’ a,b∈ C(h0,Ti)a Z T

0

(a(t)h(t) +b(t)h0(t))dt =0 pro každou funkci h∈ C1(h0,Ti)spl ˇnující

h(0) =h(T) =0. Pak funkce b má nah0,Tiderivaci a platí zde b0 =a.

(38)

Lemma 2.5

Necht’ funkceϕ∈ C(h0,Ti)je nezáporná a RT

0 ϕ(s)ds=0. Pakϕ=0nah0,Ti.

Lemma 2.6 (základní lemma variaˇcního poˇctu)

Necht’ a,b∈ C(h0,Ti)a Z T

0

(a(t)h(t) +b(t)h0(t))dt =0 pro každou funkci h∈ C1(h0,Ti)spl ˇnující

h(0) =h(T) =0. Pak funkce b má nah0,Tiderivaci a platí zde b0 =a.

(39)

Lemma 2.7

Necht’ T a F jsou jako v (P1), y,u ∈ C1(h0,Ti). Necht’

zobrazení G:C1(h0,Ti)→Rje definováno takto:

G(u) = Z T

0

F(t,y(t) +u(t),y0(t) +u0(t))dt.

Potom

δG(o,h) =

Z T

0

Fy(t,y(t),y0(t))·h(t)+Fy0(t,y(t),y0(t))·h0(t) dt pro libovolné h∈ C1(h0,Ti).

(40)

Lemma 2.7

Necht’ T a F jsou jako v (P1), y,u ∈ C1(h0,Ti). Necht’

zobrazení G:C1(h0,Ti)→Rje definováno takto:

G(u) = Z T

0

F(t,y(t) +u(t),y0(t) +u0(t))dt. Potom

δG(o,h) =

Z T

0

Fy(t,y(t),y0(t))·h(t)+Fy0(t,y(t),y0(t))·h0(t) dt pro libovolné h∈ C1(h0,Ti).

(41)

2.4 Pevný koncový ˇcas a volná koncová hodnota

Formulace úlohy (P2).

Dáno:T ∈R,T >0;A∈R;F ∈ C2(h0,Ti ×R×R)

Hledáme takovéy ∈ C1(h0,Ti),y(0) = A, že hodnota V(y) =

Z T

0

F(s,y(s),y0(s))ds je minimální (resp. maximální).

(42)

2.4 Pevný koncový ˇcas a volná koncová hodnota

Formulace úlohy (P2).

Dáno:T ∈R,T >0;A∈R;F ∈ C2(h0,Ti ×R×R) Hledáme takovéy ∈ C1(h0,Ti),y(0) = A, že hodnota

V(y) = Z T

0

F(s,y(s),y0(s))ds je minimální (resp. maximální).

(43)

V ˇeta 2.8

Necht’ y je bodem extrému pro (P2). Pak y spl ˇnuje Fy(t,y,y0) = d

dtFy0(t,y,y0), (ER1) Fy0(T,y(T),y0(T)) =0. (ER2)

(44)

2.5 Isoperimetrická úloha

Formulace úlohy (P3).

Dáno:T ∈R,T >0;A,Z ∈R;B∈R;

F ∈ C2(h0,Ti ×R×R)

Hledáme takovéy ∈ C1(h0,Ti)spl ˇnujícíy(0) = A, y(T) =Z,

Z T

0

G(t,y,y0) = B, že hodnota

V(y) = Z T

0

F(s,y(s),y0(s))ds je minimální (resp. maximální).

(45)

2.5 Isoperimetrická úloha

Formulace úlohy (P3).

Dáno:T ∈R,T >0;A,Z ∈R;B∈R;

F ∈ C2(h0,Ti ×R×R)

Hledáme takovéy ∈ C1(h0,Ti)spl ˇnujícíy(0) = A, y(T) =Z,

Z T

0

G(t,y,y0) =B, že hodnota

V(y) = Z T

0

F(s,y(s),y0(s))ds je minimální (resp. maximální).

(46)

V ˇeta 2.9

Necht’ y je bodem extrému pro (P3). Pak y spl ˇnuje bud’

Gy(t,y,y0) = d

dtGy0(t,y,y0), (I) nebo existujeλ∈R, že y spl ˇnuje

Fy(T,y(T),y0(T))−λGy(T,y(T),y0(T)) = d

dt Fy0(t,y,y0)−λGy0(t,y,y0)

. (II)

(47)

3. Posta ˇcující podmínky pro extrém

(48)

3.1 Globální extrémy

Definice

Necht’X je vektorový prostor aV :X →Rje funkcionál.

Rekneme, žeˇ V jekonkávní (resp. konvexní) naX, jestliže

∀x,y ∈X∀t ∈ h0,1i: V(tx+(1−t)y)≥tV(x)+(1−t)V(y)

(resp.

∀x,y ∈X∀t ∈ h0,1i: V(tx+(1−t)y)≤tV(x)+(1−t)V(y)).

(49)

3.1 Globální extrémy

Definice

Necht’X je vektorový prostor aV :X →Rje funkcionál.

Rekneme, žeˇ V jekonkávní (resp. konvexní) naX, jestliže

∀x,y ∈X∀t ∈ h0,1i: V(tx+(1−t)y)≥tV(x)+(1−t)V(y) (resp.

∀x,y ∈X∀t ∈ h0,1i: V(tx+(1−t)y)≤tV(x)+(1−t)V(y)).

(50)

V ˇeta 3.1

Necht’ V :X →Rje konkávní. JestližeδV(x,h) =0pro každé h∈X , pak V má v x maximum.

V ˇeta 3.2

Necht’ F ∈ C2(h0,Ti ×R×R).

(K) Necht’ pro každé t ∈ h0,Tije funkce [y,y0]7→F(t,y,y0)konkávní.

Pak je funkcionál V: C1(h0,Ti))→Rdefinovaný pˇredpisem

V :y 7→ Z T

0

F(t,y,y0)dt konkávní.

(51)

V ˇeta 3.1

Necht’ V :X →Rje konkávní. JestližeδV(x,h) =0pro každé h∈X , pak V má v x maximum.

V ˇeta 3.2

Necht’ F ∈ C2(h0,Ti ×R×R).

(K) Necht’ pro každé t ∈ h0,Tije funkce [y,y0]7→F(t,y,y0)konkávní.

Pak je funkcionál V: C1(h0,Ti))→Rdefinovaný pˇredpisem

V :y 7→ Z T

0

F(t,y,y0)dt konkávní.

(52)

V ˇeta 3.1

Necht’ V :X →Rje konkávní. JestližeδV(x,h) =0pro každé h∈X , pak V má v x maximum.

V ˇeta 3.2

Necht’ F ∈ C2(h0,Ti ×R×R).

(K) Necht’ pro každé t ∈ h0,Tije funkce [y,y0]7→F(t,y,y0)konkávní.

Pak je funkcionál V: C1(h0,Ti))→Rdefinovaný pˇredpisem

V :y 7→

Z T

0

F(t,y,y0)dt konkávní.

(53)

V ˇeta 3.3

Necht’ F v (P1) spl ˇnuje (K). Pak je (ER1) postaˇcující podmínkou pro maximum.

(54)

3.2 Postaˇcující podmínky pro lokální extrém

Definice

Normovaným lineárním prostoremrozumíme dvojici (X,||.||), kdeX je vektorový prostor (nadR) a||.||je norma naX, tj. zobrazení||.||:X → h0,+∞)spl ˇnující

∀x ∈X :||x||=0⇔x =o,

∀x ∈X∀λ ∈R: ||λx||=|λ| · ||x||,

∀x,y ∈X : ||x+y|| ≤ ||x||+||y||.

(55)

3.2 Postaˇcující podmínky pro lokální extrém

Definice

Normovaným lineárním prostoremrozumíme dvojici (X,||.||), kdeX je vektorový prostor (nadR) a||.||je norma naX, tj. zobrazení||.||:X → h0,+∞)spl ˇnující

∀x ∈X :||x||=0⇔x =o,

∀x ∈X∀λ ∈R: ||λx||=|λ| · ||x||,

∀x,y ∈X : ||x+y|| ≤ ||x||+||y||.

(56)

3.2 Postaˇcující podmínky pro lokální extrém

Definice

Normovaným lineárním prostoremrozumíme dvojici (X,||.||), kdeX je vektorový prostor (nadR) a||.||je norma naX, tj. zobrazení||.||:X → h0,+∞)spl ˇnující

∀x ∈X :||x||=0⇔x =o,

∀x ∈X∀λ ∈R: ||λx||=|λ| · ||x||,

∀x,y ∈X : ||x+y|| ≤ ||x||+||y||.

(57)

3.2 Postaˇcující podmínky pro lokální extrém

Definice

Normovaným lineárním prostoremrozumíme dvojici (X,||.||), kdeX je vektorový prostor (nadR) a||.||je norma naX, tj. zobrazení||.||:X → h0,+∞)spl ˇnující

∀x ∈X :||x||=0⇔x =o,

∀x ∈X∀λ ∈R: ||λx||=|λ| · ||x||,

∀x,y ∈X : ||x+y|| ≤ ||x||+||y||.

(58)

Definice

Necht’(X,||.||)je normovaný lineární prostor, f :X →Ra x0 ∈X. ˇRekneme, žef má v bod ˇex0

lokální maximum, jestliže existujer >0 takové, že

∀x ∈X,||x−x0||<r : f(x)≤f(x0);

ostré lokální maximum, jestliže existujer >0 takové, že

∀x ∈X,0<||x−x0||<r : f(x)<f(x0). Analogicky definujemelokální minimum aostré lokální minimum.

(59)

Definice

Necht’(X,||.||)je normovaný lineární prostor, f :X →Ra x0 ∈X. ˇRekneme, žef má v bod ˇex0

lokální maximum, jestliže existujer >0 takové, že

∀x ∈X,||x−x0||<r : f(x)≤f(x0);

ostré lokální maximum, jestliže existujer >0 takové, že

∀x ∈X,0<||x−x0||<r : f(x)<f(x0).

Analogicky definujemelokální minimum aostré lokální minimum.

(60)

Definice

Necht’(X,||.||)je normovaný lineární prostor, f :X →Ra x0 ∈X. ˇRekneme, žef má v bod ˇex0

lokální maximum, jestliže existujer >0 takové, že

∀x ∈X,||x−x0||<r : f(x)≤f(x0);

ostré lokální maximum, jestliže existujer >0 takové, že

∀x ∈X,0<||x−x0||<r : f(x)<f(x0).

Analogicky definujemelokální minimumaostré lokální minimum.

(61)

V ˇeta 3.4

Necht’ y ˇreší (ER1) v úloze (P1). Jestliže je matice Fyy(t,y(t),y0(t)) Fyy0(t,y(t),y0(t))

Fyy0(t,y(t),y0(t)) Fy0y0(t,y(t),y0(t))

negativn ˇe definitní pro každé t ∈ h0,Ti, pak y je bodem ostrého lokálního maxima.

(62)

4. Teorie optimálního ˇrízení

(63)

Definice

Rekneme, že funkceˇ f jepo ˇcástech spojitána intervalu h0,Ti, jestliže existuje d ˇelení 0=t0 <t1 <· · ·<tn=T takové, žef|(ti,ti+1)je spojitá na(ti,ti+1)pro každé

i ∈ {0, . . . ,n−1}a v krajních bodech existují vlastní limity.

Rekneme, že funkceˇ f jepo ˇcástech diferencovatelná na intervaluh0,Ti, jestliže existuje d ˇelení

0=t0 <t1 <· · ·<tn=T takové, žef|(ti,ti+1)má na(ti,ti+1) vlastní derivaci a v krajních bodech existují pˇríslušné jednostranné vlastní derivace pro každéi ∈ {0, . . . ,n−1}.

(64)

Definice

Rekneme, že funkceˇ f jepo ˇcástech spojitána intervalu h0,Ti, jestliže existuje d ˇelení 0=t0 <t1 <· · ·<tn=T takové, žef|(ti,ti+1)je spojitá na(ti,ti+1)pro každé

i ∈ {0, . . . ,n−1}a v krajních bodech existují vlastní limity.

Rekneme, že funkceˇ f jepo ˇcástech diferencovatelná na intervaluh0,Ti, jestliže existuje d ˇelení

0=t0<t1<· · ·<tn=T takové, žef|(ti,ti+1)má na(ti,ti+1) vlastní derivaci a v krajních bodech existují pˇríslušné jednostranné vlastní derivace pro každéi ∈ {0, . . . ,n−1}.

(65)

Formulace úlohy (P4)

Dáno:

T ∈R,T >0;A∈R;

F ∈ C(h0,Ti ×R×R),1F,∂2F jsou spojité; f ∈ C(h0,Ti ×R×R),1f,∂2f jsou spojité; U je omezený uzavˇrený interval.

(66)

Formulace úlohy (P4)

Dáno:

T ∈R,T >0;A∈R;

F ∈ C(h0,Ti ×R×R),1F,∂2F jsou spojité;

f ∈ C(h0,Ti ×R×R),1f,∂2f jsou spojité; U je omezený uzavˇrený interval.

(67)

Formulace úlohy (P4)

Dáno:

T ∈R,T >0;A∈R;

F ∈ C(h0,Ti ×R×R),1F,∂2F jsou spojité;

f ∈ C(h0,Ti ×R×R),1f,∂2f jsou spojité;

U je omezený uzavˇrený interval.

(68)

Formulace úlohy (P4)

Dáno:

T ∈R,T >0;A∈R;

F ∈ C(h0,Ti ×R×R),1F,∂2F jsou spojité;

f ∈ C(h0,Ti ×R×R),1f,∂2f jsou spojité;

U je omezený uzavˇrený interval.

(69)

Formulace úlohy (P4)

Dáno:

T ∈R,T >0;A∈R;

F ∈ C(h0,Ti ×R×R),1F,∂2F jsou spojité;

f ∈ C(h0,Ti ×R×R),1f,∂2f jsou spojité;

U je omezený uzavˇrený interval.

(70)

Hledámey po ˇcástech diferencovatelnou na intervalu h0,Tiaupo ˇcástech spojitou na h0,Titakové, že

y(0) = A,

y0(t) = f(t,y(t),u(t))pro každét ∈ h0,Tivyjma koneˇcné množiny,

u(t)∈ U pro každét ∈ h0,Ti, Z T

0

F(t,y(t),u(t))dt je maximální.

(71)

Hledámey po ˇcástech diferencovatelnou na intervalu h0,Tiaupo ˇcástech spojitou na h0,Titakové, že

y(0) = A,

y0(t) = f(t,y(t),u(t))pro každét ∈ h0,Tivyjma koneˇcné množiny,

u(t)∈ U pro každét ∈ h0,Ti, Z T

0

F(t,y(t),u(t))dt je maximální.

(72)

Hledámey po ˇcástech diferencovatelnou na intervalu h0,Tiaupo ˇcástech spojitou na h0,Titakové, že

y(0) = A,

y0(t) =f(t,y(t),u(t))pro každét ∈ h0,Tivyjma koneˇcné množiny,

u(t)∈ U pro každét ∈ h0,Ti, Z T

0

F(t,y(t),u(t))dt je maximální.

(73)

Hledámey po ˇcástech diferencovatelnou na intervalu h0,Tiaupo ˇcástech spojitou na h0,Titakové, že

y(0) = A,

y0(t) =f(t,y(t),u(t))pro každét ∈ h0,Tivyjma koneˇcné množiny,

u(t)∈ U pro každét ∈ h0,Ti,

Z T

0

F(t,y(t),u(t))dt je maximální.

(74)

Hledámey po ˇcástech diferencovatelnou na intervalu h0,Tiaupo ˇcástech spojitou na h0,Titakové, že

y(0) = A,

y0(t) =f(t,y(t),u(t))pro každét ∈ h0,Tivyjma koneˇcné množiny,

u(t)∈ U pro každét ∈ h0,Ti, Z T

0

F(t,y(t),u(t))dt je maximální.

(75)

V ˇeta 4.1 (Pontrjagin ˚uv princip maxima)

Necht’ u je bod maxima v úloze (P4). Pak existuje funkce t 7→λ(t), že pro H(t,y,u, λ) = F(t,y,u) +λf(t,y,u)(tzv.

hamiltonián) platí:

(PM1) pro každé t ∈ h0,Tivyjma koneˇcné množiny a pro každéu˜∈ U platí

H(t,y(t),u(t), λ(t))≥H(t,y(t),u, λ(t)),˜

(PM2) y0 = ∂H

∂λ (stavová rovnice), (PM3) λ0 =−∂H

∂y (pohybová rovnice), (PM4) λ(T) = 0(podmínka transverzality).

(76)

V ˇeta 4.1 (Pontrjagin ˚uv princip maxima)

Necht’ u je bod maxima v úloze (P4). Pak existuje funkce t 7→λ(t), že pro H(t,y,u, λ) = F(t,y,u) +λf(t,y,u)(tzv.

hamiltonián) platí:

(PM1) pro každé t ∈ h0,Tivyjma koneˇcné množiny a pro každéu˜ ∈ U platí

H(t,y(t),u(t), λ(t))≥H(t,y(t),u, λ(t)),˜

(PM2) y0 = ∂H

∂λ (stavová rovnice), (PM3) λ0 =−∂H

∂y (pohybová rovnice), (PM4) λ(T) = 0(podmínka transverzality).

(77)

V ˇeta 4.1 (Pontrjagin ˚uv princip maxima)

Necht’ u je bod maxima v úloze (P4). Pak existuje funkce t 7→λ(t), že pro H(t,y,u, λ) = F(t,y,u) +λf(t,y,u)(tzv.

hamiltonián) platí:

(PM1) pro každé t ∈ h0,Tivyjma koneˇcné množiny a pro každéu˜ ∈ U platí

H(t,y(t),u(t), λ(t))≥H(t,y(t),u, λ(t)),˜

(PM2) y0 = ∂H

∂λ (stavová rovnice),

(PM3) λ0 =−∂H

∂y (pohybová rovnice), (PM4) λ(T) = 0(podmínka transverzality).

(78)

V ˇeta 4.1 (Pontrjagin ˚uv princip maxima)

Necht’ u je bod maxima v úloze (P4). Pak existuje funkce t 7→λ(t), že pro H(t,y,u, λ) = F(t,y,u) +λf(t,y,u)(tzv.

hamiltonián) platí:

(PM1) pro každé t ∈ h0,Tivyjma koneˇcné množiny a pro každéu˜ ∈ U platí

H(t,y(t),u(t), λ(t))≥H(t,y(t),u, λ(t)),˜

(PM2) y0 = ∂H

∂λ (stavová rovnice), (PM3) λ0 =−∂H

∂y (pohybová rovnice),

(PM4) λ(T) = 0(podmínka transverzality).

(79)

V ˇeta 4.1 (Pontrjagin ˚uv princip maxima)

Necht’ u je bod maxima v úloze (P4). Pak existuje funkce t 7→λ(t), že pro H(t,y,u, λ) = F(t,y,u) +λf(t,y,u)(tzv.

hamiltonián) platí:

(PM1) pro každé t ∈ h0,Tivyjma koneˇcné množiny a pro každéu˜ ∈ U platí

H(t,y(t),u(t), λ(t))≥H(t,y(t),u, λ(t)),˜

(PM2) y0 = ∂H

∂λ (stavová rovnice), (PM3) λ0 =−∂H

∂y (pohybová rovnice), (PM4) λ(T) = 0(podmínka transverzality).

(80)

4.2 Postaˇcující podmínky

V ˇeta 4.2

Princip maxima je postaˇcující podmínkou pro extrém v úloze (P4), jestliže

F a f jsou diferencovatelné,

F a f jsou konkávní v(y,u),

bud’ f je lineární v y a v u neboλ(t)≥0pro každé t ∈ h0,Ti.

(81)

4.2 Postaˇcující podmínky

V ˇeta 4.2

Princip maxima je postaˇcující podmínkou pro extrém v úloze (P4), jestliže

F a f jsou diferencovatelné, F a f jsou konkávní v(y,u),

bud’ f je lineární v y a v u neboλ(t)≥0pro každé t ∈ h0,Ti.

(82)

4.2 Postaˇcující podmínky

V ˇeta 4.2

Princip maxima je postaˇcující podmínkou pro extrém v úloze (P4), jestliže

F a f jsou diferencovatelné, F a f jsou konkávní v(y,u),

bud’ f je lineární v y a v u neboλ(t)≥0pro každé t ∈ h0,Ti.

(83)

4.3 Problémy s více stavovými prom ˇennými

Formulace úlohy (P4’) Dáno:

T ∈R,T >0;y0Rn;

F ∈ C1(h0,Ti ×Rn×Rm); f ∈ C1(h0,Ti ×Rn×Rm,Rn); U1, . . . ,UmR.

(84)

4.3 Problémy s více stavovými prom ˇennými

Formulace úlohy (P4’) Dáno:

T ∈R,T >0;y0Rn; F ∈ C1(h0,Ti ×Rn×Rm);

f ∈ C1(h0,Ti ×Rn×Rm,Rn); U1, . . . ,UmR.

(85)

4.3 Problémy s více stavovými prom ˇennými

Formulace úlohy (P4’) Dáno:

T ∈R,T >0;y0Rn; F ∈ C1(h0,Ti ×Rn×Rm);

f ∈ C1(h0,Ti ×Rn×Rm,Rn);

U1, . . . ,UmR.

(86)

4.3 Problémy s více stavovými prom ˇennými

Formulace úlohy (P4’) Dáno:

T ∈R,T >0;y0Rn; F ∈ C1(h0,Ti ×Rn×Rm);

f ∈ C1(h0,Ti ×Rn×Rm,Rn);

U1, . . . ,UmR.

(87)

Hledámey = [y1, . . . ,yn], kde složky jsou po ˇcástech diferencovatelné funkce na intervaluh0,Ti,

a funkci u= [u1, . . . ,um]s po ˇcástech spojitými složkami nah0,Ti takové, že

y(0) = y0,

y0(t) =f(t,y(t),u(t))pro každét ∈ h0,Tivyjma koneˇcné množiny,

uj(t)∈ Uj pro každét ∈ h0,Ti,j =1, . . . ,m, Z T

0

F(t,y(t),u(t))dt je maximální.

(88)

Hledámey = [y1, . . . ,yn], kde složky jsou po ˇcástech diferencovatelné funkce na intervaluh0,Ti, a funkci u= [u1, . . . ,um]s po ˇcástech spojitými složkami nah0,Ti takové, že

y(0) = y0,

y0(t) =f(t,y(t),u(t))pro každét ∈ h0,Tivyjma koneˇcné množiny,

uj(t)∈ Uj pro každét ∈ h0,Ti,j =1, . . . ,m, Z T

0

F(t,y(t),u(t))dt je maximální.

(89)

Hledámey = [y1, . . . ,yn], kde složky jsou po ˇcástech diferencovatelné funkce na intervaluh0,Ti, a funkci u= [u1, . . . ,um]s po ˇcástech spojitými složkami nah0,Ti takové, že

y(0) = y0,

y0(t) =f(t,y(t),u(t))pro každét ∈ h0,Tivyjma koneˇcné množiny,

uj(t)∈ Uj pro každét ∈ h0,Ti,j =1, . . . ,m, Z T

0

F(t,y(t),u(t))dt je maximální.

(90)

Hledámey = [y1, . . . ,yn], kde složky jsou po ˇcástech diferencovatelné funkce na intervaluh0,Ti, a funkci u= [u1, . . . ,um]s po ˇcástech spojitými složkami nah0,Ti takové, že

y(0) = y0,

y0(t) =f(t,y(t),u(t))pro každét ∈ h0,Tivyjma koneˇcné množiny,

uj(t)∈ Uj pro každét ∈ h0,Ti,j =1, . . . ,m, Z T

0

F(t,y(t),u(t))dt je maximální.

(91)

Hledámey = [y1, . . . ,yn], kde složky jsou po ˇcástech diferencovatelné funkce na intervaluh0,Ti, a funkci u= [u1, . . . ,um]s po ˇcástech spojitými složkami nah0,Ti takové, že

y(0) = y0,

y0(t) =f(t,y(t),u(t))pro každét ∈ h0,Tivyjma koneˇcné množiny,

uj(t)∈ Uj pro každét ∈ h0,Ti,j =1, . . . ,m,

Z T

0

F(t,y(t),u(t))dt je maximální.

(92)

Hledámey = [y1, . . . ,yn], kde složky jsou po ˇcástech diferencovatelné funkce na intervaluh0,Ti, a funkci u= [u1, . . . ,um]s po ˇcástech spojitými složkami nah0,Ti takové, že

y(0) = y0,

y0(t) =f(t,y(t),u(t))pro každét ∈ h0,Tivyjma koneˇcné množiny,

uj(t)∈ Uj pro každét ∈ h0,Ti,j =1, . . . ,m, Z T

0

F(t,y(t),u(t))dt je maximální.

(93)

V ˇeta 4.3 (Pontrjagin ˚uv princip maxima pro (P4’))

Necht’ vektorová funkceu = (u1, . . . ,um)je bodem maxima v úloze (P4’).

Pak existuje vektorová funkce λ:h0,Ti →Rn, že pro hamiltonián

H(t,y,u,λ) = F(t,y,u) +λTf(t,y,u) platí:

(94)

V ˇeta 4.3 (Pontrjagin ˚uv princip maxima pro (P4’))

Necht’ vektorová funkceu = (u1, . . . ,um)je bodem maxima v úloze (P4’). Pak existuje vektorová funkce λ:h0,Ti →Rn, že pro hamiltonián

H(t,y,u,λ) = F(t,y,u) +λTf(t,y,u) platí:

(95)

(PM1) pro každét ∈ h0,Tivyjma koneˇcné množiny a pro každéu˜ ∈ U1× · · · × Um platí

H(t,y(t),u(t),λ(t))≥H(t,y(t),u,˜ λ(t)),

(PM2) yi0 = ∂H

∂λi, i =1, . . . ,n(stavová rovnice), (PM3) λ0i =−∂H

∂yi

, i =1, . . . ,n (pohybová rovnice), (PM4) λi(T) =0, i =1, . . . ,n (podmínka transverzality).

(96)

(PM1) pro každét ∈ h0,Tivyjma koneˇcné množiny a pro každéu˜ ∈ U1× · · · × Um platí

H(t,y(t),u(t),λ(t))≥H(t,y(t),u,˜ λ(t)),

(PM2) yi0 = ∂H

∂λi, i =1, . . . ,n(stavová rovnice),

(PM3) λ0i =−∂H

∂yi

, i =1, . . . ,n (pohybová rovnice), (PM4) λi(T) =0, i =1, . . . ,n (podmínka transverzality).

(97)

(PM1) pro každét ∈ h0,Tivyjma koneˇcné množiny a pro každéu˜ ∈ U1× · · · × Um platí

H(t,y(t),u(t),λ(t))≥H(t,y(t),u,˜ λ(t)),

(PM2) yi0 = ∂H

∂λi, i =1, . . . ,n(stavová rovnice), (PM3) λ0i =−∂H

∂yi

, i =1, . . . ,n (pohybová rovnice),

(PM4) λi(T) =0, i =1, . . . ,n (podmínka transverzality).

(98)

(PM1) pro každét ∈ h0,Tivyjma koneˇcné množiny a pro každéu˜ ∈ U1× · · · × Um platí

H(t,y(t),u(t),λ(t))≥H(t,y(t),u,˜ λ(t)),

(PM2) yi0 = ∂H

∂λi, i =1, . . . ,n(stavová rovnice), (PM3) λ0i =−∂H

∂yi

, i =1, . . . ,n (pohybová rovnice), (PM4) λi(T) =0, i =1, . . . ,n (podmínka transverzality).

(99)

(PM1) pro každét ∈ h0,Tivyjma koneˇcné množiny a pro každéu˜ ∈ U1× · · · × Um platí

H(t,y(t),u(t),λ(t))≥H(t,y(t),u,˜ λ(t)),

(PM2) yi0 = ∂H

∂λi, i =1, . . . ,n(stavová rovnice), (PM3) λ0i =−∂H

∂yi

, i =1, . . . ,n (pohybová rovnice), (PM4) λi(T) =0, i =1, . . . ,n (podmínka transverzality).

Odkazy

Související dokumenty

∗ Zadaný rozsah stran je platný: Zobrazit název periodika, jeho ISSN kód, ˇcíslo roˇcníku aktuálního výtisku, datum vydání tohoto roˇcníku (letopoˇcet),

Diagonální matice D vyhovující rovnici (10.1) musí mít stejná vlasní čísla jako matice A. Tato vlastní čísla tedy musí ležet na diagonále matice D. Algebraická

Kdyby ne, tak vzhledem k tomu, že B je typu (n, n), musí být na diagonále matice B nenulové prvky, protože každý následující řádek má na začátku alespoň o jednu nulu

Pak A je

[r]

[r]

Toníkovo ˇcíslo je dvakrát vˇetší než Kájovo.. Toníkovo ˇcíslo je tˇrikrát vˇetší

Každé následující písmeno pak vypisuje do políˇcka, které má s políˇckem, do kterého bylo vepsáno pˇredchozí písmeno, alespo ˇn jeden spoleˇcný bod.. V Kocourkovˇe