1. 引言
考虑如下无约束优化问题:
(1)
其中
是二次可微函数,
是一个可分离变量的非凸函数。
MCP函数是一种重要的可分离变量的函数,它有两个标量
和
,并以原点为唯一的零点:
(2)
其中
(3)
在机器学习和统计学中,为了防止模型过拟合和提高模型的泛化能力,需要通过引入凹惩罚函数来对模型进行正则化。变量选择在高维数据分析中起着非常重要的作用,在对模型正则化过程中可以通过对参数的惩罚来实现变量选择的效果,从而在一定程度上实现变量选择的功能,然而许多方法容易忽略在变量选择过程中的随机误差。2001年,Fan和Li [1] 为了解决这个问题,提出了一种新的非凸惩罚函数,即光滑裁剪绝对偏差(smoothly clipped absolute deviation, SCAD)罚函数。伴随着SCAD罚函数的出现,oracle性质成为了评价变量选择方法的标准。2010年,Zhang [2] 提出了新的非凸MCP函数来进行变量选择,并证明了MCP函数具有oracle性质。
在正则化模型中,通过引入MCP函数可以对模型的复杂度进行控制,防止模型过拟合,提高模型的泛化能力。在高维数据的情况下,MCP函数可以帮助减少模型的自由度,提高模型的稳定性和可解释性。MCP函数的参数
和a可以根据具体问题进行调节,以满足不同的需求,这使得MCP函数可以适用于各种不同的优化问题。
目前为止,由于MCP函数良好的oracle性质和它可以同时实现变量选择与参数估计,带MCP函数正则化问题被广泛应用于机器学习、统计学习、信号处理、图像处理等领域,详见参考文献 [3] [4] [5] [6] [7] 。为了计算带MCP函数正则化问题,Zhang [2] 采用了惩罚线性无偏选择(PLUS)算法,在此之后Shi等 [8] 开发了一种有效的交替方向乘子法(ADMM)求解带MCP函数的正则化问题。
为了研究模型(1)相关算法的收敛速度,我们需要研究该模型的稳定性,如强度量次正和二阶增长条件,他们已经被用来建立优化算法的收敛速度,详见文献 [9] [10] [11] 。为了刻画强度量次正则 [12] [13] [14] 和二阶增长条件 [13] [15] ,许多学者研究了图像导数,它是变分分析中的一类重要工具。Chieu等人 [13] 利用次微分的图形导数研究了有限维中二阶增长条件和次微分的强度量次正则之间的关系,Dontchev等人 [12] 给出了集值映射强度量次正则的图像导数准则,Bello-Cruz等人 [15] 利用图像导数刻画了
-正则化优化问题的二阶增长条件。到目前为止,大多数文献都是对带MCP函数正则化问题的应用与相关算法的研究,还没有文献对模型(1)次微分的强度量次正则进行研究,而强度量次正则在模型相关算法的收敛性分析中起着非常重要的作用。
受文献 [12] [13] [15] 的启发,本文研究了模型(1)次微分的强度量次正则。本文的结构如下:第二节,我们介绍了一些基本的定义和引理,并利用次微分的定义得到了MCP函数次微分的表达式。第三节,利用次微分的运算法则和集值映射的图像导数工具得到了MCP函数次微分的图像导数表达式,最后利用该表达式分别建立了模型(1)次微分强度量次正则的充分条件和充要条件。
2. 预备知识
本文中,
表示实数集,
是n维欧氏空间,
和
分别表示
中的范数和内积。
表示符号函数。函数
的定义域记作dom
。
是一个二次连续可微函数,
表示函数f的Hessian矩阵,
表示Hessian矩阵的第i行第j列元素,
表
示Hessian矩阵的第i行所有元素。对于矩阵
,定义
。对集值映
,定义F的图为
;定义集值映射F的逆映射
为:
。
2.1. 基础知识
定义1 [16] 设
。对于函数
和向量
,如果存在
和x的一个邻域V,使得
则称v是函数h在x处的邻近次梯度。所有这样的v构成的集合称为h在x处的邻近次微分,记作
。
定义2 [16] 函数
在
处的极限次微分定义如下:
其中
意味着
,
。显然
。
在本文中我们使用的都是极限次微分,简称为次微分。
定义3 [12] 设
是
中的非空子集,
。集合
在
处的切锥定义为:
定义4 [12] 对于集值映射
和点
。H在
处的图像导数是集值映射
,其图像是
在
处的切锥
:
或者说,
当且仅当存在序列
,
和
使得对于所有k,
定义5 [12] 映射
。对于点
,如果存在
和
的一个邻域V,使得
则称映射F在
处对于
是度量次正则的。如果
是
的孤立点,则称F在
处对于
是强度量次正则的。
引理1 [12] 集值映射
在
是局部闭的,H在
是强度量次正则的当且仅当
(4)
引理2 [12] 对于函数
在x处可微,集值映射
和任意的
,有
2.2. MCP函数的次微分表达式
命题1 MCP函数的次微分表达式:
(5)
证明 因为函数
是连续的可分离变量的函数,因此对于任意的
,有
因为单变量函数
除在原点外是连续可微的,因此,当
时,
当
时,我们需要先求得
时的邻近次微分。设
,则
。设
,则
如果
,我们有
。因此
,即
和
如果
,那么对于任意的
都成立。
如果
,我们有
。因此
,即
和
综合上述三种情况,我们可以得到当
时,
。再根据极限次微分的定义,我们可以得到
。因此,MCP函数的次微分表达式为:
3. 带MCP函数正则化问题次微分的强度量次正则
在本节中,我们研究问题(1)次微分的强度量次正则性。为了使用引理1来刻画
的强度量次正则性,我们需要先计算
的图像导数。
3.1.
的图像导数
为了研究
的图像导数,设
,
,定义:
命题2
是非空的当且仅当
。此外,对任意的
有
(6)
证明 取
,根据图像导数的定义,存在序列
和
使得
。下面考虑j的六种情形。
情形1.1
,即
且
。
若
且
,则
,
。此时存在足够大的常数L和
,使得当
时,
,因此这种情况是不满足的。若
且
,则
,
。当
,
,
。
若
且
,则
,
。当
时,
,因此这种情况是不满足。若
且
,则
,
。当
,
,
。
情形1.2
即
且
。
如果有一个序列
,使得
,我们可以得到
,当
,有
。对于
,有
。因为
,当
,有
即
。
若
,有
,则
,此时,
,这意味着
。若
,有
,则
,此时,
,这意味着
。
因此,当
时,有
,
和
。
对于
,我们可以分为三类情况讨论:
1)
,此时
。
2)
,此时
。
3)
,此时
。那么有
,当
,有
。
因此,当
时,有
和
。
情形1.3
,即
。
此时
。对于足够大的k,
,由(5)我们可以得到
当
时,有
;当
时,有
。这意味着
。
因此对任意的
,
。
情形1.4
,即
。
对于足够大的k,
,由(5)知此时
因此
,这意味着
。
情形1.5
,即
。
由(5)知此时
。若
,此时
,由 可得
,因此
,这意味着
;若
,此时
,由(5)可得
,因此
,这意味着
。
情形1.6
,即
。
由(5)式知此时
。若
,此时
,由(5)可得
,因此
,这意味着
;若
,此时
,由(5)可得
,因此
,这意味着
。
综合以上六个情形,我们得出
,同时也验证了(6)中的包含“
”。为了证明逆向包“
”,取
和任意的v属于(6)的右半部分,对于任意的
。我们证明了
,从而验证了
。对于任意的
,我们定义集值映射:
与证明“
”一样,我们考虑j的六种情形。
情形2.1
,即
且
。
因为
,所以有
,当k足够大时
。因此
情形2.2
,即
且
。
若
,我们有
;并且对于足够大的k,因为
,所以
。若
,我们有
。当k足够大时,因为
,所以有
综合上述两种情况,我们有
。
情形2.3
,即
。
如果
,由(5)知
。对于足够大的k,我们有
,因为
,所以
如果
,由(5)知
。对于足够大的k,我们有
,因为
,所以
情形2.4
,即
。
对于足够大的k,我们有
。因为
,所以
情形2.5
,即
。
由(4)式知此时
。若
,当k足够大时,
,
,因此
若
,当k足够大时,
,此时
,因此
综合上述两种情况,我们有
。
情形2.6
,即
。
由(5)式知此时
。若
,当k足够大时,
,
,因此
若
,当k足够大时,
,此时
,因此
综合上述两种情况,我们有
。
综上所述,当
时,总有
,故
3.2.
的强度量次正则
为了研究
的强度量次正则,设
,
,f在
处是二次可微的。定义:
其中
。
定理1 如果
,则
在
处是强度量次正则的。
证明因为
,f在
处是可微的,所以有
。由引理2可知:
设
,
,有
由命题2,当
时,
;当
时,有
或者
当
时,有
;当
时,有
;当
时,有
因为
,所以有
,这意味着
再由引理1知
在
处是强度量次正则的。
为了建立
强度量次正则的充要条件,我们引出如下记号:
定理2
在
处是强度量次正则的当且仅当
证明 因为
,f在
处是可微的,所以有
。由引理2可知:
(7)
由引理1知
在
是强度量次正则的当且仅当
,即
(8)
由命题2和(7),我们有(8)成立当且仅当
注意到
当且仅当
。
4. 结束语
在本文中,我们分别建立了带MCP函数正则化问题次微分强度量次正则的充分条件和充要条件。首先,我们根据次微分的定义得到了MCP函数的次微分的表达式,并在此基础上求得了MCP函数次微分的图像导数。最后,通过求和准则将MCP函数次微分的图像导数与正则化问题(1)次微分的图像导数建立联系,得到了正则化问题(1)次微分强度量次正则的充分条件和充要条件。后续我们将研究正则化问题(1)的二阶增长条件,并用强度量次正则或二阶增长条件研究该模型相关算法的收敛性。
基金项目
本文在重庆市自然科学基金(CSTB2022NSCQ-MSX0409, CSTB2022NSCQMSX0406, CQ-NCAM-2021-02)和重庆市教育委员会科技项目(KJZD-M202201303, KJQN202201349)资助下完成。