LINGO求解优化问题.pdfVIP

首 页 一、Lingo简介 1. 目标函数 一个函数解析式，你希望求它的最大或最小值 max 函数解析式； 或 min 函数解析式； 例 max=3*b+2*c^2; min=b^(1/3)-c*k ； Lingo的语句以 ; 号结束。 2. 运算 加(+) ，减(-) ，乘(*) ，除(/) ，乘方(x^a) 3. 变量 用字母或字母数字的组合表示 例 a ，b ，cc1 ，x1。 Lingo的变量缺省值为非负数。 4. 限制条件 一组等式或不等式。 Lingo的 ，与 ， 等价。 示例程序1 min=3*x1+5*x2; !x1,x2是非负变量; 3*x1+2*x2=36; 3*x1+5*x2=45; Lingo的注释语句用!开头用;结束。 5. 变量类型 变量类型 说 明 @bin( 变量名) ; 限制该变量为0或1。 @bnd( a,变量名, b); 限制该变量介于a ，b之间。 @free(变量名); 允许该变量为负数。 @gin( 变量名); 限制该变量为整数。 二、Lingo高级 sets语句 连续六个月的产量，可以用x1,x2,x3,x4,x5,x6表示, 但十二个月的产量用同样的方法表示就显繁琐。 Lingo可以通过sets语句设置数组功能使问题变得简单。 例 定义数组x, 有x(1),x(2),x(3),x(4)…x(12)个成员，用以表示十二个月的产量。 sets: r/1..12/:x; !r是组的类型名,x数组名; endsets; sets语句以sets开头，endsets结束。 示例程序2 sets: mat/1..4/: x; !mat是组的类型名,x数组名; endsets min=50*x(1)+20*x(2)+30*x(3)+80*x(4); 400*x(1)+200*x(2)+150*x(3)+500*x(4)=500; 3*x(1)+2*x(2)=6; 2*x(1)+2*x(2)+4*x(3)+4*x(4)=10; 2*x(1)+4*x(2)+x(3)+5*x(4)=8; data语句 有时，我们要用到常数数组，比如在 400*x(1)+200*x(2)+150*x(3)+500*x(4)=500 中，x(1), x(2), x(3), x(4)的系数分别为400, 200, 150, 500 ，此时,可用data语句。 例 定义数组a, 其中a(1) 400 ，a(2) 200 ，a(3) 150，a(4) 500。 sets: l/1..4/: a,x; endsets data: a 400 200 150 500; enddata data语句是以data开头，enddata结尾。 示例程序3 sets: l/1..4/: x, a; endsets data: a=7 2 3 9; !a(1)=7, a(2)=2, a(3)=3, a(4)=9; enddata max=x(1)*a(3)+x(2)*a(1)+x(3)*a(4)+x(4)*a(2); x(1)+x(4)-x(2)-x(3)a(1); x(4)+2*x(2)a(4); x(1)+x(3)a(1); Lingo含有一些针对数组的命令，方便了数组的使用。 @for循环语句: @for(数组类型名(i): 循环的语句); 示例程序4