《运筹学课程设计报告书》由会员分享，可在线阅读，更多相关《运筹学课程设计报告书（7页珍藏版）》请在金锄头文库上搜索。
1、运筹学课程设计报告书专业班级：_姓 名：指导教师：日 期：一课程设计的目的和意义运筹学是一门多学科的定量优化技术，为了从理论与实践的结合上，提高学 生应用运筹学方法与计算机软件的独立工作能力，本着“突出建模，结合软件， 加强应用”的指导思想，以学生自己动手为主，对一些实际题目进行构模，再运 用计算机软件进行求解，对解进行检验和评价，写出课程设计报告。二课程设计的时间本课程设计时间 1 周。三课程设计的基本任务和要求由于不同的同学选择的方向不同，因此给出如下两种要求，完成其一即可： 1 选择建模的同学：利用运筹学基本知识对所选案例建立合适的数学模 型，然后利用 winQSB、LINDO、LING
2、O 或者其它数学软件进行求解；2 选择编程的同学：根据运筹学基本原理以及所掌握的计算机语言知识， 对于运筹学中部分算法编写高级语言的具有可用性的程序软件。四课程设计的问题叙述临海市华安机械厂的潘厂长正考虑将该厂的一部分在市区的生产车间搬该 市的卫星城镇，好处是土地、房租费及排污处理费用都较便宜，但这样做会增加 车间之间的交通运输费用。该厂原在市区车间有 A、B、C、D、E 五个，计划搬迁去的卫星城镇有甲、乙 两处。规定无论留在市区或甲、乙两卫星城镇均不得多于 3 个车间。从市区搬至卫星城带来的年费用节约见表 4-24所示：表 4-24单位：万元 / 年ABCDE搬至甲1001501002005
3、0搬至乙100200150150150但搬迁后带来运输费用增加由C和D值决定，C为i和k车间之间的年运ikjlik量， D 为市区同卫星城镇间单位运量的运费，具体数据分别见表4-25和表4-26. jl表 4-25 C 值单位： t/ 年ikBCDEA0100015000B140012000C02000D700表4-26Djl值单位：元/t甲乙市区甲50140130乙5090市区50请为潘厂长提供一个决策建议方案，哪几个车间搬至卫星城镇及搬至甲还是 乙，能带来最大的经济上的好处。五模型的假设和建立设x为bool型变量，当j车间在i地时，此值为1,否则，此值为0。其中ij表示车间在甲地为x，表示
4、车间在乙地为x，表示车间在市区为x，A、B、C、1j2 j3jD、E车间在i地则用x、x、x、x、x表示。i1 i2 i3 i 4 i5设f为j车间搬至i地每年节约的费用，单位为万元每年ij此题所求为给厂带来最大的经济上的利润，而利润等于搬车间后带来的总节 约费用减去增加的总运输费用，即目标函数可以表示为：z = 工f * x - 工（D * 工为 （x * x * C ）；ij ijiji,i1 j, j1 i1,j1i=1 j=1i=1 j=1i1=1 j1=1其中，工f *x表示节约的总费用，工（x *x *C ）表示i地到j地ij iji,i1 j,j1i1,j1i=1 j=1i1=1
5、 j1=1之间发生的所有运输总量，工（D *工（x *x *C ）即表示增加的总运iji,i1 j, j1 i1,j1i=1 j=1i1=1 j1=1输费用。规定无论留在市区或甲、乙两卫星城镇均不得多于3个车间，即：工x 3，ijj=1i = 1,2,3 。另外，一共只有五个车间且每个车间只有一个，所以有：工x = 5，以及 ij i=1 j=1工 x = 1， i = 1,2,3,4,5。ji六模型求解此题用LINGO求解。程序如下：model:sets:dizhi/1.3/;chejian/1.5/; xuanze(dizhi,chejian):x,f; feiyong(chejian,c
6、hejian):C; yunfei(dizhi,dizhi):D; endsetsdata:C=0,0,1000,1500,0,0,0,1400,1200,01000,1400,0,0,2000,1500,1200,0,0,700, 0,0,2000,700,0;f=100,150,100,200,50,100,200,150,150,150, 0,0,0,0,0;D=0.005,0.014,0.0130,0.005,0.009,0,0,0.005;enddata!三个可供选择的地址;!五个车间;!未知变量,节约的费用!增加的费用;!单位运费;!增加的费用;!节约的费用;!单位运费;max=s
7、um(dizhi(i):sum(chejian(j):f(i,j)*x(i,j) !总收益等于总节约费 用减去增加的总费用;-sum(dizhi(i):!增加的总费用等于各地增加的费用之和;sum(dizhi(j):!各地增加的费用等于本地到其他地方增加的费用之和;D(i,j)*sum(chejian(i1):!本地到其他地方增加的费用等于单位运费乘以总运量;sum(chejian(j1): !总运量等于本地所有车间到他 地所有车间的运量之和;x(i,i1)*x(j,j1)*C(i1,j1);! 只有车间确实存在的时候，运量才发 生;for(dizhi(i):sum(chejian(j):x(
8、i,j)=3);!每个地方的车间总数不能超过3个;sum(dizhi(i):sum(chejian(j):x(i,j)=5;! 总共5个车间;for(chejian(i):sum(dizhi(j):x(j,i)=1);!车间有且只有一个;for(dizhi(i):for(chejian(j):bin(x(i,j);以为bool型变量; end结果如下：Local optimal solution found.710.400004Objective value:Extended solver steps:Total solver iterations:VariableValueReduced C
9、ostX( 1, 1)1.000000-52.80000X( 1, 2)0.0000000.000000X( 1, 3)0.0000000.000000X( 1, 4)1.000000-99.20000X( 1, 5)0.0000000.000000X( 2, 1)0.0000000.000000X( 2, 2)1.0000000.000000X( 2, 3)1.000000-8.800000X( 2, 4)0.0000000.000000X( 2, 5)1.000000-54.40000X( 3, 1)0.00000046.70000X( 3, 2)0.000000146.6000X( 3,
10、 3)0.00000086.00000X( 3, 4)0.00000095.80000X( 3, 5)0.00000042.10000F( 1, 1)100.00000.000000F( 1, 2)150.00000.000000F( 1, 3)100.00000.000000F( 1, 4)200.00000.000000F( 1, 5)50.000000.000000F( 2, 1)100.00000.000000F( 2, 2)200.00000.000000F( 2, 3)150.00000.000000F( 2, 4)150.00000.000000F( 2, 5)150.00000
11、.000000F( 3, 1)0.0000000.000000F( 3, 2)0.0000000.000000F( 3, 3)0.0000000.000000F( 3, 4)0.0000000.000000F( 3, 5)0.0000000.000000C( 1, 1)0.0000000.000000C(1,2)0.0000000.000000C(1,3)1000.0000.000000C(1,4)1500.0000.000000C(1,5)0.0000000.000000C(2,1)0.0000000.000000C(2,2)0.0000000.000000C(2,3)1400.0000.0
12、00000C(2,4)1200.0000.000000C(2,5)0.0000000.000000C(3,1)1000.0000.000000C(3,2)1400.0000.000000C(3,3)0.0000000.000000C(3,4)0.0000000.000000C(3,5)2000.0000.000000C(4,1)1500.0000.000000C(4,2)1200.0000.000000C(4,3)0.0000000.000000C(4,4)0.0000000.000000C(4,5)700.00000.000000C(5,1)0.0000000.000000C(5,2)0.0
13、000000.000000C(5,3)2000.0000.000000C(5,4)700.00000.000000C(5,5)0.0000000.000000D(1,1)0.5000000E-020.000000D(1,2)0.1400000E-010.000000D(1,3)0.1300000E-010.000000D(2,1)0.0000000.000000D(2,2)0.5000000E-020.000000D(2,3)0.9000000E-020.000000D(3,1)0.0000000.000000D(3,2)0.0000000.000000D(3,3)0.5000000E-020.000000RowSlack or SurplusDual Price1710.40001.00000021.0000000.00000030.00000050.8000043.0000000.00000050.000000
