管理运筹学论文---产销不平衡运输摘要运输问题是运筹学中的一个重要问题，也是物流系统优化中常见的问题，同时也是一种特殊的线性规划问题。
怎么样尽可能的在产地与销地之间减少运输成本和降低运输费用是很多运输公司热切关注的话题。
本文涉及的是一个总产量大于总销量的产销不平衡运输问题，通过对产地与销售地车辆运输的建立模型，在运用表上作业迭代法（最小元素法）求解后，再根据模型用lingo软件编写程序进行求解。
然后对结果进行分析，以及运输问题的延伸。
最后证明用lingo 解决车辆运输的可行性。
关键字：运输问题，产销不平衡，表上作业法，lingo目录一、问题的提出与分析 .................................................. 错误！未定义书签。
1.1问题提出 (3)1.2问题分析 (3)二、模型的建立与基本假设 ............... . (1)2.1模型的建立 (4)2.2基本假设....................................................................... 错误！未定义书签。
三、定义符号说明与表上作业法 (6)四、问题求解..................................................................... 错误！未定义书签。
4.1、Lingo求解模型......................................................... 错误！未定义书签。
4.2、Lingo结果 (9)五、模型结果分析与改进 (10)参考文献............................................................................. 错误！未定义书签。
一、问题的提出与分析1.1问题提出重庆有三家电子厂分别是新普，隆宇和恒华，生产的笔记本电脑将要运向北京，天津，广东，上海四个城市销售，其产量和销售量见下表：（单位：万台）表：1-1北京天津广东上海产量新普 6 2 6 7 30隆宇 4 9 5 3 25恒华8 8 1 5 21销量15 17 22 12 -问：哪种销售方案将会取得最少的运输费用，费用为多少？1.2问题分析图表数据显示产量总和为30+25+21=76万台，销量的总和为15+17+22+12=66万台，说明了此问题是一个总产量大于总销量的运输问题（76>66）。
该问题一方面要求满足北京，天津，广东，上海四个销售地的供货需求，而另一方面又要考虑新普，隆宇和恒华三个产地的运往销售地的运输费用，此外问题不但要求满足销售地分配要足，同时也要保证最大化的减少运输费用。
这里选择何种分配方案，将涉及不同的运输费用，所以其是一个典型的线性规划问题，同时也是一个总产量大于总销量的产销不平衡运输问题。
根据题目已知可以得出以下图论：二、模型的建立与基本假设2.1模型的建立假设某物品有m 个产地 A 1、A 2、…、 A m ，各产地的产量是a 1、a 2、…、a m ；有n 个销地B 1、B 2、…、B n ，各销售地销量分别为b 1、b 2、…、b n ；假定从产地A i （i=1,2，…，m ）向销售地B j （j=1，2，…，n ）运价单位物品的运价是c ij ，问这样调运这些物品才能使运费最少？设 x ij 为从产地Ai 运往销地Bj 的运输量，若各产地产量之和大于各销地销量之和，即有：∑∑==>nj jm i i ba 11则得到下列产销平衡运输量问题的模型：⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧≥===≤=∑∑∑∑=+===0,...,2,1,,...,2,1,min 11111ij j mi ij i n j ij m i nj ijij x n j b x m i a x x c z 其中，约束条件右侧常数ai 和bj ，约束条件最多有m+n-1个有效，即最多有m+n-1个基可行解。
为了能使用表上作业法，可增加一个假想的销地虚销地B n+1而由产地Ai （i=1,2，…，m ）调运到这个假想销地的物品数量的销量X i ，n+1（相当于松弛变量），实际上就地储存在Ai 。
因为就地储存没有运输，故单价为C i ，n+1=0，（i=1,2，…，m ）令假想销地的销量为：∑∑==+-=nj j m i i n b a b 111从而数学模型：⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧≥=====∑∑∑∑=+===0,...,2,1,,...,2,1,min 11111ij j mi ij i n j ij m i nj ijij x n j b x m i a x x c z2.2模型的基本假设针对该运输问题，为了方便计算，可以设新普（A1），隆宇（A2）和恒华（A3）分别销往北京（B1）、天津（B2）、广东（B3）和上海（B4）四个城市销售量为x11、x12、x13、x14、x21、x22、x23、x24、x31、x32、x33、x34。
建立以下模型： 表：1-2B1 B2 B3 B4 产量 A1 6 2 6 7 30 A2 4 9 5 3 25 A3 8 8 1 5 21 销量15172212-目标（The objective ） 最少费用：34333231242322312114131211415x x 8x 8x 3x 5x 9x 4x 7x 6x 2x 6x z Min +++++++++++==∑∑==i j j i j i x c约束条件：供应限制（The supply constrains ）⎪⎩⎪⎨⎧≤+++≤+++≤+++21x x x x 25x x x x 30x x x x 343332312423222114131211 指标约束（The damand constrains ）⎪⎪⎩⎪⎪⎨⎧=++=++=++=++12x x x 22x x x 17x x x 15x x x 342414332312322212312111 三、模型的定义符号说明与表上作业法定义符号说明：A1、A2、A3分别代表新普，隆宇和恒华生产商；B1、B2、B3、B4分别代表北京，天津，广东，上海销售地。
x11、x12、x13、x14、x21、x22、x23、x24、x31、x32、x33、x34为新普、隆宇和恒华分别销往北京、天津、广东和上海四个城市销售量。
Cij 为从产地Ai （i=1,2，…，m ）向销售地Bj （j=1，2，…，n ）运价单位物品的运价， x ij 为从产地Ai （i=1,2，…，m ）运往销地Bj （j=1，2，…，n ）的运输量。
Z 即为整个运输过程中涉及的运输费用。
Min z 则为该运输问题中的最小费用。
表上作业法（最小元素法）：最小元素法：是找出运价表中最小的元素，然后在运量表内对应的格填入允许取得的最大数值，若某行或者某列的产量或者销量已得到满足，则把运价表中该运价所在行或者列划去；找出未划去的运价中的最小数值，按此办法依次进行下去，直至得到一个基本可行解的方法。
表上作业法：是求解运输问题的一种简便而有效的方法，求解过程在运输表上进行行，这是一种迭代求解法，迭代步骤为： 步骤一：按某种规则找出一个初始基可行解。
步骤二：对进行解作最有判断，即求个非基变量的检验数，判别是否达到最优解。
如果已经是最优解，则停止计算；如果不是最优解，则进行下一步骤。
步骤三：在表上对初始方案进行改进，找出新的基可行解，再按照步骤二进行判别，直至找出最优解。
表上作业法具体求解如下： 表：1-3：步骤一：从表1-2中找出最小运价为1，故首先考虑此项，由于A3产地产量小于B3销量（21<22），故在表1-3的（A3，B3）交叉格填上21，由于A3产地产量已经饱和，故划去表1-3中的A3行得表1-4。
表：1-4步骤二：从表1-4中找出最小运价为2，故首先考虑此项，由于A1产地产量大于B2销量（30>17），故在表1-3的（A1，B2）交叉格填上17，由于B2销量已经饱和，故划去表1-4中的B2列得表1-5。
表：1-5-12 22 17 15销量21 5 0 1 21 8 0 8 0A325 312 5 0 9 0 4 13A2172307 06 12 6 A1产量 B4 B3 B2 B1步骤三：从表1-5中找出最小运价为3，故首先考虑此项，由于A2产地产量大于B4销量（25>12），故在表1-3的（A2，B4）交叉格填上12，由于B4销量已经饱和，故划去表1-5中的B2列得表1-6。
表：1-6步骤四：从表1-6中找出最小运价为4，故首先考虑此项，由于A2产地剩余产量小于B1销量（25-12=13<15），故在表1-3的（A2，B1）交叉格填上12，由于A2产地产量已经饱和，故划去表1-6中的A2行得表1-7。
表：1-7步骤五：从表1-7中找出最小运价都是6，故随机选择一项优先考虑此处选择（A1，B1），由于A1产地剩余产量大于B1剩余销量（30-17=13>15-13=2），故在表1-3的（A1，B1）交叉格填上2，由于B1销量已经饱和，故划去表1-5中的B2列。
步骤六：由于B3销地为达到饱和，故在（A1，B3）交叉格填上1，然后在其它空格位置统一填上0。
经以上步骤得到一个总产量大于总销量，且销量全部满足的调配方案。
经过计算，空格的检验数均大于零，最优方案为：x 2112131172ij 332412131211=======，其它，，，，，x x x x x x 最小费用为：∑∑===⨯+⨯+⨯+⨯+⨯+⨯==31411611213124136121762i jj ij ix cz四、问题求解4.1、lingo求解模型：LINGO模型：model:sets:origin/1..3/:a;sale/1..4/:b;routes(origin,sale):c,x;endsetsdata:a=30,25,21;b=15,17,22,12;c=6,2,6,7,4,9,5,3,8,8,1,5;enddata[OBJ]min=@sum(routes:c*x);@for(origin(i):[SUP]@sum(sale(j):x(i,j))<=a(i));@for(sale(j):[DEM]@sum(origin(i):x(i,j))=b(j));end4.2、lingo结果：Global optimal solution found.Objective value: 161.0000 Infeasibilities: 0.000000 Total solver iterations: 6Variable Value Reduced CostX( 1, 1) 2.000000 0.000000X( 1, 2) 17.00000 0.000000X( 1, 3) 1.000000 0.000000X( 1, 4) 0.000000 2.000000X( 2, 1) 13.00000 0.000000X( 2, 2) 0.000000 9.000000X( 2, 3) 0.000000 1.000000X( 2, 4) 12.00000 0.000000X( 3, 1) 0.000000 7.000000X( 3, 2) 0.000000 11.00000X( 3, 3) 21.00000 0.000000X( 3, 4) 0.000000 5.000000Row Slack or Surplus Dual PriceOBJ 161.0000 -1.000000SUP( 1) 10.00000 0.000000SUP( 2) 0.000000 2.000000SUP( 3) 0.000000 5.000000DEM( 1) 0.000000 -6.000000DEM( 2) 0.000000 -2.000000DEM( 3) 0.000000 -6.000000DEM( 4) 0.000000 -5.000000五、模型分析与改进从计算结果可以得出，新普（A1）分别销往北京（B1）、天津（B2）、广东（B3）和上海（B4）四个城市销售量为分别为2万台，17万台，1万台，0万台，剩余10万台；隆宇（A2）分别销往北京（B1）、天津（B2）、广东（B3）和上海（B4）四个城市销售量为别为13万台，0万台，0万台，12万台，剩余0万台；恒华湖北师范大学数学与统计学院运筹学论文（A3）分别销往北京（B1）、天津（B2）、广东（B3）和上海（B4）四个城市销售量为分别为0万台，0万台，21万台，0万台，剩余0万台；总费用为161个单位。
