电力系统三相短路计算的MATLAB代码(10页).doc

《电力系统三相短路计算的MATLAB代码(10页).doc》由会员分享，可在线阅读，更多相关《电力系统三相短路计算的MATLAB代码(10页).doc（10页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、-电力系统三相短路计算的MATLAB代码-第 10 页电力系统三相短路计算clear tim%打开文件dfile,pathname=uigetfile(*.m,Select Data File);if pathname = 0 error( you must select a valid data file)else lfile =length(dfile); eval(dfile(1:lfile-2);end%定义输出文件output_file=fopen(output.dat,w);%开始计时tic; %求解节点导纳矩阵，其中Ymatrix1是考虑了变比，且支路未近似的导纳矩阵；Ymatr

2、ix2是近似变比为1，但是支路未近似计算的节点导纳矩阵；Ymatrix3是近似变比为1，采取近似支路参数1的导纳矩阵；Ymatrix4是近似变比为1，采取近似支路参数2的导纳矩阵。Y = Ymatrix2(bus,line);%对故障点进行导纳修正fixY = FixY(Y,bus,fault);%求注入电流Iinj = Inode(bus,calcSettings);U = fixYIinj;%得到故障支路与其他支路电流Bcurrent = Ibranch( line,U,fault,Y );%如果发生支路三相短路，那么对应该支路的电流修正为-999999-j999999Ib = Revis

3、eBcurrent( fault,Bcurrent );%结束计时tim=toc;fprintf(n程序运行结果);fprintf(n计算完成，共用时%4.4fs,相关结果已保存在output.datn,tim);%输出结果fprintf_result(output_file, Ib);fprintf_result1(Ib);function fixY = FixY( Y,bus,fault )%对形成的导纳矩阵进行故障点的修正nb,mb=size(bus);nf,mf= size(fault);fixY = Y;%对发电机节点导纳修正for k=1:nb busType=bus(k,7);

4、if (busType=1) fixY(bus(k,1),bus(k,1) = fixY(bus(k,1),bus(k,1) + 1/1i/bus(k,8); endend%对节点短路和支路短路的导纳矩阵进行修正for k=1:nf nodeI=fault(k,1); nodeJ=fault(k,2); dis=fault(k,3); if (nodeI=0) fixY(nodeJ,nodeJ) = 999999+1i*999999; continue; end if (nodeJ=0) fixY(nodeI,nodeI) = 999999+1i*999999; continue; end i

5、f (dis=0)&(nodeI*nodeJ=0) fixY(nodeI,nodeI) = 999999+1i*999999; continue; end if (dis=1)&(nodeI*nodeJ=0) fixY(nodeJ,nodeJ) = 999999+1i*999999; continue; end if (dis=1)&(dis=0)&(nodeI*nodeJ=0) fixY(nodeI,nodeI) = fixY(nodeI,nodeI) - fixY(nodeI,nodeJ)/dis; fixY(nodeJ,nodeJ) = fixY(nodeJ,nodeJ) - fixY(

6、nodeI,nodeJ)/(1-dis); fixY(nodeI,nodeJ)=0; fixY(nodeJ,nodeI)=0; endendenda fprintf_result.mfunction output_args = fprintf_result( output_file, Ib )%将得到的短路电流输入到输出文件中n,m=size(Ib); fprintf( output_file, No. No. vector of I value of In);for k=1:n I=Ib(k,1); J=Ib(k,2); I01=real(Ib(k,3); I02=imag(Ib(k,3);

7、 I1=Ib(k,4); if(I02=0) fprintf( output_file, %3d %3d %10.6f+j%10.6f %10.6f,I,J,I01,I02,I1); end if(I02=0) fprintf(%3d %3d %10.6f+j%10.6f %10.6f,I,J,I01,I02,I1); end if(I020 Y(I,I)=Y(I,I)+Yt+Ym; Y(J,J)=Y(J,J)+Yt/(K*K); Y(I,J)=Y(I,J)-Yt/K; Y(J,I)=Y(I,J); end if K0 Y(I,I)=Y(I,I)+Yt+Ym; Y(J,J)=Y(J,J)+K*

8、K*Yt; Y(I,J)=Y(I,J)+K*Yt; Y(J,I)=Y(I,J); endendendg Ymatrix2.mfunction Y = Ymatrix2( bus,line )%考虑变压器变比近似为1，支路参数不等效nb,mb=size(bus);nl,ml=size(line);Y=zeros(nb,nb);for k=1:nl I=line(k,1); J=line(k,2); Zt=line(k,3)+1i*line(k,4); Yt=1/Zt; Ym=line(k,5)+1i*line(k,6); if J=0 Y(I,I)=Y(I,I)+Yt+Ym; Y(J,J)=Y(

9、J,J)+Yt+Ym; Y(I,J)=Y(I,J)-Yt; Y(J,I)=Y(I,J); end if J=0 Y(I,I)=Y(I,I)+Ym; endendendh Ymatrix3.mfunction Y = Ymatrix3( bus,line )%考虑变压器变比为1，采用支路参数近似1nb,mb=size(bus);nl,ml=size(line);Y=zeros(nb,nb);for k=1:nl I=line(k,1); J=line(k,2); Zt=line(k,3)+1i*line(k,4); Yt=imag(1/Zt); Ym=imag(line(k,5)+1i*line

10、(k,6); if J=0 Y(I,I)=Y(I,I)+Yt+Ym; Y(J,J)=Y(J,J)+Yt+Ym; Y(I,J)=Y(I,J)-Yt; Y(J,I)=Y(I,J); end if J=0 Y(I,I)=Y(I,I)+Ym; endendendi Ymatrix4.mfunction Y = Ymatrix4( bus,line )%变压器变比近似为1，采用支路等效参数2 nb,mb=size(bus);nl,ml=size(line);Y=zeros(nb,nb);for k=1:nl I=line(k,1); J=line(k,2); Zt=1i*line(k,4); Yt=1/Zt; Ym=1i*line(k,6); if J=0 Y(I,I)=Y(I,I)+Yt+Ym; Y(J,J)=Y(J,J)+Yt+Ym; Y(I,J)=Y(I,J)-Yt; Y(J,I)=Y(I,J); end if J=0 Y(I,I)=Y(I,I)+Ym; endendend