matlab编写的短路计算程序

%本程序利用节点阻抗矩阵进行三相短路计算
clc
n=input（‘请输入本网络节点数目n=‘）;
n1=input（‘请输入本网络支路数目n1=‘）;
B1=input（‘请输入支路参数矩阵B1=‘）;
%第一列储存支路的首节点
%第二列储存支路的末节点
%第三列储存支路的阻抗
%第四列储存支路的对地导纳
%第五列储存变压器变比，若没有变压器则为“1”
A=input（‘请输入节点参数矩阵A=‘）;
%第一列储存节点号
%第二列储存该节点对地阻抗，若没有则输入“0”
f=input（‘请输入短路节点号f=‘）;
SB=100;
Uav=input（‘请输入电压基准值Uav=‘）;
Y=zeros（n）;%节点导纳矩阵
U=zeros（n1）;%短路后各节点电压标幺值矩阵
I=zeros（n13）;%短路后各节点电流标幺值矩阵

%形成导纳矩阵
for i=1:n
    if A（i2）~=0
        p=A（i1）;
        Y（pp）=1/A（i2）;
    end
end

for i=1:n1
    p=B1（i1）;
    q=B1（i2）;
    Y（pp）=Y（pp）+1/（B1（i3）*B1（i5）^2）+B1（i4）/2;
    Y（qq）=Y（qq）+1/B1（i3）+B1（i4）/2;
    Y（pq）=-1/（B1（i3）*B1（i5））;
    Y（qp）=Y（pq）;
end
Y
Z=inv（Y）

%计算短路电流
if f<=n
    If=1/Z（ff）
else
    dip（‘节点号录入错误！‘）
end
%计算短路后各节点的电压标幺值
for i=1:n