钢管订购笔记.docxVIP

第三次培训钢管订购与运输问题要求：铺设天然气管道制造钢管一个钢厂至少生产500单位钢管（1单位钢管 == 1km钢管）钢厂记为Si，钢厂在指定期限内能生产该钢管的最大数量为个单位，钢管出厂销价1单位钢管为万元钢管运输（公路、铁路）铁路运输1000km以上每增加1至100km运价增加5万元。公路运输公路运输费用为1单位钢管每公里0.1万元（不足整公里部分按整公里计算）。钢管可由铁路、公路运往铺设地点（不只是运到点，而是管道全线连续问题离散化）。图示图中粗线表示铁路，单细线表示公路，双细线表示要铺设的管道(假设沿管道或者原来有公路，或者建有施工公路)，圆圈表示火车站，每段铁路、公路和管道旁的阿拉伯数字表示里程(单位km)。待求解问题问题可能需要的知识：需要将管道送到各自位置处-----连续的问题离散化优化问题-----查看《数学模型》Flody算法（最短路问题）----在matlab或lingo中实现或调用相关算法Lingo ---- 解决优化问题动态规划（多阶段决策过程的最优化问题）（背包问题、最短路径问题）求解思路：路程转化为运输单价求钢厂到火车站点的最小花费路（dijkstra）约束条件：钢厂生产力有限 钢厂分配给各火车站的量的比例、量可知 总运费+订购费最少Lingo使用说明注释 ----- 感叹号!（结尾加分号）每条语句结尾要加分号Lingo sets（集）的概念Lingo把相联系的对象聚合为集sets集分为原始集和派生集语法首先，要先定义集集的名字/成员名字/:属性名字;成员间、属性间用逗号或空格链接sets: 开始endsets 结束例：sets:?students/John??Jill,?Rose??Mike/:?sex,?age;?Endsets隐式罗列成员setname/member1..memberN/:?attribute_list;集成员设置数据data:!代码enddata定义派生集setname(parent_set_list)[/member_list/][:attribute_list];Lingo参考代码Lingo中使用函数Lingo中函数一律用@开头@abs(x);@sin(x); @cos(x); @tan(x);@exp(x); @log(x);@sign(x); x 0 返回-1 否则返回1@floor(x); 返回整数部分@smax(x1,x2,..,xn); 返回x1,x2,..,xn中的最大值@smin(x1,x2,..,xn); 返回x1,x2,..,xn中的最小值@bin(x); 限制x为0或1@bnd(L,x,U); L=x=U@free(x); 取消对变量x的默认下界为0的限制@gin(x); 限制x为整数@if(判断条件,结果为真时执行语句,结果为假时执行语句);注：默认情况下，lingo的变量都是非负的，用@free(x)函数取消限制Lingo逻辑运算符#not# !#eq# ==#ne# !=#gt# greater than#ge# = greaterequal#lt# less than#le# = lessequal#and# #or# ||优先级模型求解状态@status(); ?返回LINGO求解模型结束后的状态：0???Global?Optimum（全局最优）1???Infeasible（不可行）2???Unbounded（无界）3???Undetermined（不确定）4???Feasible（可行）?5???Infeasible?or?Unbounded（通常需要关闭“预处理”选项后重新求解模型，以确定模型究竟是不可行还是无界）?6???Local?Optimum（局部最优）?7???Locally?Infeasible（局部不可行，尽管可行解可能存在，但是LINGO并没有找到一个）8???Cutoff（目标函数的截断值被达到）9???Numeric?Error（求解器因在某约束中遇到无定义的算术运算而停止）????通常，如果返回值不是0、4或6时，那么解将不可信，几乎不能用。该函数仅被用在模型的数据部分来输出数据。!二次规划;model:sets:h/h1..h5/:x,c1,c2;s/s1..s4/:b;link(s,h):a;endsetsmax=@sum(h:c1*x^2)+@sum(h:c2*x);@for(h:@bnd(0,x,99));@for(s(j):@sum(h(i):a(j,i)*x(i))=b(j));data:c1=1 1 3 4 2;c2=-8 -2 -3 -1 -2;b=400 800 200 200;a=1 1 1 1 1 1 2 2 1 6 2 1 6 0 0