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

把测试文件的代码改好了。 

Signed-off-by: facat@ipso <facat@ipso>
This commit is contained in:
facat@ipso 2014-11-09 23:28:36 +08:00
parent 3a1655112a
commit 0f4d416ca0
2 changed files with 10 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
run.m
View File

@ -1,3 +1,4 @@
function [JMeasurement,SEMeasurement,SEMeasurementNOPQ,SEMeasurementM,SEMeasurementNOPQM]=run()
%% 利用先把负荷转换为电流的方法。这个方法要求知道电压量。 %% 利用先把负荷转换为电流的方法。这个方法要求知道电压量。
% %
close all close all
@ -293,12 +294,12 @@ rThreeLoad=[ phaseASpotLoadP'+1j*phaseASpotLoadQ';
phaseBSpotLoadP'+1j*phaseBSpotLoadQ'; phaseBSpotLoadP'+1j*phaseBSpotLoadQ';
phaseCSpotLoadP'+1j*phaseCSpotLoadQ'; phaseCSpotLoadP'+1j*phaseCSpotLoadQ';
]; ];
rThreeLoad=rThreeLoad(:,setxor(1:size(SEVoltpABC,2),Balance)); % rThreeLoad=rThreeLoad(:,setxor(1:size(SEVoltpABC,2),Balance));
SEThreeLoad=SEVoltpABC(:,setxor(1:size(SEVoltpABC,2),Balance)).*conj(-SEIpABC); SEThreeLoad=SEVoltpABC(:,setxor(1:size(SEVoltpABC,2),Balance)).*conj(-SEIpABC);
phaseLoadPError=real(rThreeLoad-SEThreeLoad)./real(rThreeLoad)*100; % phaseLoadPError=real(rThreeLoad-SEThreeLoad)./real(rThreeLoad)*100;
phaseLoadQError=imag(rThreeLoad-SEThreeLoad)./imag(rThreeLoad)*100; % phaseLoadQError=imag(rThreeLoad-SEThreeLoad)./imag(rThreeLoad)*100;
phaseLoadPError=reshape(phaseLoadPError,size(phaseLoadPError,1)*size(phaseLoadPError,2),1); % phaseLoadPError=reshape(phaseLoadPError,size(phaseLoadPError,1)*size(phaseLoadPError,2),1);
phaseLoadQError=reshape(phaseLoadQError,size(phaseLoadQError,1)*size(phaseLoadQError,2),1); % phaseLoadQError=reshape(phaseLoadQError,size(phaseLoadQError,1)*size(phaseLoadQError,2),1);
% figure() % figure()
% barPlot( phaseLoadPError, 10,['相对误差%'],['分布密度'],['有功负荷误差']); % barPlot( phaseLoadPError, 10,['相对误差%'],['分布密度'],['有功负荷误差']);
% figure() % figure()
@ -346,4 +347,5 @@ SEMeasurementM=sum(sum((( abs(rVoltABCV)-abs(mVoltABCV) )./abs(mVoltABCV)./sigma
SEMeasurementM=(SEMeasurementM/mCount)^.5; SEMeasurementM=(SEMeasurementM/mCount)^.5;
%没有PD，QD的 %没有PD，QD的
SEMeasurementNOPQM=sum(sum((( abs(rVoltABCV)-abs(mVoltABCV))./abs(mVoltABCV)./sigma).^2)); SEMeasurementNOPQM=sum(sum((( abs(rVoltABCV)-abs(mVoltABCV))./abs(mVoltABCV)./sigma).^2));
SEMeasurementNOPQM=(SEMeasurementNOPQM/mCount)^.5; SEMeasurementNOPQM=(SEMeasurementNOPQM/mCount)^.5;
end

3
test.m
View File

@ -8,7 +8,8 @@ SEMeasurementMSum=0;
SEMeasurementNOPQMSum=0; SEMeasurementNOPQMSum=0;
N=1000; N=1000;
for I=1:N for I=1:N
[JMeasurement,SEMeasurement,SEMeasurementNOPQ,SEMeasurementM,SEMeasurementNOPQM]=OPF(); [JMeasurement,SEMeasurement,SEMeasurementNOPQ,SEMeasurementM,SEMeasurementNOPQM]=run();
% [JMeasurement]=run();
JMeasurementSum=JMeasurement+JMeasurementSum; JMeasurementSum=JMeasurement+JMeasurementSum;
SEMeasurementSum=SEMeasurement+SEMeasurementSum; SEMeasurementSum=SEMeasurement+SEMeasurementSum;
SEMeasurementNOPQSum=SEMeasurementNOPQ+SEMeasurementNOPQSum; SEMeasurementNOPQSum=SEMeasurementNOPQ+SEMeasurementNOPQSum;