Ecuaciones no Lineales y Optimización | 6ta Parte
![](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmVrDxue83w7YemDc7ir68jCiLtP7CUhGBV6wQvetEo59x/Portada.png)
![](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmXGb9YQnbu6kDxJdRsfsepFWb5VGwuNNDRrTtk7FRjnii/Matematicas000.png)
La función vectorial
![Ecuación](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmenmbj7EFsB94BR1mZuC9DD7Bjn54PZTU4EMRds19Ytyi/Ecuacion01.png)
describe una órbita elíptica inclinada con un foco en (0, 0). Pensemos en el Sol como si estuviera situado en este punto. El parámetro φ1 es el ángulo de inclinación. Si A1 es mayor o igual que P1, entonces A1 y P1 son las distancias máximas y mínimas entre el Sol y la órbita. Sea
![Ecuación](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmarmh5B8DXPpJzmdWtKMt6c9HmeLM5AaKzZY97i4CSrHo/Ecuacion02.png)
una segunda órbita de este tipo y la cual podemos visualizar en la gráfica siguiente.
![Gráfica](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmShkv82dySJAmN4KjmWPP9kcUUVtka63WmF7j8WNTJN4r/Grafica01.png)
Gráfico de las órbitas. Elaborado con GNU Octave, por @abdulmath.
donde (A1, P1, φ1) = (10, 2, pi/8), y (A2, P2, φ2) = (10, 2, pi/8).
Nuestro objetivo es encontrar la distancia mínima desde un punto en la primera órbita (A1, P1, φ1) hasta un punto en la segunda órbita (A2, P2, φ2). Para una medida de distancia, utilizamos la siguiente función
![Ecuación](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmS3QAK2UNe3YCmaYiSqkn9qvKScxeLHx7Nb57QG9v4pB8/Ecuacion03.png)
Esta es una función de dos variables, y de acuerdo como lo vimos en la publiación anterior (Ecuaciones no Lineales y Optimización | 5ta Parte), la función S es la función objetivo. Notemos que t1 es un punto seleccionado en la órbita exterior y t2 es un punto seleccionado la órbita interior. A continuación estudiaremos la mínimización de funciones multivaluadas como S(t).
![](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmaoG5ZuMgFaYHr5CGwe3cA8GVEfH1edL4VwPT3fn1NCtb/SeparadorGMath38.png)
La implementación de S(t) depende de otra función que llamaremos orbita
, la cual mostramos a continuación:
![Script](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmVmmA3tKcFLjCZT25QqYiRqN2AZwVebjoznTWXFK4hKoJ/Script01.png)
Función
orbita
. Elaborada con GNU Octave, por @abdulmath.Con esta función podemos fácilmente generar y gráficas y puntos en las órbitas. El script siguiente gráfica dos órbitas y traza una línea que conecta el punto t1 = 3 en la primera órbita con el punto t2 = 4 en la segunda órbita.
![Script](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmabvVqEEeBz4UYeuBseUp4EQYEEun3CAPnRsexS2bHCrX/Script02.png)
Script usando para gráficas las órbitas. Elaborada con GNU Octave, por @abdulmath.
La separación entre estos dos puntos es una función de un vector de dos variables independientes y los parámetros A1, P1, φ1, A2, P2, φ2, el cual visualizamos en el siguiente funcion:
![Script](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmP9GS2jVCMTWB4fVtLjvQsiEkN22DKjhk2UzLnBFKmFGn/Script03.png)
Función
separacion
. Elaborado en GNU Octave, por @abdulmath.Observemos cómo se empaquetan los parámetros de órbita en un vector de dimensión 6, y usando el siguiente script a continuación podemos ver como usamos la función que acabamos de mostrar
separacion
![Script](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmUNJKRSRWNzfNsCMUpogCV72TCP7DvNF79xJyBkxz6W6N/Script04.png)
Script para calcular la separación de las órbitas. Elaborado en GNU Octave, por @abdulmath.
asigna a d el valor
![Ecuación](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmXiViRCjSbZNkUN2rgFUW75zTKTNGVRQdBirieG22NkLV/Ecuacion04.png)
con
![Ecuación](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmPTEH6C1DsuZzaNVktcWVMaQzPhWBtZJvwZe1PrKCPobS/Ecuacion05.png)
![](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmWc6PNwWG3TTTx9BQoTy4a8M7nVh88YMQAwVFjX6msEFD/SeparadorStem.png)
Supongamos que una persona se encuentra en un lugar nublado, y parada en una ladera sin un mapa. Si esta persona quiere dar un paso coherente con el objetivo de alcanzar de un valle, es razonable dar ese paso en la dirección más cuesta abajo. Matemáticamente, este hecho, dirigirse en la dirección del gradiente negativo.
Recordemos que el gradiente de una función
![Ecuación](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmc7vWMTHM9A46j1jWeY84AXxhmNSSCWAMoszjBhovcVCt/Ecuacion06.png)
es un vector de derivadas parciales:
![Ecuación](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmVqrs4rhg2xoWnyRrZMnXSuQPNfKwZHMUFDiAaXgjZc6x/Ecuacion07.png)
El gradiente siempre apunta en la dirección del mayor incremento y su negación apunta en la dirección de la mayor disminución. Tomando las parciales en t1 y t2 en la ecuación dada en (1), vemos que
![Ecuación](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmZBhZv3TZM3RshPcC2dns99w8VCwaUxUSbN7CbMwyuM4k/Ecuacion08.png)
Sustituyendo las definiciones de las funciones de las componentes, llegamos a la siguiente implementación del gradiente:
![Script05.png](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmdcW9QhfgEm6FhYvJk7pphcX5mq8hGV2V1iXFFVbk1xr6/Script05.png)
Función para calcular el gradiente de la función
separacion
. Elaborado con GNU Octave, por @abdulmath.La derivación no es importante. Pero lo que es importante es apreciar que la implementación de la función gradiente puede ser muy complicada y consumir mucho tiempo. Nuestro problema es pequeño (n = 2) y simple. Para aplicaciones de alta dimensión con funciones objetivas complicadas, a menudo es necesario usar algunos paquetes de diferenciación simbólica.
![](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmXffx8yjx3DgTBuid9akmfhwasc2QAEKLapRwGeorobFS/SeparadorGMath13.png)
A partir de la teoría de aproximación de Taylor multivalorada, sabemos que si λc es real y sc es un n-vector, entonces
![Ecuación](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmSaRrPvcrKBfmrAQZPDhegB89w3u6xjn2BXJvrvVkkn6T/Ecuacion09.png)
Por lo tanto, podemos esperar una disminución en el valor de f si establecemos que la dirección del paso sc sea
![Ecuación](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmTgShJHxcfic2yBiSJKPrY3SyG5ZkzVRykg3QJ4D4jYQt/Ecuacion10.png)
y el parámetro de longitud de paso λc para que sea lo suficientemente pequeño. Si el gradiente es cero, entonces este argumento se rompe. En ese caso tc es un punto crítico, pero puede no ser un mínimo local. Nuestra búsqueda del fondo del valle puede estancarse en ese lugar.
La determinación práctica de λc es el problema de la búsqueda de línea. Idealmente, nos gustaría minimizar f(tc + λcsc) sobre todo λ . Sin embargo, una solución aproximada a este problema de minimización unidimensional suele ser lo suficientemente buena y a menudo esencial si f es muy costosa de evaluar. Esto nos lleva al marco de descenso más empinado para calcular una mejora t+ al minimizador actual:
- calcular el gradiente
y sea sc = - gc. - Elegir a bajo costo λc tal que f(tc + λcsc) es suficientemente pequeña como f(tc).
- Sea t+ = tc + λcsc.
El diseño de una estrategia razonable de búsqueda de líneas es una tarea no trivial.
![](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmdpxUb7HtnNdmxthWLKfUjnPA5sind9336Fgw5iy7pfv3/SeparadorGMath12.png)
Queridos amigos y lectores, espero hayan disfrutado de esta segunda entrega de este tan apasionante e interesante tema del cálculo numérico, Ecuaciones no Lineales y Optimización | 6ta Parte apoyado con el entorno GNU Octave, en esta oportunidad pudimos estudiar algunas detalles o aspectos importantes, acerca de la minimización de funciones multivaloradas. Espero que el mismo sea de apoyo a ustedes en sus trabajos, o quizás sirva de apoyo para sus hijos, nietos, sobrinos o amigos que quieran aprender un poco más del maravilloso mundo de las matemáticas y la programación. No olviden dejar sus comentarios. Saludos y nos leemos pronto.
Si desean consultar un poco más del tema pueden usar las siguientes referencias:
- Demidovich, B. P., and I. A. Maron. Computational Mathematics Mir, Moscow, 1976.
- Björck, Åke. Numerical methods in matrix computations. Vol. 59. Cham: Springer, 2015.
- Burden, Richard L., and J. Douglas Faires. Numerical analysis. Ninth Edition. Cengage Learning. 2011.
Las imágenes, separadores y las ecuaciones fueron creadas y editadas por @abdulmath usando software libre, GNU Octave, , GIMP e Inkscape.
![](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmaTz241KQheijV5qjopsqAzLW8ZpBqCmKmjpaixyGV4Wz/STEMESPANOL06.png)
@SteemSTEM es un proyecto comunitario con el objetivo de promover y apoyar la Ciencia, la Tecnología, la Ingeniería y las Matemáticas en la blockchain Steem. @Stem-espanol es parte de esta comunidad, si desea apoyar el proyecto, puedes contribuir con contenido en español en las áreas de Ciencia, Tecnología, Ingeniería y Matemáticas, utilizando las etiquetas #steemstem y #stem-espanol.
![](https://steemitimages.com/640x0/https://cdn.steemitimages.com/DQmaSVacTwJ8v8uqJLMAKkNaLhisPAx7bi7KXVuzV1B3uQm/FinalGMath00.png)
Muchas gracias por el apoyo. Muy agradecido.
Saludos y Éxitos
Bueno debo decirte que aunque lo mio no son las matemáticas, debo felicitarte por tu gran trabajo mi amor, muy buena tu explicación. Besos y saludos.
Muchas gracias cariño, es bueno que tengas una buena apreciación de mis aportes. Y que hayas tenido el tiempo de leerla.
Saludos y un abrazo.
Woou sin duda es impresionante tu capacidad de comprender las matematicas. Gracias por darle valor a la comunidad hispana con tu contenido.
Hola @sandrag89, muchas gracias por tus comentarios. Trato cada día de dar un pequeño aporte a todos los que se interesen por tan bonitos temas (para mi gusto).
Saludos Cordiales.
Congratulations! This post has been upvoted from the communal account, @minnowsupport, by abdulmath from the Minnow Support Project. It's a witness project run by aggroed, ausbitbank, teamsteem, someguy123, neoxian, followbtcnews, and netuoso. The goal is to help Steemit grow by supporting Minnows. Please find us at the Peace, Abundance, and Liberty Network (PALnet) Discord Channel. It's a completely public and open space to all members of the Steemit community who voluntarily choose to be there.
If you would like to delegate to the Minnow Support Project you can do so by clicking on the following links: 50SP, 100SP, 250SP, 500SP, 1000SP, 5000SP.
Be sure to leave at least 50SP undelegated on your account.
This post has been voted on by the SteemSTEM curation team and voting trail in collaboration with @curie.
If you appreciate the work we are doing then consider voting both projects for witness by selecting stem.witness and curie!
For additional information please join us on the SteemSTEM discord and to get to know the rest of the community!