三相短路短路电流计算(matlab程序)

%% 三门水电站三相短路程序
%% 电气一班 刘永康 0805010114
%% 短路数据读入部分（读入数据为自定应格式）
clear
clc
format long
fid=fopen('C:\Users\john\Desktop\三门水电站程序\三门水电站设计短路部分\短路读入数据.txt','r');
tline=fgetl(fid);
node1=[];node2=[];node3=[];R=[];X1=[]; X2=[]; %初始化，或者后面加个变量a，用node1（a）,就不用初始定义了
while tline(1)~='L'
tline=fgetl(fid);
end
while tline(1)=='L'
node1=[node1 eval(tline(5))]; %平衡结点为1，PV结点为2，PQ结点为3
node2=[node2 eval(tline(13))];
R=[R eval(tline(19))];
X1=[X1 eval(tline(25:30))];
tline=fgetl(fid);
end
while tline(1)~='H'
tline=fgetl(fid);
end
while tline(1)=='H'
node3=[node3 eval(tline(5))];
X2=[X2 eval(tline(12:18))];
tline=fgetl(fid);
end
fclose(fid);
n=max(max(node1),max(node2)); %求取结点个数
%% 形成导纳矩阵（初始化）
Y=zeros(n,n);
[m1,m2]=size(node1);
for aa=1:m2
Y(node1(aa),node2(aa))=Y(node1(aa),node2(aa))-1/(R(aa)+j*X1(aa));
Y(node2(aa),node1(aa))=Y(node1(aa),node2(aa));
Y(node1(aa),node1(aa))=Y(node1(aa),node1(aa))+1/(R(aa)+j*X1(aa));
Y(node2(aa),node2(aa))=Y(node2(aa),node2(aa))+1/(R(aa)+j*X1(aa));
end
[m3,m4]=size(node3);
for aa=1:m4
Y(node3(aa),node3(aa))= Y(node3(aa),node3(aa))+1/(j*X2(aa));
end
%% 计算三相短路电流并输出
%% 有潮流计算结果可知：高压侧电压为1.0118，中压侧电压为1.0610，低压侧电压为1.0664
Z=inv(Y);
disp('高压侧：');
fprintf('程序计算的短路电流I": %10.4fKA ， 手算结果： 6.0079KA\n',abs(1.0118/Z(3,3))*100/sqrt(3)/230);
fprintf('程序计算的冲击电流Im: %10.4fKA ， 手算结果： 15.7184KA \n\n',1.85*sqrt(2)*abs(1.0118/Z(3,3))*100/sqrt(3)/230);
disp('中压侧：');
fprintf('程序计算的短路电流I": %10.4fKA ， 手算结果： 7.9154KA\n',abs(1.0610/Z(6,6))*100/sqrt(3)/115);
fprintf('程序计算的冲击电流Im: %10.4fKA , 手算结果： 20.709KA \n\n',1.85*sqrt(2)*abs(1.0610/Z(6,6))*100/sqrt(3)/115);
disp('低压侧：');
fprintf('程序计算的短路电流I": %10.4fKA ， 手算结果： 86.2745KA\n',abs(1.0664/Z(4,4))*100/sqrt(3)/10.5);
fprintf('程序计算的冲击电流Im: %10.4fKA , 手算结果： 225.7195KA \n\n',1.85*sqrt(2)*abs(1.0664/Z(4,4))*100/sqrt(3)/10.5);




下面是程序读入的txt文件，自己新建一个即可(格式模仿bpa格式)

1234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890
A node1 node2 R X1
L 1 2 0 0.0328

L 2 3 0 0.0189
L 3 4 0 0.1167
L 3 5 0 0.0583
L 4 6 0 0.075
L 5 6 0 0.0375
1234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890
B node3 X2
H 1 0.27758
H 1 0.4496
H 2 0.0756
H 4 0.20825
H 5 0.1041