不精确线性搜索算法 matlab程序,DFP算法+wolfe性非线性搜索解决无约束问题的matlab程序...

以下是两个 Matlab 程序，一个是不精确线性搜索算法，另一个是使用 DFP 算法和 Wolfe 性非线性搜索解决无约束问题的程序：

不精确线性搜索算法：

function [x, fx, n] = inexact_line_search(f, x0, d, alpha, rho, c)
% f: 目标函数
% x0: 初始点
% d: 搜索方向
% alpha: 初始步长
% rho: 步长缩小比例
% c: Armijo准则常数
% x: 最优解
% fx: 目标函数在最优解处的值
% n: 迭代次数

if nargin < 4
    alpha = 1;
end
if nargin < 5
    rho = 0.5;
end
if nargin < 6
    c = 1e-4;
end

n = 0;
while true
    x = x0 + alpha * d;
    fx0 = feval(f, x0);
    fx = feval(f, x);
    if fx <= fx0 + c * alpha * (d' * grad(f, x0))
        break;
    end
    alpha = rho * alpha;
    n = n + 1;
end

DFP 算法和 Wolfe 性非线性搜索解决无约束问题：

function [x, fx, n, alpha] = dfp_wolfe(f, x0, eps, max_iter)
% f: 目标函数
% x0: 初始点
% eps: 精度要求
% max_iter: 最大迭代次数
% x: 最优解
% fx: 目标函数在最优解处的值
% n: 迭代次数
% alpha: 步长

if nargin < 3
    eps = 1e-6;
end
if nargin < 4
    max_iter = 1000;
end

n = 0;
x = x0;
H = eye(length(x));
while true
    fx = feval(f, x);
    g = grad(f, x);
    if norm(g) < eps
        break;
    end
    d = -H * g;
    [alpha, nls_n] = wolfe(f, x, d);
    n = n + nls_n;
    s = alpha * d;
    x1 = x + s;
    y = grad(f, x1) - g;
    H = H + (s * s') / (s' * y) - (H * y * y' * H) / (y' * H * y);
    x = x1;
    n = n + 1;
    if n > max_iter
        break;
    end
end

其中，grad 函数是计算梯度的函数，wolfe 函数是 Wolfe 性非线性搜索的函数。你需要自己实现这两个函数。