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

电力系统三相短路计算 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，但是支路未近似计算的节点导纳矩阵；Ymatrix3是近似变比为1，采取近似支路参数1的导纳矩阵；Ymatrix4是近似变比为1，采取近似支路参数2的导纳矩阵。 Y = Ymatrix2(bus,line); %对故障点进行导纳修正 fixY = FixY(Y,bus,fault); %求注入电流 Iinj = Inode(bus,calcSettings); U = fixY\Iinj; %得到故障支路与其他支路电流 Bcurrent = Ibranch( line,U,fault,Y ); %如果发生支路三相短路，那么对应该支路的电流修正为-999999-j999999 Ib = ReviseBcurrent( fault,Bcurrent ); %结束计时 tim=toc; fprintf(\n程序运行结果); fprintf(\n计算完成，共用时%4.4fs,相关结果已保存在output.dat\n,tim); %输出结果 fprintf_result(output_file, Ib); fprintf_result1(Ib); 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,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) = 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(nodeI,nodeJ)/(1-dis); fixY(nodeI,nodeJ)=0; fixY(nodeJ,nodeI)=0; end end end fprintf_result.m function [ output_args ] = fprintf_result( output_file, Ib ) %将得到的短路电流输入到输出文件中 [n,m]=size(Ib); fprintf( output_file, No. No.