具有混合動(dòng)態(tài)約束的系統(tǒng)優(yōu)化調(diào)度

1、最優(yōu)控制課程Project報(bào)告PAGE 具有混合動(dòng)態(tài)約束的系統(tǒng)優(yōu)化調(diào)度火電機(jī)組子問(wèn)題求解第一小組成員： 楊勇 李培 李潔 張俊杰 目 錄第一章 引言問(wèn)題分析實(shí)驗(yàn)結(jié)果心得體會(huì)致謝具有混合動(dòng)態(tài)約束的系統(tǒng)優(yōu)化調(diào)度火電機(jī)組子問(wèn)題求解新方法引言?xún)?yōu)化的問(wèn)題由來(lái)已久。從人類(lèi)勞動(dòng)的歷史開(kāi)始，優(yōu)化的問(wèn)題就應(yīng)運(yùn)而生了。如何在客觀條件一定的情況下，使一件事做的更好，就成為人類(lèi)不懈追求的目標(biāo)。使一件事做的更好可能是花費(fèi)時(shí)間最少，等量產(chǎn)出的情況下消耗最小等等，這就關(guān)系到如何安排一系列事件的順序等問(wèn)題。隨著歷史的發(fā)展，優(yōu)化問(wèn)題的研究也在不斷的進(jìn)步。如今一方面，系統(tǒng)優(yōu)化問(wèn)題滲透到許多學(xué)科研究之中，另一方面，優(yōu)化也成為了一門(mén)

2、專(zhuān)門(mén)的學(xué)科并且有了堅(jiān)實(shí)的理論基礎(chǔ)和一系列的經(jīng)典方法。 生產(chǎn)優(yōu)化調(diào)度的目的是合理安排有限的生產(chǎn)資源，達(dá)到成本最低、能源最省等目標(biāo)。對(duì)于電力、化工等一大類(lèi)生產(chǎn)系統(tǒng)，通常包括多套不同型號(hào)的生產(chǎn)設(shè)備。優(yōu)化調(diào)度問(wèn)題除要考慮生產(chǎn)量、原料限制、合同產(chǎn)量等系統(tǒng)約束外，每個(gè)單臺(tái)設(shè)備通常具有十分復(fù)雜的離散動(dòng)態(tài)運(yùn)行約束如不連續(xù)輸出區(qū)間、開(kāi)關(guān)機(jī)約束，以及連續(xù)動(dòng)態(tài)運(yùn)行約束如生產(chǎn)能力、輸出變化限制等。本文的研究就是在這樣的背景之下進(jìn)行的。主要討論了火電機(jī)組的調(diào)度問(wèn)題。為了在滿(mǎn)足各種系統(tǒng)約束（生產(chǎn)量，合同成本，開(kāi)關(guān)機(jī)約束，爬升約束等等）的前提下是系統(tǒng)的費(fèi)用最低，能源最省，如何調(diào)度不同機(jī)組的開(kāi)關(guān)機(jī)時(shí)間，決定單個(gè)機(jī)組的開(kāi)關(guān)機(jī)序

3、列和發(fā)電功率就是我們的目標(biāo)。問(wèn)題分析1系統(tǒng)方法從數(shù)學(xué)規(guī)劃的角度來(lái)看，整個(gè)火電機(jī)組（24個(gè)機(jī)組）的調(diào)度問(wèn)題屬于NP問(wèn)題，即隨著系統(tǒng)規(guī)模的增大，計(jì)算量將呈現(xiàn)指數(shù)增長(zhǎng)的趨勢(shì)。當(dāng)系統(tǒng)規(guī)模增長(zhǎng)到一定程度時(shí)，直接求解這個(gè)問(wèn)題是不可能的或者幾乎不可能的，例如有時(shí)候，以目前的計(jì)算機(jī)運(yùn)行速度計(jì)算，將要需要幾十萬(wàn)年甚至上百萬(wàn)年數(shù)量級(jí)；即使是幾個(gè)月的時(shí)間也是不現(xiàn)實(shí)的，因?yàn)槿绻菢拥脑?huà)將會(huì)失去實(shí)際意義。近年的研究表明拉格朗日松弛法是求解具有可分結(jié)構(gòu)問(wèn)題的最有效算法之一，其主要思想是在對(duì)偶理論基礎(chǔ)上，利用拉氏乘子松弛系統(tǒng)約束，將原問(wèn)題分解成若干子問(wèn)題，形成二級(jí)優(yōu)化結(jié)構(gòu)。下級(jí)對(duì)子問(wèn)題分別求解，上級(jí)更新拉氏乘子。對(duì)偶解收斂

4、后，再應(yīng)用啟發(fā)式方法或系統(tǒng)化方法將對(duì)偶解調(diào)整為可行解。拉格朗日松弛法的主要優(yōu)點(diǎn)包括：將大規(guī)模問(wèn)題分解為若干子問(wèn)題求解從而降低了復(fù)雜度；便于靈活處理多種約束；計(jì)算時(shí)間隨問(wèn)題規(guī)模增大僅呈線(xiàn)性增長(zhǎng)；計(jì)算中得到的對(duì)偶函數(shù)值可定量評(píng)價(jià)可行調(diào)度方案的優(yōu)劣；Lagrange乘子在實(shí)際中有重要的經(jīng)濟(jì)意義。整個(gè)火電機(jī)組的調(diào)度問(wèn)題可以通過(guò)拉各朗日松弛法分解為多個(gè)子問(wèn)題求解，即各個(gè)單個(gè)機(jī)組的開(kāi)關(guān)機(jī)問(wèn)題。單個(gè)機(jī)組的開(kāi)關(guān)機(jī)問(wèn)題是一個(gè)具有混合參數(shù)（連續(xù)參數(shù)，離散參數(shù)）的最優(yōu)控制問(wèn)題。通過(guò)解耦的方法分別處理離散參數(shù)和連續(xù)參數(shù)，使兩方面的控制都達(dá)到最優(yōu)。首先應(yīng)用動(dòng)態(tài)規(guī)劃法，在每個(gè)連續(xù)運(yùn)行的時(shí)段內(nèi)，確定最優(yōu)的連續(xù)決策變量，然后

5、再利用動(dòng)態(tài)規(guī)劃法，求解最優(yōu)的離散狀態(tài)和控制。鑒于整個(gè)火電機(jī)組生產(chǎn)調(diào)度問(wèn)題復(fù)雜度和工組量的因素（各種約束可能達(dá)到上萬(wàn)個(gè)），我們?cè)谶@里只是嘗試解決單機(jī)組的優(yōu)化調(diào)度問(wèn)題，即下層的子問(wèn)題。那么這個(gè)子系統(tǒng)的輸入是拉各朗日乘子，輸出是這個(gè)機(jī)組的最小消耗成本。在這個(gè)子系統(tǒng)中要考慮發(fā)電量約束，最小開(kāi)關(guān)機(jī)約束，功率的爬升約束，冷開(kāi)機(jī)時(shí)間等多種約束。在這個(gè)子問(wèn)題中，單機(jī)組的最優(yōu)開(kāi)關(guān)機(jī)序列是一個(gè)重要的結(jié)果。是連接輸入?yún)?shù)和輸出參數(shù)的橋梁。2問(wèn)題描述考慮有I套火電機(jī)組的系統(tǒng)，調(diào)度周期為T(mén)個(gè)時(shí)段，那么整個(gè)機(jī)組的混合優(yōu)化問(wèn)題可以表示為： （1）這里和分別是機(jī)組在第時(shí)段內(nèi)的連續(xù)和離散狀態(tài)變量，為表示控制機(jī)組的離散控制變量。

6、通常對(duì)應(yīng)于機(jī)組在第時(shí)段內(nèi)的發(fā)電量，為機(jī)組在第時(shí)段內(nèi)的生產(chǎn)成本，為離散狀態(tài)轉(zhuǎn)換引起的費(fèi)用，一般僅當(dāng)與異號(hào)時(shí)此費(fèi)用非零，否則為零，它與機(jī)組目前的狀態(tài)和離散控制變量有關(guān)。問(wèn)題的約束可分為系統(tǒng)約束與設(shè)備約束兩類(lèi)，系統(tǒng)約束主要包括：(2)(3)其中為第時(shí)段的發(fā)電量要求， 表示電力系統(tǒng)調(diào)度中的系統(tǒng)備用，是設(shè)備在第時(shí)段內(nèi)發(fā)電量為時(shí)消耗的原料（不等號(hào)相反）, 是第時(shí)段的原料限制（或系統(tǒng)備用需求）。設(shè)備約束主要包括：連續(xù)狀態(tài)方程：，(4)其中為連續(xù)控制變量；離散狀態(tài)方程： (5)離散狀態(tài)方程通常表示設(shè)備開(kāi)機(jī)時(shí)間（取正值）或關(guān)機(jī)時(shí)間（取負(fù)值）的演變，比如表示第時(shí)段內(nèi)設(shè)備處于開(kāi)機(jī)狀態(tài)，且已連續(xù)開(kāi)了3個(gè)時(shí)段；連續(xù)狀態(tài)

7、約束：，(6)通常表示機(jī)組的生產(chǎn)能力和最小輸出，其中、分別為設(shè)備的最小和最大輸出量，因?yàn)槲锢砑s束和效率的考慮，機(jī)組的輸出區(qū)間可能是不連續(xù)的，即一般不為0；連續(xù)控制變量約束：，(7)表示機(jī)組輸出變化能力，其中為機(jī)組輸出變化率上限；離散控制變量約束：(8)表示對(duì)機(jī)組開(kāi)關(guān)機(jī)的限制，其中為最小開(kāi)機(jī)時(shí)間，為最小關(guān)機(jī)時(shí)間；3拉格朗日松弛接下來(lái)我們?cè)贚agrange松弛法框架下引入對(duì)偶函數(shù)如下： （9）這樣就可以松弛掉系統(tǒng)約束和產(chǎn)量約束，是對(duì)應(yīng)于系統(tǒng)約束的乘子向量，是對(duì)應(yīng)于產(chǎn)量約束的乘子向量。在此基礎(chǔ)上，連續(xù)控制變量（功率爬升約束）和離散控制變量（開(kāi)機(jī)，關(guān)機(jī)時(shí)間）的對(duì)偶問(wèn)題可描述為： （10）不難看出，(9

8、)可以分解為：其中： (11)稱(chēng)為第i個(gè)子問(wèn)題，它只與機(jī)組i有關(guān)，這樣通過(guò)解對(duì)偶問(wèn)題（11）就可以間接獲得原問(wèn)題的解。求解時(shí)的基本步驟如下：初始化乘子求解所有子問(wèn)題修正乘子，判斷解是否收斂，若未收斂，轉(zhuǎn)b，否則停止。我們的工作就是求解子問(wèn)題。4. 子問(wèn)題解決方案a)離散狀態(tài)變量和離散控制約束的優(yōu)化通過(guò)拉格朗日松弛法，我們得到了火電機(jī)組解耦以后的子問(wèn)題。在解對(duì)偶問(wèn)題時(shí)，由于系統(tǒng)約束已被松弛，所以只需考慮單機(jī)組的約束(4)-(8)。當(dāng)連續(xù)控制變量約束(7)不存在（即）時(shí)，解子問(wèn)題(12)可采用動(dòng)態(tài)規(guī)劃法直接求解。每一個(gè)離散狀態(tài)所對(duì)應(yīng)的最優(yōu)連續(xù)狀態(tài)變量是唯一的，若第時(shí)段，則相應(yīng)的最優(yōu)連續(xù)狀態(tài)變量為：

9、 (13)由于連續(xù)控制變量約束（7）的存在，我們采用按離散控制變量的變化來(lái)劃分階段，得到如下的狀態(tài)轉(zhuǎn)移圖（為描述方便，假定設(shè)備初始時(shí)刻處于關(guān)機(jī)狀態(tài)且已關(guān)夠最小關(guān)機(jī)時(shí)間，其它情況下?tīng)顟B(tài)轉(zhuǎn)移圖的建立過(guò)程完全類(lèi)似）：說(shuō)明：每個(gè)圓圈表示一個(gè)狀態(tài)，圈上的數(shù)字對(duì)應(yīng)一個(gè)時(shí)段編號(hào)。由于初始時(shí)設(shè)備處于關(guān)機(jī)狀態(tài)，所以第一階段的任務(wù)是決定何時(shí)開(kāi)機(jī)，又因?yàn)樵O(shè)備已關(guān)夠最小關(guān)機(jī)時(shí)間，所以最早的可以開(kāi)機(jī)的時(shí)段是第1時(shí)段。既然第1時(shí)段可以作為第1次開(kāi)機(jī)時(shí)段，那么其它時(shí)段都可以作為第1次開(kāi)機(jī)時(shí)段，另外，一種可能的決策是設(shè)備一直不開(kāi)機(jī)，為容納這種決策，在圖中加入一個(gè)時(shí)段，編號(hào)為，因此第1階段的所有決策含義已明確（每個(gè)決策對(duì)應(yīng)一條

10、線(xiàn)段）。接下來(lái)是第2階段，第2階段的任務(wù)是決定第一次關(guān)機(jī)操作在哪個(gè)時(shí)段進(jìn)行。對(duì)于離散狀態(tài)變量和離散控制約束，我們可以用動(dòng)態(tài)規(guī)劃的方法直接求解。動(dòng)態(tài)規(guī)劃的基本思想：動(dòng)態(tài)規(guī)劃是解決多段決策過(guò)程的極為有力的工具。有以下假定：目標(biāo)函數(shù)具有Markov性質(zhì)認(rèn)為現(xiàn)在和將來(lái)的決策不影響過(guò)去的狀態(tài)，決策和目標(biāo)。存在狀態(tài)反饋控制認(rèn)為u(k)=u(x(k),k)，即x(k)可以獲得，不論是直接測(cè)量或是重構(gòu)；根據(jù)x（k），立即作出決策u（k）。用從后往前的逆推方法，即如果某一時(shí)刻以后的控制是最優(yōu)的，那么它以后的任何時(shí)段的控制也都是最優(yōu)的，就可以得到每個(gè)時(shí)段的最優(yōu)控制u，這就是Bellman方程： 定義：Bellma

11、n方程：b)連續(xù)控制變量的優(yōu)化在計(jì)算開(kāi)機(jī)區(qū)間決策費(fèi)用時(shí)，動(dòng)態(tài)規(guī)劃的方法是一種比較有效的方法，通過(guò)求得各時(shí)段的出發(fā)費(fèi)用(Cost-to-go)函數(shù)來(lái)得到最優(yōu)的連續(xù)狀態(tài)和控制。已有文獻(xiàn)在解決此問(wèn)題時(shí)基本思想是結(jié)構(gòu)動(dòng)態(tài)規(guī)劃(6,10,14)，但由于未能充分利用問(wèn)題(1)-(8)的特性，計(jì)算過(guò)程比較復(fù)雜。通過(guò)分析連續(xù)動(dòng)態(tài)約束特性，在這里，我們使用了一種簡(jiǎn)化結(jié)構(gòu)動(dòng)態(tài)規(guī)劃新算法，它的基本思想是通過(guò)分段線(xiàn)性函數(shù)的平移和相加運(yùn)算來(lái)求得各時(shí)段的出發(fā)費(fèi)用。方法計(jì)算過(guò)程簡(jiǎn)單且計(jì)算量非常小。假定，是兩個(gè)時(shí)段編號(hào)且，再假定對(duì)應(yīng)于圖2中一個(gè)開(kāi)機(jī)區(qū)間，由(12)式知該決策的費(fèi)用為離散狀態(tài)轉(zhuǎn)換費(fèi)用(相應(yīng)于開(kāi)機(jī)區(qū)間之后的關(guān)機(jī)操

12、作的費(fèi)用，此項(xiàng)費(fèi)用容易計(jì)算)與下述連續(xù)優(yōu)化問(wèn)題的最優(yōu)值之和：(14)其中所受約束包括：(15)(16)(17) (14)-(17)是一個(gè)非線(xiàn)性規(guī)劃問(wèn)題，且只含連續(xù)變量。在很多實(shí)際問(wèn)題中，設(shè)備性能曲線(xiàn)通常都是以離散點(diǎn)方式給出，是連續(xù)分段線(xiàn)性函數(shù)。應(yīng)用簡(jiǎn)化結(jié)構(gòu)動(dòng)態(tài)規(guī)劃法對(duì)(14)-(17)求解。首先遞歸定義出發(fā)費(fèi)用函數(shù)：(18)上式的極小化受(16)、(17)約束?？紤]到（15），記(19)于是(20)由此可得問(wèn)題(14)-(17)的動(dòng)態(tài)規(guī)劃遞推算法如下：1反向遞推：初始化，置。對(duì)，由(20)求得。2正向遞推：初始化，置：，其中受 (17)約束。對(duì)，依次計(jì)算：在(19)中置為，求得(19)的最優(yōu)解

13、為，置。根據(jù)動(dòng)態(tài)規(guī)劃原理，上述算法中得到的、就是(14)-(17)的最優(yōu)解。在上述算法的第一步，要求得到出發(fā)費(fèi)用函數(shù)的表達(dá)式，這是結(jié)構(gòu)動(dòng)態(tài)規(guī)劃方法與一般動(dòng)態(tài)動(dòng)態(tài)規(guī)劃方法的最大區(qū)別。若，均為凸的分段線(xiàn)性函數(shù)，問(wèn)題可以大大簡(jiǎn)化。此時(shí)可證也是凸分段線(xiàn)性函數(shù)(6,10,12)，所以只需要確定轉(zhuǎn)折點(diǎn)，就可以確定的表達(dá)式。考慮到(14)-(17)的特點(diǎn)，提出簡(jiǎn)化結(jié)構(gòu)動(dòng)態(tài)規(guī)劃法計(jì)算，使得第二步的計(jì)算非常簡(jiǎn)單。以下是針對(duì)g凸函數(shù)的情況下所進(jìn)行的相應(yīng)的平移操作，具體證明從略。其中，情況一是g為單調(diào)增的情況，情況2是g為單調(diào)減的情況，情況3為先減后增的情況。圖3. 與的關(guān)系下面我們以求解開(kāi)機(jī)區(qū)間23，24的出發(fā)費(fèi)

14、用為例，說(shuō)明上述算法及其物理意義。首先進(jìn)行反向遞推：給出24時(shí)刻所對(duì)應(yīng)的函數(shù)（，在不影響問(wèn)題本質(zhì)的前提下，我們簡(jiǎn)化了f的表達(dá)式，忽略了乘子u），并找出g的極小值點(diǎn)，圖中是x130。從最低點(diǎn)處，平移爬升約束個(gè)單位，如下圖所示：將其與23時(shí)刻的g函數(shù)相加： 得到23時(shí)刻的出發(fā)費(fèi)用函數(shù)如下圖所示：類(lèi)似的，如果開(kāi)機(jī)區(qū)間為22，23，24，求22時(shí)刻的出發(fā)費(fèi)用，則需將23時(shí)刻的出發(fā)費(fèi)用函數(shù)進(jìn)行響應(yīng)平移，如下圖所示：然后再按照遞推公式，與22時(shí)刻的f函數(shù)相加，即可得到22時(shí)刻的出發(fā)費(fèi)用。其他開(kāi)機(jī)區(qū)段的出發(fā)費(fèi)用函數(shù)類(lèi)似。最后，根據(jù)正向遞推，求出最優(yōu)的連續(xù)狀態(tài)和控制。首先求出令出發(fā)費(fèi)用最低的（即：求出出發(fā)費(fèi)用

15、函數(shù)的最小值點(diǎn)）；根據(jù)下一時(shí)刻的出發(fā)費(fèi)用函數(shù)，決定出下一個(gè)最優(yōu)狀態(tài)；由和的出；循環(huán)進(jìn)行b)、c),直到求出最后一個(gè)開(kāi)機(jī)時(shí)段的最優(yōu)狀態(tài)。第三章 實(shí)驗(yàn)結(jié)果1實(shí)驗(yàn)（1）所用數(shù)據(jù)pmax=130; 最大功率pmin=20; 最小功率tminup=5; 最小開(kāi)機(jī)時(shí)間tmindown=5; 最小關(guān)機(jī)十幾年tcoldup=9; 冷關(guān)機(jī)時(shí)間delta=40; 爬升約束hotstart=500; 熱啟動(dòng)費(fèi)用coldstart=2*hotstart; 冷啟動(dòng)費(fèi)用l=12,23,17,17,36,17.5,18,27,15,6,8,16,22,50,10,20,37,16,8,26,17,12,22,30; 乘子向

16、量step=(pmax-pmin)/22;point=pmin:step:pmax;costofpoint=0.002*point.*point+16.6*point+700; 費(fèi)用函數(shù)2實(shí)驗(yàn)所求得的最優(yōu)狀態(tài)及最小費(fèi)用開(kāi)機(jī)序列（0：關(guān)機(jī)態(tài) 1：開(kāi)機(jī)態(tài)）a = Columns 1 through 12 0 0 0 1 1 1 1 1 0 0 0 0 Columns 13 through 24 0 1 1 1 1 1 0 0 0 0 0 1最優(yōu)狀態(tài)序列：e = Columns 1 through 12 0 0 0 100 130 130 130 130 0 0 0 0 Columns 13 thr

17、ough 24 0 130 90 130 130 90 0 0 0 0 0 130調(diào)度得到的最小值mincost = -3.6202e+003最優(yōu)狀態(tài)的圖形表示如下：3用遍歷法得到的結(jié)果為了與動(dòng)態(tài)規(guī)劃的方法進(jìn)行比較，我們還采用了遍歷的方法用遍歷方法得到的結(jié)果：開(kāi)機(jī)序列（0：關(guān)機(jī)態(tài) 1：開(kāi)機(jī)態(tài)）a = Columns 1 through 12 0 0 0 1 1 1 1 1 0 0 0 0 Columns 13 through 24 0 1 1 1 1 1 0 0 0 0 0 1最優(yōu)狀態(tài)序列：e = Columns 1 through 12 0 0 0 100 130 130 130 130 0

18、 0 0 0 Columns 13 through 24 0 130 90 130 130 90 0 0 0 0 0 130調(diào)度得到的最小值mincost = -3.6202e+003開(kāi)機(jī)序列（0：關(guān)機(jī)態(tài) 1：開(kāi)機(jī)態(tài)）a = Columns 1 through 12 0 0 0 1 1 1 1 1 0 0 0 0 Columns 13 through 24 0 1 1 1 1 1 0 0 0 0 0 1最優(yōu)狀態(tài)序列：e = Columns 1 through 12 0 0 0 100 130 130 130 130 0 0 0 0 Columns 13 through 24 0 130 90

19、130 130 90 0 0 0 0 0 130調(diào)度得到的最小值mincost = -3.6202e+003遍歷結(jié)果的圖形輸出如下：1實(shí)驗(yàn)（2）所用數(shù)據(jù)pmax=130; 最大功率pmin=20; 最小功率tminup=5; 最小開(kāi)機(jī)時(shí)間tmindown=5; 最小關(guān)機(jī)十幾年tcoldup=9; 冷關(guān)機(jī)時(shí)間delta=40; 爬升約束hotstart=560; 熱啟動(dòng)費(fèi)用coldstart=2*hotstart; 冷啟動(dòng)費(fèi)用l=12,23,17,17,36,17.5,18,27,15,6,8,16,22,50,10,20,37,16,8,26,17,12,22,30; 乘子向量step=(pm

20、ax-pmin)/22;point=pmin:step:pmax;costofpoint=0.00211*point.*point+16.5*point+680; 費(fèi)用函數(shù)2實(shí)驗(yàn)所求得的最優(yōu)狀態(tài)及最小費(fèi)用開(kāi)機(jī)序列（0：關(guān)機(jī)態(tài) 1：開(kāi)機(jī)態(tài)）a = Columns 1 through 12 0 0 0 1 1 1 1 1 0 0 0 0 Columns 13 through 24 0 1 1 1 1 1 0 0 0 0 0 1最優(yōu)狀態(tài)序列：e = Columns 1 through 12 0 0 0 120 130 130 130 130 0 0 0 0 Columns 13 through 24

21、 0 130 90 130 130 90 0 0 0 0 0 130調(diào)度得到的最小值mincost = -3.9552e+003最優(yōu)狀態(tài)的圖形表示如下：3用遍歷法得到的結(jié)果為了與動(dòng)態(tài)規(guī)劃的方法進(jìn)行比較，我們還采用了遍歷的方法用遍歷方法得到的結(jié)果：開(kāi)機(jī)序列（0：關(guān)機(jī)態(tài) 1：開(kāi)機(jī)態(tài)）a = Columns 1 through 12 0 0 0 1 1 1 1 1 0 0 0 0 Columns 13 through 24 0 1 1 1 1 1 0 0 0 0 0 1最優(yōu)狀態(tài)序列：e = Columns 1 through 12 0 0 0 120 130 130 130 130 0 0 0 0

22、Columns 13 through 24 0 130 90 130 130 90 0 0 0 0 0 130調(diào)度得到的最小值mincost = -3.9552e+003開(kāi)機(jī)序列（0：關(guān)機(jī)態(tài) 1：開(kāi)機(jī)態(tài)）a = Columns 1 through 12 0 0 0 1 1 1 1 1 0 0 0 0 Columns 13 through 24 0 1 1 1 1 1 0 0 0 0 0 1最優(yōu)狀態(tài)序列：e = Columns 1 through 12 0 0 0 100 130 130 130 130 0 0 0 0 Columns 13 through 24 0 130 90 130 130 90 0 0 0 0 0 130調(diào)度得到的最小值mincost = -3.6202e+003遍歷結(jié)果的圖形輸出如下：4兩種方法的比較