電力系統(tǒng)分析課程設(shè)計(jì)---復(fù)雜網(wǎng)絡(luò)牛頓—拉夫遜法潮流分析與計(jì)算的設(shè)計(jì)

1、<p>　　電力系統(tǒng)分析課程設(shè)計(jì)</p><p><b>　　報(bào)告書</b></p><p>　　院(部)別 </p><p>　　班 級(jí) </p><p>　　學(xué)

2、 號(hào) </p><p>　　姓 名 </p><p>　　指導(dǎo)教師 </p><p>　　時(shí) 間 2013.11.28-2013.1

3、1.29 </p><p>　　課 程 設(shè) 計(jì) 任 務(wù) 書</p><p>　　題 目 電力系統(tǒng)分析課程設(shè)計(jì) </p><p>　　學(xué) 院 </p><p>　　專 業(yè)

4、 </p><p>　　班 級(jí) </p><p>　　學(xué)生姓名 </p><p>　　學(xué) 號(hào) </p><

5、p>　　11 月 18 日至 11 月 29 日 共 2 周</p><p>　　指導(dǎo)教師(簽字) </p><p>　　院 長(簽字) </p><p>　　年 月 日</p><p>　　復(fù)雜網(wǎng)絡(luò)N-R法潮流分析與計(jì)算的設(shè)計(jì)</p>

6、<p>　　電力系統(tǒng)的潮流計(jì)算是電力系統(tǒng)分析課程基本計(jì)算的核心部分之一。它既有自身的獨(dú)立意義，又有電力系統(tǒng)規(guī)劃設(shè)計(jì)、運(yùn)行和研究的理論基礎(chǔ)，因此課程設(shè)計(jì)的重要性自不待言。</p><p><b>　　設(shè)計(jì)題目</b></p><p><b>　　1.系統(tǒng)圖的確定</b></p><p>　　選擇六節(jié)點(diǎn)、環(huán)網(wǎng)、兩電源

7、和多引出的電力系統(tǒng)，簡化電力系統(tǒng)圖如圖1所示，等值導(dǎo)納圖如圖2所示。運(yùn)用以直角坐標(biāo)表示的牛頓-拉夫遜計(jì)算如圖1所示系統(tǒng)中的潮流分布。計(jì)算精度要求各節(jié)點(diǎn)電壓的誤差或修正量不大于。 </p><p>　　圖1 電力系統(tǒng)圖</p><p>　　圖2 電力系統(tǒng)等值導(dǎo)

8、納圖</p><p>　　2.各節(jié)點(diǎn)的初值及阻抗參數(shù)</p><p>　　該系統(tǒng)中，節(jié)點(diǎn)①為平衡節(jié)點(diǎn)，保持=1.05+j0為定值，節(jié)點(diǎn)⑥為PV節(jié)點(diǎn)，其他四個(gè)節(jié)點(diǎn)都是PQ節(jié)點(diǎn)。給定的注入電壓標(biāo)幺值、線路阻抗標(biāo)幺值、輸出功率標(biāo)幺值如下表注釋。</p><p>　　表1 各節(jié)點(diǎn)電壓標(biāo)幺值參數(shù)</p><p&g

9、t;　　表2 線路、變壓器阻抗標(biāo)幺值 </p><p>　　表3 節(jié)點(diǎn)輸出功率</p><p>　　注：各PQ節(jié)點(diǎn)的電壓取1是為了方便計(jì)算和最后驗(yàn)證程序的正確性。</p><p><b>　　N-R法的求解過程</b></p><p>

10、;　　給定個(gè)節(jié)點(diǎn)電壓初始值</p><p>　　將以上電壓初始值代入下式（1）式，求出修正方程式常數(shù)項(xiàng)向量。</p><p>　　將電壓初始值代入下式（2）式，求出修正方程式中系數(shù)矩陣（雅可比矩陣）的元素（為2（n-1）階方陣）。</p><p>　　解修正方程式，求出修正量</p><p>　　方法：①高斯消元法，即將雅可比矩陣化簡為上三角

11、陣或下三角陣</p><p>　?、谇笱趴杀染仃嚨哪婢仃?lt;/p><p><b>　　修正各節(jié)點(diǎn)電壓</b></p><p>　　6、將再代入（1）式，求出</p><p>　　7、校驗(yàn)是否收斂，其收斂條件為式中，為向量中大分量的絕對(duì)值。這個(gè)收斂條件比較直觀，它可以直接顯示出最終結(jié)果的功率誤差</p>&l

12、t;p>　　8、如果收斂，就進(jìn)一步計(jì)算各段電力線路潮流和平衡節(jié)點(diǎn)功率，并打印出計(jì)算結(jié)果；如果不收斂，轉(zhuǎn)回3步進(jìn)行下一次迭代計(jì)算，直到收斂為止。</p><p><b>　　三、手算潮流計(jì)算</b></p><p>　　用圖1和上表的數(shù)據(jù)和等值網(wǎng)絡(luò)計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣</p><p><b>　　1. 節(jié)點(diǎn)導(dǎo)納矩陣</b>

13、</p><p>　　由圖1可知，該系統(tǒng)以串聯(lián)支路的導(dǎo)納標(biāo)幺值和對(duì)地并聯(lián)導(dǎo)納標(biāo)幺值得等值電路如圖2所示?？傻脠D2的等值節(jié)點(diǎn)導(dǎo)納矩陣。</p><p><b>　　對(duì)角線上的元素為：</b></p><p>　　非對(duì)角線上的元素為：</p><p><b>　　Y=</b></p>&l

14、t;p>　　2.計(jì)算各節(jié)點(diǎn)功率的修正方程的初始值（不平衡量）</p><p>　　設(shè)各節(jié)點(diǎn)電壓初始值為：</p><p>　　根據(jù)上述N-R法的求解過程，將各節(jié)點(diǎn)電壓的初始值代入（1）式和（2）式，進(jìn)行潮流計(jì)算，得節(jié)點(diǎn)功率和節(jié)點(diǎn)電壓：</p><p><b>　　雅可比矩陣：</b></p><p><b&g

15、t;　　J=</b></p><p><b>　　修正各節(jié)點(diǎn)電壓：</b></p><p>　　四、計(jì)算機(jī)算法潮流計(jì)算</p><p><b>　　導(dǎo)納矩陣=</b></p><p><b>　　time =1</b></p><p><

16、;b>　　JJ =</b></p><p><b>　　E =</b></p><p>　　1.1000 1.1712 1.0291 0.9668 1.1667 1.1000</p><p><b>　　F =</b></p><p>　　0

17、-0.1174 -0.1824 -0.3595 -0.0391 0.0358</p><p><b>　　U =</b></p><p>　　1.1000 1.1712-0.1174i 1.0291-0.1824i 0.9668-0.3595i 1.1667-0.0391i 1.1000+0.0358i</p><p&

18、gt;<b>　　dU =</b></p><p><b>　　0.1712</b></p><p><b>　　-0.1174</b></p><p><b>　　0.0291</b></p><p><b>　　-0.1824</b&g

19、t;</p><p><b>　　-0.0332</b></p><p><b>　　-0.3595</b></p><p><b>　　0.1667</b></p><p><b>　　-0.0391</b></p><p>&l

20、t;b>　　0</b></p><p><b>　　0.0358</b></p><p><b>　　PQ =</b></p><p>　　-2.0000 4.1863 -1.8000 0.1000 -1.6000 -0.3000 -3.7000 8.5726 5.0000 0</

21、p><p><b>　　time =2</b></p><p><b>　　JJ =</b></p><p><b>　　E =</b></p><p>　　1.1000 1.1208 0.9578 0.8007 1.1317 1.1008</p

22、><p><b>　　F =</b></p><p>　　0 -0.1321 -0.1605 -0.3271 -0.0785 -0.0054</p><p><b>　　U =</b></p><p>　　1.1000 1.1208-0.1321i 0.9578-0.1605

23、i 0.8007-0.3271i 1.1317-0.0785i 1.1008-0.0054i</p><p><b>　　dU =</b></p><p><b>　　-0.0504</b></p><p><b>　　-0.0147</b></p><p><b&

24、gt;　　-0.0713</b></p><p><b>　　0.0219</b></p><p><b>　　-0.1660</b></p><p><b>　　0.0324</b></p><p><b>　　-0.0350</b><

25、/p><p><b>　　-0.0394</b></p><p><b>　　0.0008</b></p><p><b>　　-0.0412</b></p><p><b>　　PQ =</b></p><p>　　-0.3576

26、-1.0319 0.0288 -0.3150 0.0666 -0.5211 0.1533 -1.8333 -0.3787 -0.0013</p><p><b>　　time =3</b></p><p><b>　　JJ =</b></p><p><b>　　E =</b><

27、/p><p>　　1.1000 1.1142 0.9438 0.7501 1.1271 1.0999</p><p><b>　　F =</b></p><p>　　0 -0.1351 -0.1560 -0.3221 -0.0889 -0.0169</p><p>&l

28、t;b>　　U =</b></p><p>　　1.1000 1.1142-0.1351i 0.9438-0.1560i 0.7501-0.3221i 1.1271-0.0889i 1.0999-0.0169i</p><p><b>　　dU =</b></p><p><b>　　-0.0065<

29、/b></p><p><b>　　-0.0030</b></p><p><b>　　-0.0140</b></p><p><b>　　0.0045</b></p><p><b>　　-0.0507</b></p><p&g

30、t;<b>　　0.0050</b></p><p><b>　　-0.0046</b></p><p><b>　　-0.0105</b></p><p><b>　　-0.0008</b></p><p><b>　　-0.0114</

31、b></p><p><b>　　PQ =</b></p><p>　　0.0000 -0.0509 -0.0099 -0.0207 -0.0154 -0.1040 0.0685 -0.0554 -0.0935 -0.0017</p><p><b>　　time =4</b></p>

32、<p><b>　　JJ =</b></p><p><b>　　E =</b></p><p>　　1.1000 1.1136 0.9424 0.7436 1.1267 1.0999</p><p><b>　　F =</b></p><

33、p>　　0 -0.1355 -0.1555 -0.3215 -0.0901 -0.0181</p><p><b>　　U =</b></p><p>　　1.1000 1.1136-0.1355i 0.9424-0.1555i 0.7436-0.3215i 1.1267-0.0901i 1.0999-0.0181i</p

34、><p><b>　　dU =</b></p><p><b>　　-0.0006</b></p><p><b>　　-0.0004</b></p><p><b>　　-0.0015</b></p><p><b>　　0

35、.0005</b></p><p><b>　　-0.0065</b></p><p><b>　　0.0006</b></p><p><b>　　-0.0005</b></p><p><b>　　-0.0012</b></p>

36、<p><b>　　-0.0001</b></p><p><b>　　-0.0012</b></p><p><b>　　PQ =</b></p><p>　　0.0002 -0.0004 -0.0006 -0.0003 -0.0024 -0.0111 0.0011 0.

37、0003 -0.0028 -0.0001</p><p><b>　　time =5</b></p><p><b>　　JJ =</b></p><p><b>　　E =</b></p><p>　　1.1000 1.1136 0.9423 0.743

38、5 1.1267 1.0999</p><p><b>　　F =</b></p><p>　　0 -0.1355 -0.1554 -0.3215 -0.0902 -0.0181</p><p><b>　　U =</b></p><p>　　1.1000

39、1.1136-0.1355i 0.9423-0.1554i 0.7435-0.3215i 1.1267-0.0902i 1.0999-0.0181i</p><p><b>　　dU =</b></p><p>　　1.0e-003 * </p><p><b>　　-0.0104</b></p&g

40、t;<p><b>　　-0.0063</b></p><p><b>　　-0.0255</b></p><p><b>　　0.0090</b></p><p><b>　　-0.1157</b></p><p><b>　　0

41、.0106</b></p><p><b>　　-0.0077</b></p><p><b>　　-0.0207</b></p><p><b>　　-0.0010</b></p><p><b>　　-0.0207</b></p>

42、;<p><b>　　PQ =</b></p><p>　　1.0e-003 *</p><p>　　0.0004 -0.0010 -0.0059 0.0003 -0.0439 -0.1922 0.0050 0.0002 -0.0289 -0.0015</p><p><b>　　time =6<

43、/b></p><p><b>　　JJ =</b></p><p><b>　　E =</b></p><p>　　1.1000 1.1136 0.9423 0.7435 1.1267 1.0999</p><p><b>　　F =</b>

44、;</p><p>　　0 -0.1355 -0.1554 -0.3215 -0.0902 -0.0181</p><p><b>　　U =</b></p><p>　　1.1000 1136-0.1355i 0.9423-0.1554i 0.7435-0.3215i 1.1267-0.0902i 1.0

45、999-0.0181i</p><p><b>　　dU =</b></p><p>　　1.0e-007 * </p><p><b>　　-0.0328</b></p><p><b>　　-0.0198</b></p><p><

46、;b>　　-0.0809</b></p><p><b>　　0.0286</b></p><p><b>　　-0.3674</b></p><p><b>　　0.0335</b></p><p><b>　　-0.0241</b>&

47、lt;/p><p><b>　　-0.0654</b></p><p><b>　　-0.0030</b></p><p><b>　　-0.0652</b></p><p><b>　　PQ =</b></p><p>　　1.0e-

48、007 *</p><p>　　0.0007 -0.0026 -0.0163 0.0023 -0.1421 -0.6090 0.0149 -0.0022 -0.0872 -0.0043</p><p><b>　　平衡節(jié)點(diǎn)功率</b></p><p><b>　　S1 =</b></p>

49、<p>　　4.7307 + 1.4446i</p><p>　　節(jié)點(diǎn)間線路上的流動(dòng)功率矩陣</p><p><b>　　S =</b></p><p>　　輸出電力線路功率總損耗</p><p><b>　　dataS =</b></p><p><b>

50、;　　個(gè)人心得</b></p><p>　　經(jīng)過兩周的課程設(shè)計(jì)讓我學(xué)到了很多以前沒有學(xué)到的東西，首先讓我明白了電力系統(tǒng)的穩(wěn)定運(yùn)行，必須經(jīng)過精密的設(shè)計(jì)和運(yùn)算。在進(jìn)行此次課程設(shè)計(jì)的過程中，加深了我對(duì)潮流計(jì)算的認(rèn)識(shí)，尤其是對(duì)牛頓拉夫遜潮流計(jì)算的求解思路有了比較透徹的理解。同時(shí)整個(gè)求解過程中用到求節(jié)點(diǎn)導(dǎo)納矩陣，求矩陣的逆陣等等，又讓我們遇到了很多問題，我們經(jīng)過查資料問同學(xué)、問老師，最后才把問題解決。使我對(duì)以前

51、學(xué)過的知識(shí)有了一次很好的復(fù)習(xí)，同時(shí)也看到了我知識(shí)的欠缺，讓我明白我應(yīng)該學(xué)的東西還有很多，更是給我更多的激勵(lì)讓我以后會(huì)更加努力地學(xué)習(xí)。此次課程設(shè)計(jì)有讓我學(xué)到了一本新的知識(shí)—Matlab的使用，起初我們經(jīng)過好大力氣才把程序設(shè)計(jì)出來，可是設(shè)計(jì)的題目里的數(shù)據(jù)又不收斂，我們反復(fù)演算了好多次，最后在老師的指導(dǎo)下改了幾個(gè)數(shù)據(jù)程序終于收斂了。</p><p>　　而且在此次課程設(shè)計(jì)中，我發(fā)現(xiàn)了自己的基礎(chǔ)只是有很多的不足。這些基礎(chǔ)

52、的缺乏給我的設(shè)計(jì)造成了不小的障礙，今后我會(huì)努力學(xué)習(xí)來彌補(bǔ)這些不足。在這個(gè)過程中，我明白了，只要用心去做，認(rèn)真去做，同學(xué)一起努力、團(tuán)結(jié)，就會(huì)克服所有困難，取得最后的勝利?？吹酵瑢W(xué)和自己設(shè)計(jì)的程序穩(wěn)定運(yùn)行，心中也有一種成功的喜悅。</p><p>　　通過此次的實(shí)習(xí)讓我發(fā)現(xiàn)手工的潮流計(jì)算和計(jì)算機(jī)潮流計(jì)算的優(yōu)缺，對(duì)于小系統(tǒng)來說手工計(jì)算很簡單，手工也很經(jīng)濟(jì)?？蓪?duì)于大系統(tǒng)來說手工計(jì)算是不可完成的任務(wù)，因此用計(jì)算機(jī)計(jì)算潮流是

53、必不可少的，只要編好程序計(jì)算機(jī)就會(huì)給你一個(gè)自己想要的滿意答案。而且方便、簡單、準(zhǔn)確。</p><p>　　此次課程設(shè)計(jì)我能成功設(shè)計(jì)完成離不開同學(xué)們的幫助、團(tuán)隊(duì)的合作以及老師細(xì)心的指導(dǎo)，在這里再次感謝老師的指導(dǎo)和同學(xué)的幫助。</p><p><b>　　附錄：</b></p><p>　　y11=1/0.03j;</p><p

54、>　　y12=-1/(1.05*0.03j);</p><p><b>　　y13=0;</b></p><p><b>　　y14=0;</b></p><p><b>　　y15=0;</b></p><p><b>　　y16=0;</b>&l

55、t;/p><p><b>　　y21=y12;</b></p><p>　　y22=0.5j+1/(0.06+0.025j)+1/(0.04+0.25j)+1/((1.05^2)*0.03j);</p><p>　　y23=-1/(0.06+0.025j);</p><p><b>　　y24=0;</b&g

56、t;</p><p>　　y25=-1/(0.04+0.25j);</p><p><b>　　y26=0;</b></p><p><b>　　y31=y13;</b></p><p><b>　　y32=y23;</b></p><p>　　y33=

57、0.5j+1/(0.06+0.025j)+1/(0.08+0.3j);</p><p>　　y34=-1/(0.08+0.3j);</p><p><b>　　y35=0; </b></p><p><b>　　y36=0;</b></p><p><b>　　y41=y14;</b

58、></p><p><b>　　y42=y24;</b></p><p><b>　　y43=y34;</b></p><p>　　y44=0.5j+1/(0.08+0.3j)+1/(0.1+0.35j);</p><p>　　y45=-1/(0.1+0.35j);</p>&l

59、t;p><b>　　y46=0;</b></p><p><b>　　y51=y15;</b></p><p><b>　　y52=y25;</b></p><p><b>　　y53=y35;</b></p><p><b>　　y54=

60、y45;</b></p><p>　　y55=0.5j+1/(0.1+0.35j)+1/(0.04+0.25j)+1/((1.05^2)*0.015j);</p><p>　　y56=-1/(1.05*0.015j);</p><p><b>　　y61=0;</b></p><p>　　y62=0;y63

61、=0;y64=0;</p><p><b>　　y65=y56;</b></p><p>　　y66=1/0.015j;%輸入導(dǎo)納參數(shù)</p><p>　　YB=[y11 y12 y13 y14 y15 y16; y21 y22 y23 y24 y25 y26;y31 y32 y33 y34 y35 y36; y41 y42 y43 y44 y

62、45 y46;y51 y52 y53 y54 y55 y56;y61 y62 y63 y64 y65 y66]%形成導(dǎo)納矩陣</p><p>　　%*********牛頓—拉夫遜法潮流計(jì)算************</p><p>　　E(2)=1.00;E(3)=1.00;E(4)=1.00;E(5)=1.00;E(6)=1.1;</p><p>　　F(2)=0;

63、F(3)=0;F(4)=0;F(5)=0; F(6)=0; %設(shè)定電壓初始值</p><p>　　G=real(YB);B=imag(YB); %提取導(dǎo)納實(shí)部,虛部分別為G,B</p><p>　　S(2)=-2-1j;S(3)=-1.8-0.4j;S(4)=-1.6-0.8j;S(5)=-3.7-1.3j; S(6)=5;%設(shè)定注入功率初始值</p><p>　　

64、P=real(S);Q=imag(S); %提取注入功率實(shí)部,虛部分別為P,Q</p><p>　　k=0;precision=0.1; %設(shè)定參數(shù)k為迭代次數(shù)初值為0及精度大于0.00001</p><p>　　N1=5;%設(shè)定除平衡節(jié)點(diǎn)外節(jié)點(diǎn)數(shù)之和</p><p>　　while precision > 0.00001</p><p&

65、gt;　　E(1)=1.1;F(1)=0;%輸入平衡節(jié)點(diǎn)電壓實(shí)部,虛部e,f</p><p>　　for m=2:N1+1</p><p>　　for n=1:N1+1</p><p>　　Pt(n)=(E(m)*(G(m,n)*E(n)-B(m,n)*F(n))+F(m)*(G(m,n)*F(n)+B(m,n)*E(n)));</p><p&g

66、t;　　Qt(n)=(F(m)*(G(m,n)*E(n)-B(m,n)*F(n))-E(m)*(G(m,n)*F(n)+B(m,n)*E(n)));%以設(shè)定值計(jì)算的注入功率</p><p><b>　　end</b></p><p>　　dP(m)=P(m)-sum(Pt);</p><p>　　dQ(m)=Q(m)-sum(Qt); %計(jì)算功

67、率誤差</p><p><b>　　end</b></p><p>　　for m=2:N1+1</p><p>　　for n=1:N1+1 </p><p>　　Ai(n)=G(m,n)*E(n)-B(m,n)*F(n); </p><p>　　Bi(n)=G(m,n)*F(n)+B(m,n)

68、*E(n);</p><p><b>　　end</b></p><p>　　N(m,m)=sum(Ai)+G(m,m)*E(m)+B(m,m)*F(m);%計(jì)算dP/de</p><p>　　H(m,m)=sum(Bi)-B(m,m)*E(m)+G(m,m)*F(m);%計(jì)算dP/df</p><p>　　L(m,m)

69、= -sum(Bi)-B(m,m)*E(m)+G(m,m)*F(m);%計(jì)算dQ/de</p><p>　　J(m,m)= sum(Ai)-G(m,m)*E(m)-B(m,m)*F(m);%計(jì)算dQ/df</p><p><b>　　end</b></p><p>　　for m=2:N1+1</p><p>　　JJ(

70、2*m-3,2*m-3)=N (m,m);</p><p>　　JJ(2*m-3,2*m-2)=H(m,m);</p><p>　　JJ(2*m-2,2*m-3)=L(m,m);</p><p>　　JJ(2*m-2,2*m-2)=J(m,m);</p><p>　　end %給雅可比矩陣的對(duì)角元素賦值</p><p>

71、;　　for m=2:N1+1</p><p>　　for n=2:N1+1</p><p><b>　　if m==n</b></p><p><b>　　else</b></p><p>　　N(m,n)=G(m,n)*E(m)+B(m,n)*F(m); %計(jì)算dP/de</p>

72、<p>　　H(m,n)=-B(m,n)*E(m)+G(m,n)*F(m); %計(jì)算dP/df</p><p>　　L(m,n)=-B(m,n)*E(m)+G(m,n)*F(m); %計(jì)算dQ/de</p><p>　　J(m,n)=-B(m,n)*F(m)-G(m,n)*E(m);%計(jì)算dQ/df</p><p>　　%************給雅可比矩

73、陣的非對(duì)角元素賦值**********</p><p>　　JJ(2*m-3,2*n-3)=N(m,n);</p><p>　　JJ(2*m-3,2*n-2)=H(m,n);</p><p>　　JJ(2*m-2,2*n-3)=L(m,n);</p><p>　　JJ(2*m-2,2*n-2)=J(m,n);</p><p

74、><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　for i=1:8</b></p><p>　　JJ(10,i)=0;%修正雅可比矩陣第十行1,到8元素

75、為0</p><p><b>　　end</b></p><p>　　JJ(10,9)=2*E(6); %修正雅可比矩陣第十行9元素為2e</p><p>　　JJ(10,10)=2*F(6);%修正雅可比矩陣第十行10元素為2f</p><p>　　for m=1:N1</p><p>　　PQ

76、(2*m-1)=dP(m+1); PQ(2*m)=dQ(m+1);</p><p>　　end%寫PQ矩陣，奇數(shù)行為dP，偶數(shù)行為dQ</p><p>　　PQ(10)=1.1^2-(E(6)^2+F(6)^2);%修正PQ矩陣第十個(gè)元素為du^2</p><p>　　dU=inv(JJ)*PQ';%求dU為JJ的逆矩陣乘以PQ'</p>

77、;<p>　　precision=max(abs(PQ)); %取精度為PQ的絕對(duì)值的最大值</p><p>　　for n=1:N1</p><p>　　E(n+1)=E(n+1)+dU(2*n-1);</p><p>　　F(n+1)=F(n+1)+dU(2*n); %修正各節(jié)點(diǎn)電壓</p><p><b>　　

78、end</b></p><p>　　for n=1:N1+1</p><p>　　U(n)=E(n)+(F(n))*j;%把電壓寫成e+jf的形式</p><p><b>　　end</b></p><p>　　k=k+1;%迭代次數(shù)加1</p><p>　　time=k,JJ, E

79、,F, U , dU,PQ %輸出k,JJ, E,F,U,du,ds</p><p><b>　　if k>9;</b></p><p><b>　　break</b></p><p><b>　　end</b></p><p><b>　　end</b&

80、gt;</p><p>　　%************功率及功率損耗的計(jì)算***********</p><p>　　for m=1:N1+1</p><p>　　I(m)=YB(1,m)*U(m);%計(jì)算平衡節(jié)點(diǎn)電流</p><p><b>　　end</b></p><p>　　S1=U(1)

81、*sum(conj(I)); %計(jì)算平衡節(jié)點(diǎn)功率</p><p>　　y10=0.05/(1.05*0.03j);</p><p>　　y20=-0.05/((1.05^2)*0.03j)+0.5j;</p><p><b>　　y30=0.5j;</b></p><p><b>　　y40=0.5j;<

82、/b></p><p>　　y50=0.5j-0.05/((1.05^2)*0.015j);</p><p>　　y60=0.05/(1.05*0.015j);</p><p>　　Y=[y10 y20 y30 y40 y50 y60]; %輸入各節(jié)點(diǎn)對(duì)地導(dǎo)納矩陣</p><p><b>　　for m=1:6</b&g

83、t;</p><p><b>　　for n=1:6</b></p><p>　　if m==n%當(dāng)m=n時(shí)</p><p>　　S(m,n)=0; %流動(dòng)節(jié)點(diǎn)功率為0</p><p><b>　　else%否則</b></p><p>　　S(m,n)=U(m)*( con

84、j(U(m)*Y(m))+( conj(U(m))- conj(U(n)))* conj(YB(m,n))); %節(jié)點(diǎn)間線路上的流動(dòng)功率矩陣</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><