Broyden法による求根

Broyden法(ブライデン法)は複数の方程式から得られる多次元の解を求める数値計算手法(球根アルゴリズム)で、セカント法を一般化した手法.セカント法が傾きで微分を近似したように、計算が複雑なヤコビ行列を一つ前のヤコビ行列の更新で実現しているのが特徴.

セカント法についてはリンク参照ください.

Broyden法の流れ 1. 初期ヤコビアン行列\(\boldsymbol{J}_0\)と解に近い二つの多次元ベクトル\(\boldsymbol{x}_0, \boldsymbol{x}_1\)、関数\(\boldsymbol{f}(x)\)を用意 2. 下記の一つ目のヤコビアン更新式に基づいてヤコビアン\(\boldsymbol{J}_n\)を算出. 3. 下記の二つ目の更新式に基づいて\(\boldsymbol{x}_{n+1}\) を算出. 4. 2.と3.を繰り返して求められる精度になるまで繰り返す. ヤコビアン更新式 $$\boldsymbol{J_n} = \boldsymbol{J_{n-1}} + \frac{ \Delta\boldsymbol{f}_n - \boldsymbol{J}_{n-1} \Delta{x_{n } }}{||\Delta x_n||^2 } \Delta{x_n}^T$$ 更新式 $$\boldsymbol{x}_{n+1} = \boldsymbol{x}_n - \boldsymbol{J}_n^{-1}\boldsymbol{f}(\boldsymbol{x }_n)$$ 別の作戦 上記ではヤコビアン\(J_n\)を更新したが、更新式で欲しいのは\(J_n^{-1}\)であるので、それを直接更新する方法も提案されている.