2017年数学建模讲座文件01第1讲lingo软件入门

1、1入：首先介绍两个简单的 LINGO 程序。例1入：首先介绍两个简单的 LINGO 程序。例z72x1 x1 x2 12x 8x 1 例 抛物面zx2y2被平面xyz1截成一椭圆，求原点到这椭圆的最短, x2 y2 zx y z z x2 y2LINGO 求解程序如下： free(x);free(y);说明：LINGOxy例1 2xixi 2xixi x2遍，free(x1); ; free(x99)。1.1 LINGO的模型窗 LINGO 程序中，每一个语句都必须要用一个英文状态下的分号结束，一个语 2 LINGO程序中“或或=0 x 的最大整数；当x=0 x 的最大整数；当x0时，返回不低

2、于 x 的最大整数。 bin(x)x01； 在默认情况下，LINGO 0，上界为+。free 取可以取消默认下界为 0 的约束。 或电子表格（EXCEL）文件中。也就是说，只能输出到TEXT函数，而不能输出到OLE和ODBC 函数。 例使用所示8运销A4953858252197433767392712395726555228143xij(i1,2,6; j 1,2,8)cij 表示产AiBjAi Bjxij(i1,2,6; j 1,2,8)cij 表示产AiBjAi Bj ij c x i1 j6djj 8eii ji1,2,6; j 6 1.5中的目标函数表达式min cijxij LING

3、Oi1 j 式中，sum LINGO 6 当于c ij i1 j如果表达式中参与运算的属性属于同一个集合，则sum语句中索引（相当于矩性的索引。 的目标函数中的属性 c 和x 属于同一个集合，因而可以表示成：81.5中的约束条件xij ei,i 1,66LINGOj 语句中的for LINGO 该集合的所有成员生成对应的约束表达式，上述for warehouses，它表示66个约束表达式；for的第二个参数是约束表达式的具体内容，此处再调用sumvendors8个成 j8 合，所以不能省略索引i，j96同样地，约束条件xij dj, j 1,8用LINGO语句表示为 6同样地，约束条件xij

4、dj, j 1,8用LINGO语句表示为 m:问题warehouses/1.6/: e; vendors/1.8/: d; links(warehouses,vendors):c,x; e=605551 4341d=35372232 413243 c=6 267 4 29 4 95 3 85 85 2 72 39 5 72 65 52 2 81 4 即把矩阵的元素逐行排成一个行向量的格式，上面LINGO语句中为了容易看清属性cLINGO函数file从外部纯文本文件读入数据，执行一次file1个计算的 LINGO 程序如下：warehouses/1.6/: e; vendors/1.8/: d;

5、 links(warehouses,vendors):c,x;min=sum(links: c*x); for(vendors(J):sum(warehouses(I):x(I,J)=d(J); ，text(sdata12.txt)=table(x!把计算结果以表格形式输出到外部纯文本文件; 645725487978text(sdata12.txt)=table(x!把计算结果以表格形式输出到外部纯文本文件; 645725487978292657393529837523153运价是（3）通过Excel文件中导入数据，也能把计算结果写入Excel文件。从 Excel 文件中导入数据的调用格式如下

6、： 1 warehouses/1. 6/: e; vendors/1.8/: d; links(warehouses,vendors):c,x;min=sum(links: c*x); 列表range_name_list也可以替换为Excel oleobject_list; object_listole warehouses/1.6/: e; vendors/1.8/: d; links(warehouses,vendors):c,x;warehouses/1.6/: e; vendors/1.8/: d; links(warehouses,vendors):c,x;): 问题的求解结果的LI

7、NGO 子模:( 例当参数a0,1, 2,3, 4；b2, 4,例当参数a0,1, 2,3, 4；b2, 4,6,7时，分别求下列的非线性规划问题min 4x3 ax 2x x1 x2 2x x s.t.x bx 解a的取值有5种可能b 的取值有 4种可(ab的取值总共有20种组合，这需要求20 个非线性规划问题。利用LINGO的子模型功能，只要编写一LINGO程序就计算的LINGO var1/1.5/:a0;!a0用于存放a的取值; var2/1.4/:b0; !b0用于存放b的取值; var3/1 2/:x;a0=01 2 3b0=24 6 obj min=4*x(1)3-a*x(1)-2

8、*x(2); !为了下面!调用目标函数和约write(a=,a0(i),b=,b0(j), 时 ， 最 优 解 x1=,x(1),x2=,x(2), 最 优 值 为 ,obj, 1.4 LINGO1.1 求解线性规划问题z 72x1 x1 x2 12x 8x 11.2 min Z120 x190z 72x1 x1 x2 12x 8x 11.2 min Z120 x190 x280 x370 x430 x5x1x2x5 x3 x4 s.t. 3x 2x 3x 2x x xj0, j,1.3A、设B 的台时、调试工序时间及每天可用于这两种家电的能力、各销售一件时的获利情况，如表 1.2 所示。问该公司应制造两种家电各多少件，使获取的利润为最大。7378857597896466937758