Das Verfahren der Lagrange-Multiplikatoren (nach Joseph-Louis Lagrange) ist in der mathematischen Optimierung eine Methode zur Lösung von Optimierungsproblemen mit Nebenbedingungen. Ein Optimierungsproblem mit Nebenbedingungen ist die Aufgabe, ein lokales Extremum einer Funktion in mehreren Veränderlichen mit einer oder mehreren Nebenbedingungen zu finden, wobei die Nebenbedingungen als Nullstellen von Funktionen definiert sind. Diese Methode führt eine neue unbekannte skalare Variable für jede Nebenbedingung ein, einen Lagrange-Multiplikator, und definiert eine Linearkombination, die die Multiplikatoren als Koeffizienten einbindet. Die Lösungen der ursprünglichen Optimierungsaufgabe können dann unter gewissen Voraussetzungen als kritische Punkte dieser sogenannten Lagrange-Funktion bestimmt werden.

Zunächst betrachten wir den zweidimensionalen Fall mit einer Nebenbedingung. Nehmen wir an, wir wollen eine Funktion ${\displaystyle f(x,y)}$ maximieren, wobei die Nebenbedingung ${\displaystyle g(x,y)=0}$ einzuhalten ist; manche Quellen verwenden stattdessen ${\displaystyle g(x,y)=c}$ mit einer Konstante ${\displaystyle c}$. Die Nebenbedingung filtert bestimmte Punkte der ${\displaystyle x}$-${\displaystyle y}$-Ebene heraus, die zusammengenommen Kurven bilden. Für unsere Betrachtung nehmen wir an, die Nebenbedingung sei so geartet, dass sie durch eine einzelne Kurve dargestellt werden kann (siehe nebenstehendes Bild, rote Kurve). Wenn wir uns auf dieser Kurve bewegen, berühren oder schneiden wir Höhenlinien von ${\displaystyle f(x,y)}$. Wir sehen nun, dass wir immer nur dann ein Maximum ${\displaystyle d}$ der Funktion ${\displaystyle f(x,y)}$ erreichen, wenn unsere Bewegung auf der Kurve ${\displaystyle g(x,y)=0}$ tangential zur Höhenlinie ${\displaystyle f(x,y)=d}$ verläuft: Andernfalls könnten wir durch Vorwärts- oder Rückwärtsbewegung auf der Kurve ${\displaystyle g(x,y)=0}$ den Funktionswert von ${\displaystyle f(x,y)}$ noch weiter vergrößern, ohne die Nebenbedingung zu verletzen.

Ein bekanntes Beispiel kann den Wetterkarten mit ihren Höhenlinien für Temperaturen und Druck entnommen werden. Die Extrema unter der Nebenbedingung treten dort auf, wo sich beim Überlagern der Karten Linien berühren. Geometrisch übersetzen wir die Tangentenbedingung, indem wir sagen, dass die Gradienten von ${\displaystyle f}$ und ${\displaystyle g}$ beim Maximum parallele Vektoren sind, wobei der Gradient von ${\displaystyle g}$ nicht verschwinden darf.

Wir suchen also Punkte ${\displaystyle (x,y)}$ mit ${\displaystyle g(x,y)=0}$, an denen ${\displaystyle \nabla _{x,y}g\neq 0}$ und

${\displaystyle \nabla _{x,y}f=-\lambda \nabla _{x,y}g}$.

Dabei wurden die folgenden Abkürzungen bzw. Definitionen für die zugehörigen Gradienten benutzt:

${\displaystyle \nabla _{x,y}f:=\left({\frac {\partial f}{\partial x}},{\frac {\partial f}{\partial y}}\right)}$

und

${\displaystyle \nabla _{x,y}g:=\left({\frac {\partial g}{\partial x}},{\frac {\partial g}{\partial y}}\right)}$

Der konstante Lagrange-Multiplikator ${\displaystyle \lambda }$ wird dabei benötigt, weil die beiden Gradienten zwar parallel sein sollen, aber als Vektoren unterschiedlich lang sein können. Um alle genannten Bedingungen zu einer Gleichung zusammenzufassen, ist es nützlich, die folgende Lagrange-Funktion zu verwenden:

${\displaystyle \Lambda (x,y,\lambda ):=f(x,y)+\lambda \cdot g(x,y)}$

Die Lösung des oben beschriebenen Optimierungsproblems mit einer Nebenbedingung entspricht jetzt einem kritischen Punkt der Lagrange-Funktion, d. h.:

${\displaystyle \nabla _{x,y,\lambda }\Lambda (x,y,\lambda )=0.}$

Die ${\displaystyle x}$- und die ${\displaystyle y}$-Komponente dieser Gleichung entsprechen dabei der Forderung nach Parallelität der zwei ursprünglichen Gradienten, die dritte Komponente ${\displaystyle \nabla _{\lambda }\Lambda (x,y,\lambda )=0}$ ist identisch mit ${\displaystyle g(x,y)=0}$.

Punkte, bei denen der Gradient der Nebenbedingung ${\displaystyle g}$ verschwindet, müssen gesondert betrachtet werden, weil das Verfahren der Lagrange-Multiplikatoren über sie keine Aussage treffen kann.

Da im Allgemeinen nicht jeder kritische Punkt der Lagrange-Funktion das ursprüngliche Optimierungsproblem löst, liefert dieses Verfahren nur eine notwendige Bedingung für die Lösung des Optimierungsproblems.

In diesem Beispiel soll die Funktion ${\displaystyle f(x,y)=x+y}$ unter der Nebenbedingung ${\displaystyle x^{2}+y^{2}=1}$ optimiert werden. Die Nebenbedingung entspricht also dem Einheitskreis. Mit Hilfe der Grafik kann das Maximum bei ${\displaystyle ({\sqrt {2}}/2,{\sqrt {2}}/2)}$ bestimmt werden. Das Minimum des Optimierungsproblems liegt bei ${\displaystyle (-{\sqrt {2}}/2,-{\sqrt {2}}/2)}$.

Zunächst überprüfen wir, an welchen Punkten des Einheitskreises der Gradient der Nebenbedingungsfunktion ${\displaystyle g(x,y)=x^{2}+y^{2}-1}$ verschwindet. Wir berechnen also

${\displaystyle \nabla _{x,y}g=(2x,2y)}$

und sehen, dass dies nur im Ursprung gleich ${\displaystyle (0,0)}$ ist. Jedoch liegt dieser Punkt nicht auf dem Einheitskreis, erfüllt also nicht die Nebenbedingung und wird somit nicht in die Liste der kritischen Punkte aufgenommen.

Um die Methode der Lagrange-Multiplikatoren anwenden zu können, sei

${\displaystyle \Lambda (x,y,\lambda )=f(x,y)+\lambda g(x,y)=x+y+\lambda x^{2}+\lambda y^{2}-\lambda }$.

Die Bedingung ${\displaystyle d\Lambda =0}$ ergibt die folgenden drei Gleichungen:

${\displaystyle {\begin{aligned}{\frac {\partial \Lambda }{\partial x}}&=1+2\lambda x&&=0,\qquad {\text{(i)}}\\{\frac {\partial \Lambda }{\partial y}}&=1+2\lambda y&&=0,\qquad {\text{(ii)}}\\{\frac {\partial \Lambda }{\partial \lambda }}&=x^{2}+y^{2}-1&&=0.\qquad {\text{(iii)}}\end{aligned}}}$

Die dritte Gleichung (iii) entspricht dabei wie immer der geforderten Nebenbedingung. Mit ${\displaystyle \lambda \neq 0}$ kann (i) nach ${\displaystyle x}$ aufgelöst werden. Dasselbe macht man für Gleichung (ii) und ${\displaystyle y}$. Man erhält somit ${\displaystyle x=-1/(2\lambda )=y}$. Wird das in (iii) eingesetzt, erhält man ${\displaystyle 2\lambda ^{2}=1}$, also ${\displaystyle \lambda =\pm 1/{\sqrt {2}}=\pm {\sqrt {2}}/2}$. Die kritischen Punkte berechnen sich damit zu ${\displaystyle ({\sqrt {2}}/2,{\sqrt {2}}/2)}$ und ${\displaystyle (-{\sqrt {2}}/2,-{\sqrt {2}}/2)}$. Die zu optimierende Funktion ${\displaystyle f}$ hat an diesen zwei Punkten die Werte ${\displaystyle {\sqrt {2}}}$, bzw. ${\displaystyle -{\sqrt {2}}}$.

Ein Grundstück soll die Form einer Ellipse mit Schwerpunkt im Ursprung und Haupt- und Nebenachse parallel zur ${\displaystyle x}$- und ${\displaystyle y}$-Achse aufweisen. Außerdem soll das Grundstück mit der Fläche

${\displaystyle A=\pi ab}$

so klein wie möglich sein. Die Ellipse soll durch einen gegebenen Punkt ${\displaystyle (x,y)}$ gehen. Diese Nebenbedingung liegt in Form der Gleichung

${\displaystyle {\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{b^{2}}}=1}$

vor, die eine Ellipse mit Zentrum im Ursprung, Hauptachse der Länge ${\displaystyle 2a}$ und Nebenachse der Länge ${\displaystyle 2b}$ beschreibt. Hauptachse und Nebenachse liegen parallel zur ${\displaystyle x}$- und ${\displaystyle y}$-Achse, weshalb die Ellipse durch die Punkte ${\displaystyle (a,0)}$ und ${\displaystyle (0,b)}$ verläuft.

Die allgemeine Lagrange-Funktion mit beliebigen Werten für ${\displaystyle x}$ und ${\displaystyle y}$ lautet

${\displaystyle \Lambda (a,b,\lambda )=\pi ab+\lambda \left({\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{b^{2}}}-1\right)=\pi ab+x^{2}{\frac {\lambda }{a^{2}}}+y^{2}{\frac {\lambda }{b^{2}}}-\lambda }$

mit ${\displaystyle \lambda \neq 0}$. Der Gradient hiervon wird auf Null gesetzt, um die kritischen Punkte zu bestimmen. Hierbei wird ${\displaystyle a>x>0}$ und ${\displaystyle b>y>0}$ vorausgesetzt. Denn für ${\displaystyle a=0}$ oder ${\displaystyle b=0}$ erhielte man eine leere Ellipse.

${\displaystyle {\frac {\partial }{\partial a}}\Lambda (a,b,\lambda )=\pi b-2x^{2}{\frac {\lambda }{a^{3}}}=0\Leftrightarrow \pi a^{3}b-2x^{2}\lambda =0}$

${\displaystyle {\frac {\partial }{\partial b}}\Lambda (a,b,\lambda )=\pi a-2y^{2}{\frac {\lambda }{b^{3}}}=0\Leftrightarrow \pi ab^{3}-2y^{2}\lambda =0}$

${\displaystyle {\frac {\partial }{\partial \lambda }}\Lambda (a,b,\lambda )={\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{b^{2}}}-1=0\Leftrightarrow x^{2}b^{2}+y^{2}a^{2}=a^{2}b^{2}}$

Die erste Gleichung wird zu ${\displaystyle \lambda ={\frac {\pi a^{3}b}{2x^{2}}}}$ umgeformt und in die zweite Gleichung eingesetzt. ${\displaystyle \pi ab^{3}-{\frac {\pi y^{2}a^{3}b}{x^{2}}}=0}$ ergibt ${\displaystyle x^{2}b^{2}=y^{2}a^{2}}$. Setzt man diesen Ausdruck in die dritte Gleichung ein, erhält man durch Rückeinsetzen die Lösungen

${\displaystyle a={\sqrt {2}}x,\ b={\sqrt {2}}y\ {\text{und}}\ \lambda ={\pi ab}={2\pi xy}\neq 0}$.

Der Gradient der Nebenbedingungsfunktion

${\displaystyle \nabla g(a,b)=\left(-{\frac {2x^{2}}{a^{3}}},-{\frac {2y^{2}}{b^{3}}}\right)}$

verschwindet für ${\displaystyle x=y=0}$. Da der Punkt ${\displaystyle (0,0)}$ nicht auf der Ellipse liegt, handelt es sich bei ${\displaystyle (a,b)=({\sqrt {2}}x,{\sqrt {2}}y)}$ um das gesuchte Minimum. Dies muss im Einzelfall grafisch überprüft werden, da die Lagrange-Multiplikatoren nur ein notwendiges Kriterium liefern.

Wir betrachten die Funktion ${\displaystyle f\colon \mathbb {R} _{0}^{+}\times \mathbb {R} _{0}^{+}\to \mathbb {R} }$ mit ${\displaystyle f(x,y)=e^{-(x+y)}}$. Untersucht man die Funktion nun auf Extrema, so kann man mithilfe des hinreichenden Kriteriums für lokale Extremstellen alle Extrema im Inneren des Definitionsbereiches bestimmen. Die Randextrema werden jedoch mithilfe des Lagrange-Multiplikator gefunden. Dabei bildet der Rand des Definitionsbereiches die Nebenbedingung. Hier sind es die beiden positiven Koordinatenachsen und der Ursprung. Wir finden also die Nebenbedingung ${\displaystyle g\colon \mathbb {R} _{0}^{+}\times \mathbb {R} _{0}^{+}\to \mathbb {R} }$ mit ${\displaystyle g(x,y)=xy=0}$.

Wir stellen zunächst die Lagrange-Funktion auf:

${\displaystyle \Lambda (x,y,\lambda )=f(x,y)+\lambda g(x,y)=e^{-(x+y)}+\lambda xy}$

Die Gleichung

${\displaystyle \nabla _{x,y,\lambda }\Lambda =0}$

führt uns auf das Gleichungssystem

${\displaystyle {\begin{aligned}{\frac {\partial \Lambda }{\partial x}}&=-e^{-(x+y)}+\lambda y&&=0,\qquad {\text{(i)}}\\{\frac {\partial \Lambda }{\partial y}}&=-e^{-(x+y)}+\lambda x&&=0,\qquad {\text{(ii)}}\\{\frac {\partial \Lambda }{\partial \lambda }}&=xy&&=0.\qquad {\text{(iii)}}\end{aligned}}}$

Die dritte Gleichung besagt, dass ${\displaystyle x=0}$ oder ${\displaystyle y=0}$. Angenommen es wäre ${\displaystyle x=0}$, dann führt dies – in die zweite Gleichung eingesetzt – auf einen Widerspruch, denn die Gleichung

${\displaystyle -e^{-y}=0}$

hat keine Lösung, da die ${\displaystyle e}$-Funktion keine Nullstellen besitzt. Analog führt man den Fall ${\displaystyle y=0}$ mit der ersten Gleichung auf einen Widerspruch. Der Lagrange-Multiplikator liefert also keine kritischen Punkte.

Jedoch haben wir nicht überprüft, an welchen Stellen der Gradient der Nebenbedingung verschwindet. Es gilt

${\displaystyle \nabla _{x,y}g=(y,x)=(0,0)\Leftrightarrow (x,y)=(0,0)}$

Im Ursprung verschwindet also der Gradient der Nebenbedingung, und dieser liegt auch auf dem Rand des Definitionsbereiches von ${\displaystyle f}$ (er erfüllt die Nebenbedingung). Wie oben beschrieben, müssen diese Punkte auch als Kandidaten für Extrema in Betracht gezogen werden. In der Tat ist ${\displaystyle f(0,0)=e^{0}=1}$ und ${\displaystyle f(x,y)=e^{-(x+y)}<1}$ für alle ${\displaystyle (x,y)\neq (0,0)}$. Der Ursprung ist also das globale Maximum der Funktion.

Wie eine Variante des Beispiels zeigt, sagt das Vorhandensein von kritischen Punkten allerdings nichts über das Vorhandensein von Extrema aus. Würde man in obigem Beispiel die Definitionsbereiche von ${\displaystyle f}$ und ${\displaystyle g}$ durch ${\displaystyle \mathbb {R} \times \mathbb {R} }$ ersetzen, so würde man zwar denselben einzigen kritischen Punkt erhalten, jedoch wäre der Ursprung kein globales (und auch kein lokales) Maximum von ${\displaystyle f}$ (z. B. divergiert die Funktion im 3. Quadranten). In der Tat besäße dieses ${\displaystyle f}$ keine lokalen Maxima oder Minima.

Es sei ${\displaystyle f}$ eine in einer offenen Teilmenge ${\displaystyle U\subseteq \mathbb {R} ^{n}}$ definierte reellwertige Funktion. Wir definieren ${\displaystyle s}$ voneinander unabhängige Nebenbedingungen ${\displaystyle g_{k}(x)=0}$, ${\displaystyle k=1,\ldots ,s}$. D. h., die Gradienten der Nebenbedingungen sind für jeden Punkt ${\displaystyle x\in U}$, mit ${\displaystyle g_{k}(x)=0}$ für alle ${\displaystyle k=1,\ldots ,s}$ linear unabhängig. Insbesondere bedeutet dies, dass keiner der Gradienten verschwindet. Sollten die Gradienten doch an einer Stelle linear abhängig sein, so wird dieser Punkt in die Liste der kritischen Punkte aufgenommen. Nun setzen wir

${\displaystyle \Lambda (x,\lambda )=f(x)+\sum _{k=1}^{s}\lambda _{k}g_{k}(x),}$

wobei ${\displaystyle x=(x_{1},\dots ,x_{n})}$ und ${\displaystyle \lambda =(\lambda _{1},\dots ,\lambda _{s})}$ ist.

Wir schauen uns nun den kritischen Punkt von ${\displaystyle \Lambda }$ an

${\displaystyle {\frac {\partial \Lambda }{\partial x_{i}}}=0,}$

was äquivalent ist zu

${\displaystyle {\frac {\partial f}{\partial x_{i}}}=-\sum _{k=1}^{s}\lambda _{k}{\frac {\partial g_{k}}{\partial x_{i}}}.}$

Wir erhalten ein lineares Gleichungssystem für die unbekannten Multiplikatoren ${\displaystyle \lambda _{1},\dotsc ,\lambda _{s}}$, dessen Koeffizientenmatrix

${\displaystyle \left({\frac {\partial g_{k}}{\partial x_{i}}}\right)_{1\leq i\leq n,1\leq k\leq s}}$

gemäß obiger Voraussetzung maximalen Spaltenrang ${\displaystyle s}$, also trivialen Nullraum besitzt. Falls es lösbar ist, sind die Lagrangemultiplikatoren durch diese Bedingung also eindeutig bestimmt; das Gleichungssystem kann aber auch innerhalb der Nebenbedingungen widersprüchlich sein und keine Lösung besitzen, siehe obiges Beispiel. Im Allgemeinen werden durch die Gleichungen des kritischen Punktes ${\displaystyle {\tfrac {\partial \Lambda }{\partial x_{i}}}=0}$ simultan diejenigen Wertekombinationen ${\displaystyle x=(x_{1},\dots ,x_{n})}$ und ${\displaystyle \lambda =(\lambda _{1},\dots ,\lambda _{s})}$ bestimmt, für die das lineare Gleichungssystem Lösungen besitzt. Das Vorliegen eines kritischen Punktes ist eine notwendige Bedingung dafür, dass ${\displaystyle f}$ auf der Menge der Punkte, welche die Nebenbedingungen erfüllen, ein Extremum hat. D. h., auch hier müssen die Extrema aus der Liste der kritischen Punkte mit anderen Mitteln herausgefiltert werden.

Man beachte, dass es deshalb insbesondere falsch ist, davon zu sprechen, die „Lagrange-Funktion zu maximieren“. Die Lagrange-Funktion ist unbeschränkt und besitzt deshalb keine globalen Extrema und kann somit nicht maximiert werden. Lediglich die kritischen Stellen der Lagrange-Funktion geben Punkte an, an denen die Zielfunktion bezüglich der Nebenbedingungen möglicherweise ein Maximum annimmt.

Dieses Verfahren liefert nur eine notwendige Bedingung für Extremstellen. Um die Extremstellen nachzuweisen und ihre Art zu bestimmen, gibt es verschiedene Kriterien. Generell wird die geränderte Hesse-Matrix gebildet und deren Determinante bzw. bestimmte Unterdeterminanten berechnet. Dieser Ansatz führt aber nicht immer zu einer Aussage. Alternativ kann man auch auf eine Visualisierung bzw. geometrische Überlegungen zurückgreifen, um die Art der Extremstelle festzustellen.

Die Bedeutung der Lagrange-Multiplikatoren in der Physik wird bei der Anwendung in der klassischen Mechanik sichtbar. Hierfür wurden sie von Lagrange um das Jahr 1777 auch eingeführt. Die Bewegungsgleichungen der klassischen Mechanik lassen sich im Lagrange-Formalismus mit Hilfe der Euler-Lagrange-Gleichung aus der Bedingung gewinnen, dass die Wirkung – bei Variation der Koordinaten und ihrer Zeitableitungen unabhängig voneinander – ein Extremum annimmt. Eine physikalische Zwangsbedingung, die die Bewegung einschränkt, erscheint als Nebenbedingung des Extremums. Der Lagrange-Multiplikator, mit dem die Zwangsbedingung in die Lagrange-Funktion eingefügt wird, steht im engen Zusammenhang zu der physikalischen Zwangskraft, mit der das durch die Bewegungsgleichung beschriebene Objekt zur Einhaltung der Zwangsbedingung gebracht wird. Das folgende Beispiel einer freien Punktmasse ${\displaystyle m}$, die sich in zwei Dimensionen auf einer Bahn mit konstantem Radius ${\displaystyle R}$ bewegt, macht dieses klar:

Lagrange-Funktion (kinetische Energie in Polarkoordinaten):

${\displaystyle L={\frac {1}{2}}m({\dot {r}}^{2}+r^{2}{\dot {\varphi }}^{2}),}$

Zwangsbedingung:

${\displaystyle g=r-R=0,}$

neue Lagrange-Funktion:

${\displaystyle L^{\prime }=L+\lambda g,}$

Euler-Lagrange-Gleichung (hier nur für die radiale Koordinate formuliert, da die Zwangsbedingung von dieser abhängt; die Winkelkoordinate ergibt die Drehimpulserhaltung für diese Bewegung):

${\displaystyle {\frac {d}{dt}}{\frac {\partial L^{\prime }}{\partial {\dot {r}}}}-{\frac {\partial L^{\prime }}{\partial r}}=0}$

${\displaystyle m{\ddot {r}}-mr{\dot {\varphi }}^{2}-\lambda =0;}$

mit ${\displaystyle {\ddot {r}}=0}$ und ${\displaystyle r=R}$ sowie ${\displaystyle {\dot {\varphi }}=\omega }$ (Winkelgeschwindigkeit) folgt

${\displaystyle \lambda =-mR\omega ^{2}.}$

Das entspricht der in Polarkoordinaten formulierten Zentripetalkraft, die die Punktmasse zur Bewegung auf eine Kreisbahn zwingt.

Die Karush-Kuhn-Tucker-Bedingungen und die Fritz-John-Bedingungen sind eine Verallgemeinerung der Lagrange-Multiplikatoren für Nebenbedingungen, die auch durch Ungleichungen beschrieben werden. Beide spielen eine wichtige Rolle in der nichtlinearen Optimierung. Für konvexe Optimierungsprobleme, bei denen die Funktionen nicht stetig differenzierbar sind, gibt es außerdem die Sattelpunktkriterien der Lagrange-Funktion.

• Otto Forster: Analysis 2. Vieweg+Teubner, Wiesbaden 2008, ISBN 978-3-8348-0575-1. S. 110ff
• Michael Sauer: Operations Research kompakt 1. Auflage, Oldenbourg, München 2009, ISBN 978-3-486-59082-1.
• Heinrich Rommelfanger: Mathematik für Wirtschaftswissenschaftler II Band 2 (2. Auflage, 1992), BI Wissenschaftsverlag, ISBN 978-3-86025-981-8. S. 238ff

• Konzeptuelle Einleitung (englisch)
• Studienarbeit Lagrangeoptimierung (inklusive Beispielrechnung und Literaturverweisen; PDF; 474 kB)