%本程序利用節點阻抗矩陣進行三相短路計算
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
????