一级中转场对应区域划分.docVIP

一级中转场对应区域划分 由于三个一级中转场地理位置的特殊性，到达华东地区各个二级三级中转场的距离不同，且各个二级三级中转场的快件的需求量不同，若按三个省行政区域完全划分，则不能使得费用最低，达到最优。 所以，我们通过数学模型和lingo求解，对一级中转场对应区域进行划分，来达到优化的目的。 在建模前，我们假设由一级中转场到各个二三级中转场的单位距离单位重量包裹的运费是一个常数，用C表示。并为了简化模型，只考虑由二三级中转场到一级中转场的收快件业务量，由此可知，这是一个多产地多销地的运输问题。但这个问题特殊处在于，每个产地即二三级中转场，只对应于一个销地即一级中转场。所以我们要增加一个权重变量x，它是一个01变量，取0表示该产地不对应于该销地，取1表示该产地对应于该销地。 建立数学模型如下:min i=1mj=1nCSijyijxij s.t. j=1nxij=1, i=1,2,…,m, (1) j=1nyij=bi, i=1,2,…,m , (2) i=1myij≤tj, j=1,2,…,n, (3) yij≥0, i=1,2,…,m, j=1,2,…,n, (4) xij∈0,1 i=1,2,…,m, j=1,2,…,n, (5) 其中，C=每万吨包裹每千米所需平均运费，单位：元/(万吨*千米) bi=二三级中转场i的年收件重量， 单位：万吨 tj=一级中转场j的年吞吐量， 单位：万吨 Sij=由二三级中转场i到一级中转场j的距离， 单位：千米 xij=权重，如二三级中转场i属于一级中转场j划分区域，则x=1；否则，x=0 yij=由二三级中转场i到一级中转场j的运量， 单位：万吨 目标函数最优表示的每年的运费最少，年运费=单位重量单位距离的平均运费*两地距离*二三级中转场的年运量*权重； (1)式是指对每个二三级中转场来说，权重之和为一； (2)式是指对每个二三级中转场来说，总的运量之和为该产地的年收件重量； (3)式是指对每个一级中转场来说，运往该地的总运量之和应小于等于该地的年吞吐量； (4)式指所有运量都是正数； (5)式指如二三级中转场i属于一级中转场j划分区域，则x=1；否则，x=0，即x是二进制数。 在上述解上述模型之前，我们要先计算各个二三级中转场到各一级中转场的距离。先通过网上的地图找出各地区经纬度坐标，如下表： 地区纬度经度上海市31°1349121°2825南京市32°337118°4749无锡市31°3358120°1811徐州市34°1217117°173常州市31°4836119°5828苏州市31°1756120°357南通市31°5848120°5340连云港34°3548119°1317淮安市33°3637119°055盐城市33°2051120°949扬州市32°2339119°2447镇江市32°1157119°2726泰州市32°2721119°5523宿迁市33°5743118°1632杭州市30°1627120°918宁波市29°526121°3238温州市30°00120°2024嘉兴市29°412119°2324绍兴市29°5945120°3510金华市29°445119°3851衢州市28°569118°5228舟山市29°598122°1224台州市28°3923121°2515丽水市28°283119°5522再用一个根据经纬度计算距离的公式，如下图所示： 再根据上诉公式用excel计算距离。先将上述公式改为excel公式如下： =6371000*ACOS(COS(RADIANS(SUM(