Por decirlo de otra manera, nuestro objetivo ahora es optimizar la función f-λg con respecto a x, y, y z sin restricciones encontrando donde ∇(f-λg)=0.
¿Y si hay más de una restricción? Supongamos ahora que tenemos que optimizar f(x,y,x) sujeto a g(x,y,z)=c y h(x,y,z)=d. Esto es lo que hacemos. Para optimizar f sujeta a g(x,y,z)=c, optimizamos f-λg sin restricción, por lo que la primera restricción se soluciona sustituyendo f por f-λg, y queda optimizar f-λg sujeta a la restricción h(x,y,z)=d. Este es un problema con una sola restricción, así que simplemente añadimos otro multiplicador de Lagrange μ, y el problema ahora es optimizar f-λg-μh sin restricción.
Veamos ahora algunos ejemplos.
El área máxima de un rectángulo
Para un perímetro dado, ¿cuál es la mayor área posible de un rectángulo con ese perímetro? Podemos formularlo como un problema de multiplicador de Lagrange. Si el ancho y el alto son x e y, entonces deseamos maximizar f(x,y)=xy para g(x,y)=2x+2y=c. El sistema de ecuaciones resultante es:
Las dos primeras ecuaciones nos dicen enseguida que x=y, por lo que el área máxima se produce cuando el rectángulo es un cuadrado. Al introducir esto en la tercera ecuación, vemos que x=y=c/4, por lo que el área máxima de un rectángulo con perímetro c es c²/16.
Un problema de Putnam (modificado)
El concurso matemático de Putnam es una competición con formato de examen que se ofrece cada diciembre a los estudiantes universitarios. Los estudiantes tienen seis horas para resolver 12 problemas extremadamente difíciles. El examen es tan difícil que, sobre una puntuación posible de 120, la nota media en la mayoría de los años es de 0 o 1. Hoy veremos una versión modificada (para simplificar y evitar problemas de copyright) del problema A3 del examen de 2018:
Hay dos problemas. El primero es que la función objetivo tiene una característica muy indeseable: es una suma de cosenos, y las sumas de funciones trigonométricas pueden ser difíciles de tratar. El segundo es que no hay forma de hacer que la función de restricción aparezca de forma clara y explícita en la función objetivo, aunque a primera vista parezca que debería ser posible. Los problemas de Putnam suelen tener «trampas» como ésta, en las que se puede perder el tiempo con un planteamiento que no conducirá a una solución si no nos detenemos a pensar en lo que estamos haciendo. Es una trampa que sería especialmente peligrosa para un concursante con experiencia en concursos de matemáticas de secundaria, donde las identidades trigonométricas suelen tener un fuerte énfasis.
Pero cuando vemos un problema que dice «maximizar f sujeto a g=0», debemos pensar inmediatamente en los multiplicadores de Lagrange. Establezcamos nuestro sistema de ecuaciones:
Ahora podemos usar nuestras identidades trigonométricas para que la restricción aparezca explícitamente. Reescribe el sistema de ecuaciones utilizando la fórmula del doble ángulo:
El resultado es:
Ahora cancela la función seno en cada lado, y suma las ecuaciones:
Y esto significa que λ=0, así que según las ecuaciones del multiplicador de Lagrange esto significa que:
Por lo tanto w, x, y, y z son todos múltiplos enteros de π/2. Para los múltiplos enteros pares de π/2, digamos que m=2k para algún entero k:
Pero si m es un entero impar, entonces cos(mπ/2)=0. Esto significa que tenemos tres opciones de cómo podemos satisfacer la ecuación de restricción. En primer lugar, podemos establecer w, x, y y z todos a múltiplos enteros impares de π/2, digamos k, l, m y n, entonces:
La primera línea está mostrando que la ecuación de la restricción se satisface, y la segunda línea es el resultado de introducir esos valores en la función objetivo, que nos da -4.
La segunda opción es que podemos establecer dos de las variables a múltiplos enteros impares de π/2, digamos k y l, y las otras dos a enteros pares, digamos 2m y 2n donde m es par y n es impar. El resultado de insertar esto en la ecuación de la restricción es:
Así que la ecuación de la restricción se satisface. Aquí está el resultado de enchufar estos en la función objetivo:
La última opción es que podemos poner todas las variables a múltiplos enteros pares de π/2, digamos 2k, 2l, 2m y 2n, donde k y l son pares y m y n son impares. El resultado de enchufar esto en la ecuación de la restricción es:
Así que la restricción se satisface. Ahora enchufamos esto en la función objetivo para obtener:
Así que la respuesta es 4, que casualmente resulta ser el valor máximo global de la función objetivo.
Vea este otro artículo que escribí para una discusión de otro problema de Putnam.
Teorema de Thomson
Ahora, en lugar de optimizar una función, vamos a utilizar el método de los multiplicadores de Lagrange para optimizar un funcional: un objeto que toma una función y devuelve un número. Las integrales definidas son un ejemplo clave. Esta demostración final mostrará cómo se puede utilizar el método de los multiplicadores de Lagrange para encontrar la función que minimiza el valor de una integral definida.
Como demostración, consideramos un hecho fundamental de la electrostática, el estudio de los sistemas sometidos a fuerzas eléctricas en equilibrio mecánico de modo que ninguna de las cargas está en movimiento. Demostraremos que cuando un conductor, un cuerpo en el que la carga eléctrica puede moverse libremente, está en equilibrio con cualquier sistema de fuerzas eléctricas, entonces no hay campo eléctrico dentro del cuerpo del conductor. Este es el Teorema de Thomson. Este teorema se suele enunciar en términos del hecho físico fundamental de que la energía de un sistema en equilibrio mecánico estable es mínima:
Teorema de Thomson: La energía electrostática de un cuerpo conductor de tamaño y forma fijos se minimiza cuando las cargas se distribuyen de forma que el potencial ϕ sea constante en el cuerpo de forma que E=-∇ϕ=0.
Sea ρ(r) la función de distribución de cargas y ϕₑₓₜ el potencial debido a cualquier fuente externa al cuerpo conductor. La energía electrostática total puede expresarse como una funcional U de ρ(r):
Son integrales de volumen sobre el cuerpo del conductor. La primera integral es una doble integral realizada primero sobre las coordenadas cebadas y luego sobre las coordenadas no cebadas. Queremos minimizar U como una funcional de ρ(r) dado que la carga total es independiente de la distribución de carga:
Afortunadamente, el método de los multiplicadores de Lagrange funciona para los funcionales con restricciones. Introduce un multiplicador indeterminado de Lagrange, y entonces en lugar de una optimización con restricciones sobre U, realiza una optimización sin restricciones sobre el siguiente funcional:
Por analogía con las derivadas de funciones, el objetivo aquí es encontrar la distribución de carga ρ tal que cuando ρ es variado por una pequeña variación δρ, el cambio en I, δI, es cero:
Esto no es del todo una derivada porque δρ es una función en lugar de un número por lo que no tiene realmente sentido hablar del límite cuando δρ llega a cero, pero la intuición básica es la misma.
Empecemos por expandir U:
Esto puede parecer malo, pero podemos simplificar esto en gran medida. Primero, como δρ es muy pequeño, (δρ)² es tan pequeño que podemos ignorarlo, así que supongamos que δρ(r)δρ(r′)=0. A continuación, cambiando el orden de integración vemos que:
También podemos ver que U aparece en la expansión. Esto nos da lo suficiente para simplificar U:
Ahora podemos expandir δI:
Ahora podemos factorizar la integral sobre δρ(r):
Para que esto sea cierto para cualquier elección de la variación δρ(r), el término entre corchetes debe ser cero, así que:
Pero esta integral es sólo el potencial en r debido a las cargas dentro del conductor por lo que el lado izquierdo de esta ecuación da el potencial total en cada punto r dentro del conductor, así que ρ(r) es la distribución que hace que el potencial total dentro del conductor sea constante. Esto completa la prueba.
Esto también sería válido si la energía electrostática del conductor se maximiza, en cuyo caso el conductor está en equilibrio inestable. En la práctica uno nunca encuentra conductores en equilibrio inestable porque el movimiento térmico aleatorio de los electrones en el conductor rápidamente, en el orden de 10-¹⁴ segundos para los conductores metálicos, conduce el sistema fuera del equilibrio inestable y en el equilibrio estable.