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

数学建模货物配送问题课程设计数学建模货物配送问题课程设计数学建模货物配送问题课程设计资料仅供参考文件编号：2022年4月数学建模货物配送问题课程设计版本号：A修改号：1页次：1.0审核：批准:发布日期：安徽工业大学—数学建模论文货物运送问题组员:班级：指导教师：侯为根 2013-7-301、问题重述一公司有二厂，分处A、B两市，另外还有4间具有存贮机构的库房，分别在P、Q、R和S市。公司出售产品给6家客户C1,C2,…,C6，由各库房或直接由工厂向客户供货。配送货物的费用由公司负担，单价见下表：表一受货者供货者A市厂B市厂P库房Q库房R库房S库房P库房----Q库房R库房S库房客户C1------------客户C2------------客户C3----客户C4--------客户C5------------客户C6--------受货者表1：单位运输费用（千元/吨）A基地B基地P库房Q库房R库房S库房P库房Q库房1——R库房S库房客户甲12——1——客户乙——客户丙——1————客户丁——注:单位元/吨;划“----”表示无供货关系.某些客户表示喜欢由某厂或某库房供货.计有:C1--------A市厂C2--------P库房C5--------Q库房C6--------R库房或S库房A市厂月供货量不能超过150千吨，B市厂月供货量不能超过200千吨。各库房的月最大流通量千吨数为表二库房PQRS流通量705010040表5：库房容量（吨）库房PQRS流量70605050各客户每月所必须满足的供货量为（单位：千吨）表三客户C1C2C3C4C5C6要求货量501040356020表3：客户需求关系（吨）客户甲乙丙丁需求货量50304030现假设可以在T市和V市建新库房，和扩大Q市的库房，而库房的个数又不能多于4个，必要时可关闭P市和S市的库房。建新库房和扩建Q市库房的费用（计入利息）摊至每月为下表所列值（万元），它们的潜在的月流通量（千吨）也列于表中表四库房月费用流通量TVQ（扩建）302520关闭P市库房月省费用1万元；关闭S市库房月省万元。涉及新库房的配送费用单价（元/吨）见下表表五供货受货ABTVTVC1----C2C3----C4----C5C62、问题分析随着经济的发展、交通网络的不断健全以及各项科技的进步。使得各个行业竞争激烈，生产商要在满足客户要求与尽量减少生产成本之间面临更复杂决策。在整个配送问题中，所有的对象有三种，一种就是厂房，它是货物的产源地分别地处Ａ、Ｂ两个市，它所生产的货物，可以直接运给客户，也可以放到库房里存放；第二种就是库房，用于存放来自于Ａ、Ｂ两个厂房的生产物以及将货物配送给它的顾客，这种库房分别位于Ｐ、Ｑ、Ｒ、Ｓ市；第三种就是客户，接收由工厂或库房提供的货物；问题一、在配送过程中，我们需要建立一个数学模型来计算如何配货公司的运输费用最低，如何配送货物，既能满足客户的要求，又能为公司节约足够的资金。当然还要考虑到增加工厂和库房的生产能力对配送费用的影响，费用单价、客户对供应货物的最低要求以及工厂和库房生产能力各微小变化对配货方案的影响等因素来进行方案设计。设计出来的方案还要能体现出公司在什么样的改进下能获得更高的经济效益。可以用数学模型来建立最优解，进而解决设计方案的建立。问题二、在问题一得基础上几乎没什么变化，A，B俩市供货量限制和客户需求量都没发生变化；改变的是库房，在T、V市新建库房，扩建Q库房，即改变了流通量，必要时刻关闭P、S库房；也就是说到底对库房做出怎样的变化，这就引进了应否关闭P，S和应否新建T，V以及应否扩大Q库房，引进零、一变量解决好此问题公司与兴建新的库房，根据实际问题条件分析下应建那些新库房Q市库房是否扩建P市和S市库房应否关闭配运费用最小的配货方案是什么根据实际情况为公司减少运费提高利润，设计出合理的配货方案。3、符号说明问题一、A、B为生产厂，P、Q、R、S为库房，C1、C2、C3、C4、C5、C6为客户。工厂向各库房和客户的供货量以及库房向客户的供货量如下两表(单位:千吨)

工厂向各库房的供应量：受货者供货者PＱＲＳＡX11X12X13X14ＢX21X22X23X24工厂和各库房向客户的供应量：受货者供货者ABPQRSC1X011X012X013X014X015X016C2X021X022X023X024X025X026C3X031X032X033X034X035X036C4X041X042X043X044X045X046C5X051X052X053X054X055X056C6X061X062X063X064X065X066模型要求公司在配货时的最小运输费用,即：min问题二、A、B给库房P、Q、R、S、T、V的货物量为：X11、X12、X13、X14、X15、X16；X21、X22、X23、X24、X25、X26；由A、B供给客户C1、C2、C3、C4、C5、C6的货物量为：y11、y12、y13、y14、y15、y16;y21、y22、y23、y24、y25、y26;由库房P、Q、R、S、T、V供给客户C1、C2、C3、C4、C5、C6的货物量为：z11、z12、z13、z14、z15、z16;z21、z22、z23、z24、z25、z26;z31、z32、z33、z34、z35、z36;z41、z42、z43、z44、z45、z46;z51、z52、z53、z53、z55、z56;z61、z62、z63、z64、z65、z66;由于最多只能用四个客房，故要确定选哪四个，即对P、Q、R、S、T、V五个库房定一个零、五变量：a1、a2、a3、a4、a5表示库房P、Q、S、T、V的0-1变量：a1为0表示关闭P库房，为1表示未关闭P库房；a2为0表示未扩建Q库房，为1表示未关闭P库房；a3为0表示关闭R库房，为1表示未关闭R库房；a4为0表示未新建T库房，为1表示扩建T库房；a5为0表示新建V库房，为1表示扩建V库房；4、模型假设（1）公司出售产品给6家客户C1,C2,…,C6，由各库房或直接由工厂向客户供货。（2）A市厂月供货量不能超过150千吨，B市厂月供货量不能超过200千吨。库房的月最大流通量保持不变，即在库房有货物剩余的情况下，月最大流通量不因此而加大。(3)某些客户表示喜欢由某厂或某库房供货.计有:C1--------A市厂C2--------P库房C5--------Q库房C6--------R库房或S库房假设顾客与库房之间不存在喜好关系。（4）在问题的解决过程中，由于这个问题只提及运输费用的问题，而不考虑公司在货物卖出时的收益问题，所以我们只对运输上的经济情况进行讨论，不管运输时各个运输路线的单价如何变化，我们的模型都能将最好的方案给出来。（5）假设可以在T市和V市建新库房，和扩大Q市的库房，而库房的个数又不能多于4个，必要时可关闭P市和S市的库房。5、模型建立在配货过程中,可以由A市厂和B市厂直接向客户直接供货,也可以把两厂的货物运到P、Q、R、S四个仓库之后再向客户供货，所以在这个模型中，我们首先把A，B看成生产地，同时又把它们作为与P、Q、R、S一样的库房来看待，并规定产地A、B不向库房A、B运送货物，在处理的时候，如果相互之间没有配送关系，我们可以认为配送货物的费用为“无穷大”，在具体运算时，我们再对“无穷大”赋予一个比较大的具体值。配送货物的费用由公司负担，单价见下表：受货者供货者A市厂B市厂P库房Q库房R库房S库房P库房----Q库房R库房S库房客户C1------------客户C2------------客户C3----客户C4--------客户C5------------客户C6--------注:单位元/吨;划“----”表示无供货关系.工厂向各库房的供应量：受货者供货者pQRSＡX11X12X13X14ＢX21X22X23X24可以得出以下结果：位于A地的厂向P、Q、R、S地库房供应货物所需运费：位于B地的厂向P、Q、R、S地库房供应货物所需运费：各库房货物的月流通量：库房PQRS流通量705010040则得之，X11<=70,X12+X22<=50,X13+X23<=100,X14+X24<=40;各客户每月所必须满足的供货量为（单位：千吨）客户C1C2C3C4C5C6要求货量501040356020工厂和各库房向客户的供应量：受货者供货者ABPQRSC1X011X012X013X014X015X016C2X021X022X023X024X025X026C3X031X032X033X034X035X036C4X041X042X043X044X045X046C5X051X052X053X054X055X056C6X061X062X063X064X065X066问题一、可以得到以下线性方程：客户C1可分别由工厂A、B以及库房Q三方供应货物，且配货所需单价分别为、、，则得为C1配货的花费为：C1客户的月需求为：;客户C2可分别库房P、Q、R三方供应货物，且配货所需单价分别为、、，则得为C配货的花费为：C2客户的月需求量为：;客户C3可分别由工厂A以及库房P、Q、R、S五方供应货物，且配货所需单价分别为、、、、，则得为C3配货的花费为：C3客户的月需求量为：；客户C4可分别由工厂A以及库房P、Q、S四方供应货物，且配货所需单价分别为、、、，则得为C4配货的花费为：C4客户的月需求量为：；客户C5可分别由库房Q、R、S三方供应货物，且配货所需单价分别为、、0。5，则得为C5配货的花费为：C5客户的月需求量为：；客户C6可分别由工厂A以及库房P、R、S四方供应货物，且配货所需单价分别为、、、，则得为C6配货的花费为：客户C6的月需求量为：;根据上面式子中的所有变量以及题给意思可知：A厂的每月最大供应量为：;B厂的每月最大供应量为：;从问题考虑，使费用最小，即为库房的存储量供完货后为0是最好的。P库房的货物量有A厂，运出的货物量为客户C2、C3、C4、C6，由约束条件如下：；Q库房的货物量有A、B厂，运出的货物量为客户C1、C2、C3、C4、C5，由约束条件如下：；R库房的货物量有A、B厂，运出的货物量为客户C2、C3、C5、C6，由约束条件如下：；S库房的货物量有A、B厂，运出的货物量为客户C3、C4、C5、C6，由约束条件如下：；总结以上的模型,可得配货的最小运输费用问题实际上为一个线性规划模型:目标函数：min=*X11+*X12++*X14+*X22+*X23+*X24+X011+*X012+X014+*X023+*X024+*X025+*X031+*X033+*X034+*X035+*X036+*X041+*X043+X044+*X046+*X054+*X055+*X056+X061+X063+*X065+*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+X044+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;X023+X033+X043+X063-X11=0；X014+X024+X034+X044+X054-X12-X22=0；X025+X035+X055+X065-X13-X23=0；X036+X046+X056+X066-X14-X24=0；则建立模型如下：min=*X11+*X12++*X14+*X22+*X23+*X24+X011+*X012+X014+*X023+*X024+*X025+*X031+*X033+*X034+*X035+*X036+*X041+*X043+X044+*X046+*X054+*X055+*X056+X061+X063+*X065+*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+X044+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;X023+X033+X043+X063-X11=0；X014+X024+X034+X044+X054-X12-X22=0；X025+X035+X055+X065-X13-X23=0；X036+X046+X056+X066-X14-X24=0；问题二建新库房和扩建Q市库房的费用（计入利息）摊至每月为下表所列值（万元），它们的潜在的月流通量（千吨）也列于表中库房月费用流通量TVQ（扩建）302520关闭P市库房月省费用1万元；关闭S市库房月省万元。涉及新库房的配送费用单价（元/吨）见下表供货受货ABTVTVC1----C2C3----C4----C5C6问题二基于问题一的内容故可得出如下：a1、a2、a3、a4、a5表示库房P、Q、S、T、V的0-1变量：a1为0表示关闭P库房，为1表示未关闭P库房；a2为0表示未扩建Q库房，为1表示未关闭P库房；a3为0表示关闭R库房，为1表示未关闭R库房；a4为0表示未新建T库房，为1表示扩建T库房；a5为0表示新建V库房，为1表示扩建V库房；则对P库房的调整费用为-1+a1;则对Q房的调整费用为*a2;则对R房的调整费用为+*a3;则对T房的调整费用为*a4;则对V房的调整费用为*a5;由题意可知目标函数为货物配送费用与库房调整费用相加得出：Min=*X11+*X12+X13+*X14+*X15+*X16+y11+*y13+2*y14+y16+*X22+*X23+*X24+*X25+*X26+2*y21+*z12+*z13+*z14+z16+z21+*z22+*z23+z24+*z25+*z32+2*z33+*z35+*z36+*z43+*z44+*z45+*z46+*z51+*z52+*z53+*z55+*z56+*z62+*z64+*z65+*z66-1+a1+*+*a3+*a4+*a5约束条件：有问题一可以轻松得出：X11+X12+X13+X14+X15+X16+y11+y13+y14+y16<=150;X22+X23+X24+X25+X26+y21<=200;最大存储量为70*a1千吨约束条件为：X11<=70*a1;Q库房的存储量约束条件为：X12+X22<=50+20*a2；R库房的存储量约束条件为：X13+X23<=100;S库房的存储量约束条件为：X14+X24<=40*a3;V房的存储量约束条件为：X16+X26<=25*a5;C1-C6客户需求量的约束条件：y11+y2+z21+z51=50;z12+z22+z32+z52+z62=10;y13+z13+z23+z33+z43+z53=40;y14+z14+z24+z44+z64=35;z25+z35+z45+z55+z65=60;y16+z16+z36+z46+z56+z66=20;由问题一可知同理有要求供给平衡z12+z13+z14+z16-X11=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;综上所述，有以下数学模型：Min=*X11+*X12+X13+*X14+*X15+*X16+y11+*y13+2*y14+y16+*X22+*X23+*X24+*X25+*X26+2*y21+*z12+*z13+*z14+z16+z21+*z22+*z23+z24+*z25+*z32+2*z33+*z35+*z36+*z43+*z44+*z45+*z46+*z51+*z52+*z53+*z55+*z56+*z62+*z64+*z65+*z66-1+a1+*+*a3+*a4+*a5X11+X12+X13+X14+X15+X16+y11+y13+y14+y16<=150;X22+X23+X24+X25+X26+y21<=200;X11<=70*a1;X12+X22<=50+20*a2；X13+X23<=100;X14+X24<=40*a3;X16+X26<=25*a5;y11+y2+z21+z51=50;z12+z22+z32+z52+z62=10;y13+z13+z23+z33+z43+z53=40;y14+z14+z24+z44+z64=35;z25+z35+z45+z55+z65=60;y16+z16+z36+z46+z56+z66=20;z12+z13+z14+z16-X11=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=*X11+*X12++*X14+*X22+*X23+*X24+X011+*X012+X014+*X023+*X024+*X025+*X031+*X033+*X034+*X035+*X036+*X041+*X043+X044+*X046+*X054+*X055+*X056+X061+X063+*X065+*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+X044+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;X023+X033+X043+X063-X11=0；X014+X024+X034+X044+X054-X12-X22=0；X025+X035+X055+X065-X13-X23=0；X036+X046+X056+X066-X14-X24=0；end问题二：Min=*X11+*X12+X13+*X14+*X15+*X16+y11+*y13+2*y14+y16+*X22+*X23+*X24+*X25+*X26+2*y21+*z12+*z13+*z14+z16+z21+*z22+*z23+z24+*z25+*z32+2*z33+*z35+*z36+*z43+*z44+*z45+*z46+*z51+*z52+*z53+*z55+*z56+*z62+*z64+*z65+*z66-1+a1+*+*a3+*a4+*a5X11+X12+X13+X14+X15+X16+y11+y13+y14+y16<=150;X22+X23+X24+X25+X26+y21<=200;X11<=70*a1;X12+X22<=50+20*a2；X13+X23<=100;X14+X24<=40*a3;X16+X26<=25*a5;y11+y2+z21+z51=50;z12+z22+z32+z52+z62=10;y13+z13+z23+z33+z43+z53=40;y14+z14+z24+z44+z64=35;z25+z35+z45+z55+z65=60;y16+z16+z36+z46+z56+z66=20;z12+z13+z14+z16-X11=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;@bin(a1);@bin(a2);@bin(a3);@bin(a4);@bin(a5);end7、结论解释{这段本来是lingo运行的结果，现把结果给删去了，用表格体现的结果}问题一：该问题的最小费用为：单位；对应的配货方案见下表：受货者供货者A市厂B市厂P库房Q库房R库房S库房P库房40Q库房50R库房15S库房40客户C150客户C210客户C340客户C435客户C551540客户C620A-P：40千吨，A-S：40千吨，A-C1：50千吨，A-C6：20千吨；B-Q：50千吨，B-R：15千吨；Q-C2：10千吨，Q-C4：35千吨，Q-C5：5千吨；P-C3：40千吨，R-C5：15千吨，S-C5：40千吨；问题二：该问题的最小费用为：单位；其具体的运输调运方案为：受货者供货者A市厂B市厂Q库房R库房S库房V库房Q库房20R库房40S库房V库房25客户C150客户C210客户C340客户C41025客户C560客户C620关闭S仓库，新建V仓库；R：40千吨，运往C1：50千吨，运往C6:20千吨；Q：20千吨；运往V：25千吨Q-C2：10千吨；运往C4：10千吨；运往C5:60千吨；S-C3：40千吨；V-C4：25千吨；在模型一中，我们发现虽然总费用最小，但是A、B两厂运送到s库房的货物数量并不确定，即配货方案并不唯一，经过分析，发现A、B运送到S的单价是相同的，所以不影响最后结果；在模型二中，我们也发现了类似的情况，因为从A厂运到C4和从A厂将货物运到P库房然后再运到C4的费用也是相同的，所以在这里的配送方案也不唯一。所以在实际中，可以结合实际情况进行调配模型优缺点模型的优点：

（1）

模型将复杂的配货过程简单和直观化，利用简单的线性规划进行建模，再结合计算机进行运算，整个过程简单而易操作。

（2）

模型具有一般性和普遍性，适合任何情况下的最优配货方案的设计。（3）

该模型的方法还可以推广大更大的经济领域中去，具有很大的可塑性。模型的不足：（1）模型中的变量和参数太多，在输入的时候可能比较麻烦。

（2）在模型的求解过程中，所给出的运货为零的方案在实际计算中并非真正为0，而是计算机在计算过程中采取了舍入的方法，这些对模型结果的影响并不大，因为数据足够小，而且在实际中，如果运送的货物很少的话，对资源反而是一种浪费。9、模型推广1．

客户满意度：

在模型二中，我们对客户喜好程度的处理是采取将近100%地满足客户的要求，而客户对某个厂房和库房的喜好也是100%，但是，在实际中并非如此。所以在处理客户的选择的时候应该引入客户对某厂房和某库房货物的满意度，即客户对该厂房或库房的喜好程度，所要求的货物在总的需求量中的比例等，来作为对货物配送中的一个主要的考虑因素。

2．

公司的收益：在整个模型中，我们只考虑了怎么使运输费用最少的问题，却没有考虑货物卖出之后能挣多少的问题，一批货物运到这一客户和运到那一客户收益有什么不同，货物的多少对收益又有什么影响，在实际中，收益才是配货的首要考虑问题。所以，在模型中加入这一因素，就可以实现对整个公司经济活动的宏观调控。10、附件（各种程序等）问题一中的lingo程序：min=*X11+*X12++*X14+*X22+*X23+*X24+X011+*X012+X014+*X023+*X024+*X025+*X031+*X033+*X034+*X035+*X036+*X041+*X043+X044+*X