1、-. z. . . . . 資料. . .2012年數(shù)學(xué)建模培訓(xùn)第二次測(cè)試論文 題目 運(yùn)輸優(yōu)化模型 *馬鵬 系院 數(shù)學(xué)系 專業(yè) 信息與計(jì)算科學(xué)、應(yīng)用數(shù)學(xué) 2012 年 8 月27日運(yùn)輸優(yōu)化模型 摘要在社會(huì)的經(jīng)濟(jì)生產(chǎn)活動(dòng)中，產(chǎn)地廠家與客戶都會(huì)想方設(shè)法合理調(diào)撥資源、降低運(yùn)輸費(fèi)用，實(shí)現(xiàn)利益最大化，完成資源優(yōu)化配置。本文在運(yùn)輸費(fèi)單價(jià)恒定，各產(chǎn)地發(fā)量一定，各客戶的需求量也一定的條件下，努力解決多個(gè)特定目標(biāo)實(shí)現(xiàn)問(wèn)題。力求最優(yōu)的運(yùn)輸方案。在確定問(wèn)題為不平衡的運(yùn)輸問(wèn)題時(shí)，先虛設(shè)一個(gè)產(chǎn)地，將問(wèn)題裝華為平衡運(yùn)輸問(wèn)題，將問(wèn)題轉(zhuǎn)化為目標(biāo)規(guī)劃問(wèn)題，按照目標(biāo)規(guī)劃問(wèn)題的建模思想逐步建立模型。本文的主要特點(diǎn)在于，將不平衡的線

2、性規(guī)劃問(wèn)題合理地轉(zhuǎn)化為目標(biāo)規(guī)劃問(wèn)題，在求解時(shí)充分利用LINGO軟件求解。關(guān)鍵詞： lingo 目標(biāo)規(guī)劃 線性規(guī)劃 運(yùn)輸優(yōu)化問(wèn)題 運(yùn)費(fèi)最少問(wèn)題重述運(yùn)輸功能是整個(gè)現(xiàn)代物流七大根本功能之一，占有很重要的地位，運(yùn)輸本錢在整個(gè)物流系統(tǒng)中所占的比重也很大，運(yùn)輸本錢的有效控制對(duì)物流總本錢的節(jié)約具有舉足輕重的作用。通過(guò)物流流程的改善能降低物流本錢，能給企業(yè)帶來(lái)難以預(yù)料的效益，影響運(yùn)輸本錢的因素是多樣化、綜合性的，這就要求對(duì)運(yùn)輸本錢的分析要采用系統(tǒng)的觀點(diǎn)，進(jìn)展綜合分析。由于影響物流運(yùn)輸本錢的因素很多，控制措施既涉及運(yùn)輸環(huán)節(jié)本身，也涉及供給鏈的整個(gè)物流流程。要想降低物流運(yùn)輸本錢，就必須運(yùn)用系統(tǒng)的觀點(diǎn)和方法，進(jìn)展綜

3、合分析，發(fā)現(xiàn)問(wèn)題，解決問(wèn)題，使物流運(yùn)輸活動(dòng)更加優(yōu)化、物流運(yùn)輸本錢更加合理化。本文把一種產(chǎn)品從產(chǎn)地一、二運(yùn)到客戶1、2、3處,產(chǎn)地的發(fā)量、客戶的收量及各產(chǎn)地到各客戶的運(yùn)輸單價(jià)。本文要解決問(wèn)題是：客戶1為重要部門，必須全部滿足需求量；滿足客戶2、3至少75%的的需求量；使總運(yùn)費(fèi)盡量少；從產(chǎn)地2到客戶1的運(yùn)量至少有1000個(gè)單位。問(wèn)題分析根據(jù)題目中所給出的條件知：有現(xiàn)成的兩個(gè)產(chǎn)地和需要產(chǎn)品的三個(gè)客戶。且兩個(gè)產(chǎn)地的產(chǎn)量不同，運(yùn)送到各個(gè)客戶的運(yùn)費(fèi)單價(jià)不同。三個(gè)客戶所需的貨物量不同。而三個(gè)客戶對(duì)兩個(gè)產(chǎn)地的總需求為2000+1500+5000=8500單位，而兩個(gè)產(chǎn)地總的發(fā)量為3000+4000=7000單

4、位,故需求量大于發(fā)量，屬于需求量和發(fā)量不平衡問(wèn)題。且提出四個(gè)不同的目標(biāo)。故使用目標(biāo)規(guī)劃實(shí)現(xiàn)建模。首先設(shè)置目標(biāo)約束的優(yōu)先級(jí)，建立目標(biāo)約束按目標(biāo)的優(yōu)先級(jí)，寫出相應(yīng)的目標(biāo)規(guī)劃模型 。再接著使用LINGO軟件實(shí)現(xiàn)模型的求解，并作出相應(yīng)結(jié)果的分析。三模型假設(shè)產(chǎn)品的運(yùn)輸過(guò)程不存在任何的導(dǎo)致產(chǎn)品發(fā)量和產(chǎn)品收量不相符的問(wèn)題。產(chǎn)品平安送到客戶處。即有：產(chǎn)品的發(fā)量就等于產(chǎn)品的收量。產(chǎn)品的運(yùn)輸單價(jià)始終恒定，不存在中途因?yàn)?種原因而導(dǎo)致產(chǎn)品的單價(jià)變化問(wèn)題。即運(yùn)費(fèi)只取決于所運(yùn)輸?shù)漠a(chǎn)品的數(shù)量。產(chǎn)地的生產(chǎn)量即發(fā)量有極限值，不可能超出本產(chǎn)地正常的生產(chǎn)*圍。客戶需求量在一定的*圍內(nèi)或或是特定的具體值。四符號(hào)說(shuō)明 基于題目及所要

5、建立的模型所要用到的變量及參數(shù)，作如下符號(hào)說(shuō)明：產(chǎn)地用表示，表示第產(chǎn)地i；表示其發(fā)量；客戶用其中j=1,2,3表示，表示客戶j;表示其需求量；用表示產(chǎn)地往客戶其中j=1,2,3處運(yùn)輸產(chǎn)品的單位費(fèi)用；用z表示總的運(yùn)輸費(fèi)用；用表示產(chǎn)地運(yùn)往客戶其中j=1,2,3處的物品數(shù)量；五模型建立由發(fā)量和需求量可知，發(fā)量小于需求量，故我們需要添加一個(gè)虛擬產(chǎn)地產(chǎn)地3，使各產(chǎn)地的總產(chǎn)量之和等于各客戶的需求量之和。使問(wèn)題為平衡的運(yùn)輸問(wèn)題。且令虛擬產(chǎn)地到各客戶的運(yùn)費(fèi)單價(jià)都為0，如表1所示：客戶1客戶2客戶3發(fā)量產(chǎn)地1104123000產(chǎn)地281034000產(chǎn)地30001500需求量200015005000 表1至此，基

6、于問(wèn)題的分析與假設(shè)，將問(wèn)題轉(zhuǎn)化為目標(biāo)規(guī)劃問(wèn)題。故分以下步驟進(jìn)展模型的建立。5.1設(shè)置目標(biāo)約束的優(yōu)先級(jí) P1：客戶1為重要部門，需求量必須全部滿足；P2：滿足其他兩個(gè)客戶至少75%的需要量；P3：使運(yùn)費(fèi)盡量少；P4：從產(chǎn)地2到客戶1的運(yùn)量至少有1000個(gè)單位。5.2建立目標(biāo)約束達(dá)不到客戶1的需求量超過(guò)客戶1的需求量達(dá)不到客戶2的需求量超過(guò)客戶2的需求量：超過(guò)客戶3的需求量：達(dá)不到33000的運(yùn)輸費(fèi)用超過(guò)33000的運(yùn)輸費(fèi)用 產(chǎn)地二達(dá)不到客戶1的需求量超過(guò)客戶1的需求量5.3求最少費(fèi)用 LINGO程序： model:sets: supply/1,2,3/:a; demand/1,2,3/:b; l

7、ink(supply,demand):c,*; endsets min=sum(link(i,j): c(i,j)*(i,j););for(demand(j):sum(supply(i): *(i,j)=b(j););for(supply(i):sum(demand(j): *(i,j)=a(i);); data: a=3000,4000,1500; b=2000,1500,5000; c=10,4,12 8,10,3 0,0,0; enddata EndLINGO求解結(jié)果：Global optimal solution found. Objective value: 33000.00 Inf

8、easibilities: 0.000000 Total solver iterations: 6 Variable Value Reduced Cost A( 1) 3000.000 0.000000 A( 2) 4000.000 0.000000 A( 3) 1500.000 0.000000 B( 1) 2000.000 0.000000 B( 2) 1500.000 0.000000 B( 3) 5000.000 0.000000 C( 1, 1) 10.00000 0.000000 C( 1, 2) 4.000000 0.000000 C( 1, 3) 12.00000 0.0000

9、00 C( 2, 1) 8.000000 0.000000 C( 2, 2) 10.00000 0.000000 C( 2, 3) 3.000000 0.000000 C( 3, 1) 0.000000 0.000000 C( 3, 2) 0.000000 0.000000 C( 3, 3) 0.000000 0.000000 *( 1, 1) 1500.000 0.000000 *( 1, 2) 1500.000 0.000000 *( 1, 3) 0.000000 2.000000 *( 2, 1) 0.000000 5.000000 *( 2, 2) 0.000000 13.00000

10、*( 2, 3) 4000.000 0.000000 *( 3, 1) 500.0000 0.000000 *( 3, 2) 0.000000 6.000000 *( 3, 3) 1000.000 0.000000 Row Slack or Surplus Dual Price 1 33000.00 -1.000000 2 0.000000 -10.00000 3 0.000000 -4.000000 4 0.000000 -10.00000 5 0.000000 0.000000 6 0.000000 7.000000 7 0.000000 10.00000我們?cè)趯?shù)據(jù)整理在一個(gè)表格中，如表

11、2所示：客戶1客戶2客戶3發(fā)量產(chǎn)地11500150003000產(chǎn)地20040004000產(chǎn)地3500010001500需求量200015005000 表2由上表可看出，最少的運(yùn)輸費(fèi)用為33000，但第一個(gè)目標(biāo)就不滿足，用戶1的需求的不到滿足。5.4按目標(biāo)的優(yōu)先級(jí)，寫出相應(yīng)的目標(biāo)規(guī)劃模型客戶1為重要部門，需求量必須全部滿足；則目標(biāo)可表示為：滿足其他兩個(gè)客戶至少75%的需要量；則目標(biāo)可表示為：從產(chǎn)地2到客戶1的運(yùn)量至少有1000個(gè)單位；則目標(biāo)可表示為：由最少費(fèi)用，可建立目標(biāo)約束為：故模型建立為：min z=模型求解 使用LINDO軟件將模型求解如下： LINGO程序： model:sets: Le

12、vel/1,2,3,4/:P,z,Goal; s_Con_Nun/1,2,3,4,5/:dplus,dminus; supply/1,2/:a; customer/1,2,3/:b; Routes(supply,customer):c,*;endsetsdata: p=,; Goal=,0; a=3000,4000; b=2000,1500,5000; c=14,4,12 8,10,3;enddatamin=sum(Level:P*z); z(1)=dminus(1) z(2)= dminus(2)+dminus(3); z(3)=dplus(4); z(4)=dminus(5);for(su

13、pply(i):sum(customer(j):*(i,j)=a(i);); *(1,1)+*(2,1)+dminus(1)-dplus(1)=2000;for(customer(j):sum(supply(i):*(i,2)+dminus(2)-dplus(2)=1500*0.75;sum(supply(i):*(i,3)+dminus(3)-dplus(3)=1500*0.75;sum(Routes:c*)+dminus(4)-dplus(4)=33000; *(2,1)+dminus(5)-dplus(5)=1000;for(Level(i)|i#lt#size(Level):bnd(0

14、,z(i),Goal(i););EndLINGO求解結(jié)果：No feasible solution found. Infeasibilities: 1500.000 Total solver iterations: 5 Variable Value Reduced Cost P( 1) 0.1000000+308 0.000000 P( 2) 0.1000000+308 0.000000 P( 3) 0.1000000+308 0.000000 P( 4) 0.1000000+308 0.000000 Z( 1) 0.000000 0.000000 Z( 2) 0.000000 0.00000

15、0 Z( 3) 13000.00 0.000000 Z( 4) 500.0000 0.000000 GOAL( 1) 0.1000000+308 0.000000 GOAL( 2) 0.1000000+308 0.000000 GOAL( 3) 0.1000000+308 0.000000 GOAL( 4) 0.000000 0.000000 DPLUS( 1) 0.000000 0.000000 DPLUS( 2) 375.0000 0.000000 DPLUS( 3) 3875.000 0.000000 DPLUS( 4) 13000.00 0.000000 DPLUS( 5) 0.000

16、000 0.1000000+308 DMINUS( 1) 0.000000 0.1000000+308 DMINUS( 2) 0.000000 0.1000000+308 DMINUS( 3) 0.000000 0.1000000+308 DMINUS( 4) 0.000000 0.1000000+308 DMINUS( 5) 500.0000 0.000000 A( 1) 3000.000 0.000000 A( 2) 4000.000 0.000000 B( 1) 2000.000 0.000000 B( 2) 1500.000 0.000000 B( 3) 5000.000 0.0000

17、00 C( 1, 1) 14.00000 0.000000 C( 1, 2) 4.000000 0.000000 C( 1, 3) 12.00000 0.000000 C( 2, 1) 8.000000 0.000000 C( 2, 2) 10.00000 0.000000 C( 2, 3) 3.000000 0.000000 *( 1, 1) 1500.000 0.000000 *( 1, 2) 1500.000 0.000000 *( 1, 3) 0.000000 0.2000000+308 *( 2, 1) 500.0000 -0.1146654+297 *( 2, 2) 0.00000

18、0 0.1300000+309 *( 2, 3) 5000.000 0.000000 Row Slack or Surplus Dual Price 1 3000.000 -1.000000 2 0.000000 -0.1000000+308 3 0.000000 -0.1000000+308 4 0.000000 -0.1000000+308 5 0.000000 -0.1000000+308 6 0.000000 Infinity 7 -1500.000 Infinity 8 0.000000 -Infinity 9 0.000000 -Infinity 10 0.000000 -Infinity 11 0.000000 0.000000 12 0.000000 0.000000 13 0.000000 0.000000 14 0.000000 0.000000 15 0.000000 0.000000 16 0.000000 0.000000