clear;?%clc
%數(shù)據(jù)加載
c=input（‘請(qǐng)輸入要計(jì)算的節(jié)點(diǎn)系統(tǒng)（6）:‘）
??load?Node6.txt;%節(jié)點(diǎn)數(shù)據(jù)
??load?Branch6.txt;%支路數(shù)據(jù)
??load?Generator6.txt;%發(fā)電機(jī)數(shù)據(jù)
??Node=Node6;
??Branch=Branch6;
??Generator=Generator6;

??
%節(jié)點(diǎn)數(shù)據(jù)處理
N=Node（:1）;%節(jié)點(diǎn)號(hào)
Type=Node（:2）;%節(jié)點(diǎn)類(lèi)型
Uamp=Node（:3）;%節(jié)點(diǎn)電壓幅值
Dlta=Node（:4）;%節(jié)點(diǎn)電壓相角
Pd=Node（:5）;%節(jié)點(diǎn)負(fù)荷有功
Qd=Node（:6）;%節(jié)點(diǎn)負(fù)荷無(wú)功
Pg=Node（:7）;%節(jié)點(diǎn)出力有功
Qg=Node（:8）;%節(jié)點(diǎn)出力無(wú)功
Umax=Node（:9）;%節(jié)點(diǎn)電壓幅值上限
Umin=Node（:10）;%節(jié)點(diǎn)電壓幅值下限
Bc=Node（:11）;%節(jié)點(diǎn)補(bǔ)償電容電納值
%支路數(shù)據(jù)處理
Nbr=Branch（:1）;%支路號(hào)
Nl=Branch（:2）;%支路首節(jié)點(diǎn)
Nr=Branch（:3）;%支路末節(jié)點(diǎn)
R=Branch（:4）;%支路電阻
X=Branch（:5）;%支路電抗
Z=R+1i*X;%支路阻抗=支路電阻+支路電抗
Bn=Branch（:6）;%支路對(duì)地電納
K=Branch（:7）;%支路變壓器變比，0表示無(wú)變壓器
Ptmax=Branch（:8）;%線(xiàn)路傳輸功率上限

%發(fā)電機(jī)數(shù)據(jù)處理
Ng=Generator（:1）;%發(fā)電機(jī)序號(hào)
Nbus=Generator（:2）;%所在母線(xiàn)號(hào)
Pumax=Generator（:3）;%發(fā)電機(jī)有功出力上界
Qumax=Generato