facat
/
pes2014-debug2-laplace
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

增加计算总的损耗. 

Signed-off-by: unknown <dmy@dmy-PC.(none)>
This commit is contained in:
unknown 2012-08-01 10:32:46 +08:00
parent 5e0f4c2520
commit aee15a5178
2 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
OPF.m
View File

@ -88,7 +88,10 @@ absPDLoad=abs( (PD(Loadi)-PDReal(Loadi))./PDReal(Loadi) );
maxPDError=max(absPDLoad);
disp('index');
Loadi(absPDLoad==maxPDError);
%% 计算线损
%% 计算总线损
totalLoss=(PGBal(Balance)-sum(PD(Loadi)))*100;
fprintf('总的损耗为%f(MW 有名值)\n',totalLoss);
%% 计算各线损
Lineloss(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Transfork0,Volt,UAngel);
toc

2
OPF_Init.m
View File

@ -1,5 +1,5 @@
function [Volt,UAngel,Init_Z,Init_W,Init_L,Init_U,Init_Y,PG,QG,RestraintCount,wG,wD,PD,PD0,randPDind,Loadi]=OPF_Init(Busnum,Balance,PG,QG,Volt,GenU,GenL,PVi,PGi,PVQU,PVQL,PD0,QD,PD)
Loadi=find(QD~=0 & PD~=0);
Loadi=find(QD~=0 | PD~=0);
%Loadi=[1:Busnum]';
RestraintCount=size(PVi,1)+size(PGi,1)+size(Loadi,1)+Busnum*1; %约束条件数
t_Bal_volt=Volt(Balance);