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

合理设置上下界和初值，程序收敛性很好。 

Signed-off-by: facat <dmy@dmy-PC.(none)>
This commit is contained in:
facat 2012-05-23 09:49:38 +08:00
parent 18dea8045f
commit db68f5f795
4 changed files with 30 additions and 21 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
FormLw.asv
View File

@ -1,8 +1,13 @@
function FormLw(Mat_G,Init_L,GenU,LinePLimt) function Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD0)
PL=GenU(:,1);%发电机有功下界
QL=GenU(:,2);%发电机无功下界 PU=1*GenU(:,2);%发电机有功上界
VoltL=0.8*ones(1,5); QU=PVQU(:,1);%发电机无功上界
t1=([PL',QL',VoltL,LinePLimt'])'; VoltU=1.2*ones(1,Busnum);
t2=Mat_G-Init_L'-t1; PDU=PD0;
Lz=t2; PDU(PD0>0)=1.2*PD0()';
PDU(PDU==0)=20;
PDU=25*ones(Busnum,1)';
t1=([PU',QU',PDU,VoltU])';
t2=Mat_G+Init_U'-t1;
Lw=t2;
end end

14
FormLw.m
View File

@ -1,12 +1,14 @@
function Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD0) function Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD0)
PU=20*GenU(:,2);%发电机有功上界 PU=1*GenU(:,2);%发电机有功上界
QU=PVQU(:,1);%发电机无功上界 QU=PVQU(:,1);%发电机无功上界
VoltU=1.5*ones(1,Busnum); VoltU=1.1*ones(1,Busnum);
%PDU=20*PD0'; PDU=PD0;
%PDU(PDU==0)=20; PDU(PD0>0)=1.2*PD0(PD0>0);
PDU=25*ones(Busnum,1)'; PDU(PDU==0)=.2;
t1=([PU',QU',PDU,VoltU])'; PDU(PDU<0)=0.8*PD0(PD0<0);
%PDU=25*ones(Busnum,1)';
t1=([PU',QU',PDU',VoltU])';
t2=Mat_G+Init_U'-t1; t2=Mat_G+Init_U'-t1;
Lw=t2; Lw=t2;
end end

14
FormLz.m
View File

@ -1,12 +1,14 @@
function Lz=FormLz(Mat_G,Init_L,GenL,Busnum,PVQL,PD0) function Lz=FormLz(Mat_G,Init_L,GenL,Busnum,PVQL,PD0)
PL=0.0001*GenL(:,2);%发电机有功下界 PL=1*GenL(:,2);%发电机有功下界
QL=PVQL(:,1);%发电机无功下界 QL=PVQL(:,1);%发电机无功下界
VoltL=0.5*ones(1,Busnum); VoltL=0.9*ones(1,Busnum);
%PDL=-0.001*PD0'; PDL=PD0;
%PDL(PD0<=0)=-20; PDL(PD0>0)=0.8*PD0(PD0>0);
PDL=-25*ones(Busnum,1)'; PDL(PD0==0)=-.2;
t1=([PL',QL',PDL,VoltL])'; PDL(PD0<0)=1.2*PD0(PD0<0);
t1=([PL',QL',PDL',VoltL])';
t2=Mat_G-Init_L'-t1; t2=Mat_G-Init_L'-t1;
Lz=t2; Lz=t2;

4
OPF_Init.m
View File

@ -24,6 +24,6 @@ PG(PGi)=(tPU+tPL)/2;
QG(PVi)=(tQU+tQL)/2; QG(PVi)=(tQU+tQL)/2;
wD=ones(Busnum,1); wD=ones(Busnum,1);
wD(Balance)=0; wD(Balance)=0;
PD=.5*PD0; PD=1*PD0;
PD(PD==0)=.2; %PD(PD==0)=.2;
end end