9. Anexos.

steinb6.txt
6 12 1
9 36 1
10 17 1
12 5 1
16 8 1
17 36 2
18 22 1
19 28 1
21 5 1
Anexos.
Como se puede observar, el nombre de la función es el mismo que el del fichero que la
alberga. Además otro detalle importante es que la función objetivo aparece con un signo -
, lo que indica que lo que se quiere es maximizar el resultado, no minimizarlo como
normalmente hace fmincon.
También aparecen algunas líneas comentadas, eso es debido a que se utiliza o la primera
o la segunda función de integridad, según se comenten unas líneas u otras.
Si el número de argumentos de salida que se le pide a la función es mayor que uno,
entonces el gradiente de la función también se calcula. Este gradiente es necesario para
conseguir mejores resultados puesto que si no se proporciona analíticamente como lo
hacemos nosotros, Matlab lo calcula numéricamente y los resultados son peores.
A la función objetivo se le pasan como parámetros de entrada:
• x – Vector de entrada.
• coste – Coste de los arcos.
• NumTerminales – Necesario para la segunda función de integridad.
• x0 – solución inicial. La solución inicial es necesaria para la correcta evolución de la
optimización. En el problema, la solución inicial se ha hallado para cada problema
especificado gracias a la implementación del algoritmo de Kruskal en Visual Basic. Se
ha aprovechado que para la primera parte del Proyecto ya se había obtenido esta primera
solución. Para ello, se ha exportado en cada ejecución del programa en Visual Basic cada
una de estas primeras soluciones a un fichero en texto plano que puede ser leído por
Matlab y por ello utilizado por fmincon.
Es necesario hacer hincapié en que la documentación de Matlab afirma que si esta
primera solución es no admisible, el programa va a buscar una primera solución
admisible que pueda ser utilizada para el arranque del algoritmo. Sin embargo debido a la
complejidad del problema, sobre todo de las restricciones, Matlab no es capaz de
encontrar una solución admisible si la primera solución proporcionada no lo es.
A continuación se muestra el fichero PrimeraSol.txt generado para el problema
steinb6.txt:
3