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

%數(shù)據(jù)加載
n=input（‘請輸入要計算的節(jié)點系統(tǒng)（514）:‘）
if?n==5
????load?Node5.txt;%節(jié)點數(shù)據(jù)
????load?Branch5.txt;%支路數(shù)據(jù)
????load?Generator5.txt;%發(fā)電機(jī)數(shù)據(jù)
????Node=Node5;
????Branch=Branch5;
????Generator=Generator5;
else?if?n==14
????load?Node14.txt;%節(jié)點數(shù)據(jù)
????load?Branch14.txt;%支路數(shù)據(jù)
????load?Generator14.txt;%發(fā)電機(jī)數(shù)據(jù)
????Node=Node14;
????Branch=Branch14;
????Generator=Generator14;
????end
end

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

%發(fā)電機(jī)數(shù)據(jù)處理
Ng=Generator（:1）;%發(fā)電機(jī)序號
Nbus=Generator（:2）;%所在母線號
Pumax=Generator（:3）;%發(fā)電機(jī)有功出力上界
Qumax=Generator（:4）;%發(fā)電機(jī)無功出力上界
Pumin=Generator（:5）;%發(fā)電機(jī)有功出力下界
Qumin=Generator（:6）;%發(fā)電機(jī)無功出力下界
a2=Generator（:7）;%燃料耗費曲線二次系數(shù)
a1=Generator（:8）;%燃料耗費曲線一次系數(shù)
a0=Generator（:9）;%燃料耗費曲線常數(shù)項

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

n=length（N）;%節(jié)點個數(shù)
ng=length（Ng）;%發(fā)電機(jī)臺數(shù)
nbr=length（Nbr）;%支路個數(shù)

x=zeros（2*（ng+n）1）;%控制變量+狀態(tài)變量
x（1:ng）=Pg（Nbus）;
x（ng+1:2*ng）=Qg（Nbus）;
x（（2*ng+2）:2:2*（ng+n））=Uamp;
x（（2*ng+1）:2:2*（ng+n）-1）=Dlta;
l=0.8*ones（2*ng+n+nbr1）;%松弛變量
u=1.1*ones（2*ng+n+nbr1）;%松弛變量

w=-1.5*ones（2*ng+n+nbr1）;%拉格朗日乘子
z=ones（2*ng+n+nbr1）;%拉格朗日乘子
y=zeros（2*n1）;%拉格朗日乘子
y（1:2:2*n-1）=1e-3;
y（2:2:2*n）=-1e-3;


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%計算不等式約束的上下限%%%%%%%%%%%%%%%%%%%%%%%%
%gmin
gmin=zeros（2*ng+n+nbr1）;
gmin（1:ng）=Pumin;
gmin（ng+1:2*ng）=Qumin;
gmin（2*ng+1:2*ng+n）=Umin;
gmin（2*ng+n+1:2*ng+n+nbr）=-Ptmax;??
%gmax
gmax=zeros（2*ng+n+nbr1）;
gmax（1:ng）=Pumax;
gmax（ng+1:2*ng）=Qumax;
gmax（2*ng+1:2*ng+n）=Umax;
gmax（2*ng+n+1:2*ng+n+nbr）=Ptmax;????????
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%形成導(dǎo)納矩陣%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Y=zeros（nn）;
%%%%%%%%%%%%%%%%%%%%計算非對角元素%%%%%%%%%%%%%%%%%%%%%
for?ii=1:nbr
????if?K（ii）==0%非變壓器支路
????????Y（Nl（ii）Nr（ii））=-1/Z（ii）;
????????Y（Nr（ii）Nl（ii））=Y（Nl（ii）Nr（ii））;
????else%變壓器支路
????????Y（Nl（ii）Nr（ii））=-1/Z（ii）/K（ii）;
????????Y（Nr（ii）Nl（ii））=?Y（Nl（ii）Nr（ii））;
????end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%計算對角元素%%%%%%%%%%%%%%%%%%%%%%
for?ii=1:n%將支路導(dǎo)納加入到對角元素中
????for?jj=1:nbr
????????if?K（jj）==0&&（Nl（jj）==ii||Nr（jj）==ii）%非變壓器支路
????????????????Y（iiii）=Y（iiii）+1/Z（jj）;
????????else?if?K（jj）~=0&&（Nl（jj）==ii||Nr（jj）==ii）%變壓器支路
????????????????Y（iiii）=Y（iiii）+1/Z（jj）/K（jj）;
????????????end
????????end
????end
end

for?ii=1:nbr%將對地電納加入到對角元素中
????if?K（ii）==0%非變壓器支路
????????Y（Nl（ii）Nl（ii））=Y（Nl（ii）Nl（ii））+1i*Bn（ii）;
????????Y（Nr（ii）Nr（ii））=Y（Nr（ii）N

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1030??2010-08-20?16:07??IPOPF\Branch14.txt

?????文件????????442??2010-08-20?11:01??IPOPF\Branch5.txt

?????文件????????484??2010-06-19?15:16??IPOPF\Generator14.txt

?????文件????????291??2010-07-16?19:12??IPOPF\Generator5.txt

?????文件???????1388??2010-08-21?07:40??IPOPF\Node14.txt

?????文件????????662??2010-08-23?07:51??IPOPF\Node5.txt

?????文件??????20064??2011-02-25?09:12??IPOPF\opf2.m

?????文件????????162??2011-05-31?16:23??IPOPF\說明.txt

?????目錄??????????0??2018-01-16?16:29??IPOPF

-----------?---------??----------?-----??----

????????????????24523????????????????????9