facat/stateestimateionyalmip-lu9-forcomparision
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

通过潮流方程把PD,QD反求出来。

Browse Source 
Signed-off-by: dmy@lab <dmy@lab.lab>
This commit is contained in:
dmy@lab
2015-03-19 13:40:41 +08:00
parent b7498c2b15
commit 2dd7ffb183
1 changed files with 5 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
Run_YALMIP2.m
View File

@@ -112,8 +112,11 @@ for I=1:1
SEVolt=x(1:length(Volt)); SEVolt=x(1:length(Volt));
SEVAngel=x(length(Volt)+1:end); SEVAngel=x(length(Volt)+1:end);
fprintf('最大偏差\n') fprintf('最大偏差\n')
PD=x(1:length(PDi)); AngleIJ=sparse(r,c,SEVAngel(r)-SEVAngel(c)-Angle,Busnum,Busnum);
QD=x(length(PDi)+1:length(PDi)+length(QDi)); zeroP=diag(SEVolt)*Y.*cos(AngleIJ)*SEVolt;
zeroQ=diag(SEVolt)*Y.*sin(AngleIJ)*SEVolt;
PD=-zeroP(PDi);
QD=-zeroQ(QDi);
maxDev1=MaxDeviation(sigma,mVolt,SEVolt,rVAngel,SEVAngel,rPD(PDi),rQD(QDi),PD,QD) maxDev1=MaxDeviation(sigma,mVolt,SEVolt,rVAngel,SEVAngel,rPD(PDi),rQD(QDi),PD,QD)
fprintf('统计偏差\n') fprintf('统计偏差\n')
statDev1=StatDeviation(sigma,rVolt,SEVolt,rVAngel,SEVAngel,rPD(PDi),rQD(QDi),PD,QD) statDev1=StatDeviation(sigma,rVolt,SEVolt,rVAngel,SEVAngel,rPD(PDi),rQD(QDi),PD,QD)