自在学

我们与你共同进步

  • 分类课程
  • 文章
  • 工作台
  • 订阅

  • 关于我们
  • 隐私政策
  • 使用条款

探索

  • 分类课程
  • 文章
  • 工作台
  • 订阅

网站信息

  • 关于我们
  • 隐私政策
  • 使用条款

加入社区

自在学学习社区微信二维码

微信扫码,交流学习

株洲市自在学教育科技有限公司© 2025 - 2026 版权所有

© 2025 - 2026 株洲市自在学教育科技有限公司 版权所有

湘公网安备43020302000292号|湘ICP备2025148919号-1
分类课程工作台文章订阅
分类课程工作台文章价格

优化与运筹基础:从模型到最优决策

  1. 01从问题到优化模型
  2. 02凸性:为什么局部信息能指向全局
  3. 03线性规划的几何与单纯形思想
  4. 04对偶与资源价值
  5. 05整数优化与离散决策
  6. 06网络优化:路径、流与割
  7. 07无约束数值优化:梯度与 Newton
  8. 08KKT 与约束最优性
  9. 09梯度算法的工程化使用
  10. 10综合决策与模型审计
正在加载课程章节内容
课程数学优化与运筹基础:从模型到最优决策梯度算法的工程化使用

梯度算法的工程化使用

同一个下降方向,步子不同,可能得到快速收敛、来回震荡,也可能直接发散。实际使用梯度算法时,算法名只占很小一部分;步长规则、变量尺度、停止条件和约束处理,才决定结果是否可信。

固定步长、精确线搜索与回溯线搜索的轨迹对照

三种步长策略

固定步长使用预先选定的 α\alphaα,实现简单但需要知道曲率尺度。精确线搜索沿给定方向求

αk=arg⁡min⁡α≥0f(xk+αdk).\alpha_k=\arg\min_{\alpha\ge0}f(x_k+\alpha d_k).αk​=argα≥0min​f(xk​+αdk​).

对二次函数和方向 dk=−gkd_k=-g_kdk​=−gk​,令一维函数导数为零可得

αk=gkTgkgkTQgk.\alpha_k=\frac{g_k^Tg_k}{g_k^TQg_k}.αk​=gkT​Qgk​gkT​gk​​.

回溯线搜索从较大的候选步长开始,不断乘以 β∈(0,1)\beta\in(0,1)β∈(0,1),直到满足 Armijo 充分下降:

f(xk+αdk)≤f(xk)+σαgkTdk,f(x_k+\alpha d_k)\le f(x_k)+\sigma\alpha g_k^Td_k,f(xk​+αdk​)≤f(xk​)+σαgkT​dk​,

其中 dkd_kdk​ 是下降方向、σ∈(0,1)\sigma\in(0,1)σ∈(0,1)。右侧允许的下降量与一阶预测成比例,过大的步长会被缩小。

Armijo 条件把实际函数值与一阶预测比较

条件数与锯齿

正定 Hessian 的条件数 κ=λmax⁡/λmin⁡\kappa=\lambda_{\max}/\lambda_{\min}κ=λmax​/λmin​ 衡量等高线的狭长程度。κ\kappaκ 大时,负梯度几乎横穿狭长谷底,下一步又从另一侧折回,导致锯齿和慢收敛。合适的线搜索能避免发散,却不能消除几何上的病态。

特征值比例改变时等高线从圆变成狭长椭圆

预条件化用变量变换或矩阵 PPP 改善尺度,使新坐标中的曲率更均衡。它没有改变原问题的含义,只改变算法看到的几何。若 PPP 选得不好,计算代价或数值误差也可能上升。

预条件化前后同一问题的等高线与轨迹

有约束时投影回可行域

若变量必须满足简单集合 CCC,投影梯度更新为

xk+1=ΠC(xk−αk∇f(xk)),x_{k+1}=\Pi_C(x_k-\alpha_k\nabla f(x_k)),xk+1​=ΠC​(xk​−αk​∇f(xk​)),

其中 ΠC(z)\Pi_C(z)ΠC​(z) 是离 zzz 最近的可行点。对盒约束 li≤xi≤uil_i\le x_i\le u_ili​≤xi​≤ui​,投影就是逐坐标截断;对概率向量单纯形,投影需要保持非负且总和为 1。投影不是把不可行解“修饰一下”,它改变了实际搜索方向,停止时应结合 KKT 残差判断。

投影梯度把越过边界的点拉回可行集合

诊断一份迭代日志

至少同时看目标值、梯度范数、步长和约束违反量。目标下降但约束违反越来越大,说明实现把投影/可行性遗漏了;梯度范数不降且步长不断缩小,可能是尺度、非凸性或错误梯度;目标值上下交替增长,优先检查步长和方向符号。

迭代日志中的目标值、梯度范数与约束违反量

练习

  1. 为什么 Armijo 条件右侧要含 gkTdkg_k^Td_kgkT​dk​?

gkTdk<0g_k^Td_k<0gkT​dk​<0 表示方向是一阶下降方向。右侧用它给出随步长缩放的一阶下降基准,实际下降不够时就拒绝该步。

  1. 条件数很大时,精确线搜索是否一定让梯度法一步到解?

不一定。精确线搜索只在当前方向上选最优步长,狭长谷底的梯度方向仍可能横向摆动;条件数大通常意味着整体收敛变慢。

  1. 盒约束投影如何处理候选值 zi<liz_i<l_izi​<li​?

取 ΠC(z)i=li\Pi_C(z)_i=l_iΠC​(z)i​=li​。它把该坐标截回可行区间;若 zi>uiz_i>u_izi​>ui​ 则截到 uiu_iui​,区间内则保持原值。

1
目标函数下降就足以说明每个迭代点都满足约束。
上一章KKT 与约束最优性下一章综合决策与模型审计