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

删掉一些注释，先吧bar图形隐掉。 

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk 2014-11-05 12:06:13 +08:00
parent 3875907756
commit 77dd602671
1 changed files with 5 additions and 29 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

34
run.m
View File

@ -78,26 +78,6 @@ while(k<=kmax && maxD> EPS)
iterPhaseASpotLoadQ=phaseASpotLoadQ+imag(SA); iterPhaseASpotLoadQ=phaseASpotLoadQ+imag(SA);
iterPhaseBSpotLoadQ=phaseBSpotLoadQ+imag(SB); iterPhaseBSpotLoadQ=phaseBSpotLoadQ+imag(SB);
iterPhaseCSpotLoadQ=phaseCSpotLoadQ+imag(SC); iterPhaseCSpotLoadQ=phaseCSpotLoadQ+imag(SC);
% %全部转换为负荷电流
% CurpA=conj((iterPhaseASpotLoadP+1j*iterPhaseASpotLoadQ)./VoltpA);
% CurpB=conj((iterPhaseBSpotLoadP+1j*iterPhaseBSpotLoadQ)./VoltpB);
% CurpC=conj((iterPhaseCSpotLoadP+1j*iterPhaseCSpotLoadQ)./VoltpC);
% %转换为序电流
% f012=Tp2f*conj([CurpA';CurpB';CurpC']);
% %把三序电流分离出来
% If0=conj(f012(1,:)');
% If1=conj(f012(2,:)');
% If2=conj(f012(3,:)');
% %试着算一下正序电流
% fsY11*V1;
% %形成负荷序电流的测量值
% mIf0=If0;
% mIf1=If1;
% mIf1(3)=-mIf1(2);
% mIf2=If2;
% %计算
% fsY11=fsY11+sparse(Balance,Balance,ones(length(Balance),1),busNum,bus
% Num);%这里要置0，置1，否则是奇异的
%%做最小二乘法 %%做最小二乘法
[dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance,busNum, ... [dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance,busNum, ...
PQi,PG,QG,QGi,iterPD,iterQD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量 PQi,PG,QG,QGi,iterPD,iterQD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量
@ -118,8 +98,6 @@ while(k<=kmax && maxD> EPS)
If0=conj(f012(1,:)'); If0=conj(f012(1,:)');
If1=conj(f012(2,:)'); If1=conj(f012(2,:)');
If2=conj(f012(3,:)'); If2=conj(f012(3,:)');
If0(Balance)=0; If0(Balance)=0;
If2(Balance)=0; If2(Balance)=0;
%Vf0=fsY0\If0; %Vf0=fsY0\If0;
@ -127,8 +105,6 @@ while(k<=kmax && maxD> EPS)
%Vf2=fsY2\If2; %Vf2=fsY2\If2;
Vf2=fsY2Q*(fsY2U\(fsY2L\(fsY2P*(fsY2R\If2)))); Vf2=fsY2Q*(fsY2U\(fsY2L\(fsY2P*(fsY2R\If2))));
fprintf('迭代时间%f\n',toc); fprintf('迭代时间%f\n',toc);
%
end end
FortiscueToc=toc; FortiscueToc=toc;
fprintf('Fortiscue法计算时间 %f\n',FortiscueToc); fprintf('Fortiscue法计算时间 %f\n',FortiscueToc);
@ -277,7 +253,7 @@ SEIpABC=full(Tf2p*conj([SEIf0';SEIf1(1:3)';SEIf2']));
%看一下差多少 %看一下差多少
VError=(full(abs(VoltpABC))-abs(SEVoltpABC))./abs(VoltpABC)*100; VError=(full(abs(VoltpABC))-abs(SEVoltpABC))./abs(VoltpABC)*100;
VError=reshape(VError,size(VError,1)*size(VError,2),1); VError=reshape(VError,size(VError,1)*size(VError,2),1);
barPlot( VError, 10,['相对误差%'],['分布密度'],['电压幅值']); % barPlot( VError, 10,['相对误差%'],['分布密度'],['电压幅值']);
% 三相负荷 % 三相负荷
rThreeLoad=[ phaseASpotLoadP'+1j*phaseASpotLoadQ'; rThreeLoad=[ phaseASpotLoadP'+1j*phaseASpotLoadQ';
phaseBSpotLoadP'+1j*phaseBSpotLoadQ'; phaseBSpotLoadP'+1j*phaseBSpotLoadQ';
@ -289,7 +265,7 @@ 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()
barPlot( phaseLoadQError, 10,'相对误差%','分布密度','无功负荷误差'); % barPlot( phaseLoadQError, 10,'相对误差%','分布密度','无功负荷误差');