1、<p>　　《電力系統(tǒng)分析基礎(chǔ)》課程設(shè)計(jì)說明書</p><p>　　題目：電力系統(tǒng)潮流(手工)計(jì)算</p><p>　　學(xué)生姓名 </p><p>　　學(xué) 號 </p><p>　　專 業(yè) 電氣工程及其自動化 </p>&l

2、t;p>　　班 級 </p><p>　　指導(dǎo)教師 </p><p>　　完成時(shí)間 2012-06-18 </p><p><b>　　目 錄</b></p><p>　　電力系統(tǒng)課程設(shè)計(jì)任務(wù)書1</p&

3、gt;<p>　　一、課程設(shè)計(jì)目的1</p><p>　　二、課程設(shè)計(jì)要求1</p><p>　　三、課程設(shè)計(jì)進(jìn)度2</p><p>　　四、設(shè)計(jì)成果要求2</p><p>　　電力系統(tǒng)課程設(shè)計(jì)說明書3</p><p>　　第一章 潮流計(jì)算簡介3</p><p>　　1

4、.1 潮流計(jì)算的概述3</p><p>　　1.2 潮流計(jì)算的發(fā)展史3</p><p>　　第二章 手工潮流計(jì)算4</p><p>　　2.1 手工潮流計(jì)算過程4</p><p>　　2.2 手工潮流計(jì)算心得4</p><p>　　第三章 潮流計(jì)算的編程4</p><p>　　3.1

5、 潮流計(jì)算流程圖4</p><p>　　3.2 牛頓-拉夫遜法程序6</p><p>　　3.3 牛頓-拉夫遜法程序結(jié)果10</p><p><b>　　思考題12</b></p><p><b>　　課程設(shè)計(jì)總結(jié)12</b></p><p><b>　　

6、參考文獻(xiàn)13</b></p><p><b>　　附錄14</b></p><p><b>　　摘 要</b></p><p>　　潮流計(jì)算是根據(jù)給定的電網(wǎng)結(jié)構(gòu)、參數(shù)和發(fā)電機(jī)、負(fù)荷等元件的運(yùn)行條件，確定電力系統(tǒng)各部分穩(wěn)態(tài)運(yùn)行狀態(tài)參數(shù)計(jì)算。通常給定的運(yùn)行條件有系統(tǒng)中各電源和負(fù)荷點(diǎn)的功率、樞紐點(diǎn)電壓、平衡點(diǎn)的

7、電壓和相位角。</p><p>　　待求的運(yùn)行狀態(tài)參量包括電網(wǎng)各母線節(jié)點(diǎn)電壓幅值和相角，以及各支路的功率分布、網(wǎng)絡(luò)的功率損耗等。</p><p>　　潮流計(jì)算是電力系統(tǒng)非常重要的分析計(jì)算，用于系統(tǒng)規(guī)劃和運(yùn)行中提出的各種問題。</p><p>　　對規(guī)劃中的電力系統(tǒng)，通過潮流計(jì)算可以預(yù)知各種負(fù)荷變化和網(wǎng)絡(luò)結(jié)構(gòu)的改變會不會危機(jī)系統(tǒng)的安全，系統(tǒng)中所有母線的電壓是否在允許的

8、范圍之內(nèi)，系統(tǒng)中各元件是否會出現(xiàn)過負(fù)荷，以及可能出現(xiàn)過負(fù)荷時(shí)應(yīng)事先采取那些預(yù)防措施等。</p><p>　　關(guān)鍵詞：潮流計(jì)算 牛頓法 P-Q分解法 電力系統(tǒng)分析</p><p>　　電力系統(tǒng)課程設(shè)計(jì)任務(wù)書</p><p><b>　　一、課程設(shè)計(jì)目的</b></p><p>　　掌握電力系統(tǒng)潮流計(jì)算的基本原理，</

9、p><p>　　掌握并能熟練運(yùn)用一門計(jì)算機(jī)語言對潮流計(jì)算進(jìn)行計(jì)算機(jī)編程計(jì)算。</p><p>　　采用計(jì)算機(jī)語言對潮流計(jì)算進(jìn)行計(jì)算機(jī)編程</p><p><b>　　二、課程設(shè)計(jì)要求</b></p><p><b>　　題目：</b></p><p>　　簡單系統(tǒng)如下圖所示，在該

10、電力系統(tǒng)中，1、2為PQ節(jié)點(diǎn)，3為P節(jié)點(diǎn)，4為平衡節(jié)點(diǎn)；標(biāo)幺值表示的各元件參數(shù)和節(jié)點(diǎn)數(shù)據(jù)如表1所示，給定電壓初值見表2。收斂系數(shù)：。</p><p>　　試求：采用以下兩種方法的該電力系統(tǒng)潮流分布，并將計(jì)算結(jié)果對比。</p><p>　　1.極坐標(biāo)下的手工計(jì)算(手寫,采用A4紙)，進(jìn)行2次迭代，計(jì)算出結(jié)果；</p><p>　　2.極坐標(biāo)下的牛頓-拉夫遜法或P-Q分

11、解法求解(程序源代碼、說明、流程圖、運(yùn)行結(jié)果，電子版)。</p><p>　　表1 標(biāo)幺值表示的網(wǎng)絡(luò)元件參數(shù)和各節(jié)點(diǎn)數(shù)據(jù)</p><p>　　表2 給定的各節(jié)點(diǎn)電壓初始值(標(biāo)幺值)參數(shù)</p><p><b>　　思考題</b></p><p>　　1）潮流計(jì)算的方法有哪些？各有何特點(diǎn)？</p><

12、;p>　　2）設(shè)計(jì)中遇到的問題和解決的辦法。</p><p><b>　　三、課程設(shè)計(jì)進(jìn)度</b></p><p>　　四、 設(shè)計(jì)成果要求 </p><p><b>　　1. 手算潮流；</b></p><p>　　2. 計(jì)算機(jī)計(jì)算潮流程序流程圖；</p><p>　　

13、3. 完整的潮流程序，并有程序注釋。包括變量的定義，子程序的調(diào)用等，計(jì)算結(jié)果；</p><p><b>　　4. 回答思考題</b></p><p>　　電力系統(tǒng)課程設(shè)計(jì)說明書</p><p>　　第一章 潮流計(jì)算簡介</p><p>　　1.1潮流計(jì)算的概述</p><p>　　為了使電力系統(tǒng)能

14、夠安全、優(yōu)質(zhì)、經(jīng)濟(jì)的運(yùn)行，在電力系統(tǒng)的設(shè)計(jì)、運(yùn)行以及研究工作中，需要針對電力系統(tǒng)的不同運(yùn)行狀態(tài)，做多種多樣的計(jì)算。指在給定電力系統(tǒng)網(wǎng)絡(luò)拓?fù)?、元件參?shù)和發(fā)電、負(fù)荷參量條件下，計(jì)算有功功率、無功功率及電壓在電力網(wǎng)中的分布。潮流計(jì)算是根據(jù)給定的電網(wǎng)結(jié)構(gòu)、參數(shù)和發(fā)電機(jī)、負(fù)荷等元件的運(yùn)行條件，確定電力系統(tǒng)各部分穩(wěn)態(tài)運(yùn)行狀態(tài)參數(shù)的計(jì)算。通常給定的運(yùn)行條件有系統(tǒng)中各電源和負(fù)荷點(diǎn)的功率、樞紐點(diǎn)電壓、平衡點(diǎn)的電壓和相位角。待求的運(yùn)行狀態(tài)參量包括電網(wǎng)各母線

15、節(jié)點(diǎn)的電壓幅值和相角，以及各支路的功率分布、網(wǎng)絡(luò)的功率損耗等。</p><p>　　1.2 潮流計(jì)算的發(fā)展史</p><p>　　近二十多年來，潮流計(jì)算的研究仍然非常活躍，但是大多數(shù)研究都是圍繞改進(jìn)牛頓法和P-Q分解法進(jìn)行的。此外，隨著人工智能理論的發(fā)展，遺傳算法、人工神經(jīng)網(wǎng)絡(luò)、模糊算法頁逐漸被引入潮流計(jì)算。但是，到目前為止這些新的模型和算法還不能取代牛頓法和P-Q分解法的地位。由于電力系

16、統(tǒng)規(guī)模的不斷擴(kuò)大，對計(jì)算機(jī)速度的要求不斷提高，計(jì)算機(jī)的并行計(jì)算技術(shù)也將在潮流計(jì)算中得到廣泛應(yīng)用，成為重要的研究領(lǐng)域。</p><p>　　在電力市場的改革不斷深入，“廠網(wǎng)分開，競價(jià)上網(wǎng)”的今天，傳統(tǒng)的潮流計(jì)算又遇到了一些新的挑戰(zhàn)。例如，在傳統(tǒng)潮流計(jì)算中，在每一步潮流迭代中所有節(jié)點(diǎn)的有功不平衡和無功不平衡都是平衡節(jié)點(diǎn)功率來平衡的，在潮流計(jì)算收斂之前，平衡節(jié)點(diǎn)的功率是無法確定的，相當(dāng)于被選為平衡節(jié)點(diǎn)的發(fā)電廠為全網(wǎng)的功

17、率平衡和潮流收斂做出了“犧牲”。這種做法在提倡“公平、公開、公正”的電力市場中受到質(zhì)疑。為體現(xiàn)市場公平的原則，在潮流計(jì)算的每一步中，讓所有調(diào)節(jié)能力的發(fā)電廠按照其裝機(jī)容量的大小來共同承擔(dān)系統(tǒng)總的不平衡功率，這樣不但能使各發(fā)電廠公平競爭，還能加快潮流計(jì)算的收斂速度。這就是所謂的動態(tài)潮流計(jì)算。另外，還有考慮負(fù)荷不確定性的區(qū)間潮流計(jì)算、考慮節(jié)點(diǎn)注入功率和負(fù)荷隨機(jī)特征的概率潮流計(jì)算、考慮負(fù)序非線性分布的三相潮流計(jì)算和考慮各次諧波在電力網(wǎng)中滲透情況

18、的諧波潮流計(jì)算等。</p><p>　　第二章 手工潮流計(jì)算</p><p>　　2.1 手工潮流計(jì)算過程</p><p><b>　　(見附錄)</b></p><p>　　2.2 手工潮流計(jì)算心得</p><p>　　通過此次的手工計(jì)算，我發(fā)現(xiàn)了很多細(xì)節(jié)問題，對公式的運(yùn)用更加熟練，同時(shí)對相關(guān)的

19、知識進(jìn)行了一次全面的梳理，比如運(yùn)用增廣矩陣求解多元方程組。由于數(shù)量太大，導(dǎo)致數(shù)據(jù)運(yùn)算過程中出現(xiàn)了不少的計(jì)算錯(cuò)誤，檢驗(yàn)數(shù)據(jù)正確性花費(fèi)了很多時(shí)間。最終通過一天的數(shù)據(jù)核算，發(fā)現(xiàn)誤差很小。</p><p>　　第三章 潮流計(jì)算的編程設(shè)計(jì)</p><p>　　3.1潮流計(jì)算流程圖</p><p><b>　　計(jì)算機(jī)計(jì)算</b></p>&l

20、t;p>　　（1）節(jié)點(diǎn)導(dǎo)納矩陣Y</p><p>　　根據(jù)題目提供的各節(jié)點(diǎn)的參數(shù)，求的節(jié)點(diǎn)導(dǎo)納矩陣</p><p>　　Yii=yio+Σyij</p><p><b>　　Yik=-yik</b></p><p>　　（2）簡化雅可比矩陣B’和B’’</p><p>　　通過上一步的導(dǎo)納矩

21、陣，形成有功迭代和無功迭代的簡化雅可比矩陣B’和B’’對雅可比矩陣進(jìn)行三角分解，形成因子表，為后面進(jìn)行修正方程計(jì)算做好準(zhǔn)備。</p><p><b>　?。?）修正和迭代</b></p><p>　　第一步，給定PQ節(jié)點(diǎn)初值和各節(jié)點(diǎn)電壓相角初值。</p><p>　　第二步，作第一次有功迭代，按公式計(jì)算節(jié)點(diǎn)有功功率不平衡量。</p>

22、<p>　　第三步，做第一次無功迭代，按公式計(jì)算無功功率不平衡量，計(jì)算時(shí)電壓相角最新的修正值。解修正方程式，可得各節(jié)點(diǎn)電壓幅值的修正量。</p><p>　　第四步，第一輪有功迭代和無功迭代便做完了。</p><p>　　第五步，按公式計(jì)算平衡結(jié)點(diǎn)功率。直到節(jié)點(diǎn)不平衡功率下降到10-5以下，迭代便可結(jié)束。</p><p>　　本次課程設(shè)計(jì)采用極坐標(biāo)下的

23、牛頓-拉夫遜計(jì)算網(wǎng)絡(luò)的潮流計(jì)算。其牛頓-拉夫遜潮流計(jì)算程序框圖如下所示。</p><p>　　3.2牛頓-拉夫遜法程序</p><p>　　>> %電力系統(tǒng)極坐標(biāo)下的牛頓-拉夫遜法潮流計(jì)算</p><p>　　disp('電力系統(tǒng)極坐標(biāo)下的牛頓-拉夫遜法潮流計(jì)算:');</p><p><b>　　cle

24、ar</b></p><p>　　n=input('請輸入結(jié)點(diǎn)數(shù)：n=');</p><p>　　n1=input('請輸入PV結(jié)點(diǎn)數(shù)：n1=');</p><p>　　n2=input('請輸入PQ結(jié)點(diǎn)數(shù)：n2=');</p><p>　　isb=input('請輸入平衡結(jié)

25、點(diǎn)：isb=');</p><p>　　pr=input('請輸入精確度：pr=');</p><p>　　K=input('請輸入變比矩陣看:K=');</p><p>　　C=input('請輸入支路阻抗矩陣：C=');</p><p>　　y=input('請輸入支路導(dǎo)納矩

26、陣：y=');</p><p>　　U=input('請輸入結(jié)點(diǎn)電壓矩陣：U=');</p><p>　　S=input('請輸入各結(jié)點(diǎn)的功率：S=');</p><p>　　Z=zeros(1,n);N=zeros(n1+n2,n2);L=zeros(n2,n2);QT1=zeros(1,n1+n2);</p>

27、<p><b>　　for m=1:n</b></p><p><b>　　for R=1:n</b></p><p>　　C(m,m)=C(m,m)+y(m,R);</p><p>　　if K(m,R)~=0</p><p>　　C(m,m)=C(m,m)+1/(C(m,R) /(

28、K(m,R) * (K(m,R)-1))) ;</p><p>　　C(R,R)=C(R,R)+1/(C(m,R)/(1-K(m,R)));</p><p>　　C(m,R)=C(m,R)/K(m,R);</p><p>　　C(R,m)=C(m,R);</p><p><b>　　end</b></p>

29、<p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　for m=1:n</b></p><p><b>　　for R=1:n</b></p><p><b>　　if m~=R&

30、lt;/b></p><p>　　Z(m)=Z(m)+1/C(m,R);</p><p><b>　　end </b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　f

31、or m=1:n</b></p><p><b>　　for R=1:n</b></p><p><b>　　if m==R</b></p><p>　　Y(m,m)=C(m,m)+Z(m);</p><p><b>　　else</b></p>&l

32、t;p>　　Y(m,R)=-1/C(m,R);</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　disp('結(jié)點(diǎn)導(dǎo)納矩陣:');</p>

33、<p><b>　　disp(Y);</b></p><p>　　disp('迭代中的雅克比矩陣:');</p><p>　　G=real(Y);</p><p>　　B=imag(Y);</p><p>　　O=angle(U);</p><p>　　U1=abs(U)

34、;</p><p><b>　　k=0;</b></p><p><b>　　PR=1;</b></p><p>　　P=real(S);</p><p>　　Q=imag(S);</p><p>　　while PR>pr</p><p>　　

35、for m=1:n2</p><p>　　UD(m)=U1(m);</p><p><b>　　end</b></p><p>　　for m=1:n1+n2</p><p><b>　　for R=1:n</b></p><p>　　PT(R)=U1(m)*U1(R)*(G

36、(m,R)*cos(O(m)-O(R))+B(m,R)*sin(O(m)-O(R)));</p><p><b>　　end</b></p><p>　　PT1(m)=sum(PT);</p><p>　　PP(m)=P(m)-PT1(m);</p><p>　　PP1(k+1,m)=PP(m);</p>

37、<p><b>　　end</b></p><p>　　for m=1:n2</p><p><b>　　for R=1:n</b></p><p>　　QT(R)=U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R))-B(m,R)*cos(O(m)-O(R)));</p><

38、p>　　end </p><p>　　QT1(m)=sum(QT); </p><p>　　QQ(m)=Q(m)-QT1(m); </p><p>　　QQ1(k+1,m)=QQ(m);</p><p><b>　　end</b></p><p>　　

39、PR1=max(abs(PP));</p><p>　　PR2=max(abs(QQ));</p><p>　　PR=max(PR1,PR2);</p><p>　　for m=1:n1+n2</p><p>　　for R=1:n1+n2</p><p><b>　　if m==R</b><

40、;/p><p>　　H(m,m)=U1(m)^2*B(m,m)+QT1(m);</p><p><b>　　else</b></p><p>　　H(m,R)=-U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R))-B(m,R)*cos(O(m)-O(R)));</p><p><b>　　end&

41、lt;/b></p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　for m=1:n1+n2</p><p>　　for R=1:n2</p><p><b>　　if m==R</b><

42、;/p><p>　　N(m,m)=-U1(m)^2*G(m,m)-PT1(m);</p><p><b>　　else</b></p><p>　　N(m,R)=-U1(m)*U1(R)*(G(m,R)*cos(O(m)-O(R))+B(m,R)*sin(O(m)-O(R)));</p><p><b>　　end

43、</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　for m=1:n2</p><p>　　for R=1:n1+n2</p><p><b>　　if m==R</b>&l

44、t;/p><p>　　J(m,m)=U1(m)^2*G(m,m)-PT1(m);</p><p><b>　　else</b></p><p>　　J(m,R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)-O(R))+B(m,R)*sin(O(m)-O(R)));</p><p><b>　　end&

45、lt;/b></p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　for m=1:n2</p><p>　　for R=1:n2</p><p><b>　　if m==R</b></p

46、><p>　　L(m,m)=U1(m)^2*B(m,m)-QT1(m);</p><p>　　else </p><p>　　L(m,R)=-U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R))-B(m,R)*cos(O(m)-O(R)));</p><p><b>　　end</b>

47、;</p><p><b>　　end</b></p><p><b>　　end </b></p><p>　　JJ=[H N;J L];</p><p><b>　　disp(JJ);</b></p><p>　　PQ=[PP';QQ&#

48、39;];</p><p>　　DA=-inv(JJ)*PQ; </p><p><b>　　DA1=DA';</b></p><p>　　for m=1:n1+n2</p><p>　　OO(m)=DA1(m);</p><p><b>　　end</b><

49、/p><p>　　for m=n:n1+n2+n2</p><p>　　UU1(m-n1-n2)=DA1(m);</p><p><b>　　end</b></p><p>　　UD2=diag(UD);</p><p>　　UU=UU1*UD2;</p><p>　　for

50、m=1:n1+n2</p><p>　　O(m)=O(m)+OO(m);</p><p><b>　　end</b></p><p>　　for m=1:n2</p><p>　　U1(m)=U1(m)+UU(m);</p><p><b>　　end</b></p&g

51、t;<p>　　for m=1:n1+n2</p><p>　　o(k+1,m)=180/pi*O(m);</p><p><b>　　end</b></p><p>　　for m=1:n2</p><p>　　u(k+1,m)=U1(m);</p><p><b>　　

52、end</b></p><p><b>　　k=k+1;</b></p><p><b>　　end</b></p><p><b>　　for m=1:n</b></p><p>　　b(m)=U1(m)*cos(O(m));</p><p&g

53、t;　　c(m)=U1(m)*sin(O(m));</p><p><b>　　end</b></p><p><b>　　U=b+i*c;</b></p><p><b>　　for R=1:n</b></p><p>　　PH1(R)=U(isb)*conj(Y(isb,R

54、))*conj(U(R));</p><p><b>　　end</b></p><p>　　PH=sum(PH1);</p><p><b>　　for m=1:n</b></p><p><b>　　for R=1:n</b></p><p><

55、;b>　　if m~=R</b></p><p>　　C1(m,R)=1/C(m,R);</p><p><b>　　else</b></p><p>　　C1(m,m)=C(m,m);</p><p><b>　　end</b></p><p><b

56、>　　end</b></p><p><b>　　end</b></p><p><b>　　for m=1:n</b></p><p><b>　　for R=1:n</b></p><p>　　if (C(m,R)~=inf)&(m~=R)<

57、/p><p>　　SS(m,R)=U1(m)^2*conj(C1(m,m))+U(m)*(conj(U(m))-conj(U(R)))*conj(C1(m,R));</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end<

58、;/b></p><p>　　disp('迭代中的△P：');disp(PP1);</p><p>　　disp('迭代中的△Q：');disp(QQ1);</p><p>　　disp('迭代中相角:');disp(o);</p><p>　　disp('迭代中電壓的模:'

59、;);disp(u);</p><p>　　disp('平衡結(jié)點(diǎn)的功率:');disp(PH);</p><p>　　disp('全部線路功率分布:');disp(SS);</p><p>　　3.3牛頓-拉夫遜法程序結(jié)果</p><p>　　表1 標(biāo)幺值表示的網(wǎng)絡(luò)元件參數(shù)和各節(jié)點(diǎn)數(shù)據(jù)</p>

60、<p>　　表2 給定的各節(jié)點(diǎn)電壓初始值(標(biāo)幺值)參數(shù)</p><p>　　電力系統(tǒng)極坐標(biāo)下的牛頓-拉夫遜法潮流計(jì)算:</p><p>　　請輸入結(jié)點(diǎn)數(shù)：n=4</p><p>　　請輸入PV結(jié)點(diǎn)數(shù)：n1=1</p><p>　　請輸入PQ結(jié)點(diǎn)數(shù)：n2=2</p><p>　　請輸入平衡結(jié)點(diǎn)：isb=4&l

61、t;/p><p>　　請輸入精確度：pr=0.00001</p><p>　　請輸入變比矩陣看:K= [0,0,0,0;0,0,0.9625,0;0,0,0,0;0,0,0,0]</p><p>　　請輸入支路阻抗矩陣：C=[0 0.04+0.12i 0.02+0.06i inf;0.04+0.12i 0 0.06+0.14i 0.00+0.05i;0.02+0.06i

62、 0.06+0.14i 0 0.04+0.1i;inf 0.00+0.05i 0.04+0.1i 0]</p><p>　　請輸入支路導(dǎo)納矩陣：y=[0 0.03i 0.0 0;0.03i 0 0 0;0.03i 0 0 0;0 0 0 0]</p><p>　　請輸入結(jié)點(diǎn)電壓矩陣：U=[1+0i 1+0i 1.00+0i 1.05+0i]</p><p>　　請輸

63、入各結(jié)點(diǎn)的功率：S=[-0.3-0.2i -0.2-0.1i 0.3 0]</p><p><b>　　結(jié)點(diǎn)導(dǎo)納矩陣:</b></p><p>　　7.5000 -22.4700i -2.5000 + 7.5000i -5.0000 +15.0000i 0 </p><p>　　-2.5000 + 7.500

64、0i 4.8959 -33.0604i -2.4892 + 5.8082i 0 +20.0000i</p><p>　　-5.0000 +15.0000i -2.4892 + 5.8082i 11.0345 -29.6252i -3.4483 + 8.6207i</p><p>　　0 0 +20.0000i -3.4483 + 8.62

65、07i 3.4483 -28.6207i</p><p>　　迭代中的雅克比矩陣:</p><p>　　-22.5000 7.5000 15.0000 -7.5000 2.5000</p><p>　　7.5000 -34.3082 5.8082 2.5000 -4.8025</p><p>　　15

66、.0000 5.8082 -29.6252 5.0000 2.4892</p><p>　　7.5000 -2.5000 -5.0000 -22.4400 7.5000</p><p>　　-2.5000 4.9892 -2.4892 7.5000 -31.8126</p><p>　　-22.7001 7

67、.7721 14.9279 -7.2118 2.5644</p><p>　　7.7564 -35.4387 5.9582 2.6117 -5.0365</p><p>　　15.0893 6.0569 -29.6252 4.7609 2.4595</p><p>　　7.8093 -2.5644 -5.2

68、449 -22.3033 7.7721</p><p>　　-2.6117 5.4427 -2.6898 7.7564 -35.3241</p><p>　　-22.6788 7.7573 14.9215 -7.2030 2.5639</p><p>　　7.7442 -35.3945 5.9528 2.6

69、033 -5.0267</p><p>　　15.0806 6.0477 -29.6252 4.7617 2.4608</p><p>　　7.8029 -2.5639 -5.2390 -22.2788 7.7573</p><p>　　-2.6033 5.4267 -2.6822 7.7442 -35.19

70、46</p><p>　　-22.6788 7.7573 14.9214 -7.2029 2.5638</p><p>　　7.7441 -35.3943 5.9527 2.6033 -5.0267</p><p>　　15.0806 6.0477 -29.6252 4.7616 2.4607</p&g

71、t;<p>　　7.8029 -2.5638 -5.2391 -22.2788 7.7573</p><p>　　-2.6033 5.4267 -2.6823 7.7441 -35.1943</p><p>　　-22.6788 7.7573 14.9214 -7.2029 2.5638</p><p&

72、gt;　　7.7442 -35.3943 5.9527 2.6033 -5.0267</p><p>　　15.0806 6.0477 -29.6252 4.7616 2.4607</p><p>　　7.8029 -2.5638 -5.2391 -22.2788 7.7573</p><p>　　-2.6033

73、 5.4267 -2.6823 7.7442 -35.1943</p><p><b>　　迭代中的△P：</b></p><p>　　-0.3000 -0.1067 0.3754</p><p>　　-0.0012 0.0031 -0.0084</p><p>　　-0.0000

74、 0.0000 0.0004</p><p>　　0.0000 -0.0000 -0.0000</p><p>　　-0.0000 -0.0000 0.0000</p><p><b>　　迭代中的△Q：</b></p><p>　　-0.1700 1.1478</p>&

75、lt;p>　　-0.0016 -0.0427</p><p>　　-0.0000 -0.0001</p><p>　　-0.0000 -0.0000</p><p>　　-0.0000 -0.0000</p><p><b>　　迭代中相角:</b></p><p>　　-0

76、.1980 -0.3723 0.7258</p><p>　　-0.2277 -0.3731 0.6838</p><p>　　-0.2265 -0.3725 0.6853</p><p>　　-0.2265 -0.3725 0.6853</p><p>　　-0.2265 -0.3725 0

77、.6853</p><p><b>　　迭代中電壓的模:</b></p><p>　　1.0007 1.0345</p><p>　　1.0002 1.0332</p><p>　　1.0002 1.0332</p><p>　　1.0002 1.0332</p>

78、;<p>　　1.0002 1.0332</p><p><b>　　平衡結(jié)點(diǎn)的功率:</b></p><p>　　0.2141 + 0.8491i</p><p><b>　　全部線路功率分布:</b></p><p>　　0 -0.0629 - 0.284

79、4i -0.2371 + 0.0544i 0 </p><p>　　-0.0341 - 0.0019i 0 -0.1245 - 0.0166i -0.2407 - 0.6106i</p><p>　　0.3354 + 0.1157i 0.1255 - 0.0432i 0

80、 0.0331 - 0.2774i</p><p>　　0 0.1411 + 0.3526i 0.0730 + 0.4965i 0 </p><p><b>　　思考題</b></p><p>　　1）.潮流計(jì)算的方法有哪些？各有何特點(diǎn)？</p><p>　

81、　2）.設(shè)計(jì)中遇到的問題和解決的辦法</p><p><b>　　答：</b></p><p>　　1牛頓-拉夫遜法及快速分解法?？焖俜纸夥ㄓ袃蓚€(gè)主要特點(diǎn)：a降壓，在潮流計(jì)算的修正方程中利用了有功功率主要與節(jié)點(diǎn)電壓相位有關(guān)，無功功率主要與節(jié)點(diǎn)電壓幅值有關(guān)的特點(diǎn)，實(shí)現(xiàn)P-Q分解。b因子表固定化，利用了線路兩端電壓相位差不大的假定，是修正方程系數(shù)矩陣元素變?yōu)槌?shù)，并且就是

82、節(jié)點(diǎn)導(dǎo)納的虛部。 </p><p>　　2設(shè)計(jì)過程中剛開始程序出錯(cuò)，運(yùn)行不出結(jié)果，后經(jīng)問老師和同學(xué)，改正了程序的錯(cuò)誤之處，得出正確的結(jié)果，設(shè)計(jì)程序過程中要認(rèn)真研究流程圖 。</p><p><b>　　課程設(shè)計(jì)總結(jié) </b></p><p>　　此次課程設(shè)計(jì)使我在潮流計(jì)算、MATLAB的使用方面均有所提高，但也暴露出了一些問題；理論知

83、識的儲備還是不足，對MATLAB的性能和特點(diǎn)還不能有一個(gè)全面的把握，對MATLAB中PSAT工具箱使用不夠熟悉，相信通過以后的學(xué)習(xí)能彌補(bǔ)這些不足，從而達(dá)到一個(gè)新的層次。</p><p>　　潮流計(jì)算是電力系統(tǒng)的最基本、最常用的分析計(jì)算。用以研究系統(tǒng)規(guī)劃和運(yùn)行中提出的各種問題。對規(guī)劃中的電力系統(tǒng)，通過潮流計(jì)算可以檢驗(yàn)所提出的電力系統(tǒng)規(guī)劃方案能否滿足各種運(yùn)行方式的要求；對運(yùn)行中的電力系統(tǒng)，</p>&l

84、t;p>　　通過潮流計(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>　　實(shí)際的情況遠(yuǎn)比我們計(jì)算的情況復(fù)雜，這讓我深刻了解了潮流計(jì)算的重要性。精準(zhǔn)的潮流計(jì)算不僅可以使電網(wǎng)處于穩(wěn)定且平衡的狀態(tài)，這對電網(wǎng)的安全運(yùn)行起到關(guān)鍵性的作用，還可以實(shí)現(xiàn)運(yùn)行

85、的經(jīng)濟(jì)性。</p><p>　　所以我認(rèn)為學(xué)好電力系統(tǒng)的關(guān)鍵在于學(xué)好潮流計(jì)算，電力系統(tǒng)的安全運(yùn)行不是兒戲，一個(gè)很小的失誤就可能造成整個(gè)電網(wǎng)的崩潰，這給國家和社會帶來的危害和損失之大是可想而知的。學(xué)好電力系統(tǒng)這門課、學(xué)好潮流計(jì)算是我們電專業(yè)學(xué)生的當(dāng)務(wù)之急。沒有堅(jiān)實(shí)的基礎(chǔ)就不可能建造出高樓大廈，所以，要想在以后的工作崗位上有所成就，就必須打好基礎(chǔ)，努力提高自己的專業(yè)素養(yǎng)。</p><p>&l

86、t;b>　　參考文獻(xiàn)</b></p><p>　　[1] 陳珩.《電力系統(tǒng)穩(wěn)態(tài)分析》.中國電力出版社.2007年.第二版；</p><p>　　[2] 張伯明、陳壽孫、嚴(yán)正.《高等電力網(wǎng)絡(luò)分析》.清華大學(xué)出版社.2007年.第二版；</p><p>　　[3] 西安交通大學(xué)等合編.《電力系統(tǒng)計(jì)算：電子數(shù)字計(jì)算機(jī)的應(yīng)用》.北京：水利電力出版社；<

87、;/p><p>　　[4] 王錫凡主編.《現(xiàn)代電力系統(tǒng)分析》.科學(xué)出版社；</p><p>　　[5] 祝淑萍等.《電力系統(tǒng)分析課程設(shè)計(jì)與綜合實(shí)驗(yàn)》.中國電力出版社.2007年3月.第一版；</p><p>　　[6] 蘇金明、王永利主編.《Matlab實(shí)用指南》.電子工業(yè)出版社。</p><p><b>　　附錄</b>&