OPTIMIZACIÓN ESTÁTICA
OPTIMIZACIÓN CON RESTRICCIÓN
1. Hay dos industrias, una incubadora y otra dedicada a la puesta de huevos, también hay dos
bienes, huevos y pollos. El precio de un huevo es 1, el precio del pollo es p. La industria
incubadora produce 𝑦1 pollos usando 𝑥1 huevos, los cuales son comprados a la industria
1/2
ponedora. Su función de producción es 𝑦1 = 4𝑥1 . La industria ponedora produce 𝑥2
huevos usando 𝑦2 pollos, los cuales son comprados a la empresa incubadora. Su función de
1/2
producción es 𝑥2 = 𝑦2 . (No tomaremos en cuenta cualquier consideración de los
retrasos, restricciones de no negatividad, la indivisibilidad de los pollos vivos o huevos, y
afirma que prevalece por todas las partes). Expresar en términos de 𝑝 el beneficio;
maximizándolo en función de las cantidades de pollos y huevos producidos o usados en
cada industria. Asumiremos que todos los huevos son producidos en la industria ponedora
y son usados como insumos por la empresa incubadora. ¿Cuál es el valor de 𝑝? ¿Cuál es la
cantidad de pollos dispuestos para el consumo?
Solución
MÉTODO DE SUSTITUCIÓN
𝑦1 = 𝑦2 + 𝑦3
INCUBADORA PONEDORA
INSUMO 𝑥1 𝑦2
PRODUCTO 𝑦1 𝑥2
FUNCIÓN 1/2 1/2
𝑦1 = 4𝑥1 𝑥2 = 𝑦2
1/2
Incubadora 𝑀𝑎𝑥 𝐵1 = 𝑝𝑦1 − 1𝑥1 → 𝑀𝑎𝑥 𝐵1 = 4𝑝𝑥1 − 1𝑥1
1/2
Incubadora 𝑀𝑎𝑥 𝐵2 = 1𝑥2 − 𝑝𝑦2 → 𝑀𝑎𝑥 𝐵2 = 𝑦2 − 𝑝𝑦2
Condiciones de primer orden (incubadora) 𝑝 > 0
𝜕𝐵1 2𝑝 1
= − 1 = 0 → 𝑥1 = 4𝑝 2 ∧ 𝑦1 = 4(4𝑝 2 )2 = 8𝑝
𝜕𝑥1 √𝑥1
Condiciones de segundo orden
𝜕𝐵1 𝑝 𝑝 1
=− =− =− 2<0
𝜕𝑥1 √43 𝑝 6 8𝑝
√𝑥13
𝑝 > 0 ∧ 𝑥1 > 0
Condiciones de primer orden (incubadora) 𝑝 > 0
𝜕𝐵2 1 1 1
= − 𝑝 = 0 → 𝑦2 = 2 ∧ 𝑥2 =
𝜕𝑦2 2√𝑦2 4𝑝 2𝑝
Condiciones de segundo orden
𝜕𝐵2 1 1
=− =− < 0 → 𝑠𝑖𝑒𝑚𝑝𝑟𝑒 𝑛𝑒𝑔𝑎𝑡𝑖𝑣𝑜
𝜕𝑦2 1
4 √𝑦23 4√ 3 6
4 𝑝
Por tanto, la función es cóncava, y los valores se optimizan para los valores de 𝑥1 , 𝑥2 . Ahora
calculamos el valor de 𝑝
𝑥1 = 𝑥2
1
= 4𝑝 2
2𝑝
8𝑝 3 − 1 = 0
(2𝑝 − 1)(4𝑝 2 + 4𝑝 + 1) = 0 → 𝑝 = 1/2
La cantidad de pollos dispuesta para el consumo
𝑦3 = 𝑦1 − 𝑦2
1
𝑦3 = 8𝑝 −
4𝑝 2
1
𝑦3 = 8 ( ) − 1 = 3
2
Multiplicadores de Lagrange
Si se tiene una función escalar 𝑓: ℝ𝑛 → ℝ restringida por un conjunto de 𝑚 ecuaciones:
𝑔1 (𝑥⃗) = 𝑏1
𝑔 (𝑥⃗) = 𝑏2
𝑔(𝑥⃗) = { 2
⋮
𝑔𝑚 (𝑥⃗) = 𝑏𝑚
Siendo 𝑛 > 𝑚, formamos la función Lagrangiana
𝑚
𝐿(𝑥⃗, 𝜆⃗) = 𝐹(𝑥⃗) + ∑ 𝜆𝑖 [𝑏𝑖 − 𝑔𝑖 (𝑥⃗)]
𝑖=1
Desarrollando la condición de primer orden para la función Lagrangiana
𝛻𝐿(𝑎⃗, ⃗⃗⃗⃗⃗ ⃗⃗
𝜆0 ) = 0
Desarrollando la condición de segundo orden de la función Lagrangiana:
0𝑚 𝐽𝑔⃗(𝑎⃗)
𝐻𝐿(𝑎⃗, ⃗⃗⃗⃗⃗
𝜆0 ) = ( )
[𝐽𝑔⃗(𝑎⃗)]𝑡 𝐻𝑓(𝑎⃗)
Se observa que se tiene una matriz orlada, por lo que se analiza su signo de los últimos (𝑛 − 𝑚)
menores principales dominantes de 𝐻, empezando con el determinante de 𝐻.
a) Si el determinante de 𝐻 tiene el mismo signo que (−1)𝑛 y si sus últimos (𝑛 − 𝑚) MPD
alternan de signo, entonces 𝐻 es DN en las restricciones dadas. Por lo tanto , ese punto
es máximo relativo
b) Si el determinante de 𝐻 y todos los últimos (𝑛 − 𝑚) MPD tienen el mismo signo que
(−1)𝑚 , entonces H es DP en las restricciones dadas. Por lo tanto , ese punto es mínimo
relativo.
USADO MULTIPLICADORES DE LAGRANGE
Solución
𝑦1 = 𝑦2 + 𝑦3
INCUBADORA PONEDORA
INSUMO 𝑥1 𝑦2
PRODUCTO 𝑦1 𝑥2
FUNCIÓN 1/2 1/2
𝑦1 = 4𝑥1 𝑥2 = 𝑦2
1° Incubadora 𝑀𝑎𝑥 𝐵1 = 𝑝𝑦1 − 1𝑥1
s.a
1/2
𝑦1 = 4𝑥1
Función Lagrangiana
1/2
𝐿1 (𝑥1 , 𝑦1 , 𝜆1 ) = 𝑝𝑦1 − 𝑥1 − 𝜆1 (𝑦1 − 4𝑥1 )
Condición de primer orden
𝜕𝐿1
𝜕𝑥1 1 1 1/2
−1 − 𝜆1 (−2𝑥−1/2
1 ) 0
𝜆1 = −1/2 = 𝑥1
2
𝜕𝐿1 2 𝑥1
∇𝐿1 (𝑥1 , 𝑦1 , 𝜆1 ) = =( 𝑝 − 𝜆1 ) = (0) →
𝜕𝑦1 1/2
𝑝 = 𝜆1
−(𝑦1 − 4𝑥1 ) 0 1/2
𝜕𝐿1 𝑦1 = 4𝑥1
(𝜕𝜆1 )
De 1 y 2
1 1/2 1/2
𝑝= 𝑥1 → 𝑥1 = 2𝑝, 𝑥1 = 4𝑝2 → 𝑒𝑛 3 𝑦1 = 4(2𝑝) = 8𝑝
2
Condición de segundo orden:
1
−
0 𝑔𝑥1 𝑔𝑦1 0 −2𝑥1 2 1
𝑔
𝐻𝐿1 (𝑥1 , 𝑦1 , 𝜆1 ) = [ 𝑥1 𝐿𝑥1 𝑥1 𝐿𝑥1 𝑦1 ] = −
1
−
3
2
𝑔𝑦1 𝐿𝑦1 𝑥1 𝐿𝑦1 𝑦1 −2𝑥 1 −1𝑥1 2 0
[ 1 0 0]
1 1
0 −2 ( ) 1 0 − 1
2𝑝 𝑝
1 1 3 = 1 1
− −1 ( ) 0 − − 3 0
𝑝 2𝑝 𝑝 8𝑝
[ 1 0 0] [ 1 0 0]
1
|𝐻3 | = > 0 𝑦𝑎 𝑞𝑢𝑒 𝑝 > 0 → 𝑀á𝑥𝑖𝑚𝑜
8𝑝 3
|𝐻3 | > 0, como representa el signo de (−1)𝑛 = (−1)2 = 1, entonces 𝐻𝐿 es definida
negativa. Por lo tanto, para
2° Ponedora 𝑀𝑎𝑥 𝐵2 = 1𝑥2 − 𝑝𝑦2
s.a
1/2
𝑥2 = 𝑦2
Función Lagrangiana
1/2
𝐿1 (𝑥2 , 𝑦2 , 𝜆2 ) = 𝑥2 − 𝑝𝑦2 − 𝜆2 (𝑥2 − 𝑦2 )
Condición de primer orden
𝜕𝐿2
𝜕𝑥2 1 − 𝜆2 𝜆2 = 1
𝜕𝐿2 1 −1/2
0 1 −1/2
∇𝐿2 (𝑥2 , 𝑦2 , 𝜆2 ) = = −𝑝 + 𝜆2 𝑦2 = (0) → 𝑝 = 𝜆2 𝑦2
𝜕𝑦2 2 2
1/2 0 1/2
𝜕𝐿2 −(𝑥2 − 𝑦2 ) 𝑥2 = 𝑦2
( )
(𝜕𝜆2 )
De 1 y 2
1 −1/2 1/2 1 1
𝑝 = 𝑦2 → 𝑦2 = → 𝑦2 = 2
2 2𝑝 4𝑝
En 2
1
𝑥2 =
2𝑝
Condición de segundo orden:
1
0 𝑔𝑥2 𝑔𝑦2 0 1 − 𝑦 −1/2
2
𝑔
𝐻𝐿2 (𝑥2 , 𝑦2 , 𝜆2 ) = [ 𝑥2 𝐿 𝑥2 𝑥2 𝐿𝑥2 𝑦2 ] = 1 0 0
𝑔𝑦2 𝐿𝑦2 𝑥2 𝐿𝑦2 𝑦2 1 −1/2 1 1 −3/2
− 𝑦 0 (− ) 𝜆2 𝑦2
[ 2 2 2 ]
1
0 1 − (2𝑝) 0 1 −𝑝
2
1 0 0 =[1 0 0 ]
1 1 −𝑝 0 −2𝑝 3
[− 2 (2𝑝) 0 − (2𝑝)3 ]
4
|𝐻3 | = 2𝑝 3 > 0 𝑦𝑎 𝑞𝑢𝑒 𝑝 > 0
|𝐻3 | > 0, como representa el signo de (−1)𝑛 = (−1)2 = 1, entonces 𝐻𝐿 es definida
negativa. Por lo tanto, para
Por tanto la función es cóncava, y los valores se optimizan para los valores de 𝑥1 , 𝑥2 . Ahora
calculamos el valor de 𝑝
𝑥1 = 𝑥2
1
= 4𝑝 2
2𝑝
8𝑝 3 − 1 = 0
(2𝑝 − 1)(4𝑝 2 + 2𝑝 + 1) = 0 → 𝑝 = 1/2
La cantidad de pollos dispuesta para el consumo
𝑦3 = 𝑦1 − 𝑦2
1
𝑦3 = 8𝑝 −
4𝑝 2
1
𝑦3 = 8 ( ) − 1 = 3
2
2. Un consumidor tiene como utilidad la función 𝑈 = ln 𝐶 + ln(24 − 𝑁), donde 𝐶 es
consumo y 𝑁 es oferta de trabajo. La restricción presupuestaria es 𝑝𝐶 = 𝑀 ̅ + 𝑤𝑁, donde
𝑝 es el precio del bien consumido, 𝑤 la tasa salarial y el ingreso del consumidor a parte del
salario.
a. Formule el problema de maximización de la utilidad sujeta a la restricción
presupuestaria, y derive las condiciones de primer orden usando el método del
multiplicador de Lagrange ignorando las restricciones no negativas.
b. Halle la función de demanda 𝐶 = 𝐶 ∗ (𝑝, 𝑤, 𝑀̅ ) y la función de oferta de trabajo 𝑁 =
𝑁 ∗ (𝑝, ̅ ) es decir, exprese 𝐶 𝑦 𝑁 en términos de 𝑝, 𝑤 𝑦 𝑀
𝑤, 𝑀 ̅ . Muestre que 𝐶 ∗ 𝑦 𝑁 ∗ y son
homogéneas de grado cero en (𝑝, 𝑤, 𝑀 ̅) .
̅ ) + ln(24 − 𝑁 ∗ (𝑝, 𝑤, 𝑀
c. Si 𝑈 ∗ = ln 𝐶 ∗ (𝑝, 𝑤, 𝑀 ̅ ) ).
𝜕𝑈 ∗ 𝜕𝑈 ∗
• Muestre que ̅
>0𝑦 < 0.
𝜕𝑀 𝜕𝑝
𝜕𝑈 ∗
• ¿Cuál es la relación entre ̅ y el multiplicador de Lagrange?
𝜕𝑀
Solución:
Parte a:
Max 𝑈(𝐶. 𝑁) = ln 𝐶 + ln(24 − 𝑁)
𝐷𝑜𝑚 = {(𝐶, 𝑁)𝜖 ℝ⁄𝑐 > 0 ∧ 24 − 𝑁 > 0}
𝐷𝑜𝑚 = {(𝐶, 𝑁)𝜖 ℝ⁄𝑐 > 0 ∧ 𝑁 < 24}
s.a
̅ + 𝑤𝑁
𝑝𝐶 = 𝑀
Función Lagrangiana
̅ − 𝑤𝑁)
𝐿(𝜆, 𝐶, 𝑁) = ln 𝐶 + ln(24 − 𝑁) − 𝜆(𝑝𝐶 − 𝑀
Condición de primer orden
𝜕𝐿 ̅ − 𝑤𝑁
𝑝𝐶 − 𝑀
𝜕𝜆 1
𝜕𝐿 0
∇𝐿(𝜆, 𝐶, 𝑁) = = − 𝜆𝑝 = (0)
𝐶
𝜕𝐶 1 0
𝜕𝐿
(− 24 − 𝑁 + 𝜆𝑤 )
(𝜕𝑁 )
De 2 y 3
1 1
=𝜆 ∧ = 𝜆 → 𝐶𝑝 = (24 − 𝑁)𝑤
𝐶𝑝 (24 − 𝑁)𝑤
Reemplazamos en 1
̅
24𝑤 − 𝑀 24𝑤 + 𝑀
̅ + 𝑤𝑁 → 𝑁∗ =
(24 − 𝑁)𝑤 = 𝑀 ∧ 𝐶∗ =
2𝑤 2𝑝
2
Y 𝜆=
24𝑤+𝑀
Condición de segundo orden:
0 𝐿𝜆𝐶 𝐿𝜆𝑁 0 𝑝 −𝑤
2
𝐻𝐿(𝜆, 𝐶, 𝑁) = [ 𝐿𝐶𝜆 𝐿𝐶𝐶 𝐿𝐶𝑁 ] = [ 𝑝 −1/𝐶 0 ]
𝐿𝑁𝜆 𝐿𝑁𝐶 𝐿𝑁𝑁 −𝑤 0 −1/(24 − 𝑁) 2
𝑝 𝑤 𝑝2 𝑤2
|𝐻3 | = −𝑝 ∗ [− ] − 𝑤 ∗ [− ] = + >0
(24 − 𝑁)2 𝐶2 (24 − 𝑁)2 𝐶 2
Se analizan los últimos (𝑛 − 𝑚) menores principales dominantes, es decir (2 − 1), el
último MPD.
|𝐻3 | > 0, como representa el signo de (−1)𝑛 = (−1)2 = 1, entonces 𝐻𝐿 es definida
negativa. Por lo tanto, para
̅
24𝑤 − 𝑀 ̅
24𝑤 + 𝑀
𝑁∗ = ∧ 𝐶∗ = , 𝑈 𝑒𝑠 𝑚á𝑥𝑖𝑚𝑜.
2𝑤 2𝑝
Parte b:
̅
24𝑤+𝑀
▪ ̅) =
𝐶 ∗(𝑝, 𝑤, 𝑀 2𝑝
Solución:
̅
24(𝑡𝑤)+𝑡𝑀 ̅)
𝑡(24𝑤+𝑀 ̅
24𝑤+𝑀
▪ 𝐶 ∗(𝑡𝑝, 𝑡𝑤, ̅̅̅̅
𝑡𝑀) = 2𝑡𝑝
= 𝑡(2𝑝) = 𝑡 0 [ 2𝑝 ]
por lo tanto, es homogénea de grado cero.
̅ ̅ ̅
▪ ̅̅̅̅) = 24(𝑡𝑤)−𝑡𝑀 = 𝑡(24𝑤−𝑀) = 𝑡 0 [24𝑤−𝑀]
𝑁 ∗ (𝑡𝑝, 𝑡𝑤, 𝑡𝑀 2𝑡𝑝 𝑡(2𝑝) 2𝑝
por lo tanto, es homogénea de grado cero.
24𝑤+𝑀̅ 24𝑤−𝑀̅
̅ ) = ln (
Parte c: Sea: 𝑈 ∗ (𝑝, 𝑤, 𝑀 ) + ln (24 − ).
2𝑝 2𝑤
Solución:
̅
48𝑤 − 24𝑤 + 𝑀
𝑈 ∗ (𝑝, 𝑤, 𝑀 ̅ ) − ln (2p) + ln (
̅ ) = ln(24𝑤 + 𝑀 )
2𝑤
̅ ) − ln (2p) + ln(24𝑤 + 𝑀
̅ ) = ln(24𝑤 + 𝑀
𝑈 ∗ (𝑝, 𝑤, 𝑀 ̅ ) − ln (2w)
Entonces:
𝜕𝑈 ∗ 2 ̅ > 0.
▪ ̅
= ̅ >
0 𝑝𝑢𝑒𝑠 𝑤 > 0, 𝑀
𝜕𝑀 24𝑤+𝑀
𝜕𝑈 ∗ 1
▪ = − <0 𝑝𝑢𝑒𝑠 𝑝 > 0
𝜕𝑝 𝑝
𝜕𝑈 ∗
▪ 𝐿𝑎 𝑟𝑒𝑙𝑎𝑐𝑖ó𝑛 𝑑𝑒 ̅
𝑦 𝜆.
𝜕𝑀
𝜕𝑈 ∗ 2
𝐷𝑒 𝑑𝑜𝑛𝑑𝑒 𝑡𝑒𝑛𝑒𝑚𝑜𝑠 𝑞𝑢𝑒: ̅
=𝜆= ̅ , son iguales.
𝜕𝑀 24𝑤+𝑀
3. Ejemplo:
𝑦
max 𝑓(𝑥, 𝑦) = ∫ (1 + 𝑡 − 𝛼𝑦 2 )(𝑥𝑡)(𝑦 − 𝑡)𝑑𝑡
0
s.a
𝑦
𝑔(𝑥, 𝑦) = ∫ (𝑥𝑡)(𝑦 − 𝑡)𝑑𝑡
0
Solución
Calculemos la función objetivo
𝑦 𝑦
∫ (1 + 𝑡 − 𝛼𝑦 )(𝑥𝑡)(𝑦 − 𝑡)𝑑𝑡 = 𝑥 ∫ [𝑡𝑦 + 𝑦𝑡2 − 𝛼𝑡𝑦3 − 𝑡2 − 𝑡3 + 𝛼𝑦2 𝑡2 ]𝑑𝑡
2
0 0
𝑦 𝑦
3) 2
(𝑦 − 𝛼𝑦3 )𝑡2 (𝑦 − 1 + 𝛼𝑦2 )𝑡3 𝑡4
2 3
𝑥 ∫ [(𝑦 − 𝛼𝑦 𝑡 + (𝑦 − 1 + 𝛼𝑦 )𝑡 − 𝑡 ]𝑑𝑡 = 𝑥 [ + − ]
0 2 3 4 0
2𝑥𝑦 3 − 2𝛼𝑥𝑦 5 + 𝑥𝑦 4
max 𝑓(𝑥, 𝑦) =
12
Calculamos la restricción
𝑦 𝑦
𝑦𝑡 2 𝑡 3 𝑥𝑦 3
𝑔(𝑥, 𝑦) = ∫ (𝑥𝑡)(𝑦 − 𝑡)𝑑𝑡 = 𝑥 [ − ] = =𝑀
0 2 3 0 6
6𝑀
Despejamos 𝑥 = y lo reemplazamos en 𝑓
𝑦3
6𝑀 3 6𝑀 6𝑀
2( ) 𝑦 − 2𝛼 ( 3 ) 𝑦 5 + ( 3 ) 𝑦 4 𝑀𝑦
𝑦3 𝑦 𝑦
𝑓(𝑦) = = 𝑀 − 𝛼𝑀𝑦 2 +
12 2
Condiciones de primer orden
𝜕𝑓 𝑀 1
= −2𝛼𝑀𝑦 + = 0 → 𝑦 ∗ =
𝜕𝑦 2 4𝛼
Condiciones de segundo orden
𝜕2𝑓
= −2𝛼𝑀 < 0
𝜕𝑦 2
1
Entonces para 𝑦 = , 𝑥 = 384𝑀𝛼 3 , la función tiene un máximo.
4𝛼
Ejercicios adicionales
4) Determine, en cada caso, los máximos y mínimos de f sobre D
a ) f (x, y) = 4xy , D = {(x, y) ∈ R2 | x9 + y16 = 1 }
2 2
b ) f (x, y) = x2 y , D = {(x, y) ∈ R2 | x2 + 2y 2 = 6}
c ) f (x, y) = 4x3 + y 2 , D = {(x, y) ∈ R2 | 2x2 + y 2 = 1}
d ) f (x, y) = xy , D = {(x, y) ∈ R2 | xa2 + yb2 = 1 }
2 2
e ) f (x, y, z) = x + 3y − z , D = {(x, y, z) ∈ R3 | x2 + y 2 + z 2 = 4}
f ) f (x, y, z) = xyz , D = {(x, y, z) ∈ R3 | x2 + 2y 2 + 3z 2 = 6}
5) Resuelva los siguientes problemas indicando el valor óptimo:
a ) máx −x2 − y 2 − z 2 s.a. xyz = 4
b ) mı́n x2 + y 2 + z 2 s.a. x + y + z = a
6) Calcule la mínima distancia del origen a la hipérbola de ecuación xy = 1.
(Sugerencia: puede considerar el cuadrado de la distancia como función objetivo).
7) Al invertir x unidades de trabajo e y unidades de capital, un fabricante de relojes
0,4 0,6
puede producir f (x, y) = 50x y relojes. Determine el número máximo de relojes
que se pueden producir con un presupuesto de $20000 si la mano de obra cuesta
$100 por unidad y el capital cuesta $200 por unidad.