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

由于之前100%情况不收敛，再次改为这个不等式上下界。 

Signed-off-by: dmy@lab <dmy@lab.lab>
This commit is contained in:
dmy@lab 2015-04-05 21:16:57 +08:00
parent 8306ac80f7
commit 920fcc1cd5
3 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
FormLw.m
View File

@ -26,8 +26,8 @@ function Lw=FormLw(Mat_G,Init_U,Loadi,rPD3P,rQD3P,mVoltABCV)
% QDU=1.0*PD(Loadi).*sqrt(1 -PF.^2)./PF;
% PDU=1.8*ones(length(Loadi),1);
% QDU=1.8*ones(length(Loadi),1);
PDU=1.20*rPD3P.*(1+unifrnd(-0.2,0.2,length(rPD3P),1));
QDU=1.20*rQD3P.*(1+unifrnd(-0.2,0.2,length(rPD3P),1));
PDU=1.25*rPD3P.*(1+unifrnd(-0.2,0.2,length(rPD3P),1));
QDU=1.25*rQD3P.*(1+unifrnd(-0.2,0.2,length(rPD3P),1));
% t1=([PDU',QDU',1*ones(length(Volt),1)'])';
t1=([PDU',QDU'])';
t2=Mat_G+Init_U'-t1;

4
FormLz.m
View File

@ -23,8 +23,8 @@ function Lz=FormLz(Mat_G,Init_L,Loadi,rPD3P,rQD3P,mVoltABCV)
% QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF;
% PDL=0*ones(length(Loadi),1);
% QDL=0*ones(length(Loadi),1);
PDL=0.8*rPD3P.*(1+unifrnd(-0.15,0.15,length(rPD3P),1));
QDL=0.8*rQD3P.*(1+unifrnd(-0.15,0.15,length(rQD3P),1));
PDL=0.75*rPD3P.*(1+unifrnd(-0.15,0.15,length(rPD3P),1));
QDL=0.75*rQD3P.*(1+unifrnd(-0.15,0.15,length(rQD3P),1));
% t1=([PDL',QDL',0.88*ones(length(Volt),1)'])';
t1=([PDL',QDL'])';
t2=Mat_G-Init_L'-t1;

6
OPF.m
View File

@ -225,10 +225,10 @@ wQD=1./(abs(mQD3P*sigma).^2);
%% 没有量测量
% noLoadi=[1,5,6,10,11];
noLoadi=[1,8,9,10,12];
% noLoadi=[1,8,9,10,12];
% noLoadi=Loadi;
% noLoadi=[1,11];
% noLoadi=[];
noLoadi=[];
noPQi3P=zeros(length(noLoadi)*3,1);
noPQi3P(1:3:end)=(noLoadi-1)*3+1;
noPQi3P(2:3:end)=(noLoadi-1)*3+2;
@ -255,7 +255,7 @@ Precision=1e-5;
CenterA=0.1;
%% 加误差
Volt=1*Vp3m;
Volt=mVoltABCV;
Volt=rVoltABCV;
UAngel=Vp3a;
while(abs(Gap)>Precision)
if KK>kmax