facat
/
stateestimateionyalmip-lu9-forcomparision
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

只约束0注入量 

Signed-off-by: facat <dugg@21cn.com>
This commit is contained in:
facat 2013-04-06 21:22:56 +08:00
parent 7b85093fb1
commit 8c79065f82
1 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
Run_YALMIP.m
View File

@ -35,6 +35,7 @@ maxdPQ=max([dPD(dPD<10);dQD(dQD<10)]);
[Volt,UAngel,Init_Z,Init_W,Init_L,Init_U,Init_Y,PG,QG,RestraintCount,wPG,wQG,wPD,wQD,PD,PD0,QD,randPDind,Loadi,notLoadi]=OPF_Init(Busnum,Balance,PG,QG,Volt,GenU,GenL,PVi,PGi,PVQU,PVQL,PD0,QD0,QD,PD); [Volt,UAngel,Init_Z,Init_W,Init_L,Init_U,Init_Y,PG,QG,RestraintCount,wPG,wQG,wPD,wQD,PD,PD0,QD,randPDind,Loadi,notLoadi]=OPF_Init(Busnum,Balance,PG,QG,Volt,GenU,GenL,PVi,PGi,PVQU,PVQL,PD0,QD0,QD,PD);
%% ¨Òå±äÁ¿ %% ¨Òå±äÁ¿
BalVolt=Volt(Balance);
Volt=sdpvar(Busnum,1); Volt=sdpvar(Busnum,1);
UAngel=sdpvar(Busnum,1); UAngel=sdpvar(Busnum,1);
% PG=sdpvar(Busnum,1); % PG=sdpvar(Busnum,1);
@ -55,13 +56,16 @@ assign(QD(:),QD0(:));
%% YALMIP²¿·Ö %% YALMIP²¿·Ö
dP=PG0-PD-diag(Volt)*Y.*cos( sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum,Busnum) )*Volt; dP=PG0-PD-diag(Volt)*Y.*cos( sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum,Busnum) )*Volt;
dQ=QG0-QD-diag(Volt)*Y.*sin( sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum,Busnum) )*Volt; dQ=QG0-QD-diag(Volt)*Y.*sin( sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum,Busnum) )*Volt;
Loadi=PD0~=0 | QD0~=0; Loadi=PD0~=0 | QD0~=0 |PG0~=0|QG0~=0;
Constraints = [%AngleIJ-sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum,Busnum)==0, ... Constraints = [%AngleIJ-sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum,Busnum)==0, ...
dP==0, ... dP(setdiff(1:Busnum,Loadi))==0, ...
dQ==0, ... dQ(setdiff(1:Busnum,Loadi))==0, ...
% dP==0, ...
% dQ==0, ...
PD(PD0==0)==0, ... PD(PD0==0)==0, ...
QD(QD0==0)==0, ... QD(QD0==0)==0, ...
0.9*ones(Busnum,1)<=Volt<=1.1*ones(Busnum,1), ... 0.9*ones(Busnum,1)<=Volt<=1.1*ones(Busnum,1), ...
Volt(Balance)==BalVolt;
]; ];
options = sdpsettings('verbose',2,'showprogress',1,'debug',0,'solver','ipopt','usex0','1'); options = sdpsettings('verbose',2,'showprogress',1,'debug',0,'solver','ipopt','usex0','1');
sol = solvesdp(Constraints,Objective,options); sol = solvesdp(Constraints,Objective,options);