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

如果不置0置1的话，潮流方程又不平衡了，还是该回来。 

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk 2014-10-19 15:58:55 +08:00
parent 4544a6d009
commit a3b1c5a30e
6 changed files with 12 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
FormLw.m
View File

@ -1,5 +1,5 @@
function Lw=FormLw(Loadi,Mat_G,Init_U)
upper=0.1*sparse(ones(length(Loadi)*1,1));
upper=0.01*sparse(ones(length(Loadi)*1,1));
Lw=Mat_G+Init_U-upper;
end

3
IPMLoop.m
View File

@ -29,7 +29,6 @@ while(abs(Gap)>0.000001)
deltH=func_deltH(busNum,fsY1,Loadi,Balance);
%% 形成不等式约束的雅克比
deltG=func_deltG(busNum,Loadi,I1r,I1i);
deltG=zeros(size(deltG));
%%
L_1Z=diag(Init_Z./Init_L);
U_1W=diag(Init_W./Init_U);
@ -41,7 +40,6 @@ while(abs(Gap)>0.000001)
ddh=0;
%% 开始构建ddg
ddg=func_ddg(busNum,Loadi,Init_Z,Init_W);
ddg=zeros(size(ddg));
%% 开始构建deltF
deltF=func_deltF(Vmeasurement,Imeasurement,busNum,fsY1,Loadi,V1r,V1i,I1r,I1i);
% deltF=0;
@ -49,7 +47,6 @@ while(abs(Gap)>0.000001)
Luu=Init_U.*Init_W+Init_u*ones(RestraintCount,1);
Lul=Init_L.*Init_Z-Init_u*ones(RestraintCount,1);
Mat_G=FormG(I1r,I1i);
Mat_G=zeros(size(Mat_G));
Mat_H=FormH(fsY1,Loadi,V1r,V1i,I1r,I1i,BalI1r,BalI1i,Balance);
Ly=Mat_H;
Lz=FormLz(Loadi,Mat_G,Init_L);

4
Modification.m
View File

@ -12,9 +12,9 @@ Init_Y=Init_Y+AlphaD*deltY;
%PG(PGi)=PG(PGi)+deltX(1:size(PGi,1));
% PG(PGi)=PG(PGi)+AlphaP*deltX(1:size(PGi,1));
V1r=V1r+AlphaP*deltX(1:Busnum);
% V1r(Balance)=Vref;
V1r(Balance)=Vref;
V1i=V1i+AlphaP*deltX(Busnum+1:2*Busnum);
% V1i(Balance)=0;
V1i(Balance)=0;
%QG(PVi)=QG(PVi)+deltX(size(PGi,1)+1:size(PVi,1)+size(PGi,1) );
% QG(PVi)=QG(PVi)+AlphaP*deltX(size(PGi,1)+1:size(PVi,1)+size(PGi,1) );
t=deltX(Busnum*2+1:ContrlCount);

14
SolveIt.m
View File

@ -11,15 +11,15 @@ yy=[LxComa;-Ly];
% t=size(PVi,1)+size(PGi,1);%电压不变
% yy(t+Balance)=0;
% aa(Balance,:)=0;
% aa(:,Balance)=0;
% aa=aa+sparse(Balance,Balance,ones(length(Balance),1),ContrlCount+2*Busnum,ContrlCount+2*Busnum);
aa(Balance,:)=0;
aa(:,Balance)=0;
aa=aa+sparse(Balance,Balance,ones(length(Balance),1),ContrlCount+2*Busnum,ContrlCount+2*Busnum);
%暂时改一下
% t=Busnum;%相角不变
% aa(t+Balance,:)=0;
% aa(:,t+Balance)=0;
% aa(t+Balance,t+Balance)=1;
t=Busnum;%相角不变
aa(t+Balance,:)=0;
aa(:,t+Balance)=0;
aa(t+Balance,t+Balance)=1;
dxdy=aa\yy;
dX=dxdy(1:ContrlCount);
dY=dxdy(ContrlCount+1:ContrlCount+2*Busnum);

2
func_deltF.m
View File

@ -11,7 +11,7 @@ deltF=[
%zeros(busNum*2,1);
sparse(Loadi,1,-2*(real(Vmeasurement)-V1r(Loadi)),busNum,1);
sparse(Loadi,1,-2*(imag(Vmeasurement)-V1i(Loadi)),busNum,1);
-0*( [real(Imeasurement);imag(Imeasurement)]-[I1r;I1i]);
-2*( [real(Imeasurement);imag(Imeasurement)]-[I1r;I1i]);
];
end

2
func_deltdeltF.m
View File

@ -4,7 +4,7 @@ deltdeltF=[
%zeros(busNum*2,busNum*2+length(Loadi)*2);
sparse(Loadi,Loadi,2,busNum,busNum*2+length(Loadi)*2);
sparse(Loadi,busNum+Loadi,2,busNum,busNum*2+length(Loadi)*2);
zeros(length(Loadi)*2,busNum*2),0*eye(length(Loadi)*2);
zeros(length(Loadi)*2,busNum*2),2*eye(length(Loadi)*2);
%sparse(Loadi,busNum*2+Loadi,2,busNum,busNum*2+length(Loadi)*2);
];
end