1

Taller
Ejercicios prácticos y evaluación
Métodos de Optimización
Cesar A. Avendaño T., Miembro, IEEE

x1 , x2 , x3 , x4 ≥ 0
Abstract—.

Index Terms—Permitividad, permeabilidad, vacío.

Luego tenemos que:
I. INTRODUCCIÓN

E 2 00 0

[ ]
el presente documento se presenta la solución a los
N
problemas de programación lineal, programación H= 0 80 0
cuadrática y programación sin restricciones, a través del uso 0 00 0
de los métodos simplex, Newton, Levenberg-Marquardt, 0 00 0
BFGS, DFP y Nelder-Mead.

II.PROGRAMACIÓN CUADRÁTICA c= [−8 −16 0 0 ]

Resuelva el siguiente problema de optimización:
A. Resuelva el siguiente problema de optimización: A= [ 11 1 10
0 01 ]
min f ( x )=−8 x 1−16 x 2 + x 21+ 4 x 22 5
Sujeto a b=[ ]
3
( x 1 + x2 ) ≤ 5 , x 1 ≤3 , x 1 ≥ 0 , x 2 ≥0 Ahora el problema se resume con las condiciones de Kuhn-
Tucker como sigue a continuación:
Ahora, para poder resolver el problema de programación
cuadrática tipo dos tenemos que: Ax=b
T
1 H x + A λ−μ=−c
f ( x )= x T H x +c T x x ,μ≥0
2
x T μ=0
Ax=b x ≥ 0
Ahora, podemos obtener una minimización del problema si
2 2 introducimos dos variables de holgura o SLACK de la
f ( x )=x +4 x +(−8 x 1−16 x2 )
1 2 siguiente manera:
Donde,
min
u +u
1 T x , μ , λ , u1 ,u 2 1 2
x H x =x 21 +4 x 22
2
Y b≥0
T c ≤0
c x=−8 x 1−16 x 2
Ax+u1 =b
Ahora,

x 1+ x2 + x 3=5 H x + A T λ−μ+u 2=−c

x 1+ x 4=3
 2

x , μ , u1 ,u 2 ≥ 0 1 1 1 0 0 0 0 0 0 0

x T μ=0
Lo que es igual a:

A 0 0 I0
1
N 0= 2
0
0
0
[ 0
0
8
0
0
0
0
0
0
0
1
0
0
0
0
0
1
1
1
0
0 0 0
1 −1 0
0 0 −1 0
0 0
1 0
0
0

0 −1 0
0
0
0
0

0 −1
]
A0 =
[ ] 1 0 0 0 0 0

[ ]
H A T −I 0 I 0 1 0 0 0 0
b 0= 0 0 1 0 0 0
b 0= [−cb ] 0
0
0
0
0
0
1
0
0
1
0
0
0 0 0 0 0 1
0
0
c0 0
1
1
[] x B 0 =B−1
0 b 0 =[ 5 3 816 0 0 ]

c B 0=[ 11 11 11 ]
'
'

c N 0 = [ x 1 x 2 x3 x 4 λ 1 λ 2 μ 1 μ 2 μ 3 μ 4 ]
'

x
λ
x0 μ
u1
u2
[]
Ahora el problema se puede escribir como:
c N 0= [ 0 0 0 0 0 0 0 0 0 0 ]

Z 0=c TB 0 x B 0=32 ; x T μ=0

'

PT =c TN 0−cTB 0 B−1
0 N0
min c 0 x0 sujeto a : A 0 x0 =b0 y x 0 ≥ 0

De lo cual se obtiene: PT =[ −4−9−1−1−3−2 1 11 1 ]

A0 =¿ '
y=B−1
0 N .2 =[ 1 0 0 8 0 0 ]

1 1 1 0 0 0 0 0 0 0 1 0 0 0Iteración
0 01

[ 1
2
0
0
0
0
0
8
0
0
0
0
0
0
0
1
0
0
0
0
0
1
1
1
0
0 0

0 0
1 0
0
1 −1 0
0 0 −1 0
0
0

0 −1 0
0
0
0
0

0 −1
0
0
0
0
0
1
0
0
0
0
0
1
0
0
0
0 0 0
0min

0 1 0
x
( )
0 0B =2
1 0 0y

0B (04 )=N

1
1 (2)
] 1 1 0 0 0 0 0 0 0
1

[ 0 0 1 0 0 0 0 0 0

]
'
b 0=[ 53 8 16 0 0 ]
N 1= 2 0 0 0 1 1 −1 0 0 0
c 0= [ 0 0 0 0 0 0 0 0 0 0 1 11 11 1 ] 0 1 0 0 1 0 0 −1 0 0
0 0 0 0 1 0 0 0 −1 0
0 0 0 0 0 1 0 0 0 −1
 3

'
1 0 0 1 0 0 c B 2=[ 01 1 0 11 ]
0
b 1= 0
0
0
0
[ 1
0
0
0
0

−1
0
1
0
0
0

x B 1=B b 0=[ 33 8 2 0 0 ]
0
0
8
0
0
0
0
0
1
0
0
0
0
0
1

'
] c N 2=[ μ11 μ 22 x 3 x 4 λ1 λ2 μ 1 μ2 μ3 μ 4 ]

c N 2=[ 11 0 0 0 0 0 0 0 0 ]

Z 0=c TB 2 x B 2=2 ; x T μ=0

'
'

1
PT =c TN 2−cTB 2 B−1
2 N2
'
c B 1 =[ 1 1 10 1 1 ]
PT =[ 4 0.625 3−1−2.375−21 0. 3751 1 ]
'
c B 1 =[ x1 μ22 x 3 x 4 λ 1 λ 2 μ1 μ 2 μ3 μ4 ] '
y= [−0.1250.125 1.25 0.125 1 ]
'
c B 1 =[ 01 0 0 0 0 0 0 0 0 ]
Iteración 3

Z 0=c TB 1 x B 1=14 ; x T μ=0 xB

PT =c TN 1−c TB 1 B−1
1 N1
min ( )
y
=0

B ( 2 )=N (5)
PT =[ −4 1.125−1−1−1.875−2 1−0. 125 1 1 ]
' 1 0 1 0 0 0 0 0 0 0

[ ]
y=B−1
1 N 1 .1= [ 1 12 0 0 0 ] 0 0 0 1 1 0 0 0 0 0
N 3= 0
Iteración 2
0 0 0 0 1 −1 0 0 0
0 1 0 0 0 0 0 −1 0 0
xB 0 0 0 0 0 0 0 0 −1 0
min ( ) y
=1 0 0 0 0 0 1 0 0 0 −1

B (1 )=N (1) 1 0 0 1 0 0

1
1
N 2= 0
0
0
0
[
0
0
0
1
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
1
1
0
0 0
0 0
0
0
1 −1 0
0 0 −1 0
0 0
1 0
0
0
0

0 −1 0
0
0
0
0
0

0 −1
] [ 1
b 3= 2
0
0
0

x B 3=B−1
1
0
0
0
0
0
1
0
0
0
0
0
8
0
0
0
0
0
1
0
0
0
0
0
1

'
]
3 b 0=[ 3 0 22 0 0 ]

'
1 0 0 1 0 0 c B 3 =[ 0 0 10 1 1 ]
1
b 2= 2
0
0
0
[
x B 2=B−1
1
0
0
0
0
0
1
0
0
0
0
0
8
0
0
0
0
0
1
0
0
0
0
0
1

'
] c N 3=[ μ11 μ 22 x 3 x 4 μ12 λ 2 μ1 μ2 μ3 μ4 ]

c N 3=[ 11 0 0 1 0 0 0 0 0 ]

Z 0=c TB 3 x B 3=2 ; x T μ=0

'
'

2 b 0= [ 3 02 2 0 0 ]
PT =c TN 3−cTB 3 B−1
3 N3
 4

B ( 6 )=N (6)
PT =[ −15 3−16 18 19−2 1−2 11 ]
1 0 0 0 0 0 0 0 0 0

[ ]
'
y= [ 0−8 8 1 8 0 ] 0 0 0 1 1 0 0 0 0 0
N 5= 0 0 0 0 0 0 −1 0 0 0
Iteración 4
0 1 0 0 0 0 0 −1 0 0
0 0 1 0 0 0 0 0 −1 0
xB
min ( )y
=0 0 0 0 0 0 1 0 0 0 −1

B (5 )=N (3) 1 0 0 1 1 0

1
0
N4= 0
0
0 [
0
0
0
1
0
0
0
0
0
1
0
1
0
0
0
0
1
0
0
0
0 0
0 0
0
0
1 −1 0
0 0 −1 0
0 0
0
0
0

0 −1 0
0
0
0
0
] [ 1
b 5= 2
0
0
0

x B 5=B−1
0
1
1
1
0
0
1
0
0
0
0
0
8
0
0
0
0
0
0
0
0
1
0
0
1

'
]
0 0 0 0 0 1 0 0 0 −1 5 b 0 =[ 3 0 22 0 0 ]

'
1 0 0 1 1 0 c B 5 =[ 0 0 10 0 0 ]
1
b 4= 2
0
0
0
[ ]
x B 4 =B−1
0
1
1
1
0
0
1
0
0
0
0
0
8
0
0
0
0
0
0
0
0
0
0
0
1

'
c N 5=[ μ11 μ 22 μ23 x 4 μ12 μ 24 μ 1 μ2 μ 3 μ 4 ]

c N 5=[ 11 1 0 11 0 0 0 0 ]

Z 0=c TB 5 x B 5=2 ; x T μ=0

'
'

4 b0= [ 3 0 2 20 0 ]
PT =c TN 5−cTB 5 B−1
5 N5
'
c B 4= [ 0 0 1 0 0 1 ]
PT =[ 11 22 3 21 0−1−1 ]
'
c N 4 =[ μ 11 μ22 μ 23 x 4 μ 12 λ 2 μ1 μ 2 μ3 μ 4 ] '
y= [ 0 0 1 0 0−1 ]
'
c N 4 =[ 1 1 10 1 0 0 0 0 0 ] Iteración 6

Z 0=c TB 4 x B 4=2 ; x T μ=0 y >0

PT =c TN 4 −c TB 4 B−1 B (3 )=N (10)

4 N4

PT =[ 11 22 3−21 0−11 ] 1 0 0 0 0 0 0 0 0 0

y= [ 0 0 1 0 0 1 ]

Iteración 5

min ( xy )=0
B
'
0
N 6= 0
0
0
0
[ 0
0
1
0
0
0
0
0
1
0
1
0
0
0
0
1
0
0
0
0
0 0 0
1 −1 0
0 0
0 1
0 0 −1 0 0
0 0
1 0
0 −1 0
0 0 0
]
 5

1 0 0 1 1 0

[
1
b 6= 2
0
0
0

x B 6 =B−1
0 0
1 0
1 0
1 0
0 −1
0
0
8
0
0
0
0
0
0
0

'
0
0
0
0
1
]
6 b 0=[ 3 0 22 0 0 ]

'
c B 6=[ 0 0 0 0 0 0 ]
'
c N 6=[ μ11 μ 22 μ23 x 4 μ12 μ 24 μ 1 μ2 μ3 μ 21 ]
Fig. 1. Tensión usada en la generación, transmisión y distribución de energía
' eléctrica en Canadá. Fuente: Table 25-10-0022-01 Installed plants, annual
c N 6=[ 1 11 0 11 0 0 0 1 ] generating capacity by type of electricity generation, consulta realizada el 3 de
octubre de 2020.

Z 0=c TB 6 x B 6=2; x T μ=0

III. CONCLUSIONES
PT =c TN 6−cTB 6 B−1
6 N6
Se puede concluir, que el uso de energía eléctrica en EEUU
supera enormemente al uso de energía en Canadá,
PT =[ 11 10 1 1 0 0 01 ] probablemente al territorio electrificado y habitado, puesto
que Canadá teniendo una extensión de terreno de 9985
PT =≥ 0 millones de km², mayor a la de EEUU, tiene aún mucho
territorio que no es habitado.
Luego los puntos óptimos serán iguales a: De igual forma se puede decir que el uso de tecnologías
renovables esta apenas en sus comienzos y el desarrollo de
3 nuevas tecnologías es necesario para poder potencializar las

0
0
[] []
x= 2 λ= 0
2
infraestructuras existentes y aumentar la capacidad de
producción de este tipo de energías, amigables con el medio
ambiente.

REFERENCIAS
0 0

0
2
[] [] []
μ= 0 μ1= 0 μ2=
2
0
0
0
Banco Mundial. (2020). Banco Mundial.

Medio
TABLA I
CONSUMO ELÉCTRICO CANADÁ VS EEUU

EEUU (Ktoe) Canadá

Transporte 638100 68195

Industria 276582 45839
Residencial 270940 35937
Comercial y publica 217594 26598
Energía no utilizada 149417 22178
Para el radio de los conductores se toma como referencia los Agricultura 21153 7253
parámetros que se observan en la siguiente tabla 1. No especificada 20343 63

Fuente: https://www.iea.org/data-and-statistics?country=USA&fuel=Energy
%20consumption&indicator=TotElecCons , consulta realizada el 3 de
Octubre de 2020.

https://datos.bancomundial.org/pais/canada

CCA. (1999). Estudio temático 1. El maíz en México: Algunas

 6

implicaciones ambientales del Tratado de Libre

Comercio de América del Norte Índice. 118.

CER. (2020). Statistics Canada.

https://doi.org/https://doi.org/10.25318/2510002201-eng

IEA. (2020). Data and statistics. https://www.iea.org/data-

and-statistics?country=CANADA&fuel=Energy
supply&indicator=TPESbySource

Cesar A. Avendaño T. Es Tecnólogo en

mecatrónica de la Institución Universitaria
Pascual Bravo, Tecnólogo en
electromecánica de la Institución
Universitaria Pascual Bravo, Ingeniero en
Electromecánica del Instituto Tecnológico
Metropolitano, Magister en Automatización
y Control Industrial del Instituto
Tecnológico Metropolitano, es investigador
del laboratorio de sistemas de control y
robótica del Instituto Tecnológico Metropolitano y es
estudiante de Doctorado en Ingeniería Mecánica y
Mecatrónica de la Universidad Nacional de Colombia, sede
Medellín.
Fue docente de las áreas de Robótica y Electrónica digital en
la misma institución, ha trabajado desde hace 14 años con
Sistemas Automáticos de Control y Robótica, en el Metro de
Medellín y entre sus publicaciones ha aplicado sistemas
embebidos para la adquisición de datos y la caracterización de
señales.