lingo求解多目标规划 例题

1、实验一、实验目得目标规划就是由线性规划发展演变而来得，线性规划考虑得就是只有一个目标函数得 问题，而实际问题中往往需要考虑多个目标函数，这些目标不仅有主次关系，而且有得还 相互矛盾。这些问题用线性规划求解就比较困难，因而提出了目标规划。熟悉目标规划模 型得建立，求解过程及结果分析。二、目标规划得一般模型设Xj(j1,2.n)就是目标规划得决策变量，共有m个约束就是国内刚性约束，可能就是等式约束，di ,di (i也可能就是不等式约束。设有 l 个柔性目标约束，其目标规划约束得偏差就是1,2,.,l)。设有q个优先级别，分别为P1, P2,. Pq。在同一个优先级Pk中，有不同得权重，分别记为w

2、kj,wkj(j1,2,.,l)。因此目标规划模型得一般数学表达式为:minqPkk1l(wkjdjwkjdj );1s、 t 、naijx jj1( ,)bi,i 1,2,.m,ncij xjdidigi,i1,2,.l,j1xjo,x 1,2,.,n,di ,di0,i1,2,.l.三、实验设备及分组实验在计算机中心机房进行，使用微型电子计算机，每人一机(一组)四、实验内容及步骤1、打开LINGO，并利用系统菜单与向导在E盘创建一个项目。目录与项目名推荐使用学生自己得学号。2、以此题为例，建立数学模型，并用说明语句进行说明，增强程序得可读性。例 2、1 :某工厂生产i、n两种产品，需要用到

3、A B, C三种设备，已知有关数据见下表。企业得经营目标不仅仅就是利润，还需要考虑多个方面:(1)力求使利润不低于 1500元;考虑到市场需求，i、n两种产品得产量比应尽量保持1： 2;设备A为贵重设备，严格禁止超时使用;设备C可以适当加班，但要控制；设备 B即要求充分利用，又尽可能不加班。在重要性上，设备 C就是设备B得3倍。In设备得生产能力/hA (h/件)2212B (h/件)4016C (h/件)0515利润元/件200300解:此题中只有设备 A就是刚性约束，其余都就是柔性约束。首先，最重要得指标就是企1: 2得比例,业得利润，将它得优先级列为第一级；其次就是i、n两种产品得产量保

4、持B得系数就是设备 C得3倍。列为第二级；再次，设备 B、C得工作时间要有所控制，列为第三级。在第三级中，设备该计划问题可用数学模型表示为：目标函数min zp 1d1p 2(d2d2)P3(3d满足约束条件2x1 2x212200x1300x2 d1d115002x1 x2 d2d204x1d3d3165x2 d4d415X1,X2,di ,di0,i1,2,3得重要性就是设备 C得3倍，因此它们得权重不一样，设备LINGO程序为:3 3d3 d4 )model : sets : ! 集合定义部分（从“ sets: ”开始，到“ endsets ”结束）：定义集合变量及其元素（含义类 似数组

5、得下标）与属性（含义类似于数组）。level/1 、 3/:p,z,goal;p , z，与goal 。其中p表示优先级， goal！level 说明得就是目标规划得优先级，有三个变量 表示相应优先级时得最优目标值。！“ 1 、 3”得意思就是从 1 到 3 得所有整数。！基本集合得定义格式为： setname/member_ list/:attribute_list;其中 setname为定义得集合名， member_list 为元素列表， attribute_list为属性列表。在“ ”中得内容，表示就是可选得项，即该项可以有也可以没有。variable/1 、 2/:x;！x 为决策变量

6、向量。h_con_num/1 、 1/:b;）。 b 表示得就是！在目标规划中，约束有两类。一类就是对资源有严格限制得，同线性规划得处理相同，用严 格得等式或者不等式约束来处理，称此约束为刚性约束（hard constraint刚性约束得资源向量。softs_con_num/1 、 4/:g,dplus,dminus; ！另一类约束就是可以不严格限制得，连同原线性规划得目标，构成柔性约束（ constraint ）。 g 表示得就是柔性约束得资源向量， dplus,dminus 就是偏差变量。在目标规 划中，用偏差变量（ deviational variables ）来表示实际值与目标值之间得

7、差异， dplus 为超出目标得差值，称为正偏差变量， dminus 为未达到目标得差值，称为负偏差变量。h_cons（h_con_num,variable）:A;！刚性约束得价值向量。s_cons（s_con_num,variable）:c;！柔性约束得价值向量。obj（level,s_con_num）:wplus,wminus;！柔性约束在不同优先级下得权重。endsets data :！数据输入部分（从“ data ：”开始，到“ enddata ”结束）：作用在于对集合得属性（数组） 输入必要得常数数据。p=? ? ?;！常数列表中得数据之间可以用“，”或者“空格”或者“回车”分开。如

8、果想在运行时才对 参数赋值， 可以在数据段使用输入语句， 但这仅用于对单个变量赋值， 而不能用于属性变量 （数 值）。输入语句格式为“变量名 =？；”。goal=? ? 0;b=12;g=1500 0 16 15;a=2 2;c=200 300 2 -1 4 0 0 5;！LINGO 中得数据就是按列赋值得，而不就是按行赋值得。wplus=0 0 0 00 1 0 00 0 3 1;wmi nus=1 0 0 00 0 3 0;en ddata min =sum(level: p*z);!目标函数(“ min= ”后面所接得表达式)就是用求与函数“sum(集合下标：关于集合属性得表达式)”得方

9、式定义得。这个函数得功能就是对语句中冒号“：”后面得表达式，按照“：” 前面得集合指定得下标(元素)进行求与。这里“sum”相当于求与符号“X”。for (level(i):z(i)=sum(s_c on_n um(j):w plus(i,j)*d plus®)+sum(s_c on_n um(j):wm in us(i,j)*dm in us(j);!约束就是用循环函数“for(集合(下标)：关于集合得属性得约束关系)”得方式定义得。后面得约束关系式都要成意思就是对冒号“：”前面得集合得每个元素(下标)，冒号“：” 立。for (h_con_num(i):sum(variable(

10、j):a(i,j)*x(j)<=b(i);for (s_con_num(i):sum(variable(j)：c(i,j)*x(j)+dm in us(i)-d plus(i)=g(i););for (level(i)|i#lt#bn d(0,z(i),goal(i);!限制 0 =z (i )< =goal (i);!这个限制条件与集合之间有一个 就是一个逻辑表达式，意思就是 于”；size (levelEnd3、下面开始用size (level):)表示集合“ | ”分开，称为过滤条件。< size ( level )。#lt#level 元素得个数。限制条件 “ i#l

11、t#size ( level ) ”就是逻辑运算符号，意思就是“小LINGO中得图标或者Solve命令编译模型，当程序运行时，会出现一个对话框，如图2、1。在作第一级目标计算时,P (1)，P (2)，P ( 3)分别输入 1，0，0,goal( 1 )与goal( 2)输入两个较大得值(例如100000)，表明这两项约束不起作用。运行状态窗口如图2、2,相应信息含义见实验一表1、 1。IVLIIKrO Runt ine InputlOiK图2、1 LINGO 得实时参数窗口LINGO Solver Status l)aok0402Solver atuEHlodelLPVari eblesot

12、>l：StateGiobsL丄 Opt1100>jectivetsibi li ty-CoRs traints(Jtal: vitlEg；iratiflnsSolverNori： 4roE otaL： anliiL 直加r：230BestGerierator Hemory Used. CK)Cbj Aounl;Stipe:nllaped ft-untime (Jih； nin : ss)Activft:OO.OC：42UpdateClose 12、2:LINGO运行状态窗口计算结果如下:Global op timal soluti on foundObjective value:T

13、otal solver iterati ons:000000VariableValueReduced CostP( 1)1、0000000、000000P( 2)0、0000000、000000P( 3)0、0000000、000000Z( 1)0、0000000、000000Z( 2)5、0000000、000000Z( 3)58、000000、000000GOAL( 1)100000、00、000000GOAL( 2)1000000、0、000000GOAL( 3)0、0000000、000000X( 1)0、0000000、000000X( 2)5、0000000、000000B( 1)

14、12、000000、000000G( 1)1500、0000、000000G( 2)0、0000000、000000G( 3)16、000000、000000G( 4)15、000000、000000DP LUS( 1)0、0000000、000000DP LUS( 2)0、0000000、000000DP LUS( 3)0、0000000、0000001DPLUS( 4)10、 000000、000000DMINUS( 1)0、0000001、 000000DMINUS( 2)5、0000000、 000000DMINUS( 3)16、000000、 000000DMINUS( 4)0、00

15、00000、 000000A( 1, 1) 2、0000000、 000000A( 1, 2) 2、0000000、 000000C( 1, 1) 200、 00000、 000000C( 1, 2) 300、 00000、 000000C( 2, 1) 2、0000000、 000000C( 2, 2) -1、0000000、000000C( 3, 1) 4、0000000、 000000C( 3, 2) 0、0000000、 000000C( 4, 1) 0、0000000、 000000C( 4, 2) 5、0000000、 000000WPLUS( 1, 1)0、0000000、 00

16、0000WPLUS( 1, 2)0、0000000、 000000WPLUS( 1, 3)0、0000000、 000000WPLUS( 1, 4)0、0000000、 000000WPLUS( 2, 1)0、0000000、 000000WPLUS( 2, 2)1、0000000、 000000WPLUS( 2, 3)0、0000000、 000000WPLUS( 2, 4)0、0000000、 000000WPLUS( 3, 1)0、0000000、 000000WPLUS( 3, 2)0、0000000、 000000WPLUS( 3, 3)3、0000000、 000000WPLUS(

17、 3, 4)1、0000000、 000000WMINUS( 1, 1)1、 0000000、000000WMINUS( 1, 2)0、 0000000、000000WMINUS( 1, 3)0、 0000000、000000WMINUS( 1, 4)0、 0000000、000000WMINUS( 2, 1)0、 0000000、000000WMINUS( 2, 2)1、 0000000、000000WMINUS( 2, 3)0、 0000000、000000WMINUS( 2, 4)0、 0000000、000000WMINUS( 3, 1)0、 0000000、000000WMINUS(

18、 3, 2)0、 0000000、000000WMINUS( 3, 3)3、 0000000、000000WMINUS( 3, 4)0、 0000000、000000Row Slack or SurplusDual Price1 0、 000000-1、 0000002 0、 000000-1、 0000003 0、 0000000、0000004 0、 0000000、0000005 2、 0000000、0000006 0、 0000000、0000007 0、 0000000、0000008 0、 0000000、0000009 0、 0000000、000000第一级得最优偏差为 0

19、，进行第二级计算。在第二级目标计算时，p（1），p（3）分别输入 0， 1，0。由于第一级得偏差为 0，因此 goal（1）得输入值为 0， goal入一个较大得值（例如100000)。计算结果如下：Global optimal solution found、Objective value:0、 000000Total solver iterations:2VariableValueReduced CostP( 1)0、0000000、000000P( 2)1、0000000、000000P( 3)0、0000000、000000Z( 1)0、0000000、000000Z( 2)0、0000

20、001、000000Z( 3)29、250000、000000GOAL( 1)0、0000000、 000000GOAL( 2)100000、00、 000000GOAL( 3)0、0000000、 000000X( 1)1、8750000、000000X( 2)3、7500000、000000B( 1)12、000000、000000G( 1)1500、0000、000000G( 2)0、0000000、000000G( 3)16、000000、000000G( 4)15、000000、000000DPLUS( 1)0、 0000000、000000DPLUS( 2)0、 0000000、0

21、00000DPLUS( 3)0、 0000000、000000DPLUS( 4)3、 7500000、000000DMINUS( 1)0、0000000、 000000DMINUS( 2)0、0000000、 000000DMINUS( 3)8、5000000、 000000DMINUS( 4)0、0000000、 000000A( 1, 1)2、0000000、 000000A( 1, 2)2、0000000、 000000C( 1, 1)200、 00000、 000000C( 1, 2)300、 00000、 000000C( 2, 1)2、0000000、 000000C( 2, 2)

22、-1、0000000、000000C( 3, 1)4、0000000、 000000C( 3, 2)0、0000000、 0000002），p2）输C( 4, 1) 0、0000000、 000000C( 4, 2) 5、0000000、 000000WPLUS( 1, 1)0、0000000、 000000WPLUS( 1, 2)0、0000000、 000000WPLUS( 1, 3)0、0000000、 000000WPLUS( 1, 4)0、0000000、 000000WPLUS( 2, 1)0、0000000、 000000WPLUS( 2, 2)1、0000000、 000000

23、WPLUS( 2, 3)0、0000000、 000000WPLUS( 2, 4)0、0000000、 000000WPLUS( 3, 1)0、0000000、 000000WPLUS( 3, 2)0、0000000、 000000WPLUS( 3, 3)3、0000000、 000000WPLUS( 3, 4)1、0000000、 000000WMINUS( 1, 1)1、 0000000、000000WMINUS( 1, 2)0、 0000000、000000WMINUS( 1, 3)0、 0000000、000000WMINUS( 1, 4)0、 0000000、000000WMINUS

24、( 2, 1)0、 0000000、000000WMINUS( 2, 2)1、 0000000、000000WMINUS( 2, 3)0、 0000000、000000WMINUS( 2, 4)0、 0000000、000000WMINUS( 3, 1)0、 0000000、000000WMINUS( 3, 2)0、 0000000、000000WMINUS( 3, 3)3、 0000000、000000WMINUS( 3, 4)0、 0000000、000000RowSlack or Surplus Dual Price、 000000-1、 000000、 000000、000000、 0

25、00000、000000、 000000、000000、 7500000、 000000、 000000000000、 000000000000、 000000000000、 000000000000第二级得最优偏差为 0 ，进行第三级计算。在第三级目标计算时，(3)分别输入 0，0，1。由于第一级、第二级得偏差均为 入值也均就是 0。计算结果如下：0，因此 goalp(1)，p( 2)，p(1)、goal (2)得输GGlobal optimal solution foundObjective value:29、00000Total solver iterations:VariableVal

26、ueReduced CostP( 1)、000000、000000P( 2)0、0000000、000000P( 3)1、0000000、000000Z( 1)0、0000000、000000Z( 2)0、000000-5、 666667Z( 3)29、000000、000000GOAL( 1)0、0000000、 000000GOAL( 2)0、0000000、 000000GOAL( 3)0、0000000、 000000X( 1)2、0000000、000000X( 2)4、0000000、000000B( 1)12、000000、000000G( 1)1500、0000、000000G

27、( 2)0、0000000、000000G( 3)16、000000、000000G( 4)15、000000、000000DPLUS( 1)100、 00000、000000DPLUS( 2)0、 0000000、000000DPLUS( 3)0、 0000006、000000DPLUS( 4)5、 0000000、000000DMINUS( 1)0、0000000、 000000DMINUS( 2)0、00000011、33333DMINUS( 3)8、0000000、 000000DMINUS( 4)0、0000001、 000000A( 1, 1)2、0000000、 000000A( 1, 2)2、0000000、 000000C( 1, 1)200、 00000、 000000C( 1, 2)300、 00000、 000000C( 2, 1)2、0000000、 000000C( 2, 2)-1、0000000、000000C( 3, 1)4、0000000、 000000C( 3, 2)0、0000000、 000000C( 4, 1)0、0000000、 000000C( 4, 2)5、0000000、 000000WPLUS( 1, 1)0、0000000、 000000W