diff --git a/FormLw.m b/FormLw.m index d16c155..9683a8d 100644 --- a/FormLw.m +++ b/FormLw.m @@ -26,8 +26,8 @@ function Lw=FormLw(Mat_G,Init_U,Loadi,rPD3P,rQD3P,mVoltABCV) % QDU=1.0*PD(Loadi).*sqrt(1 -PF.^2)./PF; % PDU=1.8*ones(length(Loadi),1); % QDU=1.8*ones(length(Loadi),1); -PDU=1.25*rPD3P.*(1+unifrnd(-0.2,0.2,length(rPD3P),1)); -QDU=1.25*rQD3P.*(1+unifrnd(-0.2,0.2,length(rPD3P),1)); +PDU=1.3*rPD3P.*(1+unifrnd(-0.15,0.15,length(rPD3P),1)); +QDU=1.3*rQD3P.*(1+unifrnd(-0.15,0.15,length(rPD3P),1)); % t1=([PDU',QDU',1*ones(length(Volt),1)'])'; t1=([PDU',QDU'])'; t2=Mat_G+Init_U'-t1; diff --git a/FormLz.m b/FormLz.m index b429825..ea56136 100644 --- a/FormLz.m +++ b/FormLz.m @@ -23,8 +23,8 @@ function Lz=FormLz(Mat_G,Init_L,Loadi,rPD3P,rQD3P,mVoltABCV) % QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF; % PDL=0*ones(length(Loadi),1); % QDL=0*ones(length(Loadi),1); -PDL=0.75*rPD3P.*(1+unifrnd(-0.15,0.15,length(rPD3P),1)); -QDL=0.75*rQD3P.*(1+unifrnd(-0.15,0.15,length(rQD3P),1)); +PDL=0.70*rPD3P.*(1+unifrnd(-0.15,0.15,length(rPD3P),1)); +QDL=0.70*rQD3P.*(1+unifrnd(-0.15,0.15,length(rQD3P),1)); % t1=([PDL',QDL',0.88*ones(length(Volt),1)'])'; t1=([PDL',QDL'])'; t2=Mat_G-Init_L'-t1; diff --git a/OPF.m b/OPF.m index cbc0be3..516ea27 100644 --- a/OPF.m +++ b/OPF.m @@ -225,10 +225,10 @@ wQD=1./(abs(mQD3P*sigma).^2); %% 没有量测量 % noLoadi=[1,5,6,10,11]; -% noLoadi=[1,8,9,10,12]; +noLoadi=[1,8,10,12]; % noLoadi=Loadi; % noLoadi=[1,11]; -noLoadi=[]; +% noLoadi=[]; noPQi3P=zeros(length(noLoadi)*3,1); noPQi3P(1:3:end)=(noLoadi-1)*3+1; noPQi3P(2:3:end)=(noLoadi-1)*3+2; @@ -255,7 +255,7 @@ Precision=1e-5; CenterA=0.1; %% 加误差 Volt=1*Vp3m; -Volt=rVoltABCV; +% Volt=mVoltABCV; UAngel=Vp3a; while(abs(Gap)>Precision) if KK>kmax diff --git a/feeder13/data1.txt b/feeder13/data1.txt index da177f9..27e0b6e 100644 --- a/feeder13/data1.txt +++ b/feeder13/data1.txt @@ -1,21 +1,18 @@ 650 4.16 0 -600 632 645 500 -602 632 633 500 -601 633 634 600 -600 645 646 300 -601 650 632 2000 -602 684 652 800 -601 632 671 2000 -600 671 684 300 -600 671 692 300 -602 684 611 300 -601 692 675 500 +600 632 645 50 +602 632 633 50 +601 633 634 60 +600 645 646 30 +601 650 632 200 +602 684 652 80 +601 632 671 200 +600 671 684 30 +600 671 692 30 +602 684 611 30 +601 692 675 50 0 -611 100 100 100 -652 200 200 200 -684 200 200 200 -692 200 200 200 +692 0 0 0 0 634 120 90 120 90 160 110 645 40 30 50 40 60 50