distributionnetwork-power2c.../FormLwstate1.m

46 lines
1006 B
Matlab

function Lw=FormLwstate1(Loadi,Mat_G,Init_U,I1measurement,dI_F,flag,guessIf1,rIf1,upper)
% 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(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);
% if flag~=2
% upper(pIr)=0.8;
% upper(nIr)=0.8;
% upper(pIi+length(Ir))=0.8;
% upper(nIi+length(Ir))=0.8;
% end
% if real(rIf1(Loadi))>upper( union( pIr,nIr) )
% warning('not bounded');
% end
%
% if imag(rIf1(Loadi))>upper(union( pIr,nIr)+length(Ir))
% warning('not bounded');
% end
%̫СµÄÊý¶¼Òª·Å¿íһЩ
% tooSmall=find(abs(Ir)<0.05);
% upper(tooSmall)=0.8;
% tooSmall=find(abs(Ii)<0.05);
% upper(tooSmall+length(Ir))=0.8;
Lw=Mat_G+Init_U-upper;
end