distributionnetwork-power2c.../FormLwstate1.m

21 lines
481 B
Matlab

function Lw=FormLwstate1(Loadi,Mat_G,Init_U,I1measurement)
Ir=real(I1measurement);
pIr=find(Ir>0);
nIr=find(Ir<0);
Ii=imag(I1measurement);
pIi=find(Ii>0);
nIi=find(Ii<0);
% lower=-0.2*sparse(ones(length(Loadi)*2,1));
upper=ones(length(Loadi)*2,1);
upper(pIr)=1.2*Ir(pIr);
upper(nIr)=0.8*Ir(nIr);
upper(pIi+length(Ir))=1.2*Ii(pIi);
upper(nIi+length(Ir))=0.8*Ii(nIi);
% upper([4,5,6])=[1;1;1];
% upper([4,5,6])=0.2*Ii(nIi);
% upper=ones(length(Ir)*2,1);
Lw=Mat_G+Init_U-upper;
end