matlab三相短路电流计算程序_电力系统三相短路计算的MATLAB代码.doc

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

电力系统三相短路计算

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 );

%如果发生支路三相短路，那