火力電廠經(jīng)營(yíng)問(wèn)題的優(yōu)化方案

.../火力電廠經(jīng)營(yíng)問(wèn)題的優(yōu)化方案摘要針對(duì)火力電廠的經(jīng)營(yíng)優(yōu)化問(wèn)題,本文建立了采購(gòu)資源安排的優(yōu)化模型,在火力電廠獲取長(zhǎng)遠(yuǎn)經(jīng)濟(jì)效益的前提下,分析解決了所給的三個(gè)問(wèn)題。對(duì)于問(wèn)題一,利用matlab軟件擬合出耗油量和柴油價(jià)格規(guī)律曲線,預(yù)測(cè)了今年每月的耗油量及每星期的油價(jià)。建立非線性規(guī)劃模型,利用預(yù)測(cè)結(jié)果計(jì)算得到具體的購(gòu)油方案,以及所用費(fèi)用,最后算出考慮支付利息是的總費(fèi)用。針對(duì)問(wèn)題二,采用遠(yuǎn)期和約來(lái)購(gòu)買柴油時(shí),利用圖論的知識(shí),把從一月到第十二個(gè)月簽訂和約的多種組合方式設(shè)定為求兩點(diǎn)之間最小距離的問(wèn)題,建立圖論——最短路徑模型。用matlab編程計(jì)算得到十二個(gè)月的遠(yuǎn)期購(gòu)油計(jì)劃：先簽訂一個(gè)月的購(gòu)油和約,到期后再簽訂兩個(gè)月的購(gòu)油和約,如此方式再簽訂六個(gè)月、一個(gè)月、兩個(gè)月的購(gòu)油和約。共耗資金70704112元。針對(duì)問(wèn)題三,分析15天最少柴油儲(chǔ)備合理性,從火電廠的經(jīng)濟(jì)效益出發(fā)考慮,分別利用問(wèn)題1、2所建立的模型,計(jì)算比較在不同的柴油儲(chǔ)備情況下,購(gòu)買柴油總費(fèi)用的變化,費(fèi)用越小,合理性就越高。最后得出政府規(guī)定的15天的最少柴油儲(chǔ)備不合理,最佳天數(shù)為13天。關(guān)鍵詞：擬合非線性規(guī)劃模型圖論1.問(wèn)題的提出1.1基本情況某火力發(fā)電廠負(fù)責(zé)一小城市生產(chǎn)和生活用電的供應(yīng)。該電廠發(fā)電機(jī)以柴油為燃料。發(fā)電廠有多臺(tái)發(fā)電機(jī),可以根據(jù)用電需求決定啟用多少臺(tái)發(fā)電機(jī)適應(yīng)點(diǎn)亮的季節(jié)性和隨機(jī)性的波動(dòng)。正常情況下柴油可以在市場(chǎng)上隨時(shí)買到并立即運(yùn)到電廠,但柴油價(jià)格隨市場(chǎng)是經(jīng)常波動(dòng)的。政府規(guī)定：發(fā)電廠必須至少有15天的滿足正常用電需求的燃油儲(chǔ)備以應(yīng)付用電需求的隨機(jī)波動(dòng)和不能按時(shí)購(gòu)進(jìn)燃油的偶然事件。〔1截止到去年底該廠前48個(gè)月的每月耗油量〔噸：1120,1180,1320,1290,1210,1350,1480,1480,1360,1190,1040,1180,1150,1260,1410,1350,1250,1490,1700,1700,1580,1330,1140,1400,1450,1500,1780,1630,1720,1780,1990,1990,1840,1620,1460,1660,1710,1800,1930,1810,1830,2180,2300,2420,2090,1910,1720,1940〔2截止到去年底前50個(gè)星期的市場(chǎng)柴油價(jià)格〔噸運(yùn)輸費(fèi),元：2494,2490,2499,2505,2515,2490,2476,2488,2504,2507,2530,2537,2550,2562,2560,2574,2640,2600,2613,2604,2616,2608,2589,2590,2589,2574,2577,2579,2574,2573,2576,2589,2578,2577,2572,2575,2568,2575,2570,2576,2573,2575,2596,2611,2629,2633,2628,2618,2622,2627〔3每購(gòu)買一次柴油所需固定費(fèi)用〔手續(xù)費(fèi)等為5萬(wàn)元。一噸柴油儲(chǔ)存一天所花費(fèi)用為0.1元,該廠從銀行貸款所用年利率〔連續(xù)復(fù)利為5%。1.2需要解決的問(wèn)題問(wèn)題1：請(qǐng)?jiān)诿磕昴瓿鯙樵搹S制定一個(gè)經(jīng)濟(jì)上盡量合理的柴油購(gòu)買與庫(kù)存方案〔即一年購(gòu)油多少噸,并根據(jù)此估算出發(fā)電廠在該年購(gòu)買柴油的費(fèi)用。請(qǐng)按所有資金不支付利息和支付利息兩種情況討論。問(wèn)題2：若采用所謂遠(yuǎn)期和約來(lái)購(gòu)買柴油,即發(fā)電廠與售油公司事先計(jì)算出一個(gè)雙方認(rèn)為的油價(jià)并簽訂一個(gè)合和約,將來(lái)購(gòu)油時(shí)就采用這個(gè)油價(jià),這樣的遠(yuǎn)期和約有期限為一個(gè)月、兩個(gè)月、三個(gè)月、四個(gè)月、六個(gè)月、一年等共六種。請(qǐng)你為發(fā)電廠選擇一種最有利的遠(yuǎn)期和約,使得費(fèi)用最省。問(wèn)題三：你認(rèn)為政府規(guī)定的15天的最少柴油儲(chǔ)備合理嗎？給出改進(jìn)意見(jiàn)。2.問(wèn)題的分析對(duì)于火力電廠的經(jīng)營(yíng)優(yōu)化問(wèn)題,這是一個(gè)帶有復(fù)雜約束條件的優(yōu)化與規(guī)劃問(wèn)題。在當(dāng)今時(shí)代,企業(yè)采購(gòu)資源規(guī)劃是企業(yè)生產(chǎn)經(jīng)營(yíng)規(guī)劃的重要組成部分。在采購(gòu)資源的安排中,不僅要注重?cái)?shù)量,更重要的是采購(gòu)質(zhì)量的提高以及如何利用好現(xiàn)有的資源獲得最大的效益,這對(duì)企業(yè)的發(fā)展至關(guān)重要。解決這個(gè)問(wèn)題就是要進(jìn)行企業(yè)預(yù)測(cè)〔特別是需求預(yù)測(cè),例如制定銷售計(jì)劃、生產(chǎn)計(jì)劃、庫(kù)存計(jì)劃以及組織計(jì)劃都可以用到這些預(yù)測(cè)結(jié)果。問(wèn)題1的分析要求每年年初為該廠制定一個(gè)經(jīng)濟(jì)上盡量合理的柴油購(gòu)買與庫(kù)存方案,這個(gè)問(wèn)題就轉(zhuǎn)化為"今年購(gòu)油多少次,每次夠有多少噸"的問(wèn)題,題目給出了前四年火電廠對(duì)柴油的需求,以及截止到去年年底前50個(gè)星期的市場(chǎng)柴油價(jià)格,很明顯,需要我們根據(jù)給出的以往的統(tǒng)計(jì)數(shù)據(jù),對(duì)今年的情況作出預(yù)測(cè)。我們對(duì)已知數(shù)據(jù)進(jìn)行擬合,運(yùn)用matlab軟件,畫出前四年的月耗油量曲線,據(jù)此擬合出今年的月耗油量曲線,進(jìn)而算出今年每月的耗油量。同理,很據(jù)往年的市場(chǎng)柴油價(jià)格,預(yù)計(jì)出今年的柴油價(jià)格在各星期的波動(dòng)情況。之后,又對(duì)采購(gòu)的次數(shù)、采購(gòu)量進(jìn)行了優(yōu)化,預(yù)測(cè)了發(fā)電廠在今年購(gòu)買柴油的費(fèi)用；發(fā)電廠購(gòu)買柴油的費(fèi)用可以表示為：購(gòu)買柴油的總費(fèi)用=〔每次購(gòu)買的數(shù)量*〔該次柴油的價(jià)格+〔每次購(gòu)買的手續(xù)費(fèi)+〔柴油的儲(chǔ)存費(fèi)用；約束條件為：〔每月購(gòu)買的柴油量+〔上月余下的庫(kù)存量—〔當(dāng)月的消費(fèi)量≥〔下個(gè)月預(yù)測(cè)消費(fèi)量的一半。II.問(wèn)題2的分析在本問(wèn)題中,因?yàn)檫h(yuǎn)期合約期限有一個(gè)月、兩個(gè)月、三個(gè)月、四個(gè)月、六個(gè)月、一年共六種情況,而一年當(dāng)中各種合約的組合方式有多種,如何使一年購(gòu)買柴油的費(fèi)用最省,其實(shí)也就是如何安排一年當(dāng)中各種合約的組織方式,使得各種合約條件下費(fèi)用的總和最小。最小費(fèi)用就是和約油價(jià)與需求量的乘積之和最小。采用遠(yuǎn)期和約來(lái)購(gòu)買柴油,考慮到與供應(yīng)商建立和約后,火電廠與供應(yīng)商之間在簽約允許的時(shí)間內(nèi),供應(yīng)商能夠隨時(shí)提供火電廠的柴油需求,這樣就不必考慮每次購(gòu)買柴油后的儲(chǔ)存問(wèn)題了。要求我們選擇一個(gè)合理的購(gòu)油價(jià)格和購(gòu)油合同期限。在決定購(gòu)油價(jià)格時(shí),要使雙方都滿意,我們?nèi)「鱾€(gè)時(shí)間段內(nèi)的平均值,這樣,最合理的購(gòu)油和約就轉(zhuǎn)化為求在一年之內(nèi)所有和約價(jià)格與需求量成績(jī)最小的那一組和約。我們利用圖論的知識(shí),把從一月到第十二個(gè)月簽訂和約的多種組合方式設(shè)定為求兩點(diǎn)之間最小距離的問(wèn)題,建立圖論——最短路徑模型。用matlab編程計(jì)算得到十二個(gè)月的遠(yuǎn)期購(gòu)油計(jì)劃。III.問(wèn)題3的分析要求對(duì)政府規(guī)定的15天的最少柴油儲(chǔ)備的合理性進(jìn)行分析評(píng)價(jià)。對(duì)于這個(gè)問(wèn)題,從火電廠的經(jīng)濟(jì)效益出發(fā)考慮,利用問(wèn)題1所建立的模型,計(jì)算比較在不同的柴油儲(chǔ)備情況下,購(gòu)買柴油總費(fèi)用的變化,費(fèi)用越小,合理性就越高；再利用問(wèn)題2所建立的模型,采用遠(yuǎn)期和約來(lái)購(gòu)買柴油的方法,在合約規(guī)定內(nèi),供應(yīng)商應(yīng)當(dāng)按照電廠的需要提供柴油。計(jì)算在不同柴油儲(chǔ)備情況下,每月的購(gòu)油量發(fā)生了怎樣的變化,比較使用不同遠(yuǎn)期和約購(gòu)油的費(fèi)用,得到合理的購(gòu)油方案。3.模型的假設(shè)〔1每年各星期的柴油市場(chǎng)價(jià)格走勢(shì)基本一致?！?為計(jì)算方便,假設(shè)一個(gè)月為四個(gè)星期,15天記為半個(gè)月?！?發(fā)電廠與售油公司簽訂的雙方都認(rèn)為合理的購(gòu)油價(jià)格,即為合同期限內(nèi)預(yù)測(cè)油價(jià)的平均值?！?發(fā)電廠與售油公司簽訂合同后,購(gòu)買柴油手續(xù)費(fèi)固定,且與購(gòu)油次數(shù)無(wú)關(guān)?！?售油公司油量充足,能滿足電廠需求?！?從銀行貸款所用的年利率保持不變,且貸款按年計(jì)算。4.符號(hào)約定表示每月的用油量；表示每月的購(gòu)油量；表示每月初擁有的油量；表示0-1變量；表示每月的平均油價(jià)；5.對(duì)問(wèn)題一的求解模型建立與求解〔1.每月耗油量的確定題目中給出了截至去年年底該火電廠前48個(gè)月的每月耗油量,我們利用matlab軟件對(duì)每年的十二個(gè)月的耗油量進(jìn)行擬合對(duì)比,發(fā)現(xiàn)擬合8次多項(xiàng)式效果最佳,所以均用8次多項(xiàng)式擬合。第一年的各月耗油量方程：擬合曲線如圖1所示圖1.第一年的耗油量擬合曲線〔2第二年的各月耗油量方程：擬合曲線如圖2所示圖2.第二年的耗油量擬合曲線〔3第三年的各月耗油量方程：擬合曲線如圖3所示圖3.第三年的耗油量擬合曲線〔4第四年的各月耗油量方程：擬合曲線如圖4所示圖4.第四年的耗油量擬合曲線由這四條曲線我們發(fā)現(xiàn),每個(gè)月的耗油規(guī)律十分相似,符合實(shí)際。又發(fā)現(xiàn),每年的最后一個(gè)數(shù)據(jù)與其次年的第一個(gè)數(shù)據(jù)十分接近,把四年中的三次交替差值取平均值作為第四年最后一個(gè)月的耗油量與今年第一個(gè)月耗油量的差值,算出其為23噸,即今年第一個(gè)月的耗油量為1963噸。將上面四個(gè)方程去掉常數(shù)項(xiàng)求平均值,再代入初始值x=1時(shí)y=1963,求的常數(shù)項(xiàng)。最后得到今年每個(gè)月耗油量的方程為：擬合曲線如圖5所示圖5.今年的耗油量擬合曲線帶入初始值求得今年每月的耗油量如表5.1表5.1.1月份123456今年的耗油量196320402218212121152297月份789101112今年的耗油量248324932330211419462150〔2.柴油價(jià)格的確定根據(jù)權(quán)威專家分析,每年柴油價(jià)格的變化都有一定的規(guī)律,即價(jià)格曲線走勢(shì)大致相同。由此,我們可以利用去年的柴油價(jià)格來(lái)預(yù)測(cè)今年柴油價(jià)格的大致走向。首先對(duì)去年的柴油價(jià)格進(jìn)行擬合。為了準(zhǔn)確擬合我們采用交叉擬合,即將50個(gè)油價(jià)數(shù)據(jù)分成四段,第一段為對(duì)前七個(gè)數(shù)據(jù),第二段為第7——17個(gè)數(shù)據(jù),第三段為第17——41個(gè)數(shù)據(jù),第四段為第41——50個(gè)數(shù)據(jù),通過(guò)擬合找到最佳擬合次數(shù),對(duì)第一段數(shù)據(jù)進(jìn)行八次擬合,第二段進(jìn)行三次擬合,第三段進(jìn)行八次擬合,第四段進(jìn)行四次擬合,得到的擬合圖如下：圖6.第一段數(shù)據(jù)八次擬合曲線圖7.第二段數(shù)據(jù)三次擬合曲線圖8.第三段數(shù)據(jù)八次擬合曲線圖9.第四段數(shù)據(jù)四次擬合曲線由于每條曲線的最后一個(gè)數(shù)據(jù)與下一條曲線的第一個(gè)數(shù)據(jù)相差非常小,所以我們將第四段曲線的最后一個(gè)數(shù)據(jù)作為今年第一個(gè)星期的油價(jià),加入原第一段曲線的擬合函數(shù),得到的函數(shù)可以近似看做今年的前七個(gè)星期的油價(jià)函數(shù),而后可以求得前七個(gè)星期的所有油價(jià),依此類推依次算出后面43個(gè)星期的所有油價(jià),得到今年每個(gè)星期的油價(jià)表,如下表5.2表5.2.1月份星期一星期二星期三星期四星期一月2627262326322638二月2648262326092630.2三月2645.12655.62663.62670.9四月2679.52691.32708.22732五月2764.62773.127312736.2六月2738.62735.527282718.7七月2710.22704.22701.32701.1八月2702.827052706.62707九月2705.82703.627012698.9十月26982698.42699.72700.7十一月2700.22704.827222741十二月2754.92760.12756.92749.5另外兩星期2745.52756.3擬合曲線得到圖10如下：圖10.預(yù)測(cè)的今年每月油價(jià)曲線圖〔3.模型的建立及求解考慮到火電廠的最大利益,和購(gòu)買柴油和庫(kù)存費(fèi)用最低以及滿足15天的庫(kù)存限制,根據(jù)每月的耗油量和每月的油價(jià),我們建立非線性微分方程模型：目標(biāo)函數(shù)：約束條件：用油量的約束條件每月初所擁有的擁有量的約束其中變量約束計(jì)算流程圖如圖11所示：N滿足各約束條件開(kāi)始輸入yi,xi,ai,RiN滿足各約束條件開(kāi)始輸入yi,xi,ai,Rii=1,2,……,12N結(jié)束輸出x、y、R結(jié)束輸出x、y、R圖11.模型一的算法流程圖通過(guò)計(jì)算我們求得不支付利息時(shí)的購(gòu)油方案：第一個(gè)月初購(gòu)油1998.5噸,第二個(gè)月初購(gòu)油24339噸,其余的時(shí)間不購(gòu)油,此時(shí)花費(fèi)的總費(fèi)用為69263700元；考慮支付利息是花費(fèi)的總費(fèi)用為72726885元。〔4.模型的分析與評(píng)價(jià)由上述計(jì)算結(jié)果我們發(fā)現(xiàn)購(gòu)油只在一二月份進(jìn)行,其他時(shí)間均不購(gòu)油,從圖10預(yù)測(cè)的今年每月油價(jià)曲線圖也可以看出第一二個(gè)月的油價(jià)比其它十個(gè)月都低,并且題目中要求必須有15天的油量?jī)?chǔ)備,所以第一個(gè)月我們必須購(gòu)油,不然就不滿足要求,而第二個(gè)月油價(jià)最低。綜合考慮油的儲(chǔ)備問(wèn)題,可以解釋為,第一個(gè)月的取油量剛好可以滿足第一個(gè)月的需求量,第二個(gè)月購(gòu)滿下十個(gè)月的總油量,所用的儲(chǔ)油費(fèi)用最小。結(jié)合圖以及計(jì)算結(jié)果的分析我們可以肯定我們的模型結(jié)果是正確的,具有一定的可靠性。6.對(duì)問(wèn)題二的求解〔圖論——最短路徑模型〔1確定各種和約油價(jià)表火電廠與售油公司采用遠(yuǎn)期和約來(lái)購(gòu)買柴油,即事先計(jì)算出一個(gè)雙方都認(rèn)為合理的油價(jià)并簽訂一個(gè)和約,將來(lái)購(gòu)油時(shí)就采用這個(gè)油價(jià),這樣的和約期限為一個(gè)月、兩個(gè)月、三個(gè)月、四個(gè)月、六個(gè)月一年共六種?，F(xiàn)在要選擇一種最有利的遠(yuǎn)期和約,使得費(fèi)用最省。我們?cè)O(shè)定雙方都認(rèn)為合理的油價(jià)為合約期內(nèi)個(gè)星期油價(jià)的平均值。對(duì)于第一種遠(yuǎn)期和約,即一個(gè)月為期限,共有十二種情況,具體油價(jià)表如表所示月份123456油價(jià)26302627.552658.22702.752751.232730.2月份789101112油價(jià)2704.22706.22702.32699.227172755.4表.一個(gè)月為期限的和約油價(jià)對(duì)于第二種遠(yuǎn)期和約,即二個(gè)月為期限,共有十一種情況,具體油價(jià)表如表所示表6.1.2.二個(gè)月為期限的和約油價(jià)月份1,22,33,44,55,66,7油價(jià)2628.782643.182680.782726.992740.712717.2月份7,88,99,1010,1111,12油價(jià)2705.22704.262700.762708.12736.18對(duì)于第三種遠(yuǎn)期和約,即三個(gè)月為期限,共有十種情況,具體油價(jià)表如表所示表6.1.3.三個(gè)月為期限的和約油價(jià)月份1,2,32,3,43,4,54,5,65,6,76,7,8油價(jià)2638.782663.032704.262728.062728.542713.53月份7,8,98,9,109,10,1110,11,12油價(jià)2704.242702.582706.182723.85對(duì)于第四種遠(yuǎn)期和約,即四個(gè)月為期限,共有九種情況,具體油價(jià)表如表所示表6.1.4.四個(gè)月為期限的和約油價(jià)月份1,2,3,42,3,4,53,4,5,64,5,6,75,6,7,86,7,8,9油價(jià)2654.782685.082710.742722.092722.962710.73月份7,8,9,108,9,10,119,10,11,12油價(jià)2702.982706.182718.47對(duì)于第五種遠(yuǎn)期和約,即六個(gè)月為期限,共有七種情況,具體油價(jià)表如表所示表6.1.5.六個(gè)月為期限的和約油價(jià)月份1,2,3,4,5,62,3,4,5,6,73,4,5,6,7,84,5,6,7,8,9油價(jià)2683.422695.792708.892716.15月份5,6,7,8,9,106,7,8,9,10,117,8,9,10,11,12油價(jià)2715.562714.042716.02以一年為期限的話,和約油價(jià)為2699.7元。模型的建立與求解按照預(yù)測(cè)出的柴油價(jià)格,我們可以計(jì)算出一年當(dāng)中所有情況下雙方都滿意的購(gòu)油價(jià)格我們把合約期限內(nèi)的柴油價(jià)格與電廠的需求的乘積看做兩個(gè)時(shí)間內(nèi)的距離,利用圖論的知識(shí)建立最短路徑模型。利用matlab編程解答程序見(jiàn)附件1考慮到15天的預(yù)備儲(chǔ)油量,計(jì)算得到十二個(gè)月的采購(gòu)計(jì)劃為：<1>,<2,3>,<4,5,6,7,8,9>,<10>,<11,12>,即年初簽訂一個(gè)月的遠(yuǎn)期購(gòu)油和約,二月初簽訂兩個(gè)月的遠(yuǎn)期購(gòu)油和約,四月初簽訂六個(gè)月的遠(yuǎn)期購(gòu)油和約,十月初簽訂一個(gè)月的遠(yuǎn)期購(gòu)油和約,十一月初簽訂兩個(gè)月的遠(yuǎn)期購(gòu)油和約,共耗費(fèi)資金70445030元。七．對(duì)問(wèn)題三的求解模型的建立與求解利用問(wèn)題一所建立的模型分析首先建立一個(gè)同問(wèn)題一的模型,在編程解答時(shí),只需要對(duì)柴油的儲(chǔ)存費(fèi)用進(jìn)行一定的約束即可,所建模型如下：目標(biāo)函數(shù)：約束條件：用油量的約束條件每月初所擁有的擁有量的約束其中變量約束如果將15天的柴油儲(chǔ)備量用單位"1"來(lái)表示,可得到如下結(jié)果：{0.5,0.6915997},{0.6,0.6919072},{0.7,0.6922146},{0.8,0.6925221},{0.9,0.6928295},{1.0,0.6931369},{1.1,0.6934444},{1.2,0.6937518},{1.3,0.6940592}由此說(shuō)明,最少柴油儲(chǔ)備天數(shù)與所消耗的費(fèi)用是線性關(guān)系,也就是說(shuō)政府規(guī)定的15天的最少柴油儲(chǔ)備是值得懷疑的。這只是定性來(lái)說(shuō)明。利用問(wèn)題二所建模的模型分析根據(jù)問(wèn)題二中所建立的模型,采用遠(yuǎn)期和約來(lái)購(gòu)買柴油的方法,在合約規(guī)定內(nèi),銷售商應(yīng)當(dāng)按照火電廠的需要提供柴油。分別對(duì)庫(kù)存量為5天,10天,12天,13天,14天,15天,20天的情況進(jìn)行分析,得出各種情況下的總費(fèi)用,就可以得到預(yù)備庫(kù)存天數(shù)與費(fèi)用的大致關(guān)系。得到的預(yù)備庫(kù)存天數(shù)與費(fèi)用的關(guān)系：5天10天12天13天14天15天20天70528520704595297040203070354220703631207044503070505020由此看出最少柴油儲(chǔ)備的天數(shù)在從第5天到第15天之間的變化是先減后增的,即最少柴油儲(chǔ)備天數(shù)在15天以下也是可行的,而且可看出使用費(fèi)用最少的儲(chǔ)備天數(shù)是13天。參考文獻(xiàn)朱道元.數(shù)學(xué)建模案例精選.北京：科學(xué)出版社,2003汪國(guó)強(qiáng).數(shù)學(xué)建模優(yōu)秀案例選編〔工科數(shù)學(xué)基地建設(shè)叢書.XX：華南理工大學(xué)出版社,1998盧開(kāi)澄.單目標(biāo)、多目標(biāo)與整數(shù)規(guī)劃.北京：清華出版社,1999錢頌迪.運(yùn)籌學(xué).北京：清華大學(xué)出版社,1990李炯生.數(shù)學(xué)競(jìng)賽中的圖論方法.XX：中國(guó)科學(xué)技術(shù)大學(xué)出版社,1996附件1.圖論模型的matlab程序w=[2629.250,2627.125,2635.417,2650.063,inf,2680.417,inf,inf,inf,inf,inf,2699.895;inf,2624.750,2638.375,2656.917,2679.375,inf,2694.997,inf,inf,inf,inf,inf;inf,inf,2652.150,2673.125,2697.583,2707.063,inf,2704.747,inf,inf,inf,inf;inf,inf,inf,2694.150,2720.375,2725.417,2722.120,inf,2707.080,inf,inf,inf;inf,inf,inf,inf,2746.750,2741.125,2731.493,2726.245,inf,2719.622,inf,inf;inf,inf,inf,inf,inf,2735.520,2723.875,2719.410,2716.558,inf,2715.372,inf;inf,inf,inf,inf,inf,inf,2712.250,2711.375,2710.251,2708.871,inf,2719.375;inf,inf,inf,inf,inf,inf,inf,2710.500,2709.250,2707.750,2711.125,inf;inf,inf,inf,inf,inf,inf,inf,inf,2708.150,2706.375,2711.333,2723.375;inf,inf,inf,inf,inf,inf,inf,inf,inf,2704.750,2713.225,2728.520;inf,inf,inf,inf,inf,inf,inf,inf,inf,inf,2721.250,2740.375;inf,inf,inf,inf,inf,inf,inf,inf,inf,inf,inf,2759.520;]n=size<w,1>;w1=w<1,:>;fori=1:nl<i>=w1<i>;z<i>=1;ends=[];s<1>=1;u=s<1>;k=1lzwhilek<nfori=1:nforj=1:kifi~=s<j>ifl<i>>l<u>+w<u,i>l<i>=l<u>+w<u,i>;z<i>=u;endendendendlzll=l;fori=1:nforj=1:kifi~=s<j>ll<i>=ll<i>;elsell<i>=inf;endendendlv=inf;fori=1:nifll<i><lvlv=ll<i>;v=i;endendlv;v;s<k+1>=v;k=k+1;u=s<k>;endlz運(yùn)行結(jié)果w=1.0e+003*Columns1through112.62922.62712.63542.6501Inf2.6804InfInfInfInfInfInf2.62482.63842.65692.6794Inf2.6950InfInfInfInfInfInf2.65222.67312.69762.7071Inf2.7047InfInfInfInfInfInf2.69422.72042.72542.7221Inf2.7071InfInfInfInfInfInf2.74682.74112.73152.7262Inf2.7196InfInfInfInfInfInf2.73552.72392.71942.7166Inf