facat
/
pes2014-wronglowervoltagebound
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

不等式约束 

Mat_G=[
    sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
    sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
    sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;
    sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    ];

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk 2014-05-22 09:46:10 +08:00
parent 6fcf53a485
commit 305f059cf6
5 changed files with 8 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
FormG.m
View File

@ -1,10 +1,10 @@
function Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi,QDbi,mPD)
function Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi,QDbi,mPD,mQD)
Mat_G=[
sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
sparse(QD(Loadi));
sparse(QD(Loadi))-0.001;
sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;
sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;
Volt'-mVolt'-bigM*Vbi-0.03;
Volt'-mVolt'+bigM*Vbi+0.03;
Vbi;

2
FormLw.m
View File

@ -24,7 +24,7 @@ QDU(indQD(3:12:end))=1.55*realQD(indQD(3:12:end));
QDU(indQD(9:12:end))=1.05*realQD(indQD(9:12:end));
% PF=0.85;
% QDU=1.0*PD(Loadi).*sqrt(1 -PF.^2)./PF;
t1=([0*PDU',100*ones(1,length(PDU)),QDU',QDU',0*VoltU,100*ones(1,length(VoltU)),1*ones(1,Busnum+length(Loadi)*2)])';
t1=([0*PDU',100*ones(1,length(PDU)),0*QDU',100*ones(1,length(QDU)),0*VoltU,100*ones(1,length(VoltU)),1*ones(1,Busnum+length(Loadi)*2)])';
t2=Mat_G+Init_U'-t1;
Lw=t2;

2
FormLz.m
View File

@ -21,7 +21,7 @@ indQD=find(realQD>0);
QDL(indQD(3:12:end))=0.95*realQD(indQD(3:12:end));
QDL(indQD(9:12:end))=0.95*realQD(indQD(9:12:end));
% QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF;
t1=([-100*ones(1,length(PDL)),0*PDL',QDL',QDL',-100*ones(1,length(VoltL)),0*VoltL,0*ones(1,Busnum+length(Loadi)*2)])';
t1=([-100*ones(1,length(PDL)),0*PDL',-100*ones(1,length(QDL)),0*QDL',-100*ones(1,length(VoltL)),0*VoltL,0*ones(1,Busnum+length(Loadi)*2)])';
t2=Mat_G-Init_L'-t1;
Lz=t2;

2
OPF.m
View File

@ -85,7 +85,7 @@ while(abs(Gap)>Precision)
%% 形成方程矩阵
Luu=Init_U'.*Init_W'+Init_u*ones(RestraintCount,1);
Lul=Init_L'.*Init_Z'-Init_u*ones(RestraintCount,1);
Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi,QDbi,mPD);
Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi,QDbi,mPD,mQD);
Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi);
Ly=Mat_H;
Lz=FormLz(Mat_G,Init_L,GenL,Busnum,PVQL,PD,RealPD,RealQD,Loadi,KK,PF);

4
func_deltG.m
View File

@ -54,8 +54,8 @@ dg7_dPDbi=sparse(length(Loadi),Busnum);
dgPDbi_dPDbi=sparse(eye(length(Loadi)));
dgQDbi_dPDbi=sparse(zeros(length(Loadi)));
%%
dg3_dQDbi=sparse(-zeros(sizeLoadi));
dg32_dQDbi=sparse(zeros(sizeLoadi));
dg3_dQDbi=sparse(-eye(sizeLoadi));
dg32_dQDbi=sparse(eye(sizeLoadi));
dg4_dQDbi=sparse(length(Loadi),length(Loadi));
dg42_dQDbi=dg4_dPDbi;
dg5_dQDbi=sparse(length(Loadi),Busnum);