一級中轉場對應區(qū)域劃分

1、一級中轉場對應區(qū)域劃分由于三個一級中轉場地理位置的特殊性，到達華東地區(qū)各個二級三級中轉場的距離不同, 且各個二級三級中轉場的快件的需求量不同，若按三個省行政區(qū)域完全劃分，則不能使得費 用最低，達到最優(yōu)。所以，我們通過數學模型和lingo求解，對一級中轉場對應區(qū)域進行劃分，來達到優(yōu)化 的目的。在建模前，我們假設由一級中轉場到各個二三級中轉場的單位距離單位重量包裹的運費 是一個常數，用C表示。并為了簡化模型，只考慮由二三級中轉場到一級中轉場的收快件 業(yè)務量，由此可知，這是一個多產地多銷地的運輸問題。但這個問題特殊處在于，每個產地 即二三級中轉場，只對應于一個銷地即一級中轉場。所以我們要增加一個權重

2、變量x,它是 一個01變量，取0表示該產地不對應于該銷地，取1表示該產地對應于該銷地。建立數學模型如下：2孕SjyijXijj =bL j=l t(4)(5)yij 2 0, i = 1,2,j = 1,2,Xjj e0,1 i = 1,2,j = 1,2,其中，C=每萬噸包裹每千米所需平均運費，單位：元/(萬噸*千米)0=二三級中轉場1的年收件重量，單位：萬噸tj級中轉場j的年吞吐量， 單位：萬噸，耳=由二三級中轉場1到一級中轉場j的距離，單位：千米Xjj= 權重，如二三級中轉場1屬于一級中轉場j劃分區(qū)域，則否則，x=0)弓=由二三級中轉場1到一級中轉場j的運量，單位：萬噸目標函數最優(yōu)表示的

3、每年的運費最少，年運費=單位重量單位距離的平均運費*兩地 距離*二三級中轉場的年運量*權重；式是指對每個二三級中轉場來說，權重之和為一；式是指對每個二三級中轉場來說，總的運量之和為該產地的年收件重量；式是指對每個一級中轉場來說，運往該地的總運量之和應小于等于該地的年吞吐式指所有運量都是正數；式指如二三級中轉場1屬于一級中轉場j劃分區(qū)域，則否則，x=0,即x是 二進制數。在上述解上述模型之前，我們要先計算各個二三級中轉場到各一級中轉場的距離。先 通過網上的地圖找出各地區(qū)經緯度坐標，如下表：地區(qū)緯度經度上海市31。1349”121。28切“南京市323,37H118。4749“無錫市31。3358

4、“12018,ir,徐州市34。1217”I。；”常州市31。4836”11958,28,1蘇州市31。1756”120357H南通市31。5848”12053,40u連云港34。3548”119。1317”淮安市3336,37H1190,55M鹽城市33o20,51H1209,49H揚州市3223,39M1192447”鎮(zhèn)江市32。1157”119O27I26U泰州市3227,21H119。55。23“宿遷市3357,43M118。1652“杭州市30。1827”1209* 18M寧波市2952,6H121O32,38,1溫州市300,0H120o20,24,1嘉興市294* 12H119。

5、23*24“紹興市2959,45MUOWIO”金華市29。445”11938,5r,衢州市2856,9H118O52,28,舟山市29059-8-12212,24“臺州市2839,23M121025,15,1麗水市28 283”119。5522”再用一個根據經緯度計算距離的公式，如下圖所示：d=lll. 12cosl/sin0AsinPB + cosOAcosOBcos(A BA A)(其中A點經度，緯度分別為入AW A, B點的經度、緯度分別為人B和d 為距離)。再根據上訴公式用excel計算距離。先將上述公式改為excel公式如下：=6371000*ACOS (COS (RADIANS (

6、SUM (1*LEFT (A2, FIND C A2)-l) , MID(A2, FIND C ，A2)+l, FIND(A2) -FIND( ，A2)-l)/60, RIGHT (LEFT (A2, LEN(A2)-1), LEN(A2)-FIND( A2)-l)/3600)-SUN(l*LEFT(Al, FI NDC A1)-1),MID(A1,FINDC =A1)+L，FIND( = A1)-FIND( , AD-D/60, RIGHT (LEFT (Al, LEN(Al) -1),LEN(A1)-FIND( Al)-1)/3600)*COS(RADIANS(SUM(1*LEFT(B2,

7、FINDC0 B2)-l), MID(B2,FINDC 0 L B2)+l, FINDC，B2)-FINDC I B2)-l)/60, RIGHT (LEFT(B2, LEN(B2)-l), LEN (B2)-FIND(。B2)-l)/ 3600)-SUM(1*LEFT(Bl,FINDC FINDC ,Bl)+1, FINDC* , B1)-FINDC B1)-1)/6O,RIGHT (LEFT (Bl, LEN (Bl)-1), LEN (Bl)-FIND(#, I Bl)-1)/3600)其中：Al二第一個地點的經度，A2二第二個地點的經度，B1二第一個地點的緯度，B2二第二 個地點的緯度

8、;A格式如例“33。3637”所示，B格式如例“119。2726”所示。計算 得到距離如下表：上海無錫杭州上海市0. 0311. 6181. 1南京市311. 60. 0249. 5無錫市135.4176.3144.6徐州市571. 1291. 8541. 1常州市178.7133.9172.0蘇州市99. 1216. 1123.6南通市105.3233.4206. 7連云港450.3286. 0491. 7淮安市380. 4174. 1392.0鹽城市276.8208. 8341.7揚州市263. 177.9249. 7鎮(zhèn)江市248.775.0227.7泰州市219.8132.7244.0宿

9、遷市467.6219. 3460. 2杭州市181. 1249. 50. 0寧波市151. 6390. 7160. 9溫州市186. 0286. 236. 8嘉興市333.7339. 0158. 6紹興市169. 1303.757. 1金華市313.7344.7144.4衢州市385.5347.5206. 0舟山市160. 6443. 7230. 4臺州市286. 3477. 8228. 4麗水市352.3418. 6202.5通過查找資料，我們得出各二三級中轉場的年收件重量如卜表所示:地區(qū)收件量 （萬噸）上海市1.50南京市0.62無錫市0.47徐州市0.28常州市0.27蘇州市0.54南通

10、市0.38連云港0.10淮安市0.13鹽城市0.31揚州市0.30鎮(zhèn)江市0.19泰州市0.29宿遷市0.08杭州市0.56寧波市0.43溫州市0.35嘉興市0.16紹興市0.20金華市0.19衢州市0.08舟山市0.05臺州市0.22麗水市0.08各一級中轉場的年吞吐量如下表所示:上海無錫杭州年吞吐量（萬噸）2.42.312由上述數據，運用ling。求解該數學模型，程序如下:MODEL：sets:M/1.24/:b;N/3,1,15/:t;LINK (M,N) : x,y, s;endsetsdata:c=2;s =0.0311. 6181. 1311. 60. 0249.5135.4176.

11、3144.6571. 1291.8541. 1178. 7133.9172.099. 1216. 1123.6105.3233.4206.7450. 3286.0491.7380.4174. 1392.0276. 8208.8341.7263. 177.9249.7248. 775.0227.7219. 8132.7244.0467.6219.3460.2181. 1249. 50.0151. 6390.7160.9186. 0286.236.8333.7339.0158.6169. 1303.757. 1313.7344.7144.4385.5347.5206.0160. 6443.7230

12、.4286. 3477.8228.4352.3418.6202.5b=l.50,0 62,0 47,0 28,0.27,0.54,0.38,0.10,0.13,0 31,0.30,0.19,0.29,0.08,0.56,043,0.35,0.16,0.20, 0 19,0 08,0.05,0.22,0.08;t=2.3,2.4, 12;endclataobjmin=sum(link:c*x*y*s);for (M (i) :sum (N (j): x (iz j) =1;);for (M (i) :sum (N (j): x (iz j)*y (i, j ) ) =b (i););for (N (j) :sum (M (i): y (iz j) t (j););for(link(iz j):bin(x(iz j););END求解得到狀態(tài)窗11如下：由輸出報告中的數值可得到運費最小的各一級中轉場對應二三級中轉場情況如卜.表:上海無錫杭州上海市100南京市010無錫市001徐州市010常州市001蘇州市001南通市100連云港010淮安市010鹽城市010揚州市010鎮(zhèn)江市010泰州市010宿遷市010杭