From e26a2f9ef064316baa55e61b400d0f20133f399a Mon Sep 17 00:00:00 2001
From: "dmy@lab" <dmy@lab.lab>
Date: Sun, 5 Apr 2015 19:31:39 +0800
Subject: [PATCH] =?UTF-8?q?1.=E4=BF=AE=E6=94=B9=E4=BA=86=E8=BF=AD=E4=BB=A3?=
 =?UTF-8?q?=E6=9C=80=E5=A4=A7=E6=AC=A1=E6=95=B0=202.=E5=85=88=E7=94=A8?=
 =?UTF-8?q?=E8=BF=99=E4=B8=AA=E7=AE=97=E4=BE=8B=E6=95=B0=E6=8D=AE=E5=92=8C?=
 =?UTF-8?q?=E5=8F=82=E6=95=B0=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>
---
 OPF.m | 13 ++++---------
 1 file changed, 4 insertions(+), 9 deletions(-)

diff --git a/OPF.m b/OPF.m
index 88de316..b1f4374 100644
--- a/OPF.m
+++ b/OPF.m
@@ -224,24 +224,19 @@ wPD=1./(abs(mPD3P*sigma).^2);
 wQD=1./(abs(mQD3P*sigma).^2);
 
 %% 没有量测量
-noLoadi=[1,5,6,10,11];
-% noLoadi=Loadi;
+% noLoadi=[1,5,6,10,11];
+noLoadi=Loadi;
 % noLoadi=[1,11];
 % noLoadi=[];
 noPQi3P=zeros(length(noLoadi)*3,1);
 noPQi3P(1:3:end)=(noLoadi-1)*3+1;
 noPQi3P(2:3:end)=(noLoadi-1)*3+2;
 noPQi3P(3:3:end)=(noLoadi-1)*3+3;
-noLoadi=noPQi3P;
+% noLoadi=noPQi3P;
 wVolt(setdiff(1:length(wVolt),Loadi))=0;%只有负荷处才有电压量测。
 wVolt( noLoadi)=0;
 wPD(ismember( Loadi,noLoadi))=0;
 wQD(ismember( Loadi,noLoadi))=0;
-
-
-% wVolt(5)=0;
-% wPD(5)=0;
-% wQD(5)=0;
 %%
 % RestraintCount=size(Loadi,1)*2+length(rVoltABCV); %约束条件数
 RestraintCount=size(Loadi,1)*2; %约束条件数
@@ -253,7 +248,7 @@ Init_Y=sparse(1,2*Busnum*3);%
 Gap=(Init_L*Init_Z'-Init_U*Init_W');
 KK=0;
 ContrlCount=size(Loadi,1)*2+Busnum*6;
-kmax=500;
+kmax=200;
 %%
 Precision=1e-5;
 CenterA=0.1;