基于matlab的兩機五節(jié)點網(wǎng)絡(luò)潮流仿真計算—牛拉法計劃書(編輯修改稿)
2025-06-09 23:41 本頁面
【文章內(nèi)容簡介】 路功率分布，PV節(jié)點無功功率和平衡節(jié)點柱入功率。以直角坐標(biāo)系形式表示①. 迭代推算式 采用直角坐標(biāo)時,節(jié)點電壓相量及復(fù)數(shù)導(dǎo)納可表示為: (310)將以上二關(guān)系式代入上式中,展開并分開實部和虛部。假定系統(tǒng)中的第1,2,m號為P—Q節(jié)點,第m+1,m+2,n1為P—V節(jié)點,根據(jù)節(jié)點性質(zhì)的不同,得到如下迭代推算式: ⑴對于PQ節(jié)點 (311)⑵對于PV節(jié)點 (312)⑶對于平衡節(jié)點 平衡節(jié)點只設(shè)一個,電壓為已知,不參見迭代,其電壓為: (313)②. 修正方程選定電壓初值及變量修正量符號之后代入式中,并將其按泰勒級數(shù)展開,略去二次方程及以后各項,得到修正方程如下: (314) ③.雅可比矩陣各元素的算式式(37)中, 雅可比矩陣中的各元素可通過對式(311)和(312), 雅可比矩陣中非對角元素為 (315)當(dāng)時,雅可比矩陣中對角元素為: (316) 牛頓—拉夫遜法的程序框圖 第四章 牛頓—拉夫遜法潮流具體計算 牛頓—拉夫遜直角坐標(biāo)潮流計算Matlab程序及運行結(jié)果、Matlab程序clccleardisp(39。節(jié)點總數(shù)為：39。)。N=5disp(39。平衡節(jié)點為：39。)。1disp(39。PQ節(jié)點為：39。)。JD=[2,3,4,5]e=[ 1 1 1 1]。f=[0 0 0 0 0]。P1=0。Q1=0。P2=。Q2=。P3=。Q3=。P4=。Q4=。P5=。Q5=。G=[,0,0。,。,0。0,。0,0,]。%形成電導(dǎo)矩陣。B=[,0,0。,。,0。0,。0,0,]。%形成電納矩陣。disp(39。節(jié)點電導(dǎo)矩陣G為:39。)。 disp(G)disp(39。節(jié)點電納矩陣B為:39。)。disp(B)k=0。for v=1:7I=[0,0。0,0。0,0。0,0。0,0]。for n=1:5 I(1,1)=I(1,1)+G(1,n)*e(n)B(1,n)*f(n)。 I(1,2)=I(1,2)+G(1,n)*f(n)+B(1,n)*e(n)。endfor n=1:5 I(2,1)=I(2,1)+G(2,n)*e(n)B(2,n)*f(n)。 I(2,2)=I(2,2)+G(2,n)*f(n)+B(2,n)*e(n)。endfor n=1:5 I(3,1)=I(3,1)+G(3,n)*e(n)B(3,n)*f(n)。 I(3,2)=I(3,2)+G(3,n)*f(n)+B(3,n)*e(n)。endfor n=1:5 I(4,1)=I(4,1)+G(4,n)*e(n)B(4,n)*f(n)。 I(4,2)=I(4,2)+G(4,n)*f(n)+B(4,n)*e(n)。endfor n=1:5 I(5,1)=I(5,1)+G(5,n)*e(n)B(5,n)*f(n)。 I(5,2)=I(5,2)+G(5,n)*f(n)+B(5,n)*e(n)。 endH=[]。N=[]。M=[]。L=[]。J=[]。P2=P2e(2)*I(2,1)f(2)*I(2,2)。 %有功功率的不平衡量Q2=Q2f(2)*I(2,1)+e(2)*I(2,2)。 %無功功率的不平衡量P3=P3e(3)*I(3,1)f(3)*I(3,2)。Q3=Q3f(3)*I(3,1)+e(3)*I(3,2)。P4=P4e(4)*I(4,1)f(4)*I(4,2)。Q4=Q4f(4)*I(4,1)+e(4)*I(4,2)。P5=P5e(5)*I(5,1)f(5)*I(5,2)。Q5=Q5f(5)*I(5,1)+e(5)*I(5,2)。for m=2:5 for n=2:5 if(m==n) H(m,m)=B(m,m)*e(m)+G(m,m)*f(m)+I(m,2)。 N(m,m)=G(m,m)*e(m)+B(m,m)*f(m)+I(m,1)。 M(m,m)=G(m,m)*e(m)B(m,m)*f(m)+I(m,1)。 L(m,m)=B(m,m)*e(m)+G(m,m)*f(m)I(m,2)。 else H(m,n)=B(m,n)*e(m)+G(m,n)*f(m)。 N(m,n)=G(m,n)*e(m)+B(m,n)*f(m)。 M(m,n)=N(m,n)。 L(m,n)=H(m,n)。 end endendJ=[H(2,2),N(2,2),H(2,3),N(2,3),H(2,4),N(2,4),H(2,5),N(2,5)。M(2,2),L(2,2),M(2,3),L(2,3),M(2,4),L(2,4),M(2,5),L(2,5)。H(3,2),N(3,2),H(3,3),N(3,3),H(3,4),N(3,4),H(3,5),N(3,5)。M(3,2),L(3,2),M(3,3),L(3,3),M(3,4),L(3,4),M(3,5),L(3,5)。H(4,2),N(4,2),H(4,3),N(4,3),H(4,4),N(4,4),H(4,5),N(4,5)。M(4,2),L(4,2),M(4,3),L(4,3),M(4,4),L(4,4),M(4,5),L(4,5)。H(5,2),N(5,2),H(5,3),N(5,3),H(5,4),N(5,4),H(5,5),N(5,5)。M(5,2),L(5,2),M(5,3),L(5,3),M(5,4),L(5,4),M(5,5),L(5,5)]。disp(39。雅克比矩陣J:39。)。 disp(J)。A=[]。C=[P2。Q2。P3。Q3。P4。Q4。P5。Q5]A=J\C。%解修正方程式disp(39。第M次修正方程的解A:39。)。disp(A)。 f(2)=f(2) +A(1,1)。e(2)=e(2) +A(2,1)。 %計算新值f(3)=f(3) +A(3,1)。
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1