「双対性理論」の版間の差分

提供: ORWiki
ナビゲーションに移動 検索に移動
3行目: 3行目:
 
 [[双対性理論]] (duality theory)は,非線形計画のみならず線形計画,多目的計画,離散凸解析などの分野で主問題とその双対問題の関係,および集合や関数の双対関係を説明する重要な基礎理論である [1, 2, 3, 4].
 
 [[双対性理論]] (duality theory)は,非線形計画のみならず線形計画,多目的計画,離散凸解析などの分野で主問題とその双対問題の関係,および集合や関数の双対関係を説明する重要な基礎理論である [1, 2, 3, 4].
  
 「双対」 (dual) と「共役」 (conjugate) は元々同義語として用いられ,数学の関数解析の分野では,ノルム空間 $X$ 上の有界線形汎関数の全体を$X$ の双対空間 (dual space) または共役空間 (conjugate space) といい,$X^{*}$ と表して,$x\in{X}$ における $x^{*}\in{X^*}$ の値を$\langle x, x^{*}\rangle$ または $x^{*}(x)$ と書く.$X$ $n$ 次元実ユークリッド空間 $\mathbf{R}^n$ の場合は,$(\mathbf{R}^n)^{*}$ $\mathbf{R}^n$ は同一視でき,$(\mathbf{R}^n)^{**}=\mathbf{R}^n$ となり,$\langle x, x^{*}\rangle$ $x$ $x^{*}$ の内積 $x^{\top}x^{*}$ となる.以下に述べる事柄は,無限次元空間に対しても拡張できるが,ここでは簡単のため $\mathbf{R}^n$ に限定して説明する.
+
 「双対」 (dual) と「共役」 (conjugate) は元々同義語として用いられ,数学の関数解析の分野では,ノルム空間 <math>X\, </math> 上の有界線形汎関数の全体を<math>X\,</math>  の双対空間 (dual space) または共役空間 (conjugate space) といい,<math>X^{*}\, </math> と表して,<math>x\in{X}\, </math> における <math>x^{*}\in{X^*}</math> の値を<math>\langle x, x^{*}\rangle</math> または <math>x^{*}(x)</math> と書く.<math>X\, </math> <math>n\, </math> 次元実ユークリッド空間 <math>\mathbf{R}^n</math> の場合は,<math>(\mathbf{R}^n)^{*}</math> <math>\mathbf{R}^n</math> は同一視でき,<math>(\mathbf{R}^n)^{**}=\mathbf{R}^n</math> となり,<math>\langle x, x^{*}\rangle</math> <math>x\, </math> <math>x^{*}\, </math> の内積 <math>x^{\top}x^{*}</math> となる.以下に述べる事柄は,無限次元空間に対しても拡張できるが,ここでは簡単のため <math>\mathbf{R}^n</math> に限定して説明する.
 空間 $\mathbf{R}^n$ 上で定義された拡張実数値関数 $f: \mathbf{R}^n\to\bar{\mathbf{R}}$ に対して(ただし,$\bar{\mathbf{R}}=\mathbf{R}\cup \{ \infty , -\infty\}$),
+
 空間 <math>\mathbf{R}^n</math> 上で定義された拡張実数値関数 <math>f: \mathbf{R}^n\to\bar{\mathbf{R}}</math> に対して(ただし,<math>\bar{\mathbf{R}}=\mathbf{R}\cup \{ \infty , -\infty\}</math>),
  
  
f^*(\xi):=\sup_{x\in{\mathbf{R}^n}} \{ \, \xi^{\top}x - f(x) \, \}
+
<math>f^*(\xi):=\sup_{x\in{\mathbf{R}^n}} \{ \, \xi^{\top}x - f(x) \, \}</math>
  
  
で定義される関数 $f^*$ $f$ の共役関数 (conjugate function) という.共役関数 $f^*$ に対して,さらにその共役関数 $f^{**}=(f^*)^{*}$ を考えることができるが,$f$ が下半連続な真凸関数のときには,$f^{**}$ $f$ に一致する.$f$ $f^*$ を対応させる写像をルジャンドル--フェンシェル変換 (Legendre-Fenchel transform) と呼ぶ.
+
で定義される関数 <math>f^*\, </math> <math>f\, </math> の共役関数 (conjugate function) という.共役関数 <math>f^*\, </math> に対して,さらにその共役関数 <math>f^{**}=(f^*)^{*}</math> を考えることができるが,<math>f\, </math> が下半連続な真凸関数のときには,<math>f^{**}\, </math> <math>f\, </math> に一致する.<math>f\, </math> <math>f^*\, </math> を対応させる写像をルジャンドル-フェンシェル変換 (Legendre-Fenchel transform) と呼ぶ.
  
 下半連続な真凸関数 $f: \mathbf{R}^n\times{\mathbf{R}^m}\to\bar{\mathbf{R}}$に対して,関数 $\varphi : \mathbf{R}^n \to \bar{\mathbf{R}}$ $\psi : \mathbf{R}^m \to \bar{\mathbf{R}}$ をそれぞれ $\varphi{(x)}:=f(x,0)$ $\psi{(y)}:=-f^{*}(0,y)$ で定義し,次の問題(P)と(D)を主問題 (primal problem) とその双対問題 (dual problem) と呼ぶ [1, 4].
+
 下半連続な真凸関数 <math>f: \mathbf{R}^n\times{\mathbf{R}^m}\to\bar{\mathbf{R}}</math>に対して,関数 <math>\varphi : \mathbf{R}^n \to \bar{\mathbf{R}}</math> <math>\psi : \mathbf{R}^m \to \bar{\mathbf{R}}</math> をそれぞれ <math>\varphi{(x)}:=f(x,0)</math> <math>\psi{(y)}:=-f^{*}(0,y)</math> で定義し,次の問題(P)と(D)を主問題 (primal problem) とその双対問題 (dual problem) と呼ぶ [1, 4].
  
  
 
\begin{center}
 
\begin{center}
 +
\end{center}
 +
 +
<math>
 
\begin{tabular}{cll}
 
\begin{tabular}{cll}
(P) & $\min_{x\in \mathbf{R}^n}$ & \hspace*{-5mm} $\varphi{(x)}$ \\
+
(P) & \min_{x\in \mathbf{R}^n} & \hspace*{-5mm} \varphi{(x)} \\
(D) & $\max_{y\in \mathbf{R}^m}$ & \hspace*{-5mm} $\psi{(y)}$
+
(D) & \max_{y\in \mathbf{R}^m} & \hspace*{-5mm} \psi{(y)}  
 
\end{tabular}
 
\end{tabular}
\end{center}
+
</math>
  
  
26行目: 29行目:
  
  
U:=\{u\in{\mathbf{R}^m}\,|\:\inf_{x\in{\mathbf{R}^n}}f(x,u)<+\infty\}
+
<math>U:=\{u\in{\mathbf{R}^m}\,|\:\inf_{x\in{\mathbf{R}^n}}f(x,u)<+\infty\}
 
\quad
 
\quad
  
V:=\{v\in{\mathbf{R}^n}\,|\:\inf_{y\in{\mathbf{R}^m}}f^{*}(v,y)<+\infty\}
+
V:=\{v\in{\mathbf{R}^n}\,|\:\inf_{y\in{\mathbf{R}^m}}f^{*}(v,y)<+\infty\}</math>
  
  
とおくと,$U$ $V$ は凸集合となる.このとき,以下が成立する.
+
とおくと,<math>U\, </math> <math>V\, </math> は凸集合となる.このとき,以下が成立する.
  
  
 
\medskip
 
\medskip
 
\noindent
 
\noindent
(i) \  $\inf_{x}\varphi{(x)}\ge\sup_{y}\psi{(y)}$ \\
+
(i) \  <math>\inf_{x}\varphi{(x)}\ge\sup_{y}\psi{(y)}</math> \\
(ii) \ $0\in{\mbox{\rm int}\,U}\; または \; 0\in{\mbox{\rm int}\,V}
+
(ii) \ <math>0\in{\mbox{\rm int}\,U}\;</math> または<math> \; 0\in{\mbox{\rm int}\,V}
         \Longrightarrow\inf_{x}\varphi{(x)}=\sup_{y}\psi{(y)}$
+
         \Longrightarrow\inf_{x}\varphi{(x)}=\sup_{y}\psi{(y)}
 +
</math>
  
  
ここで,$\mbox{\rm int}\,U$ $U$ の内部を表す.(i)を弱双対定理 (weak duality theorem),(ii)を[[双対定理]] (duality theorem) と呼び,$\inf_{x}\varphi{(x)}=\sup_{y}\psi{(y)}$ が満たされるとき,主問題(P)と双対問題(D)の間に双対性 (duality) が成立するという.(i)により,$\sup_{y}\psi{(y)}=+\infty$ なら主問題(P)は実行可能解を持たないが,$-\infty<\varphi{(\bar{x})}=\psi{(\bar{y})}<+\infty$ となる $\bar{x}$ $\bar{y}$ が存在すれば,それぞれ(P)と(D)の最適解となり,強い意味の双対性が成立する.一方,$\inf_{x}\varphi{(x)}>\sup_{y}\psi{(y)}$ となるとき,主問題と双対問題の間に[[双対性のギャップ]] (duality gap) が存在するという.
+
ここで,<math>\mbox{\rm int}\,U</math> <math>U\, </math> の内部を表す.(i)を弱双対定理 (weak duality theorem),(ii)を[[双対定理]] (duality theorem) と呼び,<math>\inf_{x}\varphi{(x)}=\sup_{y}\psi{(y)}</math> が満たされるとき,主問題(P)と双対問題(D)の間に双対性 (duality) が成立するという.(i)により,\<math>sup_{y}\psi{(y)}=+\infty</math> なら主問題(P)は実行可能解を持たないが,<math>-\infty<\varphi{(\bar{x})}=\psi{(\bar{y})}<+\infty</math> となる <math>\bar{x}</math> <math>\bar{y}</math> が存在すれば,それぞれ(P)と(D)の最適解となり,強い意味の双対性が成立する.一方,<math>\inf_{x}\varphi{(x)}>\sup_{y}\psi{(y)}</math> となるとき,主問題と双対問題の間に[[双対性のギャップ]] (duality gap) が存在するという.
  
 主問題(P)において,$f(x,u):=c^{\top}x+k(x)+h(b-Ax+u)$(ただし,$k: \mathbf{R}^n\to\bar{\mathbf{R}}$ $h: \mathbf{R}^m\to\bar{\mathbf{R}}$ は下半連続な真凸関数で$A\in{\mathbf{R}^{m\times{n}}}$, $b\in{\mathbf{R}^m}$, $c\in{\mathbf{R}^n}$ )とすると,$f^{*}(v,y)=-b^{\top}y+h^{*}(y)+k^{*}(A^{\top}y-c+v)$となり \cite{[4]},主問題(P)と双対問題(D)はそれぞれ
+
 主問題(P)において,<math>f(x,u):=c^{\top}x+k(x)+h(b-Ax+u)</math>(ただし,<math>k: \mathbf{R}^n\to\bar{\mathbf{R}}</math> <math>h: \mathbf{R}^m\to\bar{\mathbf{R}}</math> は下半連続な真凸関数で<math>A\in{\mathbf{R}^{m\times{n}}}</math>, <math>b\in{\mathbf{R}^m}</math>, <math>c\in{\mathbf{R}^n}</math> )とすると,<math>f^{*}(v,y)=-b^{\top}y+h^{*}(y)+k^{*}(A^{\top}y-c+v)</math>となり [4],主問題(P)と双対問題(D)はそれぞれ
  
  
\begin{eqnarray}
+
<math>\begin{eqnarray}
 
\mbox{\rm min}_{x\in \mathbf{R}^n} &&  
 
\mbox{\rm min}_{x\in \mathbf{R}^n} &&  
 
\hspace*{-5mm} c^{\top}x+k(x)+h(b-Ax)
 
\hspace*{-5mm} c^{\top}x+k(x)+h(b-Ax)
55行目: 59行目:
 
         \label{A-B-03+siki2}
 
         \label{A-B-03+siki2}
 
\end{eqnarray}
 
\end{eqnarray}
 +
</math>
  
  
となる.ここで$b\in\mbox{\rm int}\,(A\mbox{\rm dom}k+\mbox{\rm dom}h)$ または$c\in\mbox{\rm int}\,(A^{\top}\mbox{\rm dom}h^{*}-\mbox{\rm dom}k^{*})$ が成立すれば,(ii)により主問題(\ref{A-B-03+siki1})と双対問題(\ref{A-B-03+siki2})の間に双対性が成立する.(ただし,dom は拡張実数値関数の実効定義域を表す.)これを[[フェンシェルの双対性]] (Fenchel duality) と呼んでいる.通常は,簡略化して $c$ $b$ を零ベクトル,$-A$ を恒等写像として,新たに$f(x)$ を凸関数 $f_1(x):=k(x)$ と凹関数 $f_2(x):=-h(x)$ の差で表し,主問題 $\min_{x}\{f_1(x)-f_2(x)\}$に対して,$\max_{y}\{f_{2}^{*}(y)-f_{1}^{*}(y)\}$ をフェンシェルの双対問題 (Fenchel dual problem) と呼ぶ.ただし,$f_{2}^{*}(y):=\inf_{x\in{\mathbf{R}^n}}\{y^{\top}x-f_{2}(x)\}$.双対性は $\mbox{\rm int}\,(\mbox{\rm dom}f_1)\,\cap\,        \mbox{\rm int}\,(\mbox{\rm dom}f_2)\neq\emptyset$ のとき成立する.また,$k(x):=\sup_{s\le{0}}x^{\top}s$, $h(z):=\sup_{w\ge{0}}z^{\top}w$ とすると,(\ref{A-B-03+siki1})と(\ref{A-B-03+siki2})は線形計画の主問題と双対問題となる \cite{A-B-03-c,[4]}
+
となる.ここで<math>b\in\mbox{\rm int}\,(A\mbox{\rm dom}k+\mbox{\rm dom}h)</math> または<math>c\in\mbox{\rm int}\,(A^{\top}\mbox{\rm dom}h^{*}-\mbox{\rm dom}k^{*})</math> が成立すれば,(ii)により主問題 (1) と双対問題 (2) の間に双対性が成立する.(ただし,dom は拡張実数値関数の実効定義域を表す.)これを[[フェンシェルの双対性]] (Fenchel duality) と呼んでいる.通常は,簡略化して <math>c\, </math> <math>b\, </math> を零ベクトル,<math>-A\, </math> を恒等写像として,新たに<math>f(x)\, </math> を凸関数 <math>f_1(x):=k(x)</math> と凹関数 <math>f_2(x):=-h(x)</math> の差で表し,主問題 <math>\min_{x}\{f_1(x)-f_2(x)\}</math>に対して,<math>\max_{y}\{f_{2}^{*}(y)-f_{1}^{*}(y)\}</math> をフェンシェルの双対問題 (Fenchel dual problem) と呼ぶ.ただし,<math>f_{2}^{*}(y):=\inf_{x\in{\mathbf{R}^n}}\{y^{\top}x-f_{2}(x)\}</math>.双対性は <math>\mbox{\rm int}\,(\mbox{\rm dom}f_1)\,\cap\,        \mbox{\rm int}\,(\mbox{\rm dom}f_2)\neq\emptyset</math> のとき成立する.また,<math>k(x):=\sup_{s\le{0}}x^{\top}s, h(z):=\sup_{w\ge{0}}z^{\top}w</math> とすると,(1) と(2) は線形計画の主問題と双対問題となる [2, 4].
  
  
\begin{center}
+
<math>\begin{center}
 
\begin{tabular}{clcll}
 
\begin{tabular}{clcll}
(P$_{LP}$) & min. & $c^{\top}x$ & s. t. & $x\ge{0}, \ Ax\ge{b}.$ \\
+
(P_{LP}) & min. & c^{\top}x & s. t. & x\ge{0}, \ Ax\ge{b}. \\
(D$_{LP}$) & max. & $b^{\top}y$ & s. t. & $y\ge{0}, \ A^{\top}y\le{c}.$
+
(D_{LP}) & max. & b^{\top}y & s. t. & y\ge{0}, \ A^{\top}y\le{c}.  
 
\end{tabular}
 
\end{tabular}
 
\end{center}
 
\end{center}
 +
</math>
  
  
 
次に,[[ラグランジュ関数]] (Lagrangian function) を
 
次に,[[ラグランジュ関数]] (Lagrangian function) を
  
 +
<math>\begin{equation} \label{A-B-03+siki3}
 +
L(x,y):=\inf_{u\in{\mathbf{R}^m}}\{\, f(x,u)-y^{\top}u \,\}
 +
\end{equation}</math>
  
\begin{equation} \label{A-B-03+siki3}
 
L(x,y):=\inf_{u\in{\mathbf{R}^m}}\{\, f(x,u)-y^{\top}u \,\}
 
\end{equation}
 
  
と定義する.$-L(x,\cdot)=\left(f(x,\cdot)\right)^{*}, f(x,\cdot)=\left(-L(x,\cdot)\right)^{*}$
+
と定義する.<math>-L(x,\cdot)=\left(f(x,\cdot)\right)^{*}, f(x,\cdot)=\left(-L(x,\cdot)\right)^{*}</math>が成立しているので,
が成立しているので,
 
  
\begin{equation} \label{A-B-03+siki4}
+
<math>\begin{equation} \label{A-B-03+siki4}
 
\varphi{(x)}=\sup_{y}L(x,y), \quad\quad \psi{(y)}=\inf_{x}L(x,y)
 
\varphi{(x)}=\sup_{y}L(x,y), \quad\quad \psi{(y)}=\inf_{x}L(x,y)
\end{equation}
+
\end{equation}</math>
 +
 
  
となる.通常,$\inf_{x}L(x,\bar{y})=L(\bar{x},\bar{y})=\sup_{y}L(\bar{x},y)$ すなわち,すべての $x\in{\mathbf{R}^n}$ $y\in{\mathbf{R}^m}$ に対して$L(x,\bar{y})\ge{L(\bar{x},\bar{y})}\ge{L(\bar{x},y)}$が成り立つとき,$(\bar{x},\bar{y})$ を関数 $L$ $\mathbf{R}^{n}\times{\mathbf{R}^m}$ 上での鞍点 (saddle point) と呼ぶ.(4) により,
+
となる.通常,<math>\inf_{x}L(x,\bar{y})=L(\bar{x},\bar{y})=\sup_{y}L(\bar{x},y)</math> すなわち,すべての <math>x\in{\mathbf{R}^n}</math> <math>y\in{\mathbf{R}^m}</math> に対して<math>L(x,\bar{y})\ge{L(\bar{x},\bar{y})}\ge{L(\bar{x},y)}</math>が成り立つとき,<math>(\bar{x},\bar{y})</math> を関数 <math>L\, </math> <math>\mathbf{R}^{n}\times{\mathbf{R}^m}</math> 上での鞍点 (saddle point) と呼ぶ.(4) により,
  
  
\medskip
+
<math>\medskip
 
\noindent
 
\noindent
(iii) \ $\inf_{x}\varphi{(x)}=\inf_{x}\,[\,\sup_{y}L(x,y)\,]
+
(iii) \ \inf_{x}\varphi{(x)}=\inf_{x}\,[\,\sup_{y}L(x,y)\,]
         \ge\sup_{y}\,[\,\inf_{x}L(x,y)\,]=\sup_{y}\psi{(y)}$ \\
+
         \ge\sup_{y}\,[\,\inf_{x}L(x,y)\,]=\sup_{y}\psi{(y)} \\
(iv) \ $\varphi{(\bar{x})}=\inf_{x}\varphi{(x)}
+
(iv) \ \varphi{(\bar{x})}=\inf_{x}\varphi{(x)}
         =\sup_{y}\psi{(y)}=\psi{(\bar{y})}$
+
         =\sup_{y}\psi{(y)}=\psi{(\bar{y})}
         $\Longleftrightarrow$ $(\bar{x},\bar{y})$ $L$ の鞍点
+
         \Longleftrightarrow (\bar{x},\bar{y}) が L の鞍点
 
         \\ \hspace*{1cm}
 
         \\ \hspace*{1cm}
       $\Longleftrightarrow$ $\min_{x}\sup_{y}L(x,y)=\max_{y}\inf_{x}L(x,y)$
+
       \Longleftrightarrow \min_{x}\sup_{y}L(x,y)=\max_{y}\inf_{x}L(x,y)
         $\Longleftrightarrow$
+
         \Longleftrightarrow  
         $\varphi{(\bar{x})}=L(\bar{x},\bar{y})=\psi{(\bar{y})}$
+
         \varphi{(\bar{x})}=L(\bar{x},\bar{y})=\psi{(\bar{y})}  
  
 
\medskip
 
\medskip
 
\noindent
 
\noindent
 +
</math>
  
  
が成立する.(iv)を[[鞍点定理]] (saddle point theorem) と呼ぶ.非線形計画問題
+
が成立する.(iv) を[[鞍点定理]] (saddle point theorem) と呼ぶ.非線形計画問題
  
\mbox{(NLP)} \; \;  
+
<math>\mbox{(NLP)} \; \;  
 
\mbox{min.} \ f_{0}(x) \quad \mbox{s. t.} \; \;  
 
\mbox{min.} \ f_{0}(x) \quad \mbox{s. t.} \; \;  
 
g_{i}(x)\le{0} \ (i=1,\ldots, k), \
 
g_{i}(x)\le{0} \ (i=1,\ldots, k), \
h_{j}(x)=0 \ (j=1,\ldots,l),  
+
h_{j}(x)=0 \ (j=1,\ldots,l),</math>
 +
 
  
(ただし,$f_0,g_i,h_j$ $\mathbf{R}^n$ で定義された実数値関数,$k+l=m$)に対して,
+
(ただし,<math>f_0,g_i,h_j</math> <math>\mathbf{R}^n</math> で定義された実数値関数,<math>k+l=m</math>) に対して,
  
  
\begin{center}
+
<math>\begin{center}
 
\begin{tabular}{r@{}l}
 
\begin{tabular}{r@{}l}
$F(x)\,$&$:=(g_{1}(x),\ldots,g_{k}(x),h_{1}(x),\ldots,h_{l}(x))^{\top}$\\
+
F(x)\,&:=(g_{1}(x),\ldots,g_{k}(x),h_{1}(x),\ldots,h_{l}(x))^{\top}\\
$\theta(w)\,$&$:=\sup_{\lambda,\mu}\{\lambda^{\top}w_{1}+\mu^{\top}w_{2}\,|\,
+
\theta(w)\,&:=\sup_{\lambda,\mu}\{\lambda^{\top}w_{1}+\mu^{\top}w_{2}\,|\,
         (\lambda,\mu)\in{\mathbf{R}^{k}_{+}\times{\mathbf{R}^{l}}},w=(w_1,w_2)^{\top}\}$\\
+
         (\lambda,\mu)\in{\mathbf{R}^{k}_{+}\times{\mathbf{R}^{l}}},w=(w_1,w_2)^{\top}\}\\
$f(x,u)\,$&$:=f_{0}(x)+\theta{(F(x)+u)}$
+
f(x,u)\,&:=f_{0}(x)+\theta{(F(x)+u)}
 
\end{tabular}
 
\end{tabular}
 
\end{center}
 
\end{center}
 +
</math>
  
  
とおくと,(\ref{A-B-03+siki3})により$L(x,y)=f_{0}(x)+y^{T}F(x)-\theta^{*}(y)$ となり,$y=(\lambda,\mu)=(\lambda_{1},\ldots,\lambda_{k},\mu_{1},\ldots,\mu_{l}) \in{\mathbf{R}^{k}_{+}\times{\mathbf{R}^{l}}}$に対する問題(NLP)のラグランジュ関数は
+
とおくと,(3) により<math>L(x,y)=f_{0}(x)+y^{T}F(x)-\theta^{*}(y)</math> となり,<math>y=(\lambda,\mu)=(\lambda_{1},\ldots,\lambda_{k},\mu_{1},\ldots,\mu_{l}) \in{\mathbf{R}^{k}_{+}\times{\mathbf{R}^{l}}}</math>に対する問題(NLP)のラグランジュ関数は
  
  
\begin{equation} \label{A-B-03+siki5}
+
<math>\begin{equation}
 
L(x,\lambda,\mu)=f_{0}(x)+\sum_{i=1}^{k}\lambda_{i}g_{i}(x)
 
L(x,\lambda,\mu)=f_{0}(x)+\sum_{i=1}^{k}\lambda_{i}g_{i}(x)
 
                 +\sum_{j=1}^{l}\mu_{j}h_{j}(x)
 
                 +\sum_{j=1}^{l}\mu_{j}h_{j}(x)
\end{equation}
+
\end{equation}</math>
  
  
となる.この $(\lambda,\mu)$ をラグランジュ乗数 (Lagrange multipliers) と呼ぶ.このとき,主問題と双対問題は
+
となる.この <math>(\lambda,\mu)</math> をラグランジュ乗数 (Lagrange multipliers) と呼ぶ.このとき,主問題と双対問題は
  
  
\begin{center}
+
<math>\begin{center}
 
\begin{tabular}{cllll}
 
\begin{tabular}{cllll}
(P$_{L}$) & min.  
+
(P_{L}) & min.  
         & $\displaystyle \sup_{\lambda\ge{0},\mu} L(x, \lambda, \mu)$
+
         & \displaystyle \sup_{\lambda\ge{0},\mu} L(x, \lambda, \mu)  
         & s. t. & $\displaystyle{x \in {\mathbf{R}^n}}$, \\
+
         & s. t. & \displaystyle{x \in {\mathbf{R}^n}}, \\
(D$_{L}$) & max.
+
(D_{L}) & max.
         & $\displaystyle \inf_{x} L(x,\lambda,\mu)$
+
         & \displaystyle \inf_{x} L(x,\lambda,\mu)
 
         & s. t.  
 
         & s. t.  
         & $\displaystyle{0 \le \lambda \in {\mathbf{R}^{k}}}$, \
+
         & \displaystyle{0 \le \lambda \in {\mathbf{R}^{k}}}, \
           $\displaystyle \mu \in {\mathbf{R}^{l}}$,  
+
           \displaystyle \mu \in {\mathbf{R}^{l}},  
 
\end{tabular}
 
\end{tabular}
 
\end{center}
 
\end{center}
 +
</math>
 +
  
 +
となり,一般に問題(D_{L})をラグランジュの双対問題 (Lagrangian dual problem)と呼ぶ.鞍点定理により,<math>L\, </math> の鞍点 <math>(\bar{x},\bar{\lambda},\bar{\mu})</math> が存在すれば,つまり
  
となり,一般に問題(D$_{L}$)をラグランジュの双対問題 (Lagrangian dual problem)と呼ぶ.鞍点定理により,$L$ の鞍点 $(\bar{x},\bar{\lambda},\bar{\mu})$ が存在すれば,つまり
 
  
 +
<math>\max_{\lambda\ge{0},\mu}L(\bar{x},\lambda,\mu)=L(\bar{x},\bar{\lambda},\bar{\mu})
 +
        =\min_{x}L(x,\bar{\lambda},\bar{\mu})</math>
  
\max_{\lambda\ge{0},\mu}L(\bar{x},\lambda,\mu)=L(\bar{x},\bar{\lambda},\bar{\mu})
 
        =\min_{x}L(x,\bar{\lambda},\bar{\mu})
 
  
 +
が成立すれば,<math>\bar{x}</math> と <math>(\bar{\lambda},\bar{\mu})</math> はそれぞれ問題(P<math>_{L}</math>)と(D<math>_{L}</math>)の最適解となり最適値が一致する.これを[[ラグランジュの双対性]] (Lagrangian duality) と呼ぶ.(iv)により,<math>\min_{x}\sup_{y}L(x,y)=\max_{y}\inf_{x}L(x,y)</math> が成立すれば,この双対性が保証される.この等式に対する十分条件を述べた定理を[[ミニマックス定理]](minimax theorem) と呼ぶ [1]. 逆に,主問題の目的関数 <math>f_0\, </math> と制約関数 <math>g_i\, </math> がすべて凸で,<math>h_j\, </math> がすべてアフィン関数であるような[[凸計画問題]] (convex programming problem) においては,適当な条件のもとで,問題(P<math>_{L}</math>)の最適解 <math>\bar{x}</math> に対して,<math>\bar{\lambda}\ge{0}</math> であるようなラグランジュ乗数 <math>(\bar{\lambda},\bar{\mu})</math> が存在して,<math>(\bar{x},\bar{\lambda},\bar{\mu})</math> がラグランジュ関数 <math>L\, </math> の鞍点となる.さらに,次のような[[拡張ラグランジュ関数]](augmented Lagrangian function) に基づく双対性も考えられている [2, 3, 4].
  
が成立すれば,$\bar{x}$ と $(\bar{\lambda},\bar{\mu})$ はそれぞれ問題(P$_{L}$)と(D$_{L}$)の最適解となり最適値が一致する.これを[[ラグランジュの双対性]] (Lagrangian duality) と呼ぶ.(iv)により,$\min_{x}\sup_{y}L(x,y)=\max_{y}\inf_{x}L(x,y)$ が成立すれば,この双対性が保証される.この等式に対する十分条件を述べた定理を[[ミニマックス定理]](minimax theorem) と呼ぶ [1]. 逆に,主問題の目的関数 $f_0$ と制約関数 $g_i$ がすべて凸で,$h_j$ がすべてアフィン関数であるような[[凸計画問題]] (convex programming problem) においては,適当な条件のもとで,問題(P$_{L}$)の最適解 $\bar{x}$ に対して,$\bar{\lambda}\ge{0}$ であるようなラグランジュ乗数 $(\bar{\lambda},\bar{\mu})$ が存在して,$(\bar{x},\bar{\lambda},\bar{\mu})$ がラグランジュ関数 $L$ の鞍点となる.さらに,次のような[[拡張ラグランジュ関数]](augmented Lagrangian function) に基づく双対性も考えられている [2, 3, 4].
+
<math>\bar{L}(x,y):=\inf_{u\in{\mathbf{R}^m}}\{\, f(x,u)+r\sigma{(u)}-y^{\top}u\,\}</math>
  
\bar{L}(x,y):=\inf_{u\in{\mathbf{R}^m}}\{\, f(x,u)+r\sigma{(u)}-y^{\top}u\,\}
 
  
ただし,$r$ は正定数,$\sigma:\mathbf{R}^{m}\rightarrow\bar{\mathbf{R}}$ $u\neq{0}$ に対して $0=\sigma{(0)}<\sigma{(u)}$ を満足する下半連続な真凸関数である.関数 $\sigma$ の例としては $\sigma{(u)}:=\frac{1}{2}\|u\|^{2}$ などがある.
+
ただし,<math>r\, </math> は正定数,<math>\sigma:\mathbf{R}^{m}\rightarrow\bar{\mathbf{R}}</math> <math>u\neq{0}</math> に対して <math>0=\sigma{(0)}<\sigma{(u)}</math> を満足する下半連続な真凸関数である.関数 <math>\sigma</math> の例としては <math>\sigma{(u)}:=\frac{1}{2}\|u\|^{2}</math> などがある.
  
  

2007年6月28日 (木) 12:10時点における版

そうついせいりろん}{duality theory}

 双対性理論 (duality theory)は,非線形計画のみならず線形計画,多目的計画,離散凸解析などの分野で主問題とその双対問題の関係,および集合や関数の双対関係を説明する重要な基礎理論である [1, 2, 3, 4].

 「双対」 (dual) と「共役」 (conjugate) は元々同義語として用いられ,数学の関数解析の分野では,ノルム空間 上の有界線形汎関数の全体を の双対空間 (dual space) または共役空間 (conjugate space) といい,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle X^{*}\, } と表して,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x\in{X}\, } における の値を構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \langle x, x^{*}\rangle} または 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x^{*}(x)} と書く.構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle X\, }構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle n\, } 次元実ユークリッド空間 の場合は, は同一視でき, となり,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x^{*}\, } の内積 となる.以下に述べる事柄は,無限次元空間に対しても拡張できるが,ここでは簡単のため に限定して説明する.  空間 上で定義された拡張実数値関数 に対して(ただし,),



で定義される関数 の共役関数 (conjugate function) という.共役関数 に対して,さらにその共役関数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle f^{**}=(f^*)^{*}} を考えることができるが, が下半連続な真凸関数のときには, に一致する.構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle f\, } を対応させる写像をルジャンドル-フェンシェル変換 (Legendre-Fenchel transform) と呼ぶ.

 下半連続な真凸関数 に対して,関数 をそれぞれ 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \varphi{(x)}:=f(x,0)} で定義し,次の問題(P)と(D)を主問題 (primal problem) とその双対問題 (dual problem) と呼ぶ [1, 4].


\begin{center} \end{center}

構文解析に失敗 (不明な関数「\begin{tabular}」): {\displaystyle \begin{tabular}{cll} (P) & \min_{x\in \mathbf{R}^n} & \hspace*{-5mm} \varphi{(x)} \\ (D) & \max_{y\in \mathbf{R}^m} & \hspace*{-5mm} \psi{(y)} \end{tabular} }


また,


構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle U:=\{u\in{\mathbf{R}^m}\,|\:\inf_{x\in{\mathbf{R}^n}}f(x,u)<+\infty\} \quad V:=\{v\in{\mathbf{R}^n}\,|\:\inf_{y\in{\mathbf{R}^m}}f^{*}(v,y)<+\infty\}}


とおくと,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle V\, } は凸集合となる.このとき,以下が成立する.


\medskip \noindent (i) \ \\ (ii) \ 構文解析に失敗 (構文エラー): {\displaystyle 0\in{\mbox{\rm int}\,U}\;} または構文解析に失敗 (構文エラー): {\displaystyle \; 0\in{\mbox{\rm int}\,V} \Longrightarrow\inf_{x}\varphi{(x)}=\sup_{y}\psi{(y)} }


ここで,構文解析に失敗 (構文エラー): {\displaystyle \mbox{\rm int}\,U} の内部を表す.(i)を弱双対定理 (weak duality theorem),(ii)を双対定理 (duality theorem) と呼び,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \inf_{x}\varphi{(x)}=\sup_{y}\psi{(y)}} が満たされるとき,主問題(P)と双対問題(D)の間に双対性 (duality) が成立するという.(i)により,\ なら主問題(P)は実行可能解を持たないが, となる が存在すれば,それぞれ(P)と(D)の最適解となり,強い意味の双対性が成立する.一方, となるとき,主問題と双対問題の間に双対性のギャップ (duality gap) が存在するという.

 主問題(P)において,(ただし, は下半連続な真凸関数で, , 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle c\in{\mathbf{R}^n}} )とすると,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle f^{*}(v,y)=-b^{\top}y+h^{*}(y)+k^{*}(A^{\top}y-c+v)} となり [4],主問題(P)と双対問題(D)はそれぞれ


構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{eqnarray} \mbox{\rm min}_{x\in \mathbf{R}^n} && \hspace*{-5mm} c^{\top}x+k(x)+h(b-Ax) \label{A-B-03+siki1}\\ \mbox{\rm max}_{y\in \mathbf{R}^m} && \hspace*{-5mm} b^{\top}y-h^{*}(y)-k^{*}(A^{\top}y-c) \label{A-B-03+siki2} \end{eqnarray} }


となる.ここで構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle b\in\mbox{\rm int}\,(A\mbox{\rm dom}k+\mbox{\rm dom}h)} または構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle c\in\mbox{\rm int}\,(A^{\top}\mbox{\rm dom}h^{*}-\mbox{\rm dom}k^{*})} が成立すれば,(ii)により主問題 (1) と双対問題 (2) の間に双対性が成立する.(ただし,dom は拡張実数値関数の実効定義域を表す.)これをフェンシェルの双対性 (Fenchel duality) と呼んでいる.通常は,簡略化して 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle c\, }構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle b\, } を零ベクトル,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle -A\, } を恒等写像として,新たに構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle f(x)\, } を凸関数 と凹関数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle f_2(x):=-h(x)} の差で表し,主問題 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \min_{x}\{f_1(x)-f_2(x)\}} に対して,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \max_{y}\{f_{2}^{*}(y)-f_{1}^{*}(y)\}} をフェンシェルの双対問題 (Fenchel dual problem) と呼ぶ.ただし,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle f_{2}^{*}(y):=\inf_{x\in{\mathbf{R}^n}}\{y^{\top}x-f_{2}(x)\}} .双対性は 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mbox{\rm int}\,(\mbox{\rm dom}f_1)\,\cap\, \mbox{\rm int}\,(\mbox{\rm dom}f_2)\neq\emptyset} のとき成立する.また,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle k(x):=\sup_{s\le{0}}x^{\top}s, h(z):=\sup_{w\ge{0}}z^{\top}w} とすると,(1) と(2) は線形計画の主問題と双対問題となる [2, 4].


構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{center} \begin{tabular}{clcll} (P_{LP}) & min. & c^{\top}x & s. t. & x\ge{0}, \ Ax\ge{b}. \\ (D_{LP}) & max. & b^{\top}y & s. t. & y\ge{0}, \ A^{\top}y\le{c}. \end{tabular} \end{center} }


次に,ラグランジュ関数 (Lagrangian function) を

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{equation} \label{A-B-03+siki3} L(x,y):=\inf_{u\in{\mathbf{R}^m}}\{\, f(x,u)-y^{\top}u \,\} \end{equation}}


と定義する.構文解析に失敗 (構文エラー): {\displaystyle -L(x,\cdot)=\left(f(x,\cdot)\right)^{*}, f(x,\cdot)=\left(-L(x,\cdot)\right)^{*}} が成立しているので,

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{equation} \label{A-B-03+siki4} \varphi{(x)}=\sup_{y}L(x,y), \quad\quad \psi{(y)}=\inf_{x}L(x,y) \end{equation}}


となる.通常, すなわち,すべての 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x\in{\mathbf{R}^n}}構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y\in{\mathbf{R}^m}} に対して構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle L(x,\bar{y})\ge{L(\bar{x},\bar{y})}\ge{L(\bar{x},y)}} が成り立つとき, を関数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mathbf{R}^{n}\times{\mathbf{R}^m}} 上での鞍点 (saddle point) と呼ぶ.(4) により,


構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \medskip \noindent (iii) \ \inf_{x}\varphi{(x)}=\inf_{x}\,[\,\sup_{y}L(x,y)\,] \ge\sup_{y}\,[\,\inf_{x}L(x,y)\,]=\sup_{y}\psi{(y)} \\ (iv) \ \varphi{(\bar{x})}=\inf_{x}\varphi{(x)} =\sup_{y}\psi{(y)}=\psi{(\bar{y})} \Longleftrightarrow (\bar{x},\bar{y}) が L の鞍点 \\ \hspace*{1cm} \Longleftrightarrow \min_{x}\sup_{y}L(x,y)=\max_{y}\inf_{x}L(x,y) \Longleftrightarrow \varphi{(\bar{x})}=L(\bar{x},\bar{y})=\psi{(\bar{y})} \medskip \noindent }


が成立する.(iv) を鞍点定理 (saddle point theorem) と呼ぶ.非線形計画問題

構文解析に失敗 (構文エラー): {\displaystyle \mbox{(NLP)} \; \; \mbox{min.} \ f_{0}(x) \quad \mbox{s. t.} \; \; g_{i}(x)\le{0} \ (i=1,\ldots, k), \ h_{j}(x)=0 \ (j=1,\ldots,l),}


(ただし,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle f_0,g_i,h_j}構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mathbf{R}^n} で定義された実数値関数,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle k+l=m} ) に対して,


構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{center} \begin{tabular}{r@{}l} F(x)\,&:=(g_{1}(x),\ldots,g_{k}(x),h_{1}(x),\ldots,h_{l}(x))^{\top}\\ \theta(w)\,&:=\sup_{\lambda,\mu}\{\lambda^{\top}w_{1}+\mu^{\top}w_{2}\,|\, (\lambda,\mu)\in{\mathbf{R}^{k}_{+}\times{\mathbf{R}^{l}}},w=(w_1,w_2)^{\top}\}\\ f(x,u)\,&:=f_{0}(x)+\theta{(F(x)+u)} \end{tabular} \end{center} }


とおくと,(3) により構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle L(x,y)=f_{0}(x)+y^{T}F(x)-\theta^{*}(y)} となり,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y=(\lambda,\mu)=(\lambda_{1},\ldots,\lambda_{k},\mu_{1},\ldots,\mu_{l}) \in{\mathbf{R}^{k}_{+}\times{\mathbf{R}^{l}}}} に対する問題(NLP)のラグランジュ関数は


構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{equation} L(x,\lambda,\mu)=f_{0}(x)+\sum_{i=1}^{k}\lambda_{i}g_{i}(x) +\sum_{j=1}^{l}\mu_{j}h_{j}(x) \end{equation}}


となる.この 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (\lambda,\mu)} をラグランジュ乗数 (Lagrange multipliers) と呼ぶ.このとき,主問題と双対問題は


構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{center} \begin{tabular}{cllll} (P_{L}) & min. & \displaystyle \sup_{\lambda\ge{0},\mu} L(x, \lambda, \mu) & s. t. & \displaystyle{x \in {\mathbf{R}^n}}, \\ (D_{L}) & max. & \displaystyle \inf_{x} L(x,\lambda,\mu) & s. t. & \displaystyle{0 \le \lambda \in {\mathbf{R}^{k}}}, \ \displaystyle \mu \in {\mathbf{R}^{l}}, \end{tabular} \end{center} }


となり,一般に問題(D_{L})をラグランジュの双対問題 (Lagrangian dual problem)と呼ぶ.鞍点定理により,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle L\, } の鞍点 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (\bar{x},\bar{\lambda},\bar{\mu})} が存在すれば,つまり


構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \max_{\lambda\ge{0},\mu}L(\bar{x},\lambda,\mu)=L(\bar{x},\bar{\lambda},\bar{\mu}) =\min_{x}L(x,\bar{\lambda},\bar{\mu})}


が成立すれば,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (\bar{\lambda},\bar{\mu})} はそれぞれ問題(P構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle _{L}} )と(D構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle _{L}} )の最適解となり最適値が一致する.これをラグランジュの双対性 (Lagrangian duality) と呼ぶ.(iv)により,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \min_{x}\sup_{y}L(x,y)=\max_{y}\inf_{x}L(x,y)} が成立すれば,この双対性が保証される.この等式に対する十分条件を述べた定理をミニマックス定理(minimax theorem) と呼ぶ [1]. 逆に,主問題の目的関数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle f_0\, } と制約関数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle g_i\, } がすべて凸で,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle h_j\, } がすべてアフィン関数であるような凸計画問題 (convex programming problem) においては,適当な条件のもとで,問題(P構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle _{L}} )の最適解 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \bar{x}} に対して,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \bar{\lambda}\ge{0}} であるようなラグランジュ乗数 が存在して,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (\bar{x},\bar{\lambda},\bar{\mu})} がラグランジュ関数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle L\, } の鞍点となる.さらに,次のような拡張ラグランジュ関数(augmented Lagrangian function) に基づく双対性も考えられている [2, 3, 4].

構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \bar{L}(x,y):=\inf_{u\in{\mathbf{R}^m}}\{\, f(x,u)+r\sigma{(u)}-y^{\top}u\,\}}


ただし, は正定数,構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \sigma:\mathbf{R}^{m}\rightarrow\bar{\mathbf{R}}}構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle u\neq{0}} に対して 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle 0=\sigma{(0)}<\sigma{(u)}} を満足する下半連続な真凸関数である.関数 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \sigma} の例としては 構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替(最新ブラウザーや補助ツールに推奨): サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \sigma{(u)}:=\frac{1}{2}\|u\|^{2}} などがある.


[1] I. Ekeland and R. Temam, Convex Analysis and Variational Problems, North-Holland, Amsterdam, 1976.

[2] 福島雅夫,『非線形最適化の基礎』, 朝倉書店, 2001.

[3] 今野浩, 山下浩,『非線形計画法』, 日科技連, 1978.

[4] R.T. Rockafellar and R. J-B. Wets, Variational Analysis, Springer, Berlin, 1998.