1. 引言
模糊优化能更好解决不确定性问题,通过模糊数截集可将其转化为区间值优化问题,为复杂问题求解提供有效路径。而现有模糊优化相关研究中,最优性条件多围绕伪不变凸函数展开,如文献[1],但预不变凸场景下的区间值优化最优性条件尚未得到充分探讨,存在研究空白。本文围绕预不变凸区间值优化问题的最优性条件展开研究。首先梳理了区间的表示方法、运算规则以及序关系,明确了
差等关键概念;其次定义了E-α-预不变凸区间值函数,建立了与E-可微性的等价判定定理;最后构建了E-α-预不变凸环境下的区间值优化问题(IVOPE),说明了可行集、可行方向与有效解的概念,通过理论推导给出了该问题在E-可微条件下满足KKT假设时的最优性充分条件与必要条件,为相关优化问题的求解提供了思路。
2 预备知识
实数集
上所有有界闭区间族记为,给定
,
,定义其标准加法和数乘运算:
显然,
(Minkowski差)。
下面定义给出了区间的LU-序关系。
定义1 [2]令
,则
(1)
且
;
(2)
,且
,即
中至少有一个不等式严格成立;
(3)
且
。
定义2 [2] (
差)对于
,存在
,使得:
3. 凸模糊映射与可微性
定义3 [3]设
为
的一个非空子集,如果
,存在向量值函数
使得
,则称
为关于
的不变凸集。
定义4 [3]设
为关于
的不变凸集,如果
,存在实值函数
,有,
,则称
为关于
的预不变凸实值函数。
定义5 [3]设
为
的一个非空子集,如果
,存在向量值函数
与
,使得
,则称
为关于
的E-不变凸集。
定义6 [4]设
为关于
的E-不变凸集,如果存在实值函数
,对
有
,则称
为关于
的E-预不变凸实值函数。
定义7 [5] (E-α-预不变凸区间值函数)设
为关于
与
的E-α-不变凸集,如果存在区间值函数
,对
,有
,
则称
为关于
与
的E-α-预不变凸区间值函数。
为
上关于
与
的E-α-预不变凸区间值函数,当且仅当
与
是
上关于
与
的E-α-预不变凸区间值函数。
定理1 [5]设
是关于
与
的E-α-不变凸集,
是
上的E-α-预不变凸E-可微实值函数,则对
,有:
定理2 [5]设
是关于
与
的E-α-不变凸集,
是
上的E-α-预不变凸E-可微区间值函数
,则对任意
,有:
定义8 [6]设集合
,映射
,
是
上的区间值函数,
,如果实值函数
与
是可微的,且有:
,
,
其中:当
时,
,
,则称区间值函数
在
处E-可微。
E-可微的几何意义:区间值函数由上下界曲面围成的带状图像,在该点附近可被上下界各自的切平面所夹的线性带状区域局部近似,实现区间值函数的“局部线性化”,是实值函数切线近似的推广。
定理3 [6]设
是
上的E-α-预不变凸函数,
,则
仍为
上的E-α-预不变凸函数。
证明:由于
是
上的E-α-预不变凸区间值函数,由定义8,对
,有:
(1)
(2)
对
,将式(1)乘
、式(2)乘
后相加,得:
将右侧展开并整理,可得:
上式满足定义8中E-α-预不变凸区间值函数的判定条件,故
是
上的E-α-预不变凸区间值函数。
4. 最优性条件
是关于
与
的E-α-不变凸集,在E-可微情况下可建立如下优化问题(IVOPE)
其中
与
均为定义在
上的E-α-预不变凸区间值函数,约束条件
,由
序关系可得
且
。
和
为E-α-预不变凸实值函数,则问题(IVOPE)的E-可行集为
其中
为
的E-α-不变凸子集。
定义9 [7]设
是非空可行集,
。如果存在
,使得
,对于
成立,则称
为在
处的可行方向,
称为可行方向集,记作
若
为在
处的可行方向且
,,则称
为在
处的可行下降方向。易知,若
是(IVOPE)问题的有效解,则函数
在
处没有可行下降方向。
命题1 [8]设
是非空可行集,且
。若
在
处E-可微
,不妨令积极约束指标集为
,
则有,
定义10 [9]如果不存在可行点
,使得
,则称
是问题(IVOPE)的有效解。
KKT假设条件:若
为(IVOPE)问题的可行集且点
。其中
是定义在
上的E-α-预不变凸实值函数,并且对每个
在
处连续E-可微,则称实值函数
在
处满足KKT假设条件。
下面,给出(IVOPE)的KKT最优性充分条件。
定理4 设
是定义在
上的E-可微E-α-预不变凸区间值函数,
在
处满足KKT假设条件。如果存在拉格朗日乘子有
(3)
(4)
(5)
(6)
则
是(IVOPE)优化问题的有效解。
证明:假设
不是问题(IVOPE)的有效解,则存在
,使得
,由
可知有以下三种情况
(i)
,
(ii)
,
(iii)
。
不失一般性,仅考虑(i)
情形。由
是E-α-预不变凸的E-可微区间值函数,结合预不变凸函数的微分性质,有:
(7)
(8)
令
。因为
是E-α-不变凸子集,且
,则对
有
,即
是
处的可行方向。结合命题1,可得
,,
。
对式(5)和(6)两边同乘
,整理得:
,。记
,分下界与上界两种情况分析:
情况1:(下界
的矛盾推导)令
,
,则
,
。由[8]易知,不存在
使得
,结合式(3)的互补松弛条件,
,当
时
,与定理的KKT条件假设矛盾。
情况2:(上界
的矛盾推导)令,,则
,
。同理,不存在
使得
,结合式(3)的互补松弛条件,
,当
时
,与定理的KKT条件假设矛盾。综上,反证的假设不成立,故
是(IVOPE)优化问题的有效解。
接下来,给出(IVOPE)的KKT最优性必要条件。
定理5 设
是定义在
上的E-可微E-α-预不变凸区间值函数,
在
处满足KKT假设条件。如果
是(IVOPE)优化问题的有效解,则存在拉格朗日乘子有
(9)
(10)
(11)
(12)
证明:(反证法),假设存在着一个
,使得
因为
在
处满足KKT假设条件,且
,则存在
,使得
又因为
是K上的E-α-预不变凸实值函数,所以有
令
,那么由于
,通过E-α-不变凸性则有
所以可得
,因此
为
的一个可行方向。设
,则有
因此
为
处的一个可行方向。对于足够大的n,我们有
即
将上述式子分两种情况来讨论:
(1) 当时,显然存在拉格朗日乘子
和
使得其满足(7) (8)。
(2) 当时,显然
是可行下降方向,与
是有效解矛盾。因此,根据文献[8]可知,存在拉格朗日乘子使得
此外,令
,则
综上所述,若
是优化问题的有效解,则有(7)、(8)、(9)、(10)均成立。
5. 数值算例
考虑以下区间值目标规划
.
区间值函数
(
为
上闭区间全体),其上下界函数定义为:
显然
为非传统凸函数,其不变凸集为
,辅助函数为
,
。对
,不妨取
,则:
,由定义8可知
为E-α-预不变凸函数。
由于
在
时单调递减,
时单调递减;约束
下,
处
(下界最小)、
(上界最小),故最优解
。
先看充分条件,要求存在
,满足假设条件:
(13)
(14)
计算得:
,
,
,
,满足(12)。代入(11)得
,解得
,满足假设条件。由定理4知
是最优解。
再看必要条件:若
是最优解,则必存在
满足假设条件。已知
是最优解,易找到
满足条件,故必要条件成立。
若假设
为最优解,
,则有
,解得
,但
,故
非最优解,体现了定理5具有筛选作用。