1. 引言
本文考虑如下约束向量优化问题(constraint vector optimization,简记CVOP):
C是
中的非空闭凸子集,
是连续可微的,且
,
,
,
是连续可微的意味着它的每一个分量
都是连续的。这种类型的问题在金融、环境分析和管理学等领域都有所涉及,见文献 [1] [2] [3] 。
在求解向量优化问题的众多方法中,主要分成两类:标量化方法和下降方法。
标量化方法的基础思想是通过特定的参数将向量值问题转化为标量值问题。其中,著名的加权求和方法是通过对各个目标函数进行线性组合来得到一个单目标优化,新问题的最优解亦可看成是原始问题的最优解,如 [4] 。
向量优化的下降方法不需要提前知晓任何的参数信息,而且关于收敛性有理论的证明。它们是从单目标非线性规划的算法中拓展而来的。当向量优化无约束的时候,即如下形式:
有下面几种方法被提出:
Jörg Fliege和Benar Fux Svaiter [5] 提出解决向量优化的最速下降法,下降方向定义如下:
在解决约束向量优化问题时,目标函数是Lipschitz连续的条件下,证明了该算法收敛到一个满足某些一阶必要条件的点。
Fliege, J.等人 [6] 提出解决向量优化的牛顿法,下降方向定义如下:
文中也证明了在目标函数满足二阶连续可微的条件下,该算法生成的序列点是局部超线性收敛到最优点的。
对于带约束的向量优化问题,2004年,Graña Drummond和Iusem [7] 提出了解决向量优化的投影梯度算法(projected gradient method,简记为PGM),下降方向定义如下:
其中,
表示的是F在x处的雅可比矩阵,
,
。上面几种方法可以通过选择适合的初始点来找到帕累托最优点,并且它们都是迭代型的算法,每次迭代的时候选择的下降方向都是通过解决一个目标函数是强凸的单目标优化来得到的,这样保证了下降方向是唯一的。特别地,Graña Drummond和Iusem也就投影梯度算法的两个方面给出了收敛性分析。一个是加上平方和步长和凸假设,他们证明了算法收敛到一个弱帕累托解。另一个则是给定目标函数是凸函数的假设,但是没有额外加步长,此时,他们证明了聚点是驻点。作为一种经典而有效的方法,投影梯度算法被广泛研究和应用在各种多目标规划中。关于固定步长的投影梯度算法,也有一些重要的工作。针对非可微凸向量优化问题,Bello Cruz等人 [8] 提出投影次梯度算法,并且证明了由算法生成的序列的聚点是弱帕累托解。最近,通过使用标量化的方法,Brito等人 [9] 和Bento等人 [10] 使用带权重的次梯度算法来解决非可微凸向量优化问题,并且获得了收敛到帕累托解的结果。虽然上述方法没有使用单调线搜索技术,但是它们依然有着很强的收敛结果。
在迭代过程中,改变步长能否对算法进行改进是一个值得研究的问题。本文受此启发,引入非单调步长来得到非单调投影梯度算法,并在合适的条件下证明算法的有效性。
2. 相关定义和引理
定义2.1 集合C被称为是锥,如果对任意的
和
,都有
.进一步,集合C被称为是尖锥,如果它是锥,并且满足
.
定义2.2 令集合K是非空闭尖凸锥,偏序“
”的定义如下:
定义 2.3 集合
称为集合K的极锥,如果它满足:
定义2.4 设集合C是一个非空的集合。
1)
被称为是帕累托最优,如果不存在
使得
和
2)
被称为是弱帕累托最优,如果不存在
使得
3)
被称为是F的驻点,如果满足
其中,
。
很多情况下,(弱)帕累托解很难被证明,许多学者会考虑证明生成的序列收敛到驻点来取代证明收敛到(弱)帕累托解。根据驻点的定义,我们知道如果
不是驻点,那么,存在
使得
,这等价于
,取
时,这进一步等价于如下表达式:
又F是连续可微的函数,故
此式说明
在
处是一个下降方向,即不是驻点就存在下降方向,反之,不存在下降方向的点就是驻点。
定义2.5 映射F在H上是L-Lipschitz连续映射,当且仅当存在常数
使得映射F满足
如果
,则称F是非扩张的;如果
,则称F是压缩的。
定义2.6 (收敛速度)设序列
是
中的序列,并且满足依范数收敛到
,称序列
具有Q-线性收敛率,如果存在
使得当下标k充分大时,有
下面给出一些引理用于后面算法的收敛性证明。
引理 2.7 [11] 令C是一个锥,则对于任意的
,有
.
引理2.8 [12] 对于问题(CVOP),如果
是帕累托最优点,则x也是驻点。
3. 算法及其收敛性
3.1. 算法描述
步骤1初始化:给定初始点
,参数
,
,
,
以及非负整数M和N,令
和
;
步骤2确定寻找方向:计算
(3-1)
步骤3停止准则:如果迭代次数
或者
,停止;否则,继续步骤4;
步骤4令
和
:
(3-2)
步骤5确定步长:取
,其中
是使得下面不等式成立的最小非负整数:
(3-3)
步骤6令
,并返回到步骤2。
注3.1:当
是由算法1生成的序列时,由(3-2)可以推导出
。
首先,定义函数
如下:
(3-4)
结合极锥
和集合G的定义,集合
和
可以转化下面两个形式:
(3-5)
和
(3-6)
下面一个引理说明函数
的一些性质。
引理3.1 ( [7] Proposition 2)
的定义如(3-4)所示,下面几个性质成立:
1) 函数
具有一阶正齐次性,即
;
2) 对于所有的
,
和
都成立;
3) 对于所有的
,若
(或者
),则
(或者
);
4) 函数
是Lipschitz连续的,且Lipschitz常数
。
其次,将(3-1)中目标函数记为
:
(3-7)
其中,
,且
是使得
取值最小的一个向量。
下面的引理则陈述关于此函数的一些性质,这些性质对于算法的收敛性证明非常有用。
引理3.2 ( [7] Proposition 3) v和
的定义分别如(3-1)和(3-7)所示,下面几个结论成立:
(1)
;
(2) 当
时,下面几个条件等价:
(a) x不是F的驻点;
(b)
;
(c)
.
注3.2:(1) 由命题3.2(1)和
的定义知
,由
的定义知
,从而知
是一个下降方向;
(2) 从引理3.2(2)可以推导出
是F的驻点
.
在步骤3的停止准则中,其中一个是要求
。一般来说,参数
是任意小的,所以从注3.2(2)可知此条件意味着
是驻点。
命题3.3对于给定的点
,
定义如(3-1),则它满足下面两个不等式:
(3-8)
(3-9)
证明:第一个不等式根据
的定义和命题3.2(1)可以得到。再利用
的定义和柯西-施瓦茨不等式,有
将此与第一个不等式相结合即可得到第二个不等式。
下面的引理给出了函数
的连续性性质。
引理3.4 ( [7] , Proposition 4) 函数
在C中是连续的。
下面的引理说明极大型线搜索策略和投影梯度算法结合是可行的。
引理3.5 (1) 算法1是有意义的;
(2) 令
是算法1生成的序列,则
。
证明:(1) 显然函数
是一个强凸函数,则下降方向
可以得到且唯一。由注3.1可知(3-3)是有意义的。故结论成立。
(2) 下面用数学归纳法来证明结论。
当下标
时,由算法1的初始化步骤可知
,结论成立;
假设当下标
时,
,由下降方向
的定义知,
,而
,
以及
和C是一个闭凸集合,所以
结论成立。
3.2. 全局收敛
命题3.6 假设F连续可微有下界,序列
是由算法1生成的,则序列
是单调非增的,并且存在极限。
证明:因为
是一个下降方向,故
对所有的
均成立,将其带入到(3-3)中,有
。
当序列
是由算法1生成的序列时,根据(3-2)和
,对于向量的每一个分量i和迭代次数k,有
这也就说明了序列
是单调非增的。同时考虑到
和F是有下界的,根据单调有界定理知序列
存在极限。证毕。
引理3.7 假设F有下界且连续可微,序列
是由算法1生成的,那么,
证明:对于任意迭代次数k和
,令
为使得下面两个条件成立的整数:
(3-10)
(3-11)
利用(3-3)和
的迭代公式,有
(3-12)
对于条件(3-10),可得
,当k足够大时
取常数,再根据命题3.6,有
又因为
和
,结合(3-12)可得
(3-13)
因为
是任意取的,进一步可得
从而
(3-14)
定义
显然
,在证明结论之前,先证明下面两个结果:对任意给定的
,有
(3-15)
(3-16)
接下来用数学归纳法来证明这两个结果。
当
时,结合(3-14)与
,可得(3-15)。考虑到(3-9)和
的迭代公式,有
(3-17)
因为
,在等式两端同时对k取极限,并将(3-15)带入可得
又
是连续函数,故
后一等式的成立利用了
的定义(3-2)。因此,(3-16)成立。
假设当
时,(3-15)和(3-16)成立。利用(3-3),有
(3-18)
注意到
和(3-18)在
时成立,利用命题3.6,
,
和
,由(3-18)可以推导出
由i的任意性有
利用命题3.1,上式可以进一步变为
(3-19)
则(3-15)成立。
根据命题3.3和
的迭代公式,有
将(3-19)与上式相结合,可推导出
又因为
是连续函数,故
于是,(3-15)和(3-16)对任意给定的
都成立。
对所有的迭代次数
,有
(3-20)
注意到
,利用命题3.3和式(3-15)、(3-20),整理有
再一次用到
的连续性,得
,将此带入到不等式(3-3)中即可得到结论。证毕。
定理3.8 假设F是有下界且连续可微的。序列
是算法1生成的序列,如果序列
存在聚点,那么此聚点是F的驻点。
证明:记
是序列
的聚点,因为
以及集合C是一个闭凸集,故
。同时,序列
存在一个子序列
收敛到
,其中
是一个指标集。由引理3.7可以得到
(3-21)
此式可以拆解成两种情况:
(1)
;
(2)
。
当情况(1)成立时,由引理3.4知 是连续的,则
。又
则
再根据注3.2(2)可知
是驻点。
当情况(2)成立时,存在一个指标
使得
对所有的
和
成立。
是满足条件(3-3)的步长,故
也就是
(3-22)
而左边又可以写成如下形式:
(3-23)
又因为
,即
而K是一个锥,结合引理2.7可得
(3-24)
根据中值定理,可知存在
满足
将此带入到(3-24)中,并在两边同时除以
,整理得
(3-25)
序列
是一个有界序列,则该序列存在一个收敛的子序列,即存在指标集
使得
利用
,有
在(3-25)两端同时除以
,并对
中的k取极限,可得
考虑到
,故
,又
是一个下降方向,故
,即
,而K是一个闭凸尖锥,可得
,从而
,再利用
的定义,有
另一方面,由命题3.2(1)知
,故
,因此
是F的驻点,证毕。
3.3. 线性收敛
本小节假设以下条件成立:
(A1)
是L-Lipschitz连续的;
(A2) 存在一个正数c满足下面的不等式成立:
(A3) 令
是F的驻点,存在正常数
和
满足:
下面的引理证明了极大型和平均型非单调线搜索生成的步长是有下界的。
引理3.9 假设(A1)成立,则对于任意的迭代次数k都存在一个正常数
使得
。
证明:任取迭代次数
并将其固定,再根据
的相关定义可得
故存在
使得
(3-26)
第二个不等式的成立是依据(3-2)和注3.1和注3.3。再利用(A1)和F是连续可微的函数的条件,可得
将上式与(3-26)相结合,可得
(3-27)
考虑到
和
,再将(3-27)和(3-8)相结合可得
这也就证明了结论。
在单目标规划中,Y.H. Dai等人 [13] 证明了极大型线搜索算法的线性收敛性质,其中要求目标函数是光滑和一致凸。受此启发,将这一性质扩展到多目标规划中并建立相关的线性收敛性,依然在目标函数是连续可微和一致凸的条件下。下面是关于算法1的线性收敛性分析。
定理3.10 设条件(A1)、(A2)和(A3)成立,则由算法1生成的序列
满足:存在
使得
(3-28)
其中
表示的是序列
的聚点。
证明:(3-28)等价于下面的不等式:
(3-29)
在证明结论之前,先证明下面的不等式对所有的
都成立:
(3-30)
下面用数学归纳法来证明:
当
时,(3-30)可以依据(3-3)得到。
假设(3-30)在
处成立,考虑到
、
和
,有
(3-31)
将(3-3),(3-2)以及
的定义与上式相结合,可得
即可证明了(3-30)对所有的
都成立。在式(3-30)两边同时关于
取最大值,该式即可变成
利用假设条件(A2)和引理3.9,可推导出
其中
,再在上式两端同时减去
,可得
(3-32)
此时,将证明(3-29)转变成证明下面的不等式:
(3-33)
如果上述不等式成立,利用注3.1,有
这一不等式也就是(3-29)。
下面定义如下的两个数:
依据(3-2),有
的定义如(3-10)和(3-11).
下面分成两种情况来证明(3-33):
(1)
将此带入到不等式(3-32)中即可得到(3-33);
(2)
由假设条件(A3)和
的相关定义可以获得
再将此与假设条件相结合即可得到(3-33)。证毕。
4. 结论
本文提出一种新的算法:极大型投影梯度算法,具体而言,将解决约束向量优化问题的投影梯度算法与极大型非单调线搜索技术相结合,这样做的目的是为了摆脱单调线搜索技术陷入局部收敛的可能。在适当的条件下,证明所提算法确实不再是局部收敛的,而是全局收敛的。进一步,加强条件:目标函数各分量的梯度函数是满足Lipschitz条件的,并且和下降方向与目标函数满足某些关系。在此基础上证明了所提算法是线性收敛的。