1. 引言
最优化算法是我校数学系数学与应用数学专业和信息与计算科学专业的一门专业课。这门课包括线性规划、非线性规划、多目标规划和现代智能算法等内容 [1] [2] 。在非线性规划中,极小极大问题是一类优化模型 [3] [4] ,它的目标函数一般形式为
.
例如经典的Sylvester问题(第4页, [1] ),是求覆盖平面上m个点的最大圆盘。我们设这m个点为
,平面上任一点x到各点的距离则分别为
。若以
表示x到m个点
的最大距离,即
,则以x为圆心,以
为半径的圆盘能覆盖这m个点。于
是问题转化求解
,
这是一个无约束的极小极大问题。
对于多目标规划 [5] ,在决策时采取稳妥的保守策略,即在最坏的情况下,寻求最好的结果。按照此
策略,可以构造评价函数
,然后求解
.
极小极大问题涉及非线性规划和多目标规划两部分内容,对于学生来说,透彻理解该问题是非常有益的。问题本身相对抽象,直接套用Matlab优化工具箱中的fminimax命令却能直接得出结果,所以学这部分内容时容易囫囵吞枣,做不到融会贯通。因此在最优化算法教学中,非常有必要给出容易理解且印象深刻的极小极大问题,于是我们设计了如下的植树问题。
植树问题:某班有男同学30人,女同学20人,星期天准备去植树。根据经验,一天中,男同学平均每人挖坑20个,或栽树30棵,或给25棵树浇水,女同学平均每人挖坑10个,或栽树20棵,或给15棵树浇水。问应怎样安排,才能使植树(包括挖坑、栽树、浇水)最多?
2. 模型建立
2.1. 多目标规划模型
在该植树问题中,假设男同学中有
人挖坑,
人栽树,
人浇水,女同学中有
人挖坑,
人栽树,
人浇水。按照规则,挖坑、栽树和浇水三个步骤都完成才表示成功植树1棵,所以我们可以把每个步骤单独看成一个目标,即植树问题转化为一个三目标的多目标问题,目标函数是
,其中第一个目标是挖坑
,
第二个目标是栽树
,
第三个目标是浇水
,
接着我们给出必须要满足的约束条件,包括男同学的人数约束
,
女同学的人数约束
,
其中
都为整数。
我们希望这三个目标的每一个都是越大越好,这样才可以植树最多,并考虑到挖坑和栽树都要比浇水多,可得如下模型
模型一(多目标规划模型)
2.2. 极小极大问题模型
在模型一中,原问题的答案(最多植树的数目)实际上等于第三个目标浇水的最大值,因为只要浇上水说明前面挖坑和栽树两个步骤已经完成了,那么树也植好了。在模型一约束条件中,我们考虑挖坑和栽树都要比浇水多,这表明三个目标之间是有联系,这时我们想到只要找到三个目标的最大值中的最小值,那么它就是满足挖坑、栽树和浇水三个步骤都完成的数量,从而原问题转化为求解如下的极小极大问题。
模型二(极小极大问题模型)
3. 模型的求解
极小极大问题通常可以使用Matlab优化工具箱中的fminimax命令进行求解, 但是通过查阅Matlab帮助文档,fminimax命令能够解决的极小极大问题是有限制的,见图1。
模型二的约束条件含有整数约束,不属于fminimax能够求解的极小极大问题类型,所以此时fminimax命令无法使用。考虑到植树包括挖坑、栽树和浇水三个阶段,不妨给它们定一个共同的下界s,只要让s能取得最大值就能解决原问题。由于挖的坑、栽的树和浇的水都比s多,我们得到如下单目标优化问题。

Figure 1. The type of min-max problem which can be solved by fminimax
图1. 能够使用fminimax求解的极小极大问题类型
模型三(单目标最大值问题模型)
我们使用Lingo软件编程见图2。
得到结果见图3。
我们得到的结果是最多种植345棵树,此时男同学中17人挖坑,1人栽树,12人浇水,女同学中1人挖坑,16人栽树,3人浇水。
4. 结果的验证和分析
4.1. 结果的验证
极小极大问题和极大极小问题是可以如下相互转化的
,
我们将模型二改写成如下等价的极大极小问题
模型四(极大极小问题模型)
但是同样由于约束条件中有整数约束,Matlab优化工具箱中的fminimax命令无法使用,这里我们仍然可以考虑把这个极大极小问题转化为单目标优化问题。此时
是
,
和
的共同上界,可得如下单目标最小值优化问题。
模型五(单目标最小值问题模型)
我们在实际编程中用的目标函数是
,使用Lingo软件编程见图4。
得到结果见图5。
我们得到的结果是−345,它的相反数是345,从而最多种植345棵树。从计算结果看,男同学中有17人挖坑,有1人栽树,有12人浇水,女同学中有1人挖坑,有16人栽树,有3人浇水,这些和原结果都是相同的,原结果得到了验证。
4.2. 结果的分析
我们可以分析出模型三和模型五结果相同的原因。一方面,从目标函数来看,有
,
从而模型三和模型五的目标函数是等价的。另一方面,从约束条件来看,模型三和模型五只有前三个约束条件在形式上是不同的,我们把模型三中的前三个约束条件
,
,
,
两边乘以负一之后,就变成了模型五中的前三个约束条件
,
,
,
所以模型三和模型五的约束条件是一样的。
由上分析可知,因为目标函数等价,约束条件相同,模型三和模型五的结果肯定是相同的。
5. 总结
本文通过实例给出相应的多目标规划问题模型和极小极大问题模型,该极小极大问题含有整数约束,无法简单地应用MATLAB优化工具箱命令进行求解。我们从问题的特点出发,将其转化成一个单目标最大值问题,使用Lingo软件编程求出结果。我们给出与原问题等价的极大极小问题,通过求解其转化成的单目标最小值问题,验证了原结果,并分析了原因。文中的例子简洁明了,学生容易理解。编程简单,能提高学生的学习兴趣。在解决问题的过程中,我们将极小极大问题、多目标优化和单目标优化三部分结合在一起,能给学生很好的启迪,可以解决学生学习极小极大问题无法融会贯通的问题。
基金项目
本文工作得到上海海事大学校精品课程建设项目和院重点课程建设项目(2018)资助。