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 operador de projeção ortogonal que, como já foi apontado, pode ser computacionalmente<br />

caro de calcular.<br />

Uma versão preliminar deste algoritmo foi publicada em [60].<br />

No contexto da discretização do algoritmo contínuo, o controle ótimo de Liapunov é<br />

aplicado no cálculo do comprimento do passo a fim de minimizar a norma quadrado de<br />

um resíduo ou erro. A maior dificuldade neste processo consiste na escolha adequada<br />

desse resíduo. Efetivamente, o problema de desigualdade variacional pode ser planteado<br />

como um problema de minimização da norma de um resíduo, o qual poderia ser definido<br />

das seguintes maneiras:<br />

a) r(x) = x − x ∗<br />

b) r(x) = ‖x − x ∗ ‖<br />

c) r(x) = (x − x ∗ ) T f(x)<br />

d) r(x) = e(x) = x − Pr χ [x − f(x)]<br />

No caso c), se a função objetivo satifazer (5.14), r(x) = 0 se e somente se x = x ∗ .<br />

O desconhecimento a priori do ponto ótimo x ∗ impossibilita a utilização dos três<br />

primeiros resíduos no cálculo do comprimento do passo. No caso do quarto resíduo,<br />

seu cálculo exige a utilização da projeção ortogonal, o que pretendemos evitar aqui.<br />

Tentaremos, portanto, escolher um outro valor de resíduo que não apresente os<br />

problemas mencionados.<br />

5.8.1 Cálculo do comprimento do passo por controle ótimo de<br />

Liapunov<br />

Seja a lei de atualização das variáveis<br />

u k := − [ σ(x k )f(x k ) + D T g (x k )Λuhsgn(g(x k )) + D T h(x k )Γsgn(h(x k )) ] (5.75)<br />

Então, a discretização de (5.71) por Euler, com comprimento do passo α k pode ser<br />

escrita:<br />

x k+1 = x k + α k u k (5.76)<br />

242

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

Saved successfully!

Ooh no, something went wrong!