22节点初步短路计算

clear all;
fn=‘22节点.txt’; %数据文件

sys=DataRead(fn); %读取数据文件
Y1=sys.Y1; %正序
Y2=sys.Y2; %负序
Y0=sys.Y0; %零序
Vf0=1; %故障前短路点电压
pos=11; %短路点位置
Z1=full(inv(Y1)); %结点导纳矩阵：正序
Z2=full(inv(Y2)); %负序
Z0=full(inv(Y0)); %零序
Z1n=Z1(pos,pos); %短路点戴维宁阻抗：正序，负序，零序
Z2n=Z2(pos,pos);
Z0n=Z0(pos,pos);
%% 单相接地短路
Ifn=Vf0/(Z1n+Z2n+Z0n); %故障点故障电流：正序，负序，零序
Ifn1=Ifn;
Ifn2=Ifn;
Ifn0=Ifn;

Vfn1=Vf0-Z1nIfn1; %故障点故障电压：正序，负序，零序
Vfn2=-Z2nIfn2;
Vfn0=-Z0nIfn0;
V1=1+Z1(:,pos)(-Ifn1); %各节点电压：正序，负序，零序
V2=Z2(:,pos)(-Ifn2);
V0=Z0(:,pos)(-Ifn0);
a=exp(1j2pi/3);
T=[1 1 1;a^2 a 1;a a^2 1];
Vabc=[T*[V1,V2,V0]’]’; %结点电压

%支路电流
sys=dlmread(fn);
flag=find(sys(:,1)==0);
N=sys(1,1);
Branch=sys(flag(1)+1:flag(2)-1,1:6);

%计算各支路电流：正序，负序，零序
for K=1:size(Branch,1)
I=Branch(K,3);
J=Branch(K,4);
x1=Branch(K,5);
x0=Branch(K,6);
IS1=(V1(I)-V1(J))/(Z1(I,J));
IS2=(V2(I)-V2(J))/(Z2(I,J));
if x0~=0
IS0=(V0(I)-V0(J))/(1j*x0);
else
IS0=0;

end
BSC(K,:)=full([K,I,J,IS1,IS2,IS0]);

end

%输出结果
fprintf(’ 单 相 接 地 短 路 \n\n’);
fprintf(’=\n’);
fprintf(’ 接 地 点 故 障 分 量 |\n’);
fprintf(’=\n’);
fprintf(’ 正序电流 | 负序电流 | 零序电流 |\n’);
fprintf(’=\n’);
fprintf(’ % 4.4f+j% 4.4f % 4.4f+j%4.4f % 4.4f+j% 4.4f \n’,…
real(Ifn1),imag(Ifn1),real(Ifn2),imag(Ifn2),real(Ifn0),imag(Ifn0));
fprintf(’=\n’);
fprintf(’ 正序电压 | 负序电压 | 零序电压 |\n’);
fprintf(’=\n’);
fprintf(’ % 4.4f+j% 4.4f % 4.4f+j% 4.4f % 4.4f+j% 4.4f \n’,…
real(Vfn1),imag(Vfn1),real(Vfn2),imag(Vfn2),real(Vfn0),imag(Vfn0));
fprintf(’=\n’);
fprintf(’\n\n\n’);
fprintf(’ 各 支 路 电 流 \n\n’);
fprintf(’=====================================================================================================\n’);
fprintf(’ 序号 | From BUS | To BUS | 正 序 | 负 序 | 零 序 |\n’);
fprintf(’----------