disser1.pdf (2006 KB) - Ernst-Moritz-Arndt-Universität Greifswald
disser1.pdf (2006 KB) - Ernst-Moritz-Arndt-Universität Greifswald
disser1.pdf (2006 KB) - Ernst-Moritz-Arndt-Universität Greifswald
Sie wollen auch ein ePaper? Erhöhen Sie die Reichweite Ihrer Titel.
YUMPU macht aus Druck-PDFs automatisch weboptimierte ePaper, die Google liebt.
5.5. Ökologische Modelle mit nichtlinear auftretender Steuerung<br />
Nun lösen wir die Prozessgleichungen für konstante v k :<br />
ẋ 1 (t, v k ) = (ε 1 − s 1 v k )x 1 (t, v k ), k = 0, 1<br />
mit entsprechenden Anfangsbedingungen x 1 (t k ) = x 1k , k = 0, 1 und erhalten:<br />
x 1 (t, v k ) = e (ε 1−s 1 v k )t , t ∈ [t k , t k+1 ),<br />
mit x 1 (0, v 0 ) = 0.25. Nun berechnen wir λ 1 (t, v k ) und S 1 (t, v k ) für k = 0, 1 zu:<br />
λ 1 (t, v k ) = −<br />
p 1v k<br />
+ L 1 · e (s 1v k −ε 1 )t ,<br />
ε 1 − sv k<br />
S 1 (t, v k ) = c 1 v 2 k t + L 2 .<br />
Dabei sind L 1 , L 2 Integrationskonstanten. Auf dem Intervall [10; 20) erhalten wir mit<br />
Hilfe der Transversalitätsbedingungen<br />
λ 1 (t, v 1 ) = − p 1v 1<br />
+ p 1v 1<br />
· e (s 1v 1 −ε 1 )(t−20) ,<br />
ε 1 − sv 1 ε 1 − sv 1<br />
S 1 (t, v 1 ) = c 1 v 2 1 t − 20c 1 v 2 1 .<br />
Durch Übergangsbedingungen resultieren im Intervall [0; 10):<br />
λ 1 (t, v 0 ) = − p (<br />
1v 0 p1 v 0<br />
+ − p )<br />
1v 1<br />
· e (s 1v 0 −ε 1 )(t−10)<br />
ε 1 − sv 0 ε 1 − sv 0 ε 1 − sv 1<br />
+ p 1v 1<br />
ε 1 − sv 1<br />
· e −10·(s 1v 0 +s 1 v 1 −2ε 1) · e (s 1v 0 −ε 1 )t ,<br />
S 1 (t, v 0 ) =c 1 v 0 2 t − 10c 1 (v 0 2 + v 1 2 ).<br />
Die optimalen Steuerungen werden mit Hilfe folgender Formeln berechnet:<br />
{<br />
u ∗ 1(0) = arg max − 10c 1 (v 2 0 + v 2 1 )<br />
v 0 ∈Ω<br />
[<br />
+ − p (<br />
1v 0 p1 v 0<br />
+<br />
− p )<br />
1v 1<br />
· e (s 1v 0 −ε 1 )(−10)<br />
ε 1 − s 1 v 0 ε 1 − s 1 v 0 ε 1 − s 1 v 1<br />
+ p ] }<br />
1v 1<br />
· e −10·(s 1v 0 +s 1 v 1 −2ε 1 )<br />
· 0.25 ,<br />
ε 1 − s 1 v<br />
{ 1<br />
u ∗ 2<br />
1(10) = arg max − 10c 1 v 1<br />
v 1 ∈Ω<br />
[<br />
+ − p 1v 1<br />
+ p ]<br />
}<br />
1v 1<br />
· e (s 1v 1 −ε 1 )(−10)<br />
· 0.25 · e (ε 1−s 1 v 0 )·10<br />
.<br />
ε 1 − sv 1 ε 1 − sv 1<br />
111