18.05.2014 Views

PROJETO DE ALGORITMOS PARA RESOLU¸C˜AO DE ...

PROJETO DE ALGORITMOS PARA RESOLU¸C˜AO DE ...

PROJETO DE ALGORITMOS PARA RESOLU¸C˜AO DE ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

O método LOC será utilizado para otimizar o cálculo de α k .<br />

Seguindo a abordagem utilizada no capítulo de otimização convexa, será escolhida<br />

como função de energia<br />

E(x) := σ(x)‖f(x)‖ 2 + g T (x)Λuhsgn(g(x)) + h T (x)Γsgn(h(x)) (5.77)<br />

a qual apresenta um lado direito com descontinuidades sobre a fronteira do conjunto<br />

viável.<br />

Para fazer o cálculo do passo ótimo de iteração serão consideradas as trajetórias<br />

dentro e fora do conjunto viável.<br />

5.8.1.1 Fase de alcançar o conjunto viável<br />

Nesta fase x k /∈ χ.<br />

Fora do conjunto viável, o primeiro termo do lado direito de (5.77) se anula por ser o<br />

fator σ(x k ) = 0, resultando a função de energia nos termos de penalização. Escolhendo<br />

como resíduo nesta fase, esta função de energia:<br />

r k := E(x k ) = g T (x k )Λuhsgn(g(x k )) + h T (x k )Γsgn(h(x k ))<br />

A lei de atualização das variáveis (5.75) nesta fase é u k = −∇E(x k ), onde, similarmente<br />

ao que acontece no caso contínuo, se anulam os termos que contêm as derivadas<br />

das funções descontínuas.<br />

Calculando o resíduo na iteração seguinte, aproximando por Taylor até o termo de<br />

primeira ordem e utilizando (5.76):<br />

r k+1 ≃ r k + ∇ T E(x k )(x k+1 − x k ) = E(x k ) − α k ∇ T E(x k )∇E(x k )<br />

Nesta fase, a função de Liapunov candidata discreta será escolhida como o quadrado<br />

do resíduo a ser minimizado:<br />

V k = rk<br />

2<br />

⇒ ∆V k = V k+1 − V k = rk+1 2 − r2 k ≃<br />

−2α k E(x k )‖∇E(x k )‖ 2 2 + αk 2‖∇E(x k)‖ 4 2<br />

243

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

Saved successfully!

Ooh no, something went wrong!