Skip to content

定理Theorem

精确 ℓ1 罚函数

Exact l1 penalty

以约束违背的绝对值和正部构造有限参数即可精确的罚函数,并说明乘子阈值、非光滑最优性与局部全局边界。

形式陈述 ​

对约束优化问题

minf(x),ci(x)=0,gj(x)≤0,

定义总违约量及罚目标

v(x)=∑i|ci(x)|+∑jmax(gj(x),0),Pρ(x)=f(x)+ρv(x),ρ>0.

“精确”表示某个有限参数之后,罚问题保留原问题的最优解。必须说明是在某个点附近保留局部解,还是两个全局解集相同;这两个结论需要不同条件。

先给一个完整的全局版本。设 f,gj:Rn→R 为连续可微凸函数、ci 仿射,原问题存在最优点 x∗ 与KKT 乘子 (λ∗,ν∗),其中 ν∗≥0,Lagrangian 采用

L(x,λ,ν)=f(x)+λTc(x)+νTg(x).

令 M=max(‖λ∗‖∞,‖ν∗‖∞),空乘子块按零处理。若 ρ>M,则 Pρ 的全局最优解集恰等于原约束问题的全局最优解集。若只取 ρ=M,原最优点仍是罚目标的最小点,但可能出现额外的不可行最小点。

尖角是这一性质的核心。在 ci(x)=0 处,绝对值不能用普通导数处理。上述凸版本可用凸次梯度;非凸约束产生的复合罚项则不能直接套用凸次微分演算,应采用方向导数或另行定义的广义导数。

直觉

二次罚靠离开约束一定距离来产生反向斜率;绝对值罚在约束上就能提供一整段斜率。只要这段斜率足以抵消目标在约束法向上的推动,最优点就能准确停在约束上,不必等参数趋于无穷。

乘子衡量放松约束带来的边际收益。罚率超过这一收益时,微小违约赚到的目标改善不足以付罚金。ℓ1 违约对应乘子的 ℓ∞ 阈值,来自不等式 |λTc|≤‖λ‖∞‖c‖1;更换违约范数,也会更换相应的对偶范数。

有限罚率与精确可行
例子与边界

阈值恰好是二 ​

考虑 f(x)=(x−2)2/2、约束 x=0。原解为零,驻点式 −2+λ=0 给出 λ∗=2。罚目标为

Pρ(x)=12(x−2)2+ρ|x|.

在 x>0 分支,导数为 x−2+ρ,所以候选 x=2−ρ 只有在 ρ<2 时适用。在 x<0 分支,导数为 x−2−ρ,其零点不在负半轴。最后检查尖点:

0∈−2+ρ[−1,1]⟺ρ≥2.

因此唯一最小点是 xρ=max(2−ρ,0)。例如 ρ=1 时为一,ρ=2 和三时均为零。严格凸的平方项使本例在临界参数处仍唯一。

对比二次罚:f(x)+μx2/2 的最小点为 2/(1+μ),对所有有限 μ 都大于零。这是不同的法向斜率机制,不是求解器容差造成的区别。

为什么一般定理要求严格超过阈值 ​

取 f(x)=−x、约束 x=0。最优乘子为一。当 ρ=1 时,P1(x)=−x+|x| 在整个 [0,∞) 上都等于零,所以包含无数不可行最小点;当 ρ>1 时,两侧斜率都指向零,才只留下原解。

精确罚也不是“所有驻点都可行”。对非凸约束,罚目标仍可能有不可行驻点;对 f(x)=−x4,c(x)=x,任何有限 ρ 的罚目标在远处都趋于负无穷,虽然零附近存在精确的局部最低点。局部精确性不能替代全局可解性。

推论与应用

全局解集相同的短证明 ​

凸性和 KKT 保证 L(x,λ∗,ν∗)≥f(x∗) 对所有 x 成立。逐项比较有

ρ|ci|−λ∗,ici≥(ρ−M)|ci|,ρmax(gj,0)−ν∗,jgj≥(ρ−M)max(gj,0).

第二式在 gj≤0 时左侧为 −ν∗,jgj≥0,在 gj>0 时直接成立。相加得到

Pρ(x)≥f(x∗)+(ρ−M)v(x).

原问题的任意最优点可行,罚目标等于 f(x∗);而 ρ>M 时,任何不可行点都有严格更大的罚目标。于是两个最优解集一致。证明同时给出违约证书:若 Pρ(x)−f(x∗)≤ε,则 v(x)≤ε/(ρ−M)。数值近似最小化仍只得到近似可行,精确性定理没有消除求解误差。

非凸等式问题中的局部版本 ​

设仅有光滑等式 c(x)=0,Jc(x∗) 满行秩,(x∗,λ∗) 满足驻点式,而且

dT∇xx2L(x∗,λ∗)d>0对所有 d≠0,Jc(x∗)d=0.

这是约束切空间上的二阶充分条件。可选一个足够大的有限 a>0,使 ∇xx2L+aJcTJc 正定:罚项控制法向,原来的二阶条件控制切向。因此 x∗ 是 L(x,λ∗)+a‖c(x)‖2/2 的严格局部最小点。

当 ρ>‖λ∗‖∞ 时,在充分小的邻域内,

ρ‖c(x)‖1−λ∗Tc(x)−a2‖c(x)‖2≥0,

因为线性违约项压过二次小量。把这项加回前述严格局部最小关系,便证明 x∗ 也是 Pρ 的严格局部最小点。这里的局部结论由明确的正则性和二阶条件支撑,没有使用非凸问题的全局对偶下界。

怎样把精确罚用于算法 ​

罚率往往未知,可以用当前 QP 或 KKT 乘子估计作线索,再检查约束是否改善。固定点处的绝对值尖角不宜直接交给假设处处二阶光滑的 Newton 程序。对仿射约束,可以引入 ui≥±ci(x) 与 wj≥gj(x),wj≥0,最小化 f(x)+ρ∑iui+ρ∑jwj;这样用额外变量显式表示非光滑项。

序列二次规划常把 Pρ 用作步长验收函数,同时计入目标与违约。它在此是一种求解工具。Lasso中的 ℓ1 则是希望保留的建模惩罚;虽然某些一维子问题都出现软阈值,两者对参数和最终目标的含义不同。

参考资料
关系图谱10 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系