From 0e24eed4bfd635124277aff8336960c13c72e7fd Mon Sep 17 00:00:00 2001
From: "dugg@lab-desk" <dugg@lab-desk>
Date: Sun, 20 Jul 2014 21:08:16 +0800
Subject: [PATCH] =?UTF-8?q?1.=E4=BF=AE=E5=A4=8D=E4=BA=86=E8=BF=99=E4=B8=AA?=
 =?UTF-8?q?=E5=9C=B0=E6=96=B9=E7=9A=84bug=20lPD=3Dabs(RealPD*3*sigma);=20u?=
 =?UTF-8?q?PD=3Dabs(RealPD*3*sigma);=20lQD=3Dabs(RealQD*3*sigma);=20uQD=3D?=
 =?UTF-8?q?abs(RealQD*3*sigma);=20lVolt=3Dabs(rVolt'*3*sigma);=20uVolt=3Da?=
 =?UTF-8?q?bs(rVolt'*3*sigma);=202.=E5=8A=A0=E5=85=A5=E4=BA=86=E8=80=83?=
 =?UTF-8?q?=E8=99=91=E5=8F=91=E7=94=B5=E6=9C=BA=E7=9A=84=E6=83=85=E5=86=B5?=
 =?UTF-8?q?=EF=BC=8C=E5=85=B6=E5=AE=9E=E4=B9=9F=E5=B0=B1=E6=98=AF=E5=8A=A0?=
 =?UTF-8?q?=E5=85=A5=E8=B4=9F=E6=95=B0=E8=B4=9F=E8=8D=B7=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 | 20 +++++++++++---------
 1 file changed, 11 insertions(+), 9 deletions(-)

diff --git a/OPF.m b/OPF.m
index 6d89747..d5419ff 100644
--- a/OPF.m
+++ b/OPF.m
@@ -6,7 +6,7 @@ clear
 %%
 thesis=ForThesis(1,62);
 [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:\算例\东际911_2751267_2012-09-05\newFIle9-2.txt');
+    pf('E:\算例\17\17.csv');
 %         pf('E:\算例\柳金Ⅰ926_21671693_2012-09-06\newFIle16.txt');
 %pf('D:\Project\青秀降损项目\最小化潮流\最小潮流算例\原始\津头站津视922（3-1）_0.5_120%.txt');
 %pf('D:\Project\最小化潮流\最小潮流算例\仙海919.txt');
@@ -56,12 +56,12 @@ mVolt=rVolt.*(1+normrnd(0,sigma,length(rVolt),1))';
 mPD=PD0;
 mQD=QD0;
 %加估计上下界
-lPD=RealPD*(1-3*sigma);
-uPD=RealPD*(1+3*sigma);
-lQD=RealQD*(1-3*sigma);
-uQD=RealQD*(1+3*sigma);
-lVolt=rVolt'*(1-3*sigma);
-uVolt=rVolt'*(1+3*sigma);
+lPD=abs(RealPD*3*sigma);
+uPD=abs(RealPD*3*sigma);
+lQD=abs(RealQD*3*sigma);
+uQD=abs(RealQD*3*sigma);
+lVolt=abs(rVolt'*3*sigma);
+uVolt=abs(rVolt'*3*sigma);
 %错误数据
 %mVolt(2)=5;
 bigM=10;
@@ -148,10 +148,11 @@ while(abs(Gap)>Precision*1)
     Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi);
     Ly=Mat_H;
     if KK>8
-        eps=eps*0.2;
+        eps=eps*0.3;
         if abs(eps)<1e-6
-            eps=1e-6;
+            eps=1e-5;
         end
+        eps;
     end
     Lz=FormLz(Mat_G,Init_L,GenL,Busnum,PVQL,PD,RealPD,RealQD,Loadi,KK,PF,eps);
     Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD,PD0,QD0,Loadi,KK,PF,eps);
@@ -167,6 +168,7 @@ while(abs(Gap)>Precision*1)
     KK=KK+1;
 end
 toc
+eps
 obj=sum(Vbi)+sum(PDbi)+sum(QDbi);
 fprintf('目标函数值 %.2f\n',full(obj));