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

采用新的电压幅值和相角Absolute Mean Error 

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

16
run.m
View File

@ -1,4 +1,4 @@
function [JMeasurement,SEMeasurement,SEMeasurementNOPQ,SEMeasurementM,SEMeasurementNOPQM,isConverged]=run() function [JMeasurement,AME_Volt,AME_VAngle,isConverged]=run()
%% 利用先把负荷转换为电流的方法。这个方法要求知道电压量。 %% 利用先把负荷转换为电流的方法。这个方法要求知道电压量。
% %
close all close all
@ -407,18 +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;
%估计量质量 %估计量质量
%有PD，QD的 AME_Volt=sum(abs( abs(mVoltABCV)-abs(SEVoltpABC)));
SEMeasurement=sum(sum((( abs(rVoltABCV)-abs(SEVoltpABC) )./abs(mVoltABCV)./sigma).^2))+sum(sum(((rPD3P-PD3P)./mPD3P./sigma).^2))+sum(sum(((rQD3P-QD3P)./mQD3P./sigma).^2)); AME_VAngle=sum(abs( angle(mVoltABCV)-angle(SEVoltpABC)));
SEMeasurement=(SEMeasurement/mCount)^.5;
%没有PD，QD的
SEMeasurementNOPQ=sum(sum(((abs(rVoltABCV)-abs(SEVoltpABC))./abs(mVoltABCV)./sigma).^2));
SEMeasurementNOPQ=(SEMeasurementNOPQ/mCount)^.5;
%量测量质量
SEMeasurementM=sum(sum((( abs(rVoltABCV)-abs(mVoltABCV) )./abs(mVoltABCV)./sigma).^2))+sum(sum(((rPD3P-mPD3P)./mPD3P./sigma).^2))+sum(sum(((rQD3P-mQD3P)./mQD3P./sigma).^2));
SEMeasurementM=(SEMeasurementM/mCount)^.5;
%没有PD，QD的
SEMeasurementNOPQM=sum(sum((( abs(rVoltABCV)-abs(mVoltABCV))./abs(mVoltABCV)./sigma).^2));
SEMeasurementNOPQM=(SEMeasurementNOPQM/mCount)^.5;
%返回收敛信息 %返回收敛信息
isConverged=isConverged1*isConverged2*isConverged0; isConverged=isConverged1*isConverged2*isConverged0;
end end

22
test.m
View File

@ -1,27 +1,19 @@
close all close all
clear clear
clc clc
JMeasurementSum=0; AME_VoltSum=0;
SEMeasurementSum=0; AME_VAngleSum=0;
SEMeasurementNOPQSum=0;
SEMeasurementMSum=0;
SEMeasurementNOPQMSum=0;
N=1000; N=1000;
for I=1:N for I=1:N
[JMeasurement,SEMeasurement,SEMeasurementNOPQ,SEMeasurementM,SEMeasurementNOPQM,isConverged]=run(); [JMeasurement,AME_Volt,AME_VAngle,isConverged]=run();
if isConverged==0 if isConverged==0
I=I-1; I=I-1;
continue continue
end end
% [JMeasurement]=run(); % [JMeasurement]=run();
JMeasurementSum=JMeasurement+JMeasurementSum; AME_VoltSum=AME_VoltSum+AME_Volt;
SEMeasurementSum=SEMeasurement+SEMeasurementSum; AME_VAngleSum=AME_VAngleSum+AME_VAngle;
SEMeasurementNOPQSum=SEMeasurementNOPQ+SEMeasurementNOPQSum;
SEMeasurementMSum=SEMeasurementMSum+SEMeasurementM;
SEMeasurementNOPQMSum=SEMeasurementNOPQMSum+SEMeasurementNOPQM;
end end
JMeasurementSum=JMeasurementSum/N; JMeasurementSum=JMeasurementSum/N;
SEMeasurementSum=SEMeasurementSum/N; AME_VoltSum=AME_VoltSum/N;
SEMeasurementNOPQSum=SEMeasurementNOPQSum/N; AME_VAngleSum=AME_VAngleSum/N;
SEMeasurementMSum=SEMeasurementMSum/N;
SEMeasurementNOPQMSum=SEMeasurementNOPQMSum/N;