dlinmin.m

来自「无线传感器节点定位算法」· M 代码 · 共 22 行

%|
%| REFERENCE:  Numerical recipes in C
%|

function [p, xi, fret] = dlinmin(p, xi, func, dfunc)

TOL    = 2.0e-4;  % Tolerance passed to brent.

global pcom xicom nrfunc nrdfun;
nrfunc = func;
nrdfun = dfunc;
pcom   = p;
xicom  = xi;

ax     = 0.0;  % Initial guess for brackets.
xx     = 2.0;
[ax, xx, bx] = minBracket(ax, xx, 'f1dim');
[fret, xmin] = dbrent(ax,xx,bx,'f1dim', 'df1dim',TOL);

xi     = xi.*xmin;
p      = p + xi;