From 8c79065f827631839206ad45bf1d73714cc8fc2a Mon Sep 17 00:00:00 2001
From: facat <dugg@21cn.com>
Date: Sat, 6 Apr 2013 21:22:56 +0800
Subject: [PATCH] =?UTF-8?q?=E5=8F=AA=E7=BA=A6=E6=9D=9F0=E6=B3=A8=E5=85=A5?=
 =?UTF-8?q?=E9=87=8F?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: facat <dugg@21cn.com>
---
 Run_YALMIP.m | 10 +++++++---
 1 file changed, 7 insertions(+), 3 deletions(-)

diff --git a/Run_YALMIP.m b/Run_YALMIP.m
index 4f323e5..8809520 100644
--- a/Run_YALMIP.m
+++ b/Run_YALMIP.m
@@ -35,6 +35,7 @@ maxdPQ=max([dPD(dPD<10);dQD(dQD<10)]);
 [Volt,UAngel,Init_Z,Init_W,Init_L,Init_U,Init_Y,PG,QG,RestraintCount,wPG,wQG,wPD,wQD,PD,PD0,QD,randPDind,Loadi,notLoadi]=OPF_Init(Busnum,Balance,PG,QG,Volt,GenU,GenL,PVi,PGi,PVQU,PVQL,PD0,QD0,QD,PD);
 
 %% 定义变量
+BalVolt=Volt(Balance);
 Volt=sdpvar(Busnum,1);
 UAngel=sdpvar(Busnum,1);
 % PG=sdpvar(Busnum,1);
@@ -55,13 +56,16 @@ assign(QD(:),QD0(:));
 %% YALMIP部分
 dP=PG0-PD-diag(Volt)*Y.*cos( sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum,Busnum) )*Volt;
 dQ=QG0-QD-diag(Volt)*Y.*sin( sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum,Busnum) )*Volt;
-Loadi=PD0~=0 | QD0~=0;
+Loadi=PD0~=0 | QD0~=0 |PG0~=0|QG0~=0;
 Constraints = [%AngleIJ-sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum,Busnum)==0, ...
-    dP==0, ...
-    dQ==0, ...
+    dP(setdiff(1:Busnum,Loadi))==0, ...
+    dQ(setdiff(1:Busnum,Loadi))==0, ...
+%     dP==0, ...
+%     dQ==0, ...
      PD(PD0==0)==0, ...
      QD(QD0==0)==0, ...
     0.9*ones(Busnum,1)<=Volt<=1.1*ones(Busnum,1), ...
+    Volt(Balance)==BalVolt;
     ];
 options = sdpsettings('verbose',2,'showprogress',1,'debug',0,'solver','ipopt','usex0','1');
 sol = solvesdp(Constraints,Objective,options);