12.07.2015 Views

a variacional del estado de transición a la - Páxinas persoais - USC ...

a variacional del estado de transición a la - Páxinas persoais - USC ...

a variacional del estado de transición a la - Páxinas persoais - USC ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

A. Método simplex <strong>de</strong>scen<strong>de</strong>nte multidimensional 105por el valle”. Si se da <strong>la</strong> situación en <strong>la</strong> que el simplex “trata <strong>de</strong> pasar a través<strong><strong>de</strong>l</strong> ojo <strong>de</strong> una aguja”, se autocontrae en todas <strong>la</strong>s direcciones cercando elpunto en el que el valor <strong>de</strong> <strong>la</strong> función es más bajo. Estos movimientos básicosaparecen reflejados en <strong>la</strong> Figura A.1.El algoritmo <strong><strong>de</strong>l</strong> método simplex <strong>de</strong>scen<strong>de</strong>nte repite los siguientes pasoshasta obtener el valor mínimo <strong>de</strong> <strong>la</strong> función:1. Or<strong>de</strong>nar y retiquetar los N+1 puntos <strong>de</strong> tal manera que f(x N ) > . . . >f(x 2 ) > f(x 1 ).2. Generar un punto prueba x r por reflexión, tal que x r = x g + a · (x g −x N+1 ) dón<strong>de</strong> x g es el centroi<strong>de</strong> <strong>de</strong> hiperfase <strong>de</strong> los N mejores puntos enlos vértices <strong><strong>de</strong>l</strong> simplex..3. Si f(x 1 ) < f(x r ) < f(x N ) entonces se reemp<strong>la</strong>za x N+1 por x r .4. Si f(x r ) < f(x 1 )) se genera un nuevo punto x e mediante expansión talque x e = x r + b · (x r − x g ) .5. Si f(x e ) < f(x r ), se reemp<strong>la</strong>za un nuevo punto x N+1 por x e , en otrocaso se reemp<strong>la</strong>za x N+1 por x r .6. Si f(x r ) > f(x N ), se genera un nuevo punto x C mediante contracción,tal que x e = x g + g · (x N+1− x g ).7. Si f(x C ) < f(x N+1 ) se reemp<strong>la</strong>za x N+1 por x e .8. Si f(x C ) > f(x N+1 ) se contrae a lo <strong>la</strong>rgo <strong>de</strong> todas <strong>la</strong>s direcciones hastax 1 . Se evalúan los valores <strong>de</strong> <strong>la</strong> función objetivo en los N nuevos vértices,x i = x 1 + h · (x i − x 1 ).Geométricamente, el método simplex se pue<strong>de</strong> interpretar <strong>de</strong> <strong>la</strong> siguientemanera: partiendo <strong>de</strong> uno <strong>de</strong> los vértices <strong><strong>de</strong>l</strong> poliedro, ir <strong>de</strong>sp<strong>la</strong>zándose a uno<strong>de</strong> los vértices adyacentes mejorando el valor <strong>de</strong> <strong>la</strong> función objetivo hastaencontrar el óptimo.Los criterios <strong>de</strong> convergencia para el fin <strong><strong>de</strong>l</strong> proceso pue<strong>de</strong>n ser tan pequeñoscomo se <strong>de</strong>seen, como en cualquier método <strong>de</strong> minimización multidimensiona<strong>la</strong>unque no se tiene <strong>la</strong> opción <strong>de</strong> exigir cierto grado <strong>de</strong> precisión parauna so<strong>la</strong> <strong>de</strong> <strong>la</strong>s variables in<strong>de</strong>pendientes. A este respecto es posible terminarcuando <strong>la</strong> distancia <strong><strong>de</strong>l</strong> vector por el que se mueve en unos <strong>de</strong> los ciclos estan pequeño como un parámetro <strong>de</strong> tolerancia prefijado. Alternativamente sepodría exigir una disminución en el valor <strong>de</strong> <strong>la</strong> función en el paso final que seatan pequeño como se <strong>de</strong>see.

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

Saved successfully!

Ooh no, something went wrong!