distributionnetwork-power2c.../FormLzstate1.m

18 lines
453 B
Matlab

function Lz=FormLzstate1(Loadi,Mat_G,Init_L,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));
lower=ones(length(Loadi)*2,1);
lower(pIr)=0.9*Ir(pIr);
lower(nIr)=1.1*Ir(nIr);
lower(pIi+length(Ir))=0.9*Ii(pIi);
lower(nIi+length(Ir))=1.1*Ii(nIi);
% lower=-ones(length(Ir)*2,1);
% lower([4,5,6])=[-1;-1;-1];
Lz=Mat_G-Init_L-lower;
end