yyoshiki41/Multipoint_quasi-Newton

準ニュートン法とその改良

Multipoint quasi-Newton method

更新式

x_p(k + 1) = x_p(k) + v_p(k)
v_p(k) = λ * v_p(k - 1) + c1 * r1_p(k) * {-H_p(k) * ∇f(x_p(k))^T} + c2 * r2_p(k) * {x_g-best(k) - x_p(k)}

ベンチマーク関数

• 2n-minima

// 2変数
f(x, y) = (x^4 - 16*x^2 + 5*x) + (y^4 - 16*y^2 + 5*y)
// 多変数
f(x) = Σ (x^4 - 16*x^2 + 5*x)

• rosenbrock

// 2変数
f(x, y) = (1-x)^2 + 100*(y - x^2)^2
// 多変数
f(x_n) = Σ {(1-x_n)^2 + 100*(x_n+1 - x_n^2)^2}

直線探索

黄金分割法

ディレクトリ

.
├── 2n-minima             2n-minima 関数
│   ├── normal            普通の準ニュートン法(ステップ幅 c=0.01)
│   ├── gold              黄金分割法
│   ├── n-dimensions      n次元への拡張
│   ├── run               試行回数 m回での収束回数
│   └── function_call   function call数と目的関数値をカウント
|
└── rosenbrock            rosenbrock 関数
    ├── normal            普通の準ニュートン法(ステップ幅 c=0.01)
    ├── gold              黄金分割法
    ├── n-dimensions      n次元への拡張
    ├── run               試行回数 m回での収束回数
    └── function_call   function call数と目的関数値をカウント