From cc7312d070839ee401e00bdbc7927e1fc3d3f64b Mon Sep 17 00:00:00 2001
From: "dmy@lab" <dmy@lab.lab>
Date: Fri, 3 Apr 2015 10:06:03 +0800
Subject: [PATCH] =?UTF-8?q?1.=E4=B8=8A=E4=B8=8B=E7=95=8C=E6=94=B9=E4=B8=BA?=
 =?UTF-8?q?1.2=E5=88=B00.8=202.=E7=94=A8=E9=87=8F=E6=B5=8B=E5=80=BC?=
 =?UTF-8?q?=E7=A1=AE=E5=AE=9A=E4=B8=8A=E4=B8=8B=E7=95=8C?=
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 | 6 +++---
 FormLz.m | 6 +++---
 OPF.m    | 6 +++---
 3 files changed, 9 insertions(+), 9 deletions(-)

diff --git a/FormLw.m b/FormLw.m
index c411549..cf4f829 100644
--- a/FormLw.m
+++ b/FormLw.m
@@ -1,4 +1,4 @@
-function Lw=FormLw(Mat_G,Init_U,Loadi,rPD3P,rQD3P)
+function Lw=FormLw(Mat_G,Init_U,Loadi,mPD3P,mQD3P)
 % KK=999;
 % %PU=GenU(:,2);%发电机有功上界
 % 
@@ -26,8 +26,8 @@ function Lw=FormLw(Mat_G,Init_U,Loadi,rPD3P,rQD3P)
 % 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.1*rPD3P;
-QDU=1.1*rQD3P;
+PDU=1.2*mPD3P;
+QDU=1.2*mQD3P;
 t1=([PDU',QDU'])';
 t2=Mat_G+Init_U'-t1;
 Lw=t2;
diff --git a/FormLz.m b/FormLz.m
index 5627360..fb4d880 100644
--- a/FormLz.m
+++ b/FormLz.m
@@ -1,4 +1,4 @@
-function Lz=FormLz(Mat_G,Init_L,Loadi,rPD3P,rQD3P)
+function Lz=FormLz(Mat_G,Init_L,Loadi,mPD3P,mQD3P)
 % KK=999;
 % 
 % VoltL=(0.9)*ones(1,Busnum);
@@ -23,8 +23,8 @@ function Lz=FormLz(Mat_G,Init_L,Loadi,rPD3P,rQD3P)
 % QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF;
 % PDL=0*ones(length(Loadi),1);
 % QDL=0*ones(length(Loadi),1);
-PDL=0.9*rPD3P;
-QDL=0.9*rQD3P;
+PDL=0.8*mPD3P;
+QDL=0.8*mQD3P;
 t1=([PDL',QDL'])';
 t2=Mat_G-Init_L'-t1;
 Lz=t2;
diff --git a/OPF.m b/OPF.m
index 02c3358..df33f88 100644
--- a/OPF.m
+++ b/OPF.m
@@ -210,7 +210,7 @@ rQD3P(3:3:end)=phaseCSpotLoadQ;
 rPD3P=rPD3P(Loadi);
 rQD3P=rQD3P(Loadi);
 %量测量
-sigma=0.03;
+sigma=0.01;
 VoltSigma=(1+normrnd(0,sigma,length(rVoltABCV),1));
 mVoltABCV=rVoltABCV.*VoltSigma;
 PD3PSigma=(1+normrnd(0,sigma,length(rPD3P),1));
@@ -272,8 +272,8 @@ while(abs(Gap)>Precision)
     Mat_G=FormG(PD3P,QD3P,Loadi);
     Mat_H=FormH(Busnum,Volt,PG3P,PD3P,QG3P,QD3P,Y,UAngel,r,c,Yangle,Loadi);
     Ly=Mat_H;
-    Lz=FormLz(Mat_G,Init_L,Loadi,rPD3P,rQD3P);
-    Lw=FormLw(Mat_G,Init_U,Loadi,rPD3P,rQD3P);
+    Lz=FormLz(Mat_G,Init_L,Loadi,mPD3P,mQD3P);
+    Lw=FormLw(Mat_G,Init_U,Loadi,mPD3P,mQD3P);
     Lx=FormLx(deltF,deltH,Init_Y,deltG,Init_Z,Init_W);
     %     YY=FormYY(Lul,Lz,Ly,Luu,Lw,Lx);
     %% 开始解方程