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

为了写论文，又改为用实际值正负20%作为负荷不等式约束上下限。

This commit is contained in:
facat 2013-01-21 18:28:16 +08:00
parent 5dc81f33bb
commit 7ded138066
4 changed files with 36 additions and 26 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

33
FormLw.asv
View File

@ -1,15 +1,24 @@
function Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD0,Loadi,KK)
%KK=999;
PU=1*GenU(:,2);%发电机有功上界
QU=PVQU(:,1);%发电机无功上界
%VoltU=(1.5+1/exp(KK))*ones(1,Busnum);
VoltU=1.5*ones(1,Busnum);
PDU=PD0(Loadi);
% PDU(PDU>0)=1200*PDU(PDU>0);
% PDU(PDU<0)=-800*PDU(PDU<0);
% PDU(PDU==0)=200;
PDU=1000000*ones(length(Loadi),1);
t1=([PU',QU',PDU',VoltU])';
function Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD,PD0,QD0,Loadi,KK,PF,noDataTransCapacity)
KK=999;
%PU=GenU(:,2);%发电机有功上界
PU=5*ones(length(GenU(:,2)),1);
%QU=PVQU(:,1);%发电机无功上界
QU=5*ones(length(PVQU(:,1)),1);
VoltU=(1.1)*ones(1,Busnum);
%VoltU=10*ones(1,Busnum);
% PDU=PD0(Loadi);
PDU=noDataTransCapacity;
PDU(PDU>0)=1.200*PDU(PDU>0);
PDU(PDU<0)=0.800*PDU(PDU<0);
PDU(PDU==0)=0.400;
%PDU=10*ones(length(Loadi),1);
% QDU=QD0(Loadi);
% QDU(QDU>0)=1.200*QDU(QDU>0);
% QDU(QDU<0)=0.800*QDU(QDU<0);
% QDU(QDU==0)=0.200;
% PF=0.85;
% QDU=1.0*PD(Loadi).*sqrt(1 -PF.^2)./PF;
t1=([PU',QU',PDU',QDU',VoltU])';
t2=Mat_G+Init_U'-t1;
Lw=t2;

16
FormLw.m
View File

@ -6,18 +6,18 @@ PU=5*ones(length(GenU(:,2)),1);
QU=5*ones(length(PVQU(:,1)),1);
VoltU=(1.1)*ones(1,Busnum);
%VoltU=10*ones(1,Busnum);
% PDU=PD0(Loadi);
PDU=noDataTransCapacity;
PDU=PD0(Loadi);
% PDU=noDataTransCapacity;
PDU(PDU>0)=1.200*PDU(PDU>0);
PDU(PDU<0)=0.800*PDU(PDU<0);
PDU(PDU==0)=0.400;
%PDU=10*ones(length(Loadi),1);
% QDU=QD0(Loadi);
% QDU(QDU>0)=1.200*QDU(QDU>0);
% QDU(QDU<0)=0.800*QDU(QDU<0);
% QDU(QDU==0)=0.200;
PF=0.85;
QDU=1.0*PD(Loadi).*sqrt(1 -PF.^2)./PF;
QDU=QD0(Loadi);
QDU(QDU>0)=1.200*QDU(QDU>0);
QDU(QDU<0)=0.800*QDU(QDU<0);
QDU(QDU==0)=0.200;
% PF=0.85;
% QDU=1.0*PD(Loadi).*sqrt(1 -PF.^2)./PF;
t1=([PU',QU',PDU',QDU',VoltU])';
t2=Mat_G+Init_U'-t1;
Lw=t2;

10
FormLz.m
View File

@ -11,11 +11,11 @@ PDL(PDL>0)=0.800*PDL(PDL>0);
PDL(PDL<0)=1.200*PDL(PDL<0);
PDL(PDL==0)=-0.400;
%PDL=-10*ones(length(Loadi),1);
% QDL=QD0(Loadi);
% QDL(QDL>0)=0.800*QDL(QDL>0);
% QDL(QDL<0)=1.200*QDL(QDL<0);
% QDL(QDL==0)=-0.200;
QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF;
QDL=QD0(Loadi);
QDL(QDL>0)=0.800*QDL(QDL>0);
QDL(QDL<0)=1.200*QDL(QDL<0);
QDL(QDL==0)=-0.200;
% QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF;
t1=([PL',QL',PDL',QDL',VoltL])';
t2=Mat_G-Init_L'-t1;
Lz=t2;

3
OPF.m
View File

@ -51,7 +51,8 @@ QD0(Loadi)=QD0(Loadi).*(1+normrnd(0,0.05,length(Loadi),1));
PG0(PGi)=PG0(PGi).*(1+normrnd(0,0.01,length(PGi),1));
QG0(PVi)=QG0(PVi).*(1+normrnd(0,0.01,length(PVi),1));
%% 读变压器容量
[noDataTransNum noDataTransCapacity noDataTransPowerFactor]=ReadNoDataTrans('C:/b/东际911_2751267_2012-09-05/iPso_东际911_2751267_2012-09-05_变压器无负载.txt');
%[noDataTransNum noDataTransCapacity noDataTransPowerFactor]=ReadNoDataTrans('C:/b/东际911_2751267_2012-09-05/iPso_东际911_2751267_2012-09-05_变压器无负载.txt');
noDataTransCapacity=0;
while(abs(Gap)>Precision)
if KK>kmax
break;