电力网节点导纳矩阵计算例题与程序

1、电力网节点导纳矩阵计算例题与程序余名寰编写用计算机解算电力网潮流电压和短路电流问题首先需确定电力网的节点导纳矩阵或节点阻抗矩阵。本文通过例题介绍用网络拓扑法计算节点导纳矩阵的方法和程序，程序考虑了线路并联电容和变压器支路标么变比不为1时的影响。程序用MATLAB语言编写，线路参数均采用标么值。本文稿用。fficeword2007版编写，可供电气专业人员计算相关问题时参考。1 .用网络拓扑计算节点导纳矩阵1.1 网络拓扑矩阵：【例例图1-1是有5个节点和5条支路的网络，节点5作为基准参考点，1,2,3,4为独立节点，支路编号和方向图中已标识。Z21-0.044.J025Z23=0.08*J030

2、Z13=0.1+J035Z42=J0.015Z5340.03YQ=J025YC2=J0.50YC3=J0.25对于具有n个节点b条支路的有向图，它的关联矩阵为一个NXB的矩阵Aa：A.=aij若支路j与节点i相关，且笳头背离节点i,则1二1,若箭头指向节点则.二T，若支路j与节点i无关，则.二0,图1.1所示的有向图的关联矩阵为支路编号行编号从上到下为12345节点编号（5为参考节点）去掉第5行即为独立节点的关联矩阵。以下介绍生成网络关联矩阵的M函数义件ffm.m：%MFUNCTIONffm.m%NpisnumberofnodepointrNbisnumberofbraches%nstart-

3、thestartpointofbranches,r12ndtheendpointr%A-networkincidencematrixfunctionA=ffm(nstartrnend)globalNpNbn=iength(nstart);A=zeros(Np,Nb);fori=l:nA(nstart(i)ri)=1;A(nend(i)zi)=-l;end以例图1-1网络为例调用ffm.m文件求其关联矩阵运算以上程序可行关联矩阵midj如下:mm=10100110-100-1-10-10001000001Mmij明显与A巧是相同的。1.2 生成节点导纳矩阵程序：由网络原始矩阵计算节点导纳矩阵公式

4、Y=AYs0At(1-1)丫-节点导纳矩阵A-网络关联矩阵NA的转置矩阵Ys。-网络原始导纳矩阵000-000:000:000zbn.若网络各支路阻抗为Zb=zbl/zb2zbnZbio0zb2则ZSo=0000.,00Yso=Zso1(1-2)Y=AZso“A，(13)节点导纳计算程序以例1-1网络为例，在不计对地电容和变压器变比假定为1条件卜.，节点导纳矩阵计算程序如h：globalNpNb%Npisnumberofnodepoint,NbisnumberofbrachesrNp=5;Nb=5;%nstart-thestartpointofbranches,nend一一theendpoin

5、t,%mm-network,incidencematrixnstart=2,2rlr4,5;nend=l,3,3,2,3;nnn=ffm(nstartrnend);%zbl.theseriesimpedancesoftransmissionline%ybl,theseriesadmittancesoftransmissionlinezblr=0.04/0.08/0.1/0.0,0.0;zbli=0.25,030,035,0.015,0.03;zbl=zblr+zbli*j;ybl=zbl.A(-l);yb=diag(conj(ybl);y=mm*yb*(mm),;运算程序后可得节点导纳矩阵13

6、787-6.54171-0.6240+3.9002i-0.7547+2.6415i00-0.6240+3.900211.4539-73.6789i-0.8299+3.1120i0+66.6667i0-0.7547+2.6415i-0.8299+3.1120i1.5846-39.0869i00+33.3333i0+66.6667i0-66.6667i0+333333103333331线路对地并联电容的处理例图1-1网络中，线路，两端对地电容导纳均为j0.25,节点1,2,3自导纳相应增加j0.25,j0.5,j0.25o程序处理如下:%ybO,theshuntadmittancesoftrans

7、missionlineyb0i=0.25,0.25/0.0/0.0,0.0;ybO=O+ybOi*j;yO=diag(conj(ybO)；yg=mm*y0t(mm),;yn=diag(diag(yg);yy=yn-yg;y=y+yn;程序中ybO为各线路对地电容导纳，y0为原始导纳矩阵，yn对角矩阵对应各节点白导纳增加值，yy相应各支路两端对地并联导纳，用于计算支路电流。y为计及对地电容后的节点导纳矩阵。程序运算结果yn=0+0.2500i000000+0.5000i000000+0.2500i000000000000yy=00+0.2500i0000+0.2500i00+0.2500i000

8、+0.2500i0000000000000013787-6.2917i-0.6240+3.9002i-0.7547+2.6415i00-0.6240+3.9002i1.4539-73.1789i-0.8299+3.1120i0+66.6667i0-0.7547+2.6415i-0.8299+3.112011.584638.8369i00+333333i000+66.66671000+33.3333i0-66.6667i00033.3333i非标准变压器支路的处理当支路仃变压器并且标么变比不为1时，可将变压器化作n形等值电路，如图厂2所示。图中理想变压器变比为k:l,接入节点i,高压侧经漏抗Zi

9、与节点j相连。CK(K-1)Yt图1-2变压器等值电路等值电路图中Yt=1/ZToi侧自导变为y,.=kyt+k(k-l)yt=lj=-kytAy”二-kyt+yt=-(kT)y,。j侧自导变为yjj=kyt+(l*M=yt/与K=1时相同；变压器变比不为1时修正导纳矩阵的M函数文件fdtl.m%MFUNCTIONfdtl.m%thetransformerturnsratio,withoff-nominaltap-setting,fig1-2%bttwo-dimensionalarray,lineoneispointirlinet,woispointj%pisthetransformertur

10、nsratioFkisnumberoftransformersfunctionA,B=fdt1(a,b)globalbtpmkfori=l:mforj=l:mA(i,j)=a(i,j);B(i,j)=b(i,j);endendforil=l:kkl=bt(l.il);kZ=br2,11);t=-a(klrk2);pl=p(il)；fori=l:mforj=l:mifi=klifj=klA(i,j)=A(i,j)+(pl*pl-l)*t;elseifj=k2A(i,j)=A(i,B(i,j)=pl*(pl-1)*t;endelseifi=k2ifj=klA(i,j)=A(i,j)-(pl-1)*

11、t;B(i,j)=(1-pl)*t;endendendendend程序中bt是2行K列数组，第一行为变压器接入点i的节点编号，第二行为经漏抗后所接j节点编号，k为变压器支路数。P为一维k列数组，输入各个变压器的变比；m为矩阵a,b,A,B的阶数a, A是进行变比修正前后的节点导纳矩阵，b, B是变比修正前后支路对地并联导纳矩阵。13例图节点导纳和支路对地并联导纳计算程序例图1-1共有5个节点，5条支路，所以Np=5,Nb=5支路，为变压器支路，其i节点编号为2,j节点编号为4,变比为1：1/1.05：支路也为变压器支路，其i节点编号为3,j节点编号为5,变比为1：1/1.05：所以数组bt=2

12、,3;4,5,p=l/1.05,l/1.05bk=2;5条支路按箭头方向首末端节点分别为221,4,5和13323,所以nstart=2,24,4,5;nend=13323;输入以上数据可得以卜程序:%ThefollowingProgramfornodeimpedancescalculationisbasedonMATLAB2007clear%bus1,2,3isPQbus,bus4isPVbus.bus5isslackbus%*NUll.m,example1-1,fig1-1*globalNpNbbtpmk%Npisnumberofnodepoint,Nbisnumberofbraches,

13、%pthetransformerturnsratio,withoff-nominaltap-setting,fig1-1%btistwo-dimensionalarrayJineoneispointijinetwoispointj,%misnumberofnodepoint,kisnumberoftransformersNp=b;Nb=bt=23；4,5;p=l/1.05/l/1.05;m=5;k=2;%nstart-thestartpointofbranches,nend-theendpoint,%mm-networkincidencematrixnsta11=2,24,4,5;nend=1

14、3/3/23；mm=ffm(nstart/nend);%zbLtheseriesimpedancesoftransmissionline%ybXtheseriesadmittancesoftransmissionlinezblr=0.04z0.08/0.1/0.0/0.0;zbli，025,030,035,0.015,0.03;zbl=zblr+zbli*j;ybl=zbl.A(-l);yb=diag(conj(ybl);y=mm*yb*(mm),;%ybOztheshuntadmittancesoftransmissionlineyb0i=0.25,0.25,0.0,0.0,0.0;ybO=

15、O+ybOi*j;yO=diag(conj(ybO*)；yg=mm*yO*(mm);yn=diag(diag(yg);yy=yn-yg;y=y+yn;Y,YY=fdtl(y,yy);%YYistheshundadmittancesofthelineandtransformersateachend,%Yisnode-admittancematrixfornetwork运行结果为：Y=13787-6.29171-0.6240+3.9002i-0.7547+2.6415i00-0.6240+3.9002i1.4539-66.9808i-0.8299+3.1120i0+63.4921i0-0.7547+2.6415i-0.8299+3.1120i1.5846-35.7379i00+31.7460i00+63.4921100-66