潮流計(jì)算課程設(shè)計(jì)--基于matlab的電力系統(tǒng)潮流計(jì)算

1、<p><b>　　課程設(shè)計(jì)論文</b></p><p>　　基于MATLAB的電力系統(tǒng)潮流計(jì)算</p><p><b>　　學(xué)院:電氣工程學(xué)院</b></p><p>　　專業(yè):電氣工程及自動化</p><p><b>　　班級:</b></p>&l

2、t;p><b>　　姓名:  </b></p><p><b>　　學(xué)號:</b></p><p><b>　　目錄：</b></p><p>　　一、電力系統(tǒng)潮流計(jì)算概述4</p><p>　　二、潮流計(jì)算方法概述4</p><p>

3、;<b>　　2．1節(jié)點(diǎn)分類4</b></p><p>　　2.2牛頓—拉夫遜法概要5</p><p>　　2.2.1牛頓—拉夫遜法迭代原理5</p><p>　　2．2.1 牛頓法的框圖及求解過程6</p><p>　　三、課程設(shè)計(jì)任務(wù)8</p><p>　　3.1題目原始數(shù)據(jù)8&l

4、t;/p><p>　　3.2課程設(shè)計(jì)要求9</p><p><b>　　四、初步分析9</b></p><p>　　4.1節(jié)點(diǎn)設(shè)置及分類9</p><p>　　4.2 參數(shù)計(jì)算10</p><p>　　4.3等值電路繪制11</p><p><b>　　五、

5、潮流計(jì)算12</b></p><p>　　5.1 給定負(fù)荷下的潮流計(jì)算12</p><p>　　5.1.1 B1/B2矩陣的形成12</p><p>　　5.1.2 潮流調(diào)整及分析12</p><p>　　5.1.3 仿真比較14</p><p>　　5.2 變電所負(fù)荷變化時(shí)的潮流計(jì)算15<

6、;/p><p>　　5.2.1 4個(gè)變電所負(fù)荷同時(shí)以2%的比例增大15</p><p>　　5.2.2 4個(gè)變電所負(fù)荷成比例2%下降17</p><p>　　5.2.3 1、4號負(fù)荷下降，2、3號負(fù)荷上升19</p><p>　　5.2.4 仿真比較20</p><p>　　5.3斷線潮流計(jì)算22</

7、p><p>　　5.5.1 斷開1、3節(jié)點(diǎn)間的一條支路22</p><p>　　5.3.2 斷開1、5支路的一條線23</p><p>　　5.3.3 斷開1、7支路的一條線25</p><p>　　5.3.4 斷開2、9支路的一條線26</p><p>　　5.3.5 斷開7、9支路的一條線28</p&

8、gt;<p>　　5.3.6 仿真比較29</p><p>　　六、潮流分析總結(jié)31</p><p><b>　　七、心得體會33</b></p><p><b>　　參考文獻(xiàn)：33</b></p><p><b>　　附錄：34</b></p&

9、gt;<p><b>　　潮流計(jì)算課程設(shè)計(jì)</b></p><p><b>　　摘要</b></p><p>　　電力系統(tǒng)潮流計(jì)算是研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的一種計(jì)算，它根據(jù)給定的運(yùn)行條件及系統(tǒng)接線情況確定整個(gè)電力系統(tǒng)各部分的運(yùn)行狀態(tài)：各母線的電壓。各元件中流過的功率，系統(tǒng)的功率損耗等等。在電力系統(tǒng)規(guī)劃的設(shè)計(jì)和現(xiàn)有電力系統(tǒng)運(yùn)行方式的

10、研究中，都需要利用潮流計(jì)算來定量的分析比較供電方案或運(yùn)行方式的合理性、可靠性和經(jīng)濟(jì)性。此外，電力系統(tǒng)的潮流計(jì)算也是計(jì)算機(jī)系統(tǒng)動態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)，所以潮流計(jì)算是研究電力系統(tǒng)的一種重要的計(jì)算。</p><p>　　潮流計(jì)算在數(shù)學(xué)上是多元非線性方程組的求解問題，求解的方法有很多種，牛頓—拉夫遜法是數(shù)學(xué)上解非線性方程組的有效方法，有較好的收斂性。本文應(yīng)用了電力系統(tǒng)潮流計(jì)算仿真軟件DDTR與利用程序計(jì)算的結(jié)果進(jìn)行比較

11、，使計(jì)算的結(jié)果更加準(zhǔn)確。利用成形的程序?qū)ο到y(tǒng)中出現(xiàn)的各種情況，例如負(fù)荷的變化以及線路上所發(fā)生的變化進(jìn)行計(jì)算，并對母線上不滿足范圍的數(shù)據(jù)進(jìn)行調(diào)控，使得系統(tǒng)處于一個(gè)較穩(wěn)定的狀態(tài)。</p><p>　　關(guān)鍵詞：牛頓—拉夫遜法 MATLAB DDRTS 潮流計(jì)算 </p><p>　　一、電力系統(tǒng)潮流計(jì)算概述</p><p>　　在電力系統(tǒng)的正常運(yùn)行中，隨著用電負(fù)荷

12、的變化和系統(tǒng)運(yùn)行方式的改變，網(wǎng)絡(luò)中的損耗也將發(fā)生變化。要嚴(yán)格保證所有的用戶在任何時(shí)刻都有額定的電壓是不可能的，因此系統(tǒng)運(yùn)行中個(gè)節(jié)點(diǎn)出現(xiàn)電壓的偏移是不可避免的。為了保證電力系統(tǒng)的穩(wěn)定運(yùn)行，要進(jìn)行潮流調(diào)節(jié)。隨著電力系統(tǒng)及在線應(yīng)用的發(fā)展，計(jì)算機(jī)網(wǎng)絡(luò)已經(jīng)形成，為電力系統(tǒng)的潮流計(jì)算提供了物質(zhì)基礎(chǔ)。它的發(fā)展主要圍繞這樣幾個(gè)方面：計(jì)算方法的收斂性、可靠性；計(jì)算速度的快速性；對計(jì)算機(jī)存儲容量的要求以及計(jì)算的方便、靈活等。 </p>

13、<p>　　牛頓-拉夫遜法是數(shù)學(xué)中解決非線性方程式的典型方法，有較好的收斂性。。在解決電力系統(tǒng)潮流計(jì)算問題時(shí)，是以導(dǎo)納矩陣為基礎(chǔ)的，只要我們能在迭代過程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高潮流計(jì)算的效率。</p><p>　　二、潮流計(jì)算方法概述</p><p><b>　　2．1節(jié)點(diǎn)分類</b></p><p>　

14、　常規(guī)的電力系統(tǒng)潮流計(jì)算中一般具有三種類型的節(jié)點(diǎn)：PQ、PV及平衡節(jié)點(diǎn)。一個(gè)節(jié)點(diǎn)有四個(gè)變量，即注入有功功率、注入無功功率，電壓大小及相角。常規(guī)的潮流計(jì)算一般給定其中的二個(gè)變量：PQ節(jié)點(diǎn)（注入有功功率及無功功率），PV節(jié)點(diǎn)（注入有功功率及電壓的大小），平衡節(jié)點(diǎn)（電壓的大小及相角）。</p><p>　　第一類稱PQ節(jié)點(diǎn)：等值負(fù)荷功率PLi、QLi和等值電源功率PGi、QGi是給定的，從而注入功率P、Q是給定的，待求

15、的則是節(jié)點(diǎn)電壓的大小Ui和相位角。屬于這類節(jié)點(diǎn)的有按給定有功、無功率發(fā)電的發(fā)電廠母線和沒有其他電源的變電所母線。 </p><p>　　第二類稱PV節(jié)點(diǎn)：等值負(fù)荷和等值電源的有功功率PLi、PGi是給定的，從而注入有功功率Pi是給定的。等值負(fù)荷的無功功率QLi和節(jié)點(diǎn)電壓的大小Ui也是給定的。待求的則是等值電源的無功功率QGi，從而注入無功功率Qi和節(jié)點(diǎn)電壓的相位角。有一定無功功率儲備的發(fā)電廠和有一定無功

16、功率電源的變電所母線都可以作為PV節(jié)點(diǎn)。</p><p>　　第三類平衡節(jié)點(diǎn)：潮流計(jì)算時(shí)一般只設(shè)一個(gè)平衡節(jié)點(diǎn)。等值負(fù)荷功率PLs、QLs是給定的，節(jié)點(diǎn)電壓的大小和相位也是給定的。擔(dān)負(fù)調(diào)整系統(tǒng)頻率任務(wù)的發(fā)電廠母線往往被選作為平衡節(jié)點(diǎn)。</p><p>　　2.2牛頓—拉夫遜法概要</p><p>　　2.2.1牛頓—拉夫遜法迭代原理</p><p&

17、gt;　　已知一個(gè)變量X函數(shù)為： ，由適當(dāng)?shù)慕浦党霭l(fā)，根據(jù)：</p><p>　　反復(fù)進(jìn)行計(jì)算，當(dāng)滿足適當(dāng)?shù)氖諗織l件就是上面方程的根。這樣的方法就是所謂的牛頓—拉夫遜法。</p><p>　　這一方法還可以做下面的解釋，設(shè)第次迭代得到的解語真值之差，即的誤差為時(shí)，則：</p><p>　　把在附近對用泰勒級數(shù)展開</p><p><b

18、>　　上式省略去以后部分</b></p><p>　　的誤差可以近似由上式計(jì)算出來。</p><p>　　比較兩式，可以看出牛頓—拉夫遜法的休整量和的誤差的一次項(xiàng)相等。</p><p>　　用同樣的方法考慮，給出個(gè)變量的個(gè)方程：</p><p>　　對其近似解得修正量可以通過解下邊的方程來確定：</p><

19、;p>　　式中等號右邊的矩陣都是對于的值。這一矩陣稱為雅可比（JACOBI）矩陣。按上述得到的修正向量后，得到如下關(guān)系</p><p>　　這比更接近真實(shí)值。這一步在收斂到希望的值以前重復(fù)進(jìn)行，一般要反復(fù)計(jì)算滿足</p><p>　　為預(yù)先規(guī)定的小正數(shù)，是第n次迭代的近似值。</p><p>　　2．2.1 牛頓法的框圖及求解過程</p><

20、;p>　?。?）用牛頓法計(jì)算潮流時(shí)，有以下的步驟：</p><p>　?、俳o這各節(jié)點(diǎn)電壓初始值；</p><p>　?、趯⒁陨想妷撼跏贾荡牍?，求修正方程的常數(shù)項(xiàng)向量 ；</p><p>　?、蹖㈦妷撼跏贾翟趲肷鲜龉剑蟪鲂拚匠讨邢禂?shù)矩陣的各元素。</p><p><b>　?、芙庑拚匠淌?；&l

21、t;/b></p><p>　?、菪拚鞴?jié)點(diǎn)電壓，；</p><p>　?、迣?，在帶入方程式，求出；</p><p><b>　?、邫z驗(yàn)是否收斂，即</b></p><p>　　如果收斂，迭代到此結(jié)束，進(jìn)一步計(jì)算各線路潮流和平衡節(jié)點(diǎn)功率，并打印</p><p>　　輸出結(jié)果。如果不收斂，轉(zhuǎn)回

22、②進(jìn)行下次迭代計(jì)算，直到收斂為止。</p><p>　?。?） 牛頓—拉夫遜法計(jì)算程序框圖</p><p><b>　　三、課程設(shè)計(jì)任務(wù)</b></p><p><b>　　3.1題目原始數(shù)據(jù)</b></p><p>　　1.系統(tǒng)圖：兩個(gè)發(fā)電廠分別通過變壓器和輸電線路與四個(gè)變電所相連。</p&

23、gt;<p><b>　　2、發(fā)電廠資料：</b></p><p>　　母線1和2為發(fā)電廠高壓母線，發(fā)電廠一總裝機(jī)容量為（ 300MW ），母線3為機(jī)壓母線，機(jī)壓母線上裝機(jī)容量（100MW），最大負(fù)荷和最小負(fù)荷為50MW和20MW；發(fā)電廠二總裝機(jī)容量（200MW）。</p><p><b>　　3、變電所資料：</b></p&

24、gt;<p>　　（一）變電所1、2、3、4低壓母線的電壓等級分別為：35KV 10KV 35KV 10KV（二）變電所的負(fù)荷分別為：60MW 50MW 40MW 60MW </p><p>　　每個(gè)變電所的功率因數(shù)均為cosφ=0.85；</p><p>　　變電所1和變電所3分別配有兩臺容量為75MVA的變壓器，短路損耗414KW，短路電壓（%

25、）=16.7；變電所2和變電所4分別配有兩臺容量為63MVA的變壓器，短路損耗為245KW，短路電壓（%）=10.5；</p><p><b>　　4、輸電線路資料：</b></p><p>　　發(fā)電廠和變電所之間的輸電線路的電壓等級及長度標(biāo)于圖中，單位長度的電阻為，單位長度的電抗為，單位長度的電納為。</p><p>　　3.2課程設(shè)計(jì)要求

26、</p><p>　　對給定的網(wǎng)絡(luò)查找潮流計(jì)算所需的各元件等值參數(shù)，畫出等值電路圖。</p><p>　　輸入各支路數(shù)據(jù)，各節(jié)點(diǎn)數(shù)據(jù)利用給定的程序進(jìn)行在變電所在某一負(fù)荷情況下的潮流計(jì)算，并對計(jì)算結(jié)果進(jìn)行分析。</p><p>　　跟隨變電所負(fù)荷按一定比例發(fā)生變化，進(jìn)行潮流計(jì)算分析。</p><p>　　4個(gè)變電所的負(fù)荷同時(shí)以2%的比例增大；&

27、lt;/p><p>　　4個(gè)變電所的負(fù)荷同時(shí)以2%的比例下降</p><p>　　1和4號變電所的負(fù)荷同時(shí)以2%的比例下降，而2和3號變電所的負(fù)荷同時(shí)以2%的比例上升；</p><p>　　在不同的負(fù)荷情況下，分析潮流計(jì)算的結(jié)果，如果各母線電壓不滿足要求，進(jìn)行電壓的調(diào)整。（變電所低壓母線電壓10KV要求調(diào)整范圍在9.5-10.5之間；電壓35KV要求調(diào)整范圍在35-36

28、之間）</p><p>　　輪流斷開支路雙回線中的一條，分析潮流的分布。（幾條支路斷幾次）</p><p>　　利用DDRTS軟件，進(jìn)行繪制系統(tǒng)圖進(jìn)行上述各種情況潮流的分析，并進(jìn)行結(jié)果的比較。 </p><p>　　最終形成課程設(shè)計(jì)成品說明書。</p><p><b>　　四、初步分析</b></p>&l

29、t;p>　　4.1節(jié)點(diǎn)設(shè)置及分類</p><p>　　根據(jù)系統(tǒng)圖可知此系統(tǒng)為兩端供電網(wǎng)路，為了潮流計(jì)算的完整性，將母線1，2設(shè)為節(jié)點(diǎn)1，2，將變電所1、2、3、4的高低壓側(cè)分別設(shè)為節(jié)點(diǎn)3、4、5、6、7、8、9、10。并且，將節(jié)點(diǎn)1設(shè)為平衡節(jié)點(diǎn)，將節(jié)點(diǎn)2設(shè)為PV節(jié)點(diǎn)，其余節(jié)點(diǎn)設(shè)為PQ節(jié)點(diǎn)。則網(wǎng)絡(luò)共具有10個(gè)節(jié)點(diǎn)和9條支路。</p><p><b>　　4.2 參數(shù)計(jì)算<

30、;/b></p><p>　　根據(jù)題目原始資料，計(jì)算發(fā)電廠、變壓器及線路的有名值參數(shù)。</p><p>　?。?）運(yùn)用下列公式計(jì)算變壓器低壓側(cè)參數(shù)：</p><p><b>　　（2）計(jì)算線路參數(shù)</b></p><p>　　將參數(shù)整理成如下表格：</p><p>　?。?）變電所負(fù)荷分別為

31、：</p><p>　　變電所1 =60+j37.185 </p><p>　　變電所2 =50+j30.987</p><p>　　變電所3 =40+j24.790 </p><p>　　變電所4 =60+j37.185</p><p>　　(4）計(jì)算變壓器分接頭變比</p>

32、<p>　　變壓器有5個(gè)抽頭，電壓調(diào)節(jié)范圍為2*2.5%， 對應(yīng)的分接頭開始時(shí)設(shè)變壓器高壓側(cè)接主接頭,降壓變壓器5個(gè)分接頭時(shí)變比K以備調(diào)壓時(shí)選用。</p><p><b>　　4.3等值電路繪制</b></p><p>　　根據(jù)題目給定的參數(shù)和要求，手繪形成了網(wǎng)絡(luò)的等值電路圖,并且在圖中標(biāo)注出了各參數(shù)和節(jié)點(diǎn)編號。</p><p>

33、　　圖4.3 等值電路圖</p><p><b>　　4.4計(jì)算方法</b></p><p>　　利用牛頓拉夫遜法進(jìn)行求解，用MATLAB軟件編程，可以求解系統(tǒng)潮流分</p><p>　　布。根據(jù)題目的不同要求對參數(shù)進(jìn)行調(diào)整，通過調(diào)節(jié)變壓器變比和發(fā)電廠的電壓，求解出合理的潮流分布，最后用DDTRS進(jìn)行潮流分析，將兩者進(jìn)行比較。</p>

34、;<p><b>　　五、潮流計(jì)算</b></p><p>　　5.1 給定負(fù)荷下的潮流計(jì)算</p><p>　　5.1.1 B1/B2矩陣的形成</p><p>　　根據(jù)所求參數(shù)，以及B1、B2矩陣的含義，列寫B(tài)1、B2矩陣如下： </p><p>　　0 0 235

35、 0 1</p><p>　　200 0 235 0 3</p><p>　　0 0 220 0 2</p><p>　　0 60+j37.18466 35 0 2</p><p>　　0 0

36、220 0 2</p><p>　　0 50+j30.98722 10 0 2</p><p>　　0 0 220 0 2</p><p>　　0 40+j24.78977 35 0 2</p><p>　　0 0

37、 220 0 2</p><p>　　0 60+j37.18466 10 0 2</p><p>　　5.1.2 潮流調(diào)整及分析</p><p>　　將編好的的程序(程序見附錄)輸入MATLAB中運(yùn)行，可以得到如下運(yùn)行結(jié)果:</p><p>　　由表觀察到節(jié)點(diǎn)4,6,8,10的電

38、壓值都較正常范圍偏高，因此調(diào)節(jié)發(fā)電廠電壓和變壓器分接頭，經(jīng)過多次調(diào)整，使節(jié)點(diǎn)電壓均滿足要求，調(diào)節(jié)結(jié)果如下：</p><p><b>　　調(diào)整過程：</b></p><p><b>　　調(diào)整結(jié)果：</b></p><p>　　由上表觀察到，先后進(jìn)行發(fā)電機(jī)機(jī)端電壓調(diào)節(jié)和變壓器分接頭的調(diào)節(jié)使得節(jié)點(diǎn)4,6,8,10的節(jié)點(diǎn)電壓均在允

39、許范圍內(nèi)。</p><p>　　調(diào)節(jié)前后各支路損耗對比如下：</p><p>　　由此可以看出，隨著發(fā)電機(jī)端電壓的降低，網(wǎng)絡(luò)有功損耗逐漸增大，符合實(shí)際情況，而最終調(diào)節(jié)使得電壓在規(guī)定范圍內(nèi)。調(diào)整后的各支路的功率如下表所示：</p><p>　　支路的末端功率加上支路損耗等于該支路的首段功率，符合實(shí)際情況。</p><p>　　且得到各節(jié)點(diǎn)的功率

40、 (節(jié)點(diǎn)號從小到大排列)：</p><p>　　0.2371+1.0657i 2.0000-0.4373i -0.0000+0.0000i -0.6000-0.3718i -0.0000-0.0000i -0.5000-0.3099i -0.0000+0.0000i -0.4000-0.2479i 0.0000+0.0000i -0.6000-0.3718i</p><p

41、>　　5.1.3 仿真比較</p><p>　　由MATLAB得到電壓、功角和具體的潮流分布圖，調(diào)整好以后的數(shù)據(jù)如圖：</p><p>　　在DDTRS中繪制出電路圖，仿真并進(jìn)行調(diào)節(jié)以后可以得到潮流分布圖：</p><p>　　在用DDTRS繪制的系統(tǒng)潮流分布圖中，各個(gè)節(jié)點(diǎn)電壓在調(diào)試之后均能滿足要求，均達(dá)到系統(tǒng)電壓等級標(biāo)準(zhǔn)。</p><p&

42、gt;　　將PSAT仿真得到的潮流圖中個(gè)點(diǎn)電壓的有名值與MATLAB計(jì)算得到的各節(jié)點(diǎn)電壓進(jìn)行比較：</p><p>　　由上述比較可觀察到: 應(yīng)用DDTRS仿真和MATLAB計(jì)算的到的結(jié)果大致相同，可以認(rèn)為此潮流計(jì)算是正確的。</p><p>　　5.2 變電所負(fù)荷變化時(shí)的潮流計(jì)算 </p><p>　　5.2.1 4個(gè)變電所負(fù)荷同時(shí)以2%的比例增大</p&

43、gt;<p>　　由于各個(gè)變電所的負(fù)荷同時(shí)以2%的比例增大，所以此時(shí)的各負(fù)荷變?yōu)椋?變電所1 =61.2+j37.92835 </p><p>　　變電所2 = 51+j31.60696</p><p>　　變電所3 =40.8+j25.28557</p><p>　　變電所4 =61.2+j37.

44、92835</p><p>　　相應(yīng)的，B2矩陣變?yōu)椋?lt;/p><p>　　0 0 235 0 1</p><p>　　200 0 235 0 3</p><p>　　0 0 220 0 2</p><p>

45、;　　061.2+j37.92835 35 0 2</p><p>　　0 0 220 0 2</p><p>　　0 51+j31.60696 10 0 2</p><p>　　0 0 220 0 2</p><p>　　0

46、40.8+j25.28557 35 0 2</p><p>　　0 0 220 0 2</p><p>　　0 61.2+j37.92835 10 0 2</p><p>　　將原程序中的B2矩陣替換后，運(yùn)行并進(jìn)行調(diào)節(jié)，經(jīng)過調(diào)節(jié)后達(dá)到要求，對應(yīng)的調(diào)節(jié)方法和結(jié)果見下表：</p>

47、<p><b>　　調(diào)節(jié)方法：</b></p><p><b>　　調(diào)整結(jié)果： </b></p><p>　　由上表觀察到，調(diào)節(jié)后的電壓均在題目允許的范圍內(nèi)，對線路損耗進(jìn)行分析，統(tǒng)計(jì)每次調(diào)整后各個(gè)支路的有功損耗，記錄于下表：</p><p>　　由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目

48、要求，可認(rèn)為是合理的。具體的潮流分布如下：</p><p>　　各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號從小到大排列)：</p><p>　　0.278 +1.0872i 2.0000-0.4257i 0.0000-0.0000i -0.6120-0.3793i 0.0000+0.0000i -0.5100 0.3161i -0.000-0.0000i -0.4080-0.25

49、29i 0.0000+0.0000i -0.6120-0.3793i</p><p>　　5.2.2 4個(gè)變電所負(fù)荷成比例2%下降</p><p>　　由于各個(gè)變電所的負(fù)荷同時(shí)以2%的比例下降，所以此時(shí)的各負(fù)荷變?yōu)椋?變電所1 =58.8+j36.44097 </p><p>　　變電所2 = 49+j30.36

50、747</p><p>　　變電所3 =39.2+j24.29398</p><p>　　變電所4 =58.8+j36.44097</p><p>　　相應(yīng)的，B2矩陣變?yōu)椋?lt;/p><p>　　0 0 235 0 1</p><p>　　200 0

51、 235 0 3</p><p>　　0 0 220 0 2</p><p>　　058.8+j36.44097 35 0 2</p><p>　　0 0 220 0 2</p><p>　　0 49+j30.36747

52、 10 0 2</p><p>　　0 0 220 0 2</p><p>　　039.2+j24.29398 35 0 2</p><p>　　0 0 220 0 2</p><p>　　0 58.8+j36.44097

53、10 0 2</p><p>　　將程序中的B2矩陣替換后，運(yùn)行并進(jìn)行調(diào)節(jié)，此系統(tǒng)經(jīng)過調(diào)節(jié)后達(dá)到要求，對應(yīng)的調(diào)節(jié)方法和結(jié)果見下表：</p><p><b>　　調(diào)節(jié)方法：</b></p><p><b>　　調(diào)整結(jié)果： </b></p><p>　　由上表觀察到，調(diào)節(jié)后的電壓均在題目允許的

54、范圍內(nèi)，對線路損耗進(jìn)行分析，統(tǒng)計(jì)每次調(diào)整后各個(gè)支路的有功損耗，記錄于下表：</p><p>　　由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認(rèn)為是合理的。具體的潮流分布如下：</p><p>　　各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號從小到大排列)：</p><p>　　0.1954+1.0444i 2.0000-0.4488i -0.0000+0.

55、0000i -0.5880-0.3644i 0.0000+0.0000i -0.4900-0.3037i -0.0000+0.0000i -0.3920-0.2429i -0.0000-0.0000i - 0.5880-0.3644i</p><p>　　5.2.3 1、4號變電所負(fù)荷成2%下降，2、3號變電負(fù)荷成2%上升</p><p><b>　　此時(shí)

56、的各負(fù)荷變?yōu)椋?lt;/b></p><p>　　變電所1 =58.8+j36.44097</p><p>　　變電所2 = 51+j31.60696</p><p>　　變電所3 =40.8+j25.28557</p><p>　　變電所4 =58.8+j36.44097</p><p>

57、;　　相應(yīng)的，B2矩陣變?yōu)椋?lt;/p><p>　　0 0 235 0 1</p><p>　　200 0 235 0 3</p><p>　　0 0 220 0 2</p><p>　　058.8+j36.44097 3

58、5 0 2</p><p>　　0 0 220 0 2</p><p>　　0 51+j31.60696 10 0 2</p><p>　　0 0 220 0 2</p><p>　　040.8+j25.28557 35

59、0 2</p><p>　　0 0 220 0 2</p><p>　　0 58.8+j36.44097 10 0 2</p><p>　　將程序中的B2矩陣替換后，運(yùn)行并進(jìn)行調(diào)節(jié)，此系統(tǒng)經(jīng)過調(diào)節(jié)后達(dá)到要求，對應(yīng)的調(diào)節(jié)方法和結(jié)果見下表：</p><p><b>　　調(diào)

60、節(jié)方法：</b></p><p><b>　　調(diào)整結(jié)果： </b></p><p>　　由上表觀察到，調(diào)節(jié)后的電壓均在題目允許的范圍內(nèi)，對線路損耗進(jìn)行分析，統(tǒng)計(jì)每次調(diào)整后各個(gè)支路的有功損耗，記錄于下表：</p><p>　　由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認(rèn)為是合理的。具體的潮流分布如下：&l

61、t;/p><p>　　各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號從小到大排列)：</p><p>　　0.2314+1.0657i 2.0000-0.4427i -0.0000+0.0000i -0.5880-0.3644i 0.0000+0.0000i -0.5100-0.3161i -0.0000+0.0000i -0.4080-0.2529i -0.0000-0.0000i

62、 -0.5880-0.3644i</p><p>　　5.2.4 仿真比較</p><p>　　負(fù)荷每次發(fā)生變化時(shí)相應(yīng)的修改DDTRS中電路的參數(shù)，仿真并進(jìn)行調(diào)節(jié)，可以得到相應(yīng)的潮流分布圖，現(xiàn)于由MATLAB得到的潮流分布數(shù)據(jù)圖比較（以4個(gè)變電所負(fù)荷增加2%為例）：</p><p>　　由MATLAB得到的數(shù)據(jù)圖：</p><p>　　由DD

63、TRS得到的潮流分布圖：</p><p>　　在用DDTRS繪制的系統(tǒng)潮流分布圖中，各個(gè)節(jié)點(diǎn)電壓在調(diào)試之后均能滿足要求，均達(dá)到系統(tǒng)電壓等級標(biāo)準(zhǔn)。</p><p>　　將PSAT仿真得到的潮流圖中個(gè)點(diǎn)電壓的有名值與MATLAB計(jì)算得到的各節(jié)點(diǎn)電壓進(jìn)行比較：</p><p>　　由上述比較可觀察到: 應(yīng)用DDTRS仿真和MATLAB計(jì)算的到的結(jié)果大致相同，可以認(rèn)為此潮流

64、計(jì)算是正確的。</p><p><b>　　5.3斷線潮流計(jì)算</b></p><p>　　5.5.1 斷開1、3節(jié)點(diǎn)間的一條支路</p><p>　　當(dāng)斷開1、3節(jié)點(diǎn)間的雙回線中的一條時(shí)時(shí)，B1矩陣的改變，B1矩陣變?yōu)椋?lt;/p><p>　　1 3 17+j40.2 j2.780e-4

65、 1 0 0</p><p>　　1 5 4.25+j10.05 j2.780e-4 1 0 0</p><p>　　1 7 7.65+j18.09 j5.004e-4 1 0 0</p><p>　　2 9 8.5+j20.1 j5

66、.560e-4 1 0 0</p><p>　　3 4 0.05455+j1.65024 0 231/38.5 1 1 </p><p>　　5 6 3.375e-3+j0.10083 0 231/11 1 1</p><p>　　7 8 0.

67、05455+j1.65024 0 231/38.5 1 1</p><p>　　7 9 6.8+j16.08 j4.448e-4 1 0 0</p><p>　　9 10 3.375e-3+j0.10083 0 231/11 1 1</p><p>　

68、　用修改后的B1矩陣帶入程序，用MATLAB進(jìn)行潮流計(jì)算并進(jìn)行調(diào)整，對應(yīng)的調(diào)節(jié)方法和結(jié)果見下表：</p><p><b>　　調(diào)節(jié)方法：</b></p><p><b>　　調(diào)整結(jié)果： </b></p><p>　　由上表觀察到，調(diào)節(jié)后的電壓均在題目允許的范圍內(nèi)，對線路損耗進(jìn)行分析，統(tǒng)計(jì)每次調(diào)整后各個(gè)支路的有功損耗，記錄

69、于下表：</p><p>　　由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認(rèn)為是合理的。具體的潮流分布如下：</p><p>　　各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號從小到大排列)：</p><p>　　0.2480+1.2342i 2.0000-0.4373i -0.0000-0.0000i -0.6000-0.3718i -0.0000-

70、0.0000i -0.500 -0.3099i -0.0000+0.0000i -0.4000-0.2479i 0.0000+0.0000i -0.6000-0.3718i</p><p>　　5.3.2 斷開1、5支路的一條線</p><p>　　當(dāng)斷開1、5節(jié)點(diǎn)間的雙回線中的一條時(shí)時(shí)，B1矩陣的改變，B1矩陣變?yōu)椋?lt;/p><p>　　1

71、 3 8.5+j20.1 j5.560e-4 1 0 0</p><p>　　1 5 8.5+j20.1 j1.390e-4 1 0 0</p><p>　　1 7 7.65+j18.09 j5.004e-4 1 0 0</p>

72、<p>　　2 9 8.5+j20.1 j5.560e-4 1 0 0</p><p>　　3 4 0.05455+j1.65024 0 231/38.5 1 1 </p><p>　　5 6 3.375e-3+j0.10083 0 231/11 1

73、 1</p><p>　　7 8 0.05455+j1.65024 0 231/38.5 1 1</p><p>　　7 9 6.8+j16.08 j4.448e-4 1 0 0</p><p>　　910 3.375e-3+j0.10083 0 231/1

74、1 1 1</p><p>　　用修改后的B2矩陣帶入程序，用MATLAB進(jìn)行潮流計(jì)算并進(jìn)行調(diào)整，對應(yīng)的調(diào)節(jié)方法和結(jié)果見下表：</p><p><b>　　調(diào)節(jié)方法：</b></p><p><b>　　調(diào)整結(jié)果： </b></p><p>　　由上表觀察到，調(diào)節(jié)后的電壓均在題目允許的范圍

75、內(nèi)，對線路損耗進(jìn)行分析，統(tǒng)計(jì)每次調(diào)整后各個(gè)支路的有功損耗，記錄于下表：</p><p>　　由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認(rèn)為是合理的。具體的潮流分布如下：</p><p>　　各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號從小到大排列)：</p><p>　　0.2405+1.1448i 2.0000-0.4373i -0.0000+0.

76、0000i -0.6000-0.3718i -0.0000+0.0000i</p><p>　　-0.5000-0.3099i -0.0000+0.0000i -0.4000-0.2479i 0.0000+0.0000i -0.6000-0.3718i</p><p>　　5.3.3 斷開1、7支路的一條線</p><p>　　當(dāng)斷開1、7節(jié)點(diǎn)間的

77、雙回線中的一條時(shí)時(shí)，B1矩陣的改變，B1矩陣變</p><p>　　1 3 8.5+j20.1 j5.560e-4 1 00</p><p>　　1 5 4.25+j10.05 j2.780e-4 1 00</p><p>　　1 7 15.3+j36

78、.18 j2.502e-4 1 00</p><p>　　2 9 8.5+j20.1 j5.560e-4 1 00</p><p>　　3 4 0.05455+j1.65024 0 231/38.5 1 1 </p><p>　　5 6

79、 3.375e-3+j0.10083 0 231/11 1 1</p><p>　　7 8 0.05455+j1.65024 0 231/38.5 1 1</p><p>　　7 9 6.8+j16.08 j4.448e-4 1 00</p><

80、p>　　910 3.375e-3+j0.10083 0 231/11 1 1</p><p>　　用修改后的B2矩陣帶入程序，用MATLAB進(jìn)行潮流計(jì)算并進(jìn)行調(diào)整，對應(yīng)的調(diào)節(jié)方法和結(jié)果見下表：</p><p><b>　　調(diào)節(jié)方法：</b></p><p><b>　　調(diào)整結(jié)果：</b&

81、gt;</p><p>　　由上表觀察到，調(diào)節(jié)后的電壓均在題目允許的范圍內(nèi)，對線路損耗進(jìn)行分析，統(tǒng)計(jì)每次調(diào)整后各個(gè)支路的有功損耗，記錄于下表：</p><p>　　由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認(rèn)為是合理的。具體的潮流分布如下：</p><p>　　各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號從小到大排列)：</p><p&g

82、t;　　0.2530+1.0927i 2.0000-0.2942i -0.0000+0.0000i -0.6000-0.3718i -0.0000-0.0000i</p><p>　　-0.5000-0.3099i 0.0000+0.0000i -0.4000-0.2479i -0.0000-0.0000i -0.6000-0.3718i</p><p>　　5.3.4

83、 斷開2、9支路的一條線</p><p>　　當(dāng)斷開2、9節(jié)點(diǎn)間的雙回線中的一條時(shí)時(shí)，B1矩陣的改變，B1矩陣變?yōu)?lt;/p><p>　　1 3 8.5+j20.1 j5.560e-4 1 0 0</p><p>　　1 5 4.25+j10.05 j2.780e-4 1

84、 0 0</p><p>　　1 7 7.65+j18.09 j5.004e-4 1 0 0</p><p>　　2 9 17+j40.2 j2.780e-4 1 0 0</p><p>　　3 4 0.05455+j1.65024 0 2

85、31/38.5 1 1 </p><p>　　5 6 3.375e-3+j0.10083 0 231/11 1 1</p><p>　　7 8 0.05455+j1.65024 0 231/38.5 1 1</p><p>　　7 9 6.8+j16.08

86、 j4.448e-4 1 0 0</p><p>　　910 3.375e-3+j0.10083 0 231/11 1 1</p><p>　　用修改后的B2矩陣帶入程序，用MATLAB進(jìn)行潮流計(jì)算并進(jìn)行調(diào)整，對應(yīng)的調(diào)節(jié)方法和結(jié)果見下表：</p><p><b>　　調(diào)節(jié)方法：</b>&l

87、t;/p><p><b>　　調(diào)整結(jié)果：</b></p><p>　　由上表觀察到，調(diào)節(jié)后的電壓均在題目允許的范圍內(nèi)，對線路損耗進(jìn)行分析，統(tǒng)計(jì)每次調(diào)整后各個(gè)支路的有功損耗，記錄于下表：</p><p>　　由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認(rèn)為是合理的。具體的潮流分布如下：</p><p>

88、;　　各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號從小到大排列)：</p><p>　　0.3107+1.2647i 2.000 -0.3130i -0.0000+0.0000i -0.6000-0.3718i -0.0000-0.0000i</p><p>　　-0.5000-0.3099i -0.0000+0.0000i -0.4000-0.2479i -0.0000+0.0000

89、i -0.6000-0.3718i</p><p>　　5.3.5 斷開7、9支路的一條線</p><p>　　當(dāng)斷開2、9節(jié)點(diǎn)間的雙回線中的一條時(shí)時(shí)，B1矩陣的改變，B1矩陣變?yōu)椋?lt;/p><p>　　1 3 8.5+j20.1 j5.560e-4 1 0 0</p><p>　　1

90、 5 4.25+j10.05 j2.780e-4 1 0 0</p><p>　　1 7 7.65+j18.09 j5.004e-4 1 0 0</p><p>　　2 9 8.5+j20.1 j5.560e-4 1 0 0</p>&l

91、t;p>　　3 4 0.05455+j1.65024 0 231/38.5 1 1 </p><p>　　5 6 3.375e-3+j0.10083 0 231/11 1 1</p><p>　　7 8 0.05455+j1.65024 0 231/38.5 1

92、1</p><p>　　7 9 13.6+j32.16 j2.224e-4 1 0 0</p><p>　　910 3.375e-3+j0.10083 0 231/11 1 1</p><p>　　用修改后的B2矩陣帶入程序，用MATLAB進(jìn)行潮流計(jì)算并進(jìn)行調(diào)整，對應(yīng)的調(diào)節(jié)方法和結(jié)果見下表

93、：</p><p><b>　　調(diào)節(jié)方法：</b></p><p><b>　　調(diào)整結(jié)果：</b></p><p>　　由上表觀察到，調(diào)節(jié)后的電壓均在題目允許的范圍內(nèi)，對線路損耗進(jìn)行分析，統(tǒng)計(jì)每次調(diào)整后各個(gè)支路的有功損耗，記錄于下表：</p><p>　　由電壓和功率損耗可以觀察到，最后一次的電壓和

94、有功損耗均符合題目要求，可認(rèn)為是合理的。具體的潮流分布如下：</p><p>　　各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號從小到大排列)：</p><p>　　0.2658+1.1276i 2.0000-0.3190i -0.0000+0.0000i -0.6000-0.3718i -0.0000-0.0000i</p><p>　　-0.5000-0.3099i

95、-0.0000+0.0000i -0.4000-0.2479i -0.0000-0.0000i -0.6000-0.3718i</p><p>　　5.3.6 仿真比較</p><p>　　線路每次發(fā)生變化時(shí)相應(yīng)的修改DDTRS中電路的線路參數(shù)，仿真并進(jìn)行調(diào)節(jié)，可以得到相應(yīng)的潮流分布圖，現(xiàn)與由MATLAB得到的潮流分布數(shù)據(jù)圖比較（以斷開2、9節(jié)點(diǎn)間的支路為例）：</p>

96、<p>　　由MATLAB得到潮流數(shù)據(jù)圖：</p><p>　　由DDTRS得到潮流分布圖：</p><p>　　在用DDTRS繪制的系統(tǒng)潮流分布圖中，各個(gè)節(jié)點(diǎn)電壓在調(diào)試之后均能滿足要求，均達(dá)到系統(tǒng)電壓等級標(biāo)準(zhǔn)。</p><p>　　將PSAT仿真得到的潮流圖中個(gè)點(diǎn)電壓的有名值與MATLAB計(jì)算得到的各節(jié)點(diǎn)電壓進(jìn)行比較。</p><

97、p>　　由上述比較可觀察到: 應(yīng)用DDTRS仿真和MATLAB計(jì)算的到的結(jié)果大致相同，可以認(rèn)為此潮流計(jì)算是正確的。</p><p><b>　　六、潮流分析總結(jié)</b></p><p>　　潮流計(jì)算是電力系統(tǒng)的最基本、最常用的分析計(jì)算。用以研究系統(tǒng)規(guī)劃和運(yùn)行中提出的各種問題。對規(guī)劃中的電力系統(tǒng)，通過潮流計(jì)算可以檢驗(yàn)所提出的電力系統(tǒng)規(guī)劃方案能否滿足各種運(yùn)行方式的要

98、求；對運(yùn)行中的電力系統(tǒng)，通過潮流計(jì)算可以預(yù)知各種負(fù)荷變化和網(wǎng)絡(luò)結(jié)構(gòu)的改變會不會危及系統(tǒng)的安全，系統(tǒng)中所有母線的電壓是否在允許的范圍以內(nèi)，系統(tǒng)中各種元件(線路、變壓器等)是否會出現(xiàn)過負(fù)荷，以及可能出現(xiàn)過負(fù)荷時(shí)應(yīng)事先采取哪些預(yù)防措施等。</p><p>　　電壓是衡量電力系統(tǒng)電能質(zhì)量的標(biāo)準(zhǔn)之一。電壓過高或過低，都將對人身及其用電設(shè)備產(chǎn)生重大的影響。保證用戶的電壓接近額定值是電力系統(tǒng)調(diào)度的基本任務(wù)之一。當(dāng)系統(tǒng)的電壓偏離

99、允許值時(shí)，電力系統(tǒng)必須應(yīng)用電壓調(diào)節(jié)技術(shù)調(diào)節(jié)系統(tǒng)電壓的大小，使其維持在允許值范圍內(nèi)。</p><p>　　電力系統(tǒng)調(diào)壓措施中，常用的調(diào)壓手段有：借改變發(fā)電機(jī)機(jī)端電壓、改變變壓器的變比、無功補(bǔ)償設(shè)備的添加等。</p><p>　　1.改變發(fā)電機(jī)的機(jī)端電壓</p><p>　　對于電力系統(tǒng)，調(diào)節(jié)電壓首先要考慮改變機(jī)端電壓，其是最經(jīng)濟(jì)最直接的調(diào)壓措施。適用于出廠母線附近沒有

100、其他的用電設(shè)備網(wǎng)絡(luò)，調(diào)壓范圍在0.95—1.05Un之間。由于該調(diào)壓措施的調(diào)壓范圍比較窄及上述數(shù)據(jù)分析可知，各變電所的電壓均偏高，若降低機(jī)端電壓，在此時(shí)的負(fù)荷條件下可滿足個(gè)別變電所的電壓要求，但出廠母線電壓均降到220KV以下，是系統(tǒng)運(yùn)行在低于額定電壓等級狀況下,在系統(tǒng)的負(fù)荷發(fā)生改變時(shí),在某些幾點(diǎn)處的電壓會驟然下降,引起系統(tǒng)的癱瘓,造成大面積停電.</p><p>　　2.改變變壓器的變比</p>

101、<p>　　現(xiàn)代電力系統(tǒng)網(wǎng)絡(luò)中所常用的變壓器多為有載變壓器，能在有載條件下進(jìn)行調(diào)節(jié)，不需退出運(yùn)行就能改變變比。用變壓器調(diào)壓的前提是必須保證系統(tǒng)的無功功率充足。有上述的數(shù)據(jù)分析，當(dāng)調(diào)節(jié)機(jī)端電壓不能滿足各個(gè)變電所的電壓要求時(shí)，若個(gè)別變電所電壓偏高，則升高變壓器的變比；若個(gè)別變電所電壓偏低，則降低變壓器的變比，通?？梢赃_(dá)到調(diào)壓要求。</p><p>　　3.無功補(bǔ)償設(shè)備的添加</p><p

102、>　　無功補(bǔ)償設(shè)備的添加常見的有：串聯(lián)電抗器、并聯(lián)電容器、調(diào)相機(jī)等。對于一個(gè)網(wǎng)絡(luò)的電壓調(diào)整，常在變電所添加調(diào)相機(jī)、電容器等設(shè)備。雖然通過無功補(bǔ)償可以滿足要求，但其成本投入比較大，造價(jià)高。</p><p>　　此調(diào)壓的方案沒有充分考慮兩個(gè)變電所均衡承擔(dān)負(fù)荷的問題，導(dǎo)致發(fā)電廠二滿載運(yùn)行，而發(fā)電廠一功率過小，負(fù)荷分配不均從而導(dǎo)致網(wǎng)損略大，沒有充分契合系統(tǒng)運(yùn)行的經(jīng)濟(jì)性。</p><p>&l

103、t;b>　　七、心得體會</b></p><p>　　從本次電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)中，通過計(jì)算在不同的負(fù)荷下各支路節(jié)點(diǎn)潮流分布的變化，分析負(fù)荷變動時(shí)，對電壓，功率的影響，我們掌握了常用的調(diào)整電壓的措施，鍛煉了自己實(shí)際操作分析能力。使我在潮流計(jì)算、MATLAB的使用方面均有所提高，但也暴露出了一些問題；理論知識的儲備還是不足，對MATLAB的性能和特點(diǎn)還不能有一個(gè)全面的把握，對MATLAB的使用

104、不夠熟悉，相信通過以后的學(xué)習(xí)能彌補(bǔ)這些不足，從而達(dá)到一個(gè)新的層次。</p><p>　　在潮流計(jì)算的過程中，為了滿足調(diào)壓要求，需要一次一次的嘗試，首選的調(diào)壓措施是改變機(jī)端電壓，一般不能滿足所有的節(jié)點(diǎn)的電壓要求，這就需要調(diào)節(jié)個(gè)別的變壓器分接頭，這是兩種較為經(jīng)濟(jì)的調(diào)壓措施。這個(gè)方案充分利用了這兩種調(diào)壓措施，使電壓都滿足了要求，保證了調(diào)壓的經(jīng)濟(jì)性。在DDTRS仿真時(shí)，對參數(shù)的設(shè)置需要細(xì)心，有一點(diǎn)疏忽就會導(dǎo)致運(yùn)行結(jié)果不合

105、理，經(jīng)過調(diào)整使得兩種方法的計(jì)算結(jié)果相差甚小。但是由于考慮的不夠周全，沒有使兩個(gè)發(fā)電廠均衡的承擔(dān)負(fù)荷，使發(fā)電廠一的功率過小。</p><p>　　此次課程設(shè)計(jì)，讓我充分了解了解到精準(zhǔn)的潮流計(jì)算不僅可以節(jié)省能源并充分發(fā)揮能源的作用，還可以使電網(wǎng)處于穩(wěn)定且平穩(wěn)的狀態(tài)，這對電網(wǎng)的完全運(yùn)行起到關(guān)鍵性的作用。電力系統(tǒng)的安全運(yùn)行不是兒戲，一個(gè)很小的失誤就可能造成整個(gè)電網(wǎng)的崩潰，給國家和社會帶來極大的損失。我們應(yīng)該努力學(xué)好專業(yè)課

106、，打好專業(yè)基礎(chǔ)，在以后的工作崗位上有所成就。</p><p><b>　　參考文獻(xiàn)：</b></p><p>　　1、《電力系統(tǒng)分析基礎(chǔ)》，韋剛，中國電力出版社;</p><p>　　2、《電力系統(tǒng)穩(wěn)態(tài)分析》，陳珩，中國電力出版社;</p><p>　　3、《MATLAB實(shí)用指南》，蘇金明，王永利，電子工業(yè)出版社。<

107、;/p><p><b>　　附錄：</b></p><p>　　%本程序的功能是用牛頓——拉夫遜法進(jìn)行潮流計(jì)算</p><p>　　% B1矩陣：1、支路首端號；2、末端號；3、支路阻抗；4、支路對地電納</p><p>　　% 5、支路的變比；6、支路首端處于K側(cè)為1，1側(cè)為0</p><

108、;p>　　% B2矩陣：1、該節(jié)點(diǎn)發(fā)電機(jī)功率；2、該節(jié)點(diǎn)負(fù)荷功率；3、節(jié)點(diǎn)電壓初始值</p><p>　　% 4、PV節(jié)點(diǎn)電壓V的給定值；5、節(jié)點(diǎn)所接的無功補(bǔ)償設(shè)備的容量</p><p>　　% 6、節(jié)點(diǎn)分類標(biāo)號</p><p><b>　　clear;</b></p><p>　　

109、n=10;%input('請輸入節(jié)點(diǎn)數(shù)：n=');</p><p>　　nl=9;%input('請輸入支路數(shù)：nl=');</p><p>　　isb=1;%input('請輸入平衡母線節(jié)點(diǎn)號：isb=');</p><p>　　pr=0.00001;%input('請輸入誤差精度：pr=');<

110、;/p><p>　　B1=[1 3 8.5+20.1i 5.56e-4i 1 0 0;</p><p>　　1 5 4.25+10.05i 2.78e-4i 1 0 0;</p><p>　　1 7 7.65+18.09i 5.004e-4i 1 0 0;</p><p>　　2 9 8.5+20.1i 5.56e-4i 1 0 0;&l

111、t;/p><p>　　3 4 0.05455+1.65024i 0 231/38.5 1 1;</p><p>　　5 6 3.735e-3+0.10083i 0 231/11 1 1;</p><p>　　7 8 0.05455+1.65024i 0 231/38.5 1 1;</p><p>　　7 9 13.6+32.16i

112、2.224e-4i 1 0 0;</p><p>　　9 10 3.735e-3+0.10083i 0 231/11 1 1;</p><p>　　B2=[0 0 235 0 1;</p><p>　　200 0 235 0 3;</p><p>　　0 0 220 0 2;</p><p>　　

113、0 60+37.18466i 35 0 2;</p><p>　　0 0 220 0 2;</p><p>　　0 50+30.98722i 10 0 2;</p><p>　　0 0 220 0 2;</p><p>　　0 40+24.78977i 35 0 2;</p><p>　　0 0 220 0 2;<

114、;/p><p>　　0 60+37.18466i 10 0 2] %input('請輸入各節(jié)點(diǎn)參數(shù)形成的矩陣： B2=');</p><p>　　Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);</p><p>　　%－－－－－－－修改部分－－－

115、－－－－－－－－－</p><p><b>　　ym=0;</b></p><p>　　SB=100;UB=220;</p><p>　　%ym=input('您輸入的參數(shù)是標(biāo)么值？（若不是則輸入一個(gè)不為零的數(shù)值）');</p><p><b>　　if ym~=0</b></

116、p><p>　　%SB=input('請輸入功率基準(zhǔn)值:SB=');</p><p>　　%UB=input('請輸入電壓基準(zhǔn)值:UB=');</p><p>　　YB=SB./UB./UB;</p><p><b>　　BB1=B1;</b></p><p><

117、b>　　BB2=B2;</b></p><p>　　for i=1:nl</p><p>　　B1(i,3)=B1(i,3)*YB;</p><p>　　B1(i,4)=B1(i,4)./YB;</p><p><b>　　end</b></p><p>　　disp('

118、B1矩陣B1=');</p><p><b>　　disp(B1)</b></p><p><b>　　for i=1:n</b></p><p>　　B2(i,1)=B2(i,1)./SB;</p><p>　　B2(i,2)=B2(i,2)./SB;</p><p&g

119、t;　　B2(i,3)=B2(i,3)./UB;</p><p>　　B2(i,4)=B2(i,4)./UB;</p><p>　　B2(i,5)=B2(i,5)./SB;</p><p><b>　　end</b></p><p>　　disp('B2矩陣B2=');</p><p&