【導(dǎo)讀】第一章安吉汽車物流發(fā)展現(xiàn)狀分析···································-6-. 安吉汽車物流發(fā)展概況···········································-6-. 社會與文化環(huán)境分析·····························-11-

　　

【正文】 i的運輸任務(wù)由配送中心 u和配送中心 v協(xié)同完成 ） ，否則為 0。 若 ipqZ 為 1，則客戶 i由車輛類型 p第 q輛車完成，否則為 0。 0pqt 為車輛類型 p第 q輛車從車場出發(fā)的時間。 ipqA 為車輛類型 p第 q輛車在完成客戶 i前服務(wù)客戶的集合。 ijd 為客 戶 （ 配送中心 ） i和客戶 (配送中心 )j之間的距離。 ~ijv 為客 戶 ( 配送中心 )i和客戶 ( 配送中心 )j之間的模糊行駛速度。 pr 為車輛類型 p單位里程成本。 “安吉杯” 第四屆全國物流設(shè)計大賽 夢飛 amp?！氨税丁? 24 數(shù)學(xué)模型可描述為： ? ? ? ?,39。39。1 1 1 1 1 1 11111 1 1 121m i nm a x 0 m a x 0m a x 0m a xPPQP N N M M NP ij ijp q uv uv pq iuvP q i M j M u v iNNi i iiiiiQP N Nx j ijp qPP q i M j MNuv uv pq iuvviZ r d x d x yf e t t f X t l tA d x DA d x y? ? ? ? ? ? ? ? ? ? ???? ? ? ? ? ? ? ????????? ? ???????????? ? ? ?????????? ? ? ? ? ? ???? ? ? ?? 39。39。39。1 1 1 131 1 141 1 111110,0m a x 0m a x 0. : 1 , p, q tPPPQP M MPP q uQPNi ipqpP q iQPNi iuv ipqpP q iNNijp q jip qjjNNx jpq jip q ipqjji i pq uv pqu v ADA g z vA g y z vs t x xx x ze t t x? ? ? ?? ? ?? ? ?????????????????????????????? ? ???? ? ?? ? ? ?? ? ?? ? ?????? ?39。/ 1 ipquv uv v ipqid u ut l t z? ? ? ?? 在上述模型中，式 (1)是問題的目標(biāo)函數(shù)，前三項分別是運輸費用 (包括各配送中心自身的配送費用和配送中心之間的協(xié)同運輸費用 ) 、等待費用和延時費用，后四項是某 車型車輛最大里程及容量的懲罰函數(shù) ( Ai為懲罰因子，可取一個足夠大的 數(shù) )。 式 (2)和 (3)表示一個客戶僅能被服務(wù)一次以及一個客戶僅能被一輛車服務(wù) ； 式 ( 4) 表示車輛從客戶 j～ i的時間。 根據(jù)文獻(xiàn)［ 13， 14］可知，將帶模糊參數(shù)的數(shù)學(xué)規(guī)劃模型首先利用三角模糊數(shù)轉(zhuǎn)換成相應(yīng)的機(jī)會約束規(guī)劃模型，再由可能性理論證明進(jìn)一步清晰等價在某置信水平下對應(yīng)的確定規(guī)劃模型 。 本文將帶模糊參數(shù)的目標(biāo)函數(shù) (1)及約束條件 (4) 的 CVRP轉(zhuǎn)換為置信水平 α和 β下對應(yīng)的 CVRP確定規(guī)劃模型 如下頁 。 （ 1） （ 2） （ 3） （ 4） “安吉杯” 第四屆全國物流設(shè)計大賽 夢飛 amp?！氨税丁? 25 ? ?39。~1 1 1 1 1 1 10121 ,12,m in :. . :m a x / ( 1 ) . .( 1 ) . . 0)m a x / 1 .PipqipqQp N N M M Np x j ijpq uv uv pq iuvp q i M j M u v iNi pq uv un uvi u v Avvii uv uu v AZs t r d x d x yf X t d v vut ut ltf e t d v?????? ? ? ? ? ? ? ? ? ? ?? ??????????????????????? ? ? ???????? ? ? ???? ? ? ? ? ? ????39。39。39。12101211 1 1 121 1 1 1 131 1 14.( 1 ) . . , 0)m a x 0m a x 0m a x 0m a xPPPNn uviv v pqQP N Nx j ijpqPP q i M j MQP M M Nuv uv pq iuvPP q u v iQPNi ipqpP q ivut ut tA d x DA d x y DA g z vA????? ? ? ? ? ? ? ?? ? ? ? ?? ? ?????????? ? ? ??????????????????????????? ? ? ?? ? ? ? ?? ? ?? ?39。~1 1 1012,120( / ( 1 )( 1 ) , 1PipqQPNi iuv ipqpP q ii i uv pq uv uv uvpqu v Av v i ipqg y z v ze t t x d v vut ut lt zt????? ? ??????????? ? ? ? ? ? ? ?? ? ? ? ? ? ?? ? ?? 其中 : ( 1jut ， 2jut ， 3jut ) 和 ( 1ijv ， 2ijv ， 3ijv ) 分別是 iut~ 和 ijv~ 的三角模糊數(shù) ， 其他約束條件與 CVRP相同 。 （二） 模型求解 算法 根據(jù)上述模型粒子群的運行參數(shù)設(shè)置如下 : 最大迭代 次數(shù) max=1000，粒子群大小 P=40，以及 η1=， η2 = η3 =。針對一類 CVRP的粒子編碼方式（ 5） （ 6） （ 7） “安吉杯” 第四屆全國物流設(shè)計大賽 夢飛 amp?！氨税丁? 26 表示問題的可行解，用目標(biāo)函數(shù)直接評價該粒子的適應(yīng)值，具體的 SADPSO求解 CVRP的流程如下 圖 22所示 。 開始 設(shè)置 PSO 參數(shù)，如 max、 M、 1? 、 2? 、和 3? 等 初始化 PSO 第 iter=0 代粒子 iter=iter+1 和 m=0 M=m+1 根據(jù)位置更新公式，第 m個粒子計算自己的最新位置，對應(yīng)ＣＶＲＰ解析粒子編碼方式計算目標(biāo)值，并更新自己最佳“運輸經(jīng)驗” mM? 更新群體最佳粒子“運輸經(jīng)驗” 計算 ()x? ( ) 0x?? 混沌擾動算子 maxiter iter? 結(jié)束 否 是 是 否 否 是 圖 22 SADPSO 求 解 CVRP 的流程 “安吉杯” 第四屆全國物流設(shè)計大賽 夢飛 amp?！氨税丁? 27 簡單算例驗證 假設(shè)安吉汽車物流 兩個配送中心利用兩車型車輛 （ 1 5v? ， 2v = 6， 1D = 2D =500， 1r = 2. 9和 2r = 2. 5） 為屬于它們自己的客戶 ( if = 100， f =15)配送貨物，客戶的基本信息及客戶之間行駛速度如表 22和 表 23所示。其中置信水平 α = 及β = 。根據(jù) u， v且 iuvy =0或 u， v且 iuvy = 1配送方案可分為兩類 ： 方案 1是每個配送中心獨自承擔(dān)自己的任務(wù)配送，追求每個配送中心行使路線最優(yōu) ； 方案 2 是兩配送中心之間協(xié)同為客戶運輸，制定的配送方案目標(biāo)為它們總的運輸費用最少 。 表 22 客戶服務(wù)情況一覽表 客戶 貨運量 坐標(biāo) 客戶服務(wù)時間 時間窗 服務(wù)公司 1 （ 111,37） （ 64,65,70） 8:00— 13:00 A 2 （ 99,126） （ 79,84,94） 8:00— 11:30 A 3 （ 48,117） （ 81,87,96） 8:00— 17:00 A 4 （ 97,45） （ 74,75,81） 8:00— 11:30 A 5 （ 94,127） （ 91,99,102） 8:00— 17:00 A 6 （ 101,59） （ 66,70,76） 8:00— 12:30 A 7 （ 107,129） （ 90,92,98） 10:00— 17:00 B 8 （ 44,160） （ 50,53,55） 8:00— 16:00 B 9 1,1 （ 36,13） （ 74,80,85） 8:00— 17:00 B 10 （ 39,41） （ 68,73,77） 13:00— 15:00 B 11 （ 63,27） （ 65,69,83） 14:00— 17:00 B 12 （ 55,139） （ 72,77,80） 8： 00— 17:00 B “安吉杯” 第四屆全國物流設(shè)計大賽 夢飛 amp?！氨税丁? 28 表 23 客戶之間車輛模糊行駛速度 變 行駛速度 邊 行駛速度 （ A,2） （ 55,57,76） （ B,9） （ 41,59,69） （ 2,7） （ 53,56,67） （ 9,10） （ 44,75,76） （ 7,5） （ 43,44,55） （ 10,11） （ 45,62,77） （ 5,A） （ 56,66,75） （ 11,B） （ 52,55,76） （ B,4） （ 49,72,76） （ A,12） （ 58,69,80） （ 4,1） （ 45,61,75） （ 12,8） （ 54,72,73） （ 1,6） （ 53,54,65） （ 8,3） （ 58,59,61） （ 6,A） （ 45,46,58） （ 3,A） （ 54,59,60） 通過上述程序運行 10次求得問題的方案 1和 2的尋優(yōu)結(jié)果如表 24所示。從表24可知，程序以概率 100%找到方案 1的最優(yōu)解，而分別以 80%和 20%概率找到方案 2的最優(yōu)解及次優(yōu)解 ， 方案 1的平均配送費用和運輸里程比方案 2超出 %和%。上述方案 2的總費用和里程包括協(xié)同部分費用和里程均為 421和 。方案 1和 2的最佳配送方案詳細(xì)配送情況如表 下頁 25所示 。 表 24 10次程序運行結(jié)果 序號 方案 1 方案 2 總費用 總里程 總費用 總里程 1 2 3 4 5 6 7 8 9 10 平均 “安吉杯” 第四屆全國物流設(shè)計大賽 夢飛 amp。“彼岸” 29 表 25 方案 1和方案 2配送計劃 方案 車輛類型 行駛路線及時間表 總里程 費用 A A25A/8:00— 8:25— 9:46— 11:51 B A6413A/8:00— 9:23— 10:46— 12:17— 15:23— 17:18 1 A B91011B/8:00— 8:44— 10:31— 12:16— 13:44 B B1287B /8:00— 9:42— 11:24— 13:45— 17:01 A B416A/8:00— 8:14— 9:45— 11:17