26 lines
889 B
Mathematica
26 lines
889 B
Mathematica
|
function [ output_args ] = checkSSatisfied(Balance,phaseABCY,VoltpABC,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP,phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ )
|
|||
|
|
%%<EFBFBD><EFBFBD><EFBFBD><EFBFBD>Fortiscue<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>õ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ѹ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>Ǻ<EFBFBD><EFBFBD>Ƹ<EFBFBD><EFBFBD>ɹ<EFBFBD><EFBFBD>ʣ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>Ƿ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>Ĺ<EFBFBD><EFBFBD><EFBFBD>һ<EFBFBD>¡<EFBFBD>
|
|||
|
|
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('<EFBFBD><EFBFBD><EFBFBD>ƻصĹ<EFBFBD><EFBFBD>ʲ<EFBFBD>ƥ<EFBFBD>䡣<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>Լ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>㡣\n');
|
|||
|
|
end
|
|||
|
|
output_args=ck;
|
|||
|
|
end
|
|||
|
|
|