Numerical analysis of time discretization of optimal control problems

More documents

Recommendations

Info

Order of a RK scheme II • By induction: explicit expansion, using q ℓ=0 h q q! ykiq = yk +h s j=1 aij q−1 ℓ=0 f(ykj(h))+O(h q+1 ). • Expansion: linear combination of the elementary differentials (as for the solution of the ODE) • Global error order p iff, in the expansion of yk+1, these coefficients coincide up to order p. 35
Exercice: order 2 yk +hyk+1,1+ 1 2h2yk+1,2 = yk +h s i=1bif(yk +hyki1))+O(h 3 ), yk +hyki1+ 1 2h2yki2 = yk +h s j=1aijf(yk +hyki1))+O(h 3 ) Using yki1 = s j=1 aijf(yk))+O(h 2 ) we deduce that and so 1 2 yki2 = s j=1 aijf ′ (yk)f(yk) = cif ′ (yk)f(yk), 1 2 yk+1,2 = s bicif ′ (yk)f(yk) i=1 Since yt = y0+tf + 1 2 t2 f ′ f, the scheme is (at least) of second order iff s i=1 bici = 1 2 . 36
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 and 22: Shooting problem: discretized probl
Page 23 and 24: Error analysis • F: set of C 1 ma
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: Order of a RK scheme I • We forma
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?