facat
/
ipmthreephasese
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

1.修改了迭代最大次数 

2.先用这个算例数据和参数。

Signed-off-by: dmy@lab <dmy@lab.lab>
This commit is contained in:
dmy@lab 2015-04-05 19:31:39 +08:00
parent 109c9b3800
commit e26a2f9ef0
1 changed files with 4 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
OPF.m
View File

@ -224,24 +224,19 @@ wPD=1./(abs(mPD3P*sigma).^2);
wQD=1./(abs(mQD3P*sigma).^2); wQD=1./(abs(mQD3P*sigma).^2);
%% 没有量测量 %% 没有量测量
noLoadi=[1,5,6,10,11]; % noLoadi=[1,5,6,10,11];
% noLoadi=Loadi; noLoadi=Loadi;
% noLoadi=[1,11]; % noLoadi=[1,11];
% noLoadi=[]; % noLoadi=[];
noPQi3P=zeros(length(noLoadi)*3,1); noPQi3P=zeros(length(noLoadi)*3,1);
noPQi3P(1:3:end)=(noLoadi-1)*3+1; noPQi3P(1:3:end)=(noLoadi-1)*3+1;
noPQi3P(2:3:end)=(noLoadi-1)*3+2; noPQi3P(2:3:end)=(noLoadi-1)*3+2;
noPQi3P(3:3:end)=(noLoadi-1)*3+3; noPQi3P(3:3:end)=(noLoadi-1)*3+3;
noLoadi=noPQi3P; % noLoadi=noPQi3P;
wVolt(setdiff(1:length(wVolt),Loadi))=0;%只有负荷处才有电压量测。 wVolt(setdiff(1:length(wVolt),Loadi))=0;%只有负荷处才有电压量测。
wVolt( noLoadi)=0; wVolt( noLoadi)=0;
wPD(ismember( Loadi,noLoadi))=0; wPD(ismember( Loadi,noLoadi))=0;
wQD(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+length(rVoltABCV); %约束条件数
RestraintCount=size(Loadi,1)*2; %约束条件数 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'); Gap=(Init_L*Init_Z'-Init_U*Init_W');
KK=0; KK=0;
ContrlCount=size(Loadi,1)*2+Busnum*6; ContrlCount=size(Loadi,1)*2+Busnum*6;
kmax=500; kmax=200;
%% %%
Precision=1e-5; Precision=1e-5;
CenterA=0.1; CenterA=0.1;