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

有把互补约束全部加上，收敛次数18~19次。

Browse Source 
Mat_G=[
    sparse(PD(Loadi))-mPD(Loadi)-bigM*PDbi-0.1;
    sparse(PD(Loadi))-mPD(Loadi)+bigM*PDbi+0.1;
    sparse(QD(Loadi))-mQD(Loadi)-bigM*QDbi-0.1;
    sparse(QD(Loadi))-mQD(Loadi)+bigM*QDbi+0.1;
    Volt'-mVolt'-bigM*Vbi-0.03;
    Volt'-mVolt'+bigM*Vbi+0.03;
    Vbi;
    PDbi;
    QDbi;
    (1-Vbi).*Vbi;
   (1-PDbi).*PDbi;
    (1-QDbi).*QDbi;
    ];

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk
2014-05-25 11:35:26 +08:00
parent 721035b580
commit 240f989115
7 changed files with 14 additions and 17 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
FormH.m
View File

@@ -2,7 +2,5 @@ function Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi)
AngleIJ=sparse(r,c,UAngel(r)-UAngel(c)-Angle',Busnum,Busnum);
dP=PG-PD-diag(Volt)*Y.*cos(AngleIJ)*Volt';
dQ=QG-QD-diag(Volt)*Y.*sin(AngleIJ)*Volt';
Mat_H=[dP;dQ;];
end