__author__ = 'dmy' import numpy def getPowerflowEquation(y,Vm,Va,Pd,Qd,Pg,Qg): cmpV=Vm*numpy.exp(1j*Va)#复数电压 S=cmpV*(numpy.dot(y.toarray(),cmpV).conj()) P=numpy.real(S)+Pd-Pg Q=numpy.imag(S)+Qd-Qg return (P,Q)