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

用matlab实现电力系统潮流计算 晏鸣宇 PAGE PAGE 2 . . 电力系统三相短路计算 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,faul