1、2024/5/13 周一1运筹学运筹学OPERATIONS RESEARCH.2024/5/13 周一2第三章第三章 运输问题运输问题n运输问题的数学模型运输问题的数学模型n表上作业法表上作业法n产销不平衡的运输问题及应用产销不平衡的运输问题及应用.2024/5/13 周一31 1 运输问题的典例及数学模型运输问题的典例及数学模型一、一、引例引例某公司从三个产地某公司从三个产地 ,将产品运往四个销地将产品运往四个销地 ,各产地的产量,各销地的销量,及各产地往各销,各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。应如何调运可使运费最小?地的运费单价如表所示。应如何调运可使运费最小
2、?.2024/5/13 周一4解:解:从表中可知:总产量从表中可知:总产量 =总销量。这是一个产销平衡的总销量。这是一个产销平衡的 运输问题。假设运输问题。假设 表示从产地表示从产地 运往销地运往销地 的产的产 品数量,品数量,建立如下表格:建立如下表格:于是可建立如下的数学模型于是可建立如下的数学模型:.2024/5/13 周一5目标函数目标函数:约束条件:约束条件:产量约束产量约束销量约束销量约束2020.2024/5/13 周一6设有设有m个产地,分别为个产地,分别为 ;n 个销地,分别是个销地,分别是 ;从产地从产地 运往销地运往销地 的单位运价是的单位运价是 ,运量,运量 是产地是产
3、地 的产量;的产量;是销地是销地 的销量。的销量。二、二、一般运输问题数学模型一般运输问题数学模型则该运输问题的模型如下:则该运输问题的模型如下:.2024/5/13 周一7说明说明:当:当 时,称其为产销平衡的运输问题,时,称其为产销平衡的运输问题,否则产销不平衡。否则产销不平衡。.2024/5/13 周一8说明说明:从上述模型可以看出:从上述模型可以看出:(1)这是一个线性规划的模型;)这是一个线性规划的模型;(2)变量有)变量有mn个;个;(3)约束条件有)约束条件有 m+n 个;个;(4)系数矩阵非常稀疏;系数矩阵的秩一般为(系数矩阵非常稀疏;系数矩阵的秩一般为(m+n-1),m+n-
4、1),而非而非m+n。若直接用单纯形法求解,显然单纯形表比较庞大,于是在若直接用单纯形法求解,显然单纯形表比较庞大,于是在单纯形法的基础上创建了表上作业法求解运输问题这一特单纯形法的基础上创建了表上作业法求解运输问题这一特殊的线性规划问题殊的线性规划问题 .2024/5/13 周一9 从从第第一一节节的的运运输输问问题题的的数数学学模模型型可可知知,运运输输问问题题实实际际上上也也属属于于线线性性规规划划,但但由由于于运运输输问问题题的的特特殊殊性性(变变量量个个数数较较多多,系系数数矩矩阵阵的的特特点点),如如果果用用单单纯纯形形表表格格方方法法迭迭代代,计计算算量量很很大大。今今天天介介绍
5、绍的的 “表表上上作作业业法法”,是是针针对对运运输输问问题题的的特特殊殊求求解解方方 法法,实实 质质 还还 是是 单单 纯纯 形形 法法,但但 减减 少少 了了 计计 算算 量量。表表上上作作业业法法适适用用于于求求解解产产销销平平衡衡的的运运输输问问题题。(产产销销不不平平衡衡的的问问题题可可转转化化为为平平衡衡问问题题)2 2 运输问题的表上作业法运输问题的表上作业法.2024/5/13 周一10表上作业法表上作业法 一般步骤一般步骤:1、找出初始基本可行解;、找出初始基本可行解;2、检查各非基变量的检验数,是否达到最优性条件,若达到,则得最优、检查各非基变量的检验数,是否达到最优性条
6、件,若达到,则得最优解;否则解;否则 转第三步;转第三步;3、确定出基变量、进基变量,用闭回路方法进行调整,得到新的基可、确定出基变量、进基变量,用闭回路方法进行调整,得到新的基可 行解;行解;4、重复第二、第三步,直至得到最优解。、重复第二、第三步,直至得到最优解。.2024/5/13 周一11一、确定初始基本可行解:一、确定初始基本可行解:对于有对于有m m个产地个产地n n个销地的产销平衡问题,有个销地的产销平衡问题,有m m个关于产量个关于产量的约束方程和的约束方程和n n个关于销量的约束方程。表面上,共有个关于销量的约束方程。表面上,共有m+nm+n个个约束方程。约束方程。但由于产销
7、平衡,其模型最多只有但由于产销平衡,其模型最多只有m+n-1m+n-1个独立的约束方个独立的约束方程,所以运输问题实际上有程,所以运输问题实际上有m+n-1m+n-1个基变量个基变量。在。在mnmn的产销的产销平衡表上给出平衡表上给出m+n-1m+n-1个数字格,其相对应的调运量的值即为个数字格,其相对应的调运量的值即为基变量的值。基变量的值。.2024/5/13 周一12 那么在该例中,应有那么在该例中,应有 3+4-1=63+4-1=6个基变量。个基变量。.2024/5/13 周一131.最小元素法最小元素法 最小元素法的思想是就近供应,即对单位运价最小最小元素法的思想是就近供应,即对单位
8、运价最小的变量分配运输量。的变量分配运输量。在表上找到单位运价最小的在表上找到单位运价最小的x x2121,并使,并使x x2121取尽可能大取尽可能大的值,即的值,即x x2121=3,=3,把把A A1 1的产量改为的产量改为1 1,B B1 1的销量改为的销量改为0 0,并,并把把B B1 1列划去。在剩下的列划去。在剩下的3333矩阵中再找最小运价,同矩阵中再找最小运价,同理可得其他的基本可行解。理可得其他的基本可行解。.311310851029471.2024/5/13 周一15表中填表中填有数字的格有数字的格对应于对应于基变量基变量(取值即为格中数字),而取值即为格中数字),而空格
9、空格对应对应的是的是非基变量非基变量(取值为零)(取值为零).n在求初始基本可行解时要在求初始基本可行解时要注意注意的一个问题:的一个问题:当我们取定当我们取定x xijij的值之后,会出现的值之后,会出现A Ai i的产量与的产量与B Bj j的销量都改为零的情的销量都改为零的情况,这时只能划去况,这时只能划去A Ai i行或行或B Bj j列,但不能同时划去列,但不能同时划去A Ai i行与行与B Bj j列。列。(或者在同时划去(或者在同时划去A Ai i行与行与B Bj j列时,在该行或该列的任意空格处填加一列时,在该行或该列的任意空格处填加一个个0 0。)。)这样可以保证填过数或零的
10、格为这样可以保证填过数或零的格为m+n-1m+n-1个,即保证基变量的个数为个,即保证基变量的个数为m+n-1m+n-1个。个。.2024/5/13 周一162.Vogel法法 Vogel法的思想是法的思想是:一地的产品如果不能按照最小运一地的产品如果不能按照最小运费就近供应,就考虑次小运费,这就有差额,差额越大,费就近供应,就考虑次小运费,这就有差额,差额越大,说明不能按最小运费调运时,运费增加得越多。因而差说明不能按最小运费调运时,运费增加得越多。因而差额越大处,就应当采用最小运费调运。额越大处,就应当采用最小运费调运。.2024/5/13 周一17311310294711085.2024
11、/5/13 周一18二二、最优解的判别、最优解的判别 判别解的最优性需要:判别解的最优性需要:计算检验数。计算检验数。方法有两种方法有两种 闭回路:闭回路:是在已给出的调运方案的运输表上从一个代表是在已给出的调运方案的运输表上从一个代表非基变量的空格出发,沿水平或垂直方向前进,遇到代表基非基变量的空格出发,沿水平或垂直方向前进,遇到代表基变量的填入数字的格可转变量的填入数字的格可转9090度(当然也可以不改变方向)继度(当然也可以不改变方向)继续前进,这样继续下去,直至回到出发的那个空格,由此形续前进,这样继续下去,直至回到出发的那个空格,由此形成的封闭折线叫做成的封闭折线叫做闭回路闭回路。一
12、个空格存在唯一的闭回路。一个空格存在唯一的闭回路。1 1.闭回路法闭回路法因为任意非基向量均可表示为基向量的唯一线性组因为任意非基向量均可表示为基向量的唯一线性组合,因此对于任意空格都合,因此对于任意空格都能够找到、并且只能找到能够找到、并且只能找到唯一的唯一的一条闭回路。一条闭回路。.108531131029471108531131029471.2024/5/13 周一20108531131029471从非基变量从非基变量 出发,找到一个闭回路如上表所示。回路有四个出发,找到一个闭回路如上表所示。回路有四个顶点,除顶点,除 外,其余都为基变量。外,其余都为基变量。调整调运量:调整调运量:,运
13、费增加了,运费增加了3 3元;元;,运费减少,运费减少3 3元元 ,运费增加,运费增加2 2元;元;,运费减少,运费减少1 1元元调整后,调整后,总运费增加总运费增加:3-3+2-1=13-3+2-1=1元。元。说明如果让说明如果让 为基变量,运费就会增加,其增加值为基变量,运费就会增加,其增加值1 1作为作为 的的检检验数验数,.2024/5/13 周一21闭回路法计算检验数:闭回路法计算检验数:就是对于代表非基变量的空格就是对于代表非基变量的空格(其调运量为零),把它的调运量调整为(其调运量为零),把它的调运量调整为1 1,由于产销平衡的,由于产销平衡的要求要求,必须对这个空格的闭回路中的
14、各顶点的调运量加上或减必须对这个空格的闭回路中的各顶点的调运量加上或减少少1 1。最后计算出由这些变化给整个运输方案的总运输费带来。最后计算出由这些变化给整个运输方案的总运输费带来的变化。以这个变化的数值,作为各空格(非基变量)的检验的变化。以这个变化的数值,作为各空格(非基变量)的检验数。数。判别最优解准则:判别最优解准则:如果所有代表非基变量的空格的检验如果所有代表非基变量的空格的检验数都大于等于零,则已求得最优解;否则继续改进找出最优解。数都大于等于零,则已求得最优解;否则继续改进找出最优解。.2024/5/13 周一222.2.位势法位势法 (1 1)对运输表上的每一行赋予一个数值)对
15、运输表上的每一行赋予一个数值 ,对,对每一列赋予一个数值每一列赋予一个数值 ,称为行(列,称为行(列)位势。位势。(2 2)行(列)行(列)位势的数值是由基变量的检验数所决位势的数值是由基变量的检验数所决定的,即定的,即基变量基变量要满足:要满足:非基变量非基变量 的检验数就可以用公式的检验数就可以用公式 求出。求出。.311310851029471108531131029471.2024/5/13 周一24 我们先给我们先给u u1 1赋个任意数值,不妨设赋个任意数值,不妨设u u1 1=0=0,则从基变,则从基变量量x x1111的检验数求得的检验数求得 v v3 3=c=c1313-u-
16、u1 1=3-0=3 =3-0=3 。同理可以求得同理可以求得 v v4 4=10=10,u u2 2=-1=-1,等等见上表。,等等见上表。检验数的求法,即用公式检验数的求法,即用公式 ,如如 。311310851029471.2024/5/13 周一25位势法计算检验数:位势法计算检验数:检验数:检验数:又因为基变量的检验数为又因为基变量的检验数为0 0,于是由(,于是由(m+n-1)m+n-1)个基变个基变 量的检验数量的检验数可解出可解出 ,进而计算其他非基变量的检验,进而计算其他非基变量的检验数。数。其中其中 第第i i个分量个分量第第m+jm+j个分量个分量.2024/5/13 周
17、一26三、改进运输方案的办法三、改进运输方案的办法闭回路调整法闭回路调整法当表中的某个检验数小于零时,方案不为最优,需要调整。当表中的某个检验数小于零时,方案不为最优,需要调整。方法是:选取所有负检验数中最小的非基变量作为入基变量,方法是:选取所有负检验数中最小的非基变量作为入基变量,以求尽快实现最优。以求尽快实现最优。(1 1)确定调整量确定调整量:例:取:例:取 ,表明增加一个单位的,表明增加一个单位的 运输量,可使得总运费减少运输量,可使得总运费减少1 1。在以。在以 为出发点的闭回为出发点的闭回路中,找出所有偶数顶点的调运量:路中,找出所有偶数顶点的调运量:,则调整量则调整量 (2 2
18、)调整方法调整方法:把所有闭回路上为偶数顶点的运输量都减:把所有闭回路上为偶数顶点的运输量都减少这个值,奇数顶点的运输量都增加这个值少这个值,奇数顶点的运输量都增加这个值(见下表见下表)。.2024/5/13 周一27311310851029471311310851029471.调整运量后的新方案:调整运量后的新方案:.2024/5/13 周一29311310851029471对上表用位势法进行检验如下表,可知已达最优解。对上表用位势法进行检验如下表,可知已达最优解。.2024/5/13 周一30表上作业法表上作业法的一般步骤:的一般步骤:1 1、用最小元素法或、用最小元素法或VogelVog
19、el法确定初始基可行解;法确定初始基可行解;2 2、判断是否为最优:用闭回路法或位势法计算空格检验数,、判断是否为最优:用闭回路法或位势法计算空格检验数,若所有检验数均非负,则已得到最优解;否则进入第三步;若所有检验数均非负,则已得到最优解;否则进入第三步;3 3、从所有负检验数中选择最小者对应空格作为进基变量,从所有负检验数中选择最小者对应空格作为进基变量,从此点出发作闭回路,确定调整量从此点出发作闭回路,确定调整量 ,奇点处增加,奇点处增加 ,偶,偶点处减少点处减少 。.2024/5/13 周一31例:例:用表上作业法,求解下面的用表上作业法,求解下面的 运输问题运输问题 :解解:用最小元
20、素法确定初始基可行解,如下表所示:用最小元素法确定初始基可行解,如下表所示:.2024/5/13 周一32+-.2024/5/13 周一33+-.2024/5/13 周一34此时所有非基变量的检验数均非负,故已达最优解。此时所有非基变量的检验数均非负,故已达最优解。.2024/5/13 周一35一、产销不平衡的运输问题一、产销不平衡的运输问题例例1 1:某公司从两个产地某公司从两个产地 ,将产品运往三个销地,将产品运往三个销地 ,各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。应如何调运可使运费最小?应如何调运可使
21、运费最小?3 产销不平衡的运输问题及应用产销不平衡的运输问题及应用例例1 1:某公司从两个产地某公司从两个产地 ,将产品运往三个销地,将产品运往三个销地 ,各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。应如何调运可使运费最小?应如何调运可使运费最小?例例1 1:某公司从两个产地某公司从两个产地 ,将产品运往三个销地,将产品运往三个销地 ,各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。应如何调运可使运费最小?应如何调运可使运费最小?例
22、例1 1:某公司从两个产地某公司从两个产地 ,将产品运往三个销地,将产品运往三个销地 ,各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。应如何调运可使运费最小?应如何调运可使运费最小?例例1 1:某公司从两个产地某公司从两个产地 ,将产品运往三个销地,将产品运往三个销地 ,各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。应如何调运可使运费最小?应如何调运可使运费最小?例例1 1:某公司从两个产地某公司从两个产地 ,将产品运往三个销地,将
23、产品运往三个销地 ,各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。应如何调运可使运费最小?应如何调运可使运费最小?例例1 1:某公司从两个产地某公司从两个产地 ,将产品运往三个销地,将产品运往三个销地 ,各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。各产地的产量,各销地的销量,及各产地往各销地的运费单价如表所示。应如何调运可使运费最小?应如何调运可使运费最小?.2024/5/13 周一36易知这个问题中:总产量总销量,即易知这个问题中:总产量总销量,即这时可考虑增加一个假想产地,其产量是(总销量
24、总产量这时可考虑增加一个假想产地,其产量是(总销量总产量150)150)他到各销地的单位运费是于是得到如下的表格:他到各销地的单位运费是于是得到如下的表格:.2024/5/13 周一37例例2 2:某单位有三个区,一区、二区、三区;每年需要生活某单位有三个区,一区、二区、三区;每年需要生活用煤和取暖用煤各用煤和取暖用煤各30003000吨,吨,10001000吨,吨,20002000吨;由河北临城、吨;由河北临城、山西盂县两处煤矿负责供应,两地价格和煤质相同,两煤山西盂县两处煤矿负责供应,两地价格和煤质相同,两煤矿的供应能力分别是矿的供应能力分别是15001500吨,和吨,和40004000吨
25、。由煤矿至该单位吨。由煤矿至该单位三个区的单位运价如表所示。三个区的单位运价如表所示。由于供应能力限制,经研究决定,一区供应量可减少由于供应能力限制,经研究决定,一区供应量可减少03000300吨,吨,二区全部满足,三区不能少于二区全部满足,三区不能少于15001500吨,试求使得运费最小的运吨,试求使得运费最小的运输方案?输方案?.2024/5/13 周一38根据题意,添加虚拟产地后,可作出产销平衡的运价表:根据题意,添加虚拟产地后,可作出产销平衡的运价表:.2024/5/13 周一39例:例:设有三个化肥厂供应四个地区的化肥,假设等量的化设有三个化肥厂供应四个地区的化肥,假设等量的化肥在这
26、个地区的使用效果相同。各厂的产量、各地区的需肥在这个地区的使用效果相同。各厂的产量、各地区的需要量、单位运价如表所示。求出运费最省的调拨方案。要量、单位运价如表所示。求出运费最省的调拨方案。.2024/5/13 周一40解:无论考虑需求的上限还是下限,这都是一个产销不平衡的问题。解:无论考虑需求的上限还是下限,这都是一个产销不平衡的问题。当考虑下限时,产当考虑下限时,产销;当考虑需求上限时,产销;当考虑需求上限时,产销。销。于是可以考虑在满足最低需求的情况下,兼顾最高需求。即将每于是可以考虑在满足最低需求的情况下,兼顾最高需求。即将每 个地区的需求分为个地区的需求分为最低需求最低需求和和(最高
27、(最高-最低)需求最低)需求,最低部分必须,最低部分必须 满足,高出的部分可满足也可不满足。虽然销地满足,高出的部分可满足也可不满足。虽然销地的需求无上限,的需求无上限,但根据生产能力,最多可以给她分配但根据生产能力,最多可以给她分配6060万吨。万吨。另外若将最高需求考虑进来,则需添加虚拟产地另外若将最高需求考虑进来,则需添加虚拟产地D D,其产量应为,其产量应为 5050万吨。万吨。于是可给出如下的产销平衡及运价表。于是可给出如下的产销平衡及运价表。.2024/5/13 周一41.2024/5/13 周一42二、生产与存储问题二、生产与存储问题例例4 4:某厂按照合同规定须于当年每季度末分
28、别提供某厂按照合同规定须于当年每季度末分别提供1010,1515,2525,2020台同一规格的柴油机。已知该厂各季度的生台同一规格的柴油机。已知该厂各季度的生产能力及生产每台柴油机的成本如表所示。又如果生产产能力及生产每台柴油机的成本如表所示。又如果生产出来的柴油机当季不交货,每台每积压一个季度,存储出来的柴油机当季不交货,每台每积压一个季度,存储维护费用维护费用0.150.15万元。要求在完成合同的情况下,使得全万元。要求在完成合同的情况下,使得全年生产(存储)费用最小的决策。年生产(存储)费用最小的决策。.2024/5/13 周一43注意注意:单位运价如何确定:单位运价如何确定故设:故设
29、:表示第表示第 季度生产,用于第季度生产,用于第 季度交货的产品数量。可建立数学季度交货的产品数量。可建立数学模型。模型。.2024/5/13 周一44三、转运问题三、转运问题以本章引例为例,产品从3个产地运往4个销地。现考虑;1、各产地的产品不一定直接运往销地,可以将产品集中后再一起运;2、运往个销地的产品也可先运给其中几个,再转运给其他销地;3、除产、销地外还可以有几个中间转运站,在产地之间、销地之间、或产地与销地之间进行转运;已知单位运价表如下,问在考虑产销地之间直接运输和非直接运输的各种可能方案下,如何安排运输方案,可使总运费最小?.2024/5/13 周一45.2024/5/13 周一461、所有产地、销地、中间转运站都可作为产地或销地,于是得到扩大的运输问题,有11个产地和销地。2、建立扩大的运输问题的单位运价表,其中不可能的运输方案的单位运价用任意大正数M表示。3、中间转运站的产销量相等,而且在运费最小的方案中,不可能出现物资来回倒运现象,故转运点的产销量均为20。4、产地也作为转运站,故其产量应为原产量加上转运量,销量应为20(即转运而来的数量)。销地也是类似。于是有如下的产销平衡表,可用表上作业法求解。.2024/5/13 周一47.