Lingo的基本用法

LINGO的基本用法一 集合的基本用法集合(set)及其属性(attribute)的概念基本集合与派生集合集合名/元素列表/：属性列表；集合名（父集合列表）/元素列表/：属性列表；稠密集合与稀疏集合元素过滤法基本集合的隐式列举法：数字型1.n字符数字型Car101.Car208日期型MON.FRI月份型OCT.JAN年月型OCT2007.JAN2008二 模型结构（1） 集合段 从“sets:”到“endsets”（2） 数据输入段 从“data:”到“enddata”属性常数列表（3） 目标和约束段MIN表达式（4） 计算段 从“calc:”到“endcalc”，对原始数据的计算处理（5） 初始段 从“init:”到“endinit”，定义迭代初值用属性常数列表（6） 注释 从感叹号到分号三 函数基本数学函数ABS(X)COS(X)EXP(X)FLOOR(X)LGM(X)LOG(X)MOD(X,Y)POW(X,Y)SIGN(X)SIN(X)SMAX(list)SMIN(list)SQR(X)SQRT(X)TAN(X)其中 LGM(X) ln(X-1)!集合循环函数FORMIXMINPRODSUM用法： 集合函数名(集合名(集合索引列表)|条件：表达式组）集合操作函数ININ(集合名，集合元素名， 集合元素名)INDEXINDEX(集合名，集合元素名)WRAPWRAP(i，N)，循环计数SIZESIZE(集合名)变量定界函数BND(L,X,U)BIN(X)FREE(X)GIN(X)分别对变量取值限制：上下界，0-1值，取消非负限制，整数概率分布函数PNS(X)标准正态分布PSL(X)正态线性损失PBN(P,N,X)二项分布PHG超几何分布PTD(N,X)t分布PFD(N,D,X)F分布PPS(A,X)泊松分布PPL(A,X)泊松线性损失PCX(N,X)X平方分布RAND(seed)随机数服务系统函数PEL(A,X)PFS(A,X,C)PEB(A,X)文件输入输出函数FILE(fn)TEXT(fn)OLE结果报告函数ITERS()返回迭代次数NEWLINE(n)输出n个新行STRLEN(string)返回字符串的长度NAME(reference)返回变量名或行名WRITE用于数据段，输出变量，字符串或换行WRITEFOR是WRITE在循环情况下的推广FORMAT以格式描述符方式输出数值DUAL(varname)返回解答中变量的判别数或结束行的影子价格STATUS()返回求解后的最后状态其他函数IFIF(条件，true结果，false结果)WARNWARN(text，条件)USERUSER(用户编写的函数dll或obj文件)四 文件传输通过文本文件传输数据FILE和TEXT通过Excel文件传输数据OLE例题! 背包问题 Knapsack Problem! max z=sum(i=1n)ci xi! st. sum ai xi =b, xi=0/1;model: title 背包问题; sets: wp/w1.w8/:a,c,x; endsets data: a=1 3 4 3 3 1 5 10; c=2 9 3 8 10 6 4 10; enddata max=sum(wp:c*x); for(wp:bin(x); sum(wp:a*x)=15;end!装箱问题！=;model: title 装箱问题; sets: wp/w1.w30/:w; xz/v1.v30/:y; links(wp,xz):x; endsets data: w=0.51,0.51,0.51,0.51,0.51,0.51, 0.27,0.27,0.27,0.27,0.27,0.27, 0.26,0.26,0.26,0.26,0.26,0.26, 0.23,0.23,0.23,0.23,0.23,0.23, 0.23,0.23,0.23,0.23,0.23,0.23; enddata min=sum(xz(i):y(i); C=1; for(links:bin(x); for(xz:bin(y); for(wp(i):sum(xz(j):x(i,j)=1); for(xz(j):sum(wp(i):w(i)*x(i,j)3 ) ;z2= z1+sum(o(i):x(2)(i) - sum(o(i):he(2)(i) ;z3= z2+sum(o(i):x(3)(i) - sum(o(i):he(3)(i) ;z4= z3+sum(o(i):x(4)(i) - sum(o(i):he(4)(i) ;z5= z4+sum(o(i):x(5)(i) - sum(o(i):he(5)(i) ;研究总利润和采购与加工方案适应不同的未来市场价格应如何变化。考虑如下的价格变化方式：2月份植物油价上升x%，非植物油上升2 x%；3月份植物油价上升2x%，非植物油上升4 x%；其余月份保持这种线性的上升势头。对不同的值x（直到20），就方案的必要的变化及对利润的影响，作出全面计划。for(o(j):cun(1,j)= 500+x(1,i)-he(1,j) ;for(o(j):cun(2,j)= cun(1,j)+x(2,j)-he(2,j) ;for(o(j):cun(3,j)= cun(2,j)+x(3,j)-he(3,j) ;for(o(j):cun(4,j)= cun(3,j)+x(4,j)-he(4,j) ;for(o(j):cun(5,j)= cun(4,j)+x(5,j)-he(5,j) ;for(o(j):cun(6,j)= cun(5,j)+x(6,j)-he(6,j) ;for(month(i):he(i,1) +he(i,2)200);for(month(i):he(i,3) +he(i,4)+he(i,5)250);for(links(i,j): cun(i,j) 1000 ) ;for(o(j):cun(6,j)=500 );for(month(i) : sum(o(j): y(j)*he(i,j) 3*sum(o(j):he(i,j) );飞行管理问题在约一万米高空边长160公里的长方形区域内有若干架飞机（最多6架）在水平飞行。当一架飞机新进入该区域，地面调度中心根据记录的数据（各架飞机的位置和速度向量），要立即计算并判断是否会发生飞机相撞。如果会相撞，则应计算如何调整各架飞机飞行的方向角，以避免相撞。不相撞的标准为飞机间的距离大于8公里，所有飞机的航速均为800公里/小时，新进入的飞机与其它飞机的距离都在60公里以上。试建立一个避免碰撞的数学模型，要求方向调整角的幅度尽量小，且不超过30度，计算误差不超过0.01度。设飞行区域的四个顶点的坐标为（0, 0）、（160, 0）、（160, 160）和（0, 160）。当前飞行记录为（6号为新进入的飞机）飞机编号横坐标X纵坐标Y方向角（度）12345615085150145130014085155501500243236220.515923052! 飞行管理;model:sets:feiji/p3 . p6/:zt,zt0,dzt,zt1,xi0,yi0;juli/l1 . l3001/:L;endsetsdata:xi0=150,145,130,0;yi0=155,50,150,0;zt=220.5,159,230,52;enddatafor(juli(i):L(i)=90+(i-1)/100);for(feiji:bnd(-0.08,dzt,0.08);for(feiji:zt0=zt*3180);for(feiji:zt1=zt0+dzt);fo