clear;
%clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%初始化%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%數據加載
n=input（‘請輸入要計算的節點系統（5）:‘）
????load?Node5.txt;%節點數據
????load?Branch5.txt;%支路數據
????load?Generator5.txt;%發電機數據
????Node=Node5;
????Branch=Branch5;
????Generator=Generator5;

%節點數據處理
N=Node（:1）;%節點號
Type=Node（:2）;%節點類型
Uamp=Node（:3）;%節點電壓幅值
Dlta=Node（:4）;%節點電壓相角
Pd=Node（:5）;%節點負荷有功
Qd=Node（:6）;%節點負荷無功
Pg=Node（:7）;%節點出力有功
Qg=Node（:8）;%節點出力無功
Umax=Node（:9）;%節點電壓幅值上限
Umin=Node（:10）;%節點電壓幅值下限
Bc=Node（:11）;%節點補償電容電納值
%支路數據處理
Nbr=Branch（:1）;%支路號
Nl=Branch（:2）;%支路首節點
Nr=Branch（:3）;%支路末節點
R=Branch（:4）;%支路電阻
X=Branch（:5）;%支路電抗
Z=R+1i*X;%支路阻抗=支路電阻+支路電抗
Bn=Branch（:6）;%支路對地電納
K=Branch（:7）;%支路變壓器變比，0表示無變壓器
Ptmax=Branch（:8）;%線路傳輸功率上限

%發電機數據處理
Ng=Generator（: