電力市場的輸電阻塞管理—全國數(shù)學(xué)建模大賽

1、精選文檔電力市場的輸電堵塞管理摘要針對問題一，即對各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的爭辯，我們通過線性回歸分析的思想，建立了的數(shù)學(xué)模型。然后利用MATLAB得到具體參數(shù)，從而確定各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式。通過對置信區(qū)間的分析，檢查回歸變量對因變量的影響的顯著程度。針對問題二，在輸電堵塞發(fā)生時，堵塞費(fèi)用計(jì)算應(yīng)包括兩部分：一、通過競價取得發(fā)電權(quán)的發(fā)電容量（稱序內(nèi)容量）不能出力所引起的發(fā)電的成本提高和利潤削減的部分；二、因一些在競價中未取得發(fā)電權(quán)的發(fā)電容量（稱序外容量）要在低于對應(yīng)報(bào)價的清算價上出力所引起的額外出力的費(fèi)用。為了公正起見，第一部分的堵塞費(fèi)用不應(yīng)由網(wǎng)方全部擔(dān)當(dāng)，具

2、體折扣由買賣雙方協(xié)定，此處我們假定為50%，而其次部分的堵塞費(fèi)用由網(wǎng)方全部擔(dān)當(dāng)。針對問題三，依據(jù)題目給出的電力市場交易原則，下一個時段預(yù)報(bào)的負(fù)荷需求和各機(jī)組的段容量、段價和爬坡速率的數(shù)據(jù)，可以得到下一個時段各機(jī)組的出力安排預(yù)案。針對問題四，依據(jù)題目給出的潮流限值，有問題一得到的各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式可以檢測問題三的出力安排預(yù)案是否會引起輸電堵塞；在檢測到發(fā)生輸電堵塞時，依據(jù)各機(jī)組出力對各線路潮流值的影響比例大小，調(diào)整各機(jī)組出力安排方案，并給出與該方案相應(yīng)的堵塞費(fèi)用。針對問題五，依據(jù)對問題三和問題四的分析方法，可以得到相應(yīng)結(jié)論。關(guān)鍵詞：回歸分析 出力安排預(yù)案 堵塞費(fèi)用一、

3、問題重述設(shè)某電網(wǎng)有若干臺發(fā)電機(jī)組和若干條主要線路，每條線路上的有功潮流（輸電功率和方向）取決于電網(wǎng)結(jié)構(gòu)和各發(fā)電機(jī)組的出力。電網(wǎng)每條線路上的有功潮流的確定值有一平安限值，限值還具有肯定的相對平安裕度（即在應(yīng)急狀況下潮流確定值可以超過限值的百分比的上限）。假如各機(jī)組出力安排方案使某條線路上的有功潮流的確定值超出限值，稱為輸電堵塞。當(dāng)發(fā)生輸電堵塞時，需要爭辯如何制訂既平安又經(jīng)濟(jì)的調(diào)度方案。l 電力市場交易規(guī)章：1. 以15分鐘為一個時段組織交易，每臺機(jī)組在當(dāng)前時段開頭時刻前給出下一個時段的報(bào)價。各機(jī)組將可用出力由低到高分成至多10段報(bào)價，每個段的長度稱為段容量，每個段容量報(bào)一個價（稱為段價），段價按

4、段序數(shù)單調(diào)不減。在最低技術(shù)出力以下的報(bào)價一般為負(fù)值，表示情愿付費(fèi)維持發(fā)電以避開停機(jī)帶來更大的損失。2. 在當(dāng)前時段內(nèi)，市場交易-調(diào)度中心依據(jù)下一個時段的負(fù)荷預(yù)報(bào),每臺機(jī)組的報(bào)價、當(dāng)前出力和出力轉(zhuǎn)變速率，按段價從低到高選取各機(jī)組的段容量或其部分（見下面注釋），直到它們之和等于預(yù)報(bào)的負(fù)荷，這時每個機(jī)組被選入的段容量或其部分之和形成該時段該機(jī)組的出力安排預(yù)案（初始交易結(jié)果）。最終一個被選入的段價（最高段價）稱為該時段的清算價，該時段全部機(jī)組的全部出力均按清算價結(jié)算。注釋：（a） 每個時段的負(fù)荷預(yù)報(bào)和機(jī)組出力安排方案的參照時刻均為該時段結(jié)束時刻。（b） 機(jī)組當(dāng)前出力是對機(jī)組在當(dāng)前時段結(jié)束時刻實(shí)際出力的

5、猜測值。（c） 假設(shè)每臺機(jī)組單位時間內(nèi)能增加或削減的出力相同，該出力值稱為該機(jī)組的爬坡速率。由于機(jī)組爬坡速率的約束，可能導(dǎo)致選取它的某個段容量的部分。（d） 為了使得各機(jī)組方案出力之和等于預(yù)報(bào)的負(fù)荷需求，清算價對應(yīng)的段容量可能只選取部分。市場交易-調(diào)度中心在當(dāng)前時段內(nèi)要完成的具體操作過程如下：1、 監(jiān)控當(dāng)前時段各機(jī)組出力安排方案的執(zhí)行，調(diào)度AGC幫助服務(wù)，在此基礎(chǔ)上給出各機(jī)組的當(dāng)前出力值。2、 作出下一個時段的負(fù)荷需求預(yù)報(bào)。3、 依據(jù)電力市場交易規(guī)章得到下一個時段各機(jī)組出力安排預(yù)案。4、 計(jì)算當(dāng)執(zhí)行各機(jī)組出力安排預(yù)案時電網(wǎng)各主要線路上的有功潮流，推斷是否會消滅輸電堵塞。假如不消滅，接受各機(jī)組出

6、力安排預(yù)案；否則，依據(jù)如下原則實(shí)施堵塞管理：l 輸電堵塞管理原則：（1） 調(diào)整各機(jī)組出力安排方案使得輸電堵塞消退。（2） 假如（1）做不到，還可以使用線路的平安裕度輸電，以避開拉閘限電（強(qiáng)制削減負(fù)荷需求），但要使每條線路上潮流的確定值超過限值的百分比盡量小。（3） 假如無論怎樣安排機(jī)組出力都無法使每條線路上的潮流確定值超過限值的百分比小于相對平安裕度，則必需在用電側(cè)拉閘限電。（4） 當(dāng)轉(zhuǎn)變依據(jù)電力市場交易規(guī)章得到的各機(jī)組出力安排預(yù)案時，一些通過競價取得發(fā)電權(quán)的發(fā)電容量（稱序內(nèi)容量）不能出力；而一些在競價中未取得發(fā)電權(quán)的發(fā)電容量（稱序外容量）要在低于對應(yīng)報(bào)價的清算價上出力。因此，發(fā)電商和網(wǎng)方將產(chǎn)

7、生經(jīng)濟(jì)利益沖突。網(wǎng)方應(yīng)當(dāng)為因輸電堵塞而不能執(zhí)行初始交易結(jié)果付出代價，網(wǎng)方在結(jié)算時應(yīng)當(dāng)適當(dāng)?shù)亟o發(fā)電商以經(jīng)濟(jì)補(bǔ)償，由此引起的費(fèi)用稱之為堵塞費(fèi)用。網(wǎng)方在電網(wǎng)平安運(yùn)行的保證下應(yīng)當(dāng)同時考慮盡量削減堵塞費(fèi)用。你需要做的工作如下：1. 某電網(wǎng)有8臺發(fā)電機(jī)組，6條主要線路，表1和表2中的方案0給出了各機(jī)組的當(dāng)前出力和各線路上對應(yīng)的有功潮流值，方案132給出了圍繞方案0的一些試驗(yàn)數(shù)據(jù)，試用這些數(shù)據(jù)確定各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式。2. 設(shè)計(jì)一種簡明、合理的堵塞費(fèi)用計(jì)算規(guī)章，除考慮上述電力市場規(guī)章外，還需留意：在輸電堵塞發(fā)生時公正地對待序內(nèi)容量不能出力的部分和報(bào)價高于清算價的序外容量出力的部分。

8、 3. 假設(shè)下一個時段預(yù)報(bào)的負(fù)荷需求是982.4MW，表3、表4和表5分別給出了各機(jī)組的段容量、段價和爬坡速率的數(shù)據(jù)，試依據(jù)電力市場規(guī)章給出下一個時段各機(jī)組的出力安排預(yù)案。4. 依據(jù)表6給出的潮流限值，檢查得到的出力安排預(yù)案是否會引起輸電堵塞，并在發(fā)生輸電堵塞時，依據(jù)平安且經(jīng)濟(jì)的原則，調(diào)整各機(jī)組出力安排方案，并給出與該方案相應(yīng)的堵塞費(fèi)用。5. 假設(shè)下一個時段預(yù)報(bào)的負(fù)荷需求是1052.8MW，重復(fù)34的工作。二、問題分析我們進(jìn)行了如下的分析：Step1: 首先我們明確問題中所需做的幾個任務(wù)：1. 各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式；2. 一種簡明、合理的堵塞費(fèi)用計(jì)算規(guī)章；3. 假設(shè)下

9、一個時段預(yù)報(bào)的負(fù)荷需求是982.4MW，給出下一個時段各機(jī)組的出力安排預(yù)案；4. 檢測并調(diào)整安排預(yù)案，并計(jì)算相應(yīng)的堵塞費(fèi)用。Step2: 即對各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的爭辯，我們通過線性回歸分析的思想，建立了公式的數(shù)學(xué)模型。然后利用MATLAB得到具體參數(shù)，從而確定各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式。通過對置信區(qū)間的分析，檢查回歸變量對因變量的影響的顯著程度。Step3: 在輸電堵塞發(fā)生時，堵塞費(fèi)用計(jì)算應(yīng)包括兩部分：一、通過競價取得發(fā)電權(quán)的發(fā)電容量（稱序內(nèi)容量）不能出力所引起的發(fā)電的成本提高和利潤削減的部分；二、因一些在競價中未取得發(fā)電權(quán)的發(fā)電容量（稱序外容量）要在低于對

10、應(yīng)報(bào)價的清算價上出力所引起的額外出力的費(fèi)用。為了公正起見，第一部分的堵塞費(fèi)用不應(yīng)由網(wǎng)方全部擔(dān)當(dāng)，具體折扣由買賣雙方協(xié)定，此處我們假定為50%，而其次部分的堵塞費(fèi)用由網(wǎng)方全部擔(dān)當(dāng)。Step4: 依據(jù)題目給出的電力市場交易原則，下一個時段預(yù)報(bào)的負(fù)荷需求和各機(jī)組的段容量、段價和爬坡速率的數(shù)據(jù)，可以得到下一個時段各機(jī)組的出力安排預(yù)案。Step5: 依據(jù)題目給出的潮流限值，有問題一得到的各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式可以檢測問題三的出力安排預(yù)案是否會引起輸電堵塞；在檢測到發(fā)生輸電堵塞時，依據(jù)各機(jī)組出力對各線路潮流值的影響比例大小，調(diào)整各機(jī)組出力安排方案，并給出與該方案相應(yīng)的堵塞費(fèi)用。St

11、ep6: 依據(jù)對問題三和問題四的分析方法，可以得到相應(yīng)結(jié)論。三、符號設(shè)定 第i條線路的有功潮流 第i個機(jī)組的出力 各部分堵塞費(fèi)用的總和發(fā)電廠每度電的變動成本發(fā)電廠發(fā)電的固定成本給定的各機(jī)組的段容量出力預(yù)案中具體的各機(jī)組中的段容量給定的各機(jī)組的清算價出力預(yù)案中具體的各機(jī)組中的清算價四、模型建立與求解(一)、問題一的求解5.1 記第i條線路的有功潮流為，第i個機(jī)組的出力為，依據(jù)題目給出的32種方案的各機(jī)組的出力值及相對應(yīng)的各線路的潮流值，基于上面的分析，我們首先利用來建立的近似表達(dá)式。模型預(yù)備為了大致分析與的關(guān)系，首先利用表1和表2的數(shù)據(jù)可以做出對的散點(diǎn)圖，依此類推可分別得到對的散點(diǎn)圖。圖1 對的

12、散點(diǎn)圖模型建立從圖1可以發(fā)覺，隨著的增加，的值有比較明顯的線性增長趨勢，圖中的直線是用線性模型 （5.1.1）擬合的，（其中是誤差）同樣經(jīng)過分析，發(fā)覺隨著的增加，的值有比較明顯的線性增長趨勢，圖中的直線是用線性模型 （5.1.2）綜合上面的分析，結(jié)合上述模型簡歷如下的回歸模型（5.1.3）（5.1.3）右端的稱為回歸變量（自變量），是給定機(jī)組出力，第一條線路的潮流值的平均值，其中的參數(shù)稱為回歸系數(shù)，由表1、表2的數(shù)據(jù)估量，影響的其他因素作用都包含在隨機(jī)誤差中。假如模型選擇得合適，應(yīng)大致聽從均值為零的正態(tài)分布。依此類推，我們可以用同樣的方法得出分別對的基本模型 （5.1.4）模型求解直接利用MA

13、TLAB統(tǒng)計(jì)工具箱中的命令regress求解，使用格式為： b,bint,r,rint,stats=regress(y1(:),x0(:) x1(:) x2(:) x3(:) x4(:) x5(:) x6(:) x7(:) x8(:)其中輸入y為模型（5.1.3）中y的數(shù)據(jù)，x為對應(yīng)于回歸系數(shù)的數(shù)據(jù)矩陣，置信水平。得到模型（5.1.3）的回歸系數(shù)估量值及其置信區(qū)間（置信水平）的結(jié)果見下表（5.1）。參數(shù)參數(shù)值估量參數(shù)置信區(qū)間0110.4775 109.5421 111.412910.0826 0.0808 0.084420.0478 0.0437 0.051830.0528 0.0514 0.

14、054240.1199 0.1166 0.12315-0.0257 -0.0277 -0.02376 0.12160.1190 0.12437 0.12200.1189 0.12518 -0.0015-0.0037 0.0007表（5.1） 模型（5.1.3）的計(jì)算結(jié)果同樣，以此類推我們可以得出其他各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式： （5.1.5）yi: =0 1 2 3 4 5 6 7 8其中，y1: =110.4775 0.0826 0.0478 0.0528 0.1199 -0.0257 0.1216 0.1220 -0.0015y2: =131.3521 0.0547 0

15、.1275 -0.0001 0.0332 0.0867 -0.1127 -0.0186 0.0985y3: =-108.99281 -0.0694 0.0620 -0.1565 -0.0099 0.01247 0.0024 -0.0028 -0.2012y4: =77.6116 -0.0346 -0.1028 0.2050 -0.0209 -0.0120 0.0057 0.1452 0.0763y5: =133.1334 -0.0003 0.2428 -0.0647 -0.0412 -0.0655 0.0700 -0.0039 -0.0092y6: =120.8481 0.2376 -0.06

16、07 -0.0781 -0.0929 0.0466 -0.0003 0.1664 0.0004結(jié)果分析表（5.1）的回歸系數(shù)給出了模型（5.1.3）中的估量值，檢查他們的置信區(qū)間發(fā)覺，只有8置信區(qū)間包含零點(diǎn)，表示回歸變量8（對因變量的影響）不是太顯著，但我們?nèi)詫?保留在模型中。(二)、問題二的求解5.2 在輸電堵塞發(fā)生時，堵塞費(fèi)用（記為W）計(jì)算應(yīng)包括兩部分：一、通過競價取得發(fā)電權(quán)的發(fā)電容量（稱序內(nèi)容量）不能出力所引起的發(fā)電的成本提高（記為W1）和利潤削減的部分（記為W2）；二、因一些在競價中未取得發(fā)電權(quán)的發(fā)電容量（稱序外容量）要在低于對應(yīng)報(bào)價的清算價上出力所引起的出力的額外費(fèi)用（記為W3）。記

17、給定的各機(jī)組的段容量為，記出力預(yù)案中具體的各機(jī)組中的段容量為，記給定的各機(jī)組的清算價為，記出力預(yù)案中具體的各機(jī)組中的清算價為。模型預(yù)備為了公正起見，第一部分的堵塞費(fèi)用不應(yīng)由網(wǎng)方全部擔(dān)當(dāng)，具體折扣由買賣雙方協(xié)定，此處我們假定為50%，而其次部分的堵塞費(fèi)用由網(wǎng)方全部擔(dān)當(dāng)。即 （5.2）模型建立5.2.1 計(jì)算發(fā)電的成本提高W1發(fā)電廠的發(fā)電成本分為兩部分：固定成本：包括廠房設(shè)備、保險、材料存貯、設(shè)備運(yùn)行等，不隨發(fā)電 度數(shù)的變化而變化；平均變動成本：包括原材料、燃料、工資等，隨發(fā)電度數(shù)的變化而變化。所以由給定的各機(jī)組的段容量為，出力預(yù)案中具體的各機(jī)組中的段容量為，發(fā)電的成本提高W1的值確定為 （5.2

18、.1）5.2.2 計(jì)算發(fā)電的利潤削減W2由給定的各機(jī)組的段容量為，出力預(yù)案中實(shí)際的機(jī)組中的段容量為，故少發(fā)電為則發(fā)電的利潤削減W2的值確定為 （5.2.2）5.2.3 計(jì)算出力的額外費(fèi)用W3由給定的各機(jī)組的清算價為，出力預(yù)案中具體的各機(jī)組中的清算價為，故多出的費(fèi)用為則發(fā)電的出力的額外費(fèi)用W3的值確定為 （5.2.3）5.2.4 計(jì)算堵塞費(fèi)用W由（5.2）可得堵塞費(fèi)用的值為 （5.2.4）(三)、問題三的求解5.3 已知下一階段的負(fù)荷量，我們依據(jù)電力市場交易規(guī)章給出相應(yīng)的出力安排預(yù)案。按段價從低到高選取各機(jī)組的段容量或其部分，直到它們之和等于預(yù)報(bào)的負(fù)荷，這時每個機(jī)組被選入的段容量或其部分之和形成

19、該時段該機(jī)組的出力安排預(yù)案。以下表格中的紅色部分就是依據(jù)電力市場交易規(guī)章所取的各機(jī)組的段容量對應(yīng)的段價。機(jī)組段123456789101-50501241682102523123303634892-56001822032453003203604104953-61001521892332583083564155004-5001501702002553023253804358005-59001161461882152503103965106-60701591732052523053804055207-5001201802512603063153353485488-800153183233253283

20、303318400800因考慮到各機(jī)組的最終出力總和應(yīng)當(dāng)在最低出力和到最高出力這個范圍內(nèi)，以及最終的總負(fù)荷量的大小，所以可能有些段容量只能取其中的一部分。以下就是依據(jù)上面的兩個約束條件作出調(diào)整后的結(jié)果。由于第四機(jī)組的最高出力只能是99.5，所以第四機(jī)組的第六段只能出力9.5，這樣才能保證不超過最高出力。當(dāng)清算價為302（第四機(jī)組第六段）時，段容量之和為978.5，離目標(biāo)負(fù)荷只差3.9，所以在下一段段價為303（第八機(jī)組第七段）只要取3.9就能達(dá)到目標(biāo)值了。機(jī)組12345678總出力當(dāng)前出力120731808012512581.190最高出力1538822899.5152155102.1117最

21、低出力875813260.5989560.163第一出力70301105575955070555其次出力/段5/215/320第三出力/段15/215第四出力/段50/340/35/295第五出力/段20/310/310/415/510/320/45/320/3110第七出力/段30/68/415/530/510/415/615/615/420/5158出力和一150731809012514085110953第九出力/段6/69.5/63.9/710/529.4出力和1507918099.5125143995110982.4所以下一階段出力安排預(yù)案為：機(jī)組12345678出力150791809

22、9.5125143995110清算價為：303(4) 、問題四的求解5.4.1推斷是否有輸電堵塞現(xiàn)象 利用問題一中求出的各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式，將各機(jī)組出力的具體數(shù)據(jù)代入計(jì)算，得出各線路上的潮流值。然后與題目所給的各線路潮流限值相比較，觀看是否有輸電堵塞現(xiàn)象。計(jì)算線路一上的有功潮流：y1=110.4775+0.0826*150+0.0478*79+0.0528*180+0.1199*99.5-0.0257*125+0.1216*143.9+0.1220*95-0.0015*110y1 =173.7885 165則線路一上有堵塞現(xiàn)象。依次我們可以得出y2=131.3521

23、-0.0547*150+0.1275*79-0.0001*180+0.0332*99.5+0.0867*125-0.1127*143.9-0.0186*95+0.0985*110y2=140.1930150所以線路二上無堵塞現(xiàn)象y3=-108.99281-0.0694*150+0.0620*79-0.1565*180-0.0099*99.5+0.01247*125+0.0024*143.9-0.0028*95-0.2012*110y3 =-164.1538 160所以線路三上無堵塞現(xiàn)象y4=77.6116-0.0346*150-0.1028*79+0.2050*180-0.0209*99.5-

24、0.0120*125+0.0057*143.9+0.1452*95+0.0763*110y4 =120.6281132所以線路五上有堵塞現(xiàn)象y6=120.8481+0.2376*150-0.0607*79-0.0781*180-0.0929*99.5+0.0466*125-0.0003*143.9+0.1664*9+0.0004*110y6 =135.7147162所以線路六上無堵塞現(xiàn)象綜上所述，線路一和線路五會消滅輸電堵塞現(xiàn)象5.4.2.1調(diào)整出力安排方案 在發(fā)生輸電堵塞時，依據(jù)平安且經(jīng)濟(jì)的原則，調(diào)整各機(jī)組出力安排方案。首先我們依據(jù)各線路上有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式得出各機(jī)組出力

25、對線路潮流值的影響大小，如下表（表中1代表影響最大，8代表影響最?。C(jī)組12345678線路146537218線路251864273線路345263871線路453167824線路581453276線路615436827由于我們調(diào)整出力安排方案的原則是做盡量少的調(diào)整而能使得輸電堵塞消退，所以我們要進(jìn)行調(diào)整的首要對象是對線路一和五影響較大而對其他線路影響較小的機(jī)組。其中線路一超出范圍較大，所以優(yōu)先考慮對線路一的調(diào)整。由表可看出機(jī)組6、7的出力大小對線路一的潮流值的影響最大。所以可以使機(jī)組6、7的出力減小來達(dá)到線路一的潮流值變小。而要保持最終的出力總和不變，則應(yīng)有相應(yīng)機(jī)組的出力增加。機(jī)組5、8的變

26、化對線路一潮流值的影響較小，并且是負(fù)相關(guān)的。所以我們可以相應(yīng)的增加機(jī)組5,、8的出力來保持出力總和不變，同時還可以減小線路一的潮流值。在調(diào)整的過程中，應(yīng)時刻留意各機(jī)組的出力總和應(yīng)當(dāng)在最小出力和最大出力范圍內(nèi)，并且相對考慮一下段價的大小，盡量避開段價過高，使得堵塞費(fèi)用變大。在調(diào)整的過程中，我們發(fā)覺調(diào)整后使得線路一的潮流值小于潮流限值，相應(yīng)的線路五的潮流值也會小于潮流限值，反而線路二很簡潔就超出潮流限值的范圍。所以接下來還需要做微小的調(diào)整，使得全部線路的潮流值都在潮流限值的范圍內(nèi)，即無輸電堵塞。所以調(diào)整后的出力方案為：機(jī)組12345678調(diào)整前1507918099.512514095113.9調(diào)整

27、后15386.522890.81529560.1117最高出力1538822899.5152155102.1117最低出力875813260.5989560.1635.4.2.2 檢驗(yàn)是否存在輸電堵塞現(xiàn)象y1=110.4775+0.0826*153+0.0478*86.5+0.0528*228+0.1199*90.8-0.0257*152+0.1216*95+0.1220*60.1-0.0015*117y1 =164.9776165 y2=131.3521-0.0547*153+0.1275*86.5-0.0001*228+0.0332*90.8+0.0867*152-0.1127*95-0.

28、0186*60.1+0.0985*117y2 =149.8820150y3=-108.99281-0.0694*153+0.0620*86.5-0.1565*228-0.0099*90.8+0.01247*152+0.0024*95-0.0028*60.1-0.2012*117y3 = -172.4142160y4=77.6116-0.0346*153-0.1028*86.5+0.2050*228-0.0209*90.8-0.0120*152+0.0057*95+0.1452*60.1+0.0763*117y4 = 124.6390155y5=133.1334-0.0003*153+0.242

29、8*86.5-0.0647*228-0.0412*90.8-0.0655*152+0.0700*95-0.0039*60.1-0.0092*117y5 =130.9803132y6=120.8481+0.2376*153-0.0607*86.5-0.0781*228-0.0929*90.8+0.0466*152-0.0003*95+0.1664*60.1+0.0004*117y6 =142.8104165則線路一上堵塞現(xiàn)象依次我們可以得出y2=131.3521-0.0547*150+0.1275*81-0.0001*218.2+0.0332*99.5+0.0867*135-0.1127*150

30、-0.0186*102.1+0.0985*117y2 =141.1811 150所以線路二上無堵塞現(xiàn)象y3=-108.99281-0.0694*150+0.0620*81-0.1565*218.2-0.0099*99.5+0.01247*135+0.0024*150-0.0028*102.1-0.2012*117y3 =-171.2970 160所以線路三上無堵塞現(xiàn)象y4=77.6116-0.0346*150-0.1028*81+0.2050*218.2-0.0209*99.5-0.0120*135+0.0057*150+0.1452*102.1+0.0763*117y4 =129.7333

31、132所以線路五上有堵塞現(xiàn)象y6=120.8481+0.2376*150-0.0607*81-0.0781*218.7-0.0929*99.5+0.0466*135-0.0003*150+0.1664*102.1+0.0004*117y6 =148.5687 165 y2=131.3521-0.0547*153+0.1275*88-0.0001*228+0.0332*99.5+0.0867*152-0.1127*113.2-0.0186*102.1+0.0985*117y2 =147.5298 150 y3=-108.99281-0.0694*153+0.0620*88-0.1565*228-

32、0.0099*99.5+0.01247*152+0.0024*113.2-0.0028*102.1-0.2012*117y3 =-172.4812 160 y4=77.6116-0.0346*153-0.1028*88+0.2050*228-0.0209*99.5-0.0120*152+0.0057*113.2+0.1452*102.13+0.0763*117y4 =130.5095 132y6=120.8481+0.2376*153-0.0607*88-0.0781*228-0.0929*99.5+0.0466*152-0.0003*113.2+0.1664*102.1+0.0004*117

33、y6 =148.8944 x0=1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1; x1=120 133.02 129.63 158.77 145.32 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120;x2=73 73 73 73 73 78.596 75.45 90.487 83.848 73 73 73 73 73 73 73 73 7

34、3 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73;x3=180 180 180 180 180 180 180 180 180 231.39 198.48 212.64 190.55 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180;x4=80 80 80 80 80 80 80 80 80 80 80 80 80 75.857 65.958 87.258 97.824 80 80 80 80 80 80 80 80 80 80 80 80

35、80 80 80 80;x5=125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 150.71 141.58 132.37 156.93 125 125 125 125 125 125 125 125 125 125 125 125;x6=125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 138.88 131.21 141.71 149.29 125 125 125 125 125 125 12

36、5 125;x7=81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 60.582 70.962 64.854 75.529 81.1 81.1 81.1 81.1;x8=90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 104.84 111.22 98.092 120.44;

37、 y1=164.78 165.81 165.51 167.93 166.79 164.94 164.8 165.59 165.21 167.43 165.71 166.45 165.23 164.23 163.04 165.54 166.88 164.07 164.27 164.57 163.89 166.35 165.54 166.75 167.69 162.21 163.54 162.7 164.06 164.66 164.7 164.67 164.69; y2=140.87 140.13 140.25 138.71 139.45 141.5 141.13 143.03 142.28 14

38、0.82 140.82 140.82 140.85 140.73 140.34 141.1 141.4 143.03 142.29 141.44 143.61 139.29 140.14 138.95 138.07 141.21 141 141.14 140.94 142.27 142.94 141.56 143.84; y3=-144.25 -145.14 -144.92 -146.91 -145.92 -143.84 -144.07 -143.16 -143.49 -152.26 -147.08 -149.33 -145.82 -144.18 -144.03 -144.32 -144.34 -140.97 -142.15