Método de Lagrange en Optimización
Método de Lagrange en Optimización
2022
restricciones.
restricciones.
Teorema
Si (λ0 , x0 , y0 ) es un punto crítico de la función de Lagrange
L = f + λ φ entonces se presenta una de las alternativas siguientes:
1 (x0 , y0 ) es máximo relativo de f sujeto a la restricción φ = 0.
Teorema
Si (λ0 , x0 , y0 ) es un punto crítico de la función de Lagrange
L = f + λ φ entonces se presenta una de las alternativas siguientes:
1 (x0 , y0 ) es máximo relativo de f sujeto a la restricción φ = 0.
0
φx φy
∆(λ , x, y ) = det φx Lxx Lxy (λ , x, y )
φy Lyx Lyy
0
φx φy
∆(λ , x, y ) = det φx Lxx Lxy (λ , x, y )
φy Lyx Lyy
0
φx φy
∆(λ , x, y ) = det φx Lxx Lxy (λ , x, y )
φy Lyx Lyy
Teorema
Sea (λ0 , x0 , y0 ) un punto crítico de la función de Lagrange
L = f + λ φ y sea ∆(λ , x, y ) el Hessiano Orlado asociado a f y a φ .
Se cumple:
1 Si ∆(λ0 , x0 , y0 ) > 0 entonces (x0 , y0 ) es un máximo relativo de
f sujeto a la restricción φ = 0.
2 Si ∆(λ0 , x0 , y0 ) < 0 entonces (x0 , y0 ) es un mínimo relativo de
f sujeto a la restricción φ = 0.
Observación:
Si ∆(λ0 , x0 , y0 ) = 0 entonces el criterio no da información.
Teorema
Sea (λ0 , x0 , y0 ) un punto crítico de la función de Lagrange
L = f + λ φ y sea ∆(λ , x, y ) el Hessiano Orlado asociado a f y a φ .
Se cumple:
1 Si ∆(λ0 , x0 , y0 ) > 0 entonces (x0 , y0 ) es un máximo relativo de
f sujeto a la restricción φ = 0.
2 Si ∆(λ0 , x0 , y0 ) < 0 entonces (x0 , y0 ) es un mínimo relativo de
f sujeto a la restricción φ = 0.
Observación:
Si ∆(λ0 , x0 , y0 ) = 0 entonces el criterio no da información.
Ejemplo 1
Hallar los valores máximos y mínimos de la función
f (x, y ) = 2x 2 + y 2 ,
sujeta a la restricción
x + y = 1.
L(λ , x, y ) = 2x 2 + y 2 + λ (x + y − 1)
Ejemplo 1
Hallar los valores máximos y mínimos de la función
f (x, y ) = 2x 2 + y 2 ,
sujeta a la restricción
x + y = 1.
L(λ , x, y ) = 2x 2 + y 2 + λ (x + y − 1)
Ejemplo 1
Hallar los valores máximos y mínimos de la función
f (x, y ) = 2x 2 + y 2 ,
sujeta a la restricción
x + y = 1.
L(λ , x, y ) = 2x 2 + y 2 + λ (x + y − 1)
Ejemplo 1
Hallar los valores máximos y mínimos de la función
f (x, y ) = 2x 2 + y 2 ,
sujeta a la restricción
x + y = 1.
L(λ , x, y ) = 2x 2 + y 2 + λ (x + y − 1)
Ejemplo 1
Hallar los valores máximos y mínimos de la función
f (x, y ) = 2x 2 + y 2 ,
sujeta a la restricción
x + y = 1.
L(λ , x, y ) = 2x 2 + y 2 + λ (x + y − 1)
y = 2x (4)
En (1): x + 2x − 1 = 0 , luego x = 13 .
En (4): y = 23 .
En (2): λ = − 43
De esta manera, el único PC es − 43 , 31 , 23
y = 2x (4)
En (1): x + 2x − 1 = 0 , luego x = 13 .
En (4): y = 23 .
En (2): λ = − 43
De esta manera, el único PC es − 43 , 31 , 23
y = 2x (4)
En (1): x + 2x − 1 = 0 , luego x = 13 .
En (4): y = 23 .
En (2): λ = − 43
De esta manera, el único PC es − 43 , 31 , 23
y = 2x (4)
En (1): x + 2x − 1 = 0 , luego x = 13 .
En (4): y = 23 .
En (2): λ = − 43
De esta manera, el único PC es − 43 , 13 , 23
y = 2x (4)
En (1): x + 2x − 1 = 0 , luego x = 13 .
En (4): y = 23 .
En (2): λ = − 43
De esta manera, el único PC es − 43 , 31 , 23
y = 2x (4)
En (1): x + 2x − 1 = 0 , luego x = 13 .
En (4): y = 23 .
En (2): λ = − 43
De esta manera, el único PC es − 43 , 31 , 23
y = 2x (4)
En (1): x + 2x − 1 = 0 , luego x = 13 .
En (4): y = 23 .
En (2): λ = − 43
De esta manera, el único PC es − 43 , 31 , 23
y = 2x (4)
En (1): x + 2x − 1 = 0 , luego x = 13 .
En (4): y = 23 .
En (2): λ = − 43
De esta manera, el único PC es − 43 , 31 , 23
y = 2x (4)
En (1): x + 2x − 1 = 0 , luego x = 13 .
En (4): y = 23 .
En (2): λ = − 43
De esta manera, el único PC es − 43 , 31 , 23
Ejemplo 2
La función de producción de una empresa es
Ejemplo 2
La función de producción de una empresa es
Ejemplo 2
La función de producción de una empresa es
0 = Lλ (λ , x1 , x2 , . . . , xn ) = φ (x1 , x2 , . . . , xn )
0 = Lx1 (λ , x1 , x2 , . . . , xn ) = fx1 (x1 , . . . , xn ) + λ φx1 (x1 , . . . , xn )
.. .. ..
. . .
0 =
Lxn (λ , x1 , . . . , xn ) = fxn (x1 , . . . , xn ) + λ φxn (x1 , . . . , xn )
0 = Lλ (λ , x1 , x2 , . . . , xn ) = φ (x1 , x2 , . . . , xn )
0 = Lx1 (λ , x1 , x2 , . . . , xn ) = fx1 (x1 , . . . , xn ) + λ φx1 (x1 , . . . , xn )
.. .. ..
. . .
0 =
Lxn (λ , x1 , . . . , xn ) = fxn (x1 , . . . , xn ) + λ φxn (x1 , . . . , xn )
0 = Lλ (λ , x1 , x2 , . . . , xn ) = φ (x1 , x2 , . . . , xn )
0 = Lx1 (λ , x1 , x2 , . . . , xn ) = fx1 (x1 , . . . , xn ) + λ φx1 (x1 , . . . , xn )
.. .. ..
. . .
0 =
Lxn (λ , x1 , . . . , xn ) = fxn (x1 , . . . , xn ) + λ φxn (x1 , . . . , xn )
Denición
El Hessiano Orlado asociado a f y a φ , se dene como
0
φx1 φx2 ··· φxn
φx1 Lx1 x1 Lx1 x2 ··· Lx1 xn
∆n+1 = det
φx2 Lx2 x1 Lx2 x2 ··· Lx2 xn
.. .. .. ..
. . . .
φxn Lxn x1 Lxn x2 ··· Lxn xn
Denimos también
0
φx1 φx2
∆3 = det φx1 Lx1 x1 Lx1 x2
φx2 Lx2 x1 Lx2 x2
Denición
El Hessiano Orlado asociado a f y a φ , se dene como
0
φx1 φx2 ··· φxn
φx1 Lx1 x1 Lx1 x2 ··· Lx1 xn
∆n+1 = det
φx2 Lx2 x1 Lx2 x2 ··· Lx2 xn
.. .. .. ..
. . . .
φxn Lxn x1 Lxn x2 ··· Lxn xn
Denimos también
0
φx1 φx2
∆3 = det φx1 Lx1 x1 Lx1 x2
φx2 Lx2 x1 Lx2 x2
0
φx1 φx2 φx3
φx1 Lx1 x1 Lx1 x2 Lx1 x3
∆4 = det
φx2 Lx2 x1 Lx2 x2 Lx2 x3
φx3 Lx3 x1 Lx3 x2 Lx3 x3
Análogamente se denen ∆5 , . . . , ∆n .
Observación:
En general ∆3 , ∆4 , . . . , ∆n+1 dependen de λ , x1 , . . . , xn .
La condición de segundo orden para problemas de optimización de
varias variables con una sola restricción, viene dada por el siguiente
teorema.
0
φx1 φx2 φx3
φx1 Lx1 x1 Lx1 x2 Lx1 x3
∆4 = det
φx2 Lx2 x1 Lx2 x2 Lx2 x3
φx3 Lx3 x1 Lx3 x2 Lx3 x3
Análogamente se denen ∆5 , . . . , ∆n .
Observación:
En general ∆3 , ∆4 , . . . , ∆n+1 dependen de λ , x1 , . . . , xn .
La condición de segundo orden para problemas de optimización de
varias variables con una sola restricción, viene dada por el siguiente
teorema.
0
φx1 φx2 φx3
φx1 Lx1 x1 Lx1 x2 Lx1 x3
∆4 = det
φx2 Lx2 x1 Lx2 x2 Lx2 x3
φx3 Lx3 x1 Lx3 x2 Lx3 x3
Análogamente se denen ∆5 , . . . , ∆n .
Observación:
En general ∆3 , ∆4 , . . . , ∆n+1 dependen de λ , x1 , . . . , xn .
La condición de segundo orden para problemas de optimización de
varias variables con una sola restricción, viene dada por el siguiente
teorema.
Teorema
Si (λ0 , x10 , . . . , xn0 ) es un punto crítico de la función de Lagrange
L = f + λ φ . Con las notaciones anteriores, tenemos:
1 Si ∆3 (λ0 , x 0 , . . . , x 0 ) > 0, ∆4 (λ0 , x 0 , . . . , x 0 ) < 0, . . . ,
1 n 1 n
(−1)n ∆n+1 (λ0 , x10 , . . . , xn0 ) > 0 entonces (x10 , . . . , xn0 ) es un
máximo relativo de f sujeto a la restricción φ = 0.
2 Si ∆3 (λ0 , x 0 , . . . , x 0 ) < 0, ∆4 (λ0 , x 0 , . . . , x 0 ) < 0, . . . ,
1 n 1 n
∆n+1 (x10 , . . . , xn0 ) < 0 entonces (λ0 , x10 , . . . , xn0 ) es un mínimo
relativo de f sujeto a la restricción φ = 0.
Observación:
Con cualquier otra combinación de signos, o la presencia de algún
cero, el criterio no da información.
Teorema
Si (λ0 , x10 , . . . , xn0 ) es un punto crítico de la función de Lagrange
L = f + λ φ . Con las notaciones anteriores, tenemos:
1 Si ∆3 (λ0 , x 0 , . . . , x 0 ) > 0, ∆4 (λ0 , x 0 , . . . , x 0 ) < 0, . . . ,
1 n 1 n
(−1)n ∆n+1 (λ0 , x10 , . . . , xn0 ) > 0 entonces (x10 , . . . , xn0 ) es un
máximo relativo de f sujeto a la restricción φ = 0.
2 Si ∆3 (λ0 , x 0 , . . . , x 0 ) < 0, ∆4 (λ0 , x 0 , . . . , x 0 ) < 0, . . . ,
1 n 1 n
∆n+1 (x10 , . . . , xn0 ) < 0 entonces (λ0 , x10 , . . . , xn0 ) es un mínimo
relativo de f sujeto a la restricción φ = 0.
Observación:
Con cualquier otra combinación de signos, o la presencia de algún
cero, el criterio no da información.
Ejemplo 3.
Determine (si es que existen) los máximos y mínimos de la función
f (x, y , z) = x sujeto a la restricción x 2 + y 2 + z 2 = 1
L(λ , x, y , z) = x + λ (x 2 + y 2 + z 2 − 1)
Ejemplo 3.
Determine (si es que existen) los máximos y mínimos de la función
f (x, y , z) = x sujeto a la restricción x 2 + y 2 + z 2 = 1
L(λ , x, y , z) = x + λ (x 2 + y 2 + z 2 − 1)
Ejemplo 3.
Determine (si es que existen) los máximos y mínimos de la función
f (x, y , z) = x sujeto a la restricción x 2 + y 2 + z 2 = 1
L(λ , x, y , z) = x + λ (x 2 + y 2 + z 2 − 1)
Ejemplo 3.
Determine (si es que existen) los máximos y mínimos de la función
f (x, y , z) = x sujeto a la restricción x 2 + y 2 + z 2 = 1
L(λ , x, y , z) = x + λ (x 2 + y 2 + z 2 − 1)
Ejemplo 3.
Determine (si es que existen) los máximos y mínimos de la función
f (x, y , z) = x sujeto a la restricción x 2 + y 2 + z 2 = 1
L(λ , x, y , z) = x + λ (x 2 + y 2 + z 2 − 1)
Ejemplo 3.
Determine (si es que existen) los máximos y mínimos de la función
f (x, y , z) = x sujeto a la restricción x 2 + y 2 + z 2 = 1
L(λ , x, y , z) = x + λ (x 2 + y 2 + z 2 − 1)
Ejemplo 3.
Determine (si es que existen) los máximos y mínimos de la función
f (x, y , z) = x sujeto a la restricción x 2 + y 2 + z 2 = 1
L(λ , x, y , z) = x + λ (x 2 + y 2 + z 2 − 1)
De (9) y (10): y = 0.
De (9) y (11): z = 0.
Reemplazando estos dos últimos resultados en (8): x 2 = 1, luego
x = ±1.
En (9): Si x = 1 entonces λ = −1/2 y si x = −1 entonces λ = 1/2.
Luego los PC de L son (−1/2, 1, 0, 0) y (1/2, −1, 0, 0).
Paso 3: Clasicación del PC: El Hessiano orlado es
Lλ λ Lλ x Lλ y Lλ z
Lxλ Lxx Lxy Lxz
∆4 (λ , x, y , z) = det
Ly λ
Lyx Lyy Lyz
Lzλ Lzx Lzy Lzz
0 2x 2y 2z
2x 2λ 0 0
= det
2y 0 2λ 0
2z 0 0 2λ
De (9) y (10): y = 0.
De (9) y (11): z = 0.
Reemplazando estos dos últimos resultados en (8): x 2 = 1, luego
x = ±1.
En (9): Si x = 1 entonces λ = −1/2 y si x = −1 entonces λ = 1/2.
Luego los PC de L son (−1/2, 1, 0, 0) y (1/2, −1, 0, 0).
Paso 3: Clasicación del PC: El Hessiano orlado es
Lλ λ Lλ x Lλ y Lλ z
Lxλ Lxx Lxy Lxz
∆4 (λ , x, y , z) = det
Ly λ
Lyx Lyy Lyz
Lzλ Lzx Lzy Lzz
0 2x 2y 2z
2x 2λ 0 0
= det
2y 0 2λ 0
2z 0 0 2λ
De (9) y (10): y = 0.
De (9) y (11): z = 0.
Reemplazando estos dos últimos resultados en (8): x 2 = 1, luego
x = ±1.
En (9): Si x = 1 entonces λ = −1/2 y si x = −1 entonces λ = 1/2.
Luego los PC de L son (−1/2, 1, 0, 0) y (1/2, −1, 0, 0).
Paso 3: Clasicación del PC: El Hessiano orlado es
Lλ λ Lλ x Lλ y Lλ z
Lxλ Lxx Lxy Lxz
∆4 (λ , x, y , z) = det
Ly λ
Lyx Lyy Lyz
Lzλ Lzx Lzy Lzz
0 2x 2y 2z
2x 2λ 0 0
= det
2y 0 2λ 0
2z 0 0 2λ
De (9) y (10): y = 0.
De (9) y (11): z = 0.
Reemplazando estos dos últimos resultados en (8): x 2 = 1, luego
x = ±1.
En (9): Si x = 1 entonces λ = −1/2 y si x = −1 entonces λ = 1/2.
Luego los PC de L son (−1/2, 1, 0, 0) y (1/2, −1, 0, 0).
Paso 3: Clasicación del PC: El Hessiano orlado es
Lλ λ Lλ x Lλ y Lλ z
Lxλ Lxx Lxy Lxz
∆4 (λ , x, y , z) = det
Ly λ
Lyx Lyy Lyz
Lzλ Lzx Lzy Lzz
0 2x 2y 2z
2x 2λ 0 0
= det
2y 0 2λ 0
2z 0 0 2λ
De (9) y (10): y = 0.
De (9) y (11): z = 0.
Reemplazando estos dos últimos resultados en (8): x 2 = 1, luego
x = ±1.
En (9): Si x = 1 entonces λ = −1/2 y si x = −1 entonces λ = 1/2.
Luego los PC de L son (−1/2, 1, 0, 0) y (1/2, −1, 0, 0).
Paso 3: Clasicación del PC: El Hessiano orlado es
Lλ λ Lλ x Lλ y Lλ z
Lxλ Lxx Lxy Lxz
∆4 (λ , x, y , z) = det
Ly λ
Lyx Lyy Lyz
Lzλ Lzx Lzy Lzz
0 2x 2y 2z
2x 2λ 0 0
= det
2y 0 2λ 0
2z 0 0 2λ
De (9) y (10): y = 0.
De (9) y (11): z = 0.
Reemplazando estos dos últimos resultados en (8): x 2 = 1, luego
x = ±1.
En (9): Si x = 1 entonces λ = −1/2 y si x = −1 entonces λ = 1/2.
Luego los PC de L son (−1/2, 1, 0, 0) y (1/2, −1, 0, 0).
Paso 3: Clasicación del PC: El Hessiano orlado es
Lλ λ Lλ x Lλ y Lλ z
Lxλ Lxx Lxy Lxz
∆4 (λ , x, y , z) = det
Ly λ
Lyx Lyy Lyz
Lzλ Lzx Lzy Lzz
0 2x 2y 2z
2x 2λ 0 0
= det
2y 0 2λ 0
2z 0 0 2λ
De (9) y (10): y = 0.
De (9) y (11): z = 0.
Reemplazando estos dos últimos resultados en (8): x 2 = 1, luego
x = ±1.
En (9): Si x = 1 entonces λ = −1/2 y si x = −1 entonces λ = 1/2.
Luego los PC de L son (−1/2, 1, 0, 0) y (1/2, −1, 0, 0).
Paso 3: Clasicación del PC: El Hessiano orlado es
Lλ λ Lλ x Lλ y Lλ z
Lxλ Lxx Lxy Lxz
∆4 (λ , x, y , z) = det
Ly λ
Lyx Lyy Lyz
Lzλ Lzx Lzy Lzz
0 2x 2y 2z
2x 2λ 0 0
= det
2y 0 2λ 0
2z 0 0 2λ
0 0 0 −1
2 0 0
= −2 det 0 −1 0
0 0 −1
= −4 < 0
y
0 2 0
0 0 0 −1
2 0 0
= −2 det 0 −1 0
0 0 −1
= −4 < 0
y
0 2 0
0 0 0 −1
2 0 0
= −2 det 0 −1 0
0 0 −1
= −4 < 0
y
0 2 0
0 0 0 −1
2 0 0
= −2 det 0 −1 0
0 0 −1
= −4 < 0
y
0 2 0
0 0 0 −1
2 0 0
= −2 det 0 −1 0
0 0 −1
= −4 < 0
y
0 2 0
0 0 0 1
−2 0 0
= −(−2) det 0 1 0
0 0 1
= −4 < 0
y
0 −2 0
0 0 0 1
−2 0 0
= −(−2) det 0 1 0
0 0 1
= −4 < 0
y
0 −2 0
0 0 0 1
−2 0 0
= −(−2) det 0 1 0
0 0 1
= −4 < 0
y
0 −2 0
0 0 0 1
−2 0 0
= −(−2) det 0 1 0
0 0 1
= −4 < 0
y
0 −2 0
0 0 0 1
−2 0 0
= −(−2) det 0 1 0
0 0 1
= −4 < 0
y
0 −2 0
0 = Lλ1 = φ1 (x1 , . . . , xn )
.. .. ..
. . .
0 = Lλm = φm (x1 , . . . , xn )
m
∂ φj
0 = Lx1
= fx1 (x1 , . . . , xn ) + ∑ λj (x1 , . . . , xn )
j=1 ∂ x1
.. .. ..
. . .
m
∂ φj
0 = Lxn = fxn (x1 , x2 , . . . , xn ) + ∑ λj (x1 , . . . , xn )
j=1 ∂ xm
0 = Lλ1 = φ1 (x1 , . . . , xn )
.. .. ..
. . .
0 = Lλm = φm (x1 , . . . , xn )
m
∂ φj
0 = Lx1
= fx1 (x1 , . . . , xn ) + ∑ λj (x1 , . . . , xn )
j=1 ∂ x1
.. .. ..
. . .
m
∂ φj
0 = Lxn = fxn (x1 , x2 , . . . , xn ) + ∑ λj (x1 , . . . , xn )
j=1 ∂ xm
. . . , φm = 0.
2 x0 es un mínimo relativo de f sujeto a las restricciones φ1 = 0,
. . . , φm = 0.
3 x0 es un punto silla.
. . . , φm = 0.
2 x0 es un mínimo relativo de f sujeto a las restricciones φ1 = 0,
. . . , φm = 0.
3 x0 es un punto silla.
Denición:
El Hessiano Orlado asociado a f y a las restricciones φ1 = 0, . . . ,
φm = 0, se dene como
0 0 ∂ φ1 ∂ φ1
··· ∂ x1 ··· ∂ xn
.. .. .. ..
. . . .
0 ··· 0 ∂ φm
··· ∂ φm
∆m+n (λ , x) = det ∂ φ1 ∂ x1 ∂ xn
(λ , x)
∂ φm
∂x
1 ··· ∂ x1 Lx1 x1 ··· Lx1 xn
.. .. .. ..
. . . .
∂ φ1 ∂ φm
∂ xn ··· ∂ xn Lxn x1 ··· Lxn xn
Denimos también
0 ··· 0 ∂ φ1
∂ x1
∂ φ1
∂ x2
. .. ..
.
. . .
∆m+2 (λ , x) = det 0 ··· 0 ∂ φ1 ∂ φ1
(λ , x)
∂ x1 ∂ x2
∂ φ1 ··· ∂ φm
Lx1 x1 Lx1 x2
∂ x1 ∂ x1
∂ φ1 ∂ φm
∂ x1 ··· ∂ x1 Lx2 x1 Lx2 x2
Denimos también
0 ··· 0 ∂ φ1
∂ x1
∂ φ1
∂ x2
. .. ..
.
. . .
∆m+2 (λ , x) = det 0 ··· 0 ∂ φ1 ∂ φ1
(λ , x)
∂ x1 ∂ x2
∂ φ1 ··· ∂ φm
Lx1 x1 Lx1 x2
∂ x1 ∂ x1
∂ φ1 ∂ φm
∂ x1 ··· ∂ x1 Lx2 x1 Lx2 x2
Denimos también
0 ··· 0 ∂ φ1
∂ x1
∂ φ1
∂ x2
. .. ..
.
. . .
∆m+2 (λ , x) = det 0 ··· 0 ∂ φ1 ∂ φ1
(λ , x)
∂ x1 ∂ x2
∂ φ1 ··· ∂ φm
Lx1 x1 Lx1 x2
∂ x1 ∂ x1
∂ φ1 ∂ φm
∂ x1 ··· ∂ x1 Lx2 x1 Lx2 x2
Teorema
Si (λ0 , x0 ) es un punto crítico de la función de Lagrange
L = f + λ1 φ1 + · · · + λm φm . Con las notaciones anteriores, tenemos:
1 Si (−1)m+1 ∆m+2 (λ0 , x0 ) > 0, (−1)m+1 ∆m+3 (λ0 , x0 ) < 0,
Observación:
Con cualquier otra combinación de signos, o la presencia de algún
cero, el criterio no da información.
Teorema
Si (λ0 , x0 ) es un punto crítico de la función de Lagrange
L = f + λ1 φ1 + · · · + λm φm . Con las notaciones anteriores, tenemos:
1 Si (−1)m+1 ∆m+2 (λ0 , x0 ) > 0, (−1)m+1 ∆m+3 (λ0 , x0 ) < 0,
Observación:
Con cualquier otra combinación de signos, o la presencia de algún
cero, el criterio no da información.
∆4 (λ0 , P0 ) P0
- máximo
+ mínimo
∆4 (λ0 , P0 ) ∆5 (λ0 , P0 ) P0
- + máximo
+ + mínimo
∆4 (λ0 , P0 ) P0
- máximo
+ mínimo
∆4 (λ0 , P0 ) ∆5 (λ0 , P0 ) P0
- + máximo
+ + mínimo