【正文】 m)。 else H(m,R)=U1(m)*U1(R)*(G(m,R)*sin(O(m)O(R))B(m,R)*cos(O(m)O(R)))。 end end end for m=1:n1+n2 for R=1:n2 if m==R N(m,m)=U1(m)^2*G(m,m)PT1(m)。 else N(m,R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)O(R))+B(m,R)*sin(O(m)O(R)))。 end end end for m=1:n2 for R=1:n1+n2 if m==R J(m,m)=U1(m)^2*G(m,m)PT1(m)。 else J(m,R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)O(R))+B(m,R)*sin(O(m)O(R)))。 end end end for m=1:n2 for R=1:n2 if m==R L(m,m)=U1(m)^2*B(m,m)QT1(m)。 else L(m,R)=U1(m)*U1(R)*(G(m,R)*sin(O(m)O(R))B(m,R)*cos(O(m)O(R)))。 end end end JJ=[H N。J L]。 disp(JJ)。 PQ=[PP39。39。]。 DA=inv(JJ)*PQ。 DA1=DA39。 for m=1:n1+n2 OO(m)=DA1(m)。 end for m=n:n1+n2+n2 UU1(mn1n2)=DA1(m)。 end UD2=diag(UD)。 UU=UU1*UD2。 for m=1:n1+n2 O(m)=O(m)+OO(m)。 end for m=1:n2 U1(m)=U1(m)+UU(m)。 end for m=1:n1+n2 o(k+1,m)=180/pi*O(m)。 end for m=1:n2 u(k+1,m)=U1(m)。 end k=k+1。endfor m=1:n b(m)=U1(m)*cos(O(m))。 c(m)=U1(m)*sin(O(m))。endU=b+i*c。for R=1:n PH1(R)=U(isb)*conj(Y(isb,R))*conj(U(R))。endPH=sum(PH1)。for m=1:n for R=1:n if m~=R C1(m,R)=1/C(m,R)。 else C1(m,m)=C(m,m)。 end endendfor m=1:n for R=1:n if (C(m,R)~=inf)amp。(m~=R) SS(m,R)=U1(m)^2*conj(C1(m,m))+U(m)*(conj(U(m))conj(U(R)))*conj(C1(m,R))。 end endenddisp(39。迭代中的△P：39。)。disp(PP1)。disp(39。迭代中的△Q：39。)。disp(1)。disp(39。迭代中相角:39。)。disp(o)。disp(39。迭代中電壓的模:39。)。disp(u)。disp(39。平衡結點的功率:39。)。disp(PH)。disp(39。全部線路功率分布:39。)。disp(SS)。 程序結果請輸入結點數(shù)：n=4請輸入PV結點數(shù)：n1=1請輸入PQ結點數(shù)：n2=2請輸入平衡結點：isb=4請輸入精確度：pr=請輸入變比矩陣看:K=[0 0 0 0。0 0 0 。0 0 0 0。0 0 0 0]請輸入支路阻抗矩陣：C=[0 + + inf。 + 0 + +。+ + 0 +。inf + + 0]請輸入支路導納矩陣：y=[0 0。 0 0 0。 0 0 0。0 0 0 0]請輸入結點電壓矩陣：U=[1+0i 1+0i +0i +0i]請輸入各結點的功率：S=[ 0]結點導納矩陣: + + 0 + + 0 + + + + 0 0 + + 迭代中的雅克比矩陣: 迭代中的△P： 迭代中的△Q： 迭代中相角: 迭代中電壓的模: 平衡結點的功率: + 全部線路功率分布: 0 + 0 0 0 0 + + 0 第七章 總 結 此次課程設計首先讓我明白了要使電力系統(tǒng)運行的穩(wěn)定，必須經過精密的設計和計算。在進行課題設計的過程中，加深了我對潮流計算的認識，尤其是對牛頓拉夫遜潮流計算的求解思路有了比較透徹的理解。同時由于求解過程中用到求節(jié)點導鈉矩陣，求矩陣的逆等等，又使我對以前所學的知識有了一次很好的溫習。同時也看到了研究性學習的效果，從研究中去學習，理論結合實際，將理論運用到實際，同時在實踐中發(fā)現(xiàn)問題，然后解決問題， 通過這次課程設計我在MATLAB編程，潮流計算，WORD文檔的編輯方面均有所提高，但也暴露出了一些問題；我發(fā)現(xiàn)了自己的基礎知識有很多的不足，對MATLAB的性能和特點還不能有一個全面的把握，這些基礎的缺乏給我的設計計劃造成了不小的障礙，相信通過以后的學習能彌補這些不足，從而達到一個新的層次。 第八章 參考文獻[1] [M]. 武漢：華中理工大學出版社， [2] [M].北京：水利電力出版社， [3] 陳 [M].北京：水利電力出版社， [4] [M].北京： 水利電力出版社， [5] 于永源，楊綺雯. 電力系統(tǒng)分析（第二版）[M]. 北京：中國電力出版社，30