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
A. DIRCOL C* to the interval [0,1].) C* C* EXIT-Parameters: C* ----------- C* X........initial estimates of the state variables at TAU C* (real array) C* U........initial estimates of the control variables at TAU C* (real array) C* C* If IPHASE = 0 then: C* ================== C* ENTRY-Parameters: C* ----------- C* LU.......the number of events (integer) C* C* EXIT-Parameters: C* ----------- C* U........initial estimates of the events E(1), ..., E(M) C* (real array) C* C* If IPHASE < 0 then: C* ================== C* ENTRY-Parameters: C* ----------- C* LU.......the number LP of the control parameters (integer) C* C* EXIT-Parameters: C* ----------- C* U........initial estimates of the control parmeters C* P(1), ..., P(LP) (real array) C***************************************** C IMPLICIT NONE C INTEGER IPHASE, NX, LU C**** REAL DOUBLE PRECISION + TAU, X(NX), U(LU) C 126
A.1. Programm 1. Hauptprogramm user.f (wesentliche Teile) C-----BEGIN--PROBLEM C IF (IPHASE .GT. 0) THEN C C ---- Initial estimates of X and U at TAU C X(1) = 0.25D0 X(2) = 0.8D0 X(3) = 1.0D0 X(4) = 0D0 U(1) = 0D0 U(2) = 0D0 U(3) = 0D0 C ELSE IF (IPHASE .EQ. 0) THEN C C ---- Initial estimates of the events E C U(1) = 0.D0 U(2) = 0.1D0 U(3) = 0.D0 C ELSE IF (IPHASE .LT. 0) THEN C C ---- Initial estimates of the control parameters P C C There are none. C END IF C C------END---PROBLEM- C RETURN C -- End of subroutine USRSTV END SUBROUTINE USROBJ( NR, NX, LU, LP, ENR, XL, UL, P, FOBJ, XR, UR) C* Purpose: C* ----- C* This subroutine provides the objective of the optimal C* control problem in a Mayer type form. C* The objective of a multi-phase program is assumed to be 127
- 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 86 und 87: 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: ANHANG A DIRCOL A.1. Programm 1. Ha
- 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 *
- Seite 137 und 138: ANHANG B OC-ODE B.1. Programm 1. Be
- Seite 139 und 140: B.1. Programm 1. Berechnung mit st
- Seite 141 und 142: B.1. Programm 1. Berechnung mit st
- Seite 143 und 144: B.1. Programm 1. Berechnung mit st
- Seite 145 und 146: ANHANG C MATLAB-Programme C.1. Prog
- Seite 147 und 148: C.2. Programm 2. Berechnung der Adj
- Seite 149 und 150: C.4. Programm 4. Das Integral-Maxim
- Seite 151 und 152: C.6. Programm 6. Newton-Methode fü
- Seite 153 und 154: Danksagung An dieser Stelle möchte
- Seite 155: Versicherung Hiermit erkläre ich,
A. DIRCOL<br />
C* to the interval [0,1].)<br />
C*<br />
C* EXIT-Parameters:<br />
C* -----------<br />
C* X........initial estimates of the state variables at<br />
TAU<br />
C* (real array)<br />
C* U........initial estimates of the control variables at<br />
TAU<br />
C* (real array)<br />
C*<br />
C* If IPHASE = 0 then:<br />
C* ==================<br />
C* ENTRY-Parameters:<br />
C* -----------<br />
C* LU.......the number of events (integer)<br />
C*<br />
C* EXIT-Parameters:<br />
C* -----------<br />
C* U........initial estimates of the events E(1), ...,<br />
E(M)<br />
C* (real array)<br />
C*<br />
C* If IPHASE < 0 then:<br />
C* ==================<br />
C* ENTRY-Parameters:<br />
C* -----------<br />
C* LU.......the number LP of the control parameters<br />
(integer)<br />
C*<br />
C* EXIT-Parameters:<br />
C* -----------<br />
C* U........initial estimates of the control parmeters<br />
C* P(1), ..., P(LP) (real array)<br />
C*****************************************<br />
C<br />
IMPLICIT NONE<br />
C<br />
INTEGER IPHASE, NX, LU<br />
C**** REAL<br />
DOUBLE PRECISION<br />
+ TAU, X(NX), U(LU)<br />
C<br />
126