stateestimation-ipm-forcomparison/@SEOpti/equ.m

function [ out_arg,this ] = equ(this, x )
%EQU Summary of this function goes here
%   Detailed explanation goes here
SEVolt=x(1:this.Busnum);
SEVAngel=x(this.Busnum+1:2*this.Busnum);
cmpSEV=SEVolt.*exp(1j*SEVAngel); %<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ѹ
PQ=diag(cmpSEV)*conj(this.cmpY*cmpSEV);
% YAngle=sparse(this.r,this.c,SEVAngel(this.r)-SEVAngel(this.c)-this.Yangle,this.Busnum,this.Busnum);
% PQ1=diag(SEVolt)*(this.Y.*cos(YAngle))*SEVolt;
% PQ2=diag(SEVolt)*(this.Y.*sin(YAngle))*SEVolt;
out_arg=[real(PQ(this.zerosInjectionIndex));imag(PQ(this.zerosInjectionIndex));];
% out_arg=[PQ1(this.zerosInjectionIndex);PQ2(this.zerosInjectionIndex);];
% out_arg=[out_arg;x(this.Balance+this.Busnum)];
out_arg=[out_arg;x(this.Balance+this.Busnum)];%<EFBFBD><EFBFBD><EFBFBD>ǵ<EFBFBD><EFBFBD><EFBFBD>0
BalanceVolt=this.BalanceVolt;
out_arg=[out_arg;SEVolt(this.Balance)-BalanceVolt];%ƽ<EFBFBD><EFBFBD><EFBFBD>ڵ<EFBFBD><EFBFBD><EFBFBD>ѹ
out_arg=full(out_arg);
this.mnle=zeros(length(out_arg),1);
this.mnlrhs=this.mnle;
end