【正文】 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。平衡結(jié)點(diǎn)的功率:39。)。disp(PH)。disp(39。全部線路功率分布:39。)。disp(SS)。 程序結(jié)果請(qǐng)輸入結(jié)點(diǎn)數(shù)：n=4請(qǐng)輸入PV結(jié)點(diǎn)數(shù)：n1=1請(qǐng)輸入PQ結(jié)點(diǎn)數(shù)：n2=2請(qǐng)輸入平衡結(jié)點(diǎn)：isb=4請(qǐng)輸入精確度：pr=請(qǐng)輸入變比矩陣看:K=[0 0 0 0。0 0 0 。0 0 0 0。0 0 0 0]請(qǐng)輸入支路阻抗矩陣：C=[0 + + inf。 + 0 + +。+ + 0 +。inf + + 0]請(qǐng)輸入支路導(dǎo)納矩陣：y=[0 0。 0 0 0。 0 0 0。0 0 0 0]請(qǐng)輸入結(jié)點(diǎn)電壓矩陣：U=[1+0i 1+0i +0i +0i]請(qǐng)輸入各結(jié)點(diǎn)的功率：S=[ 0]結(jié)點(diǎn)導(dǎo)納矩陣: + + 0 + + 0 + + + + 0 0 + + 迭代中的雅克比矩陣: 迭代中的△P： 迭代中的△Q： 迭代中相角: 迭代中電壓的模: 平衡結(jié)點(diǎn)的功率: + 全部線路功率分布: 0 + 0 0 0 0 + + 0 第七章 總 結(jié) 此次課程設(shè)計(jì)首先讓我明白了要使電力系統(tǒng)運(yùn)行的穩(wěn)定，必須經(jīng)過(guò)精密的設(shè)計(jì)和計(jì)算。在進(jìn)行課題設(shè)計(jì)的過(guò)程中，加深了我對(duì)潮流計(jì)算的認(rèn)識(shí)，尤其是對(duì)牛頓拉夫遜潮流計(jì)算的求解思路有了比較透徹的理解。同時(shí)由于求解過(guò)程中用到求節(jié)點(diǎn)導(dǎo)鈉矩陣，求矩陣的逆等等，又使我對(duì)以前所學(xué)的知識(shí)有了一次很好的溫習(xí)。同時(shí)也看到了研究性學(xué)習(xí)的效果，從研究中去學(xué)習(xí)，理論結(jié)合實(shí)際，將理論運(yùn)用到實(shí)際，同時(shí)在實(shí)踐中發(fā)現(xiàn)問(wèn)題，然后解決問(wèn)題， 通過(guò)這次課程設(shè)計(jì)我在MATLAB編程，潮流計(jì)算，WORD文檔的編輯方面均有所提高，但也暴露出了一些問(wèn)題；我發(fā)現(xiàn)了自己的基礎(chǔ)知識(shí)有很多的不足，對(duì)MATLAB的性能和特點(diǎn)還不能有一個(gè)全面的把握，這些基礎(chǔ)的缺乏給我的設(shè)計(jì)計(jì)劃造成了不小的障礙，相信通過(guò)以后的學(xué)習(xí)能彌補(bǔ)這些不足，從而達(dá)到一個(gè)新的層次。 第八章 參考文獻(xiàn)[1] [M]. 武漢：華中理工大學(xué)出版社， [2] [M].北京：水利電力出版社， [3] 陳 [M].北京：水利電力出版社， [4] [M].北京： 水利電力出版社， [5] 于永源，楊綺雯. 電力系統(tǒng)分析（第二版）[M]. 北京：中國(guó)電力出版社，30