diff --git a/OPF.m b/OPF.m index 88de316..b1f4374 100644 --- a/OPF.m +++ b/OPF.m @@ -224,24 +224,19 @@ wPD=1./(abs(mPD3P*sigma).^2); wQD=1./(abs(mQD3P*sigma).^2); %% 没有量测量 -noLoadi=[1,5,6,10,11]; -% noLoadi=Loadi; +% noLoadi=[1,5,6,10,11]; +noLoadi=Loadi; % noLoadi=[1,11]; % noLoadi=[]; noPQi3P=zeros(length(noLoadi)*3,1); noPQi3P(1:3:end)=(noLoadi-1)*3+1; noPQi3P(2:3:end)=(noLoadi-1)*3+2; noPQi3P(3:3:end)=(noLoadi-1)*3+3; -noLoadi=noPQi3P; +% noLoadi=noPQi3P; wVolt(setdiff(1:length(wVolt),Loadi))=0;%只有负荷处才有电压量测。 wVolt( noLoadi)=0; wPD(ismember( Loadi,noLoadi))=0; wQD(ismember( Loadi,noLoadi))=0; - - -% wVolt(5)=0; -% wPD(5)=0; -% wQD(5)=0; %% % RestraintCount=size(Loadi,1)*2+length(rVoltABCV); %约束条件数 RestraintCount=size(Loadi,1)*2; %约束条件数 @@ -253,7 +248,7 @@ Init_Y=sparse(1,2*Busnum*3);% Gap=(Init_L*Init_Z'-Init_U*Init_W'); KK=0; ContrlCount=size(Loadi,1)*2+Busnum*6; -kmax=500; +kmax=200; %% Precision=1e-5; CenterA=0.1;