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

不等式约束是这样的：

Browse Source 
sparse(PD(Loadi));
    sparse(PD(Loadi))-0.001;
    sparse(QD(Loadi));
    sparse(QD(Loadi))-0.001;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk
2014-05-21 15:35:43 +08:00
parent 1ee3b9312c
commit 9ef19c4542
11 changed files with 87 additions and 37 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
OPF_Init.m
View File

@@ -1,7 +1,7 @@
function [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,PD0,QD0,QD,PD)
Loadi=find(QD~=0 | PD~=0);
%Loadi=[1:Busnum]';
RestraintCount=size(Loadi,1)*2+Busnum*2+Busnum; %约束条件数,放开所有QD
RestraintCount=size(Loadi,1)*4+Busnum*2+Busnum+length(Loadi)*0; %约束条件数,放开所有QD
t_Bal_volt=Volt(Balance);
Volt=sparse(1*ones(1,Busnum));
Volt(Balance)=t_Bal_volt;