電力網(wǎng)潮流電壓計(jì)算例題與matlab程序

1、1電力網(wǎng)潮流電壓計(jì)算例題與電力網(wǎng)潮流電壓計(jì)算例題與MATLAB程序程序編寫編寫佘名寰佘名寰本文介紹了兩個(gè)電力網(wǎng)潮流計(jì)算例題。一例為5個(gè)節(jié)點(diǎn)5條支路；另一例有6個(gè)節(jié)點(diǎn)7條支路，有PQ負(fù)荷節(jié)點(diǎn)也有PV發(fā)電機(jī)節(jié)點(diǎn)，變壓器支路標(biāo)么變比不為1。本文給出了完整的計(jì)算程序和計(jì)算結(jié)果。程序包括計(jì)算網(wǎng)絡(luò)節(jié)點(diǎn)導(dǎo)納矩陣和潮流電壓兩部分，程序中引用的M函數(shù)見參考文獻(xiàn)1。例題選自研究生教材，比較典型，計(jì)算結(jié)果與教材核對基本一致。本文可供電力專業(yè)師生和基層技術(shù)人員

2、潮流計(jì)算時(shí)參考?！纠纠?.3】如圖14所示，五個(gè)節(jié)點(diǎn)三條母線簡單系統(tǒng)，系統(tǒng)參數(shù)均用標(biāo)么值（pu.）表示。節(jié)點(diǎn)4、5通過變壓器與節(jié)點(diǎn)2、3相連，變比為1.05:1，節(jié)點(diǎn)1、2、3為PQ負(fù)荷節(jié)點(diǎn)，節(jié)點(diǎn)4為PV發(fā)電機(jī)節(jié)點(diǎn)，節(jié)點(diǎn)5為參考節(jié)點(diǎn)，節(jié)點(diǎn)4、5電壓均為1.05，圖中箭頭表示各支路的正方向，各個(gè)節(jié)點(diǎn)功率如表2.1所示計(jì)算各節(jié)點(diǎn)電壓、發(fā)電機(jī)無功功率和系統(tǒng)潮流分布：Z23Z13Z21YC2YC3YC1Z42Z531:1.051.05:142

3、315S3S2S1Z21=0.04J0.25Z23=0.08J0.30Z13=0.1J0.35Z42=J0.015Z53=J0.03YC1=J0.25YC2=J0.50YC3=J0.25S1=1.6J0.8S2=2J1S3=3.7J1.3P4=5.0V4=1.05V5=1.05P4=5.0V4=1.05vV5=1.05圖14[例1.3]網(wǎng)絡(luò)接線圖表11母線有功負(fù)載無功負(fù)載發(fā)電機(jī)有功發(fā)電機(jī)無功1PD1=1.6QD1=0.82PD2=2.0

4、QD2=1.03PD3=3.7QD3=1.345PG4=5.0PG5=？QG4=未定QG5=？3%yn0thenodeadmittanceofshuntcapacit%yn0i=[]%yn0=0yn0ij%ync=diag(conj(yn0))%y=yync[YYY]=fdt1(yyy)G=real(Y)B=imag(Y)%YYistheshundadmittancesofthelinetransfmersateachend%Yisno

5、deadmittancematrixfwk%u=sym([u1u2u3u4u5])delt=sym([d1d2d3d4d5])%unodevoltagemagnitudedeltanglep=[1.623.75]q=[0.81.01.3]k=0precision=1N1=4%theN1istheamountofthePQPVbusN2=1%theN2istheamountofthePVbusfm=1:N1fn=1:N11pt(n)=u(

6、m)u(n)(G(mn)cos(delt(m)delt(n))B(mn)sin(delt(m)delt(n)))endpp(m)=p(m)sum(pt)endfm=1:N1N2fn=1:N11qt(n)=u(m)u(n)(G(mn)sin(delt(m)delt(n))B(mn)cos(delt(m)delt(n)))endqq(m)=q(m)sum(qt)endJ1=jacobian(pp[d1d2d3d4u1u2u3])J2=jac

7、obian(qq[d1d2d3d4u1u2u3])J=vertcat(J1J2)uu=[1.01.01.01.051.05]dd=[00000]whileprecision0.00001u1=uu(1)u2=uu(2)u3=uu(3)u4=uu(4)u5=uu(5)d1=dd(1)d2=dd(2)d3=dd(3)d4=dd(4)d5=dd(5)fm=1:N1PP(m)=eval(pp(m))endfm=1:N1N2PP(N1m)=eva