Numerical analysis of time discretization of optimal control problems

More documents

Recommendations

Info

Partitioned Euler integrators • Consider the partitioned ODE ˙y = F(y,p); ˙p = G(y,p) • Associated partitioned Euler scheme: yk+1−yk hk = F(yk,pk+1); pk+1−pk hk = G(yk,pk+1); • Here: F(y,p) = f(Υ(y,p),y); G(y,p) = −pfy(Υ(y,p),y). • It can easily be checked that S h → S locally uniformly as well as its derivatives. 21
Error analysis • F: set of C 1 mappings R n → R n in ¯ B(¯p0,1). • Define Ξ : R n ×F → R n , (p0,F) ↦→ Ξ(p0,F) := F(p0). • Clearly Ξ is C 1 and ∂Ξ(p0,F) ∂p0 = F ′ (p0). • If invertible, we can apply the IFT (Banach space setting) to Ξ at (¯p0,S). • Conclusion: there exists a locally unique ¯p h 0 solution of S h (¯p h 0) = 0 • Error estimate: |¯p h 0 − ¯p0| = O(S h −S) = O( ¯ h). More precisely: |¯p h 0 − ¯p0| = O(S h (¯p0)|). 22
Page 1 and 2: Numerical analysis of time discreti
Page 3 and 4: This class of problems is quite lar
Page 5 and 6: The Euler discretization • N: num
Page 7 and 8: The simplest optimal time problem I
Page 9 and 10: Fuller’s problem I (work with J.
Page 11 and 12: Robbins state constrained problem I
Page 13 and 14: Robbins state constrained problem I
Page 15 and 16: Beam problem: a second-order state
Page 17 and 18: Discretization by Euler’s method
Page 19 and 20: Reduction by elimination of the con
Page 21: Shooting problem: discretized probl
Page 25 and 26: Second-order optimality conditions
Page 27 and 28: Link with the shooting formulation
Page 29 and 30: A priori parameterized control •
Page 31 and 32: Example: mid-point rule • If uk c
Page 33 and 34: Order computation “by hand” Den
Page 35 and 36: Order of a RK scheme I • We forma
Page 37 and 38: Exercice: order 2 yk +hyk+1,1+ 1 2h
Page 39 and 40: s Lagrangian (contracting yk +hk j=
Page 41 and 42: Partitionned Runge-Kutta (PRK) sche
Page 43 and 44: Conditions for order 1 to 3: dj :=
Page 45 and 46: Conditions for order 5 Table 5: Ord
Page 47 and 48: Number of conditions for each order
Page 49 and 50: Definition 1. (i) A linear mapping
Page 51 and 52: Back to the mid-point rule • The
Page 53 and 54: Synthesis of part II: unconstrained
Page 55 and 56: III: CONTROL CONSTRAINTS (Hager, Do
Page 57 and 58: Active constraints Denote the set o
Page 59 and 60: Minimization of the Hamiltonian •
Page 65 and 66: Euler discretization ⎧ Min φ(yN)
Page 67 and 68: Geometrical hypotheses: junction po
Page 69 and 70: Junction points for the discretized
Page 71 and 72: Connection with stability analysis
Page 73 and 74:
Connection with stability analysis
Page 75 and 76:
High-order RK schemes and control c
Page 77 and 78:
Abstract stability result (in view
Page 79 and 80:
Since F1 is strongly coercive, it f
Page 81 and 82:
Examples: high orders are natural
Page 83 and 84:
Academic example: first-order state
Page 85 and 86:
Academic example: first boundary ar
Page 87 and 88:
Academic example: three boundary ar
Page 89 and 90:
First-order state constraint: Min
Page 91 and 92:
Third-order state constraint: Min
Page 93 and 94:
In the sequel: we discuss only firs
Page 95 and 96:
Condition at junction point Case of
Page 97 and 98:
Homotopy path For θ ∈ [0,1]: Min
Page 99 and 100:
Proof. We have that 0 = pθ k+1fu(u
Page 101 and 102:
and similarly 0 ≥ gi(y θ k−1)
Page 103 and 104:
For ε > 0 small enough, we deduce
Page 105 and 106:
OPEN PROBLEMS • First-order state
Page 107 and 108:
Sources • J.F. Bonnans, J. Lauren
Page 109:
The End ! 108
show all

Numerical analysis of time discretization of optimal control problems

Create successful ePaper yourself

Delete template?

Save as template?