應(yīng)急運(yùn)輸調(diào)度方案設(shè)計(jì)模型v2

1、20133年中央央民族大大學(xué)數(shù)學(xué)學(xué)建模作作業(yè)論文題目目：應(yīng)急運(yùn)運(yùn)輸調(diào)度度方案設(shè)設(shè)計(jì)模型型參賽隊(duì)員員 ：姓名：吳吳極 學(xué)學(xué)院：理理學(xué)院 專(zhuān)業(yè)：統(tǒng)計(jì)學(xué)學(xué) 年級(jí)級(jí)：111級(jí)姓名：劉劉超 學(xué)學(xué)院：理理學(xué)院 專(zhuān)業(yè)：統(tǒng)計(jì)學(xué)學(xué) 年級(jí)級(jí)：111級(jí)姓名：夏夏浩 學(xué)學(xué)院：理理學(xué)院 專(zhuān)業(yè)：統(tǒng)計(jì)學(xué)學(xué) 年級(jí)級(jí)：111級(jí)應(yīng)急運(yùn)輸輸調(diào)度方方案設(shè)計(jì)計(jì)模型摘 要本題要求求我們求求出每個(gè)個(gè)企業(yè)和和儲(chǔ)備庫(kù)庫(kù)在不同同情況下下給發(fā)放放地點(diǎn)運(yùn)運(yùn)輸救災(zāi)災(zāi)物資的的最優(yōu)調(diào)調(diào)運(yùn)方案案，我們們以每個(gè)個(gè)企業(yè)和和儲(chǔ)備庫(kù)庫(kù)給每個(gè)個(gè)發(fā)放地地點(diǎn)的調(diào)調(diào)運(yùn)量作作為決策策變量，以公路路的長(zhǎng)度度和運(yùn)輸輸成本的的乘積作作為單位位運(yùn)費(fèi)（價(jià)值系系數(shù)）構(gòu)構(gòu)造目標(biāo)標(biāo)函數(shù)。所

2、求問(wèn)問(wèn)題即轉(zhuǎn)轉(zhuǎn)化為最最優(yōu)路徑徑問(wèn)題和和線性規(guī)規(guī)劃問(wèn)題題。在求解問(wèn)問(wèn)題（11）（22）（33）（44）之前前，我們們首先對(duì)對(duì)題目附附件2中中的圖進(jìn)進(jìn)行預(yù)處處理。把把公路的的交點(diǎn)看看成頂點(diǎn)點(diǎn)，每個(gè)個(gè)點(diǎn)之間間的公路路看成線線段，以以公路的的長(zhǎng)度和和運(yùn)輸成成本的乘乘積作為為一條線線段的權(quán)權(quán)重，做做出賦權(quán)權(quán)圖。利利用MAATLAAB軟件件使用FFloyyd算法法計(jì)算出出每個(gè)企企業(yè)和儲(chǔ)儲(chǔ)備庫(kù)到到每個(gè)發(fā)發(fā)放地點(diǎn)點(diǎn)的最優(yōu)優(yōu)路徑（最低單單位運(yùn)費(fèi)費(fèi)和路線線）（見(jiàn)見(jiàn)表4-3-11），解決決最優(yōu)路徑徑問(wèn)題，求出了了目標(biāo)函函數(shù)中的的價(jià)值系系數(shù)。求解問(wèn)題題（1）時(shí)，把把時(shí)間因因素放在在第一位位考慮，首先求求得最快快運(yùn)輸時(shí)

3、時(shí)間t。然后以以運(yùn)輸成成本最低低為目標(biāo)標(biāo)函數(shù)，以調(diào)運(yùn)運(yùn)量小于于等于企企業(yè)和儲(chǔ)儲(chǔ)備庫(kù)儲(chǔ)儲(chǔ)存量，接收量量介于最最低需求求量與最最大需求求量之間間等作為為約束條條件，利利用Liingoo軟件求求解此線線性規(guī)劃劃問(wèn)題的的最優(yōu)解解。由此此得到物物資的最最佳調(diào)運(yùn)運(yùn)方案，包括調(diào)調(diào)運(yùn)量和和調(diào)運(yùn)路路線（見(jiàn)見(jiàn)表4-3-22）。求解問(wèn)題題（2）時(shí)，已已知時(shí)間間t，由由實(shí)際情情況可以以修改約約束條件件，令調(diào)調(diào)運(yùn)量等等于儲(chǔ)存存量，其其他約束束條件不不變。同同樣，利利用Liingoo軟件可可以求出出一個(gè)最最優(yōu)解（見(jiàn)表44-3-3）。求解問(wèn)題題（3）時(shí)，經(jīng)經(jīng)過(guò)計(jì)算算可知企企業(yè)的生生產(chǎn)能力力不能夠夠滿足發(fā)發(fā)放地點(diǎn)點(diǎn)的實(shí)際際需

4、求，我們通通過(guò)企業(yè)業(yè)增產(chǎn)來(lái)來(lái)滿足實(shí)實(shí)際需求求。此時(shí)時(shí)需要新新增三個(gè)個(gè)變量，把問(wèn)題題（1）中的約約束條件件增加幾幾個(gè)約束束條件，利用LLinggo求解解，得到到最佳調(diào)調(diào)運(yùn)方案案（見(jiàn)表44-3-4）。求解問(wèn)題題（4）時(shí)，主主體思路路不變。由于道道路中斷斷，我們們只需要要重新利利用MAATLAAB軟件件求出最最優(yōu)路徑徑和目標(biāo)標(biāo)函數(shù)的的價(jià)值系系數(shù)（見(jiàn)見(jiàn)表4-3-55），再利利用Liingoo軟件求求解線性性規(guī)劃問(wèn)問(wèn)題即可可（見(jiàn)表4-33-6、表4-33-7、表4-33-8）。最后，我我們客觀觀地評(píng)價(jià)價(jià)了該模模型的優(yōu)優(yōu)缺點(diǎn)，并且做做出了相相應(yīng)的改改進(jìn)和推推廣。關(guān)鍵詞：最優(yōu)路徑徑 Flloydd算法 MA

5、TTLABB 線性性規(guī)劃 Linngo 一、問(wèn)題題的提出出與分析析1.1 問(wèn)題重重述在某地區(qū)區(qū)有生產(chǎn)產(chǎn)某種救救災(zāi)物質(zhì)質(zhì)的企業(yè)業(yè)有三家家，設(shè)置置物資發(fā)發(fā)放點(diǎn)八八個(gè)，儲(chǔ)儲(chǔ)備倉(cāng)庫(kù)庫(kù)兩個(gè)。在災(zāi)害害發(fā)生時(shí)時(shí)，企業(yè)業(yè)、各物物資發(fā)放放地點(diǎn)、儲(chǔ)備倉(cāng)倉(cāng)庫(kù)的庫(kù)庫(kù)存情況況，及各各發(fā)放點(diǎn)點(diǎn)的最低低需求和和實(shí)際需需求情況況見(jiàn)附件件1。企業(yè)業(yè)、發(fā)放放點(diǎn)、倉(cāng)倉(cāng)庫(kù)及道道路分布布情況見(jiàn)見(jiàn)附件22。設(shè)該該種物資資的運(yùn)輸輸成本為為高等級(jí)級(jí)公路220元/公里百件，普通公公路122元/公里百件。（1）預(yù)預(yù)案要求求盡快滿滿足各發(fā)發(fā)放點(diǎn)對(duì)對(duì)救災(zāi)物物質(zhì)的最最低需求求，并盡盡量使運(yùn)運(yùn)輸成本本降低。建立數(shù)數(shù)學(xué)模型型，給出出所需要要的時(shí)間間，物資

6、資的調(diào)運(yùn)運(yùn)方案，包括調(diào)調(diào)運(yùn)量和和調(diào)運(yùn)路路線。（2）在在20天內(nèi)內(nèi)，按均均衡配給給的原則則，各發(fā)發(fā)放點(diǎn)可可以得到到多少物物資？給給出相應(yīng)應(yīng)的調(diào)運(yùn)運(yùn)方案。（3）能能否在225天內(nèi)內(nèi)滿足各各發(fā)放點(diǎn)點(diǎn)的實(shí)際際需求？怎樣才才能滿足足各發(fā)放放點(diǎn)的實(shí)實(shí)際需求求？并給給出相應(yīng)應(yīng)的調(diào)運(yùn)運(yùn)方案。（4）在在災(zāi)害發(fā)發(fā)生時(shí)可可能造成成交通中中斷，以以中斷路路段：114-223，111-225，226-227，99-311為例，重新討討論上述述三個(gè)問(wèn)問(wèn)題。1.2 問(wèn)題分分析1.2.1 對(duì)問(wèn)題（1）的的分析要盡快滿滿足各發(fā)發(fā)放地點(diǎn)點(diǎn)對(duì)救災(zāi)災(zāi)物資的的最低需需求，由由現(xiàn)有總總庫(kù)存加加上企業(yè)業(yè)1,22,3tt天的生生產(chǎn)量大大于等

7、于于8個(gè)發(fā)發(fā)放點(diǎn)最最低需求求的不等等式，可可以解出出滿足題題意的最最小時(shí)間間t為88天。接著著在最小小時(shí)間tt=8的的情況下下，求最最小的運(yùn)運(yùn)輸費(fèi)用用，以企企業(yè)1,2,33及儲(chǔ)存存庫(kù)向88個(gè)發(fā)放放點(diǎn)運(yùn)輸輸?shù)奈镔Y資為決策策變量，建立目目標(biāo)函數(shù)數(shù).而建建立目標(biāo)標(biāo)函數(shù)需需要知道道決策變變量對(duì)應(yīng)應(yīng)的價(jià)值值系數(shù),我們根根據(jù)附件件2, 把公路路的交點(diǎn)點(diǎn)看成頂頂點(diǎn)，每每個(gè)點(diǎn)之之間的公公路看成成線段，從而把把地理圖圖轉(zhuǎn)化為為聯(lián)通無(wú)無(wú)向圖.又由題題目條件件高等級(jí)級(jí)公路和和普通公公路的運(yùn)運(yùn)輸費(fèi)用用不同，我們把把每公里里每百件件的運(yùn)費(fèi)費(fèi)和路程程的乘積作作為每條條線段的權(quán)權(quán)數(shù)，根根據(jù)圖論論知識(shí)，將求解解價(jià)值系系數(shù)的問(wèn)

8、問(wèn)題轉(zhuǎn)化化為最優(yōu)優(yōu)路徑問(wèn)問(wèn)題。再再根據(jù)賦賦權(quán)圖制制作權(quán)數(shù)數(shù)矩陣，然后利利用MAATLAAB使用用Flooyd算算法求出出企業(yè)11,2,3及儲(chǔ)儲(chǔ)存庫(kù)到到8個(gè)發(fā)發(fā)放點(diǎn)的的最優(yōu)路徑徑，由此此計(jì)算出出價(jià)值系系數(shù)，再利用用Lingoo軟件在相相關(guān)約束束下求出出目標(biāo)函函數(shù)的最最優(yōu)解，從而得得到物資資的調(diào)運(yùn)運(yùn)方案。1.2.2 對(duì)問(wèn)題（2）的的分析按照均衡衡配給的的原則，求200天后各各發(fā)放點(diǎn)點(diǎn)收到物物資的情情況以及及最佳運(yùn)運(yùn)輸方案案，同樣樣是線性性規(guī)劃問(wèn)問(wèn)題，將將模型中中的t=20，再對(duì)約約束條件件進(jìn)行修修改，使使“庫(kù)存+生產(chǎn)量量=發(fā)放放點(diǎn)接收收量”，得出出最優(yōu)答答案。1.2.3 對(duì)問(wèn)題（3）的的分析要知道

9、225天之之內(nèi)能否否滿足各各發(fā)放點(diǎn)點(diǎn)的實(shí)際際需求，即計(jì)算算25天天的庫(kù)存存和生產(chǎn)產(chǎn)量之和和，與發(fā)發(fā)放點(diǎn)的的最大需求求量進(jìn)行行比較，實(shí)際上上并沒(méi)有有達(dá)到需需求，解解決辦法法是讓企企業(yè)增產(chǎn)產(chǎn)，使之之滿足各各發(fā)放點(diǎn)點(diǎn)的最高高需求，再用線線性規(guī)劃劃模型求求出最優(yōu)優(yōu)方案即即可。1.2.4 對(duì)問(wèn)題（4）的的分析我們要知知道指定定路段中中斷后，上述建建立的數(shù)數(shù)學(xué)模型型是否可可用，也也就是說(shuō)說(shuō)只要檢檢驗(yàn)到模模型中所所選取的的路線是是否經(jīng)過(guò)過(guò)該路段段，如果果不經(jīng)過(guò)過(guò)，則中中斷路線線對(duì)模型型沒(méi)有影影響，若若經(jīng)過(guò)，可將路路段中斷斷后的圖圖采用第第一步的的方式重重新處理理計(jì)算，分別求求解出最最佳運(yùn)輸輸方案。二、基本本

10、假設(shè)2.1 假設(shè)災(zāi)災(zāi)難發(fā)生生時(shí)，企企業(yè)1,2,33只向發(fā)發(fā)放點(diǎn)運(yùn)運(yùn)送物資資，不向向儲(chǔ)備庫(kù)庫(kù)運(yùn)送物物資，而而儲(chǔ)備庫(kù)則是只只出不進(jìn)進(jìn)的向各各個(gè)發(fā)放放點(diǎn)運(yùn)送送物資。2.2 假設(shè)災(zāi)災(zāi)難發(fā)生生當(dāng)天企企業(yè)是生生產(chǎn)物資資的，即即從災(zāi)難難發(fā)生第第一天起，每每天零點(diǎn)點(diǎn)時(shí)每個(gè)個(gè)企業(yè)的的庫(kù)存量量都增加加其日生生產(chǎn)量。2.3 假設(shè)道道路的運(yùn)運(yùn)輸能力力足夠大大，沒(méi)有有運(yùn)輸限限制。2.4 假設(shè)調(diào)調(diào)運(yùn)過(guò)程程中沒(méi)有有衍生災(zāi)災(zāi)害，各各個(gè)路段段道路通通暢，無(wú)無(wú)意外發(fā)發(fā)生。2.5 不考慮慮各點(diǎn)間間的時(shí)間間，假設(shè)設(shè)所有物物資瞬時(shí)時(shí)到達(dá)。2.6 假設(shè)運(yùn)運(yùn)輸時(shí)走走高等級(jí)級(jí)公路和和普通公公路除了了費(fèi)用的的差別外外，在運(yùn)運(yùn)輸結(jié)果果上沒(méi)有有其

11、他差差別。2.7 假設(shè)發(fā)發(fā)放點(diǎn)，企業(yè)和和儲(chǔ)存庫(kù)庫(kù)與公路路的交點(diǎn)點(diǎn)處是重重合的。三、符號(hào)號(hào)說(shuō)明i=1,2,33,4,5，當(dāng)當(dāng)i=11,2,3時(shí)為為企業(yè)11,2,3，當(dāng)當(dāng)i=44,5時(shí)時(shí)為儲(chǔ)備備庫(kù)1,2j=1,2,33,4,5,66,7,8，都都為發(fā)放放點(diǎn):從i運(yùn)運(yùn)到發(fā)放放點(diǎn)j的的物資量量:從i到到j(luò)每百百件的運(yùn)運(yùn)費(fèi): 發(fā)放放地點(diǎn)jj的現(xiàn)有有庫(kù)存：發(fā)放地地點(diǎn)j的的最低需需求：發(fā)放地地點(diǎn)j的的最大需需求: i的的現(xiàn)有庫(kù)庫(kù)存，ii=1,2,33,4,5：企業(yè)ii的日生生產(chǎn)量，i=11,2,3：企業(yè)ii增產(chǎn)后后的日生生產(chǎn)量，其中，t:一個(gè)個(gè)調(diào)運(yùn)方方案的所所需時(shí)間間Z:一個(gè)個(gè)調(diào)運(yùn)方方案的總總運(yùn)費(fèi)四、模型型

12、的建立立與求解解4.1 數(shù)據(jù)處處理 將附件件2中的的公路的的交點(diǎn)看看成頂點(diǎn)點(diǎn)，每個(gè)個(gè)點(diǎn)之間間的公路路看成線線段，形形成一個(gè)個(gè)聯(lián)通無(wú)無(wú)向圖。并且以以每公里里每百件件的運(yùn)費(fèi)費(fèi)和路程程的乘積積作為權(quán)權(quán)重，做做出如下下賦權(quán)圖圖：圖4-11-1 賦賦權(quán)圖4.2 模型建建立根據(jù)題意意，建立立目標(biāo)函函數(shù)Z表表示從企企業(yè)1,2,33和儲(chǔ)備備庫(kù)1,2向88個(gè)發(fā)放放點(diǎn)運(yùn)送送物資的的總費(fèi)用用，根據(jù)各各個(gè)發(fā)放放點(diǎn)的物物資需求求量寫(xiě)出出其約束束條件，如下：目標(biāo)函數(shù)數(shù)：約束條件件：s.t.4.3 模型求求解4.3.1 求求解最優(yōu)優(yōu)路徑，確定：利用MAATLAAB軟件件由Flloydd算法(源程序序見(jiàn)附錄錄1、22)求出從

13、發(fā)出地地點(diǎn)i到到發(fā)放地地點(diǎn)j的的每百件件運(yùn)費(fèi)和和最優(yōu)路路徑，其其結(jié)果如如下表所所示：表4-33-1 各發(fā)出出地點(diǎn)與與發(fā)放地地點(diǎn)每百件件運(yùn)費(fèi)和和最優(yōu)路徑表 *企業(yè)1企業(yè)2企業(yè)3儲(chǔ)備庫(kù)11儲(chǔ)備庫(kù)22發(fā)放點(diǎn)1118448,224,226,225,115,442,2286966,411,422,28826888,334,332,339,330,229,22822772,227,440,66,411,422,28814664,330,229,228發(fā)放點(diǎn)2215000,224,226,119,118,22318884,441,442,115,118,22337440,334,11,2,7,227,226

14、,119,118,22319880,227,226,119,118,22334220,330,229,228,442,115,118,223發(fā)放點(diǎn)3340880,224,226,227,99,311,322,35536772,441,66,400,9,31,32,3514776,334,332,33528880,227,99,311,322,35521000,330,339,332,335發(fā)放點(diǎn)4423004,224,226,227,99,31118996,441,66,400,9,319000,344,322,31111004,227,99,31115224,330,339,332,331發(fā)

15、放點(diǎn)5515660,224,220,22224772,441,442,115,118,119,22238000,334,11,2,7,227,226,119,22220440,227,226,119,22240008,330,229,228,442,115,118,119,222發(fā)放點(diǎn)6634444,224,226,227,99,2,3,33630336,441,66,400,9,2,33,36617440,334,11,333,36622444,227,99,2,3,33629664,330,339,332,334,11,333,366發(fā)放點(diǎn)7725668,224,226,225,115,4

16、42,228,22914116,441,442,228,22919668,334,332,339,330,22921660,227,440,66,4,297444,300,299發(fā)放點(diǎn)8837220,224,226,227,99,311,322,38833112,441,66,400,9,31,32,3811116,334,332,33825220,227,99,311,322,38817440,330,339,332,338* 本表表中 符號(hào)表示示數(shù)組，例如118488,244,266,255,155,422,288，第一個(gè)個(gè)元素118488表示該該路徑的的每百件件運(yùn)費(fèi)，后面的的元素224,

17、226,225,115,442,228表示示該最優(yōu)優(yōu)路徑經(jīng)經(jīng)過(guò)的節(jié)節(jié)點(diǎn)序號(hào)號(hào)（包括括起點(diǎn)和和終點(diǎn)）（見(jiàn)圖4-1-1）。下同。4.3.2 求解問(wèn)問(wèn)題（11）由于要使使物資盡盡快到達(dá)達(dá)發(fā)放地地點(diǎn)，則則應(yīng)該首首先生產(chǎn)產(chǎn)出滿足足各發(fā)放放地點(diǎn)最最低需求求的物資資量，再再進(jìn)行運(yùn)運(yùn)輸規(guī)劃劃。此時(shí)時(shí)有：發(fā)放地點(diǎn)點(diǎn)最低總總需求量量=335500發(fā)放地點(diǎn)點(diǎn)，企業(yè)業(yè)與儲(chǔ)備備庫(kù)現(xiàn)儲(chǔ)儲(chǔ)存總量量=+=28840則有：解得：所以，應(yīng)應(yīng)該取最最優(yōu)調(diào)運(yùn)運(yùn)時(shí)間，再將調(diào)調(diào)運(yùn)時(shí)間間代入44.2的的線性規(guī)規(guī)劃模型型中，如如下：目標(biāo)函數(shù)數(shù)：約束條件件：s.t.利用Liingoo軟件（源源程序見(jiàn)見(jiàn)附錄33）對(duì)其其求解，得到最最優(yōu)解及及調(diào)運(yùn)

18、路路線為：表4-33-2 各發(fā)出地點(diǎn)點(diǎn)向發(fā)放放地點(diǎn)運(yùn)運(yùn)送量和和最優(yōu)路徑徑表 *企業(yè)1企業(yè)2企業(yè)3儲(chǔ)備庫(kù)11儲(chǔ)備庫(kù)22發(fā)放點(diǎn)110300（42）00160(29)發(fā)放點(diǎn)22140（26-19-18）00410（26-19-18）0發(fā)放點(diǎn)330000280(39-32)發(fā)放點(diǎn)44000320(9)0發(fā)放點(diǎn)55300（20）0000發(fā)放點(diǎn)66000260(9-22-3)0發(fā)放點(diǎn)770000470(直達(dá))發(fā)放點(diǎn)8800240（32）0290(39-32)* 本表表中（）符號(hào)表表示路徑徑，例如1440（226-119-118），括號(hào)外外1400表示該路路徑的運(yùn)運(yùn)送量，（266-199-188）表示示該路

19、徑徑經(jīng)過(guò)的的節(jié)點(diǎn)序序號(hào)（即表4-3-11中對(duì)應(yīng)應(yīng)的最優(yōu)優(yōu)路徑，不包括括起點(diǎn)和和終點(diǎn)）。下同。4.3.3 求解問(wèn)問(wèn)題（22）發(fā)放地點(diǎn)點(diǎn)最低總總需求量量=335500；發(fā)放地點(diǎn)點(diǎn)實(shí)際總總需求量量=556000；而當(dāng)時(shí)，有；即當(dāng)時(shí)，企業(yè)和和儲(chǔ)備庫(kù)庫(kù)的物資資總量肯肯定能夠夠滿足各各個(gè)發(fā)放放地點(diǎn)的的最低需需求，而而且不超超過(guò)各個(gè)個(gè)發(fā)放點(diǎn)點(diǎn)的最高高需求。從實(shí)際際情況出出發(fā)，不不可能讓讓物資堆堆積在企企業(yè)之中中，所以以按照均均勻配給給的原則則，將企企業(yè)1，2,33和儲(chǔ)備備庫(kù)1,2以及及發(fā)放點(diǎn)點(diǎn)原來(lái)所所存的所所有物資資全部發(fā)發(fā)放，得得到新的的規(guī)劃模模型如下下：目標(biāo)函數(shù)數(shù)：約束條件件：s.t. 利用Liingo

20、o軟件（源程序序見(jiàn)附錄錄4）解解出此時(shí)時(shí)的最優(yōu)優(yōu)解以及及各個(gè)發(fā)發(fā)放點(diǎn)所所得到的的物資和和調(diào)運(yùn)方方案。表4-33-3 20天內(nèi)各發(fā)出出地點(diǎn)向發(fā)發(fā)放地點(diǎn)點(diǎn)運(yùn)送量量和最優(yōu)路徑徑表 *企業(yè)1企業(yè)2企業(yè)3儲(chǔ)備庫(kù)11儲(chǔ)備庫(kù)22發(fā)放點(diǎn)110660（42）00100(29)發(fā)放點(diǎn)22480（26-19-18）00370（26-19-18）0發(fā)放點(diǎn)3300280（32）00發(fā)放點(diǎn)44000370(9)0發(fā)放點(diǎn)55440（20）0000發(fā)放點(diǎn)66000260(9-22-3)0發(fā)放點(diǎn)770000570(直達(dá))發(fā)放點(diǎn)8800200（32）0530(39-32)4.3.4 求解問(wèn)問(wèn)題（33）要知道225天之之內(nèi)能否否滿

21、足各各發(fā)放點(diǎn)點(diǎn)的實(shí)際際需求，即計(jì)算算25天天的庫(kù)存存和生產(chǎn)產(chǎn)量之和和，與發(fā)發(fā)放點(diǎn)的的實(shí)際需需求量進(jìn)進(jìn)行比較較。由題題意：發(fā)放地點(diǎn)點(diǎn)實(shí)際總總需求量量=556000；發(fā)放地點(diǎn)點(diǎn)，企業(yè)業(yè)與儲(chǔ)備備庫(kù)現(xiàn)儲(chǔ)儲(chǔ)存總量量=+=28840；25天企企業(yè)1,2,33的總產(chǎn)產(chǎn)量=22550；由于22250+28440D(ii,k)+D(k,jj) DD(i,j)=D(ii,k)+D(k,jj); ppathh(i,j)=patth(ii,k); eend eend eendend p=ssp;mp=ssp;for k=11:n if mmp=ep d=patth(mmp,eep); p=p,d; mpp=d; en

22、ddendd=D(sp,ep);pathh=p;%=ffloyyd.mm文件件結(jié)束=2. 問(wèn)問(wèn)題（11）利用用Flooyd算算法給出出系數(shù)和和各點(diǎn)間間最優(yōu)路徑徑的Maatlaab程序序%=rrun001.mm文件件開(kāi)始=loadd(llinees0.matt);%“l(fā)innes00.maat”為事先先錄入的的42*42的的權(quán)數(shù)矩矩陣，包包含每點(diǎn)點(diǎn)的 距距離*單單位路費(fèi)費(fèi) 數(shù)據(jù)據(jù)i1=11;j1=11;resuult=celll(55,8);for i=24,41,34,27,30 forr j=288,233,355,311,222,366,299,388 dd,paath=flloydd(l

23、iiness0,ii,j); %從“flooyd.m”調(diào)用ffloyyd算法法 reesulltii1,jj1=d,patth; j11=j11+1; endd i1=i1+1;jj1=11;end%=rrun001.mm文件件結(jié)束=3. 問(wèn)問(wèn)題（11）計(jì)算算最低成成本運(yùn)輸輸方案的的Linngo命命令modeel:!5發(fā)點(diǎn)點(diǎn)8收點(diǎn)運(yùn)運(yùn)輸問(wèn)題題;setss:wareehouusess/whh1.wh55/:ccapaacitty;!企業(yè)1,2,33與儲(chǔ)備備庫(kù)1,2的儲(chǔ)儲(chǔ)存量;venddorss/v11.vv8/:demmandd1,ddemaand22;!發(fā)放地地點(diǎn)的最最低需求求量和最最高需求求

24、量;linkks(wwareehouusess,veendoors):coost,vollumee;!權(quán)數(shù)矩矩陣和決決策變量量;endssetss!目標(biāo)函函數(shù);min=suum(llinkks:ccostt*voolumme);!需求約約束（儲(chǔ)儲(chǔ)存量在在最低需需求量和和最高需需求量之之間）;forr(veendoors(J): ssum(warrehoousees(II):vvoluume(I,JJ)=deemannd1(J); ssum(warrehoousees(II):vvoluume(I,JJ)=deemannd2(J); );!產(chǎn)量約約束（發(fā)發(fā)放量不不超過(guò)儲(chǔ)儲(chǔ)存量）;forr(waa

25、rehhousses(I): ssum(venndorrs(JJ):vvoluume(I,JJ)=deemannd1(J); ssum(warrehoousees(II):vvoluume(I,JJ)=40;extrra(22)=30;extrra(33)=20;!產(chǎn)量約約束（發(fā)發(fā)放量等等于儲(chǔ)存存量）;forr(waarehhousses(I): ssum(venndorrs(JJ):vvoluume(I,JJ)=cappaciity(I)+25*exttra(I); );!這里是是數(shù)據(jù);dataa: caapaccityy=1220 660 880 110000 12200; deemann

26、d=7760 8500 5880 3370 9000 4660 5570 7300; coost=18448 115000 40080 23004 115600 34444 25668 337200 6996 118844 36672 18996 224722 30036 14116 333122 26888 339844 14476 9000 440444 17740 19668 111166 22772 119800 28880 11004 220400 22244 21660 225200 14664 334200 21100 15224 440088 26676 7444 117400

27、;endddataaend最佳調(diào)運(yùn)運(yùn)方案： Gllobaal ooptiimall sooluttionn fooundd. Obbjecctivve vvaluue: 7722776000. Tootall soolveer iiterratiionss: 115 Vaariaablee VValuue Redduceed CCostt CAAPACCITYY( WWH1) 1120.00000 0.00000000 CAAPACCITYY( WWH2) 660.0000000 0.00000000 CAAPACCITYY( WWH3) 880.0000000 0.00000000 CAAP

28、ACCITYY( WWH4) 110000.0000 0.00000000 CAAPACCITYY( WWH5) 112000.0000 0.00000000 DEEMANND( V1) 7760.00000 0.00000000 DEEMANND( V2) 8850.00000 0.00000000 DEEMANND( V3) 5580.00000 0.00000000 DEEMANND( V4) 3370.00000 0.00000000 DEEMANND( V5) 9900.00000 0.00000000 DEEMANND( V6) 4460.00000 0.00000000 DEEM

29、ANND( V7) 5570.00000 0.00000000 DEEMANND( V8) 7730.00000 0.00000000 EEXTRRA( E1) 440.0000000 0.00000000 EEXTRRA( E2) 330.0000000 0.00000000 EEXTRRA( E3) 440.4400000 0.00000000 EEXTRRA( E4) 00.00000000 0.00000000 EEXTRRA( E5) 00.00000000 0.00000000 COOST( WHH1, V1) 118488.0000 0.00000000 COOST( WHH1,

30、 V2) 115000.0000 0.00000000 COOST( WHH1, V3) 440800.0000 0.00000000 COOST( WHH1, V4) 223044.0000 0.00000000 COOST( WHH1, V5) 115600.0000 0.00000000 COOST( WHH1, V6) 334444.0000 0.00000000 COOST( WHH1, V7) 225688.0000 0.00000000 COOST( WHH1, V8) 337200.0000 0.00000000 COOST( WHH2, V1) 6696.00000 0.00

31、000000 COOST( WHH2, V2) 118844.0000 0.00000000 COOST( WHH2, V3) 336722.0000 0.00000000 COOST( WHH2, V4) 118966.0000 0.00000000 COOST( WHH2, V5) 224722.0000 0.00000000 COOST( WHH2, V6) 330366.0000 0.00000000 COOST( WHH2, V7) 114166.0000 0.00000000 COOST( WHH2, V8) 333122.0000 0.00000000 COOST( WHH3,

32、V1) 226888.0000 0.00000000 COOST( WHH3, V2) 339844.0000 0.00000000 COOST( WHH3, V3) 114766.0000 0.00000000 COOST( WHH3, V4) 9900.00000 0.00000000 COOST( WHH3, V5) 440444.0000 0.00000000 COOST( WHH3, V6) 117400.0000 0.00000000 COOST( WHH3, V7) 119688.0000 0.00000000 COOST( WHH3, V8) 111166.0000 0.000

33、00000 COOST( WHH4, V1) 222722.0000 0.00000000 COOST( WHH4, V2) 119800.0000 0.00000000 COOST( WHH4, V3) 228800.0000 0.00000000 COOST( WHH4, V4) 111044.0000 0.00000000 COOST( WHH4, V5) 220400.0000 0.00000000 COOST( WHH4, V6) 222444.0000 0.00000000 COOST( WHH4, V7) 221600.0000 0.00000000 COOST( WHH4, V

34、8) 225200.0000 0.00000000 COOST( WHH5, V1) 114644.0000 0.00000000 COOST( WHH5, V2) 334200.0000 0.00000000 COOST( WHH5, V3) 221000.0000 0.00000000 COOST( WHH5, V4) 115244.0000 0.00000000 COOST( WHH5, V5) 440088.0000 0.00000000 COOST( WHH5, V6) 226766.0000 0.00000000 COOST( WHH5, V7) 7744.00000 0.0000

35、0000 COOST( WHH5, V8) 117400.0000 0.00000000 VVOLUUME( WHH1, V1) 00.00000000 15536.0000 VVOLUUME( WHH1, V2) 8800.00000 0.00000000 VVOLUUME( WHH1, V3) 00.00000000 25580.0000 VVOLUUME( WHH1, V4) 00.00000000 16680.0000 VVOLUUME( WHH1, V5) 3320.00000 0.00000000 VVOLUUME( WHH1, V6) 00.00000000 16680.0000

36、 VVOLUUME( WHH1, V7) 00.00000000 24424.0000 VVOLUUME( WHH1, V8) 00.00000000 25580.0000 VVOLUUME( WHH2, V1) 7760.00000 0.00000000 VVOLUUME( WHH2, V2) 550.0000000 0.00000000 VVOLUUME( WHH2, V3) 00.00000000 17788.0000 VVOLUUME( WHH2, V4) 00.00000000 8888.000000 VVOLUUME( WHH2, V5) 00.00000000 5228.0000

37、00 VVOLUUME( WHH2, V6) 00.00000000 8888.000000 VVOLUUME( WHH2, V7) 00.00000000 8888.000000 VVOLUUME( WHH2, V8) 00.00000000 17788.0000 VVOLUUME( WHH3, V1) 00.00000000 24400.0000 VVOLUUME( WHH3, V2) 00.00000000 25508.0000 VVOLUUME( WHH3, V3) 5580.00000 0.00000000 VVOLUUME( WHH3, V4) 00.00000000 3000.0

38、00000 VVOLUUME( WHH3, V5) 00.00000000 25508.0000 VVOLUUME( WHH3, V6) 4410.00000 0.00000000 VVOLUUME( WHH3, V7) 00.00000000 18848.0000 VVOLUUME( WHH3, V8) 1100.00000 0.00000000 VVOLUUME( WHH4, V1) 00.00000000 14480.0000 VVOLUUME( WHH4, V2) 00.00000000 0.00000000 VVOLUUME( WHH4, V3) 00.00000000 9000.0

39、00000 VVOLUUME( WHH4, V4) 3370.00000 0.00000000 VVOLUUME( WHH4, V5) 5580.00000 0.00000000 VVOLUUME( WHH4, V6) 550.0000000 0.00000000 VVOLUUME( WHH4, V7) 00.00000000 15536.0000 VVOLUUME( WHH4, V8) 00.00000000 9000.000000 VVOLUUME( WHH5, V1) 00.00000000 5552.000000 VVOLUUME( WHH5, V2) 00.00000000 1332

40、0.0000 VVOLUUME( WHH5, V3) 00.00000000 0.00000000 VVOLUUME( WHH5, V4) 00.00000000 3000.000000 VVOLUUME( WHH5, V5) 00.00000000 18848.0000 VVOLUUME( WHH5, V6) 00.00000000 3112.000000 VVOLUUME( WHH5, V7) 5570.00000 0.00000000 VVOLUUME( WHH5, V8) 6630.00000 0.00000000 Roww Sllackk orr Suurpllus DDuall P

41、rricee 11 7722776000. -1.00000000 22 00.00000000 -2888.000000 33 00.00000000 -14476.0000 44 00.00000000 -14476.0000 55 00.00000000 -6000.000000 66 00.00000000 -15536.0000 77 00.00000000 -17740.0000 88 00.00000000 -1220.000000 99 00.00000000 -11116.0000 100 00.00000000 -1226000.000 111 00.00000000 -1

42、556000.000 122 00.00000000 -6000.000000 133 00.00000000 -1002000.000 144 220.4400000 0.00000000 155 00.00000000 -244.0000000 166 00.00000000 -4008.000000 177 00.00000000 0.00000000 188 00.00000000 -5004.000000 199 00.00000000 -6224.0000006. 問(wèn)問(wèn)題（44）修改改權(quán)數(shù)矩矩陣重新新計(jì)算系系數(shù)和各各點(diǎn)間最最優(yōu)路徑徑的Maatlaab程序序%=rrun004.mm文件件開(kāi)始=loadd(llinees0.matt)