facat
/
distributionnetwork-power2current
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

修复了计算AME的一个错误。 

Signed-off-by: dmy@lab <dmy@lab.lab>
This commit is contained in:
dmy@lab 2015-04-01 16:58:55 +08:00
parent 0e591ca2f9
commit 92c16fac04
2 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
run.m
View File

@ -407,8 +407,8 @@ JMeasurement=sum(sum((( abs(mVoltABCV)-abs(SEVoltpABC) )./ abs(mVoltABCV)./sigma
Busnum=busNum; Busnum=busNum;
mCount=Busnum*3+length(Loadi)*3*2; mCount=Busnum*3+length(Loadi)*3*2;
%估计量质量 %估计量质量
AME_Volt=sum(abs( abs(mVoltABCV)-abs(SEVoltpABC))); AME_Volt=sum(sum(abs( abs(mVoltABCV)-abs(SEVoltpABC))));
AME_VAngle=sum(abs( angle(mVoltABCV)-angle(SEVoltpABC))); AME_VAngle=sum(sum(abs( angle(mVoltABCV)-angle(SEVoltpABC))));
%返回收敛信息 %返回收敛信息
isConverged=isConverged1*isConverged2*isConverged0; isConverged=isConverged1*isConverged2*isConverged0;
end end

1
test.m
View File

@ -3,6 +3,7 @@ clear
clc clc
AME_VoltSum=0; AME_VoltSum=0;
AME_VAngleSum=0; AME_VAngleSum=0;
JMeasurementSum=0;
N=1000; N=1000;
for I=1:N for I=1:N
[JMeasurement,AME_Volt,AME_VAngle,isConverged]=run(); [JMeasurement,AME_Volt,AME_VAngle,isConverged]=run();