运筹学课件第一节运输问题及其数学模型.pptxVIP

第三章 运输问题模型及其特点求解思路及相关理论求解方法——表上作业法运输问题的推广1 、产销不平衡的运输问题2 、转运问题第一页，共二十七页。?第一节 运输问题及其数学模型一、运输问题的数学模型 1、 运输问题的一般提法： 某种物资有若干产地和销地，现在需要把这种物资从各个产地运到各个销地，产量总数等于销量总数。已知各产地的产量和各销地的销量以及各产地到各销地的单位运价（或运距），问应如何组织调运，才能使总运费（或总运输量）最省？第二页，共二十七页。单位 运价 销 或运距 地产地? B1 B2 … Bn? 产 量 A1 A2 ┆ Am x11c11 x12c12 …x13c1 nX21 c21x22c22 … c2n … … …xm1cm1xm2cm2xmncm n a1 a2 ┆ am? 销 量?? b1 b2 … bn 表 有关信息单位根据具体问题选择确定。第三页，共二十七页。 2、运输问题的数学模型 设xij为从产地Ai运往销地Bj的物资数量（i=1，…m；j=1，…n），由于从Ai运出的物资总量应等于Ai的产量ai，因此xij应满足： 第四页，共二十七页。同理，运到Bj的物资总量应该等于Bj的销量bj，所以xij还应满足： 总运费为： 第五页，共二十七页。运输问题的数学模型第六页，共二十七页。m行n行二、运输问题数学模型的特点1．约束方程组的系数矩阵具有特殊的结构系数矩阵A，形式如下：第七页，共二十七页。每一列只有两个元素为1，其余元素均为0； 列向量Pij =(0,…，0，1，0，…,0,1,0,…0)T，其中两个元素1分别处于第i行和第m+j行。将该矩阵分块，特点是：前m行构成m个m×n阶矩阵，而且第k个矩阵只有第k行元素全为1，其余元素全为0(k=1，…，m）；后n行构成m个n阶单位阵2.运输问题的基变量总数是m + n -1写出增广矩阵第八页，共二十七页。前m行相加之和减去后n行相加之和结果是零向量，说明m+n个行向量线性相关，因此 的秩小于m+n； ？ 由 的第二至m+n行和前n列及 对应的列交叉处元素构成m+n-1阶方阵D 非奇异； ？因此 的秩恰好等于m+n-1，又D本身就含于A中，故A的秩也等于m+n-1 证明系数矩阵A及其增广矩阵的秩都是m+n-1 第九页，共二十七页。第十页，共二十七页。可以证明：m+n个约束方程中的任意m+n-1个都是线性无关的。第十一页，共二十七页。第二节表上作业法求解运输问题表上作业法的基本思想是:先设法给出一个初始方案,然后根据确定的判别准则对初始方案进行检查、调整、改进，直至求出最优方案，如图所示。表上作业法和单纯形法的求解思想完全一致，但是具体作法更加简捷。第十二页，共二十七页。是 确定初始方案( 初 始 基本可行解) 判定是否 最 优？结 束否最优方案 改进调整（换基迭代）运输问题求解思路图第十三页，共二十七页。 一、 初始方案的确定 1、作业表（产销平衡表） 初始方案就是初始基本可行解。 将运输问题的有关信息表和决策变量——调运量结合在一起构成“作业表”（产销平衡表）。 下表是两个产地、三个销地的运输问题作业表。 第十四页，共二十七页。表 运输问题作业表（产销平衡表） 调 销地 运 量产地 B1 B2 B3 产 量 A1 c11 X11 c12 X12 c13 X13 a1 A2 c21 X21 c22 X22 c23 X23 a2 销 量 b1 b2 b3第十五页，共二十七页。其中xij是决策变量，表示待确定的从第i个产地到第j个销地的调运量，cij为从第i个产地到第j个销地的单位运价或运距。2、确定初始方案的步骤：（1）选择一个xij，令xij= min{ai，bj}= 将具体数值填入xij在表中的位置；第十六页，共二十七页。（2）调整产销剩余数量：从ai和bj中分别减去xij的值，若ai-xij=0，则划去产地Ai所在的行，即该产地产量已全部运出无剩余，而销地Bj尚有需求缺口bj-ai；若bj-xij =0，则划去销地Bj所在的列，说明该销地需求已得到满足，而产地Ai尚有存余量ai-bj；（3）当作业表中所有的行或列均被划去，说明所有的产量均已运到各个销地，需求全部满足，xij的取值构成初始方案。否则，在作业表剩余的格子中选择下一个决策变量，返回步骤（2）。第十七页，共二十七页。 按照上述步骤产生的一组变量必定不构成闭回路，其取值非负，且总数是m+n-1个，因此构成运输问题的基本可行解。 对xij的选择采用不同的规则就形成各种不同的方法，比如每次总是在作业表剩余的格子中选择