Matly

Multiplicadors de Lagrange: restriccions d'igualtat i qualificació

Vols practicar? Aquest tema té 112 exercicis interactius de 3 tipus, de deu en deu, amb correcció immediata i explicació pas a pas.

Quan el domini no és tot l'espai sinó una superfície definida per equacions, el gradient de la funció objectiu no s'anul·la a l'òptim: només s'anul·la la seva component tangent a la restricció. Aquesta frase és tot el mètode de Lagrange; la resta és àlgebra.

La geometria

Vols minimitzar una funció subjecta a unes quantes equacions h(x)=0h(x) = 0, amb menys equacions que variables. Si la matriu jacobiana de les restriccions té rang màxim al punt, el conjunt factible és localment una superfície llisa i el seu espai tangent és el nucli d'aquesta jacobiana.

En un extrem, la derivada de la funció objectiu al llarg de qualsevol direcció tangent ha de ser zero. Això vol dir que el gradient és perpendicular a tot l'espai tangent, i per tant és una combinació lineal dels gradients de les restriccions, que són justament les direccions normals.

Teorema de Lagrange

Si el punt és un extrem local sobre el conjunt factible, les funcions són de classe C1C^{1} i els gradients de les restriccions són linealment independents al punt, aleshores existeixen multiplicadors tals que

∇f(x∗)+Dh(x∗)Tλ=0,h(x∗)=0.\nabla f(x_*) + Dh(x_*)^{T}\lambda = 0, \qquad h(x_*) = 0.

Equivalentment, el gradient del lagrangià L=f+λTh\mathcal{L} = f + \lambda^{T}h respecte de totes les variables s'anul·la.

El signe amb què escrius els multiplicadors és una convenció: canviar-lo només els canvia el signe, no els punts candidats. El que no és convenció és la hipòtesi de rang.

Un exemple complet

Extrems del producte sobre la circumferència

Busca els extrems de f(x,y)=xyf(x,y) = xy sobre la circumferència unitat.

Pas 0 — existència i regularitat. La circumferència és compacta i la funció és contínua: per Weierstrass hi ha màxim i mínim. I el gradient de la restricció, (2x,2y)(2x,2y), no s'anul·la enlloc de la circumferència, de manera que la qualificació es compleix a tot arreu. Aquestes dues comprovacions van primer.

Pas 1 — les equacions.

y+2λx=0,x+2λy=0,x2+y2=1.y + 2\lambda x = 0, \qquad x + 2\lambda y = 0, \qquad x^{2}+y^{2} = 1.

Pas 2 — resoldre amb cura. Si xx fos zero, la primera equació donaria y=0y = 0, cosa incompatible amb la restricció. Igual amb yy. Per tant cap dels dos s'anul·la.

Multiplicant la primera per xx, la segona per yy i restant surt x2=y2x^{2} = y^{2}.

Pas 3 — candidats i valors. Amb x=y=±1/2x = y = \pm 1/\sqrt2 el producte val 1/21/2; amb x=−y=±1/2x = -y = \pm 1/\sqrt2 val −1/2-1/2. Com que sabem que hi ha màxim i mínim i aquests són tots els candidats, el màxim és 1/21/2 i el mínim −1/2-1/2.

Comprovació independent. De la desigualtat 2∣xy∣≤x2+y2=12|xy|\leq x^{2}+y^{2} = 1 surt directament que el producte està entre −1/2-1/2 i 1/21/2. Coincideix.

Dos gradients paral·lels, vist a la circumferència

Busca el punt de la circumferència unitat més llunyà del punt (2,0)(2,0).

Objectiu. Maximitzar el quadrat de la distància, f(x,y)=(x−2)2+y2f(x,y) = (x-2)^{2}+y^{2}, sobre x2+y2=1x^{2}+y^{2} = 1.

Equacions. Els gradients són (2(x−2), 2y)(2(x-2),\, 2y) i (2x, 2y)(2x,\, 2y). La condició de paral·lelisme dona

2(x−2)=2λx,2y=2λy.2(x-2) = 2\lambda x, \qquad 2y = 2\lambda y.

Resolent. De la segona, o bé y=0y = 0, o bé λ=1\lambda = 1. Si λ=1\lambda = 1, la primera dona −4=0-4 = 0, impossible. Per tant y=0y = 0, i amb la restricció x=±1x = \pm 1.

Comparant valors. A (1,0)(1,0) la distància al quadrat val 11; a (−1,0)(-1,0) val 99. El punt més llunyà és (−1,0)(-1,0), a distància 33.

La comprovació geomètrica. Els dos punts són els extrems del diàmetre alineat amb (2,0)(2,0), que és exactament el que diu la intuïció: el més proper i el més llunyà estan sobre la recta que uneix el punt amb el centre. Quan Lagrange i la geometria coincideixen, pots estar tranquil.

Quan el mètode no s'aplica

Sense qualificació, Lagrange pot no tenir solució

Minimitza la funció identitat subjecta a x2=0x^{2} = 0.

L'únic punt factible és l'origen, que per tant és alhora mínim i màxim. Però la condició de Lagrange seria 1+λ⋅0=01 + \lambda\cdot 0 = 0, que no té solució per a cap multiplicador.

Què ha fallat? El gradient de la restricció s'anul·la al punt, de manera que la hipòtesi de rang no es compleix. El teorema no diu res sobre aquest punt, i per tant no hi ha cap contradicció: simplement no s'hi aplica.

La conseqüència pràctica: els punts on falla el rang s'han d'estudiar a part, amb un altre argument. Si et limites a resoldre les equacions de Lagrange, te'ls deixes.

Segon ordre amb restriccions

Amb funcions de classe C2C^{2} i regularitat, la condició necessària de mínim és que la hessiana del lagrangià respecte de les variables sigui semidefinida positiva restringida a l'espai tangent. I la positivitat estricta sobre totes les direccions tangents no nul·les és suficient per a mínim local estricte.

Dues coses que aquí es fan malament:

  • Cal la hessiana del lagrangià, no la de la funció objectiu sola. La curvatura de la restricció hi entra a través dels multiplicadors, i ignorar-la pot invertir la conclusió.
  • Cal restringir-la a l'espai tangent. Una hessiana indefinida a tot l'espai pot ser perfectament definida positiva sobre el subespai que importa.

Desigualtats: el conjunt actiu

Amb restriccions del tipus gi(x)≤0g_i(x)\leq 0, una restricció és activa al punt si s'hi compleix amb igualtat. Les inactives no limiten les variacions petites i, a efectes locals, és com si no hi fossin.

Les condicions que generalitzen Lagrange són les KKT, amb multiplicadors no negatius per a les desigualtats i la condició de complementarietat que obliga cada multiplicador a anul·lar-se si la seva restricció és inactiva. Es tracten al tema de convexitat i KKT.

Amb molt poques desigualtats, un mètode elemental és enumerar conjunts actius: per a cada subconjunt de restriccions que suposes actives, resols les equacions i després comproves que totes les desigualtats i tots els signes dels multiplicadors es compleixen. Funciona, però el nombre de combinacions creix exponencialment i no escala.

Guió per a un problema amb igualtats

  1. Existència: el conjunt factible és compacte? La funció és coerciva?
  2. Qualificació: els gradients de les restriccions són independents a tot el conjunt? Anota els punts on no.
  3. Equacions de Lagrange i resolució, vigilant els casos degenerats (què passa si una variable és zero).
  4. Afegeix els punts singulars que has apuntat al pas 2.
  5. Compara valors entre tots els candidats. Si has demostrat existència, el més gran és el màxim i el més petit el mínim.
Exercici resolt: distància mínima a una recta

Minimitza x2+y2x^{2}+y^{2} subjecte a x+y=3x+y = 3.

Equacions. Amb L=x2+y2+λ(x+y−3)\mathcal{L} = x^{2}+y^{2}+\lambda(x+y-3):

2x+λ=0,2y+λ=0,x+y=3.2x+\lambda = 0, \qquad 2y+\lambda = 0, \qquad x+y = 3.

De les dues primeres surt x=yx = y, i amb la restricció x=y=3/2x = y = 3/2. El valor és 9/29/2.

Globalitat, demostrada sense Lagrange. Amb la identitat

x2+y2=(x+y)2+(x−y)22=9+(x−y)22≥92,x^{2}+y^{2} = \frac{(x+y)^{2}+(x-y)^{2}}{2} = \frac{9 + (x-y)^{2}}{2}\geq \frac92,

i la igualtat només quan x=yx = y. Per tant és el mínim global i és únic.

Interpretació. El resultat és la distància al quadrat de l'origen a la recta, i el punt òptim és el peu de la perpendicular. El multiplicador val −3-3, i el seu mòdul és la taxa a la qual canviaria el valor òptim si mogués la constant de la restricció: aquesta lectura dels multiplicadors com a sensibilitats és el que els fa útils en aplicacions.

Exercici resolt: caixa de volum màxim

Una caixa rectangular sense tapa ha de tenir una superfície total de 1212 unitats. Quin és el volum màxim?

Amb base xx per yy i alçada zz, volem maximitzar xyzxyz amb xy+2xz+2yz=12xy + 2xz + 2yz = 12 i totes les variables positives.

Equacions de Lagrange.

yz=λ(y+2z),xz=λ(x+2z),xy=λ(2x+2y).yz = \lambda(y+2z), \qquad xz = \lambda(x+2z), \qquad xy = \lambda(2x+2y).

Restant les dues primeres surt z(y−x)=λ(y−x)z(y-x) = \lambda(y-x), i com que zz i λ\lambda no poden coincidir amb la tercera equació a menys que x=yx = y, es dedueix x=yx = y.

Substituint a la tercera, x2=4λxx^{2} = 4\lambda x, o sigui x=4λx = 4\lambda; i de la primera, z=2λz = 2\lambda, és a dir z=x/2z = x/2.

Amb la restricció: x2+2x2=12x^{2} + 2x^{2} = 12, d'on x=2x = 2, y=2y = 2 i z=1z = 1. El volum és 44.

Per què és el màxim. El conjunt factible amb variables positives no és compacte, però quan qualsevol variable tendeix a zero o a infinit el volum tendeix a zero amb la restricció activa. L'únic candidat interior és, doncs, el màxim.

Ara prova-ho tu. 112 exercicis d'aquest tema, de 3 tipus diferents, per resoldre en línia, de deu en deu: escrius la resposta, te la corregeix a l'instant i, si falles, t'explica com es fa pas a pas. Quan acabis una tanda, en tens deu més, amb dades noves.

Fer els exercicisSimulacre d'examen