【導(dǎo)讀】潮流分析是研究電力系統(tǒng)的一種最基本和最重要的計(jì)算。等人編制了實(shí)際可行的計(jì)算機(jī)潮流計(jì)算程序。力系統(tǒng)提供了極其有力的計(jì)算手段。算是計(jì)算系統(tǒng)動(dòng)態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)。運(yùn)行方式的合理性、可靠性和經(jīng)濟(jì)性。潮流有其獨(dú)特優(yōu)勢(shì)。牛頓-拉夫遜法的優(yōu)點(diǎn)是收斂速度快，若。非常精確的解，而且其迭代次數(shù)與所計(jì)算的網(wǎng)絡(luò)規(guī)?；緹o(wú)關(guān)。牛頓法的缺點(diǎn)是每次迭代的計(jì)算量和所需的內(nèi)存量較大。不過(guò)，內(nèi)存占用量及每次迭代所需的時(shí)間與程序設(shè)計(jì)技巧密切相關(guān)。所以對(duì)各節(jié)電可以采用統(tǒng)一的電壓初值。的物理特點(diǎn)，對(duì)極坐標(biāo)形式的牛頓-拉夫遜法修正方程式進(jìn)行了合理的簡(jiǎn)化。另外，MATLAB提供了一種特殊的工具：。但需注意進(jìn)行數(shù)除。于前面講的數(shù)學(xué)運(yùn)算，為有所區(qū)別，我們稱(chēng)之為數(shù)組運(yùn)算。關(guān)系”的規(guī)定，數(shù)組與常數(shù)之間的除法運(yùn)算沒(méi)有任何限制。有了“對(duì)應(yīng)元素”的規(guī)定，數(shù)組的運(yùn)算實(shí)質(zhì)上就是針對(duì)數(shù)。組內(nèi)部的每個(gè)元素進(jìn)行的。矩陣的冪運(yùn)算與數(shù)組的冪運(yùn)算有很大的區(qū)別。

　　

【正文】 =1:n if B2(i,6)= =2 L=L+1。 V(i)=V(i)DQ(L)。 end end kp=1。 K=K+1。 else kq=0。 if kp~=0 K=K+1。 end end for i=1:n Dy(K1,i)=V(i)。 end end disp(?迭代次數(shù)‘ )； disp（ K）； disp（‘每次沒(méi)有達(dá)到精度要求的有功功率個(gè)數(shù)為‘）； disp（ NP）； disp（‘每次沒(méi)有達(dá)到精度要求的無(wú)功功率個(gè)數(shù)為‘）； disp（ Nq） 。 for k=1:n E(k)=V(k)*cos(O(K))+V(k)*sin(O(k))*j。 O(k)=O(k)*180./pi。 end disp(?各節(jié)點(diǎn)的電壓標(biāo)么值 E為 （ 節(jié)點(diǎn)號(hào)從小到大排 ）：‘ )。 disp(E)。 disp(?各節(jié)點(diǎn)的電壓 U大?。ü?jié)點(diǎn)號(hào)從小到大排）：‘ )為； disp(U)。 disp(?各節(jié)點(diǎn)的電壓相角 O（節(jié)點(diǎn)號(hào)從小到大排）：‘ )為； disp(O)。 for p=1:n C(p)=0。 for q=1:n C(p)= C(p)+conj(Y(p,q))*conj(E(q))。 end S(p)=E(p)*C(p)。 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)。q=B1(i,2)。 else p=B1(i,2)。q=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))))。 disp(Si(p,q))。 end disp(‘各條支路的末端功率 Sj（順序同您輸入 B1 時(shí)的一樣）：‘）為； for i=1:nl if B1(i,6)= =0 p=B1(i,1)。q=B1(i,2)。 else p=B1(i,2)。q=B1(i,1)。 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))))。 disp(Sj(q,p))。 end disp(‘各條支路的功率損耗 DS（順序同您輸入 B1 時(shí)的一樣）：‘）為； for i=1:nl if B1(i,6)= =0 p=B1(i,1)。q=B1(i,2)。 else p=B1(i,2)。q=B1(i,1)。 end DS(i)=Si(p,q)+Sj(q,p)。 disp(DS(i))。 end for i=1:K Cs(i)=i。 for j=1:n Dy(K,j)=Dy(K1,j)。 end end disp(‘以下是每次迭代后各節(jié)點(diǎn)的電壓值（如圖所示）‘ )； plot(Cs,Dy),xlable(?迭代次數(shù)‘ )， ylable(?電壓‘ ),title(?電壓迭代次數(shù)曲線’ )； 設(shè)計(jì)驗(yàn)證 現(xiàn)以一個(gè)例子來(lái)說(shuō)明該軟件得 使用方法 : 電力系統(tǒng)接線，如 ,試求潮流分布 （ 1） n為節(jié)點(diǎn)數(shù)、 nl 為支路數(shù)、 isb為平衡母線節(jié)點(diǎn)號(hào)（固定為 1）、 pr為誤差精度。 （ 2）輸入由支路參數(shù)形成的矩陣 B1 矩陣 B1的每行是由下列參數(shù)構(gòu)成的： ① 某支路的首端號(hào) P； ② 末端號(hào) Q，且 P Q； ③ 支路的阻抗（ R +jX）； ④ 支路的對(duì)地容抗； ⑤ 支路的變比 K； ⑥ 折算到哪一側(cè)的標(biāo)志（如果支路的首端 P 處于高壓側(cè)則請(qǐng)輸入 “1”，否則請(qǐng)輸入 “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)所接的無(wú)功補(bǔ)償設(shè)備的容量； ⑥ 節(jié)點(diǎn)分類(lèi)標(biāo)號(hào) igl。 輸入數(shù)據(jù)： 請(qǐng)輸入節(jié)點(diǎn)數(shù)： n=5 請(qǐng)輸入支路數(shù)： nl=5 請(qǐng)輸入平衡母線節(jié)點(diǎn)號(hào)： isb=1 請(qǐng)輸入誤差精度： pr= 請(qǐng)輸入由支路參數(shù)形成的矩陣： B1=[1 2 0 。2 3 + 1 0。2 4 + 0 1 0。3 4 + 1 0。3 5 0 1] 請(qǐng)輸入由節(jié)點(diǎn)參數(shù)形成的矩陣： B2=[0 0 0 1。0 + 0 0 2。0 2+1i 0 0 2。0 + 0 0 2。5 0 0 3] 請(qǐng)輸入由節(jié)點(diǎn)號(hào)及其對(duì)地阻抗形成的矩陣： X=[1 0。2 0。3 0。4 0。5 0] 請(qǐng)輸入 PQ節(jié)點(diǎn)數(shù) na=3 輸入框圖 在輸入需要的數(shù)據(jù)后 ,點(diǎn)擊運(yùn)算 ,出現(xiàn)結(jié)果如下： 結(jié)果： 迭代次數(shù) 10 每次沒(méi)有達(dá)到精度要求的有功功率個(gè)數(shù)為 4 4 4 4 4 4 4 4 3 0 每次沒(méi)有達(dá)到精度要求的無(wú)功功率個(gè)數(shù)為 3 3 3 3 3 3 3 3 1 0 各節(jié)點(diǎn)的電壓標(biāo)么值 E（節(jié)點(diǎn)號(hào)從小到大排）為： + + 各節(jié)點(diǎn)的電壓 U大?。ü?jié) 點(diǎn)號(hào)從小到大排）為： 各節(jié)點(diǎn)的電壓相角 O（節(jié)點(diǎn)號(hào)從小到大排）為： 0 各節(jié)點(diǎn)的功率 S（節(jié)點(diǎn)號(hào)從小到大排）為： + + 各條支路的首端功率 Si（順序同您輸入 B1 時(shí)的一樣）為 ： + + + + + 各條支路的末端功率 Sj（順序同您輸入 B1 時(shí)的一樣）為： 各條支路的功率損耗 DS（順序同您輸入 B1 時(shí)的一樣）為： + + + + Figure Simulation diagram 經(jīng)過(guò)手動(dòng)計(jì)算，得 1至 5號(hào)節(jié)點(diǎn)功率分別為： + + 舉例驗(yàn)算結(jié)論：通過(guò)手動(dòng)運(yùn)算，發(fā)現(xiàn)運(yùn)算結(jié)果與手動(dòng)結(jié)果運(yùn)算結(jié)果一樣，故軟件計(jì)算正確。 總結(jié) 通過(guò)學(xué)習(xí)潮流計(jì)算 原理參閱一些參考資料及文獻(xiàn)，對(duì)傳統(tǒng)潮流死算有所理解，對(duì)一些新的潮流計(jì)算技術(shù)有了整體的認(rèn)識(shí)。 潮流計(jì)算是電力系統(tǒng)的一項(xiàng)基本計(jì)算，它根據(jù)給定的運(yùn)行條件及系統(tǒng)接線情況確定系統(tǒng)的運(yùn)行狀態(tài)、母線的電壓，各元件中流過(guò)的功率，系統(tǒng)的功率損耗等，進(jìn)而比較運(yùn)行方式的合理性 、 可靠性和經(jīng)濟(jì)性 。 本文分析和比較了電力系統(tǒng)潮流計(jì)算的各種方法，其中 PQ 分解法具有計(jì)算速度快、所占內(nèi)存少等優(yōu)勢(shì)被廣泛采用。本設(shè)計(jì)首先進(jìn)行了 PQ 分解法的理論研究，然后 應(yīng)用 MATLAB 軟件編制 PQ 分解法潮流程序，并通過(guò)一個(gè)電力系統(tǒng)實(shí)例對(duì)所編程序進(jìn)行驗(yàn)證，得出結(jié) 論。 參考文獻(xiàn)及資料： 《電力系統(tǒng)穩(wěn)態(tài)分析》，陳珩，中國(guó)電力出版社， 2021，第三版 《電力系統(tǒng)分析》，韓禎祥，浙江大學(xué)出版社， 2021，第三版 《電力系統(tǒng)分析》，蘇小林 閻曉霞 ，中國(guó)電力出版社 《電力系統(tǒng) 分析課程實(shí)際設(shè)計(jì)與綜合實(shí)驗(yàn) 》 ，祝書(shū)萍，中國(guó)電力出版社， 2021，第一版 清華大學(xué)等校合編 .電子數(shù)字計(jì)算機(jī)的應(yīng)用 — 電力系統(tǒng)計(jì)算 [M].北京 :水利電力出版社，1978 朱仁峰 .精通 Matlab7[M].北京：清華大學(xué)出版社， 2021 《火車(chē)發(fā)電廠、變電所二次接線技術(shù)規(guī)定（強(qiáng)電 部分）》 SDGJ─78 《 MALAB 科學(xué)運(yùn)算完整解決法案 》 ，精銳創(chuàng)組編，人民郵電大學(xué)出版社 《 MALAB 應(yīng)用集錦 》 ，林雪松 周婧 林德新 朱仁峰 .精通 Matlab7[M].北京：清華大學(xué)出版社， 2021