前面建立的UMpHMP模型是NP-hard問(wèn)題,目前沒(méi)有有效的算法。為了減少計(jì)算時(shí)間,Ernst和Krishnamoorthy(1998a,1998b)為樞紐網(wǎng)絡(luò)建立了三下標(biāo)的 數(shù)學(xué)模型,大大減少了變量和約束的個(gè)數(shù),提高了求解的效率。 這個(gè)模型不采用流量的比例作為流變量,而直接用流量為流變量,并對(duì)匯運(yùn)、 轉(zhuǎn)運(yùn)和分運(yùn)分別設(shè)置不同的變量。
令Z4為OD流的匯運(yùn)流變量,即從始發(fā)地機(jī) 場(chǎng)i到樞紐機(jī)場(chǎng)k的流量,Ya是從輪輻機(jī)場(chǎng)i運(yùn)出的轉(zhuǎn)運(yùn)流量,X是O-D對(duì)(i,j) 從樞紐機(jī)場(chǎng)l分運(yùn)到目的地j的流量,再令O是從始發(fā)地機(jī)場(chǎng)i運(yùn)出的總流量, N={1,2,…,n}是機(jī)場(chǎng)集合,則三下標(biāo)樞紐網(wǎng)絡(luò)優(yōu)化模型如下; minC[2xCaZx+22aCuYu+22oC,x, O.表示起始于機(jī)場(chǎng)i的流量,因此二W,=0;式(3-15)表示共有p個(gè)樞 紐;式(3-16)~式(3-18)是流量平衡方程,保證所有的O-D流全部由起始機(jī)場(chǎng)到目的地機(jī)場(chǎng),其中式(3-16)表示從始發(fā)機(jī)場(chǎng)i出發(fā),經(jīng)過(guò)分配連接到各樞紐機(jī)場(chǎng) 的匯運(yùn)流量之和必須等于O,式(3-17)表示O-D對(duì)(i,j)經(jīng)過(guò)各樞組機(jī)場(chǎng)l中轉(zhuǎn) 后分運(yùn)到目的地機(jī)場(chǎng)j的流量之和必須等于W,,式(3-18)表示從始發(fā)地機(jī)場(chǎng)i運(yùn) 出的流量,在樞紐機(jī)場(chǎng)k中轉(zhuǎn)的運(yùn)進(jìn)流量必須等于運(yùn)出流量;式(3-19)表示從始發(fā) 機(jī)場(chǎng)i出發(fā),經(jīng)航節(jié)i→k匯運(yùn)時(shí),k一定是樞紐城市,式(3-20)表示如果O-D對(duì)(i、 j)經(jīng)過(guò)機(jī)場(chǎng)l中轉(zhuǎn)分運(yùn)到目的地機(jī)場(chǎng)j,l一定是樞紐;式(3-21)表示流變量是非 負(fù)變量,式(3-22)要求y;是0-1型變量。 這個(gè)模型盡管采用了三下標(biāo)變量,但任一O-D對(duì)仍然至多經(jīng)過(guò)匯運(yùn)、轉(zhuǎn)運(yùn)和 分運(yùn)等三個(gè)航節(jié)完成運(yùn)輸任務(wù),因此最多兩次中轉(zhuǎn)。
變量數(shù)從四下標(biāo)的O(n1)個(gè) 減少為O(n3)個(gè)。如果n=100,那么流變量數(shù)從1億個(gè)減少到100萬(wàn)個(gè)。約束條 件數(shù)從四下標(biāo)模型的(2n3+n2+1)減少到三下標(biāo)模型的(n3+3n2+n+1),當(dāng)n= 100時(shí),四下標(biāo)模型大約有201萬(wàn)個(gè),三下標(biāo)模型大約有103萬(wàn)個(gè)??梢?,三下標(biāo) 模型的規(guī)模確實(shí)下降不少。 例3-4在例3-3同樣的條件下,采用三下標(biāo)模型對(duì)樞紐航線網(wǎng)絡(luò)進(jìn)行優(yōu)化 設(shè)計(jì)。構(gòu)建三下標(biāo)模型如下: 同樣應(yīng)用ILOG/CPLEX求解上述模型,得到的結(jié)果為:y2=y;=1,即以北京 和廣州作為樞紐;最優(yōu)的總成本是950364元,這與四下標(biāo)模型得到的結(jié)果相同;流 變量的解如下:Z2=106,Z2=354,Z2=97,Za=46,Z46=22,Z4s=231,Z2=63, 模型中一共有331個(gè)約束條件和474個(gè)變量。與四下標(biāo)模型比較可以發(fā)現(xiàn), 無(wú)論是約束條件數(shù)、變量數(shù)還是運(yùn)算時(shí)間,三下標(biāo)模型都小于四下標(biāo)模型。當(dāng)問(wèn)題 的規(guī)模進(jìn)一步擴(kuò)大后,這一點(diǎn)將會(huì)更加明顯。