diff --git a/run.m b/run.m index 8c4f9fa..bfa9e9d 100644 --- a/run.m +++ b/run.m @@ -1,4 +1,4 @@ -function [JMeasurement,SEMeasurement,SEMeasurementNOPQ,SEMeasurementM,SEMeasurementNOPQM,isConverged]=run() +function [JMeasurement,AME_Volt,AME_VAngle,isConverged]=run() %% 利用先把负荷转换为电流的方法。这个方法要求知道电压量。 % close all @@ -407,18 +407,8 @@ JMeasurement=sum(sum((( abs(mVoltABCV)-abs(SEVoltpABC) )./ abs(mVoltABCV)./sigma Busnum=busNum; mCount=Busnum*3+length(Loadi)*3*2; %估计量质量 -%有PD,QD的 -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)); -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; +AME_Volt=sum(abs( abs(mVoltABCV)-abs(SEVoltpABC))); +AME_VAngle=sum(abs( angle(mVoltABCV)-angle(SEVoltpABC))); %返回收敛信息 isConverged=isConverged1*isConverged2*isConverged0; end \ No newline at end of file diff --git a/test.m b/test.m index b34017e..c6e1982 100644 --- a/test.m +++ b/test.m @@ -1,27 +1,19 @@ close all clear clc -JMeasurementSum=0; -SEMeasurementSum=0; -SEMeasurementNOPQSum=0; -SEMeasurementMSum=0; -SEMeasurementNOPQMSum=0; +AME_VoltSum=0; +AME_VAngleSum=0; N=1000; for I=1:N - [JMeasurement,SEMeasurement,SEMeasurementNOPQ,SEMeasurementM,SEMeasurementNOPQM,isConverged]=run(); + [JMeasurement,AME_Volt,AME_VAngle,isConverged]=run(); if isConverged==0 I=I-1; continue end % [JMeasurement]=run(); - JMeasurementSum=JMeasurement+JMeasurementSum; - SEMeasurementSum=SEMeasurement+SEMeasurementSum; - SEMeasurementNOPQSum=SEMeasurementNOPQ+SEMeasurementNOPQSum; - SEMeasurementMSum=SEMeasurementMSum+SEMeasurementM; - SEMeasurementNOPQMSum=SEMeasurementNOPQMSum+SEMeasurementNOPQM; + AME_VoltSum=AME_VoltSum+AME_Volt; + AME_VAngleSum=AME_VAngleSum+AME_VAngle; end JMeasurementSum=JMeasurementSum/N; -SEMeasurementSum=SEMeasurementSum/N; -SEMeasurementNOPQSum=SEMeasurementNOPQSum/N; -SEMeasurementMSum=SEMeasurementMSum/N; -SEMeasurementNOPQMSum=SEMeasurementNOPQMSum/N; \ No newline at end of file +AME_VoltSum=AME_VoltSum/N; +AME_VAngleSum=AME_VAngleSum/N;