第5章运输与指派问题文档格式.docx
《第5章运输与指派问题文档格式.docx》由会员分享,可在线阅读,更多相关《第5章运输与指派问题文档格式.docx(29页珍藏版)》请在冰豆网上搜索。
23
【解】设xij(i=1,2,3;
j=1,2,3,4)为i个产粮地运往第j个需求地的运量(万吨),如此取得以下运输问题的数学模型:
(1)使总的运输费用最小,那么目标函数为
实际总运费等于Z乘以10000。
(2)各产粮地的供给量与运出量的平稳方程
(3)供给各需求地的供给量与需要量的平稳方程
(4)粮食的运量应大于或等于零(非负要求),即
有些问题表面上与运输问题没有多大关系,其模型的数学结构与例运输问题模型形式相同,咱们把这种模型都称为运输模型。
5.1.2模型特点
运输问题的数学模型有它的专门性。
假设有m个产地,记作A1,A2,A3,…,Am,生产某种物资,可供给的产量别离为a1,a2,…,am;
有n个销地(需求地),记作B1,B2,…,Bn,其需要量别离为b1,b2,…,bn;
供需平稳,即
。
从第i个产地到j个销地的单位物资的运费(运价)为cij,在知足各地需要的前提下,求总运输费用最小的调运方案。
设xij(i=1,2,…,m;
j=1,2,…,n)为第i个产地到第j个销地的运量,那么运输问题的数学模型为
当目标是利润时,式(5.1a)改成求最大值;
当总供给量大于总需求量时,式改成“≤”约束;
当总供给量小于总需求量时,式(5.1c)改成“≤”约束。
式()是线性计划模型,在供需平稳条件下,有m+n个等式约束mn个变量,约束条件的系数矩阵A有m+n行mn列。
目标函数由运价矩阵Cm×
n与变量矩阵Xm×
n对应元素相乘求和组成。
在第1章的线性计划模型中,假设系数矩阵Am×
n知足r(A)=m,即m个约束方程是彼此独立的,因此有m个基变量(基矩阵的秩等于m),但运输问题的数学模型的基变量数不是等于m+n,而是m+n-1。
【定理】设有m个产地n个销地且产销平稳的运输问题,那么基变量数为m+n-1。
为了在mn个变量中找出m+n-1个变量作为一组基变量,确实是要在A中找出m+n-1个线性无关的列向量,下面引用闭回路的概念寻觅这些基变量。
为一个闭回路,集合中的变量称为闭回路的极点,相邻两个变量的连线为闭回路的边。
在表5-2及表5-3中的变量组组成两个闭回路。
表5-2
x11
x12
x32
x33
A4
x41
x43
表5-3
B4
B5
x23
x25
x31
x35
x42
例如在表5-3中取变量组:
A不能组成一条闭回路,但A中包括有闭回路
,B的变量数是奇数,显然不是闭回路,也不含有闭回路;
C是一条闭回路,假设把C从头写成
不难看出C仍是一条闭回路。
假设变量组
中某一个变量是它所在行或列中显现的唯一变量,那么称那个变量是关于变量组的孤立点。
在集合A中x12位于第一行第二列,A中第二列只有变量x12,故x12是变量组A的孤立点。
变量组C中没有孤立点。
变量组B中x21和x33是孤立点。
很显然,假设一个变量组不包括任何闭回路,那么变量组必有孤立点,如变量组B。
一条闭回路中必然没有孤立点,如变量组C。
有孤立点的变量组中不必然没有闭回路,如变量组A。
【定理】假设变量组B包括有闭回路
,那么B中的变量对应的列向量线性相关。
【定理】m+n-1个变量组组成基变量的充要条件是它不包括任何闭回路。
定理告知了一个求基变量的简单方式,同时也能够判定一组变量是不是能够作为某个运输问题的基变量。
这种方式是直接在运价表中进行的,不需要在系数矩阵A中去寻觅,从而给运输问题求初始基可行解带来极大的方便
例如,m=3,n=4,将运价Cij和运量xij放在同一张表中,如表5-4所示。
表5-4
Bj
Ai
ai
C11
C12
C13
C14
a1
x13
x14
C21
C22
C23
C24
a2
x21
x22
x24
C31
C32
C33
C34
a3
x34
bj
b1
b2
b3
b4
那个运输问题的基变量数量是3+4-1=6。
变量组
中有7个变量,显然不能组成一组基变量,又如
中有6个变量,但包括有一条闭回路
,故不能组成一组基变量。
变量组
有6个变量且不含有任何闭回路,故能够组成此问题的一组基变量。
运输单纯形法
运输单纯形法(或称表上作业法)是直接在运价表上求最优解的一种方式,条件是:
问题求最小值、产销平稳和运价非负。
大体步骤是:
第一步:
求初始基行可行解(初始调运方案),经常使用的方式有最小元素法、元素差额法(Vogel近似法)、左上角法。
第二步:
求查验数并判定是不是取得最优解,经常使用求查验的方式有闭回路法和位势法,假设目标函数取最小,当非基变量的查验数λij全都非负时取得最优解(关于maxZ有λij≤0时最优),假设存在查验数λlk<
0,说明尚未达到最优,转第三步。
5.2.1初始基可行解
(1)最小元素法最小元素法的思想是就近优先输送,即最小运价Cij对应的变量xij优先赋值
,然后再在剩下的运价中取最小运价对应的变量赋值并知足约束,依次下去,直到最后取得一个初始基可行解。
(2)元素差额法(Vogel近似法)最小元素法只考虑了局部运输费用最小,对整个产销系统的总运输费用来讲可能离最优值较远,有时为了节省某一处的运费,而在其它处可能运费专门大。
元素差额法对最小元素法进行了改良,考虑到产地到销地的最小运价和次小运价之间的差额,若是差额专门大,就选最小运价处先调运,不然会增加总运费。
例如下面两种运输方案
前一种按最小元素法求得,总运费Z1=10×
8+5×
2+15×
1=105,后一种方案考虑到C11与C21之间的差额是8-2=6,若是不先调运x21,到后来就有可能x11≠0,如此会使总运费增加较大,从而先调运x21,再是x22第二是x12,这时总运费Z2=10×
5+15×
2+5×
1=85<
Z1。
基于以上方式,元素差额法求初始大体可行解的步骤是:
求出每行次小运价与最小运价之差,记为ui,i=1,2,…,m;
同时求出每列次小运价与最小运价之差,记为vj,j=1,2,…,n;
找出所有行、列差额的最大值,即L=max{ui,vi},差额L对应行或列的最小运价处优先调运;
第三步:
这时必有一列或一行调运完毕,在剩下的运价中再求最大差额,进行第二次调运,依次进行下去,直到最后全数调运完毕,就取得一个初始调运方案。
用元素差额法求得的大体可行解更接近最优解,因此也称为近似方案。
(3)左上角法
左上角法(亦称西北角法)是优先从运价表的左上角的变量赋值,当行或列分派完毕后,再在表中余下部份的左上角赋值,依次类推,直到右下角元素分派完毕。
5.2.2求查验数
求最小值的运输问题的最优判别准那么是:
所有非基变量的查验数都非负,那么运输方案最优(即为最优解)。
求查验数的方式有两种,闭回路法和位势法。
(1)闭回路法闭回路法求查验数求某一非基变量的查验数的方式是:
在大体可行解矩阵中,以该非基变量为起点,以基变量为其它极点,找一条闭回路,由起点开始,别离在极点上交替标上代数符号+、-、+、-、…-,以这些符号别离乘以相应的运价,其代数和确实是那个非基变量的查验数。
一样地,当某个非基变量xij增加一个单位时,总费用的改变量ΔZ=λij,当λij<
0时,说明总费用下降了,这时增加xij的值能够降低总费用,当λij>
0时,总费用增加,说明增加xij的值总费用也随之上升。
当所有非基变量的查验数全数大于零时,说明不能增加任何非基变量的值,即不能将非基变量换入变成基变量,不然总费用增加,这时的大体可行解确实是最优解,其费用最小。
当某个非基变量的查验数λlk<
0时,说明能够增加xlk的值,使总费用下降,即将xlk由非基变量换入成基变量,这时的大体可行解也就不是最优解,需要对运输方案进行调整。
只要求得的基变量是正确的且数量为m+n-1,那么某个非基变量的闭回路存在且唯一,因此查验数唯一。
(2)位势法位势法求查验位势法求查验数是依照对偶理论推导出来的一种方式。
运输问题的对偶模型是
加入松驰变量λij将约束化为等式
记原问题基变量XB的下标集合为I,由第二章对偶性质知,原问题xij的查验数是对偶问题的松弛变量λij,
时λij=0,因此有
方程组(5.2a)有m+n-1个方程,有m+n个未知变量ui及vj,有一个自由变量,一样地令u1=0就可取得ui及vi的一组解,再由式取得非基变量的查验数,称ui、vi为运输问题关于基变量组
的对偶解,或称位势(ui为行位势,vi为列位势)。
不同的基变量组
或自由变量的取值不同,取得不同的位势,ui及vi有无穷多组解,但对同一组基变量来讲,所求得的查验数是唯一的,并与闭回路法求得的查验数相同,这种求查验数的方式称为位势法。
5.2.3调整运量
当某个查验数小于零时,基可行解不是最优解,总运费还能够下降,这时需调整运输量,改良原运输方案,改良运输方案的步骤是:
确信进基变量,
进基;
确信出基变量,在进基变量xlk的闭回路中,标有负号的最小运量作为调整量θ,θ对应的基变量为出基变量,并打上“×
”以示作为非基变量。
第三步:
调整运量,在进基变量的闭回路中标有正号的变量加上调整量θ,标有负号的变量减去调整量θ,其余变量不变,取得一组新的基可行解,然后求所有非基变量的查验数从头查验。
以上调整运量的方式为闭回路法,既可使总运费下降,又使新的大体解可行。
在第二步确信出基变量时,当显现两个或两个以上最小运量θ,在其中任选一个作为非基变量,其它θ对应的变量仍作为基变量,运量为零,取得退化大体可行解。
【例】求以下运输问题的最优解
【解】用最小元素法求得初始大体可行解如表5-5。
表5-5
70
×
40
30
80
45
35
12
14
25
15
65
50
190
求非基变量的查验数,用闭回路法得
因为有4个查验数小于零,因此这组大体可行解不是最优解。
=-4最小,对应的非基变量x11进基。
x11的闭回路是
标负号的变量是x12、x33、x21,取运量最小值即
x33最小,x33是出基量,调整量θ=15,在x11的闭回路上x11、x32、x23别离加上15,x12、x33、x21别离减去15,而且在x33处打上记号“×
”作为非基变量,其余变量的值不变,调整后取得一组新的基可行解,见表。
表
从头求所有非基变量的查验数得
λ13=3,λ22=0,λ24=7,λ31=1,λ33=4,λ34=-1
λ34=-1<
0,说明尚未取得最优解,x34进基,在x34的闭回路中,标负号的变量x14和x32,调整量为
因此x14出基变成非基变量。
x12、x34别离加上30,x14和x32减去30,其余变量不变,取得一组新的基可行解,见表。
55
再求非基变量的查验数
λ13=3,λ14=1,λ22=0,λ24=8,λ31=1,λ33=4
所有查验数λij≥0因此取得最优解
最小运费
由λ22=0知有多重解,求另一最优解的方式是在x22的闭回路{x22,x21,x11,x12}上任意调整(维持可行),目标值不变,但只能取得一个大体最优解,其它都不是大体最优解,例如x22=5调整后取得最优解
但不是大体最优解
运输单纯形法计算进程中,运量调整后必需将所有非基变量的查验数从头求一次,由于基变量的位置改变了,原查验数也可能发生改变,或由大于等于零变成小于零。
【例】有四项工作指派给甲、乙两人完成,每人完成两项工作。
两人完成各项工作的时刻(小时)见表5-8,如何安排工作使总时刻最少。
表5-8
A
B
C
D
甲
20
乙
16
【解】设xij(i=1,2;
j=1,2,3,4)为第i人完成第j项工作的状态
数学模型为
用求解0-1计划的方式能够取得问题的最优解。
采纳本章的方式求解更为简单。
写出表5-8的平稳运输表5-9,用运输单纯形法求解取得最优表5-10。
产量
销量
最优的工作分派是:
甲完成工作C和D,乙完成工作A和B,总时刻Z=47(小时)。
例的任务分派问题能够推行到更一样的情形。
有m个人完成n项工作,每项工作一人完成,第i人完成jr项工作,j1+j2+…+jm=n,用运输单纯形法求解。
5.2.4最大值问题
当运输模型求最大值时有两种方式求解。
(1)将极大化问题转化为极小化问题。
设极大化问题的运价表为C=(Cij)m×
n,用一个较大的数M(一样令M=max{Cij})去减每一个Cij取得矩阵C'
=(C'
ij)m×
n,其中C'
ij=M-Cij≥0,将C/作为极小化问题的运价表,用表上用业法求出最优解,目标函数值为
例如,以下矩阵C是Ai到Bj的吨千米利润,运输部门如何安排运输方案使总利润最大。
用最小元素法求初始方案得
λ11=8,λ12=4,λ21=2,λ23=2全数非负,取得最优运输方案X,最大利润Z=8×
9+10×
10+6×
4=240。
(2)所有非基变量的查验数λij≤0时最优。
求初始运输方案可采纳最大元素法或西北角法。
如上例,用最大元素法取得初始运输方案
求查验数:
λ11=-8,λ12=-4,λ21=-2,λ23=-2,全数非正,取得最优解运输方案,结果与第一种方式相同。
5.2.5不平稳运输问题
当总产量与总销量不相等时,称为不平稳运输问题。
这种运输问题在实际中常常碰着,它的求解方式是将不平稳问题化为平稳问题求解。
(1)产大于销时,即有
由于总产量大于总销量,必有部份产地的产量不能全数输送完,必需当场库存,即每一个产地虚设一个仓库,库存量为xi,n+1(i=1,2,…,m),总的库存量为
bn+1作为一个虚设的销地Bn+1的销量。
各产地Ai到Bn+1的运价为零,即Ci,n+1=0,(i=1,…,m)。
那么平稳问题的数学模型为
具体求解时,只在运价表右端增加一列Bn+1,运价为零,销量为bn+1即可。
(2)当销大于产时,即
由于总销量大于总产量,故必然有些需求地不完全知足需求,这时虚设一个产地Am+1,产量为
xm+1,j是Am+1运到Bj的运量,也是Bj不能知足需要的数量。
Am+1到Bj的运价为零,即Cm+1,j=0(j=1,2,…,n),那么平稳问题的数学模型为
具体计算时,在运价表的下方增加一行Am+1,运价为零,产量为am+1即可。
上述两种情形将不等式化为等式的进程,等价于加入松驰变量x1,n+1,x2,n+1,…,xm,n+1
及xm+1,1,xm+1,2,…,xm+1,n,松弛变量在目标函数中的系数为零,因此目标函数不变。
【例】求表5-11极小化运输问题的最优解。
表5-11
60
—
11
160\180
【解】因为
,表中A2不可达B1,用一个专门大的正数M表示运价C21。
虚设一个销量为b5=180-160=20的销地B5,Ci5=0,i=1,2,3,4。
表的右边增添一列取得表5-12。
用元素差额法求初始基可行解,取得表5-13。
表5-12
M
180
表5-13
请读者自行验证,所有查验数λij≥0,表5-13的运输方案最优,最小运费
产地A4还有20个单位没有输送出去。
在例中,若是用最小元素法求初始基可行解时,因为B5列全为零,能够优先输送B5,也能够最后输送B5,不阻碍最优方案。
5.2.6需求量不确信的运输问题
【例】在表5-11给出的运输问题中,假定B1的需要量是20到60之间,B2的需要量是50到70,试求极小化问题的最优解。
【解】
(1)总产量为180,B1,…,B4的最低需求量20+50+35+45=150,这时属产大于销;
(2)B1,…,B4的最高需求是60+70+35+45=210,这时属销大于产;
(3)虚设一个产地A5,产量是210-180=30,A5的产量只能供给B1或B2。
(4)将B1与B2各分成两部份,
,
的需求量是20,
的需求量是40,
与
的需求量别离是50与20,因此
必需由A1,…,A4供给,
可由A1、…、A5供给。
(5)上述A5不能供给某需求地的运价用大M表示,A5到
的运价为零。
取得表的产销平稳表。
表5-14
A5
210
用运输单纯形法计算,取得表5-15的最优方案。
表5-15
表5-15中
是基变量,说明这组解是退化大体可行解,空格处的变量是非基变量。
B1,B2,B3,B4实际收到产品数量别离是50,50,35和45个单位。
运输单纯形法求运输问题的最优解要比用一般单纯法求解简单,另外有些问题能够用图上作业法。
运输模型的应用
运输模型是线性计划模型的一种特例,在生产实践中取得普遍的应用。
【例】DF公司在接下来的三个月内每一个月都要依照销售合同生产出两种产品。
这两种产品利用相同的设备并需要投入相同的生产能力。
每一个月可供利用的生产和存储设备都会发生转变。
因此生产能力、单位生产本钱和单位存储本钱每一个月都不相同,有必要在某些月中多生产一种或多种产品并存储起来以备需要的时候利用。
关于每一个月来讲,表中给出了在正常时刻(RegularTime,缩写为RT)和加班时刻(OverTime,缩写为OT)内能够生产这两种产品的总数,依照合同需要生产的数量,在正常时刻和加班时刻内的单位产品本钱和每件产品贮存到下一个月的贮存本钱。
两种产品的数量用“/”区分开来,产品1在“/”的左侧而产品2在“/”的右边。
月
最大生产总量
产品1/产品2
销售
单位生产成本(1000元/件)
单位储存成本(1000元/件)
RT
OT
5/3
3/5
4/4
15/16
17/15
19/17
18/20
20/18
22/