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.7 Gradients<br />

The Euler-Lagrange equations are propagated backwards to obtain the<br />

adjoint variables,<br />

⎛<br />

T<br />

⎛ ∂ f ⎞ ⎞<br />

k ∈ −<br />

λ<br />

⎜ ⎜ X ⎟ ⎟<br />

⎝ ⎝ ∂ ⎠ k ⎠<br />

X<br />

[ N −1<br />

( 1 ) n + 1 ] ⇒ λk<br />

1 : = ⎜ I + ∆t<br />

⋅ ⎜ ⎟ ⎟ ⋅ k<br />

The elements of the reduced gradient,<br />

g<br />

≡<br />

⎛ ∂ H ⎞<br />

⎜<br />

U<br />

⎟<br />

⎝ ∂ ⎠<br />

T<br />

: =<br />

⎛<br />

⎜ ∂ φ ∂ φ ⎛ ∂ f X ⎞<br />

−<br />

⎜ U X<br />

⎜<br />

X<br />

⎟<br />

∂ ∂<br />

⎝<br />

⎝ ∂ ⎠<br />

7-15<br />

−1<br />

∂ f X ⋅<br />

∂ U<br />

⎛ ∂ H ⎞<br />

⎛ ∂ f ⎞<br />

k ∈<br />

⎜<br />

λ<br />

U ⎟<br />

⎜<br />

k<br />

U ⎟<br />

⎝ ∂ ⎠<br />

⎝ ∂ k ⎠<br />

Equation 7.6-7<br />

⎞<br />

⎟<br />

⎟<br />

⎠<br />

T<br />

Equation 7.7-1<br />

k<br />

[ n ( 1 ) N −1<br />

] ⇒ ⎜ ⎟ : = ∆t<br />

⋅ ⎜ ⎟ ⋅ k<br />

∂ H<br />

∂ t<br />

N<br />

⎛ ∂ H ⎞<br />

⎜<br />

X ⎟<br />

⎝ ∂ 0 ⎠<br />

: =<br />

∂ φ<br />

∂ t<br />

N<br />

T<br />

+<br />

: =<br />

1<br />

N<br />

⋅<br />

T<br />

⎛<br />

⎜ ⎛ ∂ f 0 ⎞<br />

I + ∆t<br />

⋅<br />

⎜ ⎜<br />

X ⎟<br />

0<br />

⎝ ⎝ ∂ ⎠<br />

∑ − N 1<br />

k : = n<br />

⎛<br />

⎜<br />

⎜<br />

⎝<br />

f<br />

t<br />

k<br />

+ t<br />

n<br />

T<br />

⎛ ∂ f<br />

⋅ ⎜<br />

⎝ ∂ t<br />

⎞<br />

⎟<br />

⎟<br />

⎠<br />

k<br />

k<br />

⋅ λ<br />

⎞<br />

⎟<br />

⎠<br />

0<br />

T<br />

⎞<br />

⎟<br />

⎟<br />

⎠<br />

T<br />

Equation 7.7-2<br />

Equation 7.7-3<br />

⋅ λ<br />

k<br />

Equation 7.7-4<br />

Equation 7.7-3 is applicable when all the controls are available prior to<br />

launch and the launcher angles are optimised. If the gradient at a particular<br />

step results in a control exceeding the limit, the gradient for that control is<br />

set to zero, the control to its limit, and optimisation proceeds over the<br />

unbounded set of controls. The gradients with respect to optimiser states,<br />

controls and terminal time are provided by GRADIENT. This invokes<br />

OP_MS_DXDU for the variation in state derivatives with respect to the<br />

controls, and OP_MS_DXDX for their variation with respect to the states.<br />

The gradients are scaled in GRADIENT_SC, and the direction of steepest<br />

descent in STEEPEST_DD.

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

Saved successfully!

Ooh no, something went wrong!