Lagrange-Dualität

Die Lagrange-Dualität ist eine wichtige Dualität in der mathematischen Optimierung, die sowohl Optimalitätskriterien mittels der Karush-Kuhn-Tucker-Bedingungen oder der Lagrange-Multiplikatoren liefert als auch äquivalente Umformulierungen von Optimierungsproblemen möglich macht. Ziel ist es das ursprüngliche (primale) Problem in ein duales Problem zu überführen.

Lagrange-Funktion, duales Problem

Gegeben sei folgendes Optimierungsproblem:

${\displaystyle {\begin{aligned}{\text{Minimiere }}&f(x)&\\{\text{unter den Nebenbedingungen }}&g_{i}(x)\leq 0&i=1,\dotsc ,p\\&h_{j}(x)=0&j=1,\dotsc ,q\\&x\in X&\end{aligned}}}$

Dabei kann die abstrakte Restriktion ${\displaystyle x\in X}$ beispielsweise Forderungen wie ${\displaystyle x\in \mathbb {Z} ^{n}}$ (Ganzzahligkeit) oder ${\displaystyle x\in {\mathcal {K}}}$ für einen Kegel ${\displaystyle {\mathcal {K}}}$ aufnehmen. Die auftretenden Funktionen müssen weder konvex noch differenzierbar sein.

Die Funktion

${\displaystyle L(x,\lambda ,\mu ):=f(x)+\sum _{i=1}^{p}\lambda _{i}g_{i}(x)+\sum _{j=1}^{q}\mu _{j}h_{j}(x)}$

heißt dann die zu dem obigen Optimierungsproblem gehörende Lagrange-Funktion. Gelegentlich werden die Funktionen ${\displaystyle g_{i},h_{j}}$ sowie die Skalare ${\displaystyle \lambda _{i},\mu _{j}}$ zu Vektoren ${\displaystyle g(x)=(g_{1}(x),\dots ,g_{p}(x))^{T},h(x)=(h_{1}(x),\dots ,h_{q}(x))^{T}}$ und ${\displaystyle \lambda =(\lambda _{1},\dots ,\lambda _{p})^{T},\mu =(\mu _{1},\dots ,\mu _{q})^{T}}$ zusammengefasst. Die Lagrange-Funktion vereinfacht sich dann zu

${\displaystyle L(x,\lambda ,\mu ):=f(x)+\lambda ^{T}g(x)+\mu ^{T}h(x).}$

${\displaystyle (\lambda ,\mu )}$ werden duale Variablen oder Lagrange-Multiplikatoren genannt. Die Funktion

${\displaystyle q(\lambda ,\mu )=\inf _{x\in X}L(x,\lambda ,\mu )}$

heißt dann die duale Funktion zu dem obigen Optimierungsproblem. Das Optimierungsproblem

${\displaystyle {\begin{aligned}{\text{Maximiere }}&q(\lambda ,\mu )\\{\text{unter den Nebenbedingungen }}&\lambda \geq 0\end{aligned}}}$

heißt das duale Problem des obigen Optimierungsproblems. Dabei ist ${\displaystyle \lambda \geq 0}$ komponentenweise zu verstehen. Das ursprüngliche Problem wird dann auch als primales Problem bezeichnet. Im Allgemeinen muss die duale Funktion nicht reellwertig sein, sie kann durchaus auch den Wert ${\displaystyle -\infty }$ annehmen. Man definiert dann

${\displaystyle \operatorname {dom} _{q}:=\{(\lambda ,\mu )\in \mathbb {R} ^{p}\times \mathbb {R} ^{q}\,|\,\lambda \geq 0,q(\lambda ,\mu )>-\infty \}}$

als den wesentlichen Definitionsbereich der dualen Funktion. Die dualen Variablen ${\displaystyle (\lambda ,\mu )}$ werden dual zulässig genannt, wenn sie zulässig bezüglich des dualen Problems sind, das heißt, wenn ${\displaystyle \lambda \geq 0}$ gilt.

Beispiel

Betrachtet man als Beispiel ein lineares Optimierungsproblem in Ungleichungsform:

${\displaystyle {\begin{aligned}{\text{Minimiere }}&f(x)=c^{T}x\\{\text{unter den Nebenbedingungen }}&Ax-b\leq 0\end{aligned}}}$

Dabei ist ${\displaystyle x,c\in \mathbb {R} ^{n}\,,b\in \mathbb {R} ^{m}}$ und ${\displaystyle A\in \mathbb {R} ^{m\times n}}$. Der Vollständigkeit halber setzt man ${\displaystyle X=\mathbb {R} ^{n}}$, was in diesem Fall keine Einschränkung bedeutet. Die Lagrange-Funktion ist dann

${\displaystyle L(x,\lambda )=c^{T}x+\lambda ^{T}(Ax-b)=(c^{T}+\lambda ^{T}A)x-\lambda ^{T}b}$.

Ist ${\displaystyle c^{T}=-\lambda ^{T}A}$, so ist ${\displaystyle L(x,\lambda )}$ unabhängig von ${\displaystyle x}$ und damit ist ${\displaystyle \inf _{x\in X}L(x,\lambda )=-\lambda ^{T}b}$. Ist aber ${\displaystyle c^{T}\neq -\lambda ^{T}A}$, so ist ${\displaystyle L(x,\lambda )}$ in eine Richtung nach unten unbeschränkt und demnach gilt ${\displaystyle \inf _{x\in X}L(x,\lambda )=-\infty }$. Somit lautet die duale Funktion:

${\displaystyle q(\lambda )={\begin{cases}-\lambda ^{T}b&{\text{ falls }}c^{T}+\lambda ^{T}A=0\\-\infty &{\text{ sonst }}\end{cases}}}$

Schreibt man nun die Fallunterscheidung aus der dualen Funktion als Nebenbedingung in das duale Problem, so erhält man:

${\displaystyle {\begin{aligned}{\text{Maximiere }}&-\lambda ^{T}b\\{\text{unter den Nebenbedingungen }}&A^{T}\lambda +c=0\\&\lambda \geq 0\end{aligned}}}$

Dies ist genau ein lineares Optimierungsproblem in Standardform.

Eigenschaften des dualen Problems

• Die Definitionsmenge ${\displaystyle \operatorname {dom} _{q}}$ der dualen Funktion ist konvex.
• Die duale Funktion ist konkav. Für fixiertes ${\displaystyle {\tilde {x}}}$ ist ${\displaystyle L({\tilde {x}},\lambda ,\mu )}$ eine affine Funktion und damit ist ${\displaystyle q(\lambda ,\mu )}$ als punktweises Infimum von affinen Funktionen konkav. Somit ist das duale Problem immer ein konvexes Optimierungsproblem, unabhängig von der Struktur des primalen Problems.
• Daher sind konvexe Optimierungsprobleme eine Klasse von Problemen, deren duales Problem wieder in derselben Problemklasse liegt. Weitere solche Klassen sind die lineare Optimierungsprobleme (siehe Beispiel), die konischen Programme sowie die semidefiniten Programme und die SOCPs.

Schwache Dualität

Sei ${\displaystyle R_{P}}$ die Restriktionsmenge des primalen Problems und ${\displaystyle \operatorname {dom} _{q}}$ die Definitionsmenge des dualen Problems. Dann gilt für alle ${\displaystyle x\in R_{P}}$ und ${\displaystyle (\lambda ,\mu )\in \operatorname {dom} _{q}}$

${\displaystyle q(\lambda ,\mu )\leq f(x).}$

Der Wert ${\displaystyle f(x)-q(\lambda ,\mu )}$ heißt dann die Dualitätslücke (des zulässigen Punktes ${\displaystyle (x,\lambda ,\mu )}$). Die duale Funktion ist also immer eine untere Schranke für die Zielfunktion des primalen Problems. Somit lässt sich auch das duale Problem motivieren: Es stellt die Frage nach der größten unteren Schranke für das primale Problem, die noch zulässig ist.

Ist ${\displaystyle P=f(R_{P})}$ die Wertemenge des primalen Problems und ${\displaystyle D=q(\operatorname {dom} _{q})}$ die des dualen Problems, so gilt

${\displaystyle \sup(D)\leq \inf(P)}$.

Der Wert der dualen Optimallösung ist also stets kleiner als der Wert der primalen Optimallösung. Diese Aussage wird auch schwache Dualität genannt. Der Wert ${\displaystyle \inf(P)-\sup(D)}$ heißt dann die optimale Dualitätslücke.

Diese Aussagen folgen direkt aus

${\displaystyle {\begin{aligned}q(\lambda ,\mu )&=\inf _{x\in X}L(x,\lambda ,\mu )\leq L({\tilde {x}},\lambda ,\mu )\\&=f({\tilde {x}})+\sum _{i=1}^{p}\lambda _{i}g_{i}({\tilde {x}})+\sum _{j=1}^{q}\mu _{j}h_{j}({\tilde {x}})\leq f({\tilde {x}})\quad {\text{für alle }}(\lambda ,\mu )\in \operatorname {dom} _{q}{\text{ und }}{\tilde {x}}\in R_{P}.\end{aligned}}}$

Dabei folgt die letzte Ungleichung, da ${\displaystyle g_{i}({\tilde {x}})\leq 0}$ und ${\displaystyle h_{j}({\tilde {x}})=0}$ (Zulässigkeit von ${\displaystyle {\tilde {x}}}$) und ${\displaystyle \lambda _{i}\geq 0}$ (Zulässigkeit von ${\displaystyle \lambda }$) und damit ${\displaystyle \lambda _{i}g_{i}({\tilde {x}})\leq 0}$ und ${\displaystyle \mu _{i}h_{i}({\tilde {x}})=0}$. Da die Ungleichung für beliebige ${\displaystyle (\lambda ,\mu )\in \operatorname {dom_{q}} }$ und ${\displaystyle {\tilde {x}}\in R_{P}}$ gilt, folgen die beiden obigen Aussagen.

Starke Dualität

Unter gewissen Umständen stimmen der Optimalwert des primalen Problems und der des dualen Problems überein, die optimale Dualitätslücke ist also Null. Es gilt dann

${\displaystyle \sup(D)=\inf(P)}$.

Dieser Fall wird dann starke Dualität genannt. Gilt starke Dualität und ist der Optimalwert endlich und wird in ${\displaystyle x^{*}}$ bzw. ${\displaystyle (\lambda ^{*},\mu ^{*})}$ angenommen, so gilt:

${\displaystyle \sup(D)=\inf(P)=f(x^{*})=q(\lambda ^{*},\mu ^{*})=L(x^{*},\lambda ^{*},\mu ^{*})}$

Im Allgemeinen gilt starke Dualität nicht, sondern es müssen noch weitere Regularitätsvoraussetzungen (im Englischen constraint qualifications) an das Problem erfüllt sein. Eine der wichtigsten Voraussetzungen für konvexe Optimierungsprobleme und fast-konvexe Funktionen, unter der starke Dualität gilt, ist zum Beispiel die Slater-Bedingung.

Komplementärer Schlupf

Gilt die starke Dualität, sind ${\displaystyle x^{*}}$ und ${\displaystyle (\lambda ^{*},\mu ^{*})}$ primal bzw. dual optimal und ist der Optimalwert endlich, so gilt die complementary slackness, auch komplementärer Schlupf genannt:

${\displaystyle \lambda _{i}g_{i}(x^{*})=0{\text{ für alle }}i=1,\dotsc ,p}$

Ist der ${\displaystyle i}$-te Lagrange-Multiplikator (die ${\displaystyle i}$-te Ungleichungsrestriktion) ungleich null, so muss die ${\displaystyle i}$-te Ungleichungsrestriktion (der ${\displaystyle i}$-te Lagrange-Multiplikator) gleich null sein:

${\displaystyle {\begin{aligned}\lambda _{i}>0\implies g_{i}(x^{*})=0\\g_{i}(x^{*})<0\implies \lambda _{i}=0\end{aligned}}}$

Dies folgt aus ${\displaystyle \lambda _{i}\geq 0}$ und ${\displaystyle g_{i}(x^{*})\leq 0}$. Es muss also stets mindestens einer der beiden Faktoren null sein.

Sattelpunkte

Ein Punkt ${\displaystyle (x^{*},\lambda ^{*},\mu ^{*})}$ heißt ein Sattelpunkt der Lagrange-Funktion, wenn für alle ${\displaystyle (x,\lambda ,\mu )}$ mit ${\displaystyle \lambda \geq 0}$

${\displaystyle L(x^{*},\lambda ,\mu )\leq L(x^{*},\lambda ^{*},\mu ^{*})\leq L(x,\lambda ^{*},\mu ^{*})}$

gilt. Äquivalent dazu ist

${\displaystyle \inf _{x\in X}\sup _{(\lambda ,\mu )\in dom_{q}}L(x,\lambda ,\mu )=L(x^{*},\lambda ^{*},\mu ^{*})=\sup _{(\lambda ,\mu )\in dom_{q}}\inf _{x\in X}L(x,\lambda ,\mu ).}$

Dies bedeutet, dass ${\displaystyle (\lambda ^{*},\mu ^{*})}$ ein Maximum der Lagrange-Funktion für fixiertes ${\displaystyle x^{*}}$ und ${\displaystyle x^{*}}$ ein Minimum der Lagrange-Funktion für fixiertes ${\displaystyle (\lambda ^{*},\mu ^{*})}$ ist.

Es lässt sich zeigen, dass ${\displaystyle (x^{*},\lambda ^{*},\mu ^{*})}$ genau dann ein Sattelpunkt der Lagrange-Funktion ist, wenn ${\displaystyle x^{*}}$ primal optimal ist, ${\displaystyle (\lambda ^{*},\mu ^{*})}$ dual optimal ist und starke Dualität gilt.

Sattelpunkte spielen eine wichtige Rolle bei der Bestimmung von Optimalpunkten und schlagen eine Verbindung zu den Karush-Kuhn-Tucker-Bedingungen und den Lagrange-Multiplikatoren. Sind zum Beispiel alle beteiligten Funktionen stetig differenzierbar, so lässt sich aus dem Sattelpunktkriterium ableiten, dass an einem Optimalpunkt ${\displaystyle x^{*}}$

${\displaystyle \nabla _{x}L(x^{*},\lambda ^{*},\mu ^{*})=0}$

gelten muss, da ${\displaystyle x^{*}}$ nach der Sattelpunktcharakterisierung ${\displaystyle L(x,\lambda ^{*},\mu ^{*})}$ minimiert. Diese Forderung findet sich zum Beispiel in den Karush-Kuhn-Tucker-Bedingungen zur Bestimmung eines Optimalpunktes und als notwendiges Optimalitätskriterium wieder.

Allgemeinere Formulierungen

Formulierung für Hilberträume

Betrachtet man ein Optimierungsproblem mit verallgemeinerten Ungleichungen zwischen reellen Hilberträumen (also reellen vollständigen Vektorräumen, die mit einem Skalarprodukt versehen sind), so ist dies meist von folgender Form:

${\displaystyle {\begin{aligned}{\text{Minimiere }}&f(x)&\\{\text{unter den Nebenbedingungen }}&g_{i}(x)\preccurlyeq _{K_{i}}0&i=1,\dotsc ,p\\&h_{j}(x)=0&j=1,\dotsc ,q\\&x\in X&\end{aligned}}}$

Hierbei sind ${\displaystyle f:V_{0}\mapsto \mathbb {R} }$ die Zielfunktion, ${\displaystyle g_{i}:V_{0}\mapsto V_{i}}$ die Ungleichungsrestriktionsfunktionen und ${\displaystyle h_{j}:V_{0}\mapsto V_{j}}$ die Gleichheitsrestriktionsfunktionen. Die ${\displaystyle V_{i}}$ seien mit dem echten Kegel ${\displaystyle K_{i}}$ ausgestattet, der die verallgemeinerte Ungleichung ${\displaystyle \preccurlyeq _{K_{i}}}$ induziert. Das zum Vektorraum ${\displaystyle V_{s}}$ gehörende Skalarprodukt sei mit ${\displaystyle \langle \cdot ;\cdot \rangle _{s}}$ bezeichnet. Man setzt dann ${\displaystyle \lambda :=(\lambda _{1},\dots ,\lambda _{p})^{T},\mu :=(\mu _{1},\dots ,\mu _{q})^{T}}$. Dabei gilt ${\displaystyle \lambda _{i}\in V_{i}}$ und ${\displaystyle \mu _{j}\in V_{j}}$. Dann ist die Lagrange-Funktion von der Form

${\displaystyle L(x,\lambda ,\mu ):=f(x)+\sum _{i=1}^{p}\langle \lambda _{i};g_{i}(x)\rangle _{i}+\sum _{j=1}^{q}\langle \mu _{j};h_{j}(x)\rangle _{j}}$

und die duale Funktion lautet

${\displaystyle q(\lambda ,\mu )=\inf _{x\in X}L(x,\lambda ,\mu ).}$

Das duale Problem lautet dann:

${\displaystyle {\begin{aligned}{\text{Maximiere }}&q(\lambda ,\mu )&\\{\text{unter den Nebenbedingungen }}&\lambda _{i}\succcurlyeq _{K^{D}}&i=1,\dotsc ,p\end{aligned}}}$,

Dabei ist ${\displaystyle K^{D}}$ der duale Kegel von ${\displaystyle K}$.

Alternative Formulierungen fassen alle Ungleichungsrestriktionen und Gleichheitsrestriktionen zusammen:

${\displaystyle g(x)=(g_{1}(x),\dots ,g_{p}(x))^{T}\,,\,h(x)=(h_{1}(x),\dots ,h_{j}(x))^{T},K:=K_{1}\times \dots \times K_{p}}$

Dies führt zu einer kompakteren Schreibweise, die ohne Summenzeichen und Indizes auskommt und daher für theoretische Betrachtungen bevorzugt wird. Alternativ wird auch die Forderung nach einem echten Kegel abgeschwächt, stattdessen definiert man die Ungleichung nur durch einen Ordnungskegel, der zumindest konvex sein muss. Manchmal wir auch komplett auf Gleichheitsnebenbedingungen verzichtet, man modelliert diese dann stattdessen als Ordnungskegel mit leerem Inneren. Statt ${\displaystyle g(x)\in -K\,,\,h(x)=0}$ zu fordern, definiert man einen Kegel ${\displaystyle K'=K\times \{0\}}$. Dann gilt ${\displaystyle (g(x),h(x))\in -K'}$ genau dann, wenn ${\displaystyle g(x)\in -K\,,\,h(x)=0}$. Für alle drei Varianten sind die Lagrange-Funktion und das duale Problem in der untenstehenden Tabelle angegeben. Die duale Funktion ist stets von der Form ${\displaystyle q(\lambda ,\mu )=\inf _{x\in X}L(x,\lambda ,\mu )}$ bzw., wenn nur eine duale Variable verwendet wird, von der Form ${\displaystyle q(\lambda )=\inf _{x\in X}L(x,\lambda )}$.

Dabei ist ${\displaystyle f:V_{0}\mapsto \mathbb {R} }$ die Zielfunktion, ${\displaystyle g:V_{0}\mapsto V_{1}}$, wobei auf ${\displaystyle V_{1}}$ ein Kegel ${\displaystyle K}$ ausgezeichnet ist, der im ersten Fall ein echter Kegel ist, im zweiten und dritten Fall ein konvexer Kegel ist. Es ist ${\displaystyle h:V_{0}\mapsto V_{2}}$ und ${\displaystyle \lambda \in V_{1},\,\mu \in V_{2}}$.

Formulierung für normierte Räume

Für Probleme mit Abbildungen zwischen reellen normierten Vektorräumen ist zu beachten, dass kein Skalarprodukt definiert ist. Man wählt stattdessen die duale Paarung. Dementsprechend sind die dualen Variablen aus dem Dualraum des Vektorraumes. Ist ein Problem der Form

${\displaystyle {\begin{aligned}{\text{Minimiere }}&f(x)&\\{\text{unter den Nebenbedingungen }}&g(x)\in -K\\&h(x)=0\\&x\in X&\end{aligned}}}$

gegeben, wobei ${\displaystyle f:V_{0}\mapsto \mathbb {R} }$ die Zielfunktion, ${\displaystyle g_{i}:V_{0}\mapsto V_{1}}$ die Ungleichungsrestriktionsfunktionen und ${\displaystyle h_{j}:V_{0}\mapsto V_{2}}$ die Gleichungsrestriktionsfunktion sind. ${\displaystyle K}$ sei ein Ordnungskegel auf ${\displaystyle V_{1}}$ und es seien ${\displaystyle V_{0},V_{1},V_{2}}$ Banachräume. Die Lagrange-Funktion lautet dann

${\displaystyle L(x,u,v)=f(x)+u(g(x))+v(h(x)).}$

Dabei ist ${\displaystyle u}$ aus dem Dualraum ${\displaystyle V_{1}^{*}}$ und ${\displaystyle v}$ aus dem Dualraum ${\displaystyle V_{2}^{*}.}$ Die duale Funktion ist dann wieder

${\displaystyle q(u,v)=\inf _{x\in X}L(x,u,v)}$

und damit lautet das duale Problem:

${\displaystyle {\begin{aligned}{\text{Maximiere }}&q(u,v)&\\{\text{unter den Nebenbedingungen }}&u\in K^{*}\end{aligned}}}$

Dabei ist ${\displaystyle K^{*}}$ der duale Kegel, der in diesem Fall aber eine Teilmenge von ${\displaystyle V_{1}^{*}}$ ist. Formulierungen für alternative Problemstellungen laufen analog zu den Problemen für Abbildungen zwischen Hilberträumen. Die duale Paarung ersetzt jeweils das Skalarprodukt, die dualen Variablen befinden sich im Dualraum.

Schwache und starke Dualität

Die schwache Dualität gilt auch für die beiden allgemeineren Formulierungen. Für den Beweis in der Hilbertraumformulierung nutzt man aus, dass ${\displaystyle \langle a;b\rangle \leq 0}$ ist, wenn ${\displaystyle a\succcurlyeq _{K^{D}}0}$ und ${\displaystyle b\preccurlyeq _{K}}$ gilt (für Ordnungskegel erhält man ${\displaystyle a\in K^{D},\,b\in -K\implies \langle a;b\rangle \leq 0}$). In normierten Räumen gelten analoge Aussagen mit dem Unterschied, dass ${\displaystyle a}$ ein Element des Dualraumes ist: Gilt ${\displaystyle a\in K^{*},b\in -K}$, so folgt ${\displaystyle a(b)\leq 0}$.

Die starke Dualität wird in den allgemeineren Fällen identisch zum gewöhnlichen Fall definiert. Auch im verallgemeinerten Fall existieren Regularitätsvoraussetzungen wie die Slater-Bedingung, die starke Dualität garantieren.

Literatur

• Florian Jarre, Josef Stoer: Optimierung. Springer, Berlin 2004, ISBN 3-540-43575-1.
• Stephan Boyd, Lieven Vandenberghe: Convex Optimization. Cambridge University Press, 2004, ISBN 978-0-521-83378-3 (online).
• C. Geiger, C. Kanzow: Theorie und Numerik restringierter Optimierungsaufgaben. Springer, 2002, ISBN 3-540-42790-2.
• Johannes Jahn: Introduction to the Theory of Nonlinear Optimization. 3. Auflage. Springer, Berlin 2007, ISBN 978-3-540-49378-5.