26 lines
889 B
Matlab
26 lines
889 B
Matlab
function [ output_args ] = checkSSatisfied(Balance,phaseABCY,VoltpABC,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP,phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ )
|
||
%%利用Fortiscue方法计算得到了相电压和相角后，反推负荷功率，检查是否和给定的功率一致。
|
||
busNum=size(phaseABCY,1)/3;
|
||
VoltpA=VoltpABC(1,:);
|
||
VoltpB=VoltpABC(2,:);
|
||
VoltpC=VoltpABC(3,:);
|
||
Vp3=sparse(busNum*3,1);
|
||
Vp3(1:3:end)=VoltpA;
|
||
Vp3(2:3:end)=VoltpB;
|
||
Vp3(3:3:end)=VoltpC;
|
||
Ip3=phaseABCY*Vp3;
|
||
Sp3=Vp3.*conj(Ip3);
|
||
%VoltpABC.*conj(IpABC);
|
||
pLoadABC=sparse(length(phaseASpotLoadP)*3,1);
|
||
pLoadABC(1:3:end)=phaseASpotLoadP+1j*phaseASpotLoadQ;
|
||
pLoadABC(2:3:end)=phaseBSpotLoadP+1j*phaseBSpotLoadQ;
|
||
pLoadABC(3:3:end)=phaseCSpotLoadP+1j*phaseCSpotLoadQ;
|
||
ck=Sp3+pLoadABC;
|
||
ck(3*(Balance-1)+1:3*(Balance-1)+3)=0;
|
||
if any(real(ck)>1e-5) ||any(imag(ck)>1e-5)
|
||
fprintf('反推回的功率不匹配。潮流方程约束不满足。\n');
|
||
end
|
||
output_args=ck;
|
||
end
|
||
|