【正文】 (i,3)* B1(i,5)~2)+B1(i,4)./2。 end %求導(dǎo)納矩陣 G=real(Y)。 BI(i,i)= BI(i,i)+ B2(i,5)。 f(i)=imag(B2(i,3))。 end end for i=2:n if i= =n B(i,i)=1./B(i,i)。 for k=i+1:n for j1=i+1:n B(k,j1)=B(k,j1)B(k,i)*B(i,j1)。k=0。 else k=i+1。 end end end end ICT2=1。DET=0。 for i=1:n if i~=isb C(i)=0。 DET=abs(DP1(i))。 if i~ =n IC1=i+1。 for MZ=2:IC4 I=IC4+2MZ。L=0。 end DQ1(i)=Q(i)V(i)*C(i)。 end end end end else kp=0。L=L+1。 DET=abs(DQ1(i))。 if i= =na for LZ=2:i L=i+2LZ。 end end else IC1=i+1。 V(i)=V(i)DQ(L)。 if kp~=0 K=K+1。 O(k)=O(k)*180./pi。 disp(?各節(jié)點(diǎn)的電壓相角 O（節(jié)點(diǎn)號(hào)從小到大排）：‘ )為； disp(O)。 end disp(?各節(jié)點(diǎn)的功率 S（節(jié)點(diǎn)號(hào)從小到大排）：‘ )為； disp（ S）； disp(‘各條支路的首端功率 Si 為（順序同您輸入 B1時(shí)的一樣）：‘）為； for i=1:nl if B1(i,6)= =0 p=B1(i,1)。 end Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。 else p=B1(i,2)。 end disp(‘各條支路的功率損耗 DS（順序同您輸入 B1 時(shí)的一樣）：‘）為； for i=1:nl if B1(i,6)= =0 p=B1(i,1)。 end DS(i)=Si(p,q)+Sj(q,p)。 end end disp(‘以下是每次迭代后各節(jié)點(diǎn)的電壓值（如圖所示）‘ )； plot(Cs,Dy),xlable(?迭代次數(shù)‘ )， ylable(?電壓‘ ),title(?電壓迭代次數(shù)曲線’ )； 設(shè)計(jì)驗(yàn)證 現(xiàn)以一個(gè)例子來說明該軟件得 使用方法 : 電力系統(tǒng)接線，如 ,試求潮流分布 （ 1） n為節(jié)點(diǎn)數(shù)、 nl 為支路數(shù)、 isb為平衡母線節(jié)點(diǎn)號(hào)（固定為 1）、 pr為誤差精度。2 3 + 1 0。0 + 0 0 2。2 0。 總結(jié) 通過學(xué)習(xí)潮流計(jì)算 原理參閱一些參考資料及文獻(xiàn)，對(duì)傳統(tǒng)潮流死算有所理解，對(duì)一些新的潮流計(jì)算技術(shù)有了整體的認(rèn)識(shí)。 參考文獻(xiàn)及資料： 《電力系統(tǒng)穩(wěn)態(tài)分析》，陳珩，中國電力出版社， 2021，第三版 《電力系統(tǒng)分析》，韓禎祥，浙江大學(xué)出版社， 2021，第三版 《電力系統(tǒng)分析》，蘇小林 閻曉霞 ，中國電力出版社 《電力系統(tǒng) 分析課程實(shí)際設(shè)計(jì)與綜合實(shí)驗(yàn) 》 ，祝書萍，中國電力出版社， 2021，第一版 清華大學(xué)等校合編 .電子數(shù)字計(jì)算機(jī)的應(yīng)用 — 電力系統(tǒng)計(jì)算 [M].北京 :水利電力出版社，1978 朱仁峰 .精通 Matlab7[M].北京：清華大學(xué)出版社， 2021 《火車發(fā)電廠、變電所二次接線技術(shù)規(guī)定（強(qiáng)電 部分）》 SDGJ─78 《 MALAB 科學(xué)運(yùn)算完整解決法案 》 ，精銳創(chuàng)組編，人民郵電大學(xué)出版社 《 MALAB 應(yīng)用集錦 》 ，林雪松 周婧 林德新 朱仁峰 .精通 Matlab7[M].北京：清華大學(xué)出版社， 2021 。 本文分析和比較了電力系統(tǒng)潮流計(jì)算的各種方法，其中 PQ 分解法具有計(jì)算速度快、所占內(nèi)存少等優(yōu)勢(shì)被廣泛采用。4 0。0 + 0 0 2。3 4 + 1 0。 Figure schematic （ 3）請(qǐng)輸入各節(jié)點(diǎn)參數(shù) 形成的矩陣 B2 矩陣 B2的每行是由下列參數(shù)構(gòu)成的： ① 節(jié)點(diǎn)所接發(fā)電機(jī)的功率 SG； ② 節(jié)點(diǎn)負(fù)荷的功率 SL； ③ 節(jié)點(diǎn)電壓的初始值； ④ PU 節(jié)點(diǎn)電壓 U 的給定值； ⑤ 節(jié)點(diǎn)所接的無功補(bǔ)償設(shè)備的容量； ⑥ 節(jié)點(diǎn)分類標(biāo)號(hào) igl。 end for i=1:K Cs(i)=i。 else p=B1(i,2)。 end Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))conj(E(p)))*conj(1./(B1(i,3)*B1(i, 5))))。 end disp(‘各條支路的末端功率 Sj（順序同您輸入 B1 時(shí)的一樣）：‘）為； for i=1:nl if B1(i,6)= =0 p=B1(i,1)。 else p=B1(i,2)。 for q=1:n C(p)= C(p)+conj(Y(p,q))*conj(E(q))。 disp(E)。 end end disp(?迭代次數(shù)‘ )； disp（ K）； disp（‘每次沒有達(dá)到精度要求的有功功率個(gè)數(shù)為‘）； disp（ NP）； disp（‘每次沒有達(dá)到精度要求的無功功率個(gè)數(shù)為‘）； disp（ Nq） 。 K=K+1。 end end end L=0。 for MZ=1:IC4 I=IC4+1MZ。 if ICT3~=0 L=0。 end DQ1(i)=Q(i)V(i)*C(i)。 L=0。 DET=abs(DQ1(i))。L=L+1。 end end end end for i=2:n O(i)=O(i)DP(i)。 end else for LZ=3:i L=i+3LZ。 end end end Np(K)=ICT2。 end DP1(i)=P(i)V(i)*C(i)。 while ICT2~=0|ICT3~=0 ICT2=0。kp=1,kq=1。 end A(i,i)=1. /A(i,i)。 A(p,k)=BI(i,j1)。q=0。 for j1=IC1:n B(i,j1)= B(i,j1)./B(i,i)。 end for i=1:n if B2(i,6)= =2 V(i)=sqrt(e(i)~2+f(i)~2)。Q=imag(S)。BI=imag(Y)。 Y(p,p)= Y(p,p) +1./B1(i,3)+ B1(i,4)./2。 Y(q,p) =Y(p,q)。 else p= B1(i,2)。 p=x(i,1)。f= zeros(1,n)。 na==input(‘請(qǐng)輸入 PQ節(jié)點(diǎn)號(hào) :na=?)。 pr==input(‘請(qǐng)輸入誤差精度 :pr=’ )。由于 MATLAB語言功能強(qiáng)大、人際界面友好、編程效率高、強(qiáng)大而智能化的作圖功能，且具有編程語句簡潔、靈活、表達(dá)和運(yùn)算能力強(qiáng)等顯著特點(diǎn)。 end k=k+1。%求逆矩陣 C=max(abs(dU))。 JJ(2*m,2*n)=L(m,n)。 L(m,n)=G(m,n)*F(m)B(m,n)*E(m)。 JJ(2*m,2*m)=L(m,m)。 L(m,m)=2*B(m,m)*E(m)(sum(Bi)(B(m,m)*E(m)+G(m,m)*F(m)))。 Ai(n)=G(m,n)*E(n)B(m,n)*F(n)。 end %設(shè) P,Q的改變量為 dP,dQ dP(m)=P(m)s