matlab程序阻尼牛顿法,利用matlab实现基于阻尼牛顿法函数的极小值确定

% ******forwardBackward.m******%

function [minx,maxx] = forwardBackward (f,x0,h,eps)

%目标函数：f；

%初始点：x0；

% 初始步长：h；

% 精度：eps；

% 目标函数读取包含极值区间左端点：minx；

% 目标函数读取包含极值区间右端点：maxx；

% 处理输入三个参数时情况

if nargin == 3

eps = 1.0e-6;

end

x1 = x0;

f1 = subs(f,findsym(f),x1);

x2 = x0 + h;

f2 = subs(f,findsym(f),x2);

if f1 > f2

x1 = x2;

x3 = x2;

x2 = x1 + h;

f1 = subs(f,findsym(f),x1);

f2 = subs(f,findsym(f),x2);

while f1 > f2

h = 2 * h;

x1 = x3;

x2 = x1 + h;

f1 = subs(f,findsym(f),x1);

f2 = subs(f,findsym(f),x2);

end

minx = x1;

maxx