【Matlab学习手记】线搜索Goldstein准则

程序猿老王。

于 2018-08-11 10:03:30 发布

阅读量4.8k

点赞数 7

CC 4.0 BY-SA版权

分类专栏：算法系列 Matlab学习手记文章标签： Goldstein准则线搜索

本文链接：https://blog.csdn.net/u012366767/article/details/81563555

Matlab学习手记同时被 2 个专栏收录

80 篇文章

订阅专栏

算法系列

29 篇文章

订阅专栏

本文介绍了一种基于Goldstein准则的步长搜索算法，该算法用于确定优化过程中最佳的步长参数。通过循环迭代逐步缩小搜索范围，直到找到同时满足Goldstein两个准则的步长值。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

代码

function [lam, newxk, fk, newfk] = Opt_Goldstein(xk, dk)
a = 0;       % 搜索下界
b = 10;     % 搜索上界
lam = 1;   % 初始点
c1 = 0.25;
c2 = 0.75;
t = 2;
% 循环求解最佳的步长 
while a < b  
    % 搜索步长满足 Goldstein 第一个准则 
    if fun(xk + lam*dk) <= fun(xk) + c1*lam*gfun(xk)'*dk 
        % 搜索步长满足 Goldstein 第二个准则 
        if fun(xk + lam*dk) >= fun(xk) + c2*lam*gfun(xk)'*dk
            % 输出最佳的步长 
            newxk = xk + lam*dk;
            fk = fun(xk);
            newfk = fun(newxk);
            break; 
        % 搜索步长不满足 Goldstein 准则，继续迭代 
        else 
            a = lam; 
            lam = 0.5*(a + b); 
            if b < inf 
                lam = 0.5*(a + b); 
            else 
                lam = t*lam; 
            end 
        end 
    else 
        b = lam; 
        lam = 0.5*(a + b); 
    end 
    if abs(a - b) < 1e-3
        break;
    end
end

function f=fun(x)
f = 100*(x(1)^2 - x(2))^2 + (x(1)-1)^2;

function gf=gfun(x)
gf = [400*x(1)*(x(1)^2 - x(2)) + 2*(x(1) - 1);
        -200*(x(1)^2 - x(2))];