From 2237f178a9cfd0eb9037594891b0781e0b002598 Mon Sep 17 00:00:00 2001
From: "dugg@lab-desk" <dugg@lab-desk>
Date: Wed, 4 Feb 2015 22:16:21 +0800
Subject: [PATCH] =?UTF-8?q?=E6=8C=89=E6=97=B6=E4=BF=9D=E5=AD=98=EF=BC=8C?=
 =?UTF-8?q?=E5=A5=BD=E4=B9=A0=E6=83=AF=E3=80=82=E8=BF=98=E5=9C=A8=E5=81=9A?=
 =?UTF-8?q?=E6=B5=8B=E8=AF=95=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>
---
 OPF.m    | 4 ++--
 subOPF.m | 6 ++++--
 2 files changed, 6 insertions(+), 4 deletions(-)

diff --git a/OPF.m b/OPF.m
index 28913b0..479b3b9 100644
--- a/OPF.m
+++ b/OPF.m
@@ -10,7 +10,7 @@ badDataLocation=zeros(34,10);
 maxDVolt=0;
 maxDVAngle=0;
 % for badDataNode=1:33
-    for I=1:5
+    for I=1:1
         close
         [kmax,Precision,UAngel,Volt,Busnum,PVi,PVu,Balance,Y,Angle,P0,Q0,r,c,GB,Linei,Linej,Transfori,Transforj,GenU,GenL,GenC,PG,QG,PD,QD,CenterA,PGi,PVQU,PVQL,Liner,Linex,Lineb,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0]= ...
             pf('E:\算例\feeder33\feeder33.txt');
@@ -33,7 +33,7 @@ maxDVAngle=0;
         % mVolt=load('mVolt');
         % mVolt=mVolt.mVolt;
         
-        mVolt(18)=rVolt(18)*(1-sigma*6);
+        mVolt(31)=rVolt(31)*(1-sigma*6);
         %% 画Case A的图
         % figure('Color',[1 1 1]);
         [Busnum,Loadi,Volt,PD,QD,rVolt,UAngel,RealPD,RealQD,rUAngel,Vbi,PDbi,QDbi,plotGapA]=subOPF([],PD0,QD0,mVolt,sigma);%全部有
diff --git a/subOPF.m b/subOPF.m
index 25e78c7..41656f0 100644
--- a/subOPF.m
+++ b/subOPF.m
@@ -45,7 +45,7 @@ Gap=(Init_L*Init_Z'-Init_U*Init_W');
 KK=0;
 plotGap=zeros(1,60);
 ContrlCount=size(Loadi,1)*2+Busnum*2+Busnum+length(Loadi)*2;
-kmax=500;
+kmax=100;
 Precision=Precision/1;
 %% 加误差
 %找DG
@@ -124,7 +124,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.6;
+    bigM=0.8;
     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;
@@ -154,6 +154,7 @@ while(abs(Gap)>Precision*10)
     if KK>kmax
                     break;
     end
+%             Vbi(23)=0.1;
     plotGap(KK+1)=Gap;
     Init_u=Gap/2/RestraintCount*CenterA;
     AngleIJMat=0;
@@ -188,6 +189,7 @@ while(abs(Gap)>Precision*10)
     %% 开始解方程
 %         Vbi(Vbi>0.003)=.1;
     fprintf('迭代次数 %d Gap %f\n',KK+1,plotGap(KK+1));
+
     XX=SolveIt(deltF,deltG,Init_L,Init_Z,Init_U,Init_W,deltdeltF,ddh,ddgzw,deltH,Init_Y,Ly,Lz,ContrlCount,Lw,Lul,Luu,RestraintCount,Lx,Balance,PVi,PGi,Busnum,Loadi);
     [deltZ,deltL,deltW,deltU,deltX,deltY]=AssignXX(XX,ContrlCount,RestraintCount,Busnum);
     [Init_Z,Init_L,Init_W,Init_U,Init_Y,PG,QG,Volt,UAngel,PD,QD,Vbi,PDbi,QDbi]=Modification(Init_Z,Init_L,Init_W,Init_U,Init_Y,deltZ,deltL,deltW,deltU,deltX,deltY,PG,QG,Volt,UAngel,PVi,ContrlCount,Balance,Busnum,PGi,PD,QD,Loadi,Vbi,PDbi,QDbi);