clear;?%clc
%數據加載
c=input（‘請輸入要計算的節點系統（6）:‘）
??load?Node6.txt;%節點數據
??load?Branch6.txt;%支路數據
??load?Generator6.txt;%發電機數據
??Node=Node6;
??Branch=Branch6;
??Generator=Generator6;

??
%節點數據處理
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（:1）;%發電機序號
Nbus=Generator（:2）;%所在母線號
Pumax=Generator（:3）;%發電機有功出力上界
Qumax=Generato