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

安徽工业大学
—数学建模论文
货
物
运
送
问
题
组 员: 班 级： 指导教师：侯为根
2013-7-30
1、问题重述
一公司有二厂，分处 A、B 两市，另外还有 4 间具有存贮机构的库房，分别在 P、Q、R 和 S 市。 公司出售产品给 6 家客户 C1,C2,…,C6，由各库房或直接由工厂向客户供货。
配送货物的费用由公司负担，单价见下表：
表一
受货者
A 市厂
B 市厂
供货者 P 库房 Q 库房
R 库房
S 库房

P 库房
0.5
----
Q 库房
0.5
0.3
R 库房
1.0
0.5
S 库房
0.2
0.2
客户 C1 客户 C2 客户 C3 客户 C4 客户 C5 客户 C6
1.0 ---1.5 2.0 ---1.0
2.0 ----------------
---1.5 0.5 1.5 ---1.0
1.0 0.5 0.5 1.0 0.5 ----
---1.5 2.0 ---0.5 1.5
------0.2 1.5 0.5 1.5
注:单位元/吨;划“----”表示无供货关系. 某些客户表示喜欢由某厂或某库房供货.计有: C1-------- A 市厂 C2-------- P 库房 C5--------Q 库房 C6--------R 库房或 S 库房 A 市厂月供货量不能超过 150 千吨，B 市厂月供货量不能超过 200 千吨。各库房的月最大流通
量千吨数为 表二
库房
P
Q
R
S
流通量
70
50
100
40
各客户每月所必须满足的供货量为（单位：千吨） 表三
客户
C1
C2
C3
C4
C5
C6
要求货量
50
10
40
35
60
20
现假设可以在 T 市和 V 市建新库房，和扩大 Q 市的库房，而库房的个数又不能多于 4 个，必要时可 关闭 P 市和 S 市的库房。

建新库房和扩建 Q 市库房的费用（计入利息）摊至每月为下表所列值（万元），它们的潜在的 月流通量（千吨）也列于表中
表四
库房
月费用
流通量
T
1.2
30
V
0.4
25
Q（扩建）
0.3
20
关闭 P 市库房月省费用 1 万元；关闭 S 市库房月省 0.5 万元。 涉及新库房的配送费用单价（元/吨）见下表
表五
供货 A
受货
B
T
V
T
0.6
0.4
V
0.4
0.3
C1
1.2
----
C2
0.6
0.4
C3
0.5
----
C4
----
0.5
C5
0.3
0.6
C6
0.8
0.9
2、问题分析
随着经济的发展、交通网络的不断健全以及各项科技的进步。使得各个行业竞争激烈，生产商 要在满足客户要求与尽量减少生产成本之间面临更复杂决策。在整个配送问题中，所有的对象有三 种，一种就是厂房，它是货物的产源地分别地处Ａ、Ｂ两个市，它所生产的货物，可以直接运给客 户，也可以放到库房里存放；第二种就是库房，用于存放来自于Ａ、Ｂ两个厂房的生产物以及将货 物配送给它的顾客，这种库房分别位于Ｐ、Ｑ、Ｒ、Ｓ市；第三种就是客户，接收由工厂或库房提 供的货物；
问题一、在配送过程中，我们需要建立一个数学模型来计算如何配货公司的运输费用最低，如何配 送货物，既能满足客户的要求，又能为公司节约足够的资金。当然还要考虑到增加工厂和库房的生 产能力对配送费用的影响，费用单价、客户对供应货物的最低要求以及工厂和库房生产能力各微小 变化对配货方案的影响等因素来进行方案设计。设计出来的方案还要能体现出公司在什么样的改进

下能获得更高的经济效益。
可以用数学模型来建立最优解，进而解决设计方案的建立。
问题二、在问题一得基础上几乎没什么变化，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 供货者
Ｑ
Ｒ
Ｓ
Ａ
X11
X12
X13
X14
Ｂ
X21
X22
X23
X24
工厂和各库房向客户的供应量：
受货者 A
B
P
Q
R
S
供货者
C1
X011
X012
X013
X014
X015
X016
C2
X021
X022
X023
X024
X025
X026
C3
X031
X032
X033
X034
X035
X036
C4
X041
X042
X043
X044
X045
X046
C5
X051
X052
X053
X054
X055
X056
C6
X061
X062
X063
X064
X065
X066
模型要求公司在配货时的最小运输费用,即： 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 库房
0.5
----
Q 库房
0.5
0.3
R 库房
1.0
0.5
S 库房
0.2
0.2
客户 C1 客户 C2 客户 C3 客户 C4 客户 C5 客户 C6
1.0 ---1.5 2.0 ---1.0
2.0 ----------------
---1.5 0.5 1.5 ---1.0
1.0 0.5 0.5 1.0 0.5 ----
---1.5 2.0 ---0.5 1.5
------0.2 1.5 0.5 1.5
注:单位元/吨;划“----”表示无供货关系. 工厂向各库房的供应量：
受货者 p
Q
R
S
供货者
Ａ
X11
X12
X13
X14
Ｂ
X21
X22
X23
X24
可以得出以下结果： 位于 A 地的厂向 P、Q、R、S 地库房供应货物所需运费：

位于 B 地的厂向 P、Q、R、S 地库房供应货物所需运费： 各库房货物的月流通量：
库房
P
Q
R
S
流通量
70
50
100
40
则得之，X11<=70,X12+X22<=50,X13+X23<=100,X14+X24<=40; 各客户每月所必须满足的供货量为（单位：千吨）
客户
C1
C2
C3
C4
C5
C6
要求货量
50
10
40
35
60
20
工厂和各库房向客户的供应量：
受货者 A
B
P
Q
R
S
供货者
C1
X011
X012
X013
X014
X015
X016
C2
X021
X022
X023
X024
X025
X026
C3
X031
X032
X033
X034
X035
X036
C4
X041
X042
X043
X044
X045
X046
C5
X051
X052
X053
X054
X055
X056
C6
X061
X062
X063
X064
X065
X066
问题一、可以得到以下线性方程：
客户 C1 可分别由工厂 A、B 以及库房 Q 三方供应货物，且配货所需单价分别为 1.0、2.0、1.0，则 得为 C1 配货的花费为：
C1 客户的月需求为：
X X X
50 ;
011
012
014
客户 C2 可分别库房 P、Q、R 三方供应货物，且配货所需单价分别为 1.5、0.5、1.5，则得为 C 配 货的花费为：
C2 客户的月需求量为：
X X X
10 ;
023
024
025

客户 C3 可分别由工厂 A 以及库房 P、Q、R、S 五方供应货物，且配货所需单价分别为 1.5、0.5、 0.5、2.0、0.2，则得为 C3 配货的花费为：
C3 客户的月需求量为：
X X X X X
40 ；
031
033
034
035
036
客户 C4 可分别由工厂 A 以及库房 P、Q、S 四方供应货物，且配货所需单价分别为 2.0、1.5、1.0、 1.5，则得为 C4 配货的花费为：
C4 客户的月需求量为：
X X X X
35 ；
041
043
044
046
客户 C5 可分别由库房 Q、R、S 三方供应货物，且配货所需单价分别为 0.5、0.5、0。5，则得为 C5 配货的花费为：
C5 客户的月需求量为：
X X X
60 ；
054
055
056
客户 C6 可分别由工厂 A 以及库房 P、R、S 四方供应货物，且配货所需单价分别为 1.0、1.0、1.5、 1.5，则得为 C6 配货的花费为：
客户 C6 的月需求量为：
X X X X
20 ;
061
063
065
066
根据上面式子中的所有变量以及题给意思可知： A 厂的每月最大供应量为：
X X X X X X X X
150 ;
11
12
13
14
011
031
041
061
B 厂的每月最大供应量为：
X X X X 200 ;
22
23
24
012
从问题考虑，使费用最小，即为库房的存储量供完货后为 0 是最好的。P 库房的货物量有 A 厂，运 出的货物量为客户 C2、C3、C4、C6，由约束条件如下：
X X X X X 0；
023
033
043
063
11
Q 库房的货物量有 A、B 厂，运出的货物量为客户 C1、C2、C3、C4、C5，由约束条件如下：
X X X X X X X 0 ；
014
024
034
044
054
12
22
R 库房的货物量有 A、B 厂，运出的货物量为客户 C2、C3、C5、C6，由约束条件如下：
X X X X X X 0 ；
025
035
055
065
13
23

S 库房的货物量有 A、B 厂，运出的货物量为客户 C3、C4、C5、C6，由约束条件如下：
X X X X X X 0 ；
036
046
056
066
14
24
总结以上的模型,可得配货的最小运输费用问题实际上为一个线性规划模型: 目标函数： min=0.5*X11+0.5*X12+1.0X13+0.2*X14+0.3*X22+0.5*X23+0.2*X24+X011+2.0*X012+X014+1.5*X023+0.5*X024+1.5 *X025+1.5*X031+0.5*X033+0.5*X034+2.0*X035+0.5*X036+2.0*X041+1.5*X043+X044+1.5*X046+0.5*X054+0.5*X055+0. 5*X056+X061+X063+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+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=0.5*X11+0.5*X12+1.0X13+0.2*X14+0.3*X22+0.5*X23+0.2*X24+X011+2.0*X012+X014+1.5*X023+0.5*X024+1.5 *X025+1.5*X031+0.5*X033+0.5*X034+2.0*X035+0.5*X036+2.0*X041+1.5*X043+X044+1.5*X046+0.5*X054+0.5*X055+0. 5*X056+X061+X063+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+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 市库房的费用（计入利息）摊至每月为下表所列值（万元），它们的潜在的 月流通量（千吨）也列于表中
库房
月费用
流通量
T
1.2
30
V
0.4
25
Q（扩建）
0.3
20
关闭 P 市库房月省费用 1 万元；关闭 S 市库房月省 0.5 万元。 涉及新库房的配送费用单价（元/吨）见下表
供货 A
受货
B
T
V
T
0.6
0.4
V
0.4
0.3
C1
1.2
----
C2
0.6
0.4
C3
0.5
----
C4
----
0.5

C5
0.3
0.6
C6
0.8
0.9
问题二基于问题一的内容故可得出如下： 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 房的调整费用为 0.3*a2; 则对 R 房的调整费用为-0.5+0.5*a3; 则对 T 房的调整费用为 1.2*a4; 则对 V 房的调整费用为 0.4*a5; 由题意可知目标函数为货物配送费用与库房调整费用相加得出：
Min=0.5*X11+0.5*X12+X13+0.2*X14+0.6*X15+0.4*X16+y11+1.5*y13+2*y14+y16+0.3*X22+0.5*X23+0.2*X24+0.4* X25+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 约束条件： 有问题一可以轻松得出： 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=0.5*X11+0.5*X12+X13+0.2*X14+0.6*X15+0.4*X16+y11+1.5*y13+2*y14+y16+0.3*X22+0.5*X23+0.2*X24+0.4* X25+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 X11+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=0.5*X11+0.5*X12+1.0X13+0.2*X14+0.3*X22+0.5*X23+0.2*X24+X011+2.0*X012+X014+1.5*X023+0.5*X024+1.5 *X025+1.5*X031+0.5*X033+0.5*X034+2.0*X035+0.5*X036+2.0*X041+1.5*X043+X044+1.5*X046+0.5*X054+0.5*X055+0. 5*X056+X061+X063+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+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=0.5*X11+0.5*X12+X13+0.2*X14+0.6*X15+0.4*X16+y11+1.5*y13+2*y14+y16+0.3*X22+0.5*X23+0.2*X24+0.4* X25+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 X11+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); end 7、结论解释{这段本来是 lingo 运行的结果，现把结果给删去了，用表格体现的结果} 问题一：

该问题的最小费用为： 210.5000单位；
对应的配货方案见下表：
受货者
A 市厂
P 库房
40
B 市厂
供货者 P 库房 Q 库房
Q 库房
50
R 库房
15
S 库房
40
客户 C1
50
客户 C2
客户 C3
客户 C4
客户 C5
客户 C6
20
10 40
35 5
R 库房 15
S 库房 40
A-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千吨；
问题二：
该问题的最小费用为：106.4000单位；
其具体的运输调运方案为：
受货者
A 市厂 B 市厂
Q 库房
20
供货者 Q 库房 R 库房
R 库房
40
S 库房
V 库房
25
客户 C1
50
S 库房
V 库房

客户 C2
10
客户 C3
40
客户 C4
10
25
客户 C5
60
客户 C6
20
关闭 S 仓库，新建 V 仓库；
A-R：40 千吨，运往 C1：50 千吨，运往 C6:20 千吨；
B-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 的 费用也是相同的，所以在这里的配送方案也不唯一。所以在实际中，可以结合实际情况进行调配
8、模型优缺点
模型的优点：
（1） 模型将复杂的配货过程简单和直观化，利用简单的线性规划进行建模，再结合计算机进行 运算，整个过程简单而易操作。
（2） 模型具有一般性和普遍性，适合任何情况下的最优配货方案的设计。（3） 该模型的方法 还可以推广大更大的经济领域中去，具有很大的可塑性。
模型的不足：
（1）模型中的变量和参数太多，在输入的时候可能比较麻烦。
（2）在模型的求解过程中，所给出的运货为零的方案在实际计算中并非真正为 0，而是计算机 在计算过程中采取了舍入的方法，这些对模型结果的影响并不大，因为数据足够小，而且在实际中， 如果运送的货物很少的话，对资源反而是一种浪费。
9、模型推广
1． 客户满意度：
在模型二中，我们对客户喜好程度的处理是采取将近 100%地满足客户的要求，而客户对某个厂房 和库房的喜好也是 100%，但是，在实际中并非如此。所以在处理客户的选择的时候应该引入客户 对某厂房和某库房货物的满意度，即客户对该厂房或库房的喜好程度，所要求的货物在总的需求量 中的比例等，来作为对货物配送中的一个主要的考虑因素。

2． 公司的收益： 在整个模型中，我们只考虑了怎么使运输费用最少的问题，却没有考虑货物卖出之后能挣多少的问 题，一批货物运到这一客户和运到那一客户收益有什么不同，货物的多少对收益又有什么影响，在 实际中，收益才是配货的首要考虑问题。所以，在模型中加入这一因素，就可以实现对整个公司经 济活动的宏观调控。
10、附件（各种程序等） 问题一中的 lingo 程序： min=0.5*X11+0.5*X12+1.0X13+0.2*X14+0.3*X22+0.5*X23+0.2*X24+X011+2.0*X012+X014+1.5*X023+0.5*X024+1.5 *X025+1.5*X031+0.5*X033+0.5*X034+2.0*X035+0.5*X036+2.0*X041+1.5*X043+X044+1.5*X046+0.5*X054+0.5*X055+0. 5*X056+X061+X063+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+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=0.5*X11+0.5*X12+X13+0.2*X14+0.6*X15+0.4*X16+y11+1.5*y13+2*y14+y16+0.3*X22+0.5*X23+0.2*X24+0.4* X25+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 X11+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); end