数学建模货物配送问题课程设计

1、安徽工业大学一数学建模论文货物运送问题组员：班级：指导教师：侯为根2013-7-301、问题重述一公司有二厂，分处A、B两市，另外还有4间具有存贮机构的库房，分别在P、QR和S市公司出售产品给6家客户G,C2,C6,由各库房或直接由工厂向客户供货。配送货物的费用由公司负担，单价见下表：表一I受货者i供货者；A市厂B市厂P库房Q库房R库房S库房P库房11:0.5Q库房!0.50.3R库房1.00.51S库房110.20.2客户gi-1.02.01.0客户C21.50.51.5客户Q11.50.50.52.00.2客户C42.01.51.01.5客户C5I:0.50.50.5客户G1.01.01.

2、51.5注：单位元/吨；划“-”表示无供货关系.某些客户表示喜欢由某厂或某库房供货.计有：Ci-A市厂C2P库房C5Q库房CR库房或S库房A市厂月供货量不能超过150千吨，B市厂月供货量不能超过200千吨。各库房的月最大流通量千吨数为表二库房PQRS流通量705010040各客户每月所必须满足的供货量为（单位：千吨）表三客户GGGC6要求货量501040356020现假设可以在T市和V市建新库房，和扩大Q市的库房，而库房的个数又不能多于4个，必要时可关闭P市和S市的库房。建新库房和扩建Q市库房的费用（计入利息）摊至每月为下表所列值（万元），它们的潜在的月流通量（千吨）也列于表中表四库房月费用流

3、通量T1.230V0.425Q（扩建）0.320关闭P市库房月省费用1万元；关闭S市库房月省0.5万元涉及新库房的配送费用单价（元/吨）见下表Ci1.2G0.60.40.5QC5-0.50.30.6C60.80.92、问题分析随着经济的发展、交通网络的不断健全以及各项科技的进步。使得各个行业竞争激烈，生产商要在满足客户要求与尽量减少生产成本之间面临更复杂决策。在整个配送问题中，所有的对象有三种，一种就是厂房，它是货物的产源地分别地处A、B两个市，它所生产的货物，可以直接运给客户，也可以放到库房里存放；第二种就是库房，用于存放来自于A、B两个厂房的生产物以及将货物配送给它的顾客，这种库房分别位于

4、P、Q、R、S市；第三种就是客户，接收由工厂或库房提供的货物；问题一、在配送过程中，我们需要建立一个数学模型来计算如何配货公司的运输费用最低，如何配送货物，既能满足客户的要求，又能为公司节约足够的资金。当然还要考虑到增加工厂和库房的生产能力对配送费用的影响，费用单价、客户对供应货物的最低要求以及工厂和库房生产能力各微小变化对配货方案的影响等因素来进行方案设计。设计出来的方案还要能体现出公司在什么样的改进下能获得更高的经济效益。可以用数学模型来建立最优解，进而解决设计方案的建立问题二、在问题一得基础上几乎没什么变化，A,B俩市供货量限制和客户需求量都没发生变化；改变的是库房，在T、V市新建库房，

5、扩建Q库房，即改变了流通量，必要时刻关闭P、S库房;也就是说到底对库房做出怎样的变化，这就引进了应否关闭P,S和应否新建T,V以及应否扩大Q库房，引进零、一变量解决好此问题公司与兴建新的库房，根据实际问题条件分析下应建那些新库房？Q市库房是否扩建？P市和S市库房应否关闭？配运费用最小的配货方案是什么？根据实际情况为公司减少运费提高利润，设计出合理的配货方案。3、符号说明问题一、A、B为生产厂，P、QRS为库房，C1、C2、C3C4CSC6为客户工厂向各库房和客户的供货量以及库房向客户的供货量如下两表（单位：千吨）工厂向各库房的供应量：、受货者供货PQRSAXiX12X13Xl4BX>iX

6、22X23X>4工厂和各库房向客户的供应量:爻货者供货/，ABPQRSC1X)nX012X013X)14X015X016C2X)21X022X023X)24X025X026C3X)31Xo32X033X)34X035X036C4X41X042X043X)44X045X046C5X)51Xo52X053X)54X055X056C6X)61X062X063X)64X065X066模型要求公司在配货时的最小运输费用，即:min问题二、AB给库房P、QRST、V的货物量为:Xi、X12、X13、Xl4>X15、Xl6Ki、X22、X23、X24、X25、X26；由A、B供给客户Ci、C2C

7、3C4C5C6的货物量为:yii、yi2>y13、yi4y15、yi6；y2i、y22、y23、y24、y25、y26；由库房P、QS、T、V供给客户Ci、C2C3C4C5C6的货物量为:Zii、Z2i、Z3i、Z4iZ5iZ6iZi2、Z22、Z32、Z42、Z52、Z62、Zi3、Z23、Z33、Z43、Z53、Z63、Zi4Z24Z34Z44Z53Z64Zi5、Z25、Z35、Z45、Z55、Z65、Zi6；Z26；Z36；Z46；Z56；Z66；由于最多只能用四个客房，故要确定选哪四个，即对P、QRS、T、V五个库房定一个零、五变量：ai、a2、a3、a，、a5表示库房P、QS、

8、T、V的0-i变量：ai为0表示关闭P库房，为i表示未关闭P库房；生为0表示未扩建Q库房，为i表示未关闭P库房；a3为0表示关闭R库房，为i表示未关闭R库房；a4为0表示未新建T库房，为i表示扩建T库房；a5为0表示新建V库房，为i表示扩建V库房；4、模型假设(D公司出售产品给6家客户C,C2,C6,由各库房或直接由工厂向客户供货。(2)A市厂月供货量不能超过i50千吨，B市厂月供货量不能超过200千吨。库房的月最大流通量保持不变，即在库房有货物剩余的情况下，月最大流通量不因此而加大。(3)某些客户表示喜欢由某厂或某库房供货.计有：CA市厂C2-P库房CQ库房C6-R库房或S库房假设顾客与库房

9、之间不存在喜好关系。(4)在问题的解决过程中，由于这个问题只提及运输费用的问题，而不考虑公司在货物卖出时的收益问题，所以我们只对运输上的经济情况进行讨论，不管运输时各个运输路线的单价如何变化，我们的模型都能将最好的方案给出来。(5)假设可以在T市和V市建新库房，和扩大Q市的库房，而库房的个数又不能多于4个，必要时可关闭P市和S市的库房。5、模型建立在配货过程中，可以由A市厂和B市厂直接向客户直接供货，也可以把两厂的货物运到P、QRS四个仓库之后再向客户供货，所以在这个模型中，我们首先把A,B看成生产地，同时又把它们作为与P、Q、R、S一样的库房来看待，并规定产地A、B不向库房AB运送货物，在处

10、理的时候，如果相互之间没有配送关系，我们可以认为配送货物的费用为“无穷大”，在具体运算时，我们再对“无穷大”赋予一个比较大的具体值。配送货物的费用由公司负担，单价见下表：受货者供货者jA市厂B市厂P库房Q库房R库房S库房P库房1Q库房！1R库房S库房0.510.50.311.00.50.20.2客户G1.02.01.0客尸Gi1.50.51.5客户C31.50.50.52.00.2客户C42.01.51.01.5客户C5_:0.50.50.5客户C61.01.01.51.5注：单位元/吨；划“-”表示无供货关系工厂向各库房的供应量：可以得出以下结果：位于A地的厂向P、QRS地库房供应货物所需运

11、费:位于B地的厂向P、QRS地库房供应货物所需运费:各库房货物的月流通量:库房PQRS流通量705010040则得之，Xl1<=70,Xl2+X2<=50,Xl3+X3<=100,Xl4+%4<=40；各客户每月所必须满足的供货量为（单位：千吨）客户GGG要求货量501040356020工厂和各库房向客户的供应量:供货者ABPQRSGX)nX012X013X)14X015X016QX)21X022X023X)24X025X026GX)31X032X033X)34X035X036C4X)41X042X043X)44X045X046C5X)51X052X053X)54X05

12、5X056X)61X062Xd63X)64X065X066问题一、可以得到以下线性方程：客户C1可分别由工厂A、B以及库房Q三方供应货物，且配货所需单价分别为1.0、2.0、1.0,则得为C1配货的花费为：C1客户的月需求为：X011X012X01450；客户C2可分别库房P、QR三方供应货物，且配货所需单价分别为1.5、0.5、1.5,则得为C配货的花费为：C2客户的月需求量为:X023X024X02510；客户C3可分别由工厂A以及库房P、QRS五方供应货物，且配货所需单价分别为1.5、0.5、0.5、2.0、0.2,则得为C3配货的花费为：C3客户的月需求量为：X031X033X034X

13、035X03640;客户C4可分别由工厂A以及库房PQS四方供应货物，且配货所需单价分别为2.0、1.5、1.0、1.5,则得为C4配货的花费为：C4客户的月需求量为：X041X043X044X04635；客户C5可分别由库房QR、S三方供应货物，且配货所需单价分别为0.5、0.5、0。5,则得为C5配货的花费为：C5客户的月需求量为：X054X055X05660；客户C6可分别由工厂A以及库房PRS四方供应货物，且配货所需单价分别为1.0、1.0、1.5、1.5,则得为C6配货的花费为：客户C6的月需求量为：X061X063X065X06620；根据上面式子中的所有变量以及题给意思可知:A厂

14、的每月最大供应量为：X11X12X13X14X011X031X041X061150；B厂的每月最大供应量为:X22X23X24X012200；从问题考虑，使费用最小，即为库房的存储量供完货后为0是最好的。P库房的货物量有A厂，运出的货物量为客户C2、C3C4C6,由约束条件如下：X023X033X043X063X110；Q库房的货物量有A、B厂，运出的货物量为客户C1、C2、C3C4C5,由约束条件如下：X014X024X034X044X054X12X220；R库房的货物量有A、B厂，运出的货物量为客户C2、C3C5C6,由约束条件如下：X025X035X055X065X13X230；S库房的

15、货物量有A、B厂，运出的货物量为客户C3>C4C5C6,由约束条件如下:X036X046X056X066X14X240；总结以上的模型，可得配货的最小运输费用问题实际上为一个线性规划模型：目标函数：min=0.5*Xii+0.5*Xi2+1.0Xi3+0.2*Xm+0.3*X22+0.5*X23+0.2*X24+Xii+2.0*X012+X14+I&X°23+0.5*X024+1.5*X025+1.5*X03i+0.5*X033+0.5*X034+2.0*X035+0.5*X036+2.0*X04l+1.5*X043+%44+1.5*X046+0.5*X054+0.5*

16、X055+0.5*X056+X)6l+%63+1.5*X065+1.5*X066；约束条件：Xii<=70;Xl2+X22<=50；X13+X23<=100;Xi4+Xm<=40;Xd1i+X)12+X)14=50;X023+X)24+X)25=10;Xo31+X)33+X)34+Xd35+X)36=40;Xo4l+X)43+X)44+Xp46=35;Xd54+X)55+X)56=60;Xo6l+X)63+X)65+Xo66=20;Xll+Xl2+Xl3+Xl4+Xoi1+X)3l+X)41+Xo6i<=15O;X?2+X23+X24+Xoi2<=2OO;X

17、o23+X)33+X)43+Xo63-X11=0；Xd14+X)24+X)34+Xd44+X)54-X12-X22=0；Xo25+X)35+X)55+Xo65-X13-X23=0；X336+X)46+X)56+Xo66-X14-X24=0；则建立模型如下:min=0.5*Xn+0.5*X12+1.OX13+O2X14+0.3*X22+0.5*X23+0.2*X24+X)n+2.0*X012+X)14+1.5*X023+0.5*X024+1.5*Xo25+1.5*X031+0.5*X033+0.5*X034+2.0*X035+0.5*X036+2.0*X041+1.5*X043+%44+1.5*

18、X046+0.5*X054+0.5*XO55+O.5*Xo56+X)61+Xo63+1.5*XO65+1.5*Xo66;X11<=70;X12+X22<=50;X13+X23<=100;X14+X24<=40;X011+X012+X014=50;X)23+X)24+X)25=10；X031+X)33+X)34+X)35+X)36=40；X)4l+X)43+X)44+X046=35；X+X+X'G。；居61+%63+%65+线66=2。；Xll+Xl2+Xl3+Xl4+X)11+X)3l+X41+X)61<=150X22+X3+X4+X)12<=200

19、；Xd23+X)33+X)43+X)63-X11=)；X014+X)24+X)34+X044+X)54-X12-X22=)；Xd25+X)35+X)55+X)65-X13X23=)；Xd36+X)46+X)56+X)66-X14X24=)；问题二建新库房和扩建Q市库房的费用（计入利息）摊至每月为下表所列值（万元），它们的潜在的月流通量（千吨）也列于表中库房月费用i流通量T1.2301V10.4251Q（扩建）0.320关闭P市库房月省费用1万元；关闭S市库房月省）.5万元G|0.60.4Gi0.5C40.5C50.30.6G|0.80.9I!问题二基于问题一的内容故可得出如下：ai、a2、a3

20、、a，、a5表示库房P、QS、T、V的0-1变量：ai为0表示关闭P库房，为1表示未关闭P库房；生为0表示未扩建Q库房，为1表示未关闭P库房；a3为0表示关闭R库房，为1表示未关闭R库房；a，为0表示未新建T库房，为1表示扩建T库房；a5为0表示新建V库房，为1表示扩建V库房；则对P库房的调整费用为-1+a1；则对Q房的调整费用为0.3*a2；则对R房的调整费用为-0.5+0.5*a3；则对T房的调整费用为1.2*a4；则对V房的调整费用为0.4*a5；由题意可知目标函数为货物配送费用与库房调整费用相加得出：Min=0.5*Xn+0.5*X12+X13+02Xm+0.6*X15+0.4*X16

21、+yn+1.5*y13+2*y14+y16+0.3*X22+0.5*X23+0.2*X24+0.4*X?5+0.3*X26+2*y21+1.5*Z12+0.5*Z13+1.5*Z14+Z16+Z21+0.5*Z22+0.5*Z23+Z24+0.5*Z25+1.5*Z32+2*Z33+0.5*Z35+1.5*Z36+0.2*Z43+1.5*Z44+0.5*Z45+1.5*Z46+1.2*z51+0.6*z52+0.5*z53+0.3*z55+0.8*z56+0.4*z62+0.5*Z64+0.6*z65+0.9*z66-1+a1+0.3*a2-0.5+0.5*a3+1.2*a4+0.4*a5约束

22、条件：有问题一可以轻松得出：X11+X12+X13+X14+X15+X16+y11+y13+y14+y16V=150；X22+X23+X4+X5+X6+y21V=200；最大存储量为70*a1千吨约束条件为：Xn<=70*a1；Q库房的存储量约束条件为：乂2+又2<=50+20*&R库房的存储量约束条件为：X13+X23V=100；S库房的存储量约束条件为：X14+X24V=40*a3；V房的存储量约束条件为：Xi6+X26<=25*a5；C1-C哈户需求量的约束条件：yii+y2+Z2i+Z5i=50;Zl2+Z22+Z32+Z52+Z62=10；yi3+Zi3+Z

23、23+Z33+Z43+Z53=40；yi4+Zi4+Z24+Z44+Z64=35；Z25+Z35+Z45+Z55+Z65=60；yi6+Zi6+Z36+Z46+Z56+Z66=20；由问题一可知同理有要求供给平衡Zi2+Zi3+Zi4+Zi6-Xii=0；Z2i+Z22+Z23+Z24-Xi2-X22=0；Z32+Z33+Z35+Z36-Xi3-X23=0；Z43+Z44+Z45+Z46-Xi4-X24=0；Z5i+Z52+Z53+Z55+Z56-Xi5-X25=0；Z62+Z64+Z65+Z66-Xi6-X26=0；库房的约束条件为：ai+a3+a4+a5<=2；综上所述，有以下数学模

24、型：Min=0.5*Xii+0.5*Xi2+Xi3+0.2*Xm+0.6*Xi5+0.4*Xi6+yii+i.5*yi3+2*yi4+yi6+0.3*X22+0.5*X23+0.2*X24+0.4*X25+0.3*X26+2*y2i+i.5*Z12+0.5*Z13+1.5*Z14+Zi6+Z2i+0.5*Z22+0.5*Z23+Z24+0.5*Z25+1.5*Z32+2*Z33+0.5*Z35+1.5*Z36+0.2*Z43+1.5*Z44+0.5*Z45+1.5*Z46+1.2*Z5i+0.6*Z52+0.5*Z53+0.3*Z55+0.8*z56+0.4*z62+0.5*Z64+0.6*z

25、65+0.9*z66-1+ai+0.3*a2-0.5+0.5*a3+1.2*a4+0.4*a5Xii+Xi2+Xi3+Xi4+Xi5+Xi6+yii+yi3+yi4+yi6<=i50；X22+X23+X24+X25+X26+y2i<=200；Xii<=70*ai；Xi2+X22<=50+20*a2;Xi3+X23<=100；Xi4+X24<=40*a3；Xi6+X26<=25*a5；yii+y2+Z2i+Z5i=50；Zi2+Z22+Z32+Z52+Z62=10；yi3+Zi3+Z23+Z33+Z43+Z53=40；yi4+Zl4+Z24+Z44+Z6

26、4=35；Z25+Z35+Z45+Z55+Z65=60；yi6+Zl6+Z36+Z46+Z56+Z66=20；Zl2+Zl3+Zl4+Zl6-Xl1=0；Z21+Z22+Z23+Z24-X12-X22=0；Z32+Z33+Z35+Z36-X13-X23=0；Z43+Z44+Z45+Z46-X14-X24=0；Z51+Z52+Z53+Z55+Z56-X15-X25=0；Z62+Z64+Z65+Z66-X16-X26=0；a1+a3+a4+a5<=2；6、模型求解：问题一：min=0.5*Xn+0.5*X12+1.0X13+0.2*Xm+0.3*X22+0.5*X23+0.2*X24+Xh+

27、2.0*X012+X14+15X023+05X024+1.5*X025+1.5*X031+0.5*X033+0.5*X034+2.0*X035+0.5*X036+2.0*X041+1.5*X043+X)44+1.5*X046+0.5*X054+0.5*X055+0.5*X056+X)61+Xd63+1.5*X065+1.5*X066；X11<=70；X12+X22<=50；X13+X23<=100；X14+X24<=40；X011+X012+X014=50；X023+X024+X025=10；X031+X033+X034+X035+X036=40；X041+X043+X

28、044+X046=35；X054+X055+X056=60；X061+X063+X065+X066=20；X11+X12+X13+X14+X011+X031+X041+X061<=150；X22+X23+X24+X012<=200；X323+X)33+X)43+Xd63-X11=0；X014+X)24+X)34+X044+X)54-X12-X22=0；X325+X)35+X)55+Xd65-X13-X23=0；X036+X)46+X)56+X066-X14-X24=0；end问题二:Min=0.5*Xii+0.5*X12+X13+02Xm+0.6*X15+0.4*Xi6+yii+1

29、.5*yi3+2*yi4+yie+0.3*X22+0.5*X23+0.2*X24+0.4*X25+0.3*X26+2*y2l+1.5*Z12+0.5*z13+1.5*z14+Zi6+Z2i+0.5*Z22+0.5*Z23+Z24+0.5*Z25+1.5*z32+2*Z33+0.5*Z35+1.5*Z36+0.2*Z43+1.5*Z44+0.5*Z45+1.5*Z46+1.2*z5i+0.6*Z52+0.5*z53+0.3*z55+0.8*z56+0.4*z62+0.5*z64+0.6*z65+0.9*z66-1+ai+0.3*a2-0.5+0.5*a3+1.2*a4+0.4*a5Xii+Xi2

30、+Xi3+Xi4+Xi5+Xi6+yii+yi3+yi4+y16V=150;X?2+X23+X24+X25+X26+y21V=200；Xiiv=70*ai；Xi2+X22<=50+20*a2;Xi3+X23<=100；Xi4+X24<=40*a3；Xi6+X26<=25*a5；yii+y2+Z2i+Z5i=50；Zi2+Z22+Z32+Z52+Z62=10；yi3+Zi3+Z23+Z33+Z43+Z53=40；yi4+Zi4+Z24+Z44+Z64=35；Z25+Z35+Z45+Z55+Z65=60；yi6+Zi6+Z36+Z46+Z56+Z66=20；Zi2+Zi3+

31、Zi4+Zi6-Xi1=0；Z21+Z22+Z23+Z24-X12-X22=0Z32+Z33+Z35+Z36-X13-X,23=0；Z43+Z44+Z45+Z46-X14-X24=0Z51+Z52+Z53+Z55+Z56-X15X25=0；Z62+Z64+Z65+Z66-X16-X26=0；ai+a3+&+a5V=2；bin(ai)；bin(a2)；bin(a3)；bin(a4)；bin(a5)；end7、结论解释这段本来是lingo运行的结果，现把结果给删去了，用表格体现的结果问题一：该问题的最小费用为：210,5000单位;对应的配货方案见下表:受货者'"供货者1

32、A市广B市广P库房Q库房R库房S库房P库房！40Q库房i1501R库房；!15-S库房140客户G50客户C2!10!客户C340客户C435客户C551540i客户C61:i201A-P:40千吨，A-S:4仃吨，A-Ci：5仃吨，A-G:20千吨;B-Q:50千吨,B-R:15千吨；Q-G:1仃吨,Q-C:35千吨,Q-G:5千吨；P-C3：4仃吨,R-C5:15千吨,S-Q:40千吨；问题二：该问题的最小费用为：106,4000单位;其具体的运输调运方案为：受货者i_佟货者_A币厂B币厂Q库市R库房S吊京V库存an*!anhmiar«nam>>>ankbana

33、ara；anQ库房j20R库房40S库房i:iV库房125:1客户G50客户G10客户C340i-客户C41025i客户C560客尸C620-关闭S仓库，新建V仓库；A-R:40千吨，运往G:50千吨，运往G:20千吨；B-Q:20千吨；运往V:25千吨Q-G:10千吨；运往C4:10千吨；运往G:60千吨；S-C3：40千吨;V-G:25千吨;在模型一中，我们发现虽然总费用最小，但是A、B两厂运送到s库房的货物数量并不确定，即配货方案并不唯一，经过分析，发现A、B运送到S的单价是相同的，所以不影响最后结果；在模型二中，我们也发现了类似的情况，因为从A厂运到C4和从A厂将货物运到P库房然后再运

34、到G的费用也是相同的，所以在这里的配送方案也不唯一。所以在实际中，可以结合实际情况进行调配8、模型优缺点模型的优点：(1) 模型将复杂的配货过程简单和直观化，利用简单的线性规划进行建模，再结合计算机进行运算，整个过程简单而易操作。(2) 模型具有一般性和普遍性，适合任何情况下的最优配货方案的设计。(3)该模型的方法还可以推广大更大的经济领域中去，具有很大的可塑性。模型的不足：(1)模型中的变量和参数太多，在输入的时候可能比较麻烦。(2)在模型的求解过程中，所给出的运货为零的方案在实际计算中并非真正为0,而是计算机在计算过程中采取了舍入的方法，这些对模型结果的影响并不大，因为数据足够小，而且在实

35、际中，如果运送的货物很少的话，对资源反而是一种浪费。9、模型推广1. 客户满意度：在模型二中，我们对客户喜好程度的处理是采取将近100%fe满足客户的要求，而客户对某个厂房和库房的喜好也是100%但是，在实际中并非如此。所以在处理客户的选择的时候应该引入客户对某厂房和某库房货物的满意度，即客户对该厂房或库房的喜好程度，所要求的货物在总的需求量中的比例等，来作为对货物配送中的一个主要的考虑因素。2. 公司的收益:在整个模型中，我们只考虑了怎么使运输费用最少的问题，却没有考虑货物卖出之后能挣多少的问题，一批货物运到这一客户和运到那一客户收益有什么不同，货物的多少对收益又有什么影响，在实际中，收益才

36、是配货的首要考虑问题。所以，在模型中加入这一因素，就可以实现对整个公司经济活动的宏观调控。10、附件(各种程序等)问题一中的lingo程序：min=0.5*Xii+0.5*Xi2+1.0Xi3+0.2*X14+0.3*X22+0.5*X23+0.2*X24+X)ii+2.0*Xoi2+X)i4+1.5*X023+0.5*X024+1.5*X025+1.5*X03i+0.5*X033+0.5*X034+2.0*X035+0.5*X036+2.0*X04i+1.5*X043+%44+1.5*X046+0.5*X054+0.5*X055+0.5*X056+X)6i+X063+1.5*X065+1.5*X066；Xii<=70;