From ef0954e86dad4c528c1af46eddc69c6b1f3a67c8 Mon Sep 17 00:00:00 2001
From: "facat@lab.com" <facat@lab.com>
Date: Mon, 3 Nov 2014 17:01:21 +0800
Subject: [PATCH] =?UTF-8?q?=E5=88=A0=E9=99=A4=E4=BA=86=E6=97=A0=E7=94=A8?=
 =?UTF-8?q?=E7=9A=84=E6=B3=A8=E9=87=8A?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: facat@lab.com <facat@lab.com>
---
 OPF.m | 50 --------------------------------------------------
 1 file changed, 50 deletions(-)

diff --git a/OPF.m b/OPF.m
index 0ff0b4d..cae9621 100644
--- a/OPF.m
+++ b/OPF.m
@@ -54,39 +54,6 @@ wPD=ones(length(Loadi),1);
 wQD=ones(length(Loadi),1);
 PD0=PD3P;
 QD0=QD3P;
-% orderedPDPG=sparse(PQi,1,PD-PG,busNum,1); %节点功率按顺序排列
-% orderedQDQG=sparse(PQi,1,QD-QG,busNum,1); %节点功率按顺序排列
-
-% [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:/算例/新民Ⅰ906_2729823_2012-09-06/newFIle20.txt');
-%% 计算功率因数
-
-%     PF=sqrt(PD(Loadi).^2./(QD(Loadi).^2+PD(Loadi).^2));
-%%
-% Volt;
-% UAngel*180/3.1415926;
-%% 通过潮流计算PG
-% AngleIJ=sparse(r,c,UAngel(r)-UAngel(c)-Angle',Busnum,Busnum);
-% PGBal=PD+diag(Volt)*Y.*cos(AngleIJ)*Volt';
-% QGBal=QD+diag(Volt)*Y.*sin(AngleIJ)*Volt';
-%% 初值-即测量值
-% PG0=PG;
-% QG0=QG;
-% PD0=PD;
-% QD0=QD;
-% Volt0=Volt;
-% UAngel0=UAngel;
-% %%
-% PG0(Balance)=PGBal(Balance);
-% PG(Balance)=PGBal(Balance);
-% QG0(Balance)=QGBal(Balance);
-% QG0(PVi)=QGBal(PVi);
-% QG(PVi)=QGBal(PVi);
-% %% 真实值
-% RealPG=PG0;
-% RealQG=QG0;
-% RealPD=PD0;
-% RealQD=QD0;
 %%
 RestraintCount=size(Loadi,1)*2; %约束条件数,放开所有QD
 Init_Z=sparse(ones(1,RestraintCount));
@@ -94,25 +61,14 @@ 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*3);%与学姐一致
-
-% [Volt,UAngel,Init_Z,Init_W,Init_L,Init_U,Init_Y,PG,QG,RestraintCount,wPG,wQG,wPD,wQD,PD,PD0,QD,randPDind,Loadi]=OPF_Init(Busnum,Balance,PG,QG,Volt,GenU,GenL,PVi,PGi,PVQU,PVQL,RealPD,RealQD,QD,PD);
 Gap=(Init_L*Init_Z'-Init_U*Init_W');
 KK=0;
-% plotGap=zeros(1,60);
 ContrlCount=size(Loadi,1)*2+Busnum*6;
 kmax=30;
-%% 20120523 临时
-% QD_NON_ZERO=QD(PD==0 & QD~=0);
-% QD_NON_ZERO_IND=find(PD==0 & QD~=0);
 %%
-% Precision=Precision/1;
 Precision=1e-5;
 CenterA=0.1;
 %% 加误差
-% PD0(Loadi)=PD0(Loadi).*(1+normrnd(0,0.05,length(Loadi),1));
-% QD0(Loadi)=QD0(Loadi).*(1+normrnd(0,0.05,length(Loadi),1));
-%% 读变压器容量
-% noDataTransCapacity=0;
 Volt=Vp3m;
 UAngel=Vp3a;
 while(abs(Gap)>Precision)
@@ -162,11 +118,5 @@ while(abs(Gap)>Precision)
     KK=KK+1;
 end
 fprintf('迭代次数%d\n',KK);
-% fprintf('目标值%f\n',full(ObjectiveFun(PG,PG0,PGi,QG,QG0,PVi,PD,PD0,QD,QD0,wPG,wQG,wPD,wQD,Loadi)));
-%% 计算总线损
-% totalLoss=(sum(PG)-sum(PD(Loadi)))*100;
-% fprintf('总的损耗为%f(MW 有名值)\n',full(totalLoss));
-% fprintf('线损率为 %f\n',full(totalLoss/sum(PG)));
-
 toc