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

编制了新的节点文件。 

Signed-off-by: dmy@lab <dmy@lab.lab>
This commit is contained in:
dmy@lab 2015-04-02 15:55:54 +08:00
parent 4543616833
commit 0ad8514b6a
4 changed files with 79 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
IPMLoop.m
View File

@ -8,7 +8,6 @@ KK=0;
plotGap=zeros(1,60); plotGap=zeros(1,60);
%初始化 %初始化
%状态量为 SEPD SEQD SEVmf1 SEVaf1 %状态量为 SEPD SEQD SEVmf1 SEVaf1
state=1;%state1表示用l<=I1r<=u 这样的约束形式
RestraintCount=length(Loadi)*2; RestraintCount=length(Loadi)*2;
ContrlCount=busNum*2+length(Loadi)*2; ContrlCount=busNum*2+length(Loadi)*2;
CenterA=0.1; CenterA=0.1;
@ -18,9 +17,9 @@ Init_L=1*sparse(ones(RestraintCount,1));
Init_U=1*sparse(ones(RestraintCount,1)); Init_U=1*sparse(ones(RestraintCount,1));
Init_Y=sparse(2*busNum,1);%等式约束乘子 Init_Y=sparse(2*busNum,1);%等式约束乘子
Gap=(Init_L'*Init_Z-Init_U'*Init_W); Gap=(Init_L'*Init_Z-Init_U'*Init_W);
kmax=60; kmax=15;
while(abs(Gap)>0.00001) while(abs(Gap)>0.00001)
if KK>kmax if KK>=kmax
break; break;
end end
plotGap(KK+1)=Gap; plotGap(KK+1)=Gap;
@ -64,7 +63,7 @@ while(abs(Gap)>0.00001)
KK=KK+1; KK=KK+1;
end end
% f=sum(([real(I1measurement);imag(I1measurement)]-[I1r;I1i]).^2)+sum((real(V1measurement)-V1r(Loadi)).^2)+sum((imag(V1measurement)-V1i(Loadi)).^2); % f=sum(([real(I1measurement);imag(I1measurement)]-[I1r;I1i]).^2)+sum((real(V1measurement)-V1r(Loadi)).^2)+sum((imag(V1measurement)-V1i(Loadi)).^2);
if abs(Gap)<0.00001 if abs(Gap)<0.00001 && KK<kmax
isConverged=1; isConverged=1;
else else
isConverged=0; isConverged=0;

18
feeder13/data - 副本.txt
View File

@ -1,16 +1,16 @@
650 650
0 0
603 632 645 500 600 632 645 500
602 632 633 500 602 632 633 500
-1 633 634 0 601 633 634 0
603 645 646 300 600 645 646 300
601 650 632 2000 601 650 632 2000
607 684 652 800 602 684 652 800
601 632 671 2000 601 632 671 2000
604 671 684 300 600 671 684 300
601 671 680 1000 602 671 680 1000
-2 671 692 0 600 671 692 0
605 684 611 300 602 684 611 300
606 692 675 500 601 692 675 500
0 0

46
feeder13/data1.txt
View File

@ -1,14 +1,50 @@
650 4.16 650 4.16
0 0
600 632 645 500
602 632 633 500
601 633 634 600
600 645 646 300
601 650 632 2000 601 650 632 2000
602 632 631 1800 602 684 652 800
602 631 630 1700 601 632 671 2000
600 671 684 300
602 671 680 1000
600 671 692 300
602 684 611 300
601 692 675 500
0 0
632 0 0 0 632 0 0 0
0 0
630 80 70 60 50 40 30 634 160 110 120 90 120 90
631 40 10 11 50 30 10 645 123 102 170 125 133 110
632 120 90 120 80 90 70 646 172 150 230 132 211 160
652 128 86 30 10 80 70
675 85 10 68 60 20 12
611 52 42 85 60 70 40
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
602 671 680 1000
600 671 692 300
602 684 611 300
601 692 675 500
0
632 0 0 0
0
634 160 110 120 90 120 90
645 123 102 170 125 133 110
646 172 150 230 132 0 211 160
652 128 86 30 10 80 70
671 385 220 385 220 385 220
675 485 190 68 60 290 212
692 142 140 153 121 170 151
611 156 142 185 60 170 80
0 0

26
run.m
View File

@ -60,6 +60,13 @@ VoltpA=sparse(ones(busNum,1));
VoltpB=sparse(ones(busNum,1)).*exp(1j*-120/180*pi); VoltpB=sparse(ones(busNum,1)).*exp(1j*-120/180*pi);
VoltpC=sparse(ones(busNum,1)).*exp(1j*+120/180*pi); VoltpC=sparse(ones(busNum,1)).*exp(1j*+120/180*pi);
xs=1;
phaseASpotLoadP=phaseASpotLoadP*xs;
phaseBSpotLoadP=phaseBSpotLoadP*xs;
phaseCSpotLoadP=phaseCSpotLoadP*xs;
phaseASpotLoadQ=phaseASpotLoadQ*xs;
phaseBSpotLoadQ=phaseBSpotLoadQ*xs;
phaseCSpotLoadQ=phaseCSpotLoadQ*xs;
while(k<=kmax && maxD> EPS) while(k<=kmax && maxD> EPS)
k=k+1; k=k+1;
@ -324,6 +331,25 @@ deltmCurpC=conj((mphaseCSpotLoadP+1j*mphaseCSpotLoadQ)./mVoltpC)*0.2;
deltF012=Tp2f*conj([deltmCurpA';deltmCurpB';deltmCurpC']); deltF012=Tp2f*conj([deltmCurpA';deltmCurpB';deltmCurpC']);
dI_F=conj(deltF012)'; dI_F=conj(deltF012)';
dI_F=dI_F(Loadi,:); dI_F=dI_F(Loadi,:);
% deltmCurpA1=conj((mphaseASpotLoadP+1j*mphaseASpotLoadQ)./mVoltpA)*0.8;
% deltmCurpB1=conj((mphaseBSpotLoadP+1j*mphaseBSpotLoadQ)./mVoltpB)*0.8;
% deltmCurpC1=conj((mphaseCSpotLoadP+1j*mphaseCSpotLoadQ)./mVoltpC)*0.8;
%
% deltF012=Tp2f*conj([-deltmCurpA1'+mCurpA';-deltmCurpB1'+mCurpB';-deltmCurpC1'+mCurpC']);
% dI_F=conj(deltF012)';
% dI_F=dI_F(Loadi,:);
%
%
% deltmCurpA1=conj((mphaseASpotLoadP+1j*mphaseASpotLoadQ)./mVoltpA)*1.2;
% deltmCurpB1=conj((mphaseBSpotLoadP+1j*mphaseBSpotLoadQ)./mVoltpB)*1.2;
% deltmCurpC1=conj((mphaseCSpotLoadP+1j*mphaseCSpotLoadQ)./mVoltpC)*1.2;
%
% deltF012=Tp2f*conj([deltmCurpA1'-mCurpA';deltmCurpB1'-mCurpB';deltmCurpC1'-mCurpC']);
% dI_F=conj(deltF012)';
% dI_F=dI_F(Loadi,:);
tic tic
for II=1:3 for II=1:3