MML
  • 数学プログラミング置き場
  • 数学記号記法一覧
  • 機械学習のための数学
    • 確率論
  • 機械学習・データサイエンス
    • 1. 機械学習・データサイエンス(概要)
    • 2. 機械学習
      • 2.1. 機械学習をこれから始める人へ
      • 2.2. 回帰問題
        • 2.2.1. 線形回帰の基礎
        • 2.2.2. 線形回帰モデルと正則化
        • 2.2.3. 線形基底関数モデル
          • いろいろな基底関数
        • 2.2.4. ベイズ線形回帰モデル
        • 2.2.5. カーネル回帰
        • 2.2.6. ガウス過程回帰
        • 2.2.7. 一般化線形モデル
        • 2.2.8. 一般化加法モデル
        • 2.2.9. 分位点回帰
        • 2.2.10. 3層ニューラルネットワーク
        • 2.2.11. 深層学習
        • 2.2.12. 深層学習とガウス過程
      • 分類問題
        • 線形判別分析
        • SVM
        • カーネルSVM
      • 強化学習
    • 3. データサイエンス
      • 主成分分析
      • 白色化
      • 低ランク行列近似
      • 行列補完
      • NMF
      • スパースモデリング
        • LASSO
          • FISTA
          • ADMM
        • 全状態探索法
    • 4. ニューラルネットワーク
      • CNN
      • GAN
      • VAE
      • GNN
      • word2vec
    • 5. ベイズ推論
      • 確率
      • 確率分布
      • ちょっとだけ振り返って
      • MCMC
      • ギブスサンプリング
      • 変分推論
  • 連続最適化
    • 1. Taylor展開と関数近似
    • 2. Taylorの定理と誤差の評価
    • 3. 最適化問題
    • 4. 直線探索法
      • 4.1. 最急降下法
      • 4.2. Armijo条件・Wolfe条件
      • 4.3. ニュートン法
      • 4.4. 共役勾配法
        • 4.4.1. 原始的な共役勾配法
        • 4.4.2. PRP法
        • 4.4.3. HS法
      • 4.5. 準ニュートン法
        • 4.5.1. DFP法
        • 4.5.2. BFGS法
        • 4.5.3. L-BFGS法
        • 4.5.4. Bryoden family of methods
        • 4.5.5. Broyden法
    • 5. 確率的最適化
      • 5.1. 確率的勾配降下法(SGD)
      • 5.2. 確率的準ニュートン法(SQN)
  • 行列計算
    • QR分解
    • コレスキー分解
    • 一般化コレスキー分解
    • 固有値分解
    • 一般化固有値問題
    • 特異値分解
  • その他
    • 位相限定相関法
  • カーネル法
    • カーネル法(概要)
  • Riemannian最適化
    • Riemann多様体って?
    • Riemannian最適化(概要)
    • 概説
    • 詳説
      • 解析学
      • 線形代数学(反変性・共変性・不変量)
      • テンソル解析
      • リーマン積分と微分形式
      • 位相空間論
      • 多様体論(解析幾何学)
      • Riemann多様体
  • プログラミングサイドブック
    • はじめに
    • Step 0. 勉強の仕方
      • いつ始める?どれくらい勉強する?
      • どうやって勉強する?
      • どこから手をつけたらいい?
    • Step 1. インターネットで調べる
    • Step 2. コンピュータの仕組みを学ぶ
    • Step 3. 開発環境を整える
    • Step 4. プログラミング言語を学ぶ
    • Step 5. 次のステップへ進むには
    • Step 6. 情報科学を学ぶ
  • プログラミング
    • プログラミング(概要)
    • システムの構成
    • サーバーのセットアップ
      • OSセットアップ
        • ユーザの作成
        • SSHの設定
        • ファイアウォールの設定
      • ドメインの取得と設定(Google domain)
    • Dockerによる仮想化
    • サーバーサイド(golang/gin)
    • リバースプロキシの設置(nginx)
    • SSL/TLS(Let's Encrypt)
    • フロントエンド(Elm)
    • Authentication(Firebase)
    • DB(MySQL)
    • KVS(Redis)
  • 数学勉強ノート
    • 0. 数学勉強ノート(概要)
    • 1. 第1群(基礎1)
      • 1.1. 集合
      • 1.2. ファジィ集合
      • 1.3. 論理学
      • 1.4. 位相
      • 1.5. 代数学
        • 1.5.1. 群
        • 1.5.2. 環
        • 1.5.3. 体
      • 1.6. 多様体論
    • 2. 第2群(基礎2)
      • 2.1. 線形代数学
      • 2.2. 解析学
        • 2.2.1. 常微分
        • 2.2.2. 全微分/偏微分
        • 2.2.3. 反変・共変と不変量
        • 2.2.4. リーマン積分
        • 2.2.5. 外微分
        • 2.2.6. ベクトル解析
      • 2.3. 複素解析学
      • 2.4. 圏論
      • 2.5. 測度論
      • 2.6. 確率論
    • 3. 第3群(情報科学)
      • 3.1. 情報理論
      • 3.2. 通信符号理論
      • 3.3. 暗号理論
      • 3.4. グラフ理論
    • 4. 第4群(基礎3)
      • 4.1. 位相幾何学
      • 4.2. 微分幾何学
      • 4.3. 関数解析学
      • 4.4. 偏微分方程式論
      • 4.5. 数値解析学
        • 4.5.1. 常微分方程式
        • 4.5.2. 偏微分方程式
        • 4.5.3. 固有値問題/特異値分解
      • 4.6. 数理統計学
    • 5. 第5群(物理学)
      • 5.1. 古典力学
      • 5.2. 剛体力学
      • 5.3. 流体力学
      • 5.4. 解析力学
      • 5.5. 電磁気学
      • 5.6. 相対性理論
      • 5.7. 熱力学
      • 5.8. 統計力学
      • 5.9. 量子力学
    • 6. 第6群(電気・機械工学)
      • 6.1. 電気回路学
      • 6.2. システム制御理論
    • 7. 第7群(基礎4)
      • 7.1. 情報幾何学
      • 7.2. 確率解析
  • 書籍メモ
    • 書籍メモ(概要)
  • 論文メモ
    • 論文メモ(概要)
    • 線形代数
      • The Matrix Cookbook
      • ON THE EARLY HISTORY OF THE SINGULAR VALUE DECOMPOSITION
      • The tensor rank decomposition
    • 行列分解
      • Probabilistic Matrix Factorization
      • Probabilistic Matrix Factorization with Non-random Missing Data
      • Graph Regularized Non-negative Matrix Factorization for Data Representation
      • Alternating least squares for personalized ranking
    • 連続最適化
      • Quasi-Newton Methods: Superlinear Convergence Without Line Searches for Self-Concordant Functions
    • 確率・統計
      • Statistical inference of Langevin distributionfor directional data
      • Bayesian Inference over the Stiefel Manifold via the Givens Representation
      • あとで読む(Langevin distribution)
    • Riemannian最適化
      • A Broyden Class of Quasi-Newton Methods for Riemannian Optimization
      • Riemannian stochastic quasi-Newton algorithm with variance reduction and its convergence analysis
      • Primal-Dual Optimization Algorithms over Riemannian Manifolds: an Iteration Complexity Analysis
      • Extending FISTA to Riemannian Optimization for Sparse PCA
    • RCD空間
      • On the geometry of metric measure spaces. I
      • On the geometry of metric measure spaces. II
      • Cohomology and other analyticalaspects of RCD spaces
      • Untitled
  • 雑記帳
    • 球面幾何(楕円幾何)
    • 双曲幾何
  • なんかの部屋
    • データの表現について
Powered by GitBook
On this page
  • 直線探索法
  • 参考文献

Was this helpful?

  1. 連続最適化

4. 直線探索法

2020/1/6公開

直線探索法

直線探索法(line search method)は、現在の点x(k)∈Rnx _ {(k)} \in \mathbb{R} ^ nx(k)​∈Rnからm(k)∈Rnm _ {(k)} \in \mathbb{R} ^ nm(k)​∈Rnにより方向が定められたt∈Rt \in \mathbb{R}t∈Rを媒介変数とする直線x(k)+tm(k)x _ {(k)} + t m _ {(k)}x(k)​+tm(k)​を探索直線として、目的関数f ⁣:Rn→Rf \colon \mathbb{R} ^ n \to \mathbb{R}f:Rn→Rの探索直線上における極値への更新を繰り返すことで最適解へ近づいていく手法の総称である。探索方向m(k)m _ {(k)}m(k)​に目的関数の勾配grad⁡f(x(k))\operatorname{grad} f(x _ {(k)})gradf(x(k)​)の情報を用いる場合は勾配法(gradient descent)ともいう。

すなわち、アルゴリズムとしてはまず現在の点x(k)x _ {(k)}x(k)​から更新方向m(k)m _ {(k)}m(k)​を定め、次に

ddtf(x(k)+tm(k))=0(4.1)\frac{d}{dt} f (x _ {(k)} + t m _ {(k)}) = 0 \tag{4.1}dtd​f(x(k)​+tm(k)​)=0(4.1)

を満たすtttを見つけてt(k)t _ {(k)}t(k)​とおき、

x(k+1)=x(k)+t(k)m(k)x _ {(k+1)} = x _ {(k)} + t _ {(k)} m _ {(k)}x(k+1)​=x(k)​+t(k)​m(k)​

と更新することを繰り返すことになる。このt(k)t _ {(k)}t(k)​を更新のステップサイズ(step size)という。

現実問題として (4.1) 式を満たすtttを見つけること自体もtttに関する最適化問題になっているので、(4.1) 式を満たすtttを見つけるのは難しい場合がある。たとえば目的関数fffを評価するのに時間がかかるような問題では厳密にtttを求めるのは困難である。

ステップサイズtttを厳密に求める直線探索法を厳密直線探索法という。厳密直線探索法には

  • 黄金分割探索法

  • セカント法(割線法)

などがある。

しかし、適当なt≠0t \neq 0 t=0を選んで点x′=x(k)+tm(k)x ^ \prime = x _ {(k)} + t m _ {(k)}x′=x(k)​+tm(k)​を考えたとき、点x′x ^ \primex′における探索方向m′m ^ \primem′は一般にm(k)m _ {(k)}m(k)​とは一致しない。たとえば最急降下法では

m(k)=−grad⁡f(x(k))m _ {(k)} = - \operatorname{grad} f(x _ {(k)})m(k)​=−gradf(x(k)​)

によって探索方向を定義するが、一般に

−grad⁡f(x′)≠−grad⁡f(x(k))- \operatorname{grad} f (x ^ \prime) \neq - \operatorname{grad} f (x _ {(k)})−gradf(x′)=−gradf(x(k)​)

である。つまり (4.1) 式を満たすtttを探す間にも探索すべき方向は刻一刻と変化するから、実用上はステップサイズt(k)t _ {(k)}t(k)​を厳密に求めてもあまり意味がない。

実際、理論的にもステップサイズはあとで解説する Armijo 条件や Wolfe 条件を満たしていれば収束が保証される場合が多い。Armijo 条件を満たす点を探すには

  • バックトラック法(backtracking)

がある。

探索方向に勾配を利用すると、勾配は最適解の付近ではそのノルムが 0 に近づくという性質があるから、実用上は

  • ステップサイズを小さい定数にする

  • ステップサイズを単調減少させる

といった方法でも十分であり、これらの手法もよく用いられる。

しかしその一方で、ステップサイズを

  • 小さすぎる定数にする⇒\Rightarrow⇒収束するまでに非常に時間がかかる

  • 大きすぎる定数にする⇒\Rightarrow⇒アルゴリズムが発散/振動する

  • 単調減少のスケジュールがうまくない⇒\Rightarrow⇒収束までに時間がかかったり発散/振動したりする

といった問題が常につきまとうのも事実であり、ステップサイズのチューニングに時間がかかることもある。

ステップサイズの決定は最適化アルゴリズムの実行時間や収束にクリティカルに影響する重要なパラメータだが、大抵の場合は人間が調整するハイパーパラメータとして扱われ、勘、経験、試行錯誤に頼って決められる。ステップサイズを決定するための方法論は現代においても未熟である。

参考文献

Previous3. 最適化問題Next4.1. 最急降下法

Last updated 5 years ago

Was this helpful?

『数理計画法の基礎』

http://ocw.nagoya-u.jp/files/38/chap_3_book.pdf