この 章の 目標
凸集合・凸包・錐・ 多面体を 定義し、例と 反例を 挙げられる
閉凸集合への 射影の 存在・ 一意性と 特徴づけを 証明し、そこから 分離定理・支持超平面定理を 導ける
微分 可能な 関数の 凸性を 1 次条件・2 次条件で 判定でき、定義域が 開凸集合と いう 仮定の 意味を 説明できる
凸関数の 局所最小は 大域最小である ことを 証明し、凸性を 保つ演算で 凸性を 確かめられる
強凸性と L L L -平滑性の 同値な 特徴づけを 証明し、条件数の 意味を 説明できる
劣勾配を 計算し、最適性条件 0 ∈ ∂ f ( x ) 0 \in \partial f(x) 0 ∈ ∂ f ( x ) を 使える
前提 :第1章 、01-calculus 第4章 (1 変数の 凸関数)、 01-calculus 第7章 、02-linear-algebra 第8章 (正定値性)。2.7 節では 01-calculus 第5章 の 微分積分学の 基本定理と 第8章 の 平均値の 不等式を 使う。
第1章で 予告した「凸最適化問題では 局所最小点が 大域最小点に なる」ことを 証明し、それを 支える 理論を 整える。柱は、閉凸集合への 射影から 導かれる 分離定理(線形計画の 双対定理( 第3章 )や 強双対性( 第4章 )の 証明の 核)、凸性の 判定法、勾配法の 速さを 決める 強凸性 と L L L -平滑性 (第5章 )の 三つである。最後に、微分できない 凸関数の ための 劣勾配 と 共役関数を 導入する。
⟨ x , y ⟩ = x ⊤ y \langle x, y \rangle = x^{\top}y ⟨ x , y ⟩ = x ⊤ y 、∥ x ∥ \lVert x \rVert ∥ x ∥ は ユークリッドノルムと する。対称行列 A , B A, B A , B に ついて、 A ⪰ B A \succeq B A ⪰ B は A − B A - B A − B が 半正定値、 A ≻ B A \succ B A ≻ B は 正定値である ことを 表す。
2.1 凸集合
凸集合は 第1章 定義 1.24 で 定義した。
定義 2.1 (凸結合・凸包・錐・ 多面体) λ i ≥ 0 \lambda_i \geq 0 λ i ≥ 0 , ∑ i = 1 k λ i = 1 \sum_{i=1}^{k}\lambda_i = 1 ∑ i = 1 k λ i = 1 の とき、 ∑ i λ i x i \sum_i \lambda_ix_i ∑ i λ i x i を x 1 , … , x k x_1, \dots, x_k x 1 , … , x k の 凸結合 (convex combination) と いう。 S ⊂ R n S \subset \mathbb{R}^n S ⊂ R n の 有限個の 元の 凸結合全体を 凸包 (convex hull) conv S \operatorname{conv} S conv S と いう。「 x ∈ K x \in K x ∈ K , λ ≥ 0 \lambda \geq 0 λ ≥ 0 ならば λ x ∈ K \lambda x \in K λ x ∈ K 」を 満たす K K K を 錐 (cone)、凸な 錐を 凸錐と いう。有限個の 閉半空間の 共通部分 { x ∣ A x ≤ b } \lbrace x \mid Ax \leq b \rbrace { x ∣ A x ≤ b } を 多面体 (polyhedron) と いう。
命題 2.2 (1) 凸集合の 任意個の 共通部分は 凸である。(2) 凸集合の アフィン写像 x ↦ A x + b x \mapsto Ax + b x ↦ A x + b に よる 像と 逆像は 凸である。(3) 凸集合 C C C は、C C C の 元の 凸結合を すべて 含む。(4) conv S \operatorname{conv} S conv S は S S S を 含む 最小の 凸集合である。
証明. (1) は 明らか。(2) は A ( ( 1 − t ) x + t y ) + b = ( 1 − t ) ( A x + b ) + t ( A y + b ) A((1 - t)x + ty) + b = (1 - t)(Ax + b) + t(Ay + b) A (( 1 − t ) x + t y ) + b = ( 1 − t ) ( A x + b ) + t ( A y + b ) に よる。(3) k k k に ついての 帰納法。 λ k = 1 \lambda_k = 1 λ k = 1 なら和は x k x_k x k である。λ k < 1 \lambda_k < 1 λ k < 1 なら ∑ i = 1 k λ i x i = ( 1 − λ k ) z + λ k x k \sum_{i=1}^{k}\lambda_ix_i = (1 - \lambda_k)z + \lambda_kx_k ∑ i = 1 k λ i x i = ( 1 − λ k ) z + λ k x k , z = ∑ i = 1 k − 1 λ i 1 − λ k x i ∈ C z = \sum_{i=1}^{k-1}\frac{\lambda_i}{1 - \lambda_k}x_i \in C z = ∑ i = 1 k − 1 1 − λ k λ i x i ∈ C (帰納法の 仮定)。(4) 凸結合の 凸結合は 凸結合なので conv S \operatorname{conv} S conv S は 凸で、 S S S を 含む凸集合は (3) より conv S \operatorname{conv} S conv S を 含む。 □ \square □
例 2.3 (1) 超平面・ 閉半空間 { a ⊤ x ≤ b } \lbrace a^{\top}x \leq b \rbrace { a ⊤ x ≤ b } ・球(三角不等式に よる)・ 多面体は 凸である。線形計画の 実行可能領域(第1章 例 1.2・1.3)は 多面体である。(2) 確率単体 { w ≥ 0 ∣ ∑ i w i = 1 } = conv { e 1 , … , e n } \lbrace w \geq 0 \mid \sum_i w_i = 1 \rbrace = \operatorname{conv}\lbrace e_1, \dots, e_n \rbrace { w ≥ 0 ∣ ∑ i w i = 1 } = conv { e 1 , … , e n } は 空売りなしの ポートフォリオ全体である。(3) 非負象限 R + n \mathbb{R}^n_{+} R + n 、二次錐 { ( x , t ) ∣ ∥ x ∥ ≤ t } \lbrace (x, t) \mid \lVert x \rVert \leq t \rbrace {( x , t ) ∣ ∥ x ∥ ≤ t } 、半正定値対称行列全体(一次 不等式 v ⊤ X v ≥ 0 v^{\top}Xv \geq 0 v ⊤ X v ≥ 0 の 共通部分)は 凸錐で、線形計画・ 二次錐計画・半正定値計画の 基礎に なる。(4) 交わらない 2 つの 球の 和集合、 Z n \mathbb{Z}^n Z n 、球面は 凸でない。
凸集合 C C C の 点 x x x が、x = ( 1 − t ) y + t z x = (1 - t)y + tz x = ( 1 − t ) y + t z (y , z ∈ C y, z \in C y , z ∈ C , 0 < t < 1 0 < t < 1 0 < t < 1 )なら y = z = x y = z = x y = z = x と なる とき、 端点 (extreme point) と いう。有界な 多面体は 有限個の 端点( 頂点 )の 凸包に 等しい(ミンコフスキー–ワイルの 定理。証明は Schrijver, Theory of Linear and Integer Programming などを 参照)。
2.2 射影定理
02-linear-algebra 第7章 定理 7.17 の 正射影を、部分 空間から 閉凸集合に 一般化する。
定理 2.4 (射影定理, projection theorem)C ⊂ R n C \subset \mathbb{R}^n C ⊂ R n を 空でない 閉凸集合と する。各 x ∈ R n x \in \mathbb{R}^n x ∈ R n に 対し、 ∥ x − p ∥ = min y ∈ C ∥ x − y ∥ \lVert x - p \rVert = \min_{y \in C}\lVert x - y \rVert ∥ x − p ∥ = min y ∈ C ∥ x − y ∥ と なる p ∈ C p \in C p ∈ C が ただ 一つ 存在する。これを P C ( x ) P_C(x) P C ( x ) と 書き、 x x x の C C C への 射影と いう。 p ∈ C p \in C p ∈ C に ついて
p = P C ( x ) ⟺ ⟨ x − p , y − p ⟩ ≤ 0 ( ∀ y ∈ C ) (2.1) p = P_C(x) \iff \langle x - p, y - p \rangle \leq 0 \quad (\forall y \in C) \tag{2.1} p = P C ( x ) ⟺ ⟨ x − p , y − p ⟩ ≤ 0 ( ∀ y ∈ C ) ( 2.1 )
であり、∥ P C ( x ) − P C ( x ′ ) ∥ ≤ ∥ x − x ′ ∥ \lVert P_C(x) - P_C(x') \rVert \leq \lVert x - x' \rVert ∥ P C ( x ) − P C ( x ′ )∥ ≤ ∥ x − x ′ ∥ (非拡大性 )が 成り立つ。
証明. 存在:y ↦ ∥ x − y ∥ 2 y \mapsto \lVert x - y \rVert^2 y ↦ ∥ x − y ∥ 2 は 連続で 強圧的なので、第1章 定理 1.11 より 閉集合 C C C 上で 最小値を とる。一意性: p 1 , p 2 p_1, p_2 p 1 , p 2 が ともに 最小距離 d d d を 与えるなら、中点 m ∈ C m \in C m ∈ C に ついて 中線定理より
∥ x − m ∥ 2 = 1 2 ∥ x − p 1 ∥ 2 + 1 2 ∥ x − p 2 ∥ 2 − 1 4 ∥ p 1 − p 2 ∥ 2 = d 2 − 1 4 ∥ p 1 − p 2 ∥ 2 \lVert x - m \rVert^2 = \frac{1}{2}\lVert x - p_1 \rVert^2 + \frac{1}{2}\lVert x - p_2 \rVert^2 - \frac{1}{4}\lVert p_1 - p_2 \rVert^2 = d^2 - \frac{1}{4}\lVert p_1 - p_2 \rVert^2 ∥ x − m ∥ 2 = 2 1 ∥ x − p 1 ∥ 2 + 2 1 ∥ x − p 2 ∥ 2 − 4 1 ∥ p 1 − p 2 ∥ 2 = d 2 − 4 1 ∥ p 1 − p 2 ∥ 2
で、d d d の 最小性から p 1 = p 2 p_1 = p_2 p 1 = p 2 。(2.1) の ⇒:y ∈ C y \in C y ∈ C , 0 < t ≤ 1 0 < t \leq 1 0 < t ≤ 1 なら p + t ( y − p ) ∈ C p + t(y - p) \in C p + t ( y − p ) ∈ C なので 0 ≤ ∥ x − p − t ( y − p ) ∥ 2 − ∥ x − p ∥ 2 = − 2 t ⟨ x − p , y − p ⟩ + t 2 ∥ y − p ∥ 2 0 \leq \lVert x - p - t(y - p) \rVert^2 - \lVert x - p \rVert^2 = -2t\langle x - p, y - p \rangle + t^2\lVert y - p \rVert^2 0 ≤ ∥ x − p − t ( y − p ) ∥ 2 − ∥ x − p ∥ 2 = − 2 t ⟨ x − p , y − p ⟩ + t 2 ∥ y − p ∥ 2 。t t t で 割って t → + 0 t \to +0 t → + 0 と する。⇐: ∥ x − y ∥ 2 = ∥ x − p ∥ 2 − 2 ⟨ x − p , y − p ⟩ + ∥ y − p ∥ 2 ≥ ∥ x − p ∥ 2 \lVert x - y \rVert^2 = \lVert x - p \rVert^2 - 2\langle x - p, y - p \rangle + \lVert y - p \rVert^2 \geq \lVert x - p \rVert^2 ∥ x − y ∥ 2 = ∥ x − p ∥ 2 − 2 ⟨ x − p , y − p ⟩ + ∥ y − p ∥ 2 ≥ ∥ x − p ∥ 2 。非拡大性:p = P C ( x ) p = P_C(x) p = P C ( x ) , p ′ = P C ( x ′ ) p' = P_C(x') p ′ = P C ( x ′ ) と して (2.1) を ( x , p , y = p ′ ) (x, p, y = p') ( x , p , y = p ′ ) と ( x ′ , p ′ , y = p ) (x', p', y = p) ( x ′ , p ′ , y = p ) に 使って 足すと、 ∥ p − p ′ ∥ 2 ≤ ⟨ x − x ′ , p − p ′ ⟩ ≤ ∥ x − x ′ ∥ ∥ p − p ′ ∥ \lVert p - p' \rVert^2 \leq \langle x - x', p - p' \rangle \leq \lVert x - x' \rVert\lVert p - p' \rVert ∥ p − p ′ ∥ 2 ≤ ⟨ x − x ′ , p − p ′ ⟩ ≤ ∥ x − x ′ ∥ ∥ p − p ′ ∥ 。□ \square □
(2.1) は「x − p x - p x − p と y − p y - p y − p の なす角が 直角以上」と いう 意味である。 C C C が 閉でないと 最も 近い 点が ない ことが あり(開球と 外の 点)、凸でないと 一意性が 崩れる(球面と 中心)。
例 2.5 (射影の 計算)(1) 箱 { l ≤ y ≤ u } \lbrace l \leq y \leq u \rbrace { l ≤ y ≤ u } への 射影は 成分ごとの 切り 詰め min ( max ( x i , l i ) , u i ) \min(\max(x_i, l_i), u_i) min ( max ( x i , l i ) , u i ) ([ 0 , 1 ] 2 [0, 1]^2 [ 0 , 1 ] 2 への ( 2 , − 0.5 ) (2, -0.5) ( 2 , − 0.5 ) の 射影は ( 1 , 0 ) (1, 0) ( 1 , 0 ) )。(2) 球 { ∥ y ∥ ≤ r } \lbrace \lVert y \rVert \leq r \rbrace {∥ y ∥ ≤ r } への 射影は、 ∥ x ∥ > r \lVert x \rVert > r ∥ x ∥ > r なら r x / ∥ x ∥ rx/\lVert x \rVert r x / ∥ x ∥ 。(3) 半空間 { a ⊤ y ≤ b } \lbrace a^{\top}y \leq b \rbrace { a ⊤ y ≤ b } への 射影は x − max ( 0 , a ⊤ x − b ) ∥ a ∥ 2 a x - \frac{\max(0, a^{\top}x - b)}{\lVert a \rVert^2}a x − ∥ a ∥ 2 m a x ( 0 , a ⊤ x − b ) a (a ⊤ x > b a^{\top}x > b a ⊤ x > b なら x − p = c a x - p = ca x − p = c a , c > 0 c > 0 c > 0 , a ⊤ p = b a^{\top}p = b a ⊤ p = b で、⟨ x − p , y − p ⟩ = c ( a ⊤ y − b ) ≤ 0 \langle x - p, y - p \rangle = c(a^{\top}y - b) \leq 0 ⟨ x − p , y − p ⟩ = c ( a ⊤ y − b ) ≤ 0 )。射影が 簡単な 集合の 上では、第5章の 射影勾配法が 使える。
2.3 分離定理と 支持超平面定理
交わらない 凸集合の 間には 超平面を 引ける。これが 実行不能性や 最適性の 証明書(第3章・第4章の 双対変数)を 作る 道具に なる。仮定(閉・コンパクト・開)と 分離の 強さの 対応に 注意する。
定理 2.6 (点と 閉凸集合の 分離) C C C を 空でない 閉凸集合、 x 0 ∉ C x_0 \notin C x 0 ∈ / C と する。この とき a ≠ 0 a \neq 0 a = 0 が あって sup y ∈ C a ⊤ y < a ⊤ x 0 \sup_{y \in C} a^{\top}y < a^{\top}x_0 sup y ∈ C a ⊤ y < a ⊤ x 0 と なる。
証明. p = P C ( x 0 ) p = P_C(x_0) p = P C ( x 0 ) , a = x 0 − p ≠ 0 a = x_0 - p \neq 0 a = x 0 − p = 0 と おくと、(2.1) より y ∈ C y \in C y ∈ C で a ⊤ ( y − p ) ≤ 0 a^{\top}(y - p) \leq 0 a ⊤ ( y − p ) ≤ 0 、すな わち a ⊤ y ≤ a ⊤ p = a ⊤ x 0 − ∥ a ∥ 2 a^{\top}y \leq a^{\top}p = a^{\top}x_0 - \lVert a \rVert^2 a ⊤ y ≤ a ⊤ p = a ⊤ x 0 − ∥ a ∥ 2 。□ \square □
定理 2.7 (閉凸集合と コンパクト凸集合の 強分離) C C C を 空でない 閉凸集合、 K K K を 空でない コンパクト凸集合とし、 C ∩ K = ∅ C \cap K = \emptyset C ∩ K = ∅ と する。この とき a ≠ 0 a \neq 0 a = 0 が あって sup y ∈ C a ⊤ y < min z ∈ K a ⊤ z \sup_{y \in C} a^{\top}y < \min_{z \in K} a^{\top}z sup y ∈ C a ⊤ y < min z ∈ K a ⊤ z と なる。
証明. D = C − K = { y − z ∣ y ∈ C , z ∈ K } D = C - K = \lbrace y - z \mid y \in C, z \in K \rbrace D = C − K = { y − z ∣ y ∈ C , z ∈ K } は 凸で 0 ∉ D 0 \notin D 0 ∈ / D 。D D D は 閉である : y k − z k → w y_k - z_k \to w y k − z k → w (y k ∈ C y_k \in C y k ∈ C , z k ∈ K z_k \in K z k ∈ K )なら、K K K の コンパクト性から 部分列で z k j → z ∈ K z_{k_j} \to z \in K z k j → z ∈ K 、すると y k j → w + z ∈ C y_{k_j} \to w + z \in C y k j → w + z ∈ C で、w ∈ D w \in D w ∈ D 。定理 2.6 を D D D と 0 0 0 に 使うと、 a ≠ 0 a \neq 0 a = 0 , ε > 0 \varepsilon > 0 ε > 0 が あって すべての y ∈ C y \in C y ∈ C , z ∈ K z \in K z ∈ K で a ⊤ y ≤ a ⊤ z − ε a^{\top}y \leq a^{\top}z - \varepsilon a ⊤ y ≤ a ⊤ z − ε 。□ \square □
例 2.8 (コンパクト性は 外せない) C = { ( x , y ) ∣ x > 0 , x y ≥ 1 } C = \lbrace (x, y) \mid x > 0, \ xy \geq 1 \rbrace C = {( x , y ) ∣ x > 0 , x y ≥ 1 } (凸関数 1 / x 1/x 1/ x の エピグラフ)と K = { ( x , y ) ∣ y ≤ 0 } K = \lbrace (x, y) \mid y \leq 0 \rbrace K = {( x , y ) ∣ y ≤ 0 } は 交わらない 閉凸集合だが、 ( t , 1 / t ) ∈ C (t, 1/t) \in C ( t , 1/ t ) ∈ C と ( t , 0 ) ∈ K (t, 0) \in K ( t , 0 ) ∈ K の 距離は 0 0 0 に 近づく。 inf K a ⊤ z > − ∞ \inf_K a^{\top}z > -\infty inf K a ⊤ z > − ∞ と なるのは a = ( 0 , a 2 ) a = (0, a_2) a = ( 0 , a 2 ) , a 2 < 0 a_2 < 0 a 2 < 0 の ときだけで、その とき inf K a ⊤ z = 0 = sup C a ⊤ y \inf_K a^{\top}z = 0 = \sup_C a^{\top}y inf K a ⊤ z = 0 = sup C a ⊤ y だから、定理 2.7 の 形の 分離は できない( C C C と K K K を 入れ替えても 同様)。
閉でない 凸集合も 扱う ために、有限次元の 次の 性質を 使う。
補題 2.9 凸集合 C ⊂ R n C \subset \mathbb{R}^n C ⊂ R n に ついて、閉包の 内部は C C C の 内部に 等しい : ( C ‾ ) ∘ = C ∘ (\overline{C})^{\circ} = C^{\circ} ( C ) ∘ = C ∘ 。
証明. ⊃ \supset ⊃ は 明らか。 x 0 ∈ ( C ‾ ) ∘ x_0 \in (\overline{C})^{\circ} x 0 ∈ ( C ) ∘ とし、閉球 B ‾ ( x 0 , δ n ) ⊂ C ‾ \overline{B}(x_0, \delta\sqrt{n}) \subset \overline{C} B ( x 0 , δ n ) ⊂ C と なる δ > 0 \delta > 0 δ > 0 を とる。 v i = x 0 + δ e i v_i = x_0 + \delta e_i v i = x 0 + δ e i (1 ≤ i ≤ n 1 \leq i \leq n 1 ≤ i ≤ n )、v 0 = x 0 − δ ( e 1 + ⋯ + e n ) v_0 = x_0 - \delta(e_1 + \cdots + e_n) v 0 = x 0 − δ ( e 1 + ⋯ + e n ) は この 閉球に 属し、 x 0 = 1 n + 1 ∑ i = 0 n v i x_0 = \frac{1}{n+1}\sum_{i=0}^{n}v_i x 0 = n + 1 1 ∑ i = 0 n v i 。点 w 0 , … , w n w_0, \dots, w_n w 0 , … , w n に 対し、 ( w 0 , 1 ) , … , ( w n , 1 ) ∈ R n + 1 (w_0, 1), \dots, (w_n, 1) \in \mathbb{R}^{n+1} ( w 0 , 1 ) , … , ( w n , 1 ) ∈ R n + 1 を 列と する 正方行列を M ( w ) M(w) M ( w ) と する。 v i − v 0 = δ ( e i + e 1 + ⋯ + e n ) v_i - v_0 = \delta(e_i + e_1 + \cdots + e_n) v i − v 0 = δ ( e i + e 1 + ⋯ + e n ) (i ≥ 1 i \geq 1 i ≥ 1 )は 一次独立(係数行列 δ ( I + 11 ⊤ ) \delta(I + \mathbf{1}\mathbf{1}^{\top}) δ ( I + 1 1 ⊤ ) の 固有値は δ ( n + 1 ) \delta(n + 1) δ ( n + 1 ) と δ \delta δ )なので M ( v ) M(v) M ( v ) は 正則で、クラメルの 公式より λ ( w ) = M ( w ) − 1 ( x 0 , 1 ) \lambda(w) = M(w)^{-1}(x_0, 1) λ ( w ) = M ( w ) − 1 ( x 0 , 1 ) は v v v の 近くで w w w に ついて 連続、 λ ( v ) = ( 1 n + 1 , … , 1 n + 1 ) \lambda(v) = (\frac{1}{n+1}, \dots, \frac{1}{n+1}) λ ( v ) = ( n + 1 1 , … , n + 1 1 ) 。よって w i ∈ C w_i \in C w i ∈ C を v i v_i v i の 十分 近くにとれば、 M ( w ) M(w) M ( w ) は 正則で λ ( w ) \lambda(w) λ ( w ) の 成分は すべて 正と なり、 x 0 = ∑ i λ i ( w ) w i x_0 = \sum_i \lambda_i(w)w_i x 0 = ∑ i λ i ( w ) w i , ∑ i λ i ( w ) = 1 \sum_i \lambda_i(w) = 1 ∑ i λ i ( w ) = 1 だから x 0 ∈ C x_0 \in C x 0 ∈ C (命題 2.2 (3))。開集合 ( C ‾ ) ∘ (\overline{C})^{\circ} ( C ) ∘ が C C C に 含まれたので、 C ∘ C^{\circ} C ∘ に 含まれる。 □ \square □
定理 2.10 (支持超平面定理, supporting hyperplane theorem)C ⊂ R n C \subset \mathbb{R}^n C ⊂ R n を 空でない 凸集合、 x 0 ∉ C ∘ x_0 \notin C^{\circ} x 0 ∈ / C ∘ (たとえば C C C の 境界点)と する。この とき a ≠ 0 a \neq 0 a = 0 が あって、すべての y ∈ C y \in C y ∈ C で a ⊤ y ≤ a ⊤ x 0 a^{\top}y \leq a^{\top}x_0 a ⊤ y ≤ a ⊤ x 0 と なる。
証明. 補題 2.9 より x 0 ∉ ( C ‾ ) ∘ x_0 \notin (\overline{C})^{\circ} x 0 ∈ / ( C ) ∘ なので、C ‾ \overline{C} C に 属さない 点列 x k → x 0 x_k \to x_0 x k → x 0 が とれる。凸集合の 閉包は 凸だから( C C C の 点列 y k → y y_k \to y y k → y , z k → z z_k \to z z k → z に ついて ( 1 − t ) y k + t z k → ( 1 − t ) y + t z (1 - t)y_k + tz_k \to (1 - t)y + tz ( 1 − t ) y k + t z k → ( 1 − t ) y + t z )、定理 2.6 より ∥ a k ∥ = 1 \lVert a_k \rVert = 1 ∥ a k ∥ = 1 で a k ⊤ y < a k ⊤ x k a_k^{\top}y < a_k^{\top}x_k a k ⊤ y < a k ⊤ x k (y ∈ C ‾ y \in \overline{C} y ∈ C )と なる a k a_k a k が ある。部分列で a k → a a_k \to a a k → a (∥ a ∥ = 1 \lVert a \rVert = 1 ∥ a ∥ = 1 )とし、y ∈ C y \in C y ∈ C を 固定して 極限を とればよい。 □ \square □
定理 2.11 (分離定理, separating hyperplane theorem)C 1 , C 2 ⊂ R n C_1, C_2 \subset \mathbb{R}^n C 1 , C 2 ⊂ R n を 交わらない 空でない 凸集合と する。
a ≠ 0 a \neq 0 a = 0 が あって、すべての x ∈ C 1 x \in C_1 x ∈ C 1 , y ∈ C 2 y \in C_2 y ∈ C 2 で a ⊤ x ≤ a ⊤ y a^{\top}x \leq a^{\top}y a ⊤ x ≤ a ⊤ y と なる。
さらに C 1 C_1 C 1 が 開集合なら、 α = sup C 1 a ⊤ x \alpha = \sup_{C_1}a^{\top}x α = sup C 1 a ⊤ x に ついて、すべての x ∈ C 1 x \in C_1 x ∈ C 1 , y ∈ C 2 y \in C_2 y ∈ C 2 で a ⊤ x < α ≤ a ⊤ y a^{\top}x < \alpha \leq a^{\top}y a ⊤ x < α ≤ a ⊤ y と なる。
証明. (1) D = C 1 − C 2 D = C_1 - C_2 D = C 1 − C 2 は 凸で 0 ∉ D 0 \notin D 0 ∈ / D なので、定理 2.10 を D D D と 0 0 0 に 使う。(2) x ∈ C 1 x \in C_1 x ∈ C 1 なら 小さい ε > 0 \varepsilon > 0 ε > 0 で x + ε a ∈ C 1 x + \varepsilon a \in C_1 x + ε a ∈ C 1 だから、a ⊤ x < a ⊤ ( x + ε a ) ≤ α a^{\top}x < a^{\top}(x + \varepsilon a) \leq \alpha a ⊤ x < a ⊤ ( x + ε a ) ≤ α 。□ \square □
補足
有限次元である ことは、定理 2.10(したがって 定理 2.11 (1))の 証明の 2 か所で 使った。補題 2.9(座標ベクトル e 1 , … , e n e_1, \dots, e_n e 1 , … , e n で x 0 x_0 x 0 を 囲む)と、単位ベクトルの 列 a k a_k a k から 収束する 部分列を とる ところ(単位球面の コンパクト性)である。実際、無限次元では 定理 2.11 (1) は 成り立たない。2 乗和が 有限な 実数列の 空間 ℓ 2 \ell^2 ℓ 2 の 中で、有限個を 除く 成分が 0 0 0 の 数列全体 c 00 c_{00} c 00 (稠密な 部分 空間。 10-functional-analysis 第1章 例 1.7)と、その 外の 1 点 x 0 x_0 x 0 は 交わらない 凸集合だが、 c 00 c_{00} c 00 上で 上に 有界な 連続線形汎関数は c 00 c_{00} c 00 上で 0 0 0 、したがって 稠密性から 恒等的に 0 0 0 なので、両者を 分離する 超平面は ない。無限次元の ノルム空間では、一方が 開集合の 場合や 閉凸集合と 1 点の 場合に 分離が 成り立つ(ハーン–バナッハの 定理の 幾何形。 10-functional-analysis 第4章 定理 4.7・4.8)。射影定理の 無限次元版は ヒルベルト空間の 最近 点定理である( 10-functional-analysis 第2章 定理 2.6)。
2.4 凸関数と その 判定
定義 2.12 凸集合 C C C 上の 関数 f f f が、x ≠ y x \neq y x = y , 0 < t < 1 0 < t < 1 0 < t < 1 で f ( ( 1 − t ) x + t y ) < ( 1 − t ) f ( x ) + t f ( y ) f((1 - t)x + ty) < (1 - t)f(x) + tf(y) f (( 1 − t ) x + t y ) < ( 1 − t ) f ( x ) + t f ( y ) を 満たすとき 狭義凸 、− f -f − f が 凸の とき 凹 と いう。 epi f = { ( x , s ) ∈ C × R ∣ f ( x ) ≤ s } \operatorname{epi} f = \lbrace (x, s) \in C \times \mathbb{R} \mid f(x) \leq s \rbrace epi f = {( x , s ) ∈ C × R ∣ f ( x ) ≤ s } を エピグラフ (epigraph) と いう。
命題 2.13 C C C を 凸集合、 f : C → R f\colon C \to \mathbb{R} f : C → R と する。
f f f が 凸 ⟺ \iff ⟺ epi f \operatorname{epi} f epi f が 凸集合。
f f f が 凸なら、下位集合 { x ∈ C ∣ f ( x ) ≤ α } \lbrace x \in C \mid f(x) \leq \alpha \rbrace { x ∈ C ∣ f ( x ) ≤ α } は 凸である(逆は 成り立たない : ∣ x ∣ \sqrt{\lvert x \rvert} ∣ x ∣ )。
f f f が 凸 ⟺ \iff ⟺ 任意の x ∈ C x \in C x ∈ C , d d d に ついて、 t ↦ f ( x + t d ) t \mapsto f(x + td) t ↦ f ( x + t d ) が 区間 { t ∣ x + t d ∈ C } \lbrace t \mid x + td \in C \rbrace { t ∣ x + t d ∈ C } 上で凸。
(イェンセンの 不等式) f f f が 凸で λ i ≥ 0 \lambda_i \geq 0 λ i ≥ 0 , ∑ i λ i = 1 \sum_i \lambda_i = 1 ∑ i λ i = 1 なら f ( ∑ i λ i x i ) ≤ ∑ i λ i f ( x i ) f(\sum_i \lambda_ix_i) \leq \sum_i \lambda_if(x_i) f ( ∑ i λ i x i ) ≤ ∑ i λ i f ( x i ) 。
証明. 1 は、( x , s ) , ( y , u ) ∈ epi f (x, s), (y, u) \in \operatorname{epi} f ( x , s ) , ( y , u ) ∈ epi f に ついて f ( ( 1 − t ) x + t y ) ≤ ( 1 − t ) s + t u f((1 - t)x + ty) \leq (1 - t)s + tu f (( 1 − t ) x + t y ) ≤ ( 1 − t ) s + t u と なる ことが 凸性と 同値であることに よる。2・3 は 定義から 従う。4 は 01-calculus 第4章 定理 4.31 と 同じ 帰納法に よる。 □ \square □
定理 2.14 (1 次条件)U ⊂ R n U \subset \mathbb{R}^n U ⊂ R n を 開凸集合、 f : U → R f\colon U \to \mathbb{R} f : U → R を 微分可能と する。次は 同値である。
f f f は 凸である。
すべての x , y ∈ U x, y \in U x , y ∈ U で f ( y ) ≥ f ( x ) + ⟨ ∇ f ( x ) , y − x ⟩ f(y) \geq f(x) + \langle \nabla f(x), y - x \rangle f ( y ) ≥ f ( x ) + ⟨ ∇ f ( x ) , y − x ⟩ 。
すべての x , y ∈ U x, y \in U x , y ∈ U で ⟨ ∇ f ( x ) − ∇ f ( y ) , x − y ⟩ ≥ 0 \langle \nabla f(x) - \nabla f(y), x - y \rangle \geq 0 ⟨ ∇ f ( x ) − ∇ f ( y ) , x − y ⟩ ≥ 0 (勾配の 単調性 )。
証明. (1⇒2) 0 < t ≤ 1 0 < t \leq 1 0 < t ≤ 1 で f ( x + t ( y − x ) ) − f ( x ) ≤ t ( f ( y ) − f ( x ) ) f(x + t(y - x)) - f(x) \leq t(f(y) - f(x)) f ( x + t ( y − x )) − f ( x ) ≤ t ( f ( y ) − f ( x )) 。t t t で 割って t → + 0 t \to +0 t → + 0 と すると、左辺は 方向微分 ⟨ ∇ f ( x ) , y − x ⟩ \langle \nabla f(x), y - x \rangle ⟨ ∇ f ( x ) , y − x ⟩ に 収束する(01-calculus 第7章 命題 7.11)。(2⇒1) z = ( 1 − t ) x + t y z = (1 - t)x + ty z = ( 1 − t ) x + t y とし、2 を ( z , x ) (z, x) ( z , x ) と ( z , y ) (z, y) ( z , y ) に 使って 1 − t 1 - t 1 − t 倍と t t t 倍して 足すと、 ( 1 − t ) ( x − z ) + t ( y − z ) = 0 (1 - t)(x - z) + t(y - z) = 0 ( 1 − t ) ( x − z ) + t ( y − z ) = 0 より ( 1 − t ) f ( x ) + t f ( y ) ≥ f ( z ) (1 - t)f(x) + tf(y) \geq f(z) ( 1 − t ) f ( x ) + t f ( y ) ≥ f ( z ) 。(2⇒3) 2 を ( x , y ) (x, y) ( x , y ) と ( y , x ) (y, x) ( y , x ) に 使って 足す。(3⇒1) x t = x + t ( y − x ) x_t = x + t(y - x) x t = x + t ( y − x ) , φ ( t ) = f ( x t ) \varphi(t) = f(x_t) φ ( t ) = f ( x t ) は [ 0 , 1 ] [0, 1] [ 0 , 1 ] を 含む開区間で 微分可能で、 φ ′ ( t ) = ⟨ ∇ f ( x t ) , y − x ⟩ \varphi'(t) = \langle \nabla f(x_t), y - x \rangle φ ′ ( t ) = ⟨ ∇ f ( x t ) , y − x ⟩ 。s < t s < t s < t なら φ ′ ( t ) − φ ′ ( s ) = 1 t − s ⟨ ∇ f ( x t ) − ∇ f ( x s ) , x t − x s ⟩ ≥ 0 \varphi'(t) - \varphi'(s) = \frac{1}{t - s}\langle \nabla f(x_t) - \nabla f(x_s), x_t - x_s \rangle \geq 0 φ ′ ( t ) − φ ′ ( s ) = t − s 1 ⟨ ∇ f ( x t ) − ∇ f ( x s ) , x t − x s ⟩ ≥ 0 で、φ ′ \varphi' φ ′ が 単調増加なので φ \varphi φ は 凸(01-calculus 第4章 定理 4.30)。 φ ( t ) ≤ ( 1 − t ) φ ( 0 ) + t φ ( 1 ) \varphi(t) \leq (1 - t)\varphi(0) + t\varphi(1) φ ( t ) ≤ ( 1 − t ) φ ( 0 ) + tφ ( 1 ) が 求める 不等式である。 □ \square □
2 は「接平面が グラフ全体の 下に ある」こと、つまり 局所的な 情報 ∇ f ( x ) \nabla f(x) ∇ f ( x ) が 大域的な 下界を 与える ことを 意味し、凸最適化の 要である。
定理 2.15 (2 次条件)U ⊂ R n U \subset \mathbb{R}^n U ⊂ R n を 開凸集合、 f f f を U U U 上の C 2 C^2 C 2 級関数と する。
f f f が 凸 ⟺ \iff ⟺ すべての x ∈ U x \in U x ∈ U で ∇ 2 f ( x ) ⪰ 0 \nabla^2 f(x) \succeq 0 ∇ 2 f ( x ) ⪰ 0 。
すべての x ∈ U x \in U x ∈ U で ∇ 2 f ( x ) ≻ 0 \nabla^2 f(x) \succ 0 ∇ 2 f ( x ) ≻ 0 ならば f f f は 狭義凸である(逆は 成り立たない : x 4 x^4 x 4 )。
証明. φ ( t ) = f ( x + t d ) \varphi(t) = f(x + td) φ ( t ) = f ( x + t d ) は 0 0 0 を 含む開区間で C 2 C^2 C 2 級で、φ ′ ′ ( t ) = d ⊤ ∇ 2 f ( x + t d ) d \varphi''(t) = d^{\top}\nabla^2 f(x + td)d φ ′′ ( t ) = d ⊤ ∇ 2 f ( x + t d ) d (01-calculus 第7章 7.6 節)。(1) ⇒:φ \varphi φ は 凸なので(命題 2.13 (3)) φ ′ ′ ( 0 ) ≥ 0 \varphi''(0) \geq 0 φ ′′ ( 0 ) ≥ 0 (01-calculus 第4章 定理 4.30)。⇐:d = y − x d = y - x d = y − x と すると φ ′ ′ ≥ 0 \varphi'' \geq 0 φ ′′ ≥ 0 より φ \varphi φ は 凸で、命題 2.13 (3) より f f f は 凸。(2) d ≠ 0 d \neq 0 d = 0 なら φ ′ \varphi' φ ′ は 狭義単調増加なので、平均値の 定理より 0 < t < 1 0 < t < 1 0 < t < 1 で φ ( t ) − φ ( 0 ) t = φ ′ ( ξ ) < φ ′ ( η ) = φ ( 1 ) − φ ( t ) 1 − t \frac{\varphi(t) - \varphi(0)}{t} = \varphi'(\xi) < \varphi'(\eta) = \frac{\varphi(1) - \varphi(t)}{1 - t} t φ ( t ) − φ ( 0 ) = φ ′ ( ξ ) < φ ′ ( η ) = 1 − t φ ( 1 ) − φ ( t ) (ξ < t < η \xi < t < \eta ξ < t < η )。これは φ ( t ) < ( 1 − t ) φ ( 0 ) + t φ ( 1 ) \varphi(t) < (1 - t)\varphi(0) + t\varphi(1) φ ( t ) < ( 1 − t ) φ ( 0 ) + tφ ( 1 ) と 同値である。 □ \square □
例 2.16 (定義域が 開である ことの 意味) f ( x 1 , x 2 ) = x 1 2 − x 2 2 f(x_1, x_2) = x_1^2 - x_2^2 f ( x 1 , x 2 ) = x 1 2 − x 2 2 を 内部が 空の 凸集合 C = { ( x 1 , 0 ) } C = \lbrace (x_1, 0) \rbrace C = {( x 1 , 0 )} に 制限すると x 1 2 x_1^2 x 1 2 で 凸だが、 ∇ 2 f = diag ( 2 , − 2 ) \nabla^2 f = \operatorname{diag}(2, -2) ∇ 2 f = diag ( 2 , − 2 ) は 半正定値でない。定理 2.15 の「⇒」は、定義域が 開でないと 成り立たない。
例 2.17 (凸関数の 例)(1) アフィン関数、ノルム。(2) 二次関数 1 2 x ⊤ Q x + c ⊤ x \frac{1}{2}x^{\top}Qx + c^{\top}x 2 1 x ⊤ Q x + c ⊤ x は 凸 ⟺ Q ⪰ 0 \iff Q \succeq 0 ⟺ Q ⪰ 0 で、分散 w ⊤ Σ w w^{\top}\Sigma w w ⊤ Σ w や ∥ A x − b ∥ 2 \lVert Ax - b \rVert^2 ∥ A x − b ∥ 2 は 凸。(3) e a x e^{ax} e a x 、− log x -\log x − log x 、x log x x\log x x log x (x > 0 x > 0 x > 0 )、∣ x ∣ p \lvert x \rvert^p ∣ x ∣ p (p ≥ 1 p \geq 1 p ≥ 1 )、log-sum-exp(問題 2.3)。(4) log ( 1 + e t ) \log(1 + e^{t}) log ( 1 + e t ) は 2 階導関数 σ ( 1 − σ ) \sigma(1 - \sigma) σ ( 1 − σ ) が 正( σ \sigma σ は シグモイド関数)なので 凸で、ロジスティック回帰の 損失(第1章 例 1.6)は ∇ 2 ℓ ( x ) = ∑ i σ ( a i ⊤ x ) ( 1 − σ ( a i ⊤ x ) ) a i a i ⊤ ⪰ 0 \nabla^2\ell(x) = \sum_i \sigma(a_i^{\top}x)(1 - \sigma(a_i^{\top}x))a_ia_i^{\top} \succeq 0 ∇ 2 ℓ ( x ) = ∑ i σ ( a i ⊤ x ) ( 1 − σ ( a i ⊤ x )) a i a i ⊤ ⪰ 0 より凸。(5) x 3 x^3 x 3 、sin x \sin x sin x 、x 1 x 2 x_1x_2 x 1 x 2 は 凸でない。
2.5 凸最適化の 基本定理
定理 2.18 C C C を 凸集合、 f : C → R f\colon C \to \mathbb{R} f : C → R を 凸関数と する。
f f f の 局所最小点は 大域最小点である。
最小点の 全体は 凸集合である。
f f f が 狭義凸なら、最小点は(存在すれば)ただ 一つである。
証明. (1) x ∗ x^{\ast} x ∗ を 局所最小点とし、 f ( y ) < f ( x ∗ ) f(y) < f(x^{\ast}) f ( y ) < f ( x ∗ ) と なる y ∈ C y \in C y ∈ C が あったとする。 0 < t ≤ 1 0 < t \leq 1 0 < t ≤ 1 で x t = x ∗ + t ( y − x ∗ ) ∈ C x_t = x^{\ast} + t(y - x^{\ast}) \in C x t = x ∗ + t ( y − x ∗ ) ∈ C かつ f ( x t ) ≤ ( 1 − t ) f ( x ∗ ) + t f ( y ) < f ( x ∗ ) f(x_t) \leq (1 - t)f(x^{\ast}) + tf(y) < f(x^{\ast}) f ( x t ) ≤ ( 1 − t ) f ( x ∗ ) + t f ( y ) < f ( x ∗ ) で、t → 0 t \to 0 t → 0 で x t → x ∗ x_t \to x^{\ast} x t → x ∗ だから 局所最小性に 反する。(2) 最小値を p ∗ p^{\ast} p ∗ と すると、最小点の 全体は 下位集合 { f ≤ p ∗ } \lbrace f \leq p^{\ast} \rbrace { f ≤ p ∗ } である。(3) 最小点 x ≠ y x \neq y x = y が あれば f ( x + y 2 ) < p ∗ f(\frac{x + y}{2}) < p^{\ast} f ( 2 x + y ) < p ∗ と なり矛盾。 □ \square □
定理 2.19 (凸最適化の 最適性条件) U ⊂ R n U \subset \mathbb{R}^n U ⊂ R n を 開凸集合、 f : U → R f\colon U \to \mathbb{R} f : U → R を 微分可能な 凸関数と する。
x ∗ ∈ U x^{\ast} \in U x ∗ ∈ U が U U U 上の 最小点 ⟺ \iff ⟺ ∇ f ( x ∗ ) = 0 \nabla f(x^{\ast}) = 0 ∇ f ( x ∗ ) = 0 。
C ⊂ U C \subset U C ⊂ U を 凸集合と する。 x ∗ ∈ C x^{\ast} \in C x ∗ ∈ C が C C C 上の 最小点 ⟺ \iff ⟺ すべての y ∈ C y \in C y ∈ C で ⟨ ∇ f ( x ∗ ) , y − x ∗ ⟩ ≥ 0 \langle \nabla f(x^{\ast}), y - x^{\ast} \rangle \geq 0 ⟨ ∇ f ( x ∗ ) , y − x ∗ ⟩ ≥ 0 。
証明. (1) ⇒ は 第1章 定理 1.17。(1)・ (2) の ⇐ は 1 次条件(定理 2.14)に よる。(2) ⇒:ある y ∈ C y \in C y ∈ C で ⟨ ∇ f ( x ∗ ) , y − x ∗ ⟩ < 0 \langle \nabla f(x^{\ast}), y - x^{\ast} \rangle < 0 ⟨ ∇ f ( x ∗ ) , y − x ∗ ⟩ < 0 なら、y − x ∗ y - x^{\ast} y − x ∗ は 降下方向で(第1章 補題 1.16)、小さい t > 0 t > 0 t > 0 で x ∗ + t ( y − x ∗ ) ∈ C x^{\ast} + t(y - x^{\ast}) \in C x ∗ + t ( y − x ∗ ) ∈ C の 値は f ( x ∗ ) f(x^{\ast}) f ( x ∗ ) より 小さく、矛盾する。 □ \square □
(2) で f ( y ) = 1 2 ∥ y − x ∥ 2 f(y) = \frac{1}{2}\lVert y - x \rVert^2 f ( y ) = 2 1 ∥ y − x ∥ 2 と すると、射影の 特徴づけ (2.1) に なる。
注意
凸性が 保証するのは「局所最小 = 大域最小」と「停留点 = 最小点」であって、最小点の 存在 ではない(e x e^{x} e x や 完全分離の ロジスティック回帰(第1章 例 1.13))。存在は 強圧性などで 別に 確かめる。
2.6 凸性を 保つ演算
命題 2.20 f , f i , g f, f_i, g f , f i , g を 凸関数と する(定義域は 凸集合と する)。次の 関数は 凸である。
非負の 重みつき和 ∑ i α i f i \sum_i \alpha_if_i ∑ i α i f i (α i ≥ 0 \alpha_i \geq 0 α i ≥ 0 )。
アフィン写像との 合成 f ( A x + b ) f(Ax + b) f ( A x + b ) 。
各点ごとの 上限 sup i ∈ I f i ( x ) \sup_{i \in I} f_i(x) sup i ∈ I f i ( x ) (I I I は 任意の 集合で、上限は 有限と する)。
h : R → R h\colon \mathbb{R} \to \mathbb{R} h : R → R が 凸かつ単調非減少の ときの h ( g ( x ) ) h(g(x)) h ( g ( x )) 。
証明. 1・2 は 定義から 従う。3:各 i i i で f i ( ( 1 − t ) x + t y ) ≤ ( 1 − t ) sup j f j ( x ) + t sup j f j ( y ) f_i((1 - t)x + ty) \leq (1 - t)\sup_j f_j(x) + t\sup_j f_j(y) f i (( 1 − t ) x + t y ) ≤ ( 1 − t ) sup j f j ( x ) + t sup j f j ( y ) で、左辺の 上限を とる。4: h ( g ( ( 1 − t ) x + t y ) ) ≤ h ( ( 1 − t ) g ( x ) + t g ( y ) ) ≤ ( 1 − t ) h ( g ( x ) ) + t h ( g ( y ) ) h(g((1 - t)x + ty)) \leq h((1 - t)g(x) + tg(y)) \leq (1 - t)h(g(x)) + th(g(y)) h ( g (( 1 − t ) x + t y )) ≤ h (( 1 − t ) g ( x ) + t g ( y )) ≤ ( 1 − t ) h ( g ( x )) + t h ( g ( y )) 。□ \square □
たとえば max i ( a i ⊤ x + b i ) \max_i(a_i^{\top}x + b_i) max i ( a i ⊤ x + b i ) 、ヒンジ損失 max ( 0 , 1 − t ) \max(0, 1 - t) max ( 0 , 1 − t ) 、∥ x ∥ 1 \lVert x \rVert_1 ∥ x ∥ 1 、対称行列の 最大固有値 λ max ( X ) = max ∥ v ∥ = 1 v ⊤ X v \lambda_{\max}(X) = \max_{\lVert v \rVert = 1}v^{\top}Xv λ m a x ( X ) = max ∥ v ∥ = 1 v ⊤ X v (02-linear-algebra 第8章 命題 8.29)は 3 より、ℓ ( x ) + λ 2 ∥ x ∥ 2 \ell(x) + \frac{\lambda}{2}\lVert x \rVert^2 ℓ ( x ) + 2 λ ∥ x ∥ 2 や ラッソの 目的関数 1 2 ∥ A x − b ∥ 2 + λ ∥ x ∥ 1 \frac{1}{2}\lVert Ax - b \rVert^2 + \lambda\lVert x \rVert_1 2 1 ∥ A x − b ∥ 2 + λ ∥ x ∥ 1 は 1・2 より 凸である。凸関数の 最小値・差・積は 凸とは 限らない( min ( x 2 , ( x − 2 ) 2 ) \min(x^2, (x - 2)^2) min ( x 2 , ( x − 2 ) 2 ) 、− x 2 -x^2 − x 2 、x ⋅ x 2 x \cdot x^2 x ⋅ x 2 )。
ヒント
実務では
凸性は ヘッセ行列の 計算より、命題 2.20 の 規則の 組み合わせで 確かめる ほうが 確実である。代表的な 凸最適化の モデリング用ソフトウェアは、式が この 種の 規則で 組み立てられている ときに 限って 受け付け、凸性を 機械的に 保証する(disciplined convex programming)。規則で 説明できない 式は、変数変換で 凸に 書き直すか、凸な 近似で 置き換えられないかを 考える。
2.7 強凸性と L L L -平滑性
勾配法の 速さを 論じるには、 f f f を 上下から 二次関数で 挟む定量的な 情報が 要る。
定義 2.21 U ⊂ R n U \subset \mathbb{R}^n U ⊂ R n を 開凸集合、 f : U → R f\colon U \to \mathbb{R} f : U → R , μ > 0 \mu > 0 μ > 0 , L > 0 L > 0 L > 0 と する。
f ( x ) − μ 2 ∥ x ∥ 2 f(x) - \frac{\mu}{2}\lVert x \rVert^2 f ( x ) − 2 μ ∥ x ∥ 2 が 凸の とき、 f f f は μ \mu μ -強凸 (μ \mu μ -strongly convex) であると いう。
f f f が 微分可能で、すべての x , y ∈ U x, y \in U x , y ∈ U で ∥ ∇ f ( x ) − ∇ f ( y ) ∥ ≤ L ∥ x − y ∥ \lVert \nabla f(x) - \nabla f(y) \rVert \leq L\lVert x - y \rVert ∥ ∇ f ( x ) − ∇ f ( y )∥ ≤ L ∥ x − y ∥ の とき、 f f f は L L L -平滑 (L L L -smooth) であると いう。
定理 2.22 (強凸性の 特徴づけ) U U U を 開凸集合、 f : U → R f\colon U \to \mathbb{R} f : U → R を 微分可能と する。次は 同値である。
f f f は μ \mu μ -強凸である。
すべての x , y ∈ U x, y \in U x , y ∈ U で f ( y ) ≥ f ( x ) + ⟨ ∇ f ( x ) , y − x ⟩ + μ 2 ∥ y − x ∥ 2 f(y) \geq f(x) + \langle \nabla f(x), y - x \rangle + \frac{\mu}{2}\lVert y - x \rVert^2 f ( y ) ≥ f ( x ) + ⟨ ∇ f ( x ) , y − x ⟩ + 2 μ ∥ y − x ∥ 2 。
すべての x , y ∈ U x, y \in U x , y ∈ U で ⟨ ∇ f ( x ) − ∇ f ( y ) , x − y ⟩ ≥ μ ∥ x − y ∥ 2 \langle \nabla f(x) - \nabla f(y), x - y \rangle \geq \mu\lVert x - y \rVert^2 ⟨ ∇ f ( x ) − ∇ f ( y ) , x − y ⟩ ≥ μ ∥ x − y ∥ 2 。
f f f が C 2 C^2 C 2 級なら、これらは「すべての x ∈ U x \in U x ∈ U で ∇ 2 f ( x ) ⪰ μ I \nabla^2 f(x) \succeq \mu I ∇ 2 f ( x ) ⪰ μ I 」とも 同値である。
証明. g ( x ) = f ( x ) − μ 2 ∥ x ∥ 2 g(x) = f(x) - \frac{\mu}{2}\lVert x \rVert^2 g ( x ) = f ( x ) − 2 μ ∥ x ∥ 2 に 定理 2.14・2.15 を 適用する。 ∇ g ( x ) = ∇ f ( x ) − μ x \nabla g(x) = \nabla f(x) - \mu x ∇ g ( x ) = ∇ f ( x ) − μx , ∇ 2 g = ∇ 2 f − μ I \nabla^2 g = \nabla^2 f - \mu I ∇ 2 g = ∇ 2 f − μ I と、恒等式 ∥ y ∥ 2 − ∥ x ∥ 2 − 2 ⟨ x , y − x ⟩ = ∥ y − x ∥ 2 \lVert y \rVert^2 - \lVert x \rVert^2 - 2\langle x, y - x \rangle = \lVert y - x \rVert^2 ∥ y ∥ 2 − ∥ x ∥ 2 − 2 ⟨ x , y − x ⟩ = ∥ y − x ∥ 2 に より、 g g g の 1 次条件は 2 に、∇ g \nabla g ∇ g の 単調性は 3 に、 ∇ 2 g ⪰ 0 \nabla^2 g \succeq 0 ∇ 2 g ⪰ 0 は ∇ 2 f ⪰ μ I \nabla^2 f \succeq \mu I ∇ 2 f ⪰ μ I に 書き換わる。 □ \square □
定理 2.23 (L L L -平滑性の 特徴づけ)
U U U を 開凸集合、 f : U → R f\colon U \to \mathbb{R} f : U → R を L L L -平滑と すると、すべての x , y ∈ U x, y \in U x , y ∈ U で
f ( y ) ≤ f ( x ) + ⟨ ∇ f ( x ) , y − x ⟩ + L 2 ∥ y − x ∥ 2 (2.2) f(y) \leq f(x) + \langle \nabla f(x), y - x \rangle + \frac{L}{2}\lVert y - x \rVert^2 \tag{2.2} f ( y ) ≤ f ( x ) + ⟨ ∇ f ( x ) , y − x ⟩ + 2 L ∥ y − x ∥ 2 ( 2.2 )
が 成り立つ( 降下補題 )。f f f が C 2 C^2 C 2 級なら、L L L -平滑性は すべての x ∈ U x \in U x ∈ U で − L I ⪯ ∇ 2 f ( x ) ⪯ L I -LI \preceq \nabla^2 f(x) \preceq LI − L I ⪯ ∇ 2 f ( x ) ⪯ L I と なる ことと 同値である。
2. f : R n → R f\colon \mathbb{R}^n \to \mathbb{R} f : R n → R を 微分可能な 凸関数と すると、次は 同値である。(a) L L L -平滑。(b) すべての x , y x, y x , y で (2.2)。(c) すべての x , y x, y x , y で f ( y ) ≥ f ( x ) + ⟨ ∇ f ( x ) , y − x ⟩ + 1 2 L ∥ ∇ f ( y ) − ∇ f ( x ) ∥ 2 f(y) \geq f(x) + \langle \nabla f(x), y - x \rangle + \frac{1}{2L}\lVert \nabla f(y) - \nabla f(x) \rVert^2 f ( y ) ≥ f ( x ) + ⟨ ∇ f ( x ) , y − x ⟩ + 2 L 1 ∥ ∇ f ( y ) − ∇ f ( x ) ∥ 2 。(d) すべての x , y x, y x , y で ⟨ ∇ f ( x ) − ∇ f ( y ) , x − y ⟩ ≥ 1 L ∥ ∇ f ( x ) − ∇ f ( y ) ∥ 2 \langle \nabla f(x) - \nabla f(y), x - y \rangle \geq \frac{1}{L}\lVert \nabla f(x) - \nabla f(y) \rVert^2 ⟨ ∇ f ( x ) − ∇ f ( y ) , x − y ⟩ ≥ L 1 ∥ ∇ f ( x ) − ∇ f ( y ) ∥ 2 。
特に C 2 C^2 C 2 級の 凸関数では、 L L L -平滑性は 0 ⪯ ∇ 2 f ( x ) ⪯ L I 0 \preceq \nabla^2 f(x) \preceq LI 0 ⪯ ∇ 2 f ( x ) ⪯ L I (∀ x \forall x ∀ x )と 同値である。
証明. (1) ∇ f \nabla f ∇ f は 連続なので t ↦ f ( x + t ( y − x ) ) t \mapsto f(x + t(y - x)) t ↦ f ( x + t ( y − x )) は [ 0 , 1 ] [0, 1] [ 0 , 1 ] で C 1 C^1 C 1 級であり、微分積分学の 基本定理(01-calculus 第5章 定理 5.13)と コーシー–シュワルツの 不等式より
f ( y ) − f ( x ) − ⟨ ∇ f ( x ) , y − x ⟩ = ∫ 0 1 ⟨ ∇ f ( x + t ( y − x ) ) − ∇ f ( x ) , y − x ⟩ d t ≤ ∫ 0 1 L t ∥ y − x ∥ 2 d t = L 2 ∥ y − x ∥ 2 f(y) - f(x) - \langle \nabla f(x), y - x \rangle = \int_0^1 \langle \nabla f(x + t(y - x)) - \nabla f(x), y - x \rangle\,dt \leq \int_0^1 Lt\lVert y - x \rVert^2\,dt = \frac{L}{2}\lVert y - x \rVert^2 f ( y ) − f ( x ) − ⟨ ∇ f ( x ) , y − x ⟩ = ∫ 0 1 ⟨ ∇ f ( x + t ( y − x )) − ∇ f ( x ) , y − x ⟩ d t ≤ ∫ 0 1 L t ∥ y − x ∥ 2 d t = 2 L ∥ y − x ∥ 2
C 2 C^2 C 2 級の 場合、 L L L -平滑なら ∇ 2 f ( x ) d = lim t → 0 ( ∇ f ( x + t d ) − ∇ f ( x ) ) / t \nabla^2 f(x)d = \lim_{t \to 0}(\nabla f(x + td) - \nabla f(x))/t ∇ 2 f ( x ) d = lim t → 0 ( ∇ f ( x + t d ) − ∇ f ( x )) / t より 作用素ノルムで ∥ ∇ 2 f ( x ) ∥ ≤ L \lVert \nabla^2 f(x) \rVert \leq L ∥ ∇ 2 f ( x )∥ ≤ L 。逆に ∥ ∇ 2 f ∥ ≤ L \lVert \nabla^2 f \rVert \leq L ∥ ∇ 2 f ∥ ≤ L なら、平均値の 不等式(01-calculus 第8章 定理 8.4)を ∇ f \nabla f ∇ f に 使えばよい。対称行列の 作用素ノルムは 固有値の 絶対値の 最大値なので( 02-linear-algebra 第8章 命題 8.21・問題 8.6)、∥ ∇ 2 f ( x ) ∥ ≤ L ⟺ − L I ⪯ ∇ 2 f ( x ) ⪯ L I \lVert \nabla^2 f(x) \rVert \leq L \iff -LI \preceq \nabla^2 f(x) \preceq LI ∥ ∇ 2 f ( x )∥ ≤ L ⟺ − L I ⪯ ∇ 2 f ( x ) ⪯ L I 。
(2) (a⇒b) は (1)。(b⇒c) x x x を 固定し ϕ ( z ) = f ( z ) − ⟨ ∇ f ( x ) , z ⟩ \phi(z) = f(z) - \langle \nabla f(x), z \rangle ϕ ( z ) = f ( z ) − ⟨ ∇ f ( x ) , z ⟩ と おく。 ϕ \phi ϕ は 凸で ∇ ϕ ( x ) = 0 \nabla\phi(x) = 0 ∇ ϕ ( x ) = 0 だから x x x は ϕ \phi ϕ の 最小点で(定理 2.19)、 ϕ \phi ϕ も (2.2) を 満たす(一次関数の 項は 両辺で 打ち消し合う)。 z = y − 1 L ∇ ϕ ( y ) z = y - \frac{1}{L}\nabla\phi(y) z = y − L 1 ∇ ϕ ( y ) に (2.2) を 使うと
ϕ ( x ) ≤ ϕ ( z ) ≤ ϕ ( y ) − 1 L ∥ ∇ ϕ ( y ) ∥ 2 + 1 2 L ∥ ∇ ϕ ( y ) ∥ 2 = ϕ ( y ) − 1 2 L ∥ ∇ f ( y ) − ∇ f ( x ) ∥ 2 \phi(x) \leq \phi(z) \leq \phi(y) - \frac{1}{L}\lVert \nabla\phi(y) \rVert^2 + \frac{1}{2L}\lVert \nabla\phi(y) \rVert^2 = \phi(y) - \frac{1}{2L}\lVert \nabla f(y) - \nabla f(x) \rVert^2 ϕ ( x ) ≤ ϕ ( z ) ≤ ϕ ( y ) − L 1 ∥ ∇ ϕ ( y ) ∥ 2 + 2 L 1 ∥ ∇ ϕ ( y ) ∥ 2 = ϕ ( y ) − 2 L 1 ∥ ∇ f ( y ) − ∇ f ( x ) ∥ 2
で、書き直すと (c)。(c⇒d) (c) を ( x , y ) (x, y) ( x , y ) と ( y , x ) (y, x) ( y , x ) に ついて 足す。(d⇒a) コーシー–シュワルツの 不等式より 1 L ∥ ∇ f ( x ) − ∇ f ( y ) ∥ 2 ≤ ∥ ∇ f ( x ) − ∇ f ( y ) ∥ ∥ x − y ∥ \frac{1}{L}\lVert \nabla f(x) - \nabla f(y) \rVert^2 \leq \lVert \nabla f(x) - \nabla f(y) \rVert\lVert x - y \rVert L 1 ∥ ∇ f ( x ) − ∇ f ( y ) ∥ 2 ≤ ∥ ∇ f ( x ) − ∇ f ( y )∥ ∥ x − y ∥ 。最後の 主張は (1) と 定理 2.15 に よる。 □ \square □
凸でない f f f では、(2.2) は L 2 ∥ x ∥ 2 − f \frac{L}{2}\lVert x \rVert^2 - f 2 L ∥ x ∥ 2 − f の 凸性( C 2 C^2 C 2 級なら ∇ 2 f ⪯ L I \nabla^2 f \preceq LI ∇ 2 f ⪯ L I )と 同値で(定理 2.14)、 L L L -平滑性より 弱い( − x 2 -x^2 − x 2 は (2.2) を 任意の L > 0 L > 0 L > 0 で 満たすが、導関数の リプシッツ定数は 2 2 2 )。「∇ 2 f ⪯ L I \nabla^2 f \preceq LI ∇ 2 f ⪯ L I なら L L L -平滑」は、f f f が 凸の とき(または ∇ 2 f ⪰ − L I \nabla^2 f \succeq -LI ∇ 2 f ⪰ − L I の とき)に 限る。
系 2.24 f : R n → R f\colon \mathbb{R}^n \to \mathbb{R} f : R n → R を μ \mu μ -強凸かつ L L L -平滑と する。この とき f f f は ただ 一つの 最小点 x ∗ x^{\ast} x ∗ を もち、 μ ≤ L \mu \leq L μ ≤ L で、すべての x x x に ついて
μ 2 ∥ x − x ∗ ∥ 2 ≤ f ( x ) − f ( x ∗ ) ≤ L 2 ∥ x − x ∗ ∥ 2 , 1 2 L ∥ ∇ f ( x ) ∥ 2 ≤ f ( x ) − f ( x ∗ ) ≤ 1 2 μ ∥ ∇ f ( x ) ∥ 2 \frac{\mu}{2}\lVert x - x^{\ast} \rVert^2 \leq f(x) - f(x^{\ast}) \leq \frac{L}{2}\lVert x - x^{\ast} \rVert^2, \qquad \frac{1}{2L}\lVert \nabla f(x) \rVert^2 \leq f(x) - f(x^{\ast}) \leq \frac{1}{2\mu}\lVert \nabla f(x) \rVert^2 2 μ ∥ x − x ∗ ∥ 2 ≤ f ( x ) − f ( x ∗ ) ≤ 2 L ∥ x − x ∗ ∥ 2 , 2 L 1 ∥ ∇ f ( x ) ∥ 2 ≤ f ( x ) − f ( x ∗ ) ≤ 2 μ 1 ∥ ∇ f ( x ) ∥ 2
証明. 定理 2.22 の 2 で x = 0 x = 0 x = 0 と すると f f f は 強圧的と わかり、最小点を もつ(第1章 定理 1.11)。 f f f は 凸関数と 狭義凸関数 μ 2 ∥ x ∥ 2 \frac{\mu}{2}\lVert x \rVert^2 2 μ ∥ x ∥ 2 の 和で 狭義凸だから、最小点は 一つ(定理 2.18)。 ∇ f ( x ∗ ) = 0 \nabla f(x^{\ast}) = 0 ∇ f ( x ∗ ) = 0 を 定理 2.22 の 2 と (2.2) に 代入すると 最初の 不等式を 得て、 μ ≤ L \mu \leq L μ ≤ L も 従う。定理 2.22 の 2 の 右辺は y = x − ∇ f ( x ) / μ y = x - \nabla f(x)/\mu y = x − ∇ f ( x ) / μ で 最小値 f ( x ) − 1 2 μ ∥ ∇ f ( x ) ∥ 2 f(x) - \frac{1}{2\mu}\lVert \nabla f(x) \rVert^2 f ( x ) − 2 μ 1 ∥ ∇ f ( x ) ∥ 2 を とるので、 f ( x ∗ ) ≥ f ( x ) − 1 2 μ ∥ ∇ f ( x ) ∥ 2 f(x^{\ast}) \geq f(x) - \frac{1}{2\mu}\lVert \nabla f(x) \rVert^2 f ( x ∗ ) ≥ f ( x ) − 2 μ 1 ∥ ∇ f ( x ) ∥ 2 。左端の 不等式は、定理 2.23 (c) で x x x を x ∗ x^{\ast} x ∗ に、y y y を x x x に 置き換えて 得られる。 □ \square □
右端の 不等式は ポリャク–ロヤシェヴィチの 不等式と 呼ばれる。 κ = L / μ \kappa = L/\mu κ = L / μ (≥ 1 \geq 1 ≥ 1 )を 条件数 (condition number) と いう。 第5章 では、これらの 不等式から 勾配法の 収束の 速さを 導き、 κ \kappa κ が 大きい ほど 遅い ことを 見る。
例 2.25 (条件数の 計算)(1) 1 2 x ⊤ Q x − c ⊤ x \frac{1}{2}x^{\top}Qx - c^{\top}x 2 1 x ⊤ Q x − c ⊤ x (Q ≻ 0 Q \succ 0 Q ≻ 0 )では μ = λ min ( Q ) \mu = \lambda_{\min}(Q) μ = λ m i n ( Q ) , L = λ max ( Q ) L = \lambda_{\max}(Q) L = λ m a x ( Q ) で、等高線は 軸の 比が κ \sqrt{\kappa} κ の 楕円である。(2) x x x 座標 0 , 1 , 2 0, 1, 2 0 , 1 , 2 の データに 切片つきの 直線を あてはめる 最小二乗問題 1 2 ∥ A x − b ∥ 2 \frac{1}{2}\lVert Ax - b \rVert^2 2 1 ∥ A x − b ∥ 2 では
A = ( 1 0 1 1 1 2 ) , ∇ 2 f = A ⊤ A = ( 3 3 3 5 ) A = \begin{pmatrix} 1 & 0 \\ 1 & 1 \\ 1 & 2 \end{pmatrix}, \qquad \nabla^2 f = A^{\top}A = \begin{pmatrix} 3 & 3 \\ 3 & 5 \end{pmatrix} A = 1 1 1 0 1 2 , ∇ 2 f = A ⊤ A = ( 3 3 3 5 )
で、固有値は 4 ± 10 4 \pm \sqrt{10} 4 ± 10 、κ = ( 13 + 4 10 ) / 3 ≈ 8.55 \kappa = (13 + 4\sqrt{10})/3 \approx 8.55 κ = ( 13 + 4 10 ) /3 ≈ 8.55 。リッジ正則化 λ 2 ∥ x ∥ 2 \frac{\lambda}{2}\lVert x \rVert^2 2 λ ∥ x ∥ 2 を 加えると 固有値が すべて λ \lambda λ 増え、λ = 1 \lambda = 1 λ = 1 で κ = ( 7 + 2 10 ) / 3 ≈ 4.44 \kappa = (7 + 2\sqrt{10})/3 \approx 4.44 κ = ( 7 + 2 10 ) /3 ≈ 4.44 。(3) ロジスティック回帰の 損失は、 σ ( 1 − σ ) ≤ 1 4 \sigma(1 - \sigma) \leq \frac{1}{4} σ ( 1 − σ ) ≤ 4 1 より ∇ 2 ℓ ⪯ 1 4 A ⊤ A \nabla^2\ell \preceq \frac{1}{4}A^{\top}A ∇ 2 ℓ ⪯ 4 1 A ⊤ A (A A A は a i ⊤ a_i^{\top} a i ⊤ を 行と する 行列)で 1 4 λ max ( A ⊤ A ) \frac{1}{4}\lambda_{\max}(A^{\top}A) 4 1 λ m a x ( A ⊤ A ) -平滑だが、∣ t ∣ → ∞ \lvert t \rvert \to \infty ∣ t ∣ → ∞ で σ ( t ) ( 1 − σ ( t ) ) → 0 \sigma(t)(1 - \sigma(t)) \to 0 σ ( t ) ( 1 − σ ( t )) → 0 なので 一般には 強凸でない。 λ 2 ∥ x ∥ 2 \frac{\lambda}{2}\lVert x \rVert^2 2 λ ∥ x ∥ 2 を 加えると λ \lambda λ -強凸に なる。
ヒント
実務では
条件数は 特徴量の スケールで 大きく 変わる。例 2.25 (2) で x x x 座標を 100 , 101 , 102 100, 101, 102 100 , 101 , 102 に すると κ ≈ 1.6 × 10 8 \kappa \approx 1.6 \times 10^{8} κ ≈ 1.6 × 1 0 8 だが、平均を 引いて − 1 , 0 , 1 -1, 0, 1 − 1 , 0 , 1 に すると A ⊤ A = diag ( 3 , 2 ) A^{\top}A = \operatorname{diag}(3, 2) A ⊤ A = diag ( 3 , 2 ) , κ = 1.5 \kappa = 1.5 κ = 1.5 に なる( A A A の 列空間は 同じなので、あてはめた 直線も 同じ)。勾配法の 反復回数は 条件数に 強く 依存する(第5章)ので、特徴量の 中心化・標準化は 基本的な 前処理である。リッジ正則化も 条件数を 改善するが、解その ものを 変える。
2.8 劣勾配
ラッソの ∥ x ∥ 1 \lVert x \rVert_1 ∥ x ∥ 1 や ヒンジ損失は 折れ目で 微分できない。それでも 凸関数には、勾配の 代わりに なる「下から 支える 一次関数」が ある。
定義 2.26 (劣勾配)C C C を 凸集合、 f : C → R f\colon C \to \mathbb{R} f : C → R を 凸関数、 x ∈ C x \in C x ∈ C と する。すべての y ∈ C y \in C y ∈ C で f ( y ) ≥ f ( x ) + ⟨ g , y − x ⟩ f(y) \geq f(x) + \langle g, y - x \rangle f ( y ) ≥ f ( x ) + ⟨ g , y − x ⟩ を 満たす g g g を x x x に おける 劣勾配 (subgradient) と いい、その全体 ∂ f ( x ) \partial f(x) ∂ f ( x ) (閉凸集合である)を 劣微分 (subdifferential) と いう。
定理 2.27 U U U を 開凸集合、 f : U → R f\colon U \to \mathbb{R} f : U → R を 凸関数と する。
すべての x ∈ U x \in U x ∈ U で ∂ f ( x ) ≠ ∅ \partial f(x) \neq \emptyset ∂ f ( x ) = ∅ である。
f f f が x x x で 微分可能なら、 ∂ f ( x ) = { ∇ f ( x ) } \partial f(x) = \lbrace \nabla f(x) \rbrace ∂ f ( x ) = { ∇ f ( x )} である。
証明. (1) epi f \operatorname{epi} f epi f は 凸で、 ( x , f ( x ) − ε ) ∉ epi f (x, f(x) - \varepsilon) \notin \operatorname{epi} f ( x , f ( x ) − ε ) ∈ / epi f より ( x , f ( x ) ) (x, f(x)) ( x , f ( x )) は その 内点でない。定理 2.10 より ( a , β ) ≠ ( 0 , 0 ) (a, \beta) \neq (0, 0) ( a , β ) = ( 0 , 0 ) が あって、すべての ( y , s ) ∈ epi f (y, s) \in \operatorname{epi} f ( y , s ) ∈ epi f で a ⊤ y + β s ≤ a ⊤ x + β f ( x ) a^{\top}y + \beta s \leq a^{\top}x + \beta f(x) a ⊤ y + β s ≤ a ⊤ x + β f ( x ) 。s → ∞ s \to \infty s → ∞ と して β ≤ 0 \beta \leq 0 β ≤ 0 。β = 0 \beta = 0 β = 0 なら a ⊤ ( y − x ) ≤ 0 a^{\top}(y - x) \leq 0 a ⊤ ( y − x ) ≤ 0 (∀ y ∈ U \forall y \in U ∀ y ∈ U )で、y = x + ε a y = x + \varepsilon a y = x + ε a と すると a = 0 a = 0 a = 0 と なり矛盾。よって β < 0 \beta < 0 β < 0 で、s = f ( y ) s = f(y) s = f ( y ) と して ∣ β ∣ \lvert \beta \rvert ∣ β ∣ で 割ると a / ∣ β ∣ ∈ ∂ f ( x ) a/\lvert \beta \rvert \in \partial f(x) a / ∣ β ∣ ∈ ∂ f ( x ) 。(x x x が 定義域の 内点である ことは 外せない。開でない 区間 [ 0 , ∞ ) [0, \infty) [ 0 , ∞ ) 上の 凸関数 − x -\sqrt{x} − x は、端点 0 0 0 で 劣勾配を もたない。)(2) 定理 2.14 の (1⇒2) の 証明は x x x での 微分 可能性しか 使わないので ∇ f ( x ) ∈ ∂ f ( x ) \nabla f(x) \in \partial f(x) ∇ f ( x ) ∈ ∂ f ( x ) 。逆に g ∈ ∂ f ( x ) g \in \partial f(x) g ∈ ∂ f ( x ) なら、任意の d d d と t > 0 t > 0 t > 0 で f ( x + t d ) − f ( x ) ≥ t ⟨ g , d ⟩ f(x + td) - f(x) \geq t\langle g, d \rangle f ( x + t d ) − f ( x ) ≥ t ⟨ g , d ⟩ 。t t t で 割って t → + 0 t \to +0 t → + 0 と すると ⟨ ∇ f ( x ) − g , d ⟩ ≥ 0 \langle \nabla f(x) - g, d \rangle \geq 0 ⟨ ∇ f ( x ) − g , d ⟩ ≥ 0 で、d = g − ∇ f ( x ) d = g - \nabla f(x) d = g − ∇ f ( x ) と すれば g = ∇ f ( x ) g = \nabla f(x) g = ∇ f ( x ) 。□ \square □
例 2.28 (劣微分の 計算)(1) ∣ x ∣ \lvert x \rvert ∣ x ∣ :x ≠ 0 x \neq 0 x = 0 で { sign x } \lbrace \operatorname{sign} x \rbrace { sign x } 、x = 0 x = 0 x = 0 で [ − 1 , 1 ] [-1, 1] [ − 1 , 1 ] 。(2) ∥ x ∥ \lVert x \rVert ∥ x ∥ :x ≠ 0 x \neq 0 x = 0 で { x / ∥ x ∥ } \lbrace x/\lVert x \rVert \rbrace { x / ∥ x ∥} 、x = 0 x = 0 x = 0 で 閉単位球(コーシー–シュワルツの 不等式)。(3) ∥ x ∥ 1 \lVert x \rVert_1 ∥ x ∥ 1 :∥ y ∥ 1 − ∥ x ∥ 1 − ⟨ g , y − x ⟩ = ∑ i ( ∣ y i ∣ − ∣ x i ∣ − g i ( y i − x i ) ) \lVert y \rVert_1 - \lVert x \rVert_1 - \langle g, y - x \rangle = \sum_i(\lvert y_i \rvert - \lvert x_i \rvert - g_i(y_i - x_i)) ∥ y ∥ 1 − ∥ x ∥ 1 − ⟨ g , y − x ⟩ = ∑ i (∣ y i ∣ − ∣ x i ∣ − g i ( y i − x i )) と 成分ごとに 分かれるので、 g ∈ ∂ ∥ x ∥ 1 g \in \partial\lVert x \rVert_1 g ∈ ∂ ∥ x ∥ 1 は「x i ≠ 0 x_i \neq 0 x i = 0 なら g i = sign x i g_i = \operatorname{sign} x_i g i = sign x i 、x i = 0 x_i = 0 x i = 0 なら g i ∈ [ − 1 , 1 ] g_i \in [-1, 1] g i ∈ [ − 1 , 1 ] 」と 同値である。
定理 2.29 (劣勾配に よる 最適性条件) C C C を 凸集合、 f : C → R f\colon C \to \mathbb{R} f : C → R を 凸関数と する。 x ∗ ∈ C x^{\ast} \in C x ∗ ∈ C が 最小点である ための 必要十分条件は 0 ∈ ∂ f ( x ∗ ) 0 \in \partial f(x^{\ast}) 0 ∈ ∂ f ( x ∗ ) である。
証明. 0 ∈ ∂ f ( x ∗ ) 0 \in \partial f(x^{\ast}) 0 ∈ ∂ f ( x ∗ ) は、定義に より「すべての y ∈ C y \in C y ∈ C で f ( y ) ≥ f ( x ∗ ) f(y) \geq f(x^{\ast}) f ( y ) ≥ f ( x ∗ ) 」と 同じ ことである。 □ \square □
定理は 定義の 言い換えだが、劣微分が 計算できれば 強力である。劣勾配の 不等式を 足すと ∂ f 1 ( x ) + ∂ f 2 ( x ) ⊂ ∂ ( f 1 + f 2 ) ( x ) \partial f_1(x) + \partial f_2(x) \subset \partial(f_1 + f_2)(x) ∂ f 1 ( x ) + ∂ f 2 ( x ) ⊂ ∂ ( f 1 + f 2 ) ( x ) が わかる( f 1 , f 2 f_1, f_2 f 1 , f 2 が R n \mathbb{R}^n R n 上の 凸関数なら 等号が 成り立つ。モロー–ロッカフェラーの 定理、主張のみ。Nesterov, Lectures on Convex Optimization などを 参照)。
例 2.30 (ソフト閾値)λ > 0 \lambda > 0 λ > 0 とし、f ( x ) = 1 2 ( x − a ) 2 + λ ∣ x ∣ f(x) = \frac{1}{2}(x - a)^2 + \lambda\lvert x \rvert f ( x ) = 2 1 ( x − a ) 2 + λ ∣ x ∣ を 最小化する。 x − a + λ s = 0 x - a + \lambda s = 0 x − a + λ s = 0 と なる s ∈ ∂ ∣ ⋅ ∣ ( x ) s \in \partial\lvert \cdot \rvert(x) s ∈ ∂ ∣ ⋅ ∣ ( x ) が あれば、上の 包含より 0 ∈ ∂ f ( x ) 0 \in \partial f(x) 0 ∈ ∂ f ( x ) で x x x は 最小点である( f f f は 狭義凸なので 一つ)。 ∣ a ∣ > λ \lvert a \rvert > \lambda ∣ a ∣ > λ なら x = a − λ sign a x = a - \lambda\operatorname{sign} a x = a − λ sign a (s = sign a s = \operatorname{sign} a s = sign a )、∣ a ∣ ≤ λ \lvert a \rvert \leq \lambda ∣ a ∣ ≤ λ なら x = 0 x = 0 x = 0 (s = a / λ s = a/\lambda s = a / λ )が これを 満たすので、最小点は
S λ ( a ) = sign ( a ) max ( ∣ a ∣ − λ , 0 ) S_{\lambda}(a) = \operatorname{sign}(a)\max(\lvert a \rvert - \lambda, 0) S λ ( a ) = sign ( a ) max (∣ a ∣ − λ , 0 )
である(ソフト閾値関数 , soft-thresholding。S 2 ( 3 ) = 1 S_2(3) = 1 S 2 ( 3 ) = 1 , S 2 ( 1 ) = 0 S_2(1) = 0 S 2 ( 1 ) = 0 )。∣ a ∣ ≤ λ \lvert a \rvert \leq \lambda ∣ a ∣ ≤ λ で 解が ちょうど 0 0 0 に なることが、ラッソが 0 0 0 の 多い(スパースな)解を 与える 理由であり、第5章の 近接勾配法で 使われる。
注意
g ∈ ∂ f ( x ) g \in \partial f(x) g ∈ ∂ f ( x ) でも − g -g − g は 降下方向とは 限らない。 f ( x ) = ∣ x 1 ∣ + 2 ∣ x 2 ∣ f(x) = \lvert x_1 \rvert + 2\lvert x_2 \rvert f ( x ) = ∣ x 1 ∣ + 2 ∣ x 2 ∣ , x = ( 1 , 0 ) x = (1, 0) x = ( 1 , 0 ) では g = ( 1 , 2 ) ∈ ∂ f ( x ) g = (1, 2) \in \partial f(x) g = ( 1 , 2 ) ∈ ∂ f ( x ) (∣ y 1 ∣ + 2 ∣ y 2 ∣ ≥ y 1 + 2 y 2 \lvert y_1 \rvert + 2\lvert y_2 \rvert \geq y_1 + 2y_2 ∣ y 1 ∣ + 2 ∣ y 2 ∣ ≥ y 1 + 2 y 2 に よる)だが、 0 < s < 1 0 < s < 1 0 < s < 1 で f ( x − s g ) = 1 + 3 s > f ( x ) f(x - sg) = 1 + 3s > f(x) f ( x − s g ) = 1 + 3 s > f ( x ) 。劣勾配を 使う 反復法では 値は 単調に 減るとは 限らない。
2.9 共役関数(紹介)
空でない C C C 上の 関数 f f f に 対し、 f ∗ ( y ) = sup x ∈ C ( ⟨ y , x ⟩ − f ( x ) ) ∈ ( − ∞ , + ∞ ] f^{\ast}(y) = \sup_{x \in C}(\langle y, x \rangle - f(x)) \in (-\infty, +\infty] f ∗ ( y ) = sup x ∈ C (⟨ y , x ⟩ − f ( x )) ∈ ( − ∞ , + ∞ ] を 共役関数 (conjugate function) と いう。 f ∗ f^{\ast} f ∗ は アフィン関数の 上限なので(値 + ∞ +\infty + ∞ を 許す 意味で)凸であり、 フェンシェル–ヤングの 不等式 f ( x ) + f ∗ ( y ) ≥ ⟨ x , y ⟩ f(x) + f^{\ast}(y) \geq \langle x, y \rangle f ( x ) + f ∗ ( y ) ≥ ⟨ x , y ⟩ が 成り立つ( f f f が 凸なら、等号は y ∈ ∂ f ( x ) y \in \partial f(x) y ∈ ∂ f ( x ) の ときに 限る)。たとえば Q ≻ 0 Q \succ 0 Q ≻ 0 なら 1 2 x ⊤ Q x \frac{1}{2}x^{\top}Qx 2 1 x ⊤ Q x の 共役は 1 2 y ⊤ Q − 1 y \frac{1}{2}y^{\top}Q^{-1}y 2 1 y ⊤ Q − 1 y 、e x e^{x} e x の 共役は y > 0 y > 0 y > 0 で y log y − y y\log y - y y log y − y 、y = 0 y = 0 y = 0 で 0 0 0 、y < 0 y < 0 y < 0 で + ∞ +\infty + ∞ である。凸で エピグラフが 閉じた f f f では、C C C の 外で f = + ∞ f = +\infty f = + ∞ と みなすと f ∗ ∗ = f f^{\ast\ast} = f f ∗∗ = f が 成り立つ(フェンシェル–モローの 定理、主張のみ。Boyd–Vandenberghe, Convex Optimization 3.3 節)。共役関数は ラグランジュ双対(第4章)に 現れる。たとえば A x = b Ax = b A x = b のもとで f ( x ) f(x) f ( x ) を 最小化する とき、 f ( x ) + ν ⊤ ( A x − b ) f(x) + \nu^{\top}(Ax - b) f ( x ) + ν ⊤ ( A x − b ) の x x x に ついての 下限は − b ⊤ ν − f ∗ ( − A ⊤ ν ) -b^{\top}\nu - f^{\ast}(-A^{\top}\nu) − b ⊤ ν − f ∗ ( − A ⊤ ν ) である。
まとめ
凸集合の 共通部分・アフィン像・逆像は 凸。多面体・球・確率単体・半正定値行列の 錐は 凸である。
閉凸集合への 射影は 存在して 一意で、 ⟨ x − p , y − p ⟩ ≤ 0 \langle x - p, y - p \rangle \leq 0 ⟨ x − p , y − p ⟩ ≤ 0 で 特徴づけられ、非拡大である。
点と 閉凸集合、閉凸集合と コンパクト凸集合は 強分離できる(コンパクト性は 外せない)。有限次元では 交わらない 凸集合は 超平面で 分離でき(一方が 開なら 開集合の 側は 狭義の 不等式)、内点でない 点には 支持超 平面が ある。
開凸集合上の 微分 可能な f f f が 凸 ⟺ \iff ⟺ f ( y ) ≥ f ( x ) + ⟨ ∇ f ( x ) , y − x ⟩ f(y) \geq f(x) + \langle \nabla f(x), y - x \rangle f ( y ) ≥ f ( x ) + ⟨ ∇ f ( x ) , y − x ⟩ ⟺ \iff ⟺ 勾配が 単調、 C 2 C^2 C 2 級なら ⟺ ∇ 2 f ⪰ 0 \iff \nabla^2 f \succeq 0 ⟺ ∇ 2 f ⪰ 0 。定義域が 開である ことは 外せない。
凸関数の 局所最小は 大域最小で、微分 可能なら ∇ f ( x ∗ ) = 0 \nabla f(x^{\ast}) = 0 ∇ f ( x ∗ ) = 0 が 最小の 必要十分条件(存在は 別問題)。凸性は 凸性を 保つ演算で 確かめる。
μ \mu μ -強凸 ⟺ ∇ 2 f ⪰ μ I \iff \nabla^2 f \succeq \mu I ⟺ ∇ 2 f ⪰ μ I 、凸で L L L -平滑 ⟺ 0 ⪯ ∇ 2 f ⪯ L I \iff 0 \preceq \nabla^2 f \preceq LI ⟺ 0 ⪯ ∇ 2 f ⪯ L I (C 2 C^2 C 2 級の 場合)。条件数 L / μ L/\mu L / μ は スケーリングや 正則化で 変わる。
開凸集合上の 凸関数は 劣勾配を もち、 0 ∈ ∂ f ( x ∗ ) 0 \in \partial f(x^{\ast}) 0 ∈ ∂ f ( x ∗ ) が 最適性条件。 1 2 ( x − a ) 2 + λ ∣ x ∣ \frac{1}{2}(x - a)^2 + \lambda\lvert x \rvert 2 1 ( x − a ) 2 + λ ∣ x ∣ の 最小点は ソフト閾値 S λ ( a ) S_{\lambda}(a) S λ ( a ) である。
演習問題
問題 2.1 ★ 次の 集合は 凸か。(1) { x ∈ R 2 ∣ x 1 > 0 , x 1 x 2 ≥ 1 } \lbrace x \in \mathbb{R}^2 \mid x_1 > 0, \ x_1x_2 \geq 1 \rbrace { x ∈ R 2 ∣ x 1 > 0 , x 1 x 2 ≥ 1 } (2) { x ∈ R n ∣ ∥ x − a ∥ ≤ ∥ x − b ∥ } \lbrace x \in \mathbb{R}^n \mid \lVert x - a \rVert \leq \lVert x - b \rVert \rbrace { x ∈ R n ∣ ∥ x − a ∥ ≤ ∥ x − b ∥} (a ≠ b a \neq b a = b ) (3) { x ∈ R 2 ∣ x 1 2 − x 2 2 ≤ 1 } \lbrace x \in \mathbb{R}^2 \mid x_1^2 - x_2^2 \leq 1 \rbrace { x ∈ R 2 ∣ x 1 2 − x 2 2 ≤ 1 } (4) { x ∈ R n ∣ x ⊤ P x ≤ 1 } \lbrace x \in \mathbb{R}^n \mid x^{\top}Px \leq 1 \rbrace { x ∈ R n ∣ x ⊤ P x ≤ 1 } (P ⪰ 0 P \succeq 0 P ⪰ 0 )
解答
(1) 凸。x 1 > 0 x_1 > 0 x 1 > 0 なら x 1 x 2 ≥ 1 ⟺ x 2 ≥ 1 / x 1 x_1x_2 \geq 1 \iff x_2 \geq 1/x_1 x 1 x 2 ≥ 1 ⟺ x 2 ≥ 1/ x 1 なので、凸関数 1 / x 1/x 1/ x の エピグラフである(命題 2.13)。(2) 凸。両辺を 2 乗して 展開すると 閉半空間 2 ( b − a ) ⊤ x ≤ ∥ b ∥ 2 − ∥ a ∥ 2 2(b - a)^{\top}x \leq \lVert b \rVert^2 - \lVert a \rVert^2 2 ( b − a ) ⊤ x ≤ ∥ b ∥ 2 − ∥ a ∥ 2 に なる。(3) 凸でない。 ( 1.5 , ± 1.2 ) (1.5, \pm 1.2) ( 1.5 , ± 1.2 ) は 属する( 2.25 − 1.44 ≤ 1 2.25 - 1.44 \leq 1 2.25 − 1.44 ≤ 1 )が、中点 ( 1.5 , 0 ) (1.5, 0) ( 1.5 , 0 ) は 属さない。(4) 凸。凸関数 x ⊤ P x x^{\top}Px x ⊤ P x (例 2.17)の 下位集合である。
問題 2.2 ★ 次の 関数は 凸か。(1) f ( x , y ) = x 2 / y f(x, y) = x^2/y f ( x , y ) = x 2 / y (y > 0 y > 0 y > 0 ) (2) f ( x ) = max i x i − min i x i f(x) = \max_i x_i - \min_i x_i f ( x ) = max i x i − min i x i (3) f ( x , y ) = x y f(x, y) = xy f ( x , y ) = x y (4) f ( x ) = ∥ A x − b ∥ 1 + λ ∥ x ∥ 2 f(x) = \lVert Ax - b \rVert_1 + \lambda\lVert x \rVert^2 f ( x ) = ∥ A x − b ∥ 1 + λ ∥ x ∥ 2 (λ ≥ 0 \lambda \geq 0 λ ≥ 0 )
解答
(1) 凸。開凸集合 { y > 0 } \lbrace y > 0 \rbrace { y > 0 } 上で
∇ 2 f = ( 2 / y − 2 x / y 2 − 2 x / y 2 2 x 2 / y 3 ) = 2 y 3 ( y − x ) ( y − x ) ⪰ 0 \nabla^2 f = \begin{pmatrix} 2/y & -2x/y^2 \\ -2x/y^2 & 2x^2/y^3 \end{pmatrix} = \frac{2}{y^3}\begin{pmatrix} y \\ -x \end{pmatrix}\begin{pmatrix} y & -x \end{pmatrix} \succeq 0 ∇ 2 f = ( 2/ y − 2 x / y 2 − 2 x / y 2 2 x 2 / y 3 ) = y 3 2 ( y − x ) ( y − x ) ⪰ 0
なので 凸(定理 2.15)。(2) 凸。 max i x i \max_i x_i max i x i と − min i x i = max i ( − x i ) -\min_i x_i = \max_i(-x_i) − min i x i = max i ( − x i ) は 一次関数の 最大(命題 2.20 の 3・1)。(3) 凸でない( f ( t , − t ) = − t 2 f(t, -t) = -t^2 f ( t , − t ) = − t 2 )。(4) 凸(命題 2.20 の 1・2)。
問題 2.3 ★ ★ f ( x ) = log ∑ i = 1 n e x i f(x) = \log\sum_{i=1}^{n}e^{x_i} f ( x ) = log ∑ i = 1 n e x i に ついて ∇ 2 f ( x ) = diag ( p ) − p p ⊤ \nabla^2 f(x) = \operatorname{diag}(p) - pp^{\top} ∇ 2 f ( x ) = diag ( p ) − p p ⊤ (p i = e x i / ∑ j e x j p_i = e^{x_i}/\sum_j e^{x_j} p i = e x i / ∑ j e x j )を 確かめ、 0 ⪯ ∇ 2 f ( x ) ⪯ I 0 \preceq \nabla^2 f(x) \preceq I 0 ⪯ ∇ 2 f ( x ) ⪯ I を 示せ。これから f f f が 凸かつ 1 1 1 -平滑である ことを 結論し、 f f f が 強凸でない ことを 示せ。
解答
∂ i f = p i \partial_if = p_i ∂ i f = p i , ∂ j p i = p i δ i j − p i p j \partial_jp_i = p_i\delta_{ij} - p_ip_j ∂ j p i = p i δ ij − p i p j より ∇ 2 f = diag ( p ) − p p ⊤ \nabla^2 f = \operatorname{diag}(p) - pp^{\top} ∇ 2 f = diag ( p ) − p p ⊤ 。p i > 0 p_i > 0 p i > 0 , ∑ i p i = 1 \sum_i p_i = 1 ∑ i p i = 1 なので、コーシー–シュワルツの 不等式より ( ∑ i p i v i ) 2 = ( ∑ i p i ⋅ p i v i ) 2 ≤ ∑ i p i v i 2 (\sum_i p_iv_i)^2 = (\sum_i \sqrt{p_i} \cdot \sqrt{p_i}v_i)^2 \leq \sum_i p_iv_i^2 ( ∑ i p i v i ) 2 = ( ∑ i p i ⋅ p i v i ) 2 ≤ ∑ i p i v i 2 で、
0 ≤ v ⊤ ∇ 2 f ( x ) v = ∑ i p i v i 2 − ( ∑ i p i v i ) 2 ≤ ∑ i p i v i 2 ≤ max i v i 2 ≤ ∥ v ∥ 2 0 \leq v^{\top}\nabla^2 f(x)v = \sum_i p_iv_i^2 - \Bigl(\sum_i p_iv_i\Bigr)^2 \leq \sum_i p_iv_i^2 \leq \max_i v_i^2 \leq \lVert v \rVert^2 0 ≤ v ⊤ ∇ 2 f ( x ) v = i ∑ p i v i 2 − ( i ∑ p i v i ) 2 ≤ i ∑ p i v i 2 ≤ i max v i 2 ≤ ∥ v ∥ 2
よって 0 ⪯ ∇ 2 f ⪯ I 0 \preceq \nabla^2 f \preceq I 0 ⪯ ∇ 2 f ⪯ I で、f f f は 凸(定理 2.15)かつ 1 1 1 -平滑(定理 2.23 (2))。1 = ( 1 , … , 1 ) ⊤ \mathbf{1} = (1, \dots, 1)^{\top} 1 = ( 1 , … , 1 ) ⊤ に ついて ∇ 2 f ( x ) 1 = p − p ( p ⊤ 1 ) = 0 \nabla^2 f(x)\mathbf{1} = p - p(p^{\top}\mathbf{1}) = 0 ∇ 2 f ( x ) 1 = p − p ( p ⊤ 1 ) = 0 なので、∇ 2 f ⪰ μ I \nabla^2 f \succeq \mu I ∇ 2 f ⪰ μ I と なる μ > 0 \mu > 0 μ > 0 は なく、強凸でない(定理 2.22)。実際 f ( x + t 1 ) = f ( x ) + t f(x + t\mathbf{1}) = f(x) + t f ( x + t 1 ) = f ( x ) + t は t t t の 一次関数である。
問題 2.4 ★ ★ 次の 主張の 誤りを、反例を 挙げて 指摘せよ。(a)「損失関数の ヘッセ行列を 初期点で 計算したら 半正定値だったので、この 損失関数は 凸である」。(b)「 f f f は (2.2) を L = 1 L = 1 L = 1 で 満たすので、 ∇ f \nabla f ∇ f は 1-リプシッツである」。(c)「ロジスティック回帰の 損失は 凸なので、勾配法で 必ず 最小点に 収束する」。
解答
(a) 定理 2.15 は すべての 点 で ∇ 2 f ⪰ 0 \nabla^2 f \succeq 0 ∇ 2 f ⪰ 0 を 要求する。 f ( x ) = x 4 − x 2 f(x) = x^4 - x^2 f ( x ) = x 4 − x 2 は f ′ ′ ( 1 ) = 10 > 0 f''(1) = 10 > 0 f ′′ ( 1 ) = 10 > 0 だが f ′ ′ ( 0 ) = − 2 < 0 f''(0) = -2 < 0 f ′′ ( 0 ) = − 2 < 0 で、凸でない。(b) 凸でない 関数では (2.2) から リプシッツ性は 出ない(定理 2.23 の 後の 注意)。 f ( x ) = − x 2 f(x) = -x^2 f ( x ) = − x 2 は (2.2) を L = 1 L = 1 L = 1 で 満たすが、 f ′ ( x ) = − 2 x f'(x) = -2x f ′ ( x ) = − 2 x は 1-リプシッツでない。(c) 凸性は 最小点の 存在を 保証しない。完全に 分離できる データでは 最小点が なく(第1章 例 1.13)、反復は 発散する。 λ 2 ∥ x ∥ 2 \frac{\lambda}{2}\lVert x \rVert^2 2 λ ∥ x ∥ 2 を 加えれば 強凸に なり、最小点が ただ 一つ 存在する(系 2.24)。
問題 2.5 ★ ★ (1) ヒンジ損失 h ( t ) = max ( 0 , 1 − t ) h(t) = \max(0, 1 - t) h ( t ) = max ( 0 , 1 − t ) の 劣微分を すべての t t t で 求めよ。(2) x = ( 2 , 0 , − 1 ) x = (2, 0, -1) x = ( 2 , 0 , − 1 ) に おける ∂ ∥ x ∥ 1 \partial\lVert x \rVert_1 ∂ ∥ x ∥ 1 を 求めよ。(3) a = ( 3 , − 0.5 , − 2 ) a = (3, -0.5, -2) a = ( 3 , − 0.5 , − 2 ) に ついて f ( x ) = 1 2 ∥ x − a ∥ 2 + ∥ x ∥ 1 f(x) = \frac{1}{2}\lVert x - a \rVert^2 + \lVert x \rVert_1 f ( x ) = 2 1 ∥ x − a ∥ 2 + ∥ x ∥ 1 の 最小点を 求め、 0 ∈ ∂ f ( x ∗ ) 0 \in \partial f(x^{\ast}) 0 ∈ ∂ f ( x ∗ ) を 確かめよ。
解答
(1) t < 1 t < 1 t < 1 で { − 1 } \lbrace -1 \rbrace { − 1 } 、t > 1 t > 1 t > 1 で { 0 } \lbrace 0 \rbrace { 0 } (定理 2.27 (2))。g ∈ ∂ h ( 1 ) g \in \partial h(1) g ∈ ∂ h ( 1 ) は「すべての s s s で max ( 0 , 1 − s ) ≥ g ( s − 1 ) \max(0, 1 - s) \geq g(s - 1) max ( 0 , 1 − s ) ≥ g ( s − 1 ) 」と 同値で、 s > 1 s > 1 s > 1 から g ≤ 0 g \leq 0 g ≤ 0 、s < 1 s < 1 s < 1 から g ≥ − 1 g \geq -1 g ≥ − 1 。よって ∂ h ( 1 ) = [ − 1 , 0 ] \partial h(1) = [-1, 0] ∂ h ( 1 ) = [ − 1 , 0 ] 。
(2) 例 2.28 (3) より { ( 1 , s , − 1 ) ∣ − 1 ≤ s ≤ 1 } \lbrace (1, s, -1) \mid -1 \leq s \leq 1 \rbrace {( 1 , s , − 1 ) ∣ − 1 ≤ s ≤ 1 } 。
(3) f f f は 成分ごとの 1 2 ( x i − a i ) 2 + ∣ x i ∣ \frac{1}{2}(x_i - a_i)^2 + \lvert x_i \rvert 2 1 ( x i − a i ) 2 + ∣ x i ∣ の 和なので、例 2.30 より x ∗ = ( S 1 ( 3 ) , S 1 ( − 0.5 ) , S 1 ( − 2 ) ) = ( 2 , 0 , − 1 ) x^{\ast} = (S_1(3), S_1(-0.5), S_1(-2)) = (2, 0, -1) x ∗ = ( S 1 ( 3 ) , S 1 ( − 0.5 ) , S 1 ( − 2 )) = ( 2 , 0 , − 1 ) 。確認:x ∗ − a = ( − 1 , 0.5 , 1 ) x^{\ast} - a = (-1, 0.5, 1) x ∗ − a = ( − 1 , 0.5 , 1 ) に (2) の s = − 0.5 s = -0.5 s = − 0.5 の 元 ( 1 , − 0.5 , − 1 ) (1, -0.5, -1) ( 1 , − 0.5 , − 1 ) を 足すと 0 0 0 なので、0 ∈ ( x ∗ − a ) + ∂ ∥ x ∗ ∥ 1 ⊂ ∂ f ( x ∗ ) 0 \in (x^{\ast} - a) + \partial\lVert x^{\ast} \rVert_1 \subset \partial f(x^{\ast}) 0 ∈ ( x ∗ − a ) + ∂ ∥ x ∗ ∥ 1 ⊂ ∂ f ( x ∗ ) で、定理 2.29 より x ∗ x^{\ast} x ∗ は 最小点である。
問題 2.6 ★ ★ ロジスティック損失 f ( x ) = log ( 1 + e x ) f(x) = \log(1 + e^{x}) f ( x ) = log ( 1 + e x ) の 共役関数 f ∗ f^{\ast} f ∗ を 求め、 0 ≤ y ≤ 1 0 \leq y \leq 1 0 ≤ y ≤ 1 で f ∗ ( y ) = y log y + ( 1 − y ) log ( 1 − y ) f^{\ast}(y) = y\log y + (1 - y)\log(1 - y) f ∗ ( y ) = y log y + ( 1 − y ) log ( 1 − y ) (0 log 0 = 0 0\log 0 = 0 0 log 0 = 0 と する)、それ以外で + ∞ +\infty + ∞ と なる ことを 示せ。
解答
f ∗ ( y ) = sup x ψ ( x ) f^{\ast}(y) = \sup_x \psi(x) f ∗ ( y ) = sup x ψ ( x ) , ψ ( x ) = y x − log ( 1 + e x ) \psi(x) = yx - \log(1 + e^{x}) ψ ( x ) = y x − log ( 1 + e x ) と おく。 ψ ′ ( x ) = y − σ ( x ) \psi'(x) = y - \sigma(x) ψ ′ ( x ) = y − σ ( x ) で、σ \sigma σ は R \mathbb{R} R から ( 0 , 1 ) (0, 1) ( 0 , 1 ) への 狭義単調増加な 全単射である。 0 < y < 1 0 < y < 1 0 < y < 1 なら ψ \psi ψ は 凹で、 σ ( x ) = y \sigma(x) = y σ ( x ) = y 、すな わち x = log y 1 − y x = \log\frac{y}{1 - y} x = log 1 − y y で 最大に なり、 log ( 1 + e x ) = log 1 1 − y \log(1 + e^{x}) = \log\frac{1}{1 - y} log ( 1 + e x ) = log 1 − y 1 より
f ∗ ( y ) = y log y 1 − y + log ( 1 − y ) = y log y + ( 1 − y ) log ( 1 − y ) f^{\ast}(y) = y\log\frac{y}{1 - y} + \log(1 - y) = y\log y + (1 - y)\log(1 - y) f ∗ ( y ) = y log 1 − y y + log ( 1 − y ) = y log y + ( 1 − y ) log ( 1 − y )
y = 0 y = 0 y = 0 なら ψ ( x ) = − log ( 1 + e x ) \psi(x) = -\log(1 + e^{x}) ψ ( x ) = − log ( 1 + e x ) は 負で、 x → − ∞ x \to -\infty x → − ∞ で 0 0 0 に 近づくので f ∗ ( 0 ) = 0 f^{\ast}(0) = 0 f ∗ ( 0 ) = 0 。y = 1 y = 1 y = 1 なら ψ ( x ) = − log ( 1 + e − x ) \psi(x) = -\log(1 + e^{-x}) ψ ( x ) = − log ( 1 + e − x ) で 同様に f ∗ ( 1 ) = 0 f^{\ast}(1) = 0 f ∗ ( 1 ) = 0 。y > 1 y > 1 y > 1 なら x ≥ 0 x \geq 0 x ≥ 0 で ψ ( x ) ≥ ( y − 1 ) x − log 2 → ∞ \psi(x) \geq (y - 1)x - \log 2 \to \infty ψ ( x ) ≥ ( y − 1 ) x − log 2 → ∞ 、y < 0 y < 0 y < 0 なら x ≤ 0 x \leq 0 x ≤ 0 で ψ ( x ) ≥ y x − log 2 → ∞ \psi(x) \geq yx - \log 2 \to \infty ψ ( x ) ≥ y x − log 2 → ∞ (x → − ∞ x \to -\infty x → − ∞ )なので f ∗ ( y ) = + ∞ f^{\ast}(y) = +\infty f ∗ ( y ) = + ∞ 。− f ∗ -f^{\ast} − f ∗ は ベルヌーイ分布の エントロピーであり、ロジスティック損失(交差エントロピー)と エントロピーが 共役で 結ばれている。