内容发布更新时间 : 2024/12/25 2:21:17星期一 下面是文章的全部内容请认真阅读。

用LINGO求解整数规划

在LINGO中，输入总是以model:开始，以end结束；中间的语句之间必须以“；”分开；

LINGO不区分字母的大小写；

目标函数用MAX=…；或MIN=…；给出（注意有等号“=”）。 在LINDO中所有的函数均以“@”符号开始，如约束中@gin(x1)表示x1为整数，用@bin(x1)表示x1为0-1整数。

在现在的LINDO中，默认设置假定所有变量非负。 函数中变量的界定：

@GIN(X)：限制X为整数

@BIN(X)：限定变量X为0 或 1。

@FREE(X)：取消对x的符号限制（即可取任意实数包括负数）

@BND(L,X,U)：限制 L<= X <= U

LINGO提供了大量的标准数学函数：

@abs(x) 返回x的绝对值

@sin(x) 返回x的正弦值，x采用弧度制 @cos(x) 返回x的余弦值 @tan(x) 返回x的正切值 @exp(x) 返回常数e的x次方 @log(x) 返回x的自然对数

@lgm(x) 返回x的gamma函数的自然对数 @sign(x) 如果x<0返回-1；否则，返回1 @smax(x1,x2,…,xn) 返回x1，x2，…，xn中的最大值 @smin(x1,x2,…,xn) 返回x1，x2，…，xn中的最小值

例1：整数规划模型在LINGO中可以如下输入： model:

Max=5*x1+8*x2; ！*号不能省略

x1+x2<=6; !约束条件和目标函数可以写在model：与end之间的任何位置 5*x1<=45-9*x2;

@gin(x1);@gin(x2); ！和LINDO不同，不能写在end之后 end

运行后同样得到最优解为x1=0，x2=5，最优值为40。 例2：在线性规划中的应用

max Z =5X1+3X2+6X3, s.t.

X1 +2 X2 + X3 ≤18 2 X1 + X2 +3 X3 =16 X1 + X2 + X3 =10

X1 ,X2 ≥0 , X3 为自由变量

应用LINGO 来求解该模型，只需要在 lingo窗口中输入以下信息即可： max=5*x1+3*x2+6*x3; x1+2*x2+x3<=18; 2*x1+x2+3*x3=16; x1+x2+x3=10; @free(x3);

然后按运行按钮，得到模型最优解，具体如下： Objective value: 46.00000 Variable Value Reduced Cost x1 14.00000 0.000000 x2 0.000000 1.000000 x3 -4 .000000 0.000000

由此可知，当 x1 =14 , x2 =0 , x3 =-4 时，模型得到最优值，且最优值为 46。 说明：在利用LINGO 求解线性规划时，如自变量都为非负的话，在LINGO 中输入的信息和模型基本相同；如自变量为自由变量，可以使用函数 @free来把系统默认的非负变量定义自由变量，如实例一中的 x3。

例3、用LINGO求解整数线性规划问题：

maxz?3x1?4x2?8x3?100y1?150y2?200y3?2x1?4x2?8x3?500?2x?3x?4x?300123??x1?2x2?3x3?100?3x?5x?7x?700123??x1?200y1s.t.??x2?150y2?x3?300y3??xj?0且为整数,j?1,2,3???yj?0或1,j?1,2,31、 模型的输入

使用LINGO求解上述整数规划模型，LINGO程序如下：

MODEL:

max=3*x1+4*x2+8*x3-100*y1-150*y2-200*y3; 2*x1+4*x2+8*x3<=500; 2*x1+3*x2+4*x3<=300; x1+2*x2+3*x3<=100; 3*x1+5*x2+7*x3<=700; x1<=200*y1; x2<=150*y2; x3<=300*y3;

@GIN(x1);@GIN(x2);@GIN(x3); @BIN(y1);@BIN(y2);@BIN(y3); END 2、 执行

点击LINGO菜单下的SOLVE键，或按CTRL+S键，即可求得问题的解。

此问题的解为：x1?100,x2?0,x3?0,y1?1,y2?0,y3?0，最优值为：200。 当运用LINGO求解此问题后，系统会弹出一个名为Solution Report的文本框，其文本框中包含了求解的详细信息，如下：

Rows= 8 Vars= 6 No. integer vars= 6 ( all are linear) Nonzeros= 28 Constraint nonz= 18( 4 are +- 1) Density=0.500 Smallest and largest elements in abs value= 1.00000 700.000 No. < : 7 No. =: 0 No. > : 0, Obj=MAX, GUBs <= 3 Single cols= 0

Global optimal solution found at step: 4 Objective value: 200.0000 Branch count: 0