From db68f5f795debaa538165d51f81c0c2b9f94ae7e Mon Sep 17 00:00:00 2001
From: facat <dmy@dmy-PC.(none)>
Date: Wed, 23 May 2012 09:49:38 +0800
Subject: [PATCH] =?UTF-8?q?=E5=90=88=E7=90=86=E8=AE=BE=E7=BD=AE=E4=B8=8A?=
 =?UTF-8?q?=E4=B8=8B=E7=95=8C=E5=92=8C=E5=88=9D=E5=80=BC=EF=BC=8C=E7=A8=8B?=
 =?UTF-8?q?=E5=BA=8F=E6=94=B6=E6=95=9B=E6=80=A7=E5=BE=88=E5=A5=BD=E3=80=82?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: facat <dmy@dmy-PC.(none)>
---
 FormLw.asv | 19 ++++++++++++-------
 FormLw.m   | 14 ++++++++------
 FormLz.m   | 14 ++++++++------
 OPF_Init.m |  4 ++--
 4 files changed, 30 insertions(+), 21 deletions(-)

diff --git a/FormLw.asv b/FormLw.asv
index 35736b1..150c0a5 100644
--- a/FormLw.asv
+++ b/FormLw.asv
@@ -1,8 +1,13 @@
-function FormLw(Mat_G,Init_L,GenU,LinePLimt)
-PL=GenU(:,1);%发电机有功下界
-QL=GenU(:,2);%发电机无功下界
-VoltL=0.8*ones(1,5);
-t1=([PL',QL',VoltL,LinePLimt'])';
-t2=Mat_G-Init_L'-t1;
-Lz=t2;
+function Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD0)
+
+PU=1*GenU(:,2);%发电机有功上界
+QU=PVQU(:,1);%发电机无功上界
+VoltU=1.2*ones(1,Busnum);
+PDU=PD0;
+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
\ No newline at end of file
diff --git a/FormLw.m b/FormLw.m
index 9bce6d9..5071ee3 100644
--- a/FormLw.m
+++ b/FormLw.m
@@ -1,12 +1,14 @@
 function Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD0)
 
-PU=20*GenU(:,2);%发电机有功上界
+PU=1*GenU(:,2);%发电机有功上界
 QU=PVQU(:,1);%发电机无功上界
-VoltU=1.5*ones(1,Busnum);
-%PDU=20*PD0';
-%PDU(PDU==0)=20;
-PDU=25*ones(Busnum,1)';
-t1=([PU',QU',PDU,VoltU])';
+VoltU=1.1*ones(1,Busnum);
+PDU=PD0;
+PDU(PD0>0)=1.2*PD0(PD0>0);
+PDU(PDU==0)=.2;
+PDU(PDU<0)=0.8*PD0(PD0<0);
+%PDU=25*ones(Busnum,1)';
+t1=([PU',QU',PDU',VoltU])';
 t2=Mat_G+Init_U'-t1;
 Lw=t2;
 end
\ No newline at end of file
diff --git a/FormLz.m b/FormLz.m
index 578d717..4da82f8 100644
--- a/FormLz.m
+++ b/FormLz.m
@@ -1,12 +1,14 @@
 function Lz=FormLz(Mat_G,Init_L,GenL,Busnum,PVQL,PD0)
 
-PL=0.0001*GenL(:,2);%发电机有功下界
+PL=1*GenL(:,2);%发电机有功下界
 QL=PVQL(:,1);%发电机无功下界
-VoltL=0.5*ones(1,Busnum);
-%PDL=-0.001*PD0';
-%PDL(PD0<=0)=-20;
-PDL=-25*ones(Busnum,1)';
-t1=([PL',QL',PDL,VoltL])';
+VoltL=0.9*ones(1,Busnum);
+PDL=PD0;
+PDL(PD0>0)=0.8*PD0(PD0>0);
+PDL(PD0==0)=-.2;
+PDL(PD0<0)=1.2*PD0(PD0<0);
+
+t1=([PL',QL',PDL',VoltL])';
 
 t2=Mat_G-Init_L'-t1;
 Lz=t2;
diff --git a/OPF_Init.m b/OPF_Init.m
index df39a9c..2451b89 100644
--- a/OPF_Init.m
+++ b/OPF_Init.m
@@ -24,6 +24,6 @@ PG(PGi)=(tPU+tPL)/2;
 QG(PVi)=(tQU+tQL)/2;
 wD=ones(Busnum,1);
 wD(Balance)=0;
-PD=.5*PD0;
-PD(PD==0)=.2;
+PD=1*PD0;
+%PD(PD==0)=.2;
 end
\ No newline at end of file