lingo软件及其应用(经典)分解.ppt

%matlab_pro; a=[0.50 1.41 3.00 3.37 3.40 4.72 4.72 5.43 7.57 8.38 8.98 9.50]; b=[2.00 3.50 1.50 3.51 5.50 2.00 6.24 4.10 2.01 4.50 3.41 0.80]; a1=a-floor(a); b1=b-floor(b); plot(a1,b1,r.) axis([-1,1,-1,1]); x2=a1=(0.38-0.05)a1=(0.38+0.05) y2=b1=(0.5-0.05)a1=(0.5+0.05) x3=a1([2,4,5,10]) y3=b1([2,4,5,10]) meanx=mean(x3) meany=mean(y3) ex=abs(x3-meanx) ey=abs(y3-meany) （2）LINGO求解模型 s：水平方向网格平移量 t：铅锤方向网格平移量 x1i=xi+s:网格结点横坐标 y1i=yi+t:网格结点纵坐标 决策变量 第k个旧井到第（i,j）个网格的距离0.05 否 目标函数 约束条件 可利用的旧井最多 0.4s0.5,0.5t0.6 e( 4, 6, 5)=1 e( 2, 4, 2)=1 计算结果 e( 4, 4, 4)=1 e( 9, 5, 10)=1 !lingo_pro; sets: xi/1..11/:x,x1; yj/1..8/:y,y1; pk/1..12/:a,b; link2(xi,yj,pk):e; ff/1/:f,g; endsets data: a=0.50 1.41 3.00 3.37 3.40 4.72 4.72 5.43 7.57 8.38 8.98 9.50; b=2.0,3.50 1.50 3.51 5.50 2.00 6.24 4.10 2.01 4.50 3.41 0.80; x=0,1,2,3,4,5,6,7,8,9,10; y=0,1,2,3,4,5,6,7; enddata @for(xi(i):x1(i)=x(i)+x0); @for(yj(j):y1(j)=y(j)+y0); @for(link2(i,j,k):e(i,j,k)=@if(@smax(@abs(x1(i)-a(k)),@abs(y1(j)-b(k)))#le#0.05,1,0)); max=@sum(link2:e); @bnd(0.4,x0,0.5); @bnd(0.5,y0,0.6); end 例21 销售点设置 A1:34人 A2:29人 A5:56人 A3:42人 A4:21人 A6:18人 A7:71人 7个大学生区 分布及人数 一家出版社准备在A1-A7大学生区开设两个销售点,每个销售点只能向本区 和相邻区的大学生售书。这两个销售点应该设在何处，才能使所供应 的学生数量最大？ 设售书点设在i区供应人数为xi 决策变量 第k个售书点设在i区 第k个售书点不设在i区 i=1,2,…,7;k=1,2 目标函数 约束条件 每个售书点只能设在一个区 每个区至多设一个售书点 sets: place/a1..a7/:people_num;!people_num为各区域读书人数; number/shop1..shop2/; links(place,number):e; endsets data: people_num=105,182,126,237,124,166,110; enddata max=@sum(links(i,k):e(i,k)*people_num(i)); @for(place(i):@sum(number(k):e(i,k))=1); @for(number(k):@sum(place(i):e(i,k))=1;); @for(links:@bin(e)); 计算结果 E( A2, SHOP2)＝1 E( A4, SHOP1)＝1 例22 面试问题 (1) 四名不同专业的学生参加招聘面试，面试分三阶段面试：秘书面试、主管复试和经理 面试，在任何一个阶段均按甲、乙、丙，丁顺序，面试时间参看下表（时间单位:分） 秘书面试 主管复试 经理面试 同学甲 13 15 20 同学乙 10 20 18 同学丙 20 16 10 同学丁 8 10 15 面试时间tij 这 4 名同学约定他们全部面试完以后一起离开公司，假定现在时间是早晨 8:00 问他们最早何时能离开公司? xik：第i个学生第k个面试阶段开始时间 目标函数 约束条件 每个学生k阶段结束时间早于k+1阶段开始 对每个k，后人始晚前人结束时间 k阶段开始晚于前段结束 sets: students/s1.