From 1bce6a7520a9d48d6a50de3b9344f33173c0f486 Mon Sep 17 00:00:00 2001
From: "dugg@lab-desk" <dugg@lab-desk>
Date: Sat, 7 Feb 2015 15:38:00 +0800
Subject: [PATCH] =?UTF-8?q?=E4=B8=8A=E4=B8=80=E4=B8=AA=E6=8F=90=E4=BA=A4?=
 =?UTF-8?q?=E6=8F=90=E5=89=8D=E5=81=9C=E6=AD=A2=E8=BF=AD=E4=BB=A3=E4=BA=86?=
 =?UTF-8?q?=EF=BC=8C=E9=83=BD=E6=B2=A1=E6=9C=89=E6=94=B6=E6=95=9B=EF=BC=8C?=
 =?UTF-8?q?=E4=BA=92=E8=A1=A5=E7=BA=A6=E6=9D=9F=E4=B8=8D=E4=B8=BA0?=
 =?UTF-8?q?=EF=BC=8C=E7=8E=B0=E5=9C=A8=E8=BF=99=E4=B8=AA=E5=8F=AF=E4=BB=A5?=
 =?UTF-8?q?=E4=BA=86=EF=BC=8C=E5=8A=A0=E4=BA=86=E4=BA=92=E8=A1=A5=E5=8F=98?=
 =?UTF-8?q?=E9=87=8F=E7=9A=84=E5=88=9D=E5=A7=8B=E5=8C=96=E3=80=82?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
---
 subOPF.m | 16 +++++++++++++---
 1 file changed, 13 insertions(+), 3 deletions(-)

diff --git a/subOPF.m b/subOPF.m
index 04591da..04a2667 100644
--- a/subOPF.m
+++ b/subOPF.m
@@ -42,6 +42,9 @@ QG(PVi)=QGBal(PVi);
 
 %%
 [Volt,UAngel,Init_Z,Init_W,Init_L,Init_U,Init_Y,PG,QG,RestraintCount,wPG,wQG,wPD,wQD,PD,PD00,QD,randPDind,Loadi]=OPF_Init(Busnum,Balance,PG,QG,Volt,GenU,GenL,PVi,PGi,PVQU,PVQL,RealPD,RealQD,QD,PD);
+
+Volt=0.99*Volt;
+
 Gap=(Init_L*Init_Z'-Init_U*Init_W');
 KK=0;
 plotGap=zeros(1,60);
@@ -112,6 +115,7 @@ Vbi=sparse(0.5*ones(Busnum,1));
 PDbi=sparse(0.5*ones(length(Loadi),1));
 QDbi=sparse(0.5*ones(length(Loadi),1));
 eps=1;
+
 % 第一遍，算连续的值
 fprintf('第1次迭代，算连续量。\n');
 while(abs(Gap)>Precision*10)
@@ -140,7 +144,7 @@ while(abs(Gap)>Precision*10)
     %% 形成方程矩阵
     Luu=Init_U'.*Init_W'+Init_u*ones(RestraintCount,1);
     Lul=Init_L'.*Init_Z'-Init_u*ones(RestraintCount,1);
-    bigM=0.4;
+    bigM=0.5;
     Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi,QDbi,mPD,mQD,uPD,lPD,uQD,lQD,uVolt,lVolt);
     Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi);
     Ly=Mat_H;
@@ -165,6 +169,12 @@ Gap=1;
 % eps=0.00001;
 fprintf('\n');
 
+Init_Z=sparse(ones(1,RestraintCount));
+Init_W=sparse(-1*ones(1,RestraintCount));
+Init_L=1*sparse(ones(1,RestraintCount));
+Init_U=1*sparse(ones(1,RestraintCount));
+Init_Y=sparse(1,2*Busnum);%与学姐一致
+
 fprintf('第2次迭代，算离散量。\n');
 while(abs(Gap)>Precision*10)
     if KK>kmax
@@ -193,8 +203,8 @@ while(abs(Gap)>Precision*10)
     %% 形成方程矩阵
     Luu=Init_U'.*Init_W'+Init_u*ones(RestraintCount,1);
     Lul=Init_L'.*Init_Z'-Init_u*ones(RestraintCount,1);
-    bigM=0.4;
-    eps=Gap*0.001;
+    bigM=0.5;
+    eps=Gap*0.1;
     Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi,QDbi,mPD,mQD,uPD,lPD,uQD,lQD,uVolt,lVolt);
     Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi);
     Ly=Mat_H;