distributionnetwork-power2c.../test.m

80 lines
2.1 KiB
Matlab

close all
clear
clc
AME_VoltSum=0;
AME_VAngleSum=0;
JMeasurementSum=0;
AME_PDSum=0;
AME_QDSum=0;
AME_mVoltSum=0;
AME_mPDSum=0;
AME_mQDSum=0;
totalTime1Sum=0;
totalTime2Sum=0;
totalTime0Sum=0;
timeOfPrepare=0.00272333333;%×¼±¸Êý¾ÝµÄʱ¼ä
timeOfPrepare=0;
eV2Sum=0;
eV0Sum=0;
eI2Sum=0;
eI0Sum=0;
interation1Sum=0;
interation2Sum=0;
interation0Sum=0;
N=2000;
loopN=1;
while 1
[JMeasurement,AME_Volt,AME_VAngle,AME_PD,AME_QD,AME_mVolt,AME_mPD,AME_mQD,isConverged,totalTime1,totalTime2,totalTime0,eV2,eV0,eI2,eI0,interation1,interation2,interation0]=run();
if isConverged==0
continue
end
loopN=loopN+1;
if loopN>N
break;
end
% [JMeasurement]=run();
AME_VoltSum=AME_VoltSum+AME_Volt;
AME_VAngleSum=AME_VAngleSum+AME_VAngle;
AME_PDSum=AME_PDSum+AME_PD;
AME_QDSum=AME_QDSum+AME_QD;
AME_mVoltSum=AME_mVoltSum+AME_mVolt;
AME_mPDSum=AME_mPDSum+AME_mPD;
AME_mQDSum=AME_mQDSum+AME_mQD;
eV2Sum=eV2Sum+eV2;
eV0Sum=eV0Sum+eV0;
eI2Sum=eI2Sum+eI2;
eI0Sum=eI0Sum+eI0;
interation1Sum=interation1Sum+interation1;
interation2Sum=interation2Sum+interation2;
interation0Sum=interation0Sum+interation0;
%ʱ¼ä
totalTime1Sum=totalTime1Sum+totalTime1;
totalTime2Sum=totalTime2Sum+totalTime2;
totalTime0Sum=totalTime0Sum+totalTime0;
% totalTime1Sum=min([totalTime1Sum,totalTime1]);
% totalTime2Sum=min([totalTime2Sum,totalTime2]);
% totalTime0Sum=min([totalTime0Sum,totalTime0]);
end
JMeasurementSum=JMeasurementSum/N;
AME_VoltSum=AME_VoltSum/N;
AME_VAngleSum=AME_VAngleSum/N;
AME_PDSum=AME_PDSum/N;
AME_QDSum=AME_QDSum/N;
AME_mVoltSum=AME_mVoltSum/N;
AME_mPDSum=AME_mPDSum/N;
AME_mQDSum=AME_mQDSum/N;
eV2Sum=eV2Sum/N;
eV0Sum=eV0Sum/N;
eI2Sum=eI2Sum/N;
eI0Sum=eI0Sum/N;
% totalTime1Sum=totalTime1Sum+timeOfPrepare;
% totalTime2Sum=totalTime2Sum+timeOfPrepare;
% totalTime0Sum=totalTime0Sum+timeOfPrepare;
totalTime1Sum=totalTime1Sum/N;
totalTime2Sum=totalTime2Sum/N;
totalTime0Sum=totalTime0Sum/N;
interation1Sum=interation1Sum/N;
interation2Sum=interation2Sum/N;
interation0Sum=interation0Sum/N;