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.

3.4.1 Simulações dos algoritmos contínuos de segunda ordem<br />

para achar mínimos de funções escalares<br />

A seguir, serão realizadas simulações com os algoritmos apresentados na seção anterior.<br />

Foram escolhidas, a fim de poderem ser graficadas, funções de duas variáveis. As<br />

simulações foram realizadas com o aplicativo Simulink do programa MATLAB 6. A<br />

discretização foi realizada pelo método de Euler, com comprimento do passo fixo igual<br />

a 0.01s. e foram executados durante 10s.<br />

Em todos os casos, o valor inicial da variável de controle foi u 0 = [0 0] T . Os ganhos<br />

foram variáveis segundo a função e o algoritmo.<br />

3.4.1.1 Simulações dos algoritmos de segunda ordem contínuos com a função<br />

de Camelback invertida<br />

A função de Camelback (equação 3.8), foi invertida a fim de apresentar mínimos locais<br />

e globais.<br />

Foram escolhidas como constantes a = −2, b = 1.05, c = − 1, d = −1 e e =<br />

6<br />

0. Com estas constantes a função apresenta um mínimo global em x = [0 0] T , dois<br />

mínimos locais em x = [−1.7475 0.8737] T e x = [1.7475 − 0.8737] T , e duas selas em<br />

x = [−1.0705 0.5352] T e x = [1.0705 − 0.5352] T .<br />

Foram experimentados diversos ganhos e escolhidos aqueles que apresentaram melhores<br />

resultados, no sentido das trajetórias convergirem mais diretamente para um<br />

mínimo, apresentar menores oscilações, e, principalmente, a fim de conseguir ultrapassar<br />

os mínimos locais para se estacionarem no mínimo global. Estes foram<br />

MI1: κ = 1<br />

MI3: Γ = 8I<br />

MI4: Γ = 0.1I<br />

MI5: Γ = 0.1I, a = b = 1<br />

HBF: Γ = I<br />

DIN: a = b = 1<br />

No algortimo MI2, o ganho é (∇ 2 φ(x)) 2 , razão pela qual este ganho não pode ser<br />

escolhido.<br />

Os pontos iniciais testados foram x 0 = [−1 − 1.5] T e x 0 = [2 1] T .<br />

145

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

Saved successfully!

Ooh no, something went wrong!