牛頓拉夫遜法潮流計算matlab程序

《牛頓拉夫遜法潮流計算matlab程序》由會員上傳分享，免費在線閱讀，更多相關內容在行業(yè)資料-天天文庫。

1、%電力系統(tǒng)的潮流計算,以下程序參考文獻《電力系統(tǒng)畢業(yè)設計》中國水利電力出版社%(該文獻用極坐標下的牛頓——拉夫遜方法實現(xiàn)，在此為了與課本一致做了修改)%為了計算方便將原來的下標做以下修改:S2S3S4S5U2U3U4U5改為S1S2S3S4U1U2%U3U4,即原題的平衡點1就變?yōu)楝F(xiàn)在的平衡點5%1.形成節(jié)點導納矩陣,yb55=6.250-18.750j;yb51=-5.000+15.000j;yb52=-1.250+3.750j;yb53=0.000-0.000j;yb54=0.000-0.000j;yb15=-5.000+15.000j;yb11=10.834-32.500j;yb12=

2、-1.667+5.000j;yb13=-1.667+5.000j;yb14=-2.500+7.500j;yb25=-1.250+3.750j;yb21=-1.667+5.000j;yb22=12.917-38.750j;yb23=-10.000+30.000j;yb24=0.000-0.000j;yb35=0.000-0.000j;yb31=-1.667+5.000j;yb32=-10.000+30.000j;yb33=12.917-38.750j;yb34=-1.250+3.750j;yb45=0.000-0.000j;yb41=-2.500+7.500j;yb42=0.000-0.000

3、j;yb43=-1.250+3.750j;yb44=3.750-11.250j;YB=[yb11yb12yb13yb14yb15;yb21yb22yb23yb24yb25;yb31yb32yb33yb34yb35;yb41yb42yb43yb44yb45;yb51yb52yb53yb54yb55];%計算各節(jié)點功率的不平衡量設U=E+jF；Y=G+Bj；E(1)=1.00;E(2)=1.00;E(3)=1.00;E(4)=1.00;F(1)=0;F(2)=0;F(3)=0;F(4)=0;G=real(YB);B=imag(YB);%設S=P+Bj；S(1)=0.20+0.20i;S(2)=-

4、0.45-0.15i;S(3)=-0.40-0.05i;S(4)=-0.60-0.10i;P=real(S);Q=imag(S);k=0;precision=1;N1=4;whileprecision>0.00000001E(5)=1.06;F(5)=0;form=1:N1forn=1:N1+1%計算Pi,Qi,設Pi=Pt;Qi=QtPt(n)=(E(m)*(G(m,n)*E(n)-B(m,n)*F(n))+F(m)*(G(m,n)*F(n)+B(m,n)*E(n)));Qt(n)=(F(m)*(G(m,n)*E(n)-B(m,n)*F(n))-E(m)*(G(m,n)*F(n)+B(m,

5、n)*E(n)));end%設P,Q的改變量為dP,dQdP(m)=P(m)-sum(Pt);dQ(m)=Q(m)-sum(Qt);endform=1:N1forn=1:N1+1%計算HijNijJijLijH(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);J(m,n)=-B(m,n)*F(m)-G(m,n)*E(m);L(m,n)=G(m,n)*F(m)-B(m,n)*E(m);endendform=1:N1forn=1:N1+1Bi(n)=G(m,n)*F(n)+B(m,n)*E(n);Ai(n)=G(m,n)*

6、E(n)-B(m,n)*F(n);end%計算Hii,Nii,Jii,Lii,由公式4-44b左側公式實現(xiàn),sum(Ai)，sum(Bi)用于實現(xiàn)公式中的sigerma從j到n的求和;H(m,m)=sum(Bi)-(B(m,m)*E(m)+G(m,m)*F(m))+2*G(m,m)*F(m);N(m,m)=sum(Ai)-(G(m,m)*E(m)-B(m,m)*F(m))+2*G(m,m)*E(m);J(m,m)=-2*B(m,m)*F(m)+sum(Ai)-(G(m,m)*E(m)-B(m,m)*F(m));L(m,m)=-2*B(m,m)*E(m)-(sum(Bi)-(B(m,m)*E(

7、m)+G(m,m)*F(m)));end%設雅可比矩陣為JJ，以下語句用來實現(xiàn)雅可比矩陣中對角線上元素HNJL的排列form=1:N1JJ(2*m-1,2*m-1)=H(m,m);JJ(2*m-1,2*m)=N(m,m);JJ(2*m,2*m-1)=J(m,m);JJ(2*m,2*m)=L(m,m);end%以下語句用于實現(xiàn)雅可比矩陣非對角線上元素的排列form=1:N1forn=1:N1ifm==nelse