Lingo例子(学生版).doc

Lingo实验例子 (例子来源：《运筹学教程》（第三版）胡运权主编 清华大学出版社2007年第三版) 一、线性规划及单纯形法、灵敏度分析 例1 P28页例5 Lingo程序： max=2*x1+x2; 5*x2<=15; 6*x1+2*x2<=24; x1+x2<=5; 例2 P44页习题1.7（1）Lingo程序： model: max=2*x1-x2+2*x3; x1+x2+x3>=6; -2*x1+x3>=2; 2*x2-x3>=0; end 其余课本上的例题和习题同学们自己动手编写程序并进行调试运行，分析运行结果。 二、运输问题 例3 P82页例1 Lingo程序： model: sets: gy/g1..g3/:ai; xs/x1..x4/:dj; link(gy,xs):c,x; endsets data: ai=16,10,22; dj=8,14,12,14; c=4,12,4,11 2,10,3,9 8,5,11,6; enddata min=@sum(link(i,j):c(i,j)*x(i,j)); @for(gy(i):@sum(xs(j):x(i,j))=ai(i)); @for(xs(j):@sum(gy(i):x(i,j))=dj(j)); end 例4 P98页例5（转运）Lingo程序： model: sets: plant/x1..x5/:produce; customer/y1..y5/:require; link(plant,customer):c,x; endsets data: produce=60,90,50,50,50; require=50,50,50,80,70; c=-4,5,3,2,100 5,-1,2,100,4 3,2-3,5,5 2,100,5-3,6 100,4,5,6,-5; enddata min=@sum(link:c*x); @for(plant(i):@sum(customer(j):x(i,j))=produce(i)); @for(customer(j):@sum(plant(i):x(i,j))=require(j)); end 三、目标规划 例5 P108页例2 Lingo程序： 第一种做法分三步来完成 第一步：考虑目标：min= model: min=dminus1; 5*x1+10*x2<=60; x1-2*x2+dminus1-dplus1=0; 4*x1+4*x2+dminus2-dplus2=36; 6*x1+8*x2+dminus3-dplus3=48; end 第二步：考虑目标：min= model: min=dminus1+dplus2; 5*x1+10*x2<=60; x1-2*x2+dminus1-dplus1=0; 4*x1+4*x2+dminus2-dplus2=36; 6*x1+8*x2+dminus3-dplus3=48; end 第三步：考虑目标：min= model: min=dminus1+dplus2+dminus3; 5*x1+10*x2<=60; x1-2*x2+dminus1-dplus1=0; 4*x1+4*x2+dminus2-dplus2=36; 6*x1+8*x2+dminus3-dplus3=48; end 注意：如果想要在上述各步中进行灵敏度分析，只需按快捷键Ctrl+R就可以得到相应的结果 或者采用一个程序进行输值运算： model: sets: level/1..3/:P,z,goal; variable/1..2/:x; h_con_num/1..1/:b; s_con_num/1..3/:g,dplus,dminus; h_cons(h_con_num,variable):A; s_cons(s_con_num,variable):C; obj(level,s_con_num):wplus,wminus; endsets data: P=? ? ?; goal=? ? 0; b=60; g=0 36 48; A=5 10; C=1 -2 4 4 6 8; wplus=0 0 0 0 1 0 0 0 0; wminus=1 0 0 0 0 0 0 0 1; enddata min=@sum(level:P*z); @for(level(i):z(i)=@sum(s_con_num(j):wplus(i,j)*dplus(j))+@sum(s_con_num(j):wmi