电力系统三相短路计算地MATLAB代码.docxVIP

实用文档 电力系统三相短路计算 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-