Multiplicadores de Lagrange

Nivel UniversitarioDificultad ★★★★★Concepto⌖ Ver en el mapa

¿Qué es?

En un óptimo de ff con la restricción g=0g = 0, los gradientes son paralelos: ∇f=λ∇g\nabla f = \lambda\nabla g. El multiplicador λ\lambda mide cuánto mejoraría el óptimo si se relajara la restricción.

¿Por qué existe?

Sustituir la restricción para eliminar una variable solo funciona en casos fáciles. El método de Lagrange conserva todas las variables y añade una incógnita por restricción, convirtiendo un problema con restricciones en un sistema de ecuaciones sin ellas, y las nuevas incógnitas resultan tener significado económico y físico.

Intuición

Camina por la curva g=0g = 0 mirando los conjuntos de nivel de ff. Mientras la curva cruza conjuntos de nivel, todavía puedes bajar. En el óptimo la curva es tangente a un conjunto de nivel de ff, así que sus normales (∇f\nabla f y ∇g\nabla g) apuntan en la misma recta.

Definición formal

Si x∗x^\ast es un extremo local de ff sobre {g=0}\{g = 0\}, con f,gf, g de clase C1C^1 y ∇g(x∗)≠0\nabla g(x^\ast) \ne 0, entonces existe λ\lambda con ∇f(x∗)=λ ∇g(x∗)\nabla f(x^\ast) = \lambda\,\nabla g(x^\ast). Equivalentemente, (x∗,λ)(x^\ast, \lambda) es un punto crítico del lagrangiano

ℒ(x,λ)=f(x)−λ g(x).\mathcal L(x, \lambda) = f(x) - \lambda\,g(x).

Fórmulas

∇f(x∗)=λ ∇g(x∗),g(x∗)=0\nabla f(x^\ast) = \lambda\,\nabla g(x^\ast), \qquad g(x^\ast) = 0
∂f∗∂c=λ,g(x)=c\frac{\partial f^\ast}{\partial c} = \lambda, \qquad g(x) = c
el multiplicador es la sensibilidad del óptimo («precio sombra»)

Ejemplo

Maximiza la entropía H(p)=−∑ipiln⁡piH(p) = -\sum_i p_i\ln p_i sujeta a ∑ipi=1\sum_i p_i = 1. ∂pi(H−λ(∑p−1))=−ln⁡pi−1−λ=0\partial_{p_i}\big(H - \lambda(\sum p - 1)\big) = -\ln p_i - 1 - \lambda = 0, así que todos los pip_i son iguales: la distribución uniforme. Añade una restricción sobre la energía esperada y el mismo cálculo produce la distribución softmax o de Boltzmann pi∝e−βEip_i \propto e^{-\beta E_i}.

¿Por qué importa?

Las máquinas de vectores soporte se deducen a través de su dual lagrangiano (donde aparece el truco del núcleo). El entrenamiento con restricciones, las restricciones de equidad, la justificación de la softmax por máxima entropía y la dualidad en investigación operativa usan la misma idea. En física, la otra gran idea de Lagrange (el lagrangiano de la mecánica) gobierna la dinámica de los robots.

Aplicaciones en informática

Dónde aparece en IA

  • Máquinas de vectores soporte (SVM)★★★★★fundamentalIA y machine learning

    El dual de la SVM se obtiene con multiplicadores de Lagrange; solo importan los puntos con αi>0\alpha_i > 0 (vectores soporte).

  • Regularización★★★★★avanzadaIA y machine learning

    Entrenar con penalización L+λ∥w∥2L + \lambda\norm w^2 es el lagrangiano de entrenar con una restricción de norma ∥w∥≤r\norm w \le r.

¿Dónde se utiliza?

Temas de informática a los que se llega desde aquí, con la cadena de ideas que lleva a ellos:

Qué depende de él

Ejercicios

1Cálculo directo

Maximiza f(x,y)=xyf(x,y) = xy sujeta a x+y=10x + y = 10 con un multiplicador de Lagrange.

Solución

(y,x)=λ(1,1)(y, x) = \lambda(1, 1) da x=y=λx = y = \lambda, y x+y=10x + y = 10 da x=y=5x = y = 5, f=25f = 25, λ=5\lambda = 5.

2IA

Demuestra que la distribución que maximiza la entropía con energía media fija ∑ipiEi=Eˉ\sum_i p_i E_i = \bar E tiene la forma pi∝e−βEip_i \propto e^{-\beta E_i}.

Solución

∂pi[−∑pln⁡p−λ(∑p−1)−β(∑pE−Eˉ)]=−ln⁡pi−1−λ−βEi=0\partial_{p_i}\big[-\sum p\ln p - \lambda(\sum p - 1) - \beta(\sum pE - \bar E)\big] = -\ln p_i - 1 - \lambda - \beta E_i = 0, así que pi=e−1−λe−βEip_i = e^{-1-\lambda}e^{-\beta E_i}: una softmax de −βE-\beta E.

↑ ↓ para navegar · ↵ · Esc