14.02.2013 Views

Thesis - Leigh Moody.pdf - Bad Request - Cranfield University

Thesis - Leigh Moody.pdf - Bad Request - Cranfield University

Thesis - Leigh Moody.pdf - Bad Request - Cranfield University

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Chapter 7 / Missile Trajectory Optimisation<br />

_ _<br />

7.6 TPBVP Formulation<br />

Minimising the cost function by adjusting the control sequence is replaced<br />

by the simpler problem of minimising a scalar Hamiltonian that depends on<br />

boundary conditions only. Thus the constrained TPBVP is reduced to the<br />

unconstrained minimisation of a Hamiltonian function. Pontryagin theory is<br />

used to reformulate the open-loop solution into a classical control law by<br />

relating the controls to state observer parameters, in this case using Equation<br />

7.3-4. The Hamiltonian is formed by combining the cost function, dynamic<br />

constraints using Lagrangian multipliers (λ), and inequality constraints<br />

using co-state variables (µ),<br />

H<br />

: =<br />

φ + λ<br />

T<br />

⋅ f<br />

7-14<br />

X<br />

+ µ<br />

T<br />

⋅<br />

( h + s ⊗ s − c )<br />

Equation 7.6-1<br />

The Kuhn-Tucker 1 st order conditions to be satisfied on an optimum<br />

trajectory,<br />

∂<br />

∂ X<br />

∂<br />

∂ U<br />

H T<br />

T<br />

: =<br />

: =<br />

∂ φ<br />

+<br />

∂ X<br />

λ<br />

∂ f X ⋅<br />

∂ X<br />

H T<br />

T<br />

∂<br />

∂ s<br />

∂ φ<br />

+<br />

∂ U<br />

∂ H<br />

∂ λ<br />

λ<br />

: =<br />

∂ f X ⋅<br />

∂ U<br />

f<br />

X<br />

H T<br />

: =<br />

0<br />

⇒<br />

µ<br />

+<br />

+<br />

: =<br />

µ<br />

µ<br />

0<br />

∂ h<br />

⋅<br />

∂ X<br />

∂ h<br />

⋅<br />

∂ U<br />

: =<br />

: =<br />

( c − h ) : = 0<br />

0<br />

Equation 7.6-2<br />

0<br />

Equation 7.6-3<br />

Equation 7.6-4<br />

Equation 7.6-5<br />

The optimiser states are propagated using Euler integration over the<br />

remaining controls [n(1)N-1], control limits applied, and the Hamiltonian<br />

Jacobians determined for each control. For proof of principle the discrete<br />

equations that follow involve only the dynamic equality constraints.<br />

Modifying the discrete equations in Vorley [V.4] to deal with the remaining<br />

controls, the transversality conditions at impact,<br />

λ N − 1<br />

( ) T<br />

∂ φ X<br />

: = ∂<br />

N<br />

Equation 7.6-6

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!