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
5. Numerische Methoden für verschiedene Aufgaben der Optimalsteuerung u app t k t k+1 t x app x(t k+1/2) t k t k+1/2 t k+1 t Abbildung 5.2.: Stückweise lineare Approximation der Steuerung (oben) und stückweise kubische Approximation der Trajektorie (unten) Nutzt man die Stetigkeitsbedingungen und die Differenzierbarkeit an den Punkten t k , k = 1, N − 1, so bekommt man die Unbekannten c q (k) als Funktionen von x(t k ), x(t k+1 ), u(t k ) und u(t k+1 ): c 0 (k) = x(t k ), c 1 (k) = h k f k , c 2 (k) = 3[x(t k+1 ) − x(t k )] − h k [f k+1 + 2f k ], c 3 (k) = −2[x(t k+1 ) − x(t k )] + h k [f k+1 + f k ]. Dabei ist f k = f(x(t k ), u(t k ), t k ), k = 1, . . . , N − 1. 86
5.2. Direkte Methoden Weiterhin muss die Kollokationsbedingung: f(x app (t k+ 1 2 ), u app (t k+ 1 2 ), t k+ 1 2 ) − ẋ app (t k+ 1 ) = 0 2 für k = 1, . . . , N − 1 gelten. Hier wird der Ausdruck ẋ app (t k+ 1 ) durch 2 ẋ app (t k+ 1 ) = c 1 (k) + 2c 2 (k) ( ) t 2 h k h 2 k+ 1 − t k + 3c 3 (k) ( ) 2 t 2 k h 3 k+ 1 − t k 2 k = c 1 (k) h k = c 1 (k) h k + 2c 2 (k) h 2 k + c 2 (k) h k · hk 2 + 3c 3 (k) · h2 k h 3 k 4 + 3c 3 (k) 4h k = h k f k h k + 3[x(t k+1) − x(t k )] − h k [f k+1 + 2f k ] h k + −6[x(t k+1) − x(t k )] + 3h k [f k+1 + f k ] 4h k = f k + 3[x(t k+1) − x(t k )] 2h k − [f k+1 + 2f k ] + 3 4 · [f k+1 + f k ] = 3[x(t k+1) − x(t k )] − 1 2h k 4 · [f k+1 + f k ] ersetzt. Daraus resultiert das endlich-dimensionale nichtlineare Optimierungsproblem: Minimiere J(Y ) = J(x(t N )), mit Y = (x(t 1 ), u(t 1 ), . . . , x(t N ), u(t N )) T ∈ R N(n+m) , und Kollokationsbedingung: f(x app (t k+ 1 2 ), u app (t k+ 1 2 für k = 1, . . . , N − 1. ), t k+ 1 ) − 2 ( 3[x(tk+1 ) − x(t k )] − 1 ) 2h k 4 · [f k+1 + f k ] = 0, Werden bei einer Optimalsteuerungsaufgabe zusätzliche Nebenbedingungen gestellt, so werden diese Bedingungen an den diskreten Punkten t k , 1, . . . , N − 1 ausgewertet. Die Eingabe erfolgt über FORTRAN-Subroutinen USRSTV und USROBJ, die Dimensionen des Problems und jeweilige Zustands- und Steuerbeschränkungen werden durch DATDIM und DATLIM ermittelt (genauere Beschreibung im Anhang A.3 zu finden). 5.2.2. Direktes Schießverfahren Das direkte Schießverfahren diskretisiert nur die Steuerung und erzeugt mit der folgenden numerischen Integration der Differentialgleichung ein nichtlineares 87
- Seite 36 und 37: 2. Biologische, ökologische und wi
- Seite 38 und 39: 3. Mathematische Modelle in der Bio
- Seite 40 und 41: 3. Mathematische Modelle in der Bio
- Seite 42 und 43: 3. Mathematische Modelle in der Bio
- Seite 44 und 45: 3. Mathematische Modelle in der Bio
- Seite 46 und 47: 3. Mathematische Modelle in der Bio
- Seite 48 und 49: 3. Mathematische Modelle in der Bio
- Seite 50 und 51: 4. Grundlagen der Optimalsteuerung
- Seite 52 und 53: 4. Grundlagen der Optimalsteuerung
- Seite 54 und 55: 4. Grundlagen der Optimalsteuerung
- Seite 56 und 57: 4. Grundlagen der Optimalsteuerung
- Seite 58 und 59: 4. Grundlagen der Optimalsteuerung
- Seite 60 und 61: 4. Grundlagen der Optimalsteuerung
- Seite 62 und 63: 4. Grundlagen der Optimalsteuerung
- Seite 64 und 65: 4. Grundlagen der Optimalsteuerung
- Seite 66 und 67: 4. Grundlagen der Optimalsteuerung
- Seite 68 und 69: 4. Grundlagen der Optimalsteuerung
- Seite 70 und 71: 4. Grundlagen der Optimalsteuerung
- Seite 72 und 73: 4. Grundlagen der Optimalsteuerung
- Seite 74 und 75: 4. Grundlagen der Optimalsteuerung
- Seite 76 und 77: 4. Grundlagen der Optimalsteuerung
- Seite 78 und 79: 4. Grundlagen der Optimalsteuerung
- Seite 80 und 81: 4. Grundlagen der Optimalsteuerung
- Seite 82 und 83: 4. Grundlagen der Optimalsteuerung
- Seite 84 und 85: 5. Numerische Methoden für verschi
- Seite 88 und 89: 5. Numerische Methoden für verschi
- Seite 90 und 91: 5. Numerische Methoden für verschi
- Seite 92 und 93: 5. Numerische Methoden für verschi
- Seite 94 und 95: 5. Numerische Methoden für verschi
- Seite 96 und 97: 5. Numerische Methoden für verschi
- Seite 98 und 99: 5. Numerische Methoden für verschi
- Seite 100 und 101: 5. Numerische Methoden für verschi
- Seite 102 und 103: 5. Numerische Methoden für verschi
- Seite 104 und 105: 5. Numerische Methoden für verschi
- Seite 106 und 107: 5. Numerische Methoden für verschi
- Seite 108 und 109: 5. Numerische Methoden für verschi
- Seite 110 und 111: 5. Numerische Methoden für verschi
- Seite 112 und 113: 5. Numerische Methoden für verschi
- Seite 114 und 115: 5. Numerische Methoden für verschi
- Seite 116 und 117: 6. Zusammenfassung und Ausblick in
- Seite 118 und 119: 6. Zusammenfassung und Ausblick die
- Seite 120 und 121: Literaturverzeichnis [15] M. Herman
- Seite 122 und 123: Literaturverzeichnis [54] http://ww
- Seite 125 und 126: ANHANG A DIRCOL A.1. Programm 1. Ha
- Seite 127 und 128: A.1. Programm 1. Hauptprogramm user
- Seite 129 und 130: A.1. Programm 1. Hauptprogramm user
- Seite 131 und 132: A.2. Programm 2. Programm DATDIM *
- Seite 133 und 134: A.3. Programm 3. Programm DATLIM 0
- Seite 135 und 136: A.3. Programm 3. Programm DATLIM *
5. Numerische Methoden für verschiedene Aufgaben der Optimalsteuerung<br />
u app<br />
t k<br />
t k+1<br />
t<br />
x app<br />
x(t k+1/2)<br />
t k<br />
t k+1/2<br />
t k+1<br />
t<br />
Abbildung 5.2.: Stückweise lineare Approximation der Steuerung (oben) und stückweise<br />
kubische Approximation der Trajektorie (unten)<br />
Nutzt man die Stetigkeitsbedingungen und die Differenzierbarkeit an den Punkten<br />
t k , k = 1, N − 1, so bekommt man die Unbekannten c q (k) als Funktionen von x(t k ),<br />
x(t k+1 ), u(t k ) und u(t k+1 ):<br />
c 0 (k) = x(t k ),<br />
c 1 (k) = h k f k ,<br />
c 2 (k) = 3[x(t k+1 ) − x(t k )] − h k [f k+1 + 2f k ],<br />
c 3 (k) = −2[x(t k+1 ) − x(t k )] + h k [f k+1 + f k ].<br />
Dabei ist f k = f(x(t k ), u(t k ), t k ), k = 1, . . . , N − 1.<br />
86