机器学习笔记-牛顿法搜索方向的相关证明

关于牛顿法搜索方向的相关证明

  关于之前介绍牛顿法中提出的:牛顿法的搜索方向不一定就是下降方向
  既然说牛顿法的搜索方向可以不是下降方向,那我们就必须知道在满足什么条件下牛顿法的搜索方向才是下降方向呢?下面先看一个定理:
f f f R n R^{n} Rn-> R 1 R^{1} R1在点 x ∈ R n x\in{R^n} xRn可微。若存在 p ∈ R n p\in R^{n} pRn,使:
∇ f ( x ) T p < 0 \nabla f(x) ^{T} p<0 f(x)Tp<0
则向量 p p p f f f在点 x x x处的下降方向。
  在牛顿法中我们知道搜索方向为 − ∇ 2 f ( x ) − 1 ∇ f ( x ) -\nabla^{2}f(x)^{-1}\nabla f(x) 2f(x)1f(x),所以就可以认为
p = − ∇ 2 f ( x ) − 1 ∇ f ( x ) p=-\nabla^{2}f(x)^{-1}\nabla f(x) p=2f(x)1f(x),只需要证明下面的式子成立即可:
− ∇ f ( x ) T ∇ 2 f ( x ) − 1 ∇ f ( x ) < 0 -\nabla f(x) ^{T}\nabla^{2}f(x)^{-1}\nabla f(x)<0 f(x)T2f(x)1f(x)<0
即等价于证明:
∇ f ( x ) T ∇ 2 f ( x ) − 1 ∇ f ( x ) > 0 \nabla f(x) ^{T}\nabla^{2}f(x)^{-1}\nabla f(x)>0 f(x)T2f(x)1f(x)>0
而上述证明又可以化简为:
∇ 2 f ( x ) − 1 > 0 \nabla^{2}f(x)^{-1}>0 2f(x)1>0
  这其实就是黑塞矩阵正定的条件。也就是需要满足当前点的黑塞矩阵是正定的才能保证牛顿法的迭代方向一定是沿着函数值下降的方向。
  而牛顿法在迭代过程中是不一定就能保证每次迭代所求得的黑塞矩阵都是正定的,所以说牛顿法的迭代方向不一定就是函数值的下降方向。

  • 1
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是可以直接在 Matlab 中运行的代码,使用牛顿求解无约束优化问题: ```matlab % 定义目标函数 f(x) f = @(x) 100*(x(1)^2 - x(2))^2 + (x(1) - 1)^2; % 定义目标函数的一阶和二阶导数 df = @(x) [400*x(1)*(x(1)^2 - x(2)) + 2*(x(1) - 1); -200*(x(1)^2 - x(2))]; d2f = @(x) [1200*x(1)^2 - 400*x(2) + 2, -400*x(1); -400*x(1), 200]; % 设置初始点和容差 x0 = [0; 0]; tol = 1e-6; % 牛顿主循环 while true % 计算当前点的一阶和二阶导数 dfx = df(x0); d2fx = d2f(x0); % 计算搜索方向 d d = -d2fx \ dfx; % 计算步长 alpha alpha = 1; % 更新 x0 的值 x = x0 + alpha*d; % 判断是否满足终止条件 if norm(x - x0) < tol break; end % 更新 x0 的值 x0 = x; end % 输出最优解和最优值 fprintf('最优解为:x1 = %.4f, x2 = %.4f\n', x(1), x(2)); fprintf('最优值为:%.4f\n', f(x)); ``` 在上面的代码中,我们首先定义了目标函数 $f(x)$,以及它的一阶和二阶导数。然后,我们设置了初始点和容差,并进入了牛顿的主循环。在循环中,我们先计算当前点的一阶和二阶导数,然后计算搜索方向和步长,更新 $x_0$ 的值,并判断是否满足终止条件。最后,我们输出最优解和最优值。 请注意,牛顿在计算二阶导数的时候需要保证目标函数的二阶导数存在,并且是正定的。如果二阶导数不存在或不是正定的,牛顿可能会出现收敛问题。在这种情况下,可以考虑使用其他的优化算,比如拟牛顿

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值