From 98502be3112afe010dbf51680eca9509683c2035 Mon Sep 17 00:00:00 2001
From: "dmy@lab" <dmy@lab.lab>
Date: Sun, 5 Apr 2015 20:22:38 +0800
Subject: [PATCH] =?UTF-8?q?=E7=94=A8=E8=BF=99=E4=B8=AA=E6=B2=A1=E6=9C=89?=
 =?UTF-8?q?=E9=87=8F=E6=B5=8B=E9=87=8F=E7=9A=84=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 | 9 +++++----
 1 file changed, 5 insertions(+), 4 deletions(-)

diff --git a/OPF.m b/OPF.m
index b1f4374..71e521b 100644
--- a/OPF.m
+++ b/OPF.m
@@ -225,7 +225,8 @@ wQD=1./(abs(mQD3P*sigma).^2);
 
 %% 没有量测量
 % noLoadi=[1,5,6,10,11];
-noLoadi=Loadi;
+noLoadi=[1,8,9,10,12];
+% noLoadi=Loadi;
 % noLoadi=[1,11];
 % noLoadi=[];
 noPQi3P=zeros(length(noLoadi)*3,1);
@@ -234,9 +235,9 @@ noPQi3P(2:3:end)=(noLoadi-1)*3+2;
 noPQi3P(3:3:end)=(noLoadi-1)*3+3;
 % noLoadi=noPQi3P;
 wVolt(setdiff(1:length(wVolt),Loadi))=0;%只有负荷处才有电压量测。
-wVolt( noLoadi)=0;
-wPD(ismember( Loadi,noLoadi))=0;
-wQD(ismember( Loadi,noLoadi))=0;
+wVolt( noPQi3P)=0;
+wPD(ismember( Loadi,noPQi3P))=0;
+wQD(ismember( Loadi,noPQi3P))=0;
 %%
 % RestraintCount=size(Loadi,1)*2+length(rVoltABCV); %约束条件数
 RestraintCount=size(Loadi,1)*2; %约束条件数