• 検索結果がありません。

拡張実数値DC最適化問題のラグランジュ型双対性に対する制約想定の考察 (非線形解析学と凸解析学の研究)

N/A
N/A
Protected

Academic year: 2021

シェア "拡張実数値DC最適化問題のラグランジュ型双対性に対する制約想定の考察 (非線形解析学と凸解析学の研究)"

Copied!
5
0
0

読み込み中.... (全文を見る)

全文

(1)165. 拡張実数値 DC 最適化問題の. ラグランジュ型双対性に対する 制約想定の考察 島根大学 大学院総合理工学研究科 村上 卓見 (Takumi Murakami) Interdisciplinary Graduate School of Science and Engineering, Shimane University. 島根大学 大学院総合理工学研究科 角田 侑也 (Yuya Sumida) Interdisciplinary Graduate School of Science and Engineering, Shimane University. 島根大学総合理工学部 黒岩 大史 (Daishi Kuroiwa) Major in Interdisciplinary Science and Engineering, Shimane University. 1. はじめに. 本講究録では次のようなDC 最適化問題に対する [9] における結果を紹介する: (P). minimize subject to. ただし,. I. は. 0. f_{0}(x)-g_{0}(x) f_{i}(x)-g_{i}(x)\leqq 0, \forall i\in I,. を含まない任意の添え字集合とし f_{i},. 真凸関数とする.. I. g_{i}. : \mathbb{R}^{n}arrow \mathbb{R}\cup\{+\infty\}(i\in\{0\}\cup I) は閉. が有限集合である時のこの問題のラグランジュ型双対性に対する制約. 想定の先行研究としては Mart_{1}'nez ‐Legaz and Volle [3], Harada and Kuroiwa [7] などが 挙げられる.これらの研究において,Harada and Kuroiwa [7] の結果は f_{i}(i=1, \ldots, m) が全て実数値関数の時に限り,Martínez‐Legaz and Volle [3] の結果の拡張になっている. 本講究録では,Kuroiwa, Otani and Okano [8] が示した結果に注目し 0\cdot(+\infty)=+\infty と 定めることで, f_{i}(i=1, \ldots, m) が拡張実数値関数の場合においても Harada and Kuroiwa [7] と Martl’nez‐Legaz and Volle [3] の結果の拡張となる,問題 (P) に対するラグランジュ 型双対定理を導き出した.. 2. 準備. まずは準備として,関数 f : \mathbb{R}^{n}arrow \mathbb{R}\cup\{+\infty\} 域をそれぞれ次のように定義する:. に対して, f のエピグラフと f の実行定義. epi f= { (x, r)\in \mathbb{R}^{n}\cross \mathbb{R}|x\in dom f, f(x)\leq r }, dom f=\{x\in \mathbb{R}^{n}|f(x)<+\infty\}. f のエピグラフ epif が凸集合,閉集合,非空のとき, f はそれぞれ凸関数,閉関数,真関数 であるという. f の共役関数を. f^{*}(y)= \sup_{x\in \mathbb{R}^{n} \{\langle x, y\}-f(x)\}.

(2) 166 と定義すると, f^{*} は常に閉凸関数であり, f が真凸関数ならば f^{*} は閉真凸関数であり,また f が閉真凸関数ならば f=f^{**} であることが知られている.また明らかに f(x)+f^{*}(y)\geq\{y, x\rangle. (the Young‐Fenchel inequality) が成立する.任意の. x\in. domf に対して, f の. x. における. 劣微分を次のように定義する:. \partial f(x)=\{x^{*}\in \mathbb{R}^{n}|\{x^{*}, y-x\}\leq f(y)-f(x), \forall y \in \mathbb{R}^{n}\}. このとき,次の同値関係が成立する:. f(x)+f^{*}(y)=\{y, x\}\Leftrightarrow y\in\partial f(x). .. また問題 (P) について,制約集合を S=\{x\in \mathbb{R}^{n}|f_{i}(x)-9i(x)\leqq 0, \forall i\in I\}. とおき,. y_{0}, y_{1}. , . . . , y_{m}\in \mathbb{R}^{n} に対して,(P) の子問題を (P(y_{0}, (y_{i}))). minimize. subject to. f_{0}(x)-\{x, y_{0}\rangle+g_{0}^{*}(y_{0}) f_{\dot{i}}(x)-\{x, y_{i}\}+g_{i}^{*}(y_{\dot{i}})\leqq 0,. \forall i\in I.. とおき,その制約集合を. S(y_{i})=\{x\in \mathbb{R}^{n}|f_{i}(x)-\{x, y_{i}\}+g_{i}^{*}(y_{i})\leqq 0, \forall i\in I\} とおく.. 次に Martínez‐Legaz and Volle [3] と Harada and Kuroiwa [7] の先行研究を紹介する。 これらの先行研究では,上記のようにDC 最適化問題 (P) を凸最適化問題 (P ( y_{0} , (yi)) ) に 分割することが本質である。. 定理2.1 (Martínez‐Legaz Volle [3]) I=\{1, . . .m\}, f_{i,g_{i}} : \mathbb{R}^{n}arrow \mathbb{R}\cup\{+\infty\}(i\in I\cup\{0\}) を凸関数, S=\{x\in \mathbb{R}^{n}|f_{i}(x)\leqq 0, \forall i\in I\}\neq\emptyset とし, g_{i} が S 上で劣微分可 能,. 上で g_{0}=g_{0}^{**} と仮定する.任意の (x_{1}^{*}, . x_{m}^{*}) \in\prod_{1\leq i\leq m}\{g_{\dot{i}}^{*}-f_{\dot{i}}^ {*}\leq 0\} に対して f_{i}(\overline{x})-\{\overline{x}, x_{i}^{*}\}+g_{i}^{*}(x_{i}^{*})<0, \forall i\in I をみたす \overline{x}\in domf_{0} が存在するとき,次の等式が成り S. 立つ.. x\inSx^{*}\indomg_{0}^*g_{i}^*(x_{\doti}^{*)-f_i}^{*(x_ '}^{*) \leq0\lambd R_{+}^m\dot{\imath}nf\{_0}(x)-g_{0}(x)\=inf_{1\leqi m}\dot{\imath}nf\max_{\in}{\begin{ar y}{l m g_{0}^*(x{})+\Sigma\lmbda_{i}g ^{*}(x_i^{*}) i=1 m -(f_{0}+\Sigma\lmbda_{i}f )^{*}(x '+\Sigma\lmbda_{i}x ^{*}) i=1 \end{ar y}\. 定理2.2 (Martínez‐Legaz Volle [3]) I=\{1, m\}, f_{ig_{i}} : \mathbb{R}^{n}arrow \mathbb{R}\cup\{+\infty\}(i\in I\cup\{0\}) を凸関数, S=\{x\in \mathbb{R}^{n}|f_{i}(x)\leqq 0, \forall i\in I\}\neq\emptyset, \Omega=\{(x_{1}^{*}, x_{m}^{*})\in(\mathbb{R}^{n})^{m}| \partial g_{1}^{*}(x_{1}^{*})\cap \cap\partial g_{m}^{*}(x_{m}^{*})\neq\emptyset\} とし,gi が S 上で劣微分可能, S 上で g_{0}=g_{0}^{**} と仮定す る.任意の (x_{1}^{*}, \ldots, x_{m}^{*})\in\Omega に対して f_{i}(\overline{x})-\langle baTx, x_{i}^{*}\rangle+g_{i}^{*}(x_{\dot{i}}^{*})<0, \forall i\in I をみたす \overline{x}\in domf_{0} が存在するとき,次の等式が成り立つ.. \inf_{x\in S}\{f_{0}(x)-g_{0}(x)\}=\inf_{(x^{*},x_{1}^{*},\ldots,x_{m}^{*})\in domg_{0}^{*}\cros \Omega\lambda}\max_{\in \mathbb{R}_{+}^{m}. \{begin{ary}l g_{0}^*(x{})+\sum_{i=1}^m\labd_{i}g ^{*}(x_i^{*}) -(f_{0}+\sum_{i=1}^m\labd_{i}f )^{*}(x +\sum_{i=1}^m\labd_{i} x\dot{i}^*)\ end{ary}\. ..

(3) 167 定理2.3 (R. Harada, D. Kuroiwa, [7]) f_{i,g_{i}} : \mathbb{R}^{n}arrow \mathbb{R}(i\in\{0\}\cup I) を凸関数と し, S\neq\emptyset, \bigcup_{x\in 8}\partial g_{0}(x)\subseteq D_{0}, \bigcup_{x\in S}(\prod_{i=1}^{m}\partial g_{i}(x))\subseteq D とする.任意の (y_{1}, \ldots, y_{m})\in D \cap\prod_{i=1}^{m}domg_{i}^{*} に対して \bullet. \bullet. S(y_{1}, \ldots, y_{m})\neq\emptyset,. cone co. \bigcup_{i=1}^{m. (epi f_{\dot{i}}^{*}-(y_{i}, g_{i}^{*}(y_{i})) ) +\{0\}\cross[0, +\infty ) が閉集合. ならば,次の等式が成り立つ:. x \in S\dot{ \imath} nf(f_{0}(x)-g_{0}(x) =(y0,\ldots,y_{m})\in D_{0}\cros D\lambda_{\dot{i} \geq 0x\in \mathb {R}^{n}\dot{ \imath} nf\max\dot{ \imath} nf 注意2.1. \bullet. \{ begin{ar y}{l f_{0}(x)-\{x y_{0}\+ g_{0}^*(y_{0}) +\sum_{i=1}^{m}\lambda_{i}(f_{i}(x)-\langlex,y_{i}\rangle+ g_{i}^*(y_ {i}) \end{ar y}\. .. f_{i}(i=1, \ldots, m) が実数値関数の時,定理2.3は定理2.1と定理2.2の拡. 張になっている. \bullet. f_{i}(i=1, \ldots, m) が拡張実数値関数の時は,定理2.3は定理2.1と定理2.2の拡張に なっているとは限らない.. 次に,本講究録の主結果を導くために重要な役割を果たす Kuroiwa, Otani and Okano [8] の凸最適化問題に関する結果を紹介する。. lemma 2.1 (Kuroiwa, Otani and Okano [8]) g_{i} : \mathbb{R}^{n}arrow \mathbb{R}\cup\{+\infty\}(i\in I) を閉真 凸関数とし, S=\{x\in \mathbb{R}^{n}|g_{i}(x)\leq 0, \forall i\in I\}\neq\emptyset とする.任意の i\in I に対して, g_{i} が連 続になるような点が S に存在するとき,次は同値である: 1. cone co. 2.. \bigcup_{i\in I}(epig_{i}^{*}\cup epi\delta_{domg_{i}}^{*}). は閉集合. domf \neq\emptyset かつ epif^{*}+epi\delta_{S}^{*} が閉集合となるような任意の閉真凸関数 f:\mathbb{R}^{n}arrow \mathbb{R}\cup\{+\infty\} に対して, S\cap. \inf_{x\inS}f(x)=I_{0}:fin\dot{\imath}te\max\dot{\imath}nfI_{0}\subset q Ix\in\mathb {R}^{n}\{f(x)+\sum_{i\nI_{0}\lambda_{i}g_{i}(x)\}. \lambda\in \mathb {R}_{+}^{I_{0}. ただし, 0\cdot(+\infty)=+\infty.. lemma 2.2 (Kuroiwa, Otani and Okano [8]) g_{i} : \mathbb{R}^{n}arrow \mathbb{R}\cup\{+\infty\}(i\in I) を閉真 凸関数とし, S=\{x\in \mathbb{R}^{n}|g_{i}(x)\leq 0, \forall i\in I\}\neq\emptyset とする.任意の i\in I に対して, g_{i} が連 続になるような点が S に存在するとき,もし Slater condition が成り立つならば,つまり, 任意の i\in I に対して, g_{i}(x_{0})<0 となるような x_{0}\in \mathbb{R}^{n} が存在するならば, cone co. \bigcup_{i\in I}(epig_{\dot{i} ^{*}\cup epi\delta_{domg_{i} ). は閉集合となる..

(4) 168 3. 主結果. 定理3.1 (Kuroiwa, Murakami, Sumida, [9]) f_{i,g_{i}} : \mathbb{R}^{n}arrow \mathbb{R}\cup\{+\infty\}(i\in\{0\}UI) を閉真凸関数とし, \bigcup_{x\in 8}\partial g_{0}(x)\subseteq D_{0}, \bigcup_{x\in S}\prod_{i\in I}\partial g_{i}(x)\subseteq D とする.任意の i\in I に対 して, g_{i} が S で劣微分可能で,任意の (y_{i}) \in D\cap\prod_{i\in I} dom薩に対して, \bullet. \bullet. \bullet. 任意の j\in I に対して, f_{\dot{j} が S(y_{i}) で連続,. S(y_{i})\cap domf_{0}\neq\emptyset かつ epi f_{0}^{*}+epi\delta_{S(y_{i})}^{*} が閉集合, cone co. \bigcup_{i\in I} ((epif_{\dot{i}}^{*}- (y_{i}, g_{i}^{*}(yi))) Uepi\delta_{domf_{i}}^{*}) が閉集合. ならば. x\inS\dot{\imath}nf\{_0}(x)-g_{0}(x)\}=\inf\max\dot{\imath}nf(y0,(y_{i} )\inD_{0}\cros DI_{0}\subset qIx\in\mathb {R}^{n}I_{0}:finte\lambda\in \mathb {R}_{+}^{I_0} \{ begin{ar y}{l f_{0}(x)-\{x y_{0}\+ g_{0}^{*(y_{0}) +\Sigma \lambda_{i}(f_{i}(x)-\{x,y_{i}\+g_{i}^*(y_{i}) i\nI_{0} \end{ar y}\. .. ただし, 0\cdot(+\infty)=+\infty.. 注意3.1任意の (y_{i})\in D\cap H_{i=1}^{m}dom鋪と任意の j\in\{1, m\} に対して,乃が S(y_{i}) で 連続となるような拡張実数値関数 f_{i}(i=1, \ldots, m) に対しては,定理3.1は定理2.1と定 理22の拡張になっている.. 注意3.2定理3.1は定理2.3の拡張になっている. 例3.1 次のような制約関数を持つ DC 最適化問題を考える.. f_{i}(x)=\{ begin{ar ay}{l} -x1 X<1, 0 X\in[-1, ], x-1 x>1 \end{ar ay}. かつ. g_{1}(x)=0.. このとき, f_{1}^{*}(x^{*})=|x^{*}|+\delta_{[-1,1]}(x^{*}), g_{1}^{*}(x^{*})=\delta_{\{0\}}(x^{*}) となり,cone. (epif_{1}^{*}\cup epi_{domf1}^{*}) は. 閉集合となる.よって定理3.1の仮定をみたすことがわかる.一方で 0\in\{g_{1}^{*}-f_{1}^{*}\leq 0\} か. つ,任意の x\in \mathbb{R}^{n} に対して, f_{1}(x)-\langle x, 0\rangle+g_{1}^{*}(0)=f_{1}(x)-\langle x, 0\rangle+\delta_{\{0\}}^{*}(0)=f_{1}(x)\geq 0 となることから,定理2.1と定理2.2の仮定をみたさないことがわかる. 例3.2次のような制約関数を持つ DC 最適化問題を考える.. f_{i}(x)=\{\begin{ar ay}{l } -1 X\in[-1, ], +\infty X\not\in[-1, ] \end{ar ay}. かつ. 91 (x)=0.. このとき, f_{1}^{*}(y_{1})=\Vert y_{1}\Vert+1, g_{1}^{*}(y_{1})=\delta_{\{0\}}(y_{1}) , dom f_{1}=[-1,1], \delta_{domf1}^{*}(y_{1})=\Vert y_{1}\Vert とな り,cone (epif_{1}^{*}\cup epi\delta_{domf1}) は閉集合となる.よって定理3.1の仮定をみたすことがわか. る.一方で cone epi f_{1}^{*}+\{0\}\cross[0, +\infty ) は閉集合とならないことから,定理2.3の仮定をみ たさないことがわかる..

(5) 169 参考文献 [1] V. Jeyakumar, A.M.Rubinov, B.M.Glover, Y.Ishizuka. Inequality Systems and Global optimization. J. Math. Anal. Appl. 202 (1996), no. 3, 900‐919. [2] B. Lemaire. Duality in reverse convex optimization. SIAM J. Optim. 8 (1998), no. 4, 1029‐1037.. [3] J.‐E. Martı’nez‐Legaz, M. Volle.Duality in DC programming: the case of several DC consrtaints.J. Math. Anal. Appl., 237 (1999), pp. 657‐671. [4] M.A. Goberna, V. Jeyakumar, M.A. López.Necessary and sufficient constraint qual‐ ifications for solvability of systems of infinite convex inequalities Nonlinear Anal.,. 68 (2008), pp. 1184‐1194. [5] V. Jeyakumar,G.Y.Li.New dual constraint qualification characterizing zero duality gaps of convex programs and semidefinite programs.Nonlinear Anal. 71 (2009), no. 12, e2239‐e2249.. [6] Y. Fujiwara, D. Kuroiwa. Lagrange duality in canonical DC programming. J. Math. Anal. Appl. 408 (2013), no. 2, 476‐483. [7] R. Harada, D. Kuroiwa. Lagrange‐type duality in DC programming. J. Math. Anal. Appl. 418 (2014), no. 1, 415‐424. [8] D. Kuroiwa, H. Otani, K. Okano. A necessary and sufficient constraint qualifica‐ tion for the Lagrange‐duality including the Slater condition in extended real‐valued convex optimization problems, preprint.. [9] D. Kuroiwa, T. Murakami, Y. Sumida. A constraint qualification for a Lagrange‐ type duality of extended real valued DC optimization problems, preprint.

(6)

参照

関連したドキュメント

東京大学 大学院情報理工学系研究科 数理情報学専攻. [email protected]

Max-flow min-cut theorem and faster algorithms in a circular disk failure model, INFOCOM 2014...

ポートフォリオ最適化問題の改良代理制約法による対話型解法 仲川 勇二 関西大学 * 伊佐田 百合子 関西学院大学 井垣 伸子

Research Institute for Mathematical Sciences, Kyoto University...

東北大学大学院医学系研究科の運動学分野門間陽樹講師、早稲田大学の川上

【対応者】 :David M Ingram 教授(エディンバラ大学工学部 エネルギーシステム研究所). Alistair G。L。 Borthwick

話題提供者: 河﨑佳子 神戸大学大学院 人間発達環境学研究科 話題提供者: 酒井邦嘉# 東京大学大学院 総合文化研究科 話題提供者: 武居渡 金沢大学

向井 康夫 : 東北大学大学院 生命科学研究科 助教 牧野 渡 : 東北大学大学院 生命科学研究科 助教 占部 城太郎 :