distributionnetwork-power2current/FormLwstate1.m

function Lw=FormLwstate1(Loadi,Mat_G,Init_U,I1measurement,dI_F,flag)

Ir=real(I1measurement);
pIr=find(Ir>0);
nIr=find(Ir<0);
Ii=imag(I1measurement);
pIi=find(Ii>0);
nIi=find(Ii<0);

% Loadi=[1,4,5,6,8,10];

% lower=-0.2*sparse(ones(length(Loadi)*2,1));
% upper=ones(length(Loadi)*2,1);
% upper(pIr)=1.4*Ir(pIr);
% upper(nIr)=0.6*Ir(nIr);
% upper(pIi+length(Ir))=1.4*Ii(pIi);
% upper(nIi+length(Ir))=0.6*Ii(nIi);
upper=ones(length(Loadi)*2,1);
upper(pIr)=abs(real(dI_F(pIr,flag)))+Ir(pIr);
upper(nIr)=abs(real(dI_F(nIr,flag)))+Ir(nIr);
upper(pIi+length(Ir))=abs(imag(dI_F(pIi,flag)))+Ii(pIi);
upper(nIi+length(Ir))=abs(imag(dI_F(nIi,flag)))+Ii(nIi);


%̫С<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>Ҫ<EFBFBD>ſ<EFBFBD>һЩ
% tooSmall=find(abs(Ir)<0.0005);
% upper(tooSmall)=0.2;
% tooSmall=find(abs(Ii)<0.0005);
% upper(tooSmall+length(Ir))=0.2;

Lw=Mat_G+Init_U-upper;
end