数学において、 ミラー降下法は 微分可能関数 の 局所的最小値 を見つけるための 反復的な 最適化 アルゴリズム です 。
勾配降下法 や 乗法重み などのアルゴリズムを一般化します 。
歴史
ミラー降下法は、 1983年に ネミロフスキー とユディンによって最初に提案されました 。[1]
モチベーション
微分可能関数に 学習率の列を適用した 勾配降下法 では、 の局所的最小値の 推測から始めて、 次のような
列を考える。
(
η
ん
)
ん
≥
0
{\displaystyle (\eta _{n})_{n\geq 0}}
ふ
{\displaystyle F}
x
0
{\displaystyle \mathbf {x} _{0}}
ふ
、
{\displaystyle F,}
x
0
、
x
1
、
x
2
、
…
{\displaystyle \mathbf {x} _{0},\mathbf {x} _{1},\mathbf {x} _{2},\ldots }
x
ん
+
1
=
x
ん
−
η
ん
∇
ふ
(
x
ん
)
、
ん
≥
0.
{\displaystyle \mathbf {x} _{n+1}=\mathbf {x} _{n}-\eta _{n}\nabla F(\mathbf {x} _{n}),\ n\geq 0.}
これを次のように言い換えることができる。
x
ん
+
1
=
引数
分
x
(
ふ
(
x
ん
)
+
∇
ふ
(
x
ん
)
T
(
x
−
x
ん
)
+
1
2
η
ん
‖
x
−
x
ん
‖
2
)
{\displaystyle \mathbf {x} _{n+1}=\arg \min _{\mathbf {x} }\left(F(\mathbf {x} _{n})+\nabla F(\mathbf { x} _{n})^{T}(\mathbf {x} -\mathbf {x} _{n})+{\frac {1}{2\eta _{n}}}\|\mathbf {x} -\mathbf {x} _{n}\|^{2}\right)}
言い換えると、 近接項 を追加することで、 での 1 次近似を に最小化します 。
x
ん
+
1
{\displaystyle \mathbf {x} _{n+1}}
ふ
{\displaystyle F}
x
ん
{\displaystyle \mathbf {x} _{n}}
‖
x
−
x
ん
‖
2
{\displaystyle \|\mathbf {x} -\mathbf {x} _{n}\|^{2}}
この二乗ユークリッド距離項はブレグマン距離 の特別な例です。他のブレグマン距離を使用すると 、特定の形状での最適化により適した ヘッジ などの他のアルゴリズムが生成されます。 [2] [3]
凸集合 上で最適化する 凸関数 と、 上の何らかのノルムが与えられます 。
ふ
{\displaystyle f}
け
⊂
R
ん
{\displaystyle K\subset \mathbb {R} ^{n}}
‖
⋅
‖
{\displaystyle \|\cdot \|}
R
ん
{\displaystyle \mathbb {R} ^{n}}
また、 微分可能な凸関数 も与えられます。これは 、 与えられたノルムに関して 強凸です。これは 距離生成関数 と呼ばれ、その 勾配は ミラー マップ として知られています 。
h
:
R
ん
→
R
{\displaystyle h\colon \mathbb {R} ^{n}\to \mathbb {R} }
α
{\displaystyle \alpha}
∇
h
:
R
ん
→
R
ん
{\displaystyle \nabla h\colon \mathbb {R} ^{n}\to \mathbb {R} ^{n}}
初期 から始めて 、Mirror Descent の各反復で次のようになります。
x
0
∈
け
{\displaystyle x_{0}\in K}
双対空間へのマップ:
θ
t
←
∇
h
(
x
t
)
{\displaystyle \theta _{t}\leftarrow \nabla h(x_{t})}
勾配ステップを使用してデュアル空間で更新します。
θ
t
+
1
←
θ
t
−
η
t
∇
ふ
(
x
t
)
{\displaystyle \theta _{t+1}\leftarrow \theta _{t}-\eta _{t}\nabla f(x_{t})}
原始空間にマップし直す:
x
t
+
1
′
←
(
∇
h
)
−
1
(
θ
t
+
1
)
{\displaystyle x'_{t+1}\leftarrow (\nabla h)^{-1}(\theta _{t+1})}
実行可能領域 に投影し直します 。 ここで、 Bregman ダイバージェンス となります 。
け
{\displaystyle K}
x
t
+
1
←
1つの
r
グ
分
x
∈
け
だ
h
(
x
|
|
x
t
+
1
′
)
{\displaystyle x_{t+1}\leftarrow \mathrm {arg} \min _{x\in K}D_{h}(x||x'_{t+1})}
だ
h
{\displaystyle D_{h}}
拡張機能
オンライン最適化 設定におけるミラー降下法は 、オンラインミラー降下法(OMD)として知られています。 [4]
参照
参考文献
^ アルカディ・ネミロフスキー、デイビッド・ユディン。最適化における問題の複雑さと方法の効率。ジョン・ワイリー・アンド・サンズ、1983年
^ Nemirovski, Arkadi (2012) チュートリアル: 大規模な決定論的および確率的凸最適化のためのミラー降下アルゴリズム。https://www2.isye.gatech.edu/~nemirovs/COLT2012Tut.pdf
^ 「ミラー降下アルゴリズム」。tlienart.github.io . 2022年7月10日 閲覧 。
^ Fang, Huang; Harvey, Nicholas JA; Portella, Victor S.; Friedlander, Michael P. (2021-09-03). 「 オンラインミラー降下法とデュアル平均化:動的なケースでのペースの維持」。arXiv : 2006.02585 [cs.LG]。