電力市場的輸電阻塞模型 帶matlab程序

1、電力市場的輸電阻塞模型 摘要: 本文針對第一問，建立了多元線性均值回歸模型，用于描述各線路上的有功潮流與各發(fā)電機組出力的關系。模型中利用均值定理，運用matlab軟件求得解，此模型的誤差較小，較符合實際情況。 對于第二問，我們引入基本補償價概念，利用限上和限下的方法設計阻塞費用計算規(guī)則，在阻塞費用計算規(guī)則里面我們認為序外容量在低于對應報價的清算價上出力的時候，網(wǎng)方已經(jīng)按清算價給予了發(fā)電商一部分的費用，這部分費用不作為經(jīng)濟補償，即其是不包含在阻塞費用中，對于序外容量的經(jīng)濟補償我們就用對應段容量的報價和清算價之差作為補償價格來計算阻塞費用。 對于第三問，在考慮爬坡速率和負荷需求的前提下，根據(jù)電力市

2、場規(guī)則，考慮最小購電成本為目標函數(shù)，利用最優(yōu)化的方法，我們得到了各機組的出力分配預案（如下表所示）。機組12345678出力1507918099.512514095113.9此時各機組所對應的潮流值（如下表所示）線路123456潮流值173.315796.8095-150.5051101.8151136.8380168.5319對于第四問，利用第一問求得的結果，檢驗第三問各機組的出力分配預案是否會引起輸電阻塞，如果引起輸電阻塞，根據(jù)輸電阻塞管理原則并考慮經(jīng)濟原則，用計算機編程搜索調(diào)整預案，得出阻塞費用為405.8004元，新的方案見下表機組12345678出力1538822890.177515

3、29560.1116.1225此時各線路所對應的潮流值（如下表所示）線路123456潮流值165100.6408-154.8553105.5524131.5169159.5725 對于第五問，重復三，四問的工作，我們得出方案（如下表所示），阻塞費用為94.0377元。機組12345678出力1538822899.515215560.3117此時各線路對應的潮流值（如下表所示）線路123456潮流值173.461994.1143-154.9968104.9631135.3452160.4779問題的重述：電網(wǎng)公司在組織交易，調(diào)度和配送時,必須遵循電網(wǎng)“安全第一”的原則，同時制定一個電力市場交易規(guī)

4、則，按照購電費用最小的經(jīng)濟目標來運作。市場交易-調(diào)度中心根據(jù)負荷預報和交易規(guī)則制定滿足電網(wǎng)安全運行的調(diào)度計劃各發(fā)電機組的處理(發(fā)電功率)分配方案；在執(zhí)行調(diào)度計劃的過程中，還需實時調(diào)度承擔AGC(自動發(fā)電控制)輔助服務的機組出力，以跟蹤電網(wǎng)中實時變化的負荷。設某電網(wǎng)有若干臺發(fā)電機組和若干條主要線路，每條線路上的有功潮流（輸電功率和方向）取決于電網(wǎng)結構和各發(fā)電機組的出力。電網(wǎng)每條線路上的有功潮流的絕對值有一安全限值，限值還具有一定的相對安全裕度（即在應急情況下潮流絕對值可以超過限值的百分比的上限）。如果各機組出力分配方案使某條線路上的有功潮流的絕對值超出限值，稱為輸電阻塞。當發(fā)生輸電阻塞時，需要研

5、究如何制訂既安全又經(jīng)濟的調(diào)度計劃。 電力市場交易規(guī)則：1. 以15分鐘為一個時段組織交易，每臺機組在當前時段開始時刻前給出下一個時段的報價。各機組將可用出力由低到高分成至多10段報價，每個段的長度稱為段容量，每個段容量報一個價（稱為段價），段價按段序數(shù)單調(diào)不減。在最低技術出力以下的報價一般為負值，表示愿意付費維持發(fā)電以避免停機帶來更大的損失。2. 在當前時段內(nèi)，市場交易-調(diào)度中心根據(jù)下一個時段的負荷預報,每臺機組的報價、當前出力和出力改變速率，按段價從低到高選取各機組的段容量或其部分（見下面注釋），直到它們之和等于預報的負荷，這時每個機組被選入的段容量或其部分之和形成該時段該機組的出力分配預案

6、（初始交易結果）。最后一個被選入的段價（最高段價）稱為該時段的清算價，該時段全部機組的所有出力均按清算價結算。注釋：（a） 每個時段的負荷預報和機組出力分配計劃的參照時刻均為該時段結束時刻。（b） 機組當前出力是對機組在當前時段結束時刻實際出力的預測值。（c） 假設每臺機組單位時間內(nèi)能增加或減少的出力相同，該出力值稱為該機組的爬坡速率。由于機組爬坡速率的約束，可能導致選取它的某個段容量的部分。（d） 為了使得各機組計劃出力之和等于預報的負荷需求，清算價對應的段容量可能只選取部分。市場交易-調(diào)度中心在當前時段內(nèi)要完成的具體操作過程如下：1、 監(jiān)控當前時段各機組出力分配方案的執(zhí)行，調(diào)度AGC輔助服

7、務，在此基礎上給出各機組的當前出力值。2、 作出下一個時段的負荷需求預報。3、 根據(jù)電力市場交易規(guī)則得到下一個時段各機組出力分配預案。4、 計算當執(zhí)行各機組出力分配預案時電網(wǎng)各主要線路上的有功潮流，判斷是否會出現(xiàn)輸電阻塞。如果不出現(xiàn)，接受各機組出力分配預案；否則，按照如下原則實施阻塞管理：輸電阻塞管理原則：（1） 調(diào)整各機組出力分配方案使得輸電阻塞消除。（2） 如果（1）做不到，還可以使用線路的安全裕度輸電，以避免拉閘限電（強制減少負荷需求），但要使每條線路上潮流的絕對值超過限值的百分比盡量小。（3） 如果無論怎樣分配機組出力都無法使每條線路上的潮流絕對值超過限值的百分比小于相對安全裕度，則必

8、須在用電側拉閘限電。（4） 當改變根據(jù)電力市場交易規(guī)則得到的各機組出力分配預案時，一些通過競價取得發(fā)電權的發(fā)電容量（稱序內(nèi)容量）不能出力；而一些在競價中未取得發(fā)電權的發(fā)電容量（稱序外容量）要在低于對應報價的清算價上出力。因此，發(fā)電商和網(wǎng)方將產(chǎn)生經(jīng)濟利益沖突。網(wǎng)方應該為因輸電阻塞而不能執(zhí)行初始交易結果付出代價，網(wǎng)方在結算時應該適當?shù)亟o發(fā)電商以經(jīng)濟補償，由此引起的費用稱之為阻塞費用。網(wǎng)方在電網(wǎng)安全運行的保證下應當同時考慮盡量減少阻塞費用?，F(xiàn)要解決的問題如下：（表格見附錄1）1. 某電網(wǎng)有8臺發(fā)電機組，6條主要線路，表1和表2中的方案0給出了各機組的當前出力和各線路上對應的有功潮流值，方案132給出

9、了圍繞方案0的一些實驗數(shù)據(jù)，試用這些數(shù)據(jù)確定各線路上有功潮流關于各發(fā)電機組出力的近似表達式。2. 設計一種簡明、合理的阻塞費用計算規(guī)則，除考慮上述電力市場規(guī)則外，還需注意：在輸電阻塞發(fā)生時公平地對待序內(nèi)容量不能出力的部分和報價高于清算價的序外容量出力的部分。3. 假設下一個時段預報的負荷需求是982.4MW，表3、表4和表5分別給出了各機組的段容量、段價和爬坡速率的數(shù)據(jù)，試按照電力市場規(guī)則給出下一個時段各機組的出力分配預案。4. 按照表6給出的潮流限值，檢查得到的出力分配預案是否會引起輸電阻塞，并在發(fā)生輸電阻塞時，根據(jù)安全且經(jīng)濟的原則，調(diào)整各機組出力分配方案，并給出與該方案相應的阻塞費用。5.

10、 假設下一個時段預報的負荷需求是1052.8MW，重復34的工作。問題的分析：1 初步分析題意，根據(jù)表1和表2所給的數(shù)據(jù)，我們利用數(shù)理統(tǒng)計知識和線性回歸分析方法，擬建立多元線性均值回歸模型，推算出各線路上有功潮流關于各發(fā)電機組出力的近似表達式；2 設計簡明，合理的阻塞費用的計算規(guī)則，關鍵在于體現(xiàn)網(wǎng)方是怎樣充分補償發(fā)電廠的損失的，在這里，我們引入基本補償價的概念，利用限上和限下的方法得出阻塞費用的計算規(guī)則；3 僅僅考慮電力市場交易規(guī)則和各機組的爬坡速率的影響，要得出下一時段各機組的出力分配預案，就是在各機組的總段容量不超出爬坡速率所限制的極限值的前提條件下，選取價格最低的全部機組的段容量總和，使

11、得這個總和等于負荷需求； 4 用第一問所得出來的有功潮流的表達式檢查預案是否會引起輸電阻塞，要是會發(fā)生輸電阻塞，則根據(jù)輸電阻塞管理原則并考慮電力市場規(guī)則的情況下調(diào)整預案得出一個新的方案，并根據(jù)設計的阻塞費用計算規(guī)則求出新方案對于預案的阻塞費用。模型的假設：1 各線路上的有功潮流只受到這個電網(wǎng)里面八個發(fā)電機組所發(fā)的電量的影響，與外界的電網(wǎng)無關，即電網(wǎng)結構不會改變；2 發(fā)電機組出力是按照爬坡速率一直增加或減少到符合出力分配方案的要求后，出力就是一直恒定不變的；3 發(fā)電機出力的增量是由低價格的容量增加到高價格的容量的，而出力的減量就是由高價格的容量減少到低價格的容量；符號的說明： ： 第i條線路的潮

12、流值（i1，2，6）； ： 序外容量要出力的段容量或其部分； ： 序內(nèi)容量不能出力的段容量或其部分； ： 第i個機組的出力（i1，2，8）； ： 當前方案0中第i個機組的出力（i1，2，8）； ： 第i個機組在單位時間內(nèi)的功率增量（i1，2，8）； ： 第i個機組的爬坡速率（i1，2，8）； ： 一個時間段（15分鐘）； ： 對應段容量/對應段容量所在機組的爬坡速率； ： 對應段容量的前面容量總和變化到方案0的出力所需的時間；模型的建立：問題1 多元線性均值回歸模型 根據(jù)題中所給統(tǒng)計數(shù)據(jù)，我們用回歸分析方法，考慮線性分析模型，討論得到：運用全回歸分析雖然計算簡單，但是最終的計算結果誤差比較大。

13、所以我們建立多元線性均值回歸模型。則有： （1）對（1）式的算法如下： （2） （3） 由表1和表2的132方案中知，為稀疏矩陣，元素只有一項i不為零（此項機組表現(xiàn)為相對0方案的變動），且同類的有4組，剛好8個機組的單獨變動都有4組數(shù)，很明顯我們可以均值就可求出i組的k的列向量，運用循環(huán)即可全部求得8組k的列向量，合并得k值。（程序附錄2和3）系數(shù)矩陣k觀察以上求出的系數(shù)k可以發(fā)現(xiàn)：k值比較小，而且相鄰k值的變化也比較小。 求c時（程序見附錄4），利用以上求得k值，有： （4） 與對應于132方案中的某一方案，由此運用循環(huán)當求36組，求均值可得。常數(shù)c所以寫出y和p的近似表達式為：對上述方法進

14、行誤差分析，由求出的表達式返回32組方案與原方案對比，得出y的相對誤差圖（如圖一所示）： 如圖一可知，此模型的誤差較小。問題2對于第二問，考慮的是當阻塞發(fā)生的時候網(wǎng)方所要補償給發(fā)電商的費用，這個費用稱為阻塞費用S，一共包括兩個部分：一部分是序內(nèi)容量不能出力的部分，決定這部分阻塞費用的主要因素是不能出力的機組停發(fā)或少發(fā)電的度數(shù)以及少發(fā)一度電所需的補償價格；另一部分是報價高于清算價的序外容量出力的部分，決定這部分阻塞費用的主要因素是序外容量機組多發(fā)電的度數(shù)以及多發(fā)一度電所需的補償價格。模型一：限上和限下的方法根據(jù)上面的思想，我們應用限上和限下的方法1來計算阻塞費用，規(guī)則如下：1 當出現(xiàn)輸電阻塞時，

15、在競價中未取得發(fā)電權的發(fā)電容量（序外容量）按報價由低到高的順序安排機組發(fā)電(限上),在競價中本已經(jīng)取得發(fā)電權的發(fā)電容量（序內(nèi)容量）按報價由高到低的順序安排機組停發(fā)或少發(fā)電(限下)。2 在競價中未取得發(fā)電權的發(fā)電容量（序外容量）是在低于其對應報價的清算 價上出力的，這樣就會使得發(fā)電商利益虧損而網(wǎng)方從中盈利，網(wǎng)方應該對這 些發(fā)電商給與經(jīng)濟補償，這里補償我們就用對應段容量的報價和清算價之差作為補償價格；而那些在競價中本已經(jīng)取得發(fā)電權的發(fā)電容量（序內(nèi)容量）卻不能出力的發(fā)電商，由于不能出力或部分不能出力，使得這些發(fā)電商由于網(wǎng)方線路原因造成經(jīng)濟利益的損失，這時網(wǎng)方也應該給與這些發(fā)電商經(jīng)濟補償，這里我們就以

16、清算價和其不能出力的對應的段容量的報價之差作為補償價格。根據(jù)計算規(guī)則，我們可以列出阻塞費用S的式子： （5）其中 清算價對應段容量的報價對應段容量的報價清算價模型二：模型二是對于模型一所進行的改進，在模型一中可以發(fā)現(xiàn)：對于序內(nèi)容量，當某個的段容量的報價等于它的清算價的時候，此時0，這表明這個段容量是得不到經(jīng)濟補償?shù)?。由于每一個發(fā)電廠發(fā)的每度電都存在一個利潤空間，這個利潤是包含在報價里面的，按照模型一的結論表明發(fā)電廠按清算價作為報價的段容量的利潤得不到經(jīng)濟補償，顯然這是不合理的。為了解決這個問題，我們在這里引入一個概念：基本補償價B，根據(jù)文獻資料 2 ，我們定 （6），于是我們就以清算價和其不能

17、出力的對應的段容量的報價之差再加上基本補償價作為序內(nèi)容量不能出力的發(fā)電商的經(jīng)濟補償。于是式子（5）變?yōu)椋?（7）問題3由題意可知: （8）在下一個時間段的結束時刻各機組的出力大小為 （9）根據(jù)式子（8）（9），結合表1中方案0各機組的數(shù)據(jù)和表5中各機組的爬坡速率，相對于方案0，我們可以算出在下一個時間段的結束時刻各機組所允許能選取的出力范圍（如下表31所示）表31 各機組所允許能選取的出力范圍機組12345678最大出力1538822899.5152155102.1117最少出力875813260.5989560.163根據(jù)電力市場規(guī)則，我們對附錄1中的表3，表4進行處理。（按照段價由低到高排

18、列如下表格群所示）表格群32 價格-800-610-607-590-560-505-500-500機組83651174段容量70110957570705055價格00000116120124機組12356571段容量00050151550價格146150152153159168170173機組54386146段容量054001001020價格180182183188189200203205機組72853426段容量520201501080價格210215233233245250251252機組15382571段容量0153001501530價格252253255258260283300302機

19、組68437824段容量1520100100610價格303305306308310312315318機組86735178段容量20101020100510價格320325330335348356360363機組24177321段容量21501034000價格380380396400405410415435機組64586234段容量20010150000價格489495500510520548800800機組12356784段容量408401010251按照以上表格，僅僅考慮電力市場規(guī)則，而不考慮各機組所允許的出力范圍的時候，選取到第7個表格中的303價格位的段容量就能夠滿足負荷需求（這個價格

20、位之前的段容量全部選取），這時候我們得到各個機組的出力分別為：表33機組12345678出力1507918010012514095113.4由于受到爬坡速率的限制，第四機組的出力最大值為99.5，不能達到100，將第四組不能達到的0.5段容量按照電力市場規(guī)則選取還未選取的且對應段容量價格最低的，由此，我們得出下一個時段各機組的出力分配預案為（如下表34所示）：表34 下一時段各機組的出力分配預案（負荷需求為982.4MW時）機組12345678出力1507918099.512514095113.9問題4 將上題所得的各機組的出力分配預案代入第一問的模型二中求出的y和p的關系式中，得到各線路的潮

21、流值如下表41所示：表41 各線路的潮流值線路123456潮流值173.315796.8095-150.5051101.8151136.8380168.5319上表41與附錄1中的表6相比較可以發(fā)現(xiàn)：線路1，5，6會引起輸電阻塞。 圖 二根據(jù)算法（如圖二所示）編寫程序（見附錄6），由計算機搜索得到：改變后的出力分配方案如下表42所示：表42 各機組的出力分配方案機組12345678出力15388228901116.1225此時各線路的潮流值如下表43所示：表43 改變方案后各線路的潮流值線路123456潮流值165100.6408-154.8553105.552413

22、1.5169159.5725并得到對應的阻塞費用為405.8004元。問題5 利用第三問和第四問同樣的方法，得出在下一個時段預報的負荷需求是1052.8MW時，各機組的出力分配預案如下（如下表51所示）表51 下一時段各機組的出力分配預案（負荷需求為1052.8MW時）機組12345678出力15081218.299.5135150102.1117此時各線路的潮流值如下表52所示表52 各線路的潮流值線路123456潮流值177.260894.6307-155.7527110.0612134.8521167.0870調(diào)整之后各機組的出力分配方案為：（如表53所示） 表53 調(diào)整后的各機組出力分

23、配方案機組12345678出力1538822899.515215560.3117各線路的潮流值為：（如表54所示） 表54 調(diào)整后的各線路的潮流值線路123456潮流值173.461994.1143-154.9968104.9631135.3452160.4779從表54可以知道，第1，5條線路都是超過限值的，但是都在安全裕度內(nèi)，第一條線路超出限值的百分比為5.128，第五條線路超出限值的百分比為2.5342。并得出相應的阻塞費用為：94.0377元。模型的優(yōu)缺點和改進方向：優(yōu)缺點：本論文中的平均值分析法能在一定范圍內(nèi)很好地吻合數(shù)據(jù)，經(jīng)誤差分析亦可得相對誤差0.0001左右，而且機組系數(shù)比較小

24、，能很好地說明此電網(wǎng)有比較好的穩(wěn)定性，不會因一個機組的小變化而有大的改動，但是，此模型只適應機組出力都較大的情況下，在所有機組都比較小時，此模型有較大誤差，但此范圍并不在此題的計算范圍內(nèi)，此模型已經(jīng)能非常很好地解決一般的電網(wǎng)調(diào)動。另外，本文中的計算較多應用matlab編程，對大量的數(shù)據(jù)計算提供很好的橋梁，針對每一個問題都有一個程序來解決，為以后的利用提供方便。橋梁改進方向：由于給出數(shù)據(jù)中有突變數(shù)據(jù)，所以在計算得出的相對誤差中有6個是大于0.05的，當踢除這些值代以估算的均值，能更完美地完成數(shù)據(jù)的擬合。對于賠償算法中有發(fā)電商的隱藏利潤，此利潤應該賠償過加入更加合理準確，對此中的專業(yè)估計，是改進的

25、重點。參考文獻：1 柯進,管霖 電力市場下的輸電阻塞管理技術 電力系統(tǒng)自動化 第26卷 第14期 第四頁 2002年2 西部水利水電投資項目介紹29/xbshj/xbbshd/xbbshd_indx.htm 訪問時間：2004-9-19 韓忠誠 呂伯雄 金昌市招標項目 訪問時間：2004-9-19附錄1：表1 各機組出力方案 （單位：兆瓦，記作MW）方案機組123 4 5 6780120731808012512581.1901133.02731808012512581.1902129.63731808012512581.1903158.7773180801

26、2512581.1904145.32731808012512581.190512078.5961808012512581.190612075.451808012512581.190712090.4871808012512581.190812083.8481808012512581.190912073231.398012512581.1901012073198.488012512581.1901112073212.648012512581.1901212073190.558012512581.190131207318075.85712512581.190141207318065.95812512

27、581.190151207318087.25812512581.190161207318097.82412512581.190171207318080150.7112581.190181207318080141.5812581.190191207318080132.3712581.190201207318080156.9312581.190211207318080125138.8881.190221207318080125131.2181.190231207318080125141.7181.190241207318080125149.2981.19025120731808012512560.

28、5829026120731808012512570.9629027120731808012512564.8549028120731808012512575.5299029120731808012512581.1104.8430120731808012512581.1111.2231120731808012512581.198.09232120731808012512581.1120.44表2 各線路的潮流值（各方案與表1相對應，單位：MW）方案線路1234560164.78140.87-144.25119.09135.44157.691165.81140.13-145.14118.63135.

29、37160.762165.51140.25-144.92118.7135.33159.983167.93138.71-146.91117.72135.41166.814166.79139.45-145.92118.13135.41163.645164.94141.5-143.84118.43136.72157.226164.8141.13-144.07118.82136.02157.57165.59143.03-143.16117.24139.66156.598165.21142.28-143.49117.96137.98156.969167.43140.82-152.26129.58132.

30、04153.610165.71140.82-147.08122.85134.21156.2311166.45140.82-149.33125.75133.28155.0912165.23140.85-145.82121.16134.75156.7713164.23140.73-144.18119.12135.57157.214163.04140.34-144.03119.31135.97156.3115165.54141.1-144.32118.84135.06158.2616166.88141.4-144.34118.67134.67159.2817164.07143.03-140.9711

31、8.75133.75158.8318164.27142.29-142.15118.85134.27158.3719164.57141.44-143.3119134.88158.0120163.89143.61-140.25118.64133.28159.1221166.35139.29-144.2119.1136.33157.5922165.54140.14-144.19119.09135.81157.6723166.75138.95-144.17119.15136.55157.5924167.69138.07-144.14119.19137.11157.6525162.21141.21-14

32、4.13116.03135.5154.2626163.54141-144.16117.56135.44155.9327162.7141.14-144.21116.74135.4154.8828164.06140.94-144.18118.24135.4156.6829164.66142.27-147.2120.21135.28157.6530164.7142.94-148.45120.68135.16157.6331164.67141.56-145.88119.68135.29157.6132164.69143.84-150.34121.34135.12157.64表3 各機組的段容量 （單位

33、：MW）機組段1234567891017005000300004023002081562008311004003002040040455510101010150015755150151501010106950102001510200107501551510105103287002002002010155表4 各機組的段價（單位：元/兆瓦小時，記作元/MWh）機組段123456789101-50501241682102523123303634892-56001822032453003203604104953-61001521892332583083564155004-50015017020025

34、53023253804358005-59001161461882152503103965106-60701591732052523053804055207-5001201802512603063153353485488-800153183233253283303318400800表5 各機組的爬坡速率 （單位：MW/分鐘）機組12345678速率1.821.41.8表6 各線路的潮流限值（單位：MW）和相對安全裕度線路123456限值165150160155132162安全裕度13%18%9%11%15%14%附錄2function k=zhen(fa8,r6) %求k初值

35、 p0=164.78 140.87 -144.25 119.09 135.44 157.69;pj=120 73 180 80 125 125 81.1 90;j=0;for i=1:8 j=fa8(i)-pj(i)+j;endpi=r6-p0;k=pi./j;附錄3function ok=zhen1(fa8,r6) %求k系數(shù)clcfor i=1:32 a=fa8(i,:);b=r6(i,:); k(i,:)=zhen(a,b);endoptions=;for i=1:8 for j=1:6 ok(i,j)=mean(k(4*i,j),k(4*i-1,j),k(4*i-2,j),k(4*i-

36、3,j); end end附錄4function c=zhen4(ok,fa8,r6) %求c值for k=1:32 if k=11 c(k,:)=c(k-1,:); end for j=1:6 c(k,j)=r6(k,j)-fa8(k,:)*ok(:,j); endendc=mean(c);附錄5function pout=zhen5(pin) %求潮流值ok=0.08284 0.04828 0.05297 0.11993 -0.02544 0.12201 0.12158 -0.00123;. -0.05456 0.12785 -0.00003 0.03328 -0.08685 -0.112

37、44 -0.01893 -0.09873;. -0.06954 0.06105 -0.15662 -0.00992 0.12449 0.00212 0.00251 -0.20139;. -0.03446 -0.10241 0.20516 -0.02083 -0.01183 -0.00595 0.14492 -0.07655;. 0.00053 0.24329 -0.06455 -0.04113 -0.06522 0.07034 -0.00426 -0.00891;. 0.23781 -0.06017 -0.07787 0.09298 0.04690 0.00008 0.16593 0.0006

38、9;c=110.2965 131.2289 -108.873 77.4817 132.97447 120.6633;pout=ok*pin'+c'附錄6function f=zhenfun2(x)x0=120 73 180 80 125 125 81.1 90;x1=150 81 218.2 99.5 135 150 102.1 117;v=2.2 1 3.2 1.3 1.8 2 1.4 1.8;duan=70, 0 ,50 ,0, 0, 30 ,0, 0 ,0, 40;. 30, 0, 20 ,8 ,15 ,6 ,2, 0, 0, 8;. 110 ,0 ,40 ,0, 30

39、,0 ,20, 40, 0, 40;. 55 ,5 ,10, 10, 10, 10, 15 ,0 ,0 ,1;. 75, 5 ,15 ,0 ,15 ,15, 0,10 ,10 ,10;. 95, 0, 10, 20, 0, 15 ,10 ,20 ,0 ,10;. 50 ,15 ,5 ,15 ,10, 10 ,5 ,10 ,3 ,2;. 70, 0, 20, 0 ,20, 0 ,20, 10, 15, 5; price=-505 0 124 168 210 252 312 330 363 489;. -560 0 182 203 245 300 320 360 410 495;. -610 0

40、152 189 233 258 308 356 415 500;. -500 150 170 200 255 302 325 380 435 800;. -590 0 116 146 188 215 250 310 396 510;. -607 0 159 173 205 252 305 380 405 520;. -500 120 180 251 260 306 315 335 348 548;. -800 153 183 233 253 283 303 318 400 800; n0=6 6 5 6 6 6 5 10;b=0 0 0 0 0 0 0 0 ;p=0 0 0 0 0 0 0 0

41、 ;for i=1:8 if x(i)-x0(i)<0 xi=x(i)-x0(i); t=0; if xi<0 for n=n0(i):-1:1 b(i)=duan(i,n)+b(i); if b(i)+xi>=0 b(i)=-duan(i,n)+b(i); break; end end for ni=n0(i):-1:n if ni=n if duan(i,ni)=0 p(i)=-duan(i,ni)2/v(i)/2*(303-price(i,ni)+p(i); t=t+duan(i,ni)/v(i); end else if xi<=x0(i)-x1(i) p(i)

42、=-(xi+b(i)2/v(i)/2*(303-price(i,ni)+(0.25-t-(x1(i)-x0(i)/v(i)*(xi+b(i)*(303-price(i,ni)+p(i); else p(i)=-(xi+b(i)2/v(i)/2*(303-price(i,ni)+(0.25-t-(x1(i)+xi-x0(i)/v(i)*(xi+b(i)*(303-price(i,ni)+p(i); end end end else for n=n0(i):10 b(i)=duan(i,n)+b(i); if b(i)-xi>=0 b(i)=-duan(i,n)+b(i); break; end end for ni=n0(i):n if ni=n if duan(i,ni)=0 p(i)=duan(i,ni)2/v(i)/2*(price(i,ni)-303)+p(i); t=t+duan(i,ni)/v(i); end else p(i)=(xi-b(i)2/v(i)/2*(price(i,ni)-303)+(0.25-t)*(xi-b(i)*(303-price(i,ni)+p(i); end endend elseif x(i)-x1(i)<0 xi=x(i)-x1(i); t=0; if xi<0 for n=n0(i):-1:1 b(