畢業(yè)設(shè)計-matlab電力系統(tǒng)pq潮流計算程序設(shè)計-資料下載頁
2025-08-06 11:39本頁面
【正文】 (i)=0。L=L+1。 for k=1:n C(i)= C(i) +V(k)*(G(i,k)*sin(O(i)O(k))BI(i,k)*cos(O(i)O(k)))。 end DQ1(i)=Q(i)V(i)*C(i)。 DQ(L)= DQ1(i)./V(i)。 DET=abs(DQ1(i))。endendendendNq(K)=ICT3。if ICT3~=0 L=0。 for i=1:naDQ(i)=A(i,i)*DQ(i)。if i= =na for LZ=2:i L=i+2LZ。 IC4=L1。 for MZ=1:IC4 I=IC4+1MZ。 DQ(I)=DQ(I)A(I,L)*DQ(L)。 endendelse IC1=i+1。 for k=IC1:na DQ(k)=DQ(k)A(k,i)*DQ(i)。endendendL=0。for i=1:nif B2(i,6)= =2 L=L+1。 V(i)=V(i)DQ(L)。endend kp=1。K=K+1。else kq=0。 if kp~=0K=K+1。 endendfor i=1:n Dy(K1,i)=V(i)。endenddisp(‘迭代次數(shù)‘)；disp（K）；disp（‘每次沒有達到精度要求的有功功率個數(shù)為‘）；disp（NP）；disp（‘每次沒有達到精度要求的無功功率個數(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。enddisp(‘各節(jié)點的電壓標(biāo)么值E為（節(jié)點號從小到大排）：‘)。disp(E)。disp(‘各節(jié)點的電壓U大?。ü?jié)點號從小到大排）：‘)為；disp(U)。disp(‘各節(jié)點的電壓相角O（節(jié)點號從小到大排）：‘)為；disp(O)。for p=1:n C(p)=0。 for q=1:nC(p)= C(p)+conj(Y(p,q))*conj(E(q))。 end S(p)=E(p)*C(p)。enddisp(‘各節(jié)點的功率S（節(jié)點號從小到大排）：‘)為；disp（S）；disp(‘各條支路的首端功率Si為（順序同您輸入B1時的一樣）：‘）為；for i=1:nl if B1(i,6)= =0p=B1(i,1)。q=B1(i,2)。else p=B1(i,2)。q=B1(i,1)。 endSi(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))。enddisp(‘各條支路的末端功率Sj（順序同您輸入B1時的一樣）：‘）為；for i=1:nl if B1(i,6)= =0p=B1(i,1)。q=B1(i,2)。else p=B1(i,2)。q=B1(i,1)。 endSj(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))。enddisp(‘各條支路的功率損耗DS（順序同您輸入B1時的一樣）：‘）為；for i=1:nlif B1(i,6)= =0p=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))。endfor i=1:K Cs(i)=i。 for j=1:nDy(K,j)=Dy(K1,j)。 endenddisp(‘以下是每次迭代后各節(jié)點的電壓值（如圖所示）‘)；plot(Cs,Dy),xlable(‘迭代次數(shù)‘)，ylable(‘電壓‘),title(‘電壓迭代次數(shù)曲線’)；現(xiàn)以一個例子來說明該軟件得使用方法:電力系統(tǒng)接線，,試求潮流分布（1）n為節(jié)點數(shù)、nl為支路數(shù)、isb為平衡母線節(jié)點號（固定為 1）、pr為誤差精度。（2）輸入由支路參數(shù)形成的矩陣 B1矩陣 B1的每行是由下列參數(shù)構(gòu)成的：① 某支路的首端號 P；② 末端號 Q，且 P Q；③ 支路的阻抗（R +jX）；④ 支路的對地容抗；⑤ 支路的變比 K；⑥ 折算到哪一側(cè)的標(biāo)志（如果支路的首端 P 處于高壓側(cè)則請輸入“1”，否則請輸入“0”）。圖 Figure schematic（3）請輸入各節(jié)點參數(shù)形成的矩陣 B2矩陣 B2的每行是由下列參數(shù)構(gòu)成的：① 節(jié)點所接發(fā)電機的功率 SG；② 節(jié)點負荷的功率 SL；③ 節(jié)點電壓的初始值；④ PU 節(jié)點電壓 U 的給定值；⑤ 節(jié)點所接的無功補償設(shè)備的容量；⑥ 節(jié)點分類標(biāo)號 igl。輸入數(shù)據(jù)：請輸入節(jié)點數(shù)：n=5請輸入支路數(shù)：nl=5請輸入平衡母線節(jié)點號：isb=1請輸入誤差精度：pr=請輸入由支路參數(shù)形成的矩陣：B1=[1 2 0 。2 3 + 1 0。2 4 + 0 1 0。3 4 + 1 0。3 5 0 1]請輸入由節(jié)點參數(shù)形成的矩陣：B2=[0 0 0 1。0 + 0 0 2。0 2+1i 0 0 2。0 + 0 0 2。5 0 0 3]請輸入由節(jié)點號及其對地阻抗形成的矩陣：X=[1 0。2 0。3 0。4 0。5 0]請輸入PQ節(jié)點數(shù)na=3: 輸入框圖 Input Block Diagram在輸入需要的數(shù)據(jù)后,點擊運算,出現(xiàn)結(jié)果如下：結(jié)果：迭代次數(shù)10每次沒有達到精度要求的有功功率個數(shù)為4 4 4 4 4 4 4 4 3 0每次沒有達到精度要求的無功功率個數(shù)為3 3 3 3 3 3 3 3 1 0各節(jié)點的電壓標(biāo)么值E（節(jié)點號從小到大排）為： + +各節(jié)點的電壓U大?。ü?jié)點號從小到大排）為： 各節(jié)點的電壓相角O（節(jié)點號從小到大排）為：0 各節(jié)點的功率S（節(jié)點號從小到大排）為：+ +各條支路的首端功率Si（順序同您輸入B1時的一樣）為： + + + + +各條支路的末端功率Sj（順序同您輸入B1時的一樣）為： 各條支路的功率損耗DS（順序同您輸入B1時的一樣）為： + + + +：Figure Simulation diagram經(jīng)過手動計算，得1至5號節(jié)點功率分別為：++舉例驗算結(jié)論：通過手動運算，發(fā)現(xiàn)運算結(jié)果與手動結(jié)果運算結(jié)果一樣，故軟件計算正確。5 總結(jié)通過學(xué)習(xí)潮流計算原理參閱一些參考資料及文獻，對傳統(tǒng)潮流死算有所理解，對一些新的潮流計算技術(shù)有了整體的認識。潮流計算是電力系統(tǒng)的一項基本計算，它根據(jù)給定的運行條件及系統(tǒng)接線情況確定系統(tǒng)的運行狀態(tài)、母線的電壓，各元件中流過的功率，系統(tǒng)的功率損耗等，進而比較運行方式的合理性、可靠性和經(jīng)濟性。本文分析和比較了電力系統(tǒng)潮流計算的各種方法，其中PQ分解法具有計算速度快、所占內(nèi)存少等優(yōu)勢被廣泛采用。本設(shè)計首先進行了PQ分解法的理論研究，然后應(yīng)用 MATLAB軟件編制PQ分解法潮流程序，并通過一個電力系統(tǒng)實例對所編程序進行驗證，得出結(jié)論。
黨政相關(guān)相關(guān)推薦
鄂ICP備17016276號-1