stateestimation-ipm/TransReactivePower.m

13 lines
535 B
Mathematica
Raw Normal View History

2013-04-22 21:56:16 +08:00
function [ output_args ] = TransReactivePower( newwordParameter,cmpY,Volt,VAngel )
%TRANSREAVTIVEPOWER Summary of this function goes here
% Detailed explanation goes here
transI=newwordParameter.trans.transI;
transJ=newwordParameter.trans.transJ;
cmpSEV=Volt.*exp(1j*VAngel); %<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ѹ
% cmpY=sparse(r,c,Y.*(1j*YAngel),length(Volt),length(Volt));
cmpSETransI=(cmpSEV(transI)-cmpSEV(transJ)).*cmpY(sub2ind(size(cmpY),transI,transJ));
SETransQ=imag((cmpSEV(transI)-cmpSEV(transJ)).*conj(cmpSETransI));
output_args=full(SETransQ);
2013-04-22 21:56:16 +08:00
end