%%单相接地短路计算
function [U_single_A,I_single_A]=single_PSC(f,k,node,line_num,Z1,Z2,Z0,line_parameter)
%故障点电流
If1=1/(Z1(f,f)+Z2(f,f)+Z0(f,f)+3*k);
If2=If1;
If0=If1;
%各节点短路电压
U_single1=1-Z1(:,f)*If1;
U_single2=-Z2(:,f)*If2;
U_single0=-Z0(:,f)*If0;
%各支路的短路电流
I_single1=sparse(node,node);
I_single2=sparse(node,node);
I_single0=sparse(node,node);
for n=1:line_num
I_single1=I_single1+sparse(line_parameter(n,2),line_parameter(n,3),(U_single1(line_parameter(n,2))-U_single1(line_parameter(n,3)))/Z1(line_parameter(n,2),line_parameter(n,3)),node,node);
I_single2=I_single2+sparse(line_parameter(n,2),line_parameter(n,3),(U_single2(line_parameter(n,2))-U_single2(line_parameter(n,3)))/Z2(line_parameter(n,2),line_parameter(n,3)),node,node);
I_single0=I_single0+sparse(line_parameter(n,2),line_parameter(n,3),(U_single0(line_parameter(n,2))-U_single0(line_parameter(n,3)))/Z0(line_parameter(n,2),line_parameter(n,3)),node,node);
end
%计算各相数据
%故障点短路电流
If_A=If1+If2+If0;
If_B=If1*(cos(pi/3)+1i*sin(pi/3))^2+If2*(cos(pi/3)+1i*sin(pi/3))+If0;
If_C=If1*(cos(pi/3)+1i*sin(pi/3))+If2*(cos(pi/3)+1i*sin(pi/3))^2+If0;
%节点短路电压
U_single_A=U_single1+U_single2+U_single0;
U_single_B=U_single1*(cos(pi/3)+1i*sin(pi/3))^2+U_single2*(cos(pi/3)+1i*sin(pi/3))+U_single0;
U_single_C=U_single1*(cos(pi/3)+1i*sin(pi/3))+U_single2*(cos(pi/3)+1i*sin(pi/3))^2+U_single0;
%支路短路电流
I_single_A=I_single1+I_single2+I_single0;
I_single_B=I_single1*(cos(pi/3)+1i*sin(pi/3))^2+I_single2*(cos(pi/3)+1i*sin(pi/3))+I_single0;
I_single_C=I_single1*(cos(pi/3)+1i*sin(pi/3))+I_single2*(cos(pi/3)+1i*sin(pi/3))^2+I_single0;
%%输出参数
disp('单相接地短路');
disp('故障点A相电流：');
disp(If_A);
disp('故障点B相电流：');
disp(If_B);
disp('故障点C相电流：');
disp(If_C);
disp('各节点的A相电压：');
disp(U_single_A);
disp('各节点的B相电压：');
disp(U_single_B);
disp('各节点的C相电压：');
disp(U_single_C);
disp('各支路的A相电流：');
disp(I_single_A);
disp('各支路的B相电流：');
disp(I_single_B);
disp('各支路的C相电流：');
disp(I_single_C);
end
