PNLS 82
PNLS 82
PNLS 82
3. INTRODUCCIÓN
Con esta investigación daremos a conocer varios de los aspectos que concierne la programación No
lineal y algunas de sus ramas. Empezando por definir su concepto, de donde proviene, cuál es su
importancia en ámbitos económicos, sociales y humanos. Además, se mencionan los problemas
clásicos que resuelve, que uso se le ha dado a lo largo del tiempo en las áreas de ingeniería,
administración entre otras, se enfatizará con ejercicios propuestos para ver cada uno de los métodos a
fondo para realizar este tipo de Programación y también las aplicaciones donde se abordará como la
programación No lineal fue utilizada con el fin de solucionar las problemáticas planteadas. Finalmente
mencionaremos los resultados obtenidos durante el desarrollo de la investigación.
3.1.1. Objetivos
3.1.1.1. Objetivo General
Mostrar como la programación no lineal separable, usando diferentes métodos como Base restringida,
Método de gradiente o método de secuenciales no restringidos además de las diferentes herramientas
digitales, sirve para la solución de problemas que involucran decisiones.
3.1.2. Justificación
Se busca dar a entender la importancia que tiene la programación no lineal separable en la ingeniería,
en el desarrollo económico y la administración abordando las distintas aplicaciones y ejercicios en las
cuales esta ha sido la principal herramienta de solución. Centraremos el documento en ejercicios típicos
y propuestos además de aplicaciones principales, los cuales explicaremos en detalle dando a entender
su contexto social real, justificación y método de solución a través de modelos de programación no
lineal separable.
La primera vez que se escuchó sobre Programación no lineal fue en el Año 1951 por el Matemático
Canadiense Albert W. Tucker y el Matemático Americano Harold W Kuhn los cuales venían
trabajando desde hace varios Años juntos con el campo de la teoría de juegos, en dicho año crearon
las condiciones de Kuhn – Tucker el cual es considerado como el primer método de optimización
para un problema No lineal. Este está basado en el proceso matemático conocido como multiplicador
de LaGrange y que solo incluía algunos problemas de programación lineal ya que para poder
3
resolverse por este método debía de tener unas condiciones básicas y así solucionarse. En el año 1951
estos dos matemáticos Sacan el libro “ Contributions to the Theory of Games ” donde explican y resaltan
los fundamentos de este método y más adelante en 1953 sacan el volumen 2 de su libro donde
explican casos especiales de su método expuesto, cabe resaltar que antes de la publicación del libro
ya se había difundido este método de programación no lineal a base de conferencias.
Sin embargo, investigadores encontraron años después que los fundamentos por los cuales se basaron
estos dos matemáticos ya había sido anunciados antes en 1939 por el matemático Americano William
Karush en su tesis de maestría llamada “ Minima of functions of several variables with inequalities as side
conditions” lo cual genero muchas controversias en la comunidad científica por parte de estos dos
Matemáticos, por lo que el método pionero fue re bautizado como Condiciones de Karush – Kuhn –
Tucker.
En el año 1956 sale la primera rama importante de programación no lineal denominada Programación
no lineal convexa siendo sus pioneros la Matemática francés Marguerite Frank y el Matemático
Americano Phill Wolfe con su método denominado Algoritmo de Frank – Wolfe o método del
gradiente que se publicó en un libro recopilatorio de investigaciones de la Universidad de Pricenton
llamado “Naval Research Logistics Quarterly Vol 3” y su Artículo tenia por nombre ”An algorithm for
quadratic programming” un dato curioso es que Frank y Wolfe formaron parte del proyecto logístico de
Princeton quienes dirigían los señores Kuhn y Tucker los padres de esta programación.
A partir de aquí muchos autores y matemáticos empezaron profundizar con el tema de programación no lineal
dejando de lado un poco a la programación lineal que era muy popular al comienzo de la década de los 50, en
1962 el investigador Americano William W. Cooper en colaboración de A charnes publica su libro “ Programming
with linear fractional functionals” El cual estaba enfocado a un nuevo tipo de Programación con fraccionarios, y
aunque su trabajo estuvo enfocado a programación lineal, se vieron aplicaciones bastante aceptables para
trabajar con la programación no lineal ya que casos exactos donde se requería transformar la función a lineal se
resolvían gracias a los métodos propuestos por Cooper.
y con estas bases nacen las demás ramas que constituyen esta programación y nacieron todo tipo de
investigaciones por parte de muchos autores e incluso de trabajos posteriores de los pioneros y primeros
matemáticos de este tipo de programación contribuyeron a la construcción de esta.
4
3.1.3.2 Desarrollo
La optimización No Lineal es una de las áreas más importantes de la matemática aplicada moderna,
con aplicaciones en campos de ingeniería y economía a finanzas, estadísticas, ciencias de gestión y
medicina. Existen varias posibilidades para la naturaleza del conjunto de restricciones, o función a
optimizar, dependiendo de estas características salen ramas principales donde muchos autores han
profundizado y realizado Algoritmos para resolver cada uno de estas.
La suposición adicional es que todas las funciones f(x) y g(x) son funciones separables. Una función
separable es una función en la que cada término incluye una sola variable, por lo que la función se
puede separar en una suma de funciones de variables individuales.
La Función Objetivo y las restricciones deben ser netamente separable y además La función Objetivo
puede tener exponentes En Caso de Maximización debe ser Cóncava, en caso de Minimización debe
ser Convexa
5
Las restricciones deben ser NO lineales y puede ser Convexa y La escancia de estos problemas es el de
eliminar las restricciones especiales, para convertir el problema lineal Su solución Óptima debe
satisfacer en teoría a la restricción especial.
3.1.3.3 Conceptualización
La programación no lineal en términos generales es el proceso de resolver un problema de
optimización definido por un sistema de igualdades y desigualdades, colectivamente denominadas
restricciones, sobre un conjunto de variables reales desconocidas, junto con una función objetivo a
maximizar o minimizar, donde algunas de las limitaciones o la función objetivo son de orden 2 en
adelante.
Programación no lineal sin restricciones Es el proceso de optimización de una función objetivo con
respecto a algunas variables en presencia de restricciones sobre esas variables. La función objetivo es
una función de coste o de energía que debe minimizarse, o una función de recompensa o función de
utilidad, que debe maximizarse. Las restricciones pueden ser restricciones duras que establezcan
condiciones para las variables que se requiera satisfacer o restricciones blandas que tienen algunos
valores variables que se penalizan en la función objetivo si, y en base a la extensión en que las
condiciones en las variables no son satisfechas.
6
Metodo Gradiente reduce el calculo de un minimo local a una secuencia de problemas de búsqueda
lineal modificando la búsqueda del gradiente evitando la entrada de la trayectoria de búsqueda a la
frontera de restricción
Metodo de Base Restringida consiste en identificar las funciones que dependen únicamente de una sola
variable, determinar las restricciones como una suma de funciones que dependen de una sola variable,
separar las restricciones que pueden ser lineales o separables, para obtener finalmente la solución
óptima con el método simplex normal.
Metodo de Charnes y Cooper consiste en realizar una transformada a partir del cambio de variable con
el propósito de convertir un problema no lineal en un problema lineal equivalente
Metodo Frank & Wolf está condicionado a problemas netamente convexos los cuales se solucionan a
travez de aproximaciones por medio de la matriz hessiana con el objetivo de Resolver el problema con
dicha aproximación, como un problema de PL con holguras complementarias
En 1994 Francisco Jubete sugiere el método Polar para resolución de problemas convexosde
una fomra más práctica y rápida basados en el teorema de farkas y Khun & Tucker
M. Bell & Yasuunori Iida en 1997 publicaron su libro “Transportation Network Analysis” hicieron un
amplio análisis a los métodos convexos y tomaron las aplicaciones de redes recopilándolas.
3.1.3.7 Mentefacto
MATLAB – Mathworks
Excel Solver – Microsoft
CPlex- IBM
AMPL Solver – AMPL
Gams – Gams
9
CONOPT – MPL
CVXOP – Phyton
LINGO – Lindo Systems
OPT ++ - Sandia Software
GLPK – GNU proyect
Knitro – Artelys
FSQP – pyopt
JOptimizer – Apache Maven
AMMS Optimization – AMMS
WinGULF – Martos B.
ENUNCIADO:
SOLUCION:
Pasos Procedimiento
Se separan las funciones,
tanto en la F.O. como en las
restricciones.
10
X1=5;
11
X2=5;
Z=200:
Tabla 3.1 Problema Clásico Resuelto 1
Ejercicio No. 2:
ENUNCIADO:
SOLUCIÓN:
Modelo linealizado
3. Solucionar el modelo linealizado usando el método de las 2 fases pero con el criterio de base
restringida:
Primera fase
13
Segunda fase
En
términos de las variables originales:
x1 = 1
x2 = 0.1081
FUENTE:
https://www.me.utexas.edu/~jensen/ORMM/supplements/methods/nlpmethod/S1_separable.pdf
Ejercicio No. 3:
ENUNCIADO:
14
SOLUCIÓN:
El modelo no es separable por la suma al cuadrado de la función objetivo. Para cambiar esta situación
se hace la siguiente sustitución:
3. Solucionar el modelo linealizado. Como las funciones separables son cóncavas se puede omitir la
restricción de base restringida, luego se puede aplicar el simplex sin lo que conlleva la anterior
condición.
x1 = 2, x2 = 3 x3 = 5
FUENTE:
http://web.mit.edu/15.053/www/AMP-Chapter-13.pdf
Ejercicio No. 4:
ENUNCIADO
16
Paso 0: Separar la Función Objetivo y las restricciones gij ( x j ), donde i corresponde al número de la
restricción y j a la variable asociada.
Función Objetivo:
f 1 ( x 1 )=5 x 1−x 12
f 2 ( x 2 )=3 x 2−x 22
Restricciones:
g12 ( x 2 ) =x2
Paso 1:
Paso 2: Reemplazo de variables no lineales por suma de variables lineales (t), para ello se genera un
par de tablas para cada función separada:
Paso 3: Plantear el nuevo modelo lineal, en el cuál se unen funciones objetivo y restricciones
planteadas para cada función:
Sujeto a:
t ik , S 1 , S2 ≥0
Cj 4 6 2 2 0 -4 0 0 0 0
V.B. t 12 t 13 t 22 t 23 t 24 t 25 S1 S2 t 11 t 21 bj
S1 0 2 32 1 2 3 4 1 0 0 0 32
S2 0 1 2 2 8 18 32 0 1 0 0 32
t 11 0 1 1 0 0 0 0 0 0 1 0 1
18
t 21 0 0 0 1 1 1 1 0 0 0 1 1
Zj 0 0 0 0 0 0 0 0 0 0 0
C j−Z j 4 6 2 2 0 -4 0 0 0 0
El candidato a entrar (t 13), puede hacerlo debido a que el criterio de adyacencia se mira con respecto a
las otras variables lineales t ik ,
Cj 4 6 2 2 0 -4 0 0 0 0
V.B. t 12 t 13 t 22 t 23 t 24 t 25 S1 S2 t 11 t 21 bj
S1 0 -30 0 1 2 3 4 1 0 -32 0 0
S2 0 -1 0 2 8 18 32 0 1 -2 0 30
t 13 6 1 1 0 0 0 0 0 0 1 0 1
t 21 0 0 0 1 1 1 1 0 0 0 1 1
Zj 6 6 0 0 0 0 0 0 6 0 6
C j−Z j -2 0 2 2 0 -4 0 0 0 0
Cj 4 6 2 2 0 -4 0 0 0 0
V.B. t 12 t 13 t 22 t 23 t 24 t 25 S1 S2 t 11 t 21 bj
t 22 2 -30 0 1 2 3 4 1 0 -32 0 0
S2 0 59 0 0 4 12 24 -2 1 62 0 30
t 13 6 1 1 0 0 0 0 0 0 1 0 1
t 21 0 30 0 0 -1 -2 -3 -1 0 32 1 1
Zj -54 6 2 4 6 8 2 0 -58 0 6
C j−Z j 58 0 0 -2 -6 -12 -2 0 58 0
Cj 4 6 2 2 0 -4 0 0 0 0
V.B. t 12 t 13 t 22 t 23 t 24 t 25 S1 S2 t 11 t 21 bj
t 22 2 0 0 1 1 1 1 0 0 0 1 1
19
1∗1 2∗29 59
x 1= + = 1,96667
30 30 30
x 2=1∗1=1
5∗59 59 2
Z=
30
−
30 ( )
+ 3∗1−12=7,96556
REFERENCIAS
20
Ejercicio No. 5:
PROBLEMA:
MaxZ=3 x 13 +3 x2
S . A . :3 x 12 +6 x 2 ≤ 30
x1 , x2 ≥ 0
Pasos Procedimiento
Separamos tanto la f 1 ( x 1 )=3 x 13
función objetivo como
f 2 ( x 2 )=3 x 2
las restricciones
g11 ( x 1 )=3 x 12
g12 ( x 2 ) =6 x2
x 2=√ 10=3.1622
Encontramos f ( t ) yg ( t ) , |
teniendo en cuenta que
k t1k f 1(t 1 k ) g11 ( x 1 )
t 1 k va hasta el numero
1 0 0 0
siguiente de la
2 1 3 3
aproximación de la
3 2 24 12
variable
4 3 81 27
21
5 4 192 48
f 1 ( x 1 )=0t 11 +3 t 12+24 t 13 +81 t 14 +192 t 15
t 11 ,t 12 , t 13 , t 14 , t 15 , x2 ≥ 0
t 11 ,t 12 , t 13 , t 14 , t 15 , x2 ≥ 0
Obtenemos el primer +
tablero del Simplex
cj 3 3 24 81 0 0 ---
92
VB x2 t 12 t 13 t 14 t 15 s1 t 11 b
s1 ∨0 6 3 12 27 48 1 0 30
t 11 ∨ 0 0 1 1 1 1 0 1 1
zj 0 0 0 0 0 0 0 0
c j−z j 3 3 24 81 192 0 0 ----
1/48
zj 24 12 48 108 192 4 0 120
c j−z j -21 -9 0 -27 0 0 0 ----
Teniendo como
resultado:
3 5
t 11= ,t 15=
8 8
Reemplazamos en los
valores obtenidos en
BIBLIOGRAFIA
REFERENCIA ELECTRONICA
http://www.virtual.unal.edu.co/cursos/sedes/manizales/4060015/Lecciones/Capitulo
%20V/pseparable.htm
Ejercicio No. 6:
ENUNCIADO:
MinZ=X 12−4 X 1 −X 2
Sujetoa:2 x 12+8 x 22 ≤ 30
x1 , x2 ≥ 0
23
Pasos
Procedimiento
Paso 1
Paso 2
0 x 1+ 8 x 22=30 ∴ x2=1.9364
Paso 3
24
1 0 0 0 0 0 0
2 1 1 -3 -1 2 8
3 2 2 -4 -2 8 32
Paso 4
S.a.
t 11 +t 12+t 13=1
Paso 5
Cj 0 3 4 0 1 2 0 bj
V.B. t 11 t 12 t 13 t 21 t 22 t 23 S1
S1 0 0 2 8 0 8 32 1 30
t 11 0 1 1 1 0 0 0 0 1
t 21 0 0 0 0 1 1 1 0 1
Zj 0 0 0 0 0 0 0
Z j−C j 0 -3 -4 0 -1 -2 0
25
Paso 6
Se desarrolla el método Simplex, haciendo cero los valores por encima y por debajo del pivote,
finalizando por la escogencia del nuevo pivote.
Cj 0 3 4 0 1 2 0 bj
V.B. t 11 t 12 t 13 t 21 t 22 t 23 S1
S1 0 -8 -6 0 0 8 32 1 22
t 13 4 1 1 1 0 0 0 0 1
t 21 0 0 0 0 1 1 1 0 1
Zj 4 4 4 0 0 4 0
Z j−C j 4 1 0 0 -1 2 0
Paso 7
Se desarrollará el método simplex hasta finalizar al tener todos los valores C i−Z i mayores a cero,
siempre teniendo en cuenta las restricciones de adyacencia para la escogencia del pivote.
Cj 0 3 4 0 1 2 0 bj
V.B. t 11 t 12 t 13 t 21 t 22 t 23 S1
t 23 2 -1/4 -3/16 0 -1/4 0 1 1/32 7/16
t 13 4 1 1 1 0 0 0 0 1
t 22 1 1/4 3/16 0 5/4 1 0 -1/32 9/16
Zj 15/4 61/16 4 3/4 1 2 1/32
Z j−C j 15/4 13/16 0 3/4 0 0 1/32
Tabla 3. 6 Problema Clásico Resuelto 6
26
Solución al problema:
t 13=1
t 23=7/16
t 22=9/16
Paso 8
Se reemplazan los valores obtenidos en las columnas t 1 k y t 2 k para hallar los valores de x 1 y x 2
respectivamente, luego se reemplazan estos valores en la función objetivo para hallar la solución
REFERENCIAS
https://es.scribd.com/doc/59899398/Introduccion-a-la-Programacion-No-Lineal
27
Ejercicio No. 7:
Maximizar:
2 x1 x 1+3 x 2 x 2 ≤6
x1 , x2 ≥ 0
f2(X2) = 0t21+1t22+2t23
g11(X1) = 0t11+2t12+8t13
g12(X2) = 0t21+3t22+12t23
28
2 t 12 +8 t 13 +3 t 22+12 t 23 + s1=6 ,
Proponer el nuevo problema
t 11 +t 12+t 13=1
equivalente.
t 21+t 22+ t 23=1
t ij ≥ 0.i , j=1,2,3
Con la solución:
Con la solución:
Cj 1 0 1 2 0 0 0
V.B. t12 t13 t22 t23 s1 t11 t21 bj
s1 0 0 6 0 9 1 -2 -3 1
29
t12 1 1 1 0 0 0 1 0 1
t22 1 0 0 1 1 0 0 1 1
Zj 1 1 1 1 0 1 1 2
Cj-Zj 0 -1 0 1 0 -1 -1
Con la solución:
Referencias
Ejercicio No. 8:
Enunciado:
MaxZ=3 x 1 + x 2+ 2 x 12
s. a: 3 x 12 +5 x 2 ≤ 24
x1 , x2 ≥ 0
Solución:
Pasos Procedimiento
F.O:
f 1 ( x 1 )=3 x 1 +2 x 12
Paso Inicial. Antes de
empezar, se debe separar la f 2 ( x 2 )=x 2
función objetivo y las
Restricciones:
restricciones.
Paso 2. Encontrar f ( t ) y g ( t )
.
k t 1k f 1(t 1 k ) g11 ( t 1 k )
1 0 0 0
2 1 5 3
3 2 14 12
31
4 3 27 27
Cj 1 5 14 27 0 0
V.B. x2 t 12 t 13 t 14 s1 t 11 bj
s1 0 5 3 12 27 1 0 24
t 11 0 0 1 1 1 0 1 1
Zj 0 0 0 0 0 0 0
Paso 4. Solucionar la nueva
C j−Z j 1 5 14 27 0 0
función por el Método
Simplex Modificado,
teniendo en cuenta el En la primera iteración el candidato a entrar es t 14 y a salir es s1 pero
criterio de adyacencia para no puede hacerlo ya que no es adyacente a t 11, así que la variable que
la entrada de una variable a entra es t 13 y sale t 11.
la tabla.
Cj 1 5 14 27 0 0
V.B. x2 t 12 t 13 t 14 s1 t 11 bj
s1 0 5 -9 0 15 1 -12 1
t 13 14 0 1 1 1 0 1 1
Zj 0 14 14 14 0 14 14
32
C j−Z j 1 -9 0 13 0 -14
Cj 1 5 14 27 0 0
V.B. x2 t 12 t 13 t 14 s1 t 11 bj
2
t 14 1/3 −3/5 0 1 1/15 −4/5 4 /5
7
1
t 13 −1/ 3 8 /5 1 0 −1/15 9 /5 1/5
4
Zj 13/3 31/5 14 27 13/15 18/5 122/5
C j−Z j −10/3 −6 /5 0 0 0 −18/5
122 1 4
Z= ; x 2=0 ; t 11 =0 ; t 12=0 ; t 13= ; t 14=
5 5 5
Paso 5. Deducir la solución Según lo que el resultado del método simplex, tenemos que:
del problema original.
x 1=0 t 11 +1t 12+ 2t 13+ 3t 14
x 1=2 ( 15 )+3( 45 )
x 1=2.8
Z=3 x1 + x 2 +2 x 12
2
Z=3 ( 2.8 ) +0+2 ( 2.8 )
33
Z=24.08
122
Z= ; x 1=2.8; x 2=0
5
REFERENCIAS ELECTRÓNICAS
Ejercicio No. 9:
ENUNCIADO:
s.a.
PROCEDIMIENTO:
2. Los intervalos de las variables de decisión estimados con las restricciones son [0, 3] para x1 y [0, 4]
para x2. De esta manera:
Para i = 1:
Para i = 2:
4. La solución es:
REFERENCIAS
ENUNCIADO:
1. Maximizarz=x 1−3 x 32
Sujetoa:
2 x1 +2 x 22 ≤ 32
x1 , x2 ≥ 0
PROCEDIMIENTO
Tomando en cuenta la restricción se halla una solución para x 2, que según el modelo
(derecha) es la variable no lineal, dando como resultado x 2 = 4
f ( x 1 ) =x1
k t2k f2(t2k) g12(t2k)
f ( x 2 ) =−3 x32 1 0 0 0
2 1 -3 2
g11 =2 x 1
3 2 -24 8
g12=2 x 22 4 3 -81 18
5 4 -192 32
32
x 22=
2
x 2=√ 16=4
Se resuelve el problema por simplex (se divide la fila 1 en 2) teniendo en cuenta el criterio
de adyacencia y se obtiene como tabla final (derecha)
X1 = 16
X2 = 0
Z = 16
REFERENCIAS
Una fábrica de textiles desea aumentar sus ganancias modificando el número de camisetas y
pantalones producidos por una maquina al día, la ganancia que se obtiene son 2 por camiseta
y 1 por pantalón, pero por cada venta se pierde el valor del cuadrado del número de las
camisetas vendidas, para que la producción pueda ser fabricada no puede sobrepasarse que
2(camisetas producidas)2 +3(pantalones producidos)2 sea mayor a 6, pues este es el valor que
regula la producción de la máquina.
Maximizar:
2 x1 x 1+3 x 2 x 2 ≤6
x1 , x2 ≥ 0
39
f1(X1) = 0t11+1t12+0t13
f2(X2) = 0t21+1t22+2t23
g11(X1) = 0t11+2t12+8t13
g12(X2) = 0t21+3t22+12t23
Max . z=t 12+ t 22 +2 t 23 , sujetoa :
2 t 12 +8 t 13 +3 t 22+12 t 23 + s1=6 ,
Proponer el nuevo problema
t 11 +t 12+t 13=1
equivalente.
t 21+t 22+ t 23=1
t ij ≥ 0.i , j=1,2,3
40
Cj 1 0 1 2 0 0 0
s1 0 2 8 3 12 1 0 0 6
t11 0 1 1 0 0 0 1 0 1
t21 0 0 0 1 1 0 0 1 1
Zj 0 0 0 0 0 0 0 0
Cj-Z 1 0 1 2 0 0 0
Con la solución:
Cj 1 0 1 2 0 0 0
Resolver mediante el método V.B. t12 t13 t22 t23 s1 t11 t21 bj
simplex s1 0 0 6 3 12 1 -2 0 4
t12 1 1 1 0 0 0 1 0 1
t21 0 0 0 1 1 0 0 1 1
Zj 1 1 0 0 0 1 0 1
Cj-Zj 0 -1 1 2 0 -1 0
Con la solución:
Cj 1 0 1 2 0 0 0
s1 0 0 6 0 9 1 -2 -3 1
41
t12 1 1 1 0 0 0 1 0 1
t22 1 0 0 1 1 0 0 1 1
Zj 1 1 1 1 0 1 1 2
Cj-Zj 0 -1 0 1 0 -1 -1
Con la solución:
Respuesta:
Se deben fabricar 1 pantalón y 1 camiseta por maquina para obtener una ganancia de 2.