31.01.2014 Aufrufe

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

MEHR ANZEIGEN
WENIGER ANZEIGEN

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

Hurra! Ihre Datei wurde hochgeladen und ist bereit für die Veröffentlichung.

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!