數(shù)學(xué)建模論文電力市場(chǎng)輸電阻塞管理數(shù)學(xué)模型

1、電力市場(chǎng)輸電阻塞管理數(shù)學(xué)模型隊(duì)員： 摘要：本文通過(guò)求出各線(xiàn)路有功潮流與各機(jī)組的出力之間的關(guān)系，并設(shè)計(jì)合理的阻塞費(fèi)用計(jì)算規(guī)則，建立了電力市場(chǎng)的輸電阻塞管理模型。我們運(yùn)用matlab數(shù)學(xué)軟件、線(xiàn)性回歸模型對(duì)表中數(shù)據(jù)進(jìn)行擬合，并對(duì)數(shù)據(jù)進(jìn)行了檢驗(yàn)，初步得到近似的線(xiàn)性表達(dá)式。根據(jù)負(fù)荷需求，有功潮流限值，安全裕度，各機(jī)組的爬坡速率，阻塞費(fèi)用最少的條件，建立了雙目標(biāo)規(guī)劃模型得到最優(yōu)的處理方案。關(guān)鍵詞：電力市場(chǎng) 有功潮流 輸電阻塞管理 多元線(xiàn)性回歸一、 問(wèn)題重述 我國(guó)電力市場(chǎng)初期是發(fā)電側(cè)電力市場(chǎng)，采取交易與調(diào)度一體化的模式。電網(wǎng)公司在組織交易、調(diào)度和配送時(shí)，必須遵循電網(wǎng)“安全第一”的原則，同時(shí)要制訂一個(gè)電力市

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

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

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

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

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

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

8、一個(gè)時(shí)段預(yù)報(bào)的負(fù)荷需求是1052.8mw，重復(fù)34的工作。二、問(wèn)題分析問(wèn)題要求在滿(mǎn)足”安全第一“和盡量避免輸電阻塞的情況下，尋求經(jīng)濟(jì)的各機(jī)組的出力分配方案，使得購(gòu)電費(fèi)用最小。我們的目標(biāo)是使得購(gòu)電費(fèi)用最小，約束是各線(xiàn)路的有功潮流限值與各機(jī)組的出力分配方案可能引起的輸電阻塞。以及爬坡速度對(duì)下一時(shí)段出力范圍的限制。故：1. 針對(duì)問(wèn)題1根據(jù)表1表2數(shù)據(jù)，我們運(yùn)用matlab軟件對(duì)數(shù)據(jù)進(jìn)行描點(diǎn)，發(fā)現(xiàn)他們屬于線(xiàn)性關(guān)系，并進(jìn)行一次線(xiàn)性擬合，處理后可得到各線(xiàn)路有功潮流值與各機(jī)組出力的近似表達(dá)式。2. 設(shè)計(jì)一種簡(jiǎn)明合理的阻塞費(fèi)用計(jì)算規(guī)則。首先，由于出現(xiàn)了輸電阻塞，通過(guò)改變市場(chǎng)交易規(guī)則而調(diào)整各機(jī)組的出力分配方案，

9、這樣便出現(xiàn)了部分序內(nèi)容量不能出力的補(bǔ)償和序外容量低于清算價(jià)計(jì)費(fèi)的補(bǔ)償。然后，我們要盡量減少阻塞費(fèi)用。這里我們通過(guò)比較，以序內(nèi)不能出力部分的報(bào)價(jià)和序外出力的報(bào)價(jià)為參考，設(shè)計(jì)了一種簡(jiǎn)明合理的阻塞費(fèi)用計(jì)算方法。3在滿(mǎn)足負(fù)荷需求的條件下，考慮到各機(jī)組爬坡速率的限制以及盡量使得費(fèi)用最少，依次取費(fèi)用最少的段，直到滿(mǎn)足負(fù)荷需求，得到使得各機(jī)組出力之和滿(mǎn)足負(fù)荷的出力分配方案。4考慮到各路線(xiàn)的有功潮流限值，在保證安全的前提下，盡量使得阻塞費(fèi)用最小，列出目標(biāo)函數(shù)和約束條件的線(xiàn)性規(guī)劃，用matlab進(jìn)行求解。得出調(diào)整后的出力方案。5.重復(fù)3、4三基本假設(shè)1.各線(xiàn)路的有功潮流限值是不會(huì)改變的。即電網(wǎng)的結(jié)構(gòu)是不會(huì)改變的

10、。2.每一時(shí)段，發(fā)點(diǎn)上的各機(jī)組出力是穩(wěn)定的。3.不考慮電的輸送過(guò)程中的損失。即輸電過(guò)程是理想狀態(tài)。4.市場(chǎng)交易-調(diào)度中心有權(quán)要求控制區(qū)內(nèi)的各發(fā)電機(jī)組無(wú)條件地執(zhí)行agc控制信號(hào)。四符號(hào)說(shuō)明-各機(jī)組的出力-部分機(jī)組調(diào)整后序內(nèi)容量出力部分-部分機(jī)組調(diào)整后序外容量出力部分-各線(xiàn)路有功潮流值-各線(xiàn)路調(diào)整后的潮流值-各線(xiàn)路有功潮流最大值-各線(xiàn)路的安全裕度-各線(xiàn)路潮流值超過(guò)限值的百分比-清算價(jià) -各段報(bào)價(jià)-各機(jī)組段價(jià)與清算價(jià)之差w-總的阻塞費(fèi)用-各機(jī)組的爬坡速率q-負(fù)荷需求z-支付的總費(fèi)用（i=1,26；s=1,28）五模型建立1. 我們?cè)O(shè)各線(xiàn)路潮流值為(i=1,26) 各機(jī)組出力為(s=1,28)則 =(

11、)+()+.+()+因?yàn)橛^察實(shí)驗(yàn)數(shù)據(jù)，可以看出每四個(gè)數(shù)據(jù)分別取了其中一個(gè)機(jī)組的出力不同值，而其它七個(gè)機(jī)組的出力值都保持不變，我們把實(shí)驗(yàn)數(shù)據(jù)分成8組，每組中只有一個(gè)機(jī)組的出力變動(dòng)。則我們?cè)O(shè)=()+ . =()+ (其中為常數(shù))分別找到與的關(guān)系用matlab分別對(duì) 的8組數(shù)據(jù)描點(diǎn)有 對(duì)與 的數(shù)據(jù)有圖 125130135140145150155160165.5166166.5167167.516874767880828486889092164.8164.9165165.1165.2165.3165.4165.5165.6165.7 1901952002052102152202252302351651

12、65.5166166.5167167.5 65707580859095100163163.5164164.5165165.5166166.5167 其他組圖略去， 發(fā)現(xiàn)均成線(xiàn)性關(guān)系。用matlab 計(jì)算出各系數(shù)（見(jiàn)附錄）。可得到個(gè)潮流值與各組出力的關(guān)系。2. 方法一：設(shè)當(dāng)前各機(jī)組的出力為(s=1,28)，清算價(jià)為。當(dāng)前該出力方案使得某條線(xiàn)路出現(xiàn)輸電阻塞，所以必須進(jìn)行調(diào)整，調(diào)整會(huì)使得有些機(jī)組出力增加（即出現(xiàn)序外容量），有些機(jī)組出力減少（序內(nèi)容量）。設(shè)出力減少的那部分機(jī)組在減少出力后的出力為，通過(guò)競(jìng)價(jià)取得發(fā)電權(quán)的不能出力的那部分容量的最高段價(jià)為（）。出力增加的那部分機(jī)組在出力增加后的出力為；在競(jìng)價(jià)

13、中未取得發(fā)電權(quán)卻出力的那部分容量的最高段價(jià)為()。則：序內(nèi)容量補(bǔ)償費(fèi)用為：（ ）)k序外容量補(bǔ)償費(fèi)用為：（ ）k當(dāng) =0、 =0時(shí)，即表示第s機(jī)組的出力沒(méi)有進(jìn)行調(diào)整。其中k表示的是一個(gè)系數(shù)，范圍是0k1方法二：我們以序內(nèi)和序外報(bào)價(jià)與清算價(jià)的差額為補(bǔ)償基準(zhǔn)，即考慮調(diào)整后發(fā)電商損失的利潤(rùn)問(wèn)題：對(duì)于序內(nèi)容量，由于方案的調(diào)整，使得一些機(jī)組的出力值減少，減少部分的獲利值消失。為解決這部分沖突，網(wǎng)方賠償該機(jī)組應(yīng)得的獲利值，有，補(bǔ)償費(fèi)用=（清算價(jià)調(diào)整前報(bào)價(jià)）調(diào)整量；對(duì)于序外容量，方案調(diào)整后一些機(jī)組由于出力增加，其報(bào)價(jià)也隨之增加，但報(bào)價(jià)在其清算價(jià)以上，故也存在利潤(rùn)損失，網(wǎng)方應(yīng)對(duì)這部分出力所帶來(lái)的損失給予補(bǔ)償。

14、有，補(bǔ)償費(fèi)用=（調(diào)整后報(bào)價(jià)清算價(jià)）調(diào)整量。由此得到：序內(nèi)容量的補(bǔ)償費(fèi)用： （ ）)（） 序外容量的補(bǔ)償費(fèi)用： （ ）（）當(dāng)（ ）)（）=0，（ ）（）=0時(shí)，表示該機(jī)組沒(méi)有進(jìn)行調(diào)整，或者，調(diào)整過(guò)后不存在利潤(rùn)損失。因而，網(wǎng)方對(duì)此無(wú)需補(bǔ)償。故總的阻塞費(fèi)用為：w=（計(jì)算時(shí)我們選取的第二種方法）3. 由表5 各機(jī)組的爬坡速率 （單位：mw/分鐘）機(jī)組12345678速率2.213.21.31.821.41.8設(shè)當(dāng)前時(shí)段的出力為爬坡速率為 則由于爬坡速率的限制，可以算出下一時(shí)段各機(jī)組出力范圍：即下一時(shí)段各機(jī)組出力范圍為機(jī)組12345678當(dāng)前出力120731808012512581.190下一時(shí)段出力范

15、圍87153588813222860.599.5981529515560.1102.163117同時(shí)，要滿(mǎn)足負(fù)荷需求，即q=再次，盡量滿(mǎn)足網(wǎng)方支付的費(fèi)用最少。 則有目標(biāo)函數(shù) minz=qs.t.q=982.4mw=因?yàn)閦=q ，所以總費(fèi)用取決于出清價(jià)，又價(jià)格單調(diào)不減，因此我們依次從段價(jià)由低到高選取，并滿(mǎn)足爬坡速率的限制，各得到費(fèi)用最小的各機(jī)組的出力：機(jī)組12345678出力1507918099.512514095113.9此時(shí)=303 4. 由第1問(wèn)所解得的潮流值與各機(jī)組出力的近似表達(dá)式可算出各線(xiàn)路的潮流值：線(xiàn)路123456潮流值173.4186141.0771-150.848128.1975

16、136.9146185.2873根據(jù)表6各線(xiàn)路的潮流限值約束：線(xiàn)路123456限值165150160155132162安全裕度13%18%9%11%15%14%最大值186.45177174.4172.05151.8184.68可得到第3問(wèn)解答中的1，5.，6線(xiàn)路出現(xiàn)阻塞所以必須進(jìn)行調(diào)整，調(diào)整時(shí)有目標(biāo)1阻塞費(fèi)用最小，目標(biāo)2盡量不超過(guò)安全裕度或超過(guò)安全裕度盡量少，有目標(biāo)函數(shù)： w= （因?yàn)橥?fù)所以序內(nèi)和序外都為正，沒(méi)有分別列出，其中要分別列出可能的把它變?yōu)橐阎?s.t. 如果由以上幾個(gè)不等式組可以求出解即為調(diào)整的最優(yōu)解，如果不能得到解，那么就應(yīng)該進(jìn)行安全裕度輸電或者電測(cè)拉閘限電。則有目標(biāo)函數(shù)

17、：w= s.t. ) 首先將安全裕度前的系數(shù)k取為100%，如果有功潮流值不滿(mǎn)足，就用電測(cè)拉閘限電；如果滿(mǎn)足再將安全裕度前的系數(shù)減半，如果滿(mǎn)足再減半，如果不滿(mǎn)足再將安全裕度前的系數(shù)增加減半之后的一半，這樣一直進(jìn)行下去，最終得出最后選中的那個(gè)安全裕度將使得調(diào)整之后的有功潮流值超出限值的百分比最小。經(jīng)過(guò)調(diào)整可得新的各機(jī)組的出力方案機(jī)組12345678出力1508822880.41529670.1117此時(shí)各線(xiàn)路所對(duì)應(yīng)的潮流值（如下表所示）線(xiàn)路123456潮流值165.0605149.7603-155.0151131.2020131.9054171.1415相應(yīng)的阻塞費(fèi)用w=4602.15元5.若下

18、一時(shí)段的預(yù)報(bào)負(fù)荷需求是1052.8 mw，同樣由于爬坡速率以及購(gòu)電費(fèi)用最小的限制（具體數(shù)據(jù)已在第三題的解答中給出），在根據(jù)市場(chǎng)交易規(guī)則得出下一時(shí)段各機(jī)組的出力分配方案：機(jī)組12345678出力15081218.299.5135150102.1117同樣，根據(jù)第一題所得潮流與各機(jī)組出力關(guān)系可得到此時(shí)各線(xiàn)路的有功潮流值為線(xiàn)路123456潮流值177.3874141.2228-156.1046137.5563134.8907184.9970.當(dāng)需求的負(fù)荷為1052.8mw時(shí)，得出的各機(jī)組的出力分配方案經(jīng)計(jì)算使線(xiàn)路1,5,6上的有功潮流值超出其限值，所以必須進(jìn)行調(diào)整.調(diào)整時(shí)有目標(biāo)1阻塞費(fèi)用最小，目標(biāo)2

19、盡量不超過(guò)安全裕度或超過(guò)安全裕度盡量少，有目標(biāo)函數(shù)： w= （因?yàn)橥?fù)所以序內(nèi)和序外都為正，沒(méi)有分別列出，其中要分別列出可能的把它變?yōu)橐阎?s.t. 如果由以上幾個(gè)不等式組可以求出解即為調(diào)整的最優(yōu)解，如果不能得到解，那么就應(yīng)該進(jìn)行安全裕度輸電或者電測(cè)拉閘限電。則有目標(biāo)函數(shù)：w= s.t. ) 首先將安全裕度前的系數(shù)k取為100%，如果有功潮流值不滿(mǎn)足，就用電測(cè)拉閘限電；如果滿(mǎn)足再將安全裕度前的系數(shù)減半，如果滿(mǎn)足再減半，如果不滿(mǎn)足再將安全裕度前的系數(shù)增加減半之后的一半，這樣一直進(jìn)行下去，最終得出最后選中的那個(gè)安全裕度將使得調(diào)整之后的有功潮流值超出限值的百分比最小。 經(jīng)過(guò)調(diào)整可得新的各機(jī)組的出

20、力方案：機(jī)組12345678出力149.78822893.6150.5137.888.2117此時(shí)各線(xiàn)路所對(duì)應(yīng)的潮流值（如下表所示）線(xiàn)路123456潮流值173.9252144.9515-155.223136.0269134.3795180.1813相應(yīng)的阻塞費(fèi)用w=1578.25元六模型評(píng)價(jià)優(yōu)點(diǎn)：1.充分考慮到輸電的安全和經(jīng)濟(jì)的情況，得到較優(yōu)的分配方案。2利用matlab軟件，得到的誤差較小。3建立合理的線(xiàn)性規(guī)劃模型，得到較優(yōu)的方案。缺點(diǎn)：1. 在計(jì)算阻塞費(fèi)用時(shí)，不了解一些細(xì)節(jié)，因此賠償方案和阻塞費(fèi)用不一定最優(yōu)且最合理，費(fèi)用精確度不很高。2考慮的因素不完全，與實(shí)際有一定差距（老師，因?yàn)槲覀兪?/p>

21、第一次做，很多東西完全不了解，不得不承認(rèn)我們有參考別人的優(yōu)秀論文，希望原諒，謝謝?。└戒洠河?jì)算有功潮流關(guān)于各發(fā)電機(jī)組出力的近似表達(dá)式matlab源程序x11x=133.02 129.63 158.77 145.32y=165.81 165.51 167.93 166.79x = 133.0200 129.6300 158.7700 145.3200y = 165.8100 165.5100 167.9300 166.7900 polyfit(x,y,1)ans =0.0826 154.8091x12 x=78.596 75.45 90.487 83.848y=164.94 164.8 165.

22、59 165.21x = 78.5960 75.4500 90.4870 83.8480y = 164.9400 164.8000 165.5900 165.2100 polyfit(x,y,1)ans =0.0528 160.8006x13 x=231.39 198.48 212.64 190.55y=167.43 165.71 166.45 165.23x = 231.3900 198.4800 212.6400 190.5500y = 167.4300 165.7100 166.4500 165.2300 polyfit(x,y,1)ans =0.0535 155.0693x14x=75

23、.857 65.958 87.258 97.824y=164.23 163.04 165.54 166.88x = 75.8570 65.9580 87.2580 97.8240y = 164.2300 163.0400 165.5400 166.8800 polyfit(x,y,1)ans = 0.1199 155.1246x15 x=150.71 141.58 132.37 156.93y=164.07 164.27 164.57 163.89x = 150.7100 141.5800 132.3700 156.9300y = 164.0700 164.2700 164.5700 163.8900 polyfit(x,y,1)ans = -0.