《三相短路短路电流计算matlab程序文档.pdf》由会员分享，可在线阅读，更多相关《三相短路短路电流计算matlab程序文档.pdf（25页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、 电力系统三相短路计算 a main.m 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 是考虑了变比，且支路未近似的导纳矩阵；Ymatrix2 是近似变比为 1，但是支
2、路未近似计算的节点导纳矩阵；Ymatrix3 是近似变 比为 1，采取近似支路参数 1 的导纳矩阵；Ymatrix4 是近似变比为 1，采取近似支路参 数 2 的导纳矩阵。Y=Ymatrix2(bus,line);%对故障点进行导纳修正 fixY=FixY(Y,bus,fault);%求注入电流 Iinj=Inode(bus,calcSettings);%得到故障支路与其他支路电流 Bcurrent=Ibranch(line,U,fault,Y);%如果发生支路三相短路，那么对应该支路的电流修正为-999999-j999999 Ib=ReviseBcurrent(fault,Bcurrent)
3、;%结束计时 tim=toc;fprintf(程序运行结果 );fprintf(计算完成，共用时%4.4fs,相关结果已保存在 output.dat,tim);%输出结果 fprintf_result(output_file,Ib);fprintf_result1(Ib);b FixY.m function fixY=FixY(Y,bus,fault)%对形成的导纳矩阵进行故障点的修正 nb,mb=size(bus);nf,mf=size(fault);fixY=Y;%对发电机节点导纳修正 for k=1:nb busType=bus(k,7);if(busType=1)fixY(bus(k,
4、1),bus(k,1)=fixY(bus(k,1),bus(k,1)+1/1i/bus(k,8);end end%对节点短路和支路短路的导纳矩阵进行修正 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 if(dis=0)&(nodeI*nodeJ=0)fixY(nodeI,nodeI)=99
5、9999+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(nodeI,nodeJ)/(1-dis);fixY(nodeI,nodeJ)=0;fixY(nodeJ,nodeI)=0;end en
6、d end c fprintf_result.m function output_args =fprintf_result(output_file,Ib)%将得到的短路电流输入到输出文件中 n,m=size(Ib);fprintf(output_file,No.No.vector of I value of I );for k=1:n I=Ib(k,1);J=Ib(k,2);I01=real(Ib(k,3);I02=imag(Ib(k,3);I1=Ib(k,4);if(I02=0)fprintf(output_file,%3d%3d%10.6f+j%10.6f%10.6f,I,J,I01,I0
7、2,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*K*Yt;Y(I,J)=Y(I,J)+K*Yt;Y(J,I)=Y(I,J);end end end i Ymatrix2.m function Y=Ymatrix2(bus,line)%考虑变
8、压器变比近似为 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(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;end end end j Ymatrix3.m function
9、 Y=Ymatrix3(bus,line)%考虑变压器变比为 1，采用支路参数近似 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=imag(1/Zt);Ym=imag(line(k,5)+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)
10、+Ym;end end end k Ymatrix4.m function 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;end end end