From 6e6e5e34d872d84fb0c04899ec321ba2c5943c9d Mon Sep 17 00:00:00 2001
From: unknown <dmy@dmy-PC.(none)>
Date: Fri, 16 Nov 2012 22:23:43 +0800
Subject: [PATCH] =?UTF-8?q?=E7=94=A8=E5=8A=9F=E7=8E=87=E5=9B=A0=E6=95=B0?=
 =?UTF-8?q?=E5=AE=9AQD=E7=9A=84=E4=B8=8A=E9=99=90=E3=80=82?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: unknown <dmy@dmy-PC.(none)>
---
 FormLw.m |  3 ++-
 OPF.m    | 11 ++++++++---
 2 files changed, 10 insertions(+), 4 deletions(-)

diff --git a/FormLw.m b/FormLw.m
index d2778a2..d2dc593 100644
--- a/FormLw.m
+++ b/FormLw.m
@@ -1,4 +1,4 @@
-function Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD0,QD0,Loadi,KK)
+function Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD,PD0,QD0,Loadi,KK,PF)
 KK=999;
 %PU=GenU(:,2);%发电机有功上界
 PU=5*ones(length(GenU(:,2)),1);
@@ -15,6 +15,7 @@ QDU=QD0(Loadi);
 QDU(QDU>0)=1.200*QDU(QDU>0);
 QDU(QDU<0)=0.800*QDU(QDU<0);
 QDU(QDU==0)=0.200;
+%QDU=1.2*PD(Loadi).*sqrt(1 -PF.^2)./PF;
 t1=([PU',QU',PDU',QDU',VoltU])';
 t2=Mat_G+Init_U'-t1;
 Lw=t2;
diff --git a/OPF.m b/OPF.m
index 6ba6041..68a0ab3 100644
--- a/OPF.m
+++ b/OPF.m
@@ -2,14 +2,16 @@ tic
 clc
 clear
 [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,Transfork0]= ...
-pf('ieee118.dat');
+pf('ieee30.dat');
 %pf('D:\Project\青秀降损项目\最小化潮流\最小潮流算例\原始\津头站津视922（3-1）_0.5_120%.txt'); 
 %pf('D:\Project\最小化潮流\最小潮流算例\仙海919.txt'); 
 %pf('c:/file31.txt'); 
 
 
 %% 计算功率因数
-%atan(PD(QD~=0 | PD~=0)./QD(QD~=0 | PD~=0));
+Loadi=QD~=0 | PD~=0;
+PF=sqrt(PD(Loadi).^2./(QD(Loadi).^2+PD(Loadi).^2))
+%% 
 Volt;
 UAngel*180/3.1415926;
 %% 通过潮流计算PG
@@ -22,6 +24,7 @@ QG0=QG;
 PD0=PD;
 QD0=QD;
 PDReal=PD;%真值
+QDReal=QD;%真值
 %PD0(12)=PD0(12)+0.001;
 %%
 PG0(Balance)=PGBal(Balance);
@@ -71,7 +74,7 @@ while(abs(Gap)>Precision)
     Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,QD_NON_ZERO,QD_NON_ZERO_IND,Loadi);
     Ly=Mat_H;
     Lz=FormLz(Mat_G,Init_L,GenL,Busnum,PVQL,PD0,QD0,Loadi,KK);
-    Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD0,QD0,Loadi,KK);
+    Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD,PD0,QD0,Loadi,KK,PF);
     Lx=FormLx(deltF,deltH,Init_Y,deltG,Init_Z,Init_W);
     YY=FormYY(Lul,Lz,Ly,Luu,Lw,Lx);
     %% 开始解方程
@@ -93,6 +96,8 @@ DrawGap(plotGap);
 % absPDLoad=abs(  (PD(Loadi)-PD0(Loadi))./PD0(Loadi)  );
 absPDLoad=abs(  (PD(Loadi)-PDReal(Loadi))./PDReal(Loadi)  );
 maxPDError=max(absPDLoad)
+absQDLoad=abs(  (QD(Loadi)-QDReal(Loadi))./QDReal(Loadi)  );
+maxQDError=max(absQDLoad)
 disp('index');
 Loadi(absPDLoad==maxPDError);
 %% 计算总线损