LINGO大规模规划求解.docx

文档介绍：lingo大规模规划求解首先,让我们先看看一个非常简单的规划例子在LINGO软件中实现过程:目标函数:y(min)=2*x1+x2-3*x3+5约束条件:x1+x2-3*x3≤10x1-2*x2≥5x1,x2,x3∈N求解上面目标函数的最小值,我们在lingo中可编写如下代码:model:MIN=2*X1+X2-3*X3+5;X1+X2-3*X3<=10;X1-2*X2>=******@GIN(X1);!整数约束;***@GIN(X2);***@GIN(X3);END可以看出,LINGO语言和数学专业语言很接近,很容易表示约束条件和目标函数。可是对于规模很大的约束条件,难道我们也必须这样一条一条的输入吗,显然这样做是一件非常困难和繁琐的事,lingo语言又是如何表示约束条件规模巨大的规划问题呢,带着这样的疑问,让我们一步一步得看下面的内容:一、集合域在数学中集合的定义如下:集合:具有某种相同属性的对象放在一起,就形成了一个集合,集合中的每一个对象称作该集合的元素。在大规模的优化问题中,集合是必然存在的,比如在平板车建模中,各种规格集装箱就可以看成一个集合,在货物配送问题中154个城市可以看成一个集合。在lingo语言中,将某些对象看成一个集合便可以很方便地对集合中的每一个元素进行统一处理。集合域必须在模型的约束引用集合之前定义。集合域用关键字“sets”开始,“endsets”结束。集合分类:基本集合定义统一语法格式:setname[/member-list/][:attribute-list];集合名/对象名1对象名2…对象名n/:对象属性;集合定义的几种方法:setsrow/1..20/:d1,d2,…dn;!集合名/对象名/:对象属性;endsets派生集合定义方法:setsrow/1..20/;col/1..100/;page/1..50/;link():k1,k2…kn;trd/():t1,t2,…tn;endsetsK1可以表示某个省的某个城市的人口。t1可以表示某个省的某个城市某个人的收入。二、数据域:数据域是优化问题中已知得对象的属性值,例如:人的身高,体重;车辆的载重,行驶速度。数据域以关键字“data”开头,“enddata”结束。数据域可以出现在模型中的任何地方。data:d1=…;k1=…;t1=…;enddata数据域的未知数值有时只想为一个集的部分成员的某个属性指定值,而让其余成员的该属性保持未知,以便让LINGO去求出它们的最优值。在数据声明中输入两个相连的逗号表示该位置对应的集成员的属性值未知。两个逗号间可以有空格。sets:years/1..5/:capacity;endsetsdata:capacity=,34,20,,;enddata属性capacity的第2个和第3个值分别为34和20,其余的未知。实时数据处理在某些情况,对于模型中的某些数据并不是定值。譬如模型中有一个通货膨胀率的参数,我们想在2%至6%范围内,对不同的值求解模型,来观察模型的结果对通货膨胀的依赖有多么敏感。我们把这种情况称为实时数据处理(whatifanalysis)。LINGO有一个特征可方便地做到这件事。在本该放数的地方输入一个问号(?)。data:interest_rate,inflation_rate=.085?;enddata每一次求解模型时,LINGO