facat
/
distributionnetwork-power2current
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

有点问题 

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk 2014-11-11 11:37:18 +08:00
parent 3db4877bb7
commit 00a4cfad44
4 changed files with 11 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
FormLwstate1.m
View File

@ -8,10 +8,10 @@ pIi=find(Ii>0);
nIi=find(Ii<0); nIi=find(Ii<0);
% lower=-0.2*sparse(ones(length(Loadi)*2,1)); % lower=-0.2*sparse(ones(length(Loadi)*2,1));
upper=ones(length(Loadi)*2,1); upper=ones(length(Loadi)*2,1);
upper(pIr)=1.5*Ir(pIr); upper(pIr)=1.2*Ir(pIr);
upper(nIr)=0.5*Ir(nIr); upper(nIr)=0.8*Ir(nIr);
upper(pIi+length(Ir))=1.5*Ii(pIi); upper(pIi+length(Ir))=1.2*Ii(pIi);
upper(nIi+length(Ir))=0.5*Ii(nIi); upper(nIi+length(Ir))=0.8*Ii(nIi);
% upper([4,5,6])=[1;1;1]; % upper([4,5,6])=[1;1;1];
% upper([4,5,6])=0.2*Ii(nIi); % upper([4,5,6])=0.2*Ii(nIi);

8
FormLzstate1.m
View File

@ -7,10 +7,10 @@ pIi=find(Ii>0);
nIi=find(Ii<0); nIi=find(Ii<0);
% lower=-0.2*sparse(ones(length(Loadi)*2,1)); % lower=-0.2*sparse(ones(length(Loadi)*2,1));
lower=ones(length(Loadi)*2,1); lower=ones(length(Loadi)*2,1);
lower(pIr)=0.5*Ir(pIr); lower(pIr)=0.8*Ir(pIr);
lower(nIr)=1.5*Ir(nIr); lower(nIr)=1.2*Ir(nIr);
lower(pIi+length(Ir))=0.5*Ii(pIi); lower(pIi+length(Ir))=0.8*Ii(pIi);
lower(nIi+length(Ir))=1.5*Ii(nIi); lower(nIi+length(Ir))=1.2*Ii(nIi);
% lower=-ones(length(Ir)*2,1); % lower=-ones(length(Ir)*2,1);
% lower([4,5,6])=[-1;-1;-1]; % lower([4,5,6])=[-1;-1;-1];

4
IPMLoop.m
View File

@ -43,7 +43,7 @@ while(abs(Gap)>0.00001)
% U_1W=diag(Init_W./Init_U); % U_1W=diag(Init_W./Init_U);
%% 形成海森阵 %% 形成海森阵
deltdeltF=func_deltdeltF(busNum,fsY1,Loadi,wV1r,wV1i,wI1r,wI1i,V1measurement,V1r,V1i); deltdeltF=func_deltdeltF(busNum,fsY1,Loadi,wV1r,wV1i,wI1r,wI1i,V1measurement,V1r,V1i);
% deltdeltF=0; deltdeltF=0;
%% 形成ddHy %% 形成ddHy
% ddh=func_ddh(busNum,Loadi,Init_Z,Init_W); % ddh=func_ddh(busNum,Loadi,Init_Z,Init_W);
ddh=0; ddh=0;
@ -55,7 +55,7 @@ while(abs(Gap)>0.00001)
end end
%% 开始构建deltF %% 开始构建deltF
deltF=func_deltF(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,busNum,fsY1,Loadi,V1r,V1i,I1r,I1i); deltF=func_deltF(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,busNum,fsY1,Loadi,V1r,V1i,I1r,I1i);
% deltF=0; deltF=0;
%% %%
Luu=Init_U.*Init_W+Init_u*ones(RestraintCount,1); Luu=Init_U.*Init_W+Init_u*ones(RestraintCount,1);
Lul=Init_L.*Init_Z-Init_u*ones(RestraintCount,1); Lul=Init_L.*Init_Z-Init_u*ones(RestraintCount,1);

1
run.m
View File

@ -261,6 +261,7 @@ I0measurement=mIf0(Loadi);
wV0r=abs(real(V0measurement)).*sigma; wV0r=abs(real(V0measurement)).*sigma;
wV0i=abs(imag(V0measurement)).*sigma; wV0i=abs(imag(V0measurement)).*sigma;
wI0r=abs(real(I0measurement)).*sigma; wI0r=abs(real(I0measurement)).*sigma;
wI0r(abs(wI0r)<1e-5)=1e10;
wI0i=abs(imag(I0measurement)).*sigma; wI0i=abs(imag(I0measurement)).*sigma;
% matlabpool local 3 % matlabpool local 3
tic tic