数学建模讲座优化建模与lingo优化软件.pptVIP

例1 min f(x1,x2)=-2x1-6x2+x12-2x1x2+2x22 s.t. x1+x2≤2 -x1+2x2≤2 x1≥0, x2≥0 Lingo 基 础 Lingo 程序 MIN=-2*X1-6*X2+X1*X1-2*X1*X2+2*X1*X1; X1+X2=2; -X1+2*X2=2; Lingo 基 础 计算结果 Objective value: -9.777778 X1 = 0.6666667 X2 = 1.333333  Lingo 基 础 例2 Lingo 基 础 Lingo程序 min=@exp(x1)*(4*x1*x1+2*x2*x2+4*x1*x2+2*x2+1); x1+x2=0; 1.5+x1*x2-x1-x2=0; -x1*x2-10=0; @free(x1); @free(x2); Lingo 基 础 计算结果 Objective value: 5.276848 X1 = 1.224745 X2 = -1.224745 这是局部最优值 Lingo 基 础 Lingo菜单中Options，选Global Solver ， 在Use Global Solver 中打钩，再次计算。 Lingo 基 础 例三 Lingo 基 础 例三中,x,y ,e是一维数组，长度为2，a,b,d 也是一维数组，长度为6，c是二维数组（6×2矩阵） Lingo中如何定义它们？ Lingo 基 础 Lingo 程序 example3_2 sets: demand/1..6/:a,b,d; supply/1..2/:x,y,e; link(demand,supply):c; endsets data: a=1.25,8.75,0.5,5.75,3,7.25; b=1.25,0.75,4.75,5,6.5,7.75; d=3,5,4,7,6,11; e=20,20; enddata Lingo 基 础 init: x,y=5,1,2,7; endinit min=@sum(link(i,j):c(i,j)*((x(j)-a(i))^2+(y(j)-b(i))^2)^(1/2) ); @for(demand(i): @sum(supply(j):c(i,j)) =d(i);); @for(supply(i): @sum(demand(j):c(j,i)) =e(i); ); @for(supply:@bnd(0.5,X,8.75); @bnd(0.75,Y,7.75); ); END Lingo 基 础 计算结果为 Objective value:89.88349 模型的集部分 LINGO有两种类型的集：原始集(primitive　set)和派生集(derived set)。 一个原始集是由一些最基本的对象组成的。 一个派生集是用一个或多个其它集来定义的，也就是说，它的成员来自于其它已存在的集。 定义一个原始集，用下面的语法： setname[/member_list/][:attribute_list]; 注意：用“[]”表示该部分内容可选。 Member_list是集成员列表。如果集成员放在集定义中，那么对它们可采取显式罗列和隐式罗列两种方式。如果集成员不放在集定义中，那么可以在随后的数据部分定义它们 模型的数据部分 数据部分以关键字“data:”开始，以关键字“enddata”结束。在这里，可以指定集成员、集的属性。 数据部分的未知数值 有时只想为一个集的部分成员的某个属性指定值，而让其余成员的该属性保持未知，以便让LINGO去求出它们的最优值。在数据声明中输入两个相连的逗号表示该位置对应的集成员的属性值未知。两个逗号间可以有空格 数据部分的未知数值 sets: years/1..5/: capacity; endsets data: capacity = ,34,20,,; enddata 属性capacity的第2个和第3个值分别为34和20，其余的未知