From 920fcc1cd5045a84ab5530c5450db3bbb83f1c8b Mon Sep 17 00:00:00 2001
From: "dmy@lab" <dmy@lab.lab>
Date: Sun, 5 Apr 2015 21:16:57 +0800
Subject: [PATCH] =?UTF-8?q?=E7=94=B1=E4=BA=8E=E4=B9=8B=E5=89=8D100%?=
 =?UTF-8?q?=E6=83=85=E5=86=B5=E4=B8=8D=E6=94=B6=E6=95=9B=EF=BC=8C=E5=86=8D?=
 =?UTF-8?q?=E6=AC=A1=E6=94=B9=E4=B8=BA=E8=BF=99=E4=B8=AA=E4=B8=8D=E7=AD=89?=
 =?UTF-8?q?=E5=BC=8F=E4=B8=8A=E4=B8=8B=E7=95=8C=E3=80=82?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: dmy@lab <dmy@lab.lab>
---
 FormLw.m | 4 ++--
 FormLz.m | 4 ++--
 OPF.m    | 6 +++---
 3 files changed, 7 insertions(+), 7 deletions(-)

diff --git a/FormLw.m b/FormLw.m
index edab99c..d16c155 100644
--- a/FormLw.m
+++ b/FormLw.m
@@ -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;
diff --git a/FormLz.m b/FormLz.m
index f1eaaaa..b429825 100644
--- a/FormLz.m
+++ b/FormLz.m
@@ -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;
diff --git a/OPF.m b/OPF.m
index 0eb42ef..cbc0be3 100644
--- a/OPF.m
+++ b/OPF.m
@@ -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