应用LINGO软件求解运输问题

20122013学年第 一 学期合肥学院数理系实验报告课程名称： 运筹学 实验项目： 应用LINGO软件求解运输问题 实验类别：综合性 设计性 验证性 专业班级： 姓 名： 学 号： 实验地点： 实验时间： 指导教师： 成 绩： 一.实验目的1、学会使用LINGO软件求解运输问题的步骤与方法。2、掌握使用LINGO对运输问题的求解功能，并对结果进行分析。二.实验内容1.已知某企业有甲、乙、丙三个分厂生产一种产品，其产量分别为7、9、7个单位，需运往A、B、C、D四个门市部，各门市部需要量分别为3、5、7、8个单位。已知单位运价如下表。试确定运输计划使总运费最少。ABCD甲12131011乙10121410丙141115122.现在要在五个工人中确定四个人来分别完成四项工作中的一项工作。由于每个工人的技术特长不同，他们完成各项工作所需的工时也不同。每个工人完成各项工作所需工时如下表所示，试找出一个工作分配方案，使总工时最小。工作工人ABCD943746565475752310674三. 模型建立1.由题设知，总产量为：7+9+7=23个单位，总销量为：3+5+7+8=23个单位，所以这是一个产销平衡的运输问题。设代表从第个产地运往第个销地的数量，为总运费。表示第个产地的产量，表示第个销地的销量表示从第个产地运往第个销地的单位产品运输费用。则该问题的数学模型为：2. 设0-1变量，则该问题的数学模型为：四. 模型求解（含经调试后正确的源程序）1、编写程序1-1.m如下：model:sets:warehouses/wh1.wh3/: capacity;vendors/v1.v4/: demand;links(warehouses,vendors): cost, volume;endsetsdata: capacity=7 9 7； demand=3 5 7 8； cost= 12 13 10 11 10 12 14 10 14 11 15 12；enddatamin=sum(links(I,J): cost(I,J)*volume(I,J);for(vendors(J): sum(warehouses(I): volume(I,J)=demand(J);for(warehouses(I): sum(vendors(J): volume(I,J)=capacity(I); end 2、编写程序2-1.m如下：model:sets: workers/w1.w5/; jobs/j1.j4/; links(workers,jobs): cost,volume;Endsetsdata: cost=9 4 3 7 4 6 5 6 5 4 7 5 7 5 2 3 10 6 7 4;enddatamin=sum(links: cost*volume); for(workers(I): sum(jobs(J): volume(I,J)=1); for(jobs(J): sum(workers(I): volume(I,J)=1); for(links(i,j): bin(volume(i,j);End五结果分析1、运行结果： Global optimal solution found. Objective value: 239.0000 Infeasibilities: 0. Total solver iterations: 6 Variable Value Reduced Cost CAPACITY( WH1) 7. 0. CAPACITY( WH2) 9. 0. CAPACITY( WH3) 7. 0. DEMAND( V1) 3. 0. DEMAND( V2) 5. 0. DEMAND( V3) 7. 0. DEMAND( V4) 8. 0. COST( WH1, V1) 12.00000 0. COST( WH1, V2) 13.00000 0. COST( WH1, V3) 10.00000 0. COST( WH1, V4) 11.00000 0. COST( WH2, V1) 10.00000 0. COST( WH2, V2) 12.00000 0. COST( WH2, V3) 14.00000 0. COST( WH2, V4) 10.00000 0. COST( WH3, V1) 14.00000 0. COST( WH3, V2) 11.00000 0. COST( WH3, V3) 15.00000 0. COST( WH3, V4) 12.00000 0. VOLUME( WH1, V1) 0. 1. VOLUME( WH1, V2) 0. 3. VOLUME( WH1, V3) 7. 0. VOLUME( WH1, V4) 0. 0. VOLUME( WH2, V1) 3. 0. VOLUME( WH2, V2) 0. 3. VOLUME( WH2, V3) 0. 5. VOLUME( WH2, V4) 6. 0. VOLUME( WH3, V1) 0. 2. VOLUME( WH3, V2) 5. 0. VOLUME( WH3, V3) 0. 4. VOLUME( WH3, V4) 2. 0. Row Slack or Surplus Dual Price 1 239.0000 -1. 2 0. -12.00000 3 0. -11.00000 4 0. -11.00000 5 0. -12.00000 6 0. 1. 7 0. 2. 8 0. 0.所以，最优调运方案为：甲C：7单位；甲D：0单位；乙A：3单位；乙D：6单位；丙B：5单位；丙D：2单位。最少总运费为：239。2、运行结果：Global optimal solution found. Objective value: 14.00000 Objective bound: 14.00000 Infeasibilities: 0. Extended solver steps: 0 Total solver iterations: 0 Variable Value Reduced Cost COST( W1, J1) 9. 0. COST( W1, J2) 4. 0. COST( W1, J3) 3. 0. COST( W1, J4) 7. 0. COST( W2, J1) 4. 0. COST( W2, J2) 6. 0. COST( W2, J3) 5. 0. COST( W2, J4) 6. 0. COST( W3, J1) 5. 0. COST( W3, J2) 4. 0. COST( W3, J3) 7. 0. COST( W3, J4) 5. 0. COST( W4, J1) 7. 0. COST( W4, J2) 5. 0. COST( W4, J3) 2. 0. COST( W4, J4) 3. 0. COST( W5, J1) 10.00000 0. COST( W5, J2) 6. 0. COST( W5, J3) 7. 0. COST( W5, J4) 4. 0. VOLUME( W1, J1) 0. 9. VOLUME( W1, J2) 0. 4. VOLUME( W1, J3) 1. 3. VOLUME( W1, J4) 0. 7. VOLUME( W2, J1) 1. 4. VOLUME( W2, J2) 0. 6. VOLUME( W2, J3) 0. 5. VOLUME( W2, J4) 0. 6. VOLUME( W3, J1) 0. 5. VOLUME( W3, J2) 1. 4. VOLUME( W3, J3) 0. 7. VOLUME( W3, J4) 0. 5. VOLUME( W4, J1) 0. 7. VOLUME( W4, J2) 0. 5. VOLUME( W4, J3) 0. 2. VOLUME( W4, J4) 1. 3. VOLUME( W5, J1) 0. 10.00000 VOLUME( W5, J2) 0. 6. VOLUME( W5, J3) 0. 7. VOLUME( W5, J4) 0. 4. Row Slack or Surplus Dual Price 1 14.00000 -1. 2 0. 0. 3 0. 0. 4 0. 0. 5 0. 0. 6 1. 0. 7 0. 0. 8 0. 0. 9