Matlab学习手记——线搜索Wolfe准则

最新推荐文章于 2024-05-15 11:43:35 发布

程序猿老王。

最新推荐文章于 2024-05-15 11:43:35 发布

阅读量1w

点赞数 9

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

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

80 篇文章 74 订阅

订阅专栏

算法系列

29 篇文章 9 订阅

订阅专栏

测试：

>> [alpha, newxk, fk, newfk] = Opt_Wolfe([0 0]', [1 0]')
alpha =
    0.1250
newxk =
    0.1250
         0
fk =
     1
newfk =
    0.7900

代码如下：

function [alpha, newxk, fk, newfk] = Opt_Wolfe(xk, dk)
rho = 0.25; 
sigma = 0.75;
alpha = 1; 
a = 0; 
b = 1; 
while a < b
    if fun(xk + alpha*dk) <= fun(xk) + rho*alpha*gfun(xk)'*dk
        if gfun(xk + alpha*dk)'*dk >= sigma*gfun(xk)'*dk
            break;
        else
            a = alpha;
            alpha = 0.5*(a + b);
        end
    else
        b = alpha;
        alpha = (alpha+a)/2;
    end
end
newxk = xk + alpha*dk;
fk = fun(xk);
newfk = fun(newxk);

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))];