Mnfb Mnfb_r Mngb Mngb_r Mnhb Mnhb_r

Mnfb() Sub Mnfb ( N As  Long, X() As  Double, B() As  Double, F As  LongPtr, Info As  Long, Optional Itsum As  LongPtr = NullPtr, Optional Info2 As  Long, Optional NFcall As  Long, Optional NGcall As  Long, Optional Niter As  Long, Optional Fval As  Double, Optional Rtol As  Double = -1, Optional Atol As  Double = -1, Optional MaxFcall As  Long = 0, Optional MaxIter As  Long = 0, Optional Tuner1 As  Double = -1, Optional Xctol As  Double = -1, Optional Xftol As  Double = -1, Optional Lmax0 As  Double = -1, Optional Lmaxs As  Double = -1, Optional Sctol As  Double = -1, Optional Bias As  Double = -1, Optional Eta0 As  Double = -1  ) 多変数非線形関数の最小点 (信頼領域法) (単純制約付き) (微分係数不要) 目的 本ルーチンは非線形関数 f(x1, x2, ..., xn) (2階連続的微分可能な実数関数) の局所的最小点(xs1, xs2, ..., xsn) を次の制約条件のもとで求める(b0iとb1iはそれぞれ上限値と下限値). b0i <= xsi <= b1i (i = 1 〜 n) 目的関数の微分係数は, これを解析的に求めるユーザー定義サブルーチンを使わずに, 有限差分近似により求められる. ヘッセ行列はセカント法(BFGS法)により求められる. 探索アルゴリズムはダブルドッグレッグ信頼領域法を使用する. 引数 [in] N 問題の次数(パラメータ数). (N > 0) [in,out] X() 配列 X(LX - 1) (LX >= N) [in] 解ベクトルの初期推定値. [out] 求められた解ベクトル. [in] B() 配列 B(LB1 - 1, LB2 - 1) (LB1 = 2, LB2 >= N) 解ベクトル X の境界制約条件.   B(0, I) <= X(I) <= B(1, I) (I = 0 〜 N-1) [in] F 目的関数 f(x1, x2, ..., xn) を求めるユーザー定義サブルーチンで, 次のように定義すること. Sub F(N As Long, X() As Double, Nf As Long, Fval As Double) 与えられたNおよびX()から関数値を求めFvalに設定する. Nfは呼び出しカウンタである. 与えられたX()が計算可能範囲外の場合, Nf = 0 に設定して戻ること. それ以外の変数を変更してはならない. End Sub [out] Info = 0: 正常終了. (サブコードをInfo2に返す) = -1: パラメータ N の誤り. (N < 1) = -2: パラメータ X() の誤り. = -3: パラメータ B() の誤り. = 7: 特異収束. (近傍のヘッセ行列が特異になった) = 8: 誤収束. (誤った点での収束と思われる. 目標精度が小さすぎる可能性がある) = 9: 関数評価回数の最大値を超えた. = 10: 反復回数の最大値を超えた. = 63: X()の初期点においてF(X)を求めることができない. = 65: X()において微分係数を求めることができない. [in] Itsum (省略可) 計算の途中経過を出力するためのユーザー定義サブルーチンで, 次のように定義すること. (省略時 = NullPtr) アドレスを与えると(NullPtr以外であれば)反復ごとにサブルーチンが呼び出される. Sub Itsum(N As Long, X() As Double, NIter As Long, Nf As Long, Ng As Long, Fval As Double) 以下の情報が渡されるのでこれを任意の形式で出力する. N: 変数の数. X(): 解ベクトルの現在の推定値. NIter: 何回目の反復かを示す. Nf: Fの呼び出し回数 (微分係数計算のための呼び出しを除く). Ng: Fの呼び出し回数 (微分係数計算のための呼び出し回数). Fval: X()におけるFの値. End Sub 引数の値を変更してはならない. [out] Info2 (省略可) Info = 0 のときのサブコード. = 1: X値の収束条件を満たした. = 2: 関数値相対収束条件を満たした. = 3: X値および関数値相対収束条件の両方を満たした. = 4: 関数値絶対収束条件を満たした. [out] NFcall (省略可) 関数Fの呼び出し回数(微分係数計算のための呼び出しを除く). [out] NGcall (省略可) 関数Fの呼び出し回数(微分係数計算のための呼び出し回数). [out] Niter (省略可) 反復回数. [out] Fval (省略可) 求められた解ベクトルX()における関数値. [in] Rtol (省略可) 関数値の目標相対精度. (Eps <= Rtol <= 0.1) (省略時 = 1e-10) (Epsはマシンイプシロン) (Rtol < Eps または Rtol > 0.1 であれば省略時の既定値とみなす) [in] Atol (省略可) 関数値の目標絶対精度. (省略時 = 1e-20) (Atol < 0 であれば省略時の既定値とみなす) [in] MaxFcall (省略可) 関数Fの呼び出し回数の最大値. (省略時 = 200) (MaxFcall <= 0 であれば省略時の既定値とみなす) [in] MaxIter (省略可) 反復回数の最大値. (省略時 = 150) (MaxIter <= 0 であれば省略時の既定値とみなす) [in] Tuner1 (省略可) 誤収束の判定パラメータ. (0 <= Tuner1 <= 0.5) (省略時 = 0.1) (Tuner1 < 0 または Tuner1 > 0.5 であれば省略時の既定値とみなす) [in] Xctol (省略可) X値の収束判定しきい値. (0 <= Xctol <= 1) (省略時 = Eps^(1/2)) (Xctol < 0 または Xctol > 1 であれば省略時の既定値とみなす) [in] Xftol (省略可) 誤収束の判定しきい値. (0 <= Xftol <= 1) (省略時 = 100*Eps) (Xftol < 0 または Xftol > 1 であれば省略時の既定値とみなす) [in] Lmax0 (省略可) スケーリングされた一番始めのステップ長の最大幅. (Lmax0 > 0) (省略時 = 1) (Lmax0 <= 0 であれば省略時の既定値とみなす) [in] Lmaxs (省略可) [in] Sctol (省略可) LmaxsおよびSctolは特異収束の判定パラメータである. (Lmaxs > 0) (0 <= Sctol <= 1) (省略時: Lmaxs = 1, Sctol = 1e-10) ステップ長 Lmaxs において Sctol*abs(f)よりも関数値が減少しなくなった (fは現在の反復の出発点における関数値). (Lmaxs <= 0 であれば省略時の既定値とみなす) (Sctol < 0 であれば省略時の既定値とみなす) [in] Bias (省略可) ドッグレッグ信頼領域法で使用するパラメータ. (0 <= Bias <= 1) (省略時 = 0.8) (Bias < 0 または Bias > 1 であれば省略時の既定値とみなす) [in] Eta0 (省略可) Fの推定計算精度. (真値 = 計算値*(1 + E)とすると Abs(E) <= Eta0) (Eps <= Eta0 <= 1) (省略時 = 1000*Eps) (Epsはマシンイプシロン) (Eta0 < Eps または Eta0 > 1 であれば省略時の既定値とみなす) 出典 netlib/port 使用例 次の関数の最小点を求める(ローゼンブロック関数). f(x1, x2) = 100(x2 - x1^2)^2 + (1 - x1)^2 初期値は, (x1, x2) = (-1.2, 1) とする. また, 解ベクトル X の境界制約条件を -2 <= x1 <= 2, 0 <= x2 <= 2 とする. Sub FMnfb(N As Long, X() As Double, Nf As Long, F As Double) F = 100 * (X(1) - X(0) ^ 2) ^ 2 + (1 - X(0)) ^ 2 End Sub Sub Ex_Mnfb() Const N As Long = 2 Dim X(N - 1) As Double, B(1, N - 1) As Double, Info As Long X(0) = -1.2: X(1) = 1 B(0, 0) = -2: B(1, 0) = 2 B(0, 1) = 0: B(1, 1) = 2 Call Mnfb(N, X(), B(), AddressOf FMnfb, Info) Debug.Print "X1, X2 =", X(0), X(1) Debug.Print "Info =", Info End Sub 実行結果 X1, X2 = 1.00000000015727 1.00000000031223 Info = 0