facat/threephasese-onlyenquationconstraint
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

提交第一版

Browse Source 
Signed-off-by: dmy@lab <dmy@lab.lab>
This commit is contained in:
dmy@lab
2015-05-08 22:02:19 +08:00
commit 627ba1eb72
210 changed files with 4639 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
.gitignore vendored Normal file
View File

@@ -0,0 +1 @@
*.asv

18
AssignXX.m Normal file
View File

@@ -0,0 +1,18 @@
function [deltX,deltY]=AssignXX(XX,ContrlCount,Loadi,Balance,Busnum)
% deltX=XX(1:14);
% deltY=XX(15:24);
% deltZ=XX(25:38);
% deltW=XX(39:52);
% deltL=XX(53:66);
% deltU=XX(67:80);
deltX=XX(1:ContrlCount);
t=size(Loadi,1)*2;
deltX(t+(Balance-1)*3+1)=0;
deltX(t+(Balance-1)*3+2)=0;
deltX(t+(Balance-1)*3+3)=0;
t=size(Loadi,1)*2+Busnum*3;
deltX(t+(Balance-1)*3+1)=0;
deltX(t+(Balance-1)*3+2)=0;
deltX(t+(Balance-1)*3+3)=0;
deltY=XX(ContrlCount+1:end);
end

25
CalPGQG.m Normal file
View File

@@ -0,0 +1,25 @@
function [ output_args ] = CalPGQG(Balance,phaseABCY,VoltpABC,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP,phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ )
%计算PG和QG
busNum=size(phaseABCY,1)/3;
VoltpA=VoltpABC(1,:);
VoltpB=VoltpABC(2,:);
VoltpC=VoltpABC(3,:);
Vp3=sparse(busNum*3,1);
Vp3(1:3:end)=VoltpA;
Vp3(2:3:end)=VoltpB;
Vp3(3:3:end)=VoltpC;
Ip3=phaseABCY*Vp3;
Sp3=Vp3.*conj(Ip3);
%VoltpABC.*conj(IpABC);
pLoadABC=sparse(length(phaseASpotLoadP)*3,1);
pLoadABC(1:3:end)=phaseASpotLoadP+1j*phaseASpotLoadQ;
pLoadABC(2:3:end)=phaseBSpotLoadP+1j*phaseBSpotLoadQ;
pLoadABC(3:3:end)=phaseCSpotLoadP+1j*phaseCSpotLoadQ;
ck=Sp3+pLoadABC;
% ck(3*(Balance-1)+1:3*(Balance-1)+3)=0;
% if any(abs(ck)>1e-5)
% fprintf('反推回的功率不匹配。潮流方程约束不满足。\n');
% end
output_args=ck;
end

9
FormAlphaD.m Normal file
View File

@@ -0,0 +1,9 @@
function AlphaD=FormAlphaD(Init_Z,deltZ,Init_W,deltW)
tdeltZinx=find(deltZ<0);
tdeltWinx=find(deltW>0);
t1=-Init_Z(tdeltZinx)./deltZ(tdeltZinx)';
t2=-Init_W(tdeltWinx)./deltW(tdeltWinx)';
t3=[t1,t2];
t4=min(t3);
AlphaD=0.9995*min([t4 1]);
end

9
FormAlphaP.m Normal file
View File

@@ -0,0 +1,9 @@
function AlphaP=FormAlphaP(Init_L,deltL,Init_U,deltU)
tdeltLinx=find(deltL<0);
tdeltUinx=find(deltU<0);
t1=-Init_L(tdeltLinx)./deltL(tdeltLinx)';
t2=-Init_U(tdeltUinx)./deltU(tdeltUinx)';
t3=[t1,t2];
t4=min(t3);
AlphaP=0.9995*min([t4 1]);
end

8
FormG.m Normal file
View File

@@ -0,0 +1,8 @@
function Mat_G=FormG(Volt,PD,QD,Loadi)
Mat_G=[
sparse(PD);
sparse(QD);
% Volt;
];
end

14
FormH.m Normal file
View File

@@ -0,0 +1,14 @@
function Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi)
%%
%QDcos=textread('D:\Project\×îС»¯³±Á÷\×îС³±Á÷ËãÀý\Ïɺ£919PDQDglys.txt');
%QD(QD~=0)=PD(QD~=0)./tan(QDcos);
%QD(QD_NON_ZERO_IND)=QD_NON_ZERO;
%%
%PD(Loadi)=QD(Loadi)./tan(acos(0.98));
AngleIJ=sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum*3,Busnum*3);
dP=PG-sparse(Loadi,1,PD,Busnum*3,1)-diag(Volt)*Y.*cos(AngleIJ)*Volt;
dQ=QG-sparse(Loadi,1,QD,Busnum*3,1)-diag(Volt)*Y.*sin(AngleIJ)*Volt;
Mat_H=[dP;dQ;];
end

38
FormLw.m Normal file
View File

@@ -0,0 +1,38 @@
function Lw=FormLw(Mat_G,Init_U,Loadi,rPD3P,rQD3P,mVoltABCV,PDU,QDU)
% KK=999;
% %PU=GenU(:,2);%·¢µç»úÓй¦ÉϽç
%
% VoltU=(1.1)*ones(1,Busnum);
% %VoltU=10*ones(1,Busnum);
% PDU=PD0(Loadi);
% % PDU=noDataTransCapacity;
% PDU(PDU>0)=1.200*PDU(PDU>0);
% PDU(PDU<0)=0.800*PDU(PDU<0);
% PDU(PDU==0)=0.400;
% realPD=PD0(Loadi);
% indPD=find(realPD>0);
% PDU(indPD(3:12:end))=1.55*realPD(indPD(3:12:end));
% PDU(indPD(9:12:end))=1.05*realPD(indPD(9:12:end));
% %PDU=10*ones(length(Loadi),1);
% QDU=QD0(Loadi);
% QDU(QDU>0)=1.200*QDU(QDU>0);
% QDU(QDU<0)=0.800*QDU(QDU<0);
% QDU(QDU==0)=0.400;
% realQD=QD0(Loadi);
% indQD=find(realQD>0);
% QDU(indQD(3:12:end))=1.55*realQD(indQD(3:12:end));
% QDU(indQD(9:12:end))=1.05*realQD(indQD(9:12:end));
% PF=0.85;
% 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.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;
Lw=t2;
end

4
FormLx.m Normal file
View File

@@ -0,0 +1,4 @@
function Lx=FormLx(deltF,deltH,Init_Y)
t1=deltF-deltH*Init_Y';
Lx=t1;
end

4
FormLxComa.m Normal file
View File

@@ -0,0 +1,4 @@
function LxComa=FormLxComa(Lx)
%t1=deltG*(Init_Z'+Init_W');%%
LxComa=Lx;
end

35
FormLz.m Normal file
View File

@@ -0,0 +1,35 @@
function Lz=FormLz(Mat_G,Init_L,Loadi,rPD3P,rQD3P,mVoltABCV,PDL,QDL)
% KK=999;
%
% VoltL=(0.9)*ones(1,Busnum);
% %VoltL=-10*ones(1,Busnum);
% PDL=PD0(Loadi);
% PDL(PDL>0)=0.800*PDL(PDL>0);
% PDL(PDL<0)=1.200*PDL(PDL<0);
% PDL(PDL==0)=-0.400;
% realPD=PD0(Loadi);
% indPD=find(realPD>0);
% PDL(indPD(3:12:end))=0.95*realPD(indPD(3:12:end));
% PDL(indPD(9:12:end))=0.45*realPD(indPD(9:12:end));
% %PDL=-10*ones(length(Loadi),1);
% QDL=QD0(Loadi);
% QDL(QDL>0)=0.800*QDL(QDL>0);
% QDL(QDL<0)=1.200*QDL(QDL<0);
% QDL(QDL==0)=-0.400;
% realQD=QD0(Loadi);
% indQD=find(realQD>0);
% QDL(indQD(3:12:end))=0.95*realQD(indQD(3:12:end));
% QDL(indQD(9:12:end))=0.95*realQD(indQD(9:12:end));
% QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF;
% PDL=0*ones(length(Loadi),1);
% QDL=0*ones(length(Loadi),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;
Lz=t2;
end

10
FormYY.m Normal file
View File

@@ -0,0 +1,10 @@
function YY=FormYY(Lul,Lz,Ly,Luu,Lw,Lx)
YY=[
Lx;
-Ly;
-Lz;
-Lw;
-Lul;
-Luu;
];
end

96
IEEE30.dat Normal file
View File

@@ -0,0 +1,96 @@
30 41 100.0 28 0.1
1.e-5 2
1 1
0
1 1 2 0.0192 0.0575 0.0264
2 1 3 0.0452 0.1852 0.0204
3 2 4 0.0570 0.1737 0.0184
4 3 4 0.0132 0.0379 0.0042
5 2 5 0.0472 0.1983 0.0209
6 2 6 0.0581 0.1763 0.0187
7 4 6 0.0119 0.0414 0.0045
8 5 7 0.0460 0.1160 0.0102
9 6 7 0.0267 0.0820 0.0085
10 6 8 0.0120 0.0420 0.0045
13 9 11 0.0 0.2080 0.0
15 12 13 0.0 0.1400 0.0
16 12 14 0.1231 0.2559 0.0
17 12 15 0.0662 0.1304 0.0
18 12 16 0.945 0.1987 0.0
19 14 15 0.2210 0.1997 0.0
20 16 17 0.0824 0.1923 0.0
21 15 18 0.1070 0.2185 0.0
22 18 19 0.0639 0.1292 0.0
23 19 20 0.0340 0.0680 0.0
24 10 20 0.0936 0.2090 0.0
25 10 17 0.0324 0.0845 0.0
26 10 21 0.0348 0.0749 0.0
27 10 22 0.0727 0.1499 0.0
28 21 22 0.0116 0.0236 0.0
29 15 23 0.1000 0.2020 0.0
30 22 24 0.1150 0.1790 0.0
31 23 24 0.1320 0.2700 0.0
32 24 25 0.1885 0.3292 0.0
33 25 26 0.2554 0.3800 0.0
34 25 27 0.1093 0.2087 0.0
36 27 29 0.2198 0.4153 0.0
37 27 30 0.3202 0.6027 0.0
38 29 30 0.2399 0.4533 0.0
39 8 28 0.0636 0.2000 0.0214
40 6 28 0.0169 0.0599 0.0065
41 9 10 0.0 0.1100 0.0
0
10 0.19
24 0.043
0
1 9 6 0.0 0.2080 0.978 0.9 1.1
2 6 10 0.0 0.5560 0.969 0.9 1.1
3 12 4 0.0 0.2560 0.932 0.9 1.1
4 28 27 0.0 0.3960 0.968 0.9 1.1
0
1 20. 0. 0. 0.
2 57.56 2.43 21.7 12.7
3 0. 0. 2.4 1.2
4 0. 0. 7.6 1.6
5 24.56 22.25 94.2 19.
6 0. 0. 0. 0.
7 0. 0. 22.8 10.9
8 35 37.27 30. 30.
9 0. 0. 0. 0.
10 0. 0. 5.8 2.
11 17.93 17.61 0. 0.
12 0. 0. 11.2 7.5
13 16.91 24.69 0. 0.
14 0. 0. 6.2 1.6
15 0. 0. 8.2 2.5
16 0. 0. 3.5 1.8
17 0. 0. 9. 5.8
18 0. 0. 3.2 .9
19 0. 0. 9.5 3.4
20 0. 0. 2.2 .7
21 0. 0. 17.5 11.2
22 0. 0. 0. 0.
23 0. 0. 3.2 1.6
24 0. 0. 8.7 6.7
25 0. 0. 0. 0.
26 0. 0. 3.5 2.3
27 0. 0. 0. 0.
28 0. 0. 0. 0.
29 0. 0. 2.4 .9
30 0. 0. 10.6 1.9
0
1 1.060 -50 50.
2 1.045 -40. 60.
5 1.010 -40. 40.
8 1.010 -10. 40.
11 1.082 -6. 24.
13 1.071 -6. 24.
0
1 10. 2.0 2.0 10. 60.
2 10. 1.5 2.4 10. 60.
5 20. 1.8 0.8 10. 150.
8 10. 1.0 1.2 10. 120.
11 20. 1.8 0.8 10. 150.
13 10. 1.5 2.0 10. 60.
0
0

824
IEEE300.dat Normal file
View File

@@ -0,0 +1,824 @@
300 409 100. 28 0.1
1.e-5 4
1 38
0
1 269 291 0.00080 0.00348 0.00000
2 226 271 0.05558 0.24666 0.00000
3 226 300 0.05559 0.24666 0.00000
4 227 225 0.03811 0.21648 0.00000
5 225 228 0.05370 0.07026 0.00000
6 228 229 1.10680 0.95278 0.00000
7 271 300 0.05580 0.24666 0.00000
8 300 144 0.07378 0.06352 0.00000
9 144 270 0.03832 0.02894 0.00000
10 227 68 0.23552 0.99036 0.00000
11 146 147 0.00100 0.00600 0.00000
12 230 71 0.00100 0.00900 0.00000
13 230 148 0.00600 0.02700 0.05400
14 292 272 0.00000 0.00300 0.00000
15 292 150 0.00800 0.06900 0.13900
16 292 104 0.00100 0.00700 0.00000
17 70 149 0.00200 0.01900 1.12700
18 147 72 0.00600 0.02900 0.01800
19 272 231 0.00100 0.00900 0.07000
20 272 98 0.00100 0.00700 0.01400
21 148 273 0.01300 0.05950 0.03300
22 148 75 0.01300 0.04200 0.08100
23 72 273 0.00600 0.02700 0.01300
24 273 74 0.00800 0.03400 0.01800
25 231 233 0.00200 0.01500 0.11800
26 74 232 0.00600 0.03400 0.01600
27 75 286 0.01400 0.04200 0.09700
28 286 297 0.06500 0.24800 0.12100
29 286 165 0.09900 0.24800 0.03500
30 286 166 0.09600 0.36300 0.04800
31 149 274 0.00200 0.02200 1.28000
32 150 233 0.00200 0.01800 0.03600
33 150 163 0.01300 0.08000 0.15100
34 232 77 0.01600 0.03300 0.01500
35 232 79 0.06900 0.18600 0.09800
36 233 235 0.00400 0.03400 0.28000
37 77 234 0.05200 0.11100 0.05000
38 234 78 0.01900 0.03900 0.01800
39 235 14 0.00700 0.06800 0.13400
40 78 151 0.03600 0.07100 0.03400
41 151 79 0.04500 0.12000 0.06500
42 151 15 0.04300 0.13000 0.01400
43 236 80 0.00000 0.06300 0.00000
44 236 238 0.00250 0.01200 0.01300
45 236 152 0.00600 0.02900 0.02000
46 236 287 0.00700 0.04300 0.02600
47 80 274 0.00100 0.00800 0.04200
48 237 245 0.01200 0.06000 0.00800
49 237 161 0.00600 0.01400 0.00200
50 237 293 0.01000 0.02900 0.00300
51 81 164 0.00400 0.02700 0.04300
52 297 238 0.00800 0.04700 0.00800
53 297 152 0.02200 0.06400 0.00700
54 297 287 0.01000 0.03600 0.02000
55 297 241 0.01700 0.08100 0.04800
56 297 165 0.10200 0.25400 0.03300
57 297 166 0.04700 0.12700 0.01600
58 238 287 0.00800 0.03700 0.02000
59 238 239 0.03200 0.08700 0.04000
60 82 274 0.00060 0.00640 0.40400
61 152 155 0.02600 0.15400 0.02200
62 287 274 0.00000 0.02900 0.00000
63 287 241 0.06500 0.19100 0.02000
64 287 156 0.03100 0.08900 0.03600
65 274 153 0.00200 0.01400 0.80600
66 239 275 0.02600 0.07200 0.03500
67 239 155 0.09500 0.26200 0.03200
68 239 84 0.01300 0.03900 0.01600
69 275 154 0.02700 0.08400 0.03900
70 275 157 0.02800 0.08400 0.03700
71 240 87 0.00700 0.04100 0.31200
72 240 246 0.00900 0.05400 0.41100
73 153 248 0.00500 0.04200 0.69000
74 154 277 0.05200 0.14500 0.07300
75 154 94 0.04300 0.11800 0.01300
76 155 173 0.02500 0.06200 0.00700
77 241 156 0.03100 0.09400 0.04300
78 156 83 0.03700 0.10900 0.04900
79 83 242 0.02700 0.08000 0.03600
80 84 157 0.02500 0.07300 0.03500
81 157 242 0.03500 0.10300 0.04700
82 242 243 0.06500 0.16900 0.08200
83 243 85 0.04600 0.08000 0.03600
84 243 159 0.15900 0.53700 0.07100
85 85 86 0.00900 0.02600 0.00500
86 86 158 0.00200 0.01300 0.01500
87 87 276 0.00900 0.06500 0.48500
88 276 88 0.01600 0.10500 0.20300
89 276 101 0.00100 0.00700 0.01300
90 159 19 0.02650 0.17200 0.02600
91 160 298 0.05100 0.23200 0.02800
92 160 247 0.05100 0.15700 0.02300
93 89 244 0.03200 0.10000 0.06200
94 89 20 0.02000 0.12340 0.02800
95 244 245 0.03600 0.13100 0.06800
96 244 277 0.03400 0.09900 0.04700
97 245 293 0.01800 0.08700 0.01100
98 245 21 0.02560 0.19300 0.00000
99 277 161 0.02100 0.05700 0.03000
100 277 247 0.01800 0.05200 0.01800
101 246 164 0.00400 0.02700 0.05000
102 246 23 0.02860 0.20130 0.37900
103 161 293 0.01600 0.04300 0.00400
104 293 162 0.00100 0.00600 0.00700
105 293 90 0.01400 0.07000 0.03800
106 293 22 0.08910 0.26760 0.02900
107 293 24 0.07820 0.21270 0.02200
108 162 247 0.00600 0.02200 0.01100
109 162 1 0.00000 0.03600 0.00000
110 247 298 0.09900 0.37500 0.05100
111 90 298 0.02200 0.10700 0.05800
112 248 205 0.00350 0.03300 0.53000
113 248 206 0.00350 0.03300 0.53000
114 91 249 0.00800 0.06400 0.12800
115 249 163 0.01200 0.09300 0.18300
116 249 17 0.00600 0.04800 0.09200
117 165 167 0.04700 0.11900 0.01400
118 166 168 0.03200 0.17400 0.02400
119 167 169 0.10000 0.25300 0.03100
120 167 278 0.02200 0.07700 0.03900
121 168 171 0.01900 0.14400 0.01700
122 168 250 0.01700 0.09200 0.01200
123 169 278 0.27800 0.42700 0.04300
124 278 170 0.02200 0.05300 0.00700
125 278 280 0.03800 0.09200 0.01200
126 278 171 0.04800 0.12200 0.01500
127 92 170 0.02400 0.06400 0.00700
128 92 280 0.03400 0.12100 0.01500
129 279 173 0.05300 0.13500 0.01700
130 279 174 0.00200 0.00400 0.00200
131 279 251 0.04500 0.35400 0.04400
132 279 252 0.05000 0.17400 0.02200
133 170 280 0.01600 0.03800 0.00400
134 280 172 0.04300 0.06400 0.02700
135 171 250 0.01900 0.06200 0.00800
136 172 174 0.07600 0.13000 0.04400
137 172 16 0.04400 0.12400 0.01500
138 250 173 0.01200 0.08800 0.01100
139 250 252 0.15700 0.40000 0.04700
140 174 18 0.07400 0.20800 0.02600
141 251 252 0.07000 0.18400 0.02100
142 251 94 0.10000 0.27400 0.03100
143 251 175 0.10900 0.39300 0.03600
144 252 93 0.14200 0.40400 0.05000
145 93 175 0.01700 0.04200 0.00600
146 95 256 0.00360 0.01990 0.00400
147 96 255 0.00200 0.10490 0.00100
148 97 253 0.00010 0.00180 0.01700
149 253 254 0.00000 0.02710 0.00000
150 253 142 0.00000 0.61630 0.00000
151 142 255 0.00000 -0.36970 0.00000
152 253 176 0.00220 0.29150 0.00000
153 254 255 0.00000 0.03390 0.00000
154 254 176 0.00000 0.05820 0.00000
155 256 177 0.08080 0.23440 0.02900
156 256 179 0.09650 0.36690 0.05400
157 177 178 0.03600 0.10760 0.11700
158 177 179 0.04760 0.14140 0.14900
159 179 294 0.00060 0.01970 0.00000
160 294 257 0.00590 0.04050 0.25000
161 294 181 0.01150 0.11060 0.18500
162 294 182 0.01980 0.16880 0.32100
163 294 191 0.00500 0.05000 0.33000
164 294 192 0.00770 0.05380 0.33500
165 294 196 0.01650 0.11570 0.17100
166 257 180 0.00590 0.05770 0.09500
167 257 183 0.00490 0.03360 0.20800
168 257 195 0.00590 0.05770 0.09500
169 180 299 0.00780 0.07730 0.12600
170 180 288 0.00260 0.01930 0.03000
171 181 299 0.00760 0.07520 0.12200
172 181 288 0.00210 0.01860 0.03000
173 299 182 0.00160 0.01640 0.02600
174 299 105 0.00170 0.01650 0.02600
175 299 115 0.00790 0.07930 0.12700
176 299 195 0.00780 0.07840 0.12500
177 288 295 0.00170 0.01170 0.28900
178 288 195 0.00260 0.01930 0.03000
179 288 196 0.00210 0.01860 0.03000
180 288 2 0.00020 0.01010 0.00000
181 183 99 0.00430 0.02930 0.18000
182 183 121 0.00390 0.03810 0.25800
183 99 184 0.00910 0.06230 0.38500
184 184 295 0.01250 0.08900 0.54000
185 184 106 0.00560 0.03900 0.95300
186 295 296 0.00150 0.01140 0.28400
187 295 201 0.00050 0.00340 0.02100
188 295 122 0.00070 0.01510 0.12600
189 295 262 0.00050 0.00340 0.02100
190 185 197 0.05620 0.22480 0.08100
191 296 186 0.01200 0.08360 0.12300
192 296 187 0.01520 0.11320 0.68400
193 296 282 0.04680 0.33690 0.51900
194 296 258 0.04300 0.30310 0.46300
195 296 102 0.04890 0.34920 0.53800
196 296 119 0.00130 0.00890 0.11900
197 186 258 0.02910 0.22670 0.34200
198 187 281 0.00600 0.05700 0.76700
199 281 282 0.00750 0.07730 0.11900
200 281 103 0.01270 0.09090 0.13500
201 282 258 0.00850 0.05880 0.08700
202 282 103 0.02180 0.15110 0.22300
203 258 102 0.00730 0.05040 0.07400
204 188 261 0.05230 0.15260 0.07400
205 188 200 0.13710 0.39190 0.07600
206 106 189 0.01370 0.09570 0.14100
207 189 110 0.00550 0.02880 0.19000
208 107 108 0.17460 0.31610 0.04000
209 107 120 0.08040 0.30540 0.04500
210 190 110 0.01100 0.05680 0.38800
211 191 193 0.00080 0.00980 0.06900
212 192 193 0.00290 0.02850 0.19000
213 192 109 0.00660 0.04480 0.27700
214 111 194 0.00240 0.03260 0.23600
215 111 113 0.00180 0.02450 1.66200
216 112 194 0.00440 0.05140 3.59700
217 113 114 0.00020 0.01230 0.00000
218 115 196 0.00180 0.01780 0.02900
219 197 259 0.06690 0.48430 0.06300
220 197 198 0.05580 0.22100 0.03100
221 259 198 0.08070 0.33310 0.04900
222 259 260 0.07390 0.30710 0.04300
223 259 199 0.17990 0.50170 0.06900
224 260 199 0.09040 0.36260 0.04800
225 260 200 0.07700 0.30920 0.05400
226 199 117 0.02510 0.08290 0.04700
227 117 261 0.02220 0.08470 0.05000
228 261 200 0.04980 0.18550 0.02900
229 261 118 0.00610 0.02900 0.08400
230 201 100 0.00040 0.02020 0.00000
231 201 123 0.00040 0.00830 0.11500
232 121 3 0.00250 0.02450 0.16400
233 122 262 0.00070 0.00860 0.11500
234 123 262 0.00070 0.00860 0.11500
235 262 100 0.00040 0.02020 0.00000
236 202 212 0.03300 0.09500 0.00000
237 202 131 0.04600 0.06900 0.00000
238 203 290 0.00040 0.00220 6.20000
239 203 138 0.00000 0.02750 0.00000
240 124 125 0.00300 0.04800 0.00000
241 125 218 0.00200 0.00900 0.00000
242 204 210 0.04500 0.06300 0.00000
243 204 212 0.04800 0.12700 0.00000
244 205 284 0.00310 0.02860 0.50000
245 205 25 0.00240 0.03550 0.36000
246 206 284 0.00310 0.02860 0.50000
247 263 207 0.01400 0.04000 0.00400
248 263 283 0.03000 0.08100 0.01000
249 207 289 0.01000 0.06000 0.00900
250 207 298 0.01500 0.04000 0.00600
251 289 128 0.33200 0.68800 0.00000
252 289 129 0.00900 0.04600 0.02500
253 289 283 0.02000 0.07300 0.00800
254 289 298 0.03400 0.10900 0.03200
255 126 208 0.07600 0.13500 0.00900
256 126 283 0.04000 0.10200 0.00500
257 208 283 0.08100 0.12800 0.01400
258 127 209 0.12400 0.18300 0.00000
259 129 298 0.01000 0.05900 0.00800
260 209 210 0.04600 0.06800 0.00000
261 210 211 0.30200 0.44600 0.00000
262 211 130 0.07300 0.09300 0.00000
263 211 212 0.24000 0.42100 0.00000
264 213 215 0.01390 0.07780 0.08600
265 214 215 0.00170 0.01850 0.02000
266 214 222 0.00150 0.01080 0.00200
267 215 132 0.00450 0.02490 0.02600
268 132 264 0.00400 0.04970 0.01800
269 264 216 0.00000 0.04560 0.00000
270 264 284 0.00050 0.01770 0.02000
271 264 265 0.00270 0.03950 0.83200
272 284 285 0.00030 0.00180 5.20000
273 265 216 0.00370 0.04840 0.43000
274 265 133 0.00100 0.02950 0.50300
275 265 221 0.00160 0.00460 0.40200
276 133 134 0.00030 0.00130 1.00000
277 217 218 0.01000 0.06400 0.48000
278 217 135 0.00190 0.00810 0.86000
279 218 124 0.00100 0.06100 0.00000
280 135 290 0.00050 0.02120 0.00000
281 219 220 0.00190 0.00870 1.28000
282 219 290 0.00260 0.09170 0.00000
283 219 266 0.00130 0.02880 0.81000
284 220 203 0.00000 0.06260 0.00000
285 290 136 0.00020 0.00690 1.36400
286 290 285 0.00010 0.00060 3.57000
287 136 8 0.00170 0.04850 0.00000
288 266 137 0.00020 0.02590 0.14400
289 266 285 0.00060 0.02720 0.00000
290 137 221 0.00020 0.00060 0.80000
291 138 13 0.00030 0.00430 0.00900
292 222 267 0.00820 0.08510 0.00000
293 222 268 0.01120 0.07230 0.00000
294 139 140 0.01270 0.03550 0.00000
295 139 267 0.03260 0.18040 0.00000
296 140 223 0.01950 0.05510 0.00000
297 267 223 0.01570 0.07320 0.00000
298 267 268 0.03600 0.21190 0.00000
299 223 268 0.02680 0.12850 0.00000
300 268 224 0.04280 0.12150 0.00000
301 224 141 0.03510 0.10040 0.00000
302 141 12 0.06160 0.18570 0.00000
0
97 3.250
255 0.550
107 0.345
194 -2.120
114 -1.030
259 0.530
200 0.450
203 -1.500
290 -3.000
221 -1.500
138 -1.400
224 0.456
300 0.024
54 0.017
0
1 297 269 0.0001 0.0005 1.0082 0.9043 1.1043
2 269 226 0.0244 0.4368 0.9668 0.9391 1.1478
3 269 227 0.0362 0.6490 0.9796 0.9391 1.1478
4 291 62 0.0158 0.3749 1.0435 0.9391 1.1478
5 291 63 0.0158 0.3749 0.9391 0.9391 1.1478
6 291 145 0.0160 0.3805 1.0435 0.9391 1.1478
7 291 64 0.0000 0.1520 1.0435 0.9391 1.1000
8 291 65 0.0000 0.8000 1.0435 0.9391 1.1000
9 228 47 0.4436 2.8152 1.0000 0.9391 1.1000
10 225 48 0.5075 3.2202 1.0000 0.9391 1.1000
11 229 49 0.6669 3.9440 1.0000 0.9391 1.1000
12 229 50 0.6113 3.6152 1.0000 0.9391 1.1000
13 271 66 0.4412 2.9668 1.0000 0.9391 1.1000
14 271 67 0.3079 2.0570 1.0000 0.9391 1.1000
15 300 51 0.7363 4.6724 1.0000 0.9391 1.1000
16 300 52 0.7698 4.8846 1.0000 0.9391 1.1000
17 300 53 0.7573 4.8056 1.0000 0.9391 1.1000
18 270 59 0.3661 2.4560 1.0000 0.9391 1.1000
19 270 60 1.0593 5.4536 1.0000 0.9391 1.1000
20 270 61 0.1567 1.6994 1.0000 0.9000 1.1000
21 300 54 0.1301 1.3912 1.0000 0.9391 1.1000
22 300 55 0.5448 3.4572 1.0000 0.9391 1.1000
23 300 56 0.1543 1.6729 1.0000 0.9391 1.1000
24 300 57 0.3849 2.5712 1.0000 0.9391 1.1000
25 300 58 0.4412 2.9668 1.0000 0.9391 1.1000
26 145 69 0.0000 0.7500 0.9583 0.9391 1.1000
27 4 214 0.0025 0.0380 1.0000 0.9391 1.1000
28 5 285 0.0014 0.0514 1.0000 0.9391 1.1000
29 6 290 0.0009 0.0472 1.0000 0.9391 1.1000
30 11 285 0.0005 0.0154 1.0000 0.9391 1.1000
31 292 146 0.0000 0.0520 0.9470 0.9000 1.1000
32 292 230 0.0000 0.0520 0.9560 0.9000 1.1000
33 292 70 0.0000 0.0050 0.9710 0.9000 1.1000
34 272 147 0.0000 0.0390 0.9480 0.9000 1.1000
35 272 71 0.0000 0.0390 0.9590 0.9000 1.1000
36 73 273 0.0000 0.0890 1.0460 0.9000 1.1000
37 231 73 0.0000 0.0530 0.9850 0.9000 1.1000
38 286 76 0.0194 0.0311 0.9561 0.9000 1.1000
39 149 286 0.0010 0.0380 0.9710 0.9000 1.1000
40 233 232 0.0000 0.0140 0.9520 0.9000 1.1000
41 235 234 0.0000 0.0640 0.9430 0.9000 1.1000
42 81 237 0.0000 0.0470 1.0100 0.9000 1.1000
43 240 275 0.0000 0.0200 1.0080 0.9000 1.1000
44 240 153 0.0000 0.0210 1.0000 0.9000 1.1000
45 276 158 0.0000 0.0590 0.9750 0.9000 1.1000
46 159 88 0.0000 0.0380 1.0170 0.9000 1.1000
47 277 246 0.0000 0.0244 1.0000 0.9000 1.1000
48 248 164 0.0000 0.0200 1.0000 0.9000 1.1000
49 91 279 0.0000 0.0480 1.0000 0.9000 1.1000
50 249 280 0.0000 0.0480 1.0000 0.9000 1.1000
51 163 169 0.0000 0.0460 1.0150 0.9000 1.1000
52 175 130 0.0000 0.1490 0.9670 0.9000 1.1000
53 96 178 0.0052 0.0174 1.0100 0.9000 1.1000
54 176 95 0.0000 0.0280 1.0500 0.9000 1.1000
55 256 191 0.0005 0.0195 1.0000 0.9000 1.1000
56 299 98 0.0000 0.0180 1.0522 0.9000 1.1000
57 299 104 0.0000 0.0140 1.0522 0.9000 1.1000
58 182 116 0.0010 0.0402 1.0500 0.9000 1.1000
59 186 198 0.0024 0.0603 0.9750 0.9000 1.1000
60 187 260 0.0024 0.0498 1.0000 0.9000 1.1000
61 281 101 0.0000 0.0833 1.0350 0.9000 1.1000
62 281 188 0.0013 0.0371 0.9565 0.9000 1.1000
63 282 118 0.0005 0.0182 1.0000 0.9000 1.1000
64 105 116 0.0010 0.0392 1.0500 0.9000 1.1000
65 189 120 0.0027 0.0639 1.0730 0.9000 1.1000
66 190 108 0.0008 0.0256 1.0500 0.9000 1.1000
67 193 97 0.0000 0.0160 1.0506 0.9000 1.1000
68 109 178 0.0012 0.0396 0.9750 0.9000 1.1000
69 112 295 0.0013 0.0384 0.9800 0.9000 1.1000
70 194 190 0.0009 0.0231 0.9560 0.9000 1.1000
71 119 185 0.0003 0.0131 1.0500 0.9000 1.1000
72 202 283 0.0000 0.2520 1.0300 0.9000 1.1000
73 204 263 0.0000 0.2370 1.0300 0.9000 1.1000
74 206 213 0.0008 0.0366 0.9850 0.9000 1.1000
75 208 224 0.0000 0.2200 1.0000 0.9000 1.1000
76 127 160 0.0000 0.0980 1.0300 0.9000 1.1000
77 128 298 0.0000 0.1280 1.0100 0.9000 1.1000
78 209 143 0.0200 0.2040 1.0500 0.9000 1.1000
79 131 289 0.0260 0.2110 1.0300 0.9000 1.1000
80 298 213 0.0030 0.0122 1.0000 0.9000 1.1000
81 216 284 0.0030 0.0122 0.9700 0.9000 1.1000
82 134 217 0.0012 0.0195 1.0000 0.9000 1.1000
83 220 7 0.0010 0.0332 1.0200 0.9000 1.1000
84 266 9 0.0005 0.0160 1.0700 0.9000 1.1000
85 221 10 0.0005 0.0160 1.0200 0.9000 1.1000
86 263 143 0.0001 0.0200 1.0000 0.9000 1.1000
87 254 26 0.0010 0.0230 1.0223 0.9000 1.1000
88 255 27 0.0000 0.0230 0.9284 0.9000 1.1000
89 29 230 0.0010 0.0146 1.0000 0.9000 1.1000
90 30 292 0.0000 0.0105 1.0000 0.9000 1.1000
91 41 158 0.0000 0.0238 1.0000 0.9000 1.1000
92 42 276 0.0000 0.0321 0.9500 0.9000 1.1000
93 46 114 0.0000 0.0154 1.0000 0.9000 1.1000
94 35 235 0.0000 0.0289 1.0000 0.9000 1.1000
95 28 146 0.0000 0.0195 1.0000 0.9000 1.1000
96 44 299 0.0000 0.0193 1.0000 0.9000 1.1000
97 31 273 0.0000 0.0192 1.0000 0.9000 1.1000
98 34 234 0.0000 0.0230 1.0000 0.9000 1.1000
99 38 241 0.0000 0.0124 1.0000 0.9000 1.1000
100 45 185 0.0000 0.0167 1.0000 0.9000 1.1000
101 32 231 0.0000 0.0312 1.0000 0.9000 1.1000
102 33 76 0.0000 0.0165 0.9420 0.9000 1.1000
103 36 82 0.0000 0.0316 0.9650 0.9000 1.1000
104 40 243 0.0000 0.0535 0.9500 0.9000 1.1000
105 37 275 0.0000 0.1818 0.9420 0.9000 1.1000
106 39 242 0.0000 0.1961 0.9420 0.9000 1.1000
107 43 244 0.0000 0.0690 0.9565 0.9000 1.1000
0
146 0.00 0.00 90.00 49.00
230 0.00 0.00 56.00 15.00
292 0.00 0.00 20.00 0.00
70 0.00 0.00 0.00 0.00
147 0.00 0.00 353.00 130.00
71 0.00 0.00 120.00 41.00
272 0.00 0.00 0.00 0.00
148 -5.00 0.00 58.00 14.00
72 0.00 0.00 96.00 43.00
73 -5.00 0.00 148.00 33.00
273 0.00 0.00 83.00 21.00
231 0.00 0.00 0.00 0.00
74 0.00 0.00 58.00 10.00
75 0.00 0.00 160.00 60.00
286 0.00 0.00 126.70 23.00
149 0.00 0.00 0.00 0.00
76 0.00 0.00 561.00 220.00
150 0.00 0.00 0.00 0.00
232 -10.00 0.00 595.00 120.00
233 0.00 0.00 77.00 1.00
77 0.00 0.00 81.00 23.00
234 0.00 0.00 21.00 7.00
235 0.00 0.00 0.00 0.00
78 0.00 0.00 45.00 12.00
151 0.00 0.00 28.00 9.00
79 0.00 0.00 69.00 13.00
236 0.00 0.00 55.00 6.00
80 0.00 0.00 0.00 0.00
237 0.00 0.00 0.00 0.00
81 0.00 0.00 0.00 0.00
297 0.00 0.00 85.00 32.00
238 0.00 0.00 155.00 18.00
82 0.00 0.00 0.00 0.00
152 0.00 0.00 46.00 -21.00
287 0.00 0.00 86.00 0.00
274 0.00 0.00 0.00 0.00
239 0.00 0.00 39.00 9.00
275 0.00 0.00 195.00 29.00
240 0.00 0.00 0.00 0.00
153 0.00 0.00 0.00 0.00
154 0.00 0.00 58.00 11.80
155 0.00 0.00 41.00 19.00
241 0.00 0.00 92.00 26.00
156 0.00 0.00 -5.00 5.00
83 0.00 0.00 61.00 28.00
84 0.00 0.00 69.00 3.00
157 0.00 0.00 10.00 1.00
242 0.00 0.00 22.00 10.00
243 0.00 0.00 98.00 20.00
85 0.00 0.00 14.00 1.00
86 0.00 0.00 218.00 106.00
87 0.00 0.00 0.00 0.00
158 0.00 0.00 227.00 110.00
276 0.00 0.00 0.00 0.00
159 0.00 0.00 70.00 30.00
88 0.00 0.00 0.00 0.00
160 0.00 0.00 0.00 0.00
89 0.00 0.00 56.00 20.00
244 0.00 0.00 116.00 38.00
245 0.00 0.00 57.00 19.00
277 0.00 0.00 224.00 71.00
246 0.00 0.00 0.00 0.00
161 0.00 0.00 208.00 107.00
293 0.00 0.00 74.00 28.00
162 0.00 0.00 0.00 0.00
247 0.00 0.00 48.00 14.00
90 0.00 0.00 28.00 7.00
248 0.00 0.00 0.00 0.00
1 375.00 0.00 37.00 13.00
91 0.00 0.00 0.00 0.00
249 0.00 0.00 0.00 0.00
163 0.00 0.00 0.00 0.00
164 0.00 0.00 0.00 0.00
165 0.00 0.00 44.20 0.00
166 0.00 0.00 66.00 0.00
167 155.00 0.00 17.40 0.00
168 290.00 0.00 15.80 0.00
169 0.00 0.00 60.30 0.00
278 0.00 0.00 39.90 0.00
92 68.00 0.00 66.70 0.00
279 0.00 0.00 83.50 0.00
170 0.00 0.00 0.00 0.00
280 0.00 0.00 77.80 0.00
171 0.00 0.00 32.00 0.00
172 0.00 0.00 8.60 0.00
250 0.00 0.00 49.60 0.00
173 0.00 0.00 4.60 0.00
174 117.00 0.00 112.10 0.00
251 0.00 0.00 30.70 0.00
252 0.00 0.00 63.00 0.00
93 0.00 0.00 19.60 0.00
94 0.00 0.00 26.20 0.00
175 0.00 0.00 18.20 0.00
95 0.00 0.00 0.00 0.00
96 0.00 0.00 0.00 0.00
97 0.00 0.00 0.00 0.00
253 0.00 0.00 14.10 650.00
254 1930.00 0.00 0.00 0.00
255 0.00 0.00 777.00 215.00
176 0.00 0.00 535.00 55.00
256 0.00 0.00 229.10 11.80
177 0.00 0.00 78.00 1.40
178 240.00 0.00 276.40 59.30
179 0.00 0.00 514.80 82.70
294 0.00 0.00 57.90 5.10
257 0.00 0.00 380.80 37.00
180 0.00 0.00 0.00 0.00
181 0.00 0.00 0.00 0.00
299 0.00 0.00 0.00 0.00
98 0.00 0.00 0.00 0.00
182 0.00 0.00 0.00 0.00
288 0.00 0.00 0.00 0.00
183 0.00 0.00 0.00 0.00
99 0.00 0.00 169.20 41.60
184 0.00 0.00 55.20 18.20
295 0.00 0.00 273.60 99.80
100 -192.50 0.00 826.70 135.20
185 0.00 0.00 595.00 83.30
296 0.00 0.00 387.70 114.70
186 281.00 0.00 145.00 58.00
187 0.00 0.00 56.50 24.50
281 696.00 0.00 89.50 35.50
101 0.00 0.00 0.00 0.00
282 0.00 0.00 24.00 14.00
258 84.00 0.00 0.00 0.00
102 217.00 0.00 0.00 0.00
188 0.00 0.00 63.00 25.00
103 103.00 0.00 0.00 0.00
104 0.00 0.00 0.00 0.00
105 0.00 0.00 0.00 0.00
106 372.00 0.00 17.00 9.00
189 216.00 0.00 0.00 0.00
107 0.00 0.00 70.00 5.00
190 0.00 0.00 200.00 50.00
108 0.00 0.00 75.00 50.00
191 0.00 0.00 123.50 -24.30
192 0.00 0.00 0.00 0.00
193 0.00 0.00 33.00 16.50
109 0.00 0.00 0.00 0.00
110 0.00 0.00 35.00 15.00
111 0.00 0.00 85.00 24.00
112 0.00 0.00 0.00 0.40
194 0.00 0.00 0.00 0.00
113 0.00 0.00 0.00 0.00
114 0.00 0.00 0.00 0.00
115 0.00 0.00 299.90 95.70
195 0.00 0.00 0.00 0.00
196 0.00 0.00 0.00 0.00
116 205.00 0.00 481.80 205.00
2 0.00 0.00 763.60 291.10
197 0.00 0.00 26.50 0.00
259 0.00 0.00 163.50 43.00
198 0.00 0.00 0.00 0.00
260 0.00 0.00 176.00 83.00
199 228.00 0.00 5.00 4.00
117 84.00 0.00 28.00 12.00
261 0.00 0.00 427.40 173.60
200 0.00 0.00 74.00 29.00
118 0.00 0.00 69.50 49.30
201 0.00 0.00 73.40 0.00
119 0.00 0.00 240.70 89.00
120 0.00 0.00 40.00 4.00
121 0.00 0.00 136.80 16.60
3 200.00 0.00 0.00 0.00
122 1200.00 0.00 59.80 24.30
123 1200.00 0.00 59.80 24.30
262 0.00 0.00 182.60 43.60
202 0.00 0.00 7.00 2.00
203 475.00 0.00 0.00 0.00
124 1973.00 0.00 489.00 53.00
125 0.00 0.00 800.00 72.00
204 0.00 0.00 0.00 0.00
205 0.00 0.00 0.00 0.00
206 0.00 0.00 0.00 0.00
263 0.00 0.00 10.00 3.00
207 0.00 0.00 43.00 14.00
289 424.00 0.00 64.00 21.00
126 0.00 0.00 35.00 12.00
208 0.00 0.00 27.00 12.00
127 0.00 0.00 41.00 14.00
128 0.00 0.00 38.00 13.00
129 0.00 0.00 42.00 14.00
209 0.00 0.00 72.00 24.00
210 0.00 0.00 0.00 -5.00
211 0.00 0.00 12.00 2.00
130 0.00 0.00 -21.00 -14.20
212 0.00 0.00 7.00 2.00
131 0.00 0.00 38.00 13.00
283 0.00 0.00 0.00 0.00
298 0.00 0.00 96.00 7.00
213 0.00 0.00 0.00 0.00
4 272.00 0.00 0.00 0.00
214 0.00 0.00 22.00 16.00
215 0.00 0.00 47.00 26.00
132 0.00 0.00 176.00 105.00
264 0.00 0.00 100.00 75.00
216 0.00 0.00 131.00 96.00
284 0.00 0.00 0.00 0.00
265 100.00 0.00 285.00 100.00
133 450.00 0.00 171.00 70.00
5 250.00 0.00 328.00 188.00
134 0.00 0.00 428.00 232.00
217 0.00 0.00 173.00 99.00
218 0.00 0.00 410.00 40.00
135 0.00 0.00 0.00 0.00
6 303.00 0.00 538.00 369.00
219 0.00 0.00 223.00 148.00
220 0.00 0.00 96.00 46.00
7 345.00 0.00 0.00 0.00
290 0.00 0.00 159.00 107.00
136 0.00 0.00 448.00 143.00
8 300.00 0.00 404.00 212.00
266 0.00 0.00 572.00 244.00
137 0.00 0.00 269.00 157.00
9 600.00 0.00 0.00 0.00
285 0.00 0.00 0.00 0.00
221 250.00 0.00 255.00 149.00
10 550.00 0.00 0.00 0.00
138 0.00 0.00 0.00 0.00
11 575.43 0.00 0.00 0.00
222 170.00 0.00 0.00 0.00
139 84.00 0.00 8.00 3.00
140 0.00 0.00 0.00 0.00
267 0.00 0.00 61.00 30.00
223 0.00 0.00 77.00 33.00
268 0.00 0.00 61.00 30.00
224 0.00 0.00 29.00 14.00
141 0.00 0.00 29.00 14.00
12 0.00 0.00 -23.00 -17.00
13 0.00 0.00 -33.10 -29.40
14 0.00 0.00 115.80 -24.00
15 0.00 0.00 2.40 -12.60
16 0.00 0.00 2.40 -3.90
17 0.00 0.00 -14.90 26.50
18 0.00 0.00 24.70 -1.20
19 0.00 0.00 145.30 -34.90
20 0.00 0.00 28.10 -20.50
21 0.00 0.00 14.00 2.50
22 0.00 0.00 -11.10 -1.40
23 0.00 0.00 50.50 17.40
24 0.00 0.00 29.60 0.60
25 0.00 0.00 -113.70 76.70
26 0.00 0.00 100.31 29.17
27 0.00 0.00 -100.00 34.17
142 0.00 0.00 0.00 0.00
143 0.00 0.00 0.00 0.00
28 467.00 0.00 0.00 0.00
29 623.00 0.00 0.00 0.00
30 1210.00 0.00 0.00 0.00
31 234.00 0.00 0.00 0.00
32 372.00 0.00 0.00 0.00
33 330.00 0.00 0.00 0.00
34 185.00 0.00 0.00 0.00
35 410.00 0.00 0.00 0.00
36 500.00 0.00 0.00 0.00
37 37.00 0.00 0.00 0.00
38 0.00 0.00 0.00 0.00
39 45.00 0.00 0.00 0.00
40 165.00 0.00 0.00 0.00
41 400.00 0.00 0.00 0.00
42 400.00 0.00 0.00 0.00
43 116.00 0.00 0.00 0.00
44 1292.00 0.00 0.00 0.00
45 700.00 0.00 0.00 0.00
46 553.00 0.00 0.00 0.00
269 0.00 0.00 0.00 0.00
225 -4.20 0.00 0.00 0.00
300 0.00 0.00 2.71 0.94
270 0.00 0.00 0.86 0.28
291 0.00 0.00 0.00 0.00
226 0.00 0.00 0.00 0.00
271 0.00 0.00 0.00 0.00
227 0.00 0.00 0.00 0.00
228 0.00 0.00 4.75 1.56
47 0.00 0.00 1.53 0.53
229 0.00 0.00 0.00 0.00
48 0.00 0.00 1.35 0.47
49 0.00 0.00 0.45 0.16
50 0.00 0.00 0.45 0.16
51 0.00 0.00 1.84 0.64
52 0.00 0.00 1.39 0.48
53 0.00 0.00 1.89 0.65
54 0.00 0.00 1.55 0.54
55 0.00 0.00 1.66 0.58
56 0.00 0.00 3.03 1.00
57 0.00 0.00 1.86 0.64
58 0.00 0.00 2.58 0.89
59 0.00 0.00 1.01 0.35
60 0.00 0.00 0.81 0.28
61 0.00 0.00 1.60 0.52
144 0.00 0.00 0.00 0.00
62 -35.81 0.00 0.00 0.00
63 0.00 0.00 30.00 23.00
145 -26.48 0.00 0.00 0.00
64 50.00 0.00 0.00 0.00
65 8.00 0.00 0.00 0.00
66 0.00 0.00 1.02 0.35
67 0.00 0.00 1.02 0.35
68 0.00 0.00 3.80 1.25
69 0.00 0.00 1.19 0.41
0
148 1.0153 -10.00 10.00
73 1.0205 -20.00 20.00
232 1.0010 -20.00 20.00
159 0.9583 -25.00 25.00
161 0.9632 12.00 35.00
1 1.0250 -240.00 240.00
167 1.0520 -11.00 96.00
168 1.0520 -153.00 153.00
92 1.0000 -30.00 56.00
174 0.9900 -24.00 77.00
254 1.0435 -500.00 1500.00
178 1.0233 -60.00 120.00
179 1.0103 -25.00 200.00
100 1.0550 -125.00 350.00
186 1.0510 -50.00 75.00
281 1.0435 -100.00 300.00
258 1.0528 -15.00 35.00
102 1.0528 -50.00 100.00
103 1.0735 -25.00 50.00
106 1.0535 -50.00 175.00
189 1.0435 -50.00 90.00
108 0.9630 -10.00 15.00
116 0.9290 -40.00 90.00
2 0.9829 -50.00 150.00
199 1.0522 -45.00 90.00
117 1.0077 -15.00 35.00
3 1.0522 -50.00 80.00
122 1.0650 -100.00 400.00
123 1.0650 -100.00 400.00
203 1.0551 -300.00 300.00
124 1.0435 -1000.00 1000.00
289 1.0150 -260.00 260.00
4 1.0100 -150.00 150.00
265 1.0080 -60.00 60.00
133 1.0000 -320.00 320.00
5 1.0500 -300.00 300.00
6 1.0000 -300.00 300.00
7 1.0400 -250.00 250.00
8 1.0000 -500.00 500.00
9 1.0165 -300.00 300.00
221 1.0100 -200.00 200.00
10 1.0000 -400.00 400.00
11 1.0500 -600.00 600.00
222 0.9930 40.00 100.00
139 1.0100 40.00 80.00
28 1.0507 -210.00 210.00
29 1.0507 -280.00 280.00
30 1.0323 -420.00 420.00
31 1.0145 -100.00 100.00
32 1.0145 -224.00 224.00
33 1.0507 0.00 350.00
34 1.0507 0.00 120.00
35 1.0290 -224.00 224.00
36 1.0500 -200.00 200.00
37 1.0145 0.00 42.00
38 1.0507 -500.00 500.00
39 0.9967 0.00 25.00
40 1.0212 -90.00 90.00
41 1.0145 -150.00 150.00
42 1.0017 0.00 150.00
43 0.9893 0.00 87.00
44 1.0507 -100.00 600.00
45 1.0507 -125.00 325.00
46 1.0145 -200.00 300.00
225 0.9945 -2.00 2.00
62 1.0000 -17.35 17.35
145 1.0000 -12.80 12.83
64 1.0000 -38.00 38.00
65 1.0000 -6.00 6.00
0
5 0.2000 0.2000 0.1375 210.0000 600.0000
4 0.3000 0.1750 0.1750 200.0000 580.0000
9 0.3000 0.1250 0.1000 200.0000 800.0000
11 0.3500 0.2250 0.1283 300.0000 800.0000
28 0.2200 0.1430 0.2500 300.0000 600.0000
29 0.7500 0.1250 0.2500 300.0000 800.0000
30 0.5400 0.1950 0.2520 600.0000 1600.0000
31 0.3800 0.2200 0.4390 200.0000 400.0000
32 0.3600 0.1250 0.6350 200.0000 500.0000
33 0.9000 0.1300 0.0250 200.0000 500.0000
34 0.8300 0.2300 0.0730 100.0000 300.0000
35 0.4400 0.1430 0.3120 320.0000 600.0000
36 0.1200 0.1400 0.6650 400.0000 800.0000
38 0.5400 0.1150 0.1020 600.0000 1600.0000
40 0.6600 0.1550 0.2650 100.0000 350.0000
41 0.8200 0.1600 0.7000 250.0000 600.0000
42 0.4400 0.1450 0.1050 250.0000 600.0000
43 0.3500 0.1270 0.4500 80.0000 300.0000
44 0.5400 0.1250 0.1220 600.0000 1600.0000
45 0.3800 0.2000 0.1390 500.0000 900.0000
46 0.3600 0.1250 0.2350 400.0000 800.0000
0
0

23
IEEE4.dat Normal file
View File

@@ -0,0 +1,23 @@
4 4 1. 18 .1
1.e-5 2
1 4
0
1 1 2 .1 .4 0.01528
2 1 4 .12 .5 0.0192
3 2 4 .08 .4 0.01413
0
0
1 1 3 .0 .3 0.90909 .9 1.15
0
1 0. 0. 0.3 0.18
2 0. 0. 0.55 0.13
3 0.5 0. 0. 0.
4 0.36 0.26 0. 0.
0
3 1.1 -0.1 0.6
4 1.05 -0.6 0.6
0
3 44.4 351. 50. 0.3 1.2
4 40.6 389. 50. 0.3 1.2
0
0

25
Jacobi.m Normal file
View File

@@ -0,0 +1,25 @@
function [jaco]=Jacobi(Balance,busNum,QGi,Volt,YdotSinVolt,YdotCosVolt,diag_Volt_YdotSin,diag_Volt_YdotCos)
diag_YdotSinVolt_=diag(YdotSinVolt);
diag_YdotCosVolt_=diag(YdotCosVolt);
dPdTyta=diag_Volt_YdotSin*diag(Volt)-diag_YdotSinVolt_*diag(Volt); % 简化第三次
dQdTyta=-diag_Volt_YdotCos*diag(Volt)+diag_YdotCosVolt_*diag(Volt);%dQ/dThyta
dPdV=diag_YdotCosVolt_+diag_Volt_YdotCos;%dP/dV
dQdV=diag_YdotSinVolt_+diag_Volt_YdotSin;%dQ/dV
%% 平衡节点相角不变
dPdTyta(Balance,:)=0;
dPdTyta(:,Balance)=0;
dPdTyta=dPdTyta+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum);
dQdTyta(:,Balance)=0;
dPdV(Balance,:)=0;
%% PV节点电压不变
dQdV(QGi,:)=0;
dQdV(:,QGi)=0;
dQdV=dQdV+sparse(QGi,QGi,ones(length(QGi),1),busNum,busNum);
dQdTyta(QGi,:)=0;
dPdV(:,QGi)=0;
% jaco=[
% dPdV,dPdTyta;
% dQdV,dQdTyta;
% ];
jaco=cat(1,cat(2,dPdV,dPdTyta),cat(2,dQdV,dQdTyta));
end

18
Modification.m Normal file
View File

@@ -0,0 +1,18 @@
function [Init_Y,PG,QG,Volt,UAngel,PD,QD]=Modification(Init_Y,deltX,deltY,PG,QG,Volt,UAngel,ContrlCount,Balance,Busnum,PD,QD,Loadi)
Init_Y=Init_Y+deltY';
t=deltX(1:size(Loadi,1)*2);
PD=PD+t(1:length(Loadi));
QD=QD+t(length(Loadi)+1:length(Loadi)*2);
t=deltX(size(Loadi,1)*2+1:ContrlCount)';
t(Busnum*3+(Balance-1)*3+1)=0;
t(Busnum*3+(Balance-1)*3+2)=0;
t(Busnum*3+(Balance-1)*3+3)=0;
%Volt=Volt+AlphaP*t(2:2:2*Busnum);Ôİʱ¸ÄÒ»ÏÂ20111227
%UAngel=UAngel+AlphaP*t(1:2:2*Busnum);Ôİʱ¸ÄÒ»ÏÂ20111227
% balVolt=Volt( (Balance-1)*3+1 );
Volt=Volt+t(1:Busnum*3)';
Volt( (Balance-1)*3+1 )=1;
Volt( (Balance-1)*3+2 )=1;
Volt( (Balance-1)*3+3 )=1;
UAngel=UAngel+t(Busnum*3+1:2*Busnum*3)';
end

5
Modify.m Normal file
View File

@@ -0,0 +1,5 @@
function [Volt, Vangle]=Modify(Volt,Vangle,dV,dVangle,optimalu)
%Volt=Volt+Volt.*dV;
Volt=Volt+optimalu*dV;
Vangle=Vangle+optimalu*dVangle;
end

330
OPF.m Normal file
View File

@@ -0,0 +1,330 @@
function [JMeasurement,AME_Volt,AME_VAngle,AME_PD,AME_QD,AME_mVolt,AME_mPD,AME_mQD,isConverged,totalTime]=OPF()
tic
clc
clear
lineZ=readLineZ('feeder104\lineParameter.txt');
[ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ...
phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY ...
cap]=dataRead(lineZ,'feeder104\data1.txt');
% phaseASpotLoadP(phaseASpotLoadP==0)=0.002;
% phaseBSpotLoadP(phaseBSpotLoadP==0)=0.002;
% phaseCSpotLoadP(phaseCSpotLoadP==0)=0.002;
% phaseASpotLoadQ(phaseASpotLoadQ==0)=0.002;
% phaseBSpotLoadQ(phaseBSpotLoadQ==0)=0.002;
% phaseCSpotLoadQ(phaseCSpotLoadQ==0)=0.002;
%% 潮流计算begin
a=exp(1j*2*pi/3);
Tp2f=1/3*[1 1 1;
1 a a^2;
1 a^2 a];
Tp2f=sparse(Tp2f);
Tf2p=inv(Tp2f);
fsY1amp=abs(fsY1);
[r,c,fsY1ang]=find(fsY1);
fsY1ang=angle(fsY1ang);
Pabc=phaseASpotLoadP+phaseBSpotLoadP+phaseCSpotLoadP;
Qabc=phaseASpotLoadQ+phaseBSpotLoadQ+phaseCSpotLoadQ;
busNum=length(phaseASpotLoadP);
%给序电压赋初值
Vmf1=sparse(ones(busNum,1));
Vaf1=sparse(zeros(busNum,1));
%先求解正序的
PQi=nodeNum;
PG=sparse(busNum,1);
QG=sparse(busNum,1);
QGi=[Balance];
PD=Pabc/3;
QD=Qabc/3;
Loadi=find(PD~=0);
maxD=100000;% 最大不平衡量
EPS=1e-6;
k=0;
kmax=20;
fsY11=fsY1;
fsY00=fsY0;
fsY22=fsY2;
Vf2=sparse(busNum,1);
If2=sparse(busNum,1);
Vf0=sparse(busNum,1);
If0=sparse(busNum,1);
%准备序矩阵
%平衡节点置0置1
fsY2(Balance,:)=0;
fsY2(:,Balance)=0;
fsY2=fsY2+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum);
%平衡节点置0置1
fsY0(Balance,:)=0;
fsY0(:,Balance)=0;
fsY0=fsY0+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum);
%%LU分解
[fsY0L,fsY0U,fsY0P,fsY0Q,fsY0R]=lu(fsY0);
[fsY2L,fsY2U,fsY2P,fsY2Q,fsY2R]=lu(fsY2);
%算初始补偿功率
tic
VoltpA=sparse(ones(busNum,1));
VoltpB=sparse(ones(busNum,1)).*exp(1j*-120/180*pi);
VoltpC=sparse(ones(busNum,1)).*exp(1j*+120/180*pi);
while(k<=kmax+10 && maxD> EPS)
k=k+1;
%把补偿电容看作负荷
SA=VoltpA.*conj(VoltpA.*sparse(cap.capNode,1,1j*cap.capB(:,1),busNum,1));
SB=VoltpB.*conj(VoltpB.*sparse(cap.capNode,1,1j*cap.capB(:,2),busNum,1));
SC=VoltpC.*conj(VoltpC.*sparse(cap.capNode,1,1j*cap.capB(:,3),busNum,1));
%先不要电容，已经加到序的导纳矩阵中了。20150405 By 杜孟远
SA=0;
SB=0;
SC=0;
iterPD=PD+real(SA+SB+SC)/3;
iterQD=QD+imag(SA+SB+SC)/3;
iterPhaseASpotLoadP=phaseASpotLoadP+real(SA);
iterPhaseBSpotLoadP=phaseBSpotLoadP+real(SB);
iterPhaseCSpotLoadP=phaseCSpotLoadP+real(SC);
iterPhaseASpotLoadQ=phaseASpotLoadQ+imag(SA);
iterPhaseBSpotLoadQ=phaseBSpotLoadQ+imag(SB);
iterPhaseCSpotLoadQ=phaseCSpotLoadQ+imag(SC);
[dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance,busNum, ...
PQi,PG,QG,QGi,iterPD,iterQD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量
maxD=max(abs([dP;dQ;]));
jaco=Jacobi(Balance,busNum,QGi,Vmf1,YdotSinVolt,YdotCosVolt,diag_Volt_YdotSin,diag_Volt_YdotCos);%雅克比矩阵
[dV, dVangle]=Solv(busNum,jaco,dP,dQ);%解出修正量
[Vmf1, Vaf1]=Modify(Vmf1,Vaf1,dV,dVangle,1);
fprintf('第 %d 次迭代， 最大不平衡量为 %f\n',k,full(maxD));
%转换为三相电压
VoltpABC=Tp2f\conj([ Vf0'; (Vmf1.*exp(1j*Vaf1))'; Vf2']);%用Tp2f\ 代替Tf2p*
VoltpA=conj(VoltpABC(1,:)');
CurpA=-conj((iterPhaseASpotLoadP+1j*iterPhaseASpotLoadQ)./VoltpA);
VoltpB=conj(VoltpABC(2,:)');
CurpB=-conj((iterPhaseBSpotLoadP+1j*iterPhaseBSpotLoadQ)./VoltpB);
VoltpC=conj(VoltpABC(3,:)');
CurpC=-conj((iterPhaseCSpotLoadP+1j*iterPhaseCSpotLoadQ)./VoltpC);
f012=Tp2f*conj([CurpA';CurpB';CurpC']);
If0=conj(f012(1,:)');
If1=conj(f012(2,:)');
If2=conj(f012(3,:)');
If0(Balance)=0;
If2(Balance)=0;
%Vf0=fsY0\If0;
Vf0=fsY0Q*(fsY0U\(fsY0L\(fsY0P*(fsY0R\If0))));
%Vf2=fsY2\If2;
Vf2=fsY2Q*(fsY2U\(fsY2L\(fsY2P*(fsY2R\If2))));
fprintf('迭代时间%f\n',toc);
%
end
FortiscueToc=toc;
fprintf('Fortiscue法计算时间 %f\n',FortiscueToc);
Vf1=Vmf1.*exp(1j*Vaf1);
(Vf0.*conj(fsY00*Vf0)+Vf1.*conj(fsY11*Vf1)+Vf2.*conj(fsY22*Vf2))*3;%包含补偿电容的功率
conj(Tf2p*[If0(2);If1(2);If2(2)]).*(Tf2p*[Vf0(2);Vf1(2);Vf2(2)]);
IpABC=Tf2p*conj([If0';If1';If2']);
%转换回三相电压
VoltpABC=Tf2p*conj([ Vf0'; Vf1'; Vf2']);
disp([' A B C'])
full(abs(VoltpABC'))
fprintf('节点号对应\n');
disp([setIJ,nodeNum ])
%%检查反推回去的功率是否满足
ub=checkSSatisfied(Balance,phaseABCY,VoltpABC, ...
phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP, ...
phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ );
PGQG=CalPGQG(Balance,phaseABCY,VoltpABC,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP,phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ );
fprintf('最大不平衡量为%f\n\n',full(max(abs(ub))))
%% 潮流计算end
busNum=length(nodeNum);
Busnum=busNum;
% PQi=setxor(nodeNum,Balance);
QGi=[Balance];
% fprintf('开始牛顿法迭代\n');
[r,c,GB]=find(phaseABCY);
Y=abs(phaseABCY);
Yangle=angle(GB);
Vp3=sparse(ones(busNum*3,1));%给电压赋初值
Vp3(2:3:end)=Vp3(2:3:end)*exp(1j*-120/180*pi);
Vp3(3:3:end)=Vp3(3:3:end)*exp(1j*+120/180*pi);
PQi3P=zeros(length(Loadi)*3,1);
PQi3P(1:3:end)=(Loadi-1)*3+1;
PQi3P(2:3:end)=(Loadi-1)*3+2;
PQi3P(3:3:end)=(Loadi-1)*3+3;
Loadi=PQi3P;
PD3P=sparse(Busnum*3,1);
QD3P=sparse(Busnum*3,1);
PD3P(1:3:end)=phaseASpotLoadP*0.9;
PD3P(2:3:end)=phaseBSpotLoadP*0.9;
PD3P(3:3:end)=phaseCSpotLoadP*0.9;
QD3P(1:3:end)=phaseASpotLoadQ*0.9;
QD3P(2:3:end)=phaseBSpotLoadQ*0.9;
QD3P(3:3:end)=phaseCSpotLoadQ*0.9;
PD3P=PD3P(Loadi);
QD3P=QD3P(Loadi);
QGi3P=zeros(length(QGi)*3,1);
QGi3P(1:3:end)=(QGi-1)*3+1;
QGi3P(2:3:end)=(QGi-1)*3+2;
QGi3P(3:3:end)=(QGi-1)*3+3;
Vp3m=abs(Vp3);
Vp3a=angle(Vp3);
Balance3P=zeros(length(Balance)*3,1);
Balance3P(1:3:end)=(Balance-1)*3+1;
Balance3P(2:3:end)=(Balance-1)*3+2;
Balance3P(3:3:end)=(Balance-1)*3+3;
% PGA=sum(PD3P(1:3:end))*1.03;
% PGB=sum(PD3P(2:3:end))*1.03;
% PGC=sum(PD3P(3:3:end))*1.03;
% QGA=sum(QD3P(1:3:end))*1.03;
% QGB=sum(QD3P(2:3:end))*1.03;
% QGC=sum(QD3P(3:3:end))*1.03;
% PG3P=sparse( (Balance-1)*3+1:(Balance-1)*3+3,1,[PGA,PGB,PGC],Busnum*3,1);
% QG3P=sparse( (Balance-1)*3+1:(Balance-1)*3+3,1,[QGA,QGB,QGC],Busnum*3,1);
PG3P=real(PGQG);
QG3P=imag(PGQG);
Vp3a((Balance-1)*3+1)=0;
Vp3a((Balance-1)*3+2)=-120/180*pi;
Vp3a((Balance-1)*3+3)=+120/180*pi;
%准备量测量和方差
%真实值
%三相电压幅值
rVoltABCV=zeros(busNum*3,1);
rVoltABCV(1:3:end)=abs(VoltpABC(1,:));
rVoltABCV(2:3:end)=abs(VoltpABC(2,:));
rVoltABCV(3:3:end)=abs(VoltpABC(3,:));
%三相电压相角
rVoltABCA=zeros(busNum*3,1);
rVoltABCA(1:3:end)=angle(VoltpABC(1,:));
rVoltABCA(2:3:end)=angle(VoltpABC(2,:));
rVoltABCA(3:3:end)=angle(VoltpABC(3,:));
%三相负荷
rPD3P=zeros(busNum*3,1);
rPD3P(1:3:end)=phaseASpotLoadP;
rPD3P(2:3:end)=phaseBSpotLoadP;
rPD3P(3:3:end)=phaseCSpotLoadP;
rQD3P=zeros(busNum*3,1);
rQD3P(1:3:end)=phaseASpotLoadQ;
rQD3P(2:3:end)=phaseBSpotLoadQ;
rQD3P(3:3:end)=phaseCSpotLoadQ;
rPD3P=rPD3P(Loadi);
rQD3P=rQD3P(Loadi);
%量测量
sigma=0.03;
VoltSigma=(1+normrnd(0,sigma,length(rVoltABCV),1));
mVoltABCV=rVoltABCV.*VoltSigma;
PD3PSigma=(1+normrnd(0,sigma,length(rPD3P),1));
mPD3P=rPD3P.*PD3PSigma;
QD3PSigma=(1+normrnd(0,sigma,length(rQD3P),1));
mQD3P=rQD3P.*QD3PSigma;
%量测方差
wVolt=1./(abs(mVoltABCV*sigma).^2);
wPD=1./(abs(mPD3P*sigma).^2);
wQD=1./(abs(mQD3P*sigma).^2);
%% 没有量测量
% noLoadi=[1,5,6,10,11];
noLoadi=[1,8,10,12];
% noLoadi=Loadi;
% noLoadi=[1,11];
noLoadi=[];
noPQi3P=zeros(length(noLoadi)*3,1);
noPQi3P(1:3:end)=(noLoadi-1)*3+1;
noPQi3P(2:3:end)=(noLoadi-1)*3+2;
noPQi3P(3:3:end)=(noLoadi-1)*3+3;
% noPQi3P=Loadi;
wVolt(setdiff(1:length(wVolt),Loadi))=0;%只有负荷处才有电压量测。
wVolt( noPQi3P)=0;
wPD(ismember( Loadi,noPQi3P))=0;
wQD(ismember( Loadi,noPQi3P))=0;
%%
% RestraintCount=size(Loadi,1)*2+length(rVoltABCV); %约束条件数
% RestraintCount=size(Loadi,1)*2; %约束条件数
Init_Y=sparse(1,2*Busnum*3);%与学姐一致
KK=0;
ContrlCount=size(Loadi,1)*2+Busnum*6;
kmax=200;
%%
Precision=1e-5;
CenterA=0.1;
%% 加误差
Volt=1*Vp3m;
% Volt=mVoltABCV;
UAngel=Vp3a;
maxD=100;
tic
while(maxD>Precision)
if KK>kmax
break;
end
%% 开始计算OPF
%% 形成等式约束的雅克比
deltH=func_deltH(Busnum,Volt,Y,UAngel,r,c,Yangle,Loadi);
%% 形成不等式约束的雅克比
% deltG=func_deltG(Busnum,Loadi);
%%
%% 形成海森阵
deltdeltF=func_deltdeltF(wVolt,wPD,wQD,ContrlCount,Loadi);
%% 形成ddHy
ddh=func_ddh(Volt,Init_Y,Busnum,Y,UAngel,r,c,Yangle,Loadi,ContrlCount);
%% 开始构建ddg
%% 开始构建deltF
deltF=func_deltF(wVolt,wPD,wQD,mPD3P,PD3P,QD3P,mQD3P,Volt,mVoltABCV,Busnum,Loadi);
%% 形成方程矩阵
% Luu=Init_U'.*Init_W'+Init_u*ones(RestraintCount,1);
% Lul=Init_L'.*Init_Z'-Init_u*ones(RestraintCount,1);
% Mat_G=FormG(Volt,PD3P,QD3P,Loadi);
Mat_H=FormH(Busnum,Volt,PG3P,PD3P,QG3P,QD3P,Y,UAngel,r,c,Yangle,Loadi);
Ly=Mat_H;
% if isSetBound==0
% 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));
% 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));
% isSetBound=1;
% end
% Lz=FormLz(Mat_G,Init_L,Loadi,rPD3P,rQD3P,rVoltABCV,PDL,QDL);
% Lw=FormLw(Mat_G,Init_U,Loadi,rPD3P,rQD3P,rVoltABCV,PDU,QDU);
Lx=FormLx(deltF,deltH,Init_Y);
% YY=FormYY(Lul,Lz,Ly,Luu,Lw,Lx);
%% 开始解方程
% fprintf('迭代次数 %d Gap %f\n',KK+1,plotGap(KK+1));
XX=SolveIt(deltdeltF,ddh,deltH,Init_Y,Ly,ContrlCount,Lx,Balance,Busnum,Loadi);
%%取各分量
[deltX,deltY]=AssignXX(XX,ContrlCount,Loadi,Balance,Busnum);
[Init_Y,PG,QG,Volt,UAngel,PD3P,QD3P]=Modification(Init_Y,deltX,deltY,PG,QG,Volt,UAngel,ContrlCount,Balance,Busnum,PD3P,QD3P,Loadi);
maxD=max(abs([deltX;deltY]));
fprintf('%f\n',full(maxD));
KK=KK+1;
end
totalTime=toc
(rVoltABCV-Volt)./rVoltABCV*100;
%% 计算统计量
%目标函数均值
JMeasurement=sum(((mVoltABCV-Volt)./mVoltABCV./sigma).^2)+sum(((mPD3P-PD3P)./mPD3P./sigma).^2)+sum(((mQD3P-QD3P)./mQD3P./sigma).^2);
%估计误差统计
%量测量数量
Busnum=busNum;
mCount=Busnum*3+length(Loadi)*3*2;
%估计量质量
AME_Volt=sum(sum(abs( abs(rVoltABCV)-abs(Volt))));
AME_VAngle=sum(sum(abs( rVoltABCA-UAngel)));
AME_PD=sum(sum(abs(rPD3P-PD3P)));
AME_QD=sum(sum(abs(rQD3P-QD3P)));
%计算与量测值的
AME_mVolt=sum(sum(abs( mVoltABCV-rVoltABCV)));
AME_mPD=sum(sum(abs(rPD3P-mPD3P)));
AME_mQD=sum(sum(abs(rQD3P-mQD3P)));
%返回收敛信息
isConverged=1;
if KK>=kmax
isConverged=0;
end
if abs(maxD)>Precision
isConverged=0;
end
fprintf('迭代次数%d\n',KK);
toc
end

30
OPF_Init.m Normal file
View File

@@ -0,0 +1,30 @@
function [Volt,UAngel,Init_Z,Init_W,Init_L,Init_U,Init_Y,PG,QG,RestraintCount,wPG,wQG,wPD,wQD,PD,PD0,QD,randPDind,Loadi]=OPF_Init(Busnum,Balance,PG,QG,Volt,GenU,GenL,PVi,PGi,PVQU,PVQL,PD0,QD0,QD,PD)
Loadi=find(QD~=0 | PD~=0);
%Loadi=[1:Busnum]';
t_Bal_volt=Volt(Balance);
Volt=sparse(1*ones(1,Busnum));
Volt(Balance)=t_Bal_volt;
UAngel=sparse(1,Busnum);
tPU=sparse(GenU(:,2));% 发电机有功上限
tQU=sparse(PVQU(:,1));% 无功上限
tPL=sparse(GenL(:,2));% 发电机有功下限
tQL=sparse(PVQL(:,1));% 无功下限
wPG=0;
wQG=0;
%randInt=randperm(size(Loadi,1));
%randPDind=randInt(1:10);
randPDind=0;
wPD=zeros(size(Loadi,1),1);
wPD(1:2:end)=0;
wQD=zeros(size(Loadi,1),1);
wQD(1:2:end)=0;
%wD(randPDind)=0;%一些负荷不约束
%wD(7)=0;
% wD(11)=0;
PD=0.8*PD0;
%powerFacter=0.98;
%QD=PD*sqrt((1-powerFacter^2)/powerFacter^2);
QD=.8*QD0;
end

1
README.md Normal file
View File

@@ -0,0 +1 @@
把原来内点法的程序给改了，改成只有等式约束的普通WLS方法。

16
Solv.m Normal file
View File

@@ -0,0 +1,16 @@
function [dV, dVangle]=Solv(busNum,jaco,dP,dQ)
%y=klu (jaco, '\', full(-[dP;dQ]));
y=jaco\(-[dP;dQ]);
% [L,U] = luinc(jaco,1e-3); %luinc(A,'0')也可以试一下，是一种完全不同的ILU
% tol=1e-5; %残量的精度要求
% restart=30; % 30-50之间吧，不要过小
% maxit=100; %看情况，如果不收敛就适当调大
% [x,flag]=gmres(jaco,-[dP;dQ],restart,tol,maxit);
% y=sparse(x);
y=sparse(y);
dV=y(1:busNum);
dVangle=y(busNum+1:end);
end

46
SolveIt.m Normal file
View File

@@ -0,0 +1,46 @@
function XX=SolveIt(deltdeltF,ddh,deltH,Init_Y,Ly,ContrlCount,Lx,Balance,Busnum,Loadi)
LxComa=FormLxComa(Lx);%Lx=deltF-deltH*Init_Y';
Lx=LxComa;
H=-deltdeltF+ddh;
aa=[
H,deltH;
deltH',zeros(size(Init_Y,2));
];
yy=[Lx;-Ly];
%% ƽºâ½Úµãµçѹ²»±ä
t=size(Loadi,1)*2;
aa(t+(Balance-1)*3+1,:)=0;
aa(t+(Balance-1)*3+2,:)=0;
aa(t+(Balance-1)*3+3,:)=0;
aa(:,t+(Balance-1)*3+1)=0;
aa(:,t+(Balance-1)*3+2)=0;
aa(:,t+(Balance-1)*3+3)=0;
%aa(t+Balance,t+Balance)=1;
aa=aa+sparse(t+(Balance-1)*3+1,t+(Balance-1)*3+1,ones(length(Balance),1),ContrlCount+2*Busnum*3,ContrlCount+2*Busnum*3);
aa=aa+sparse(t+(Balance-1)*3+2,t+(Balance-1)*3+2,ones(length(Balance),1),ContrlCount+2*Busnum*3,ContrlCount+2*Busnum*3);
aa=aa+sparse(t+(Balance-1)*3+3,t+(Balance-1)*3+3,ones(length(Balance),1),ContrlCount+2*Busnum*3,ContrlCount+2*Busnum*3);
%%
t=size(Loadi,1)*2+Busnum*3;
aa(t+(Balance-1)*3+1,:)=0;
aa(t+(Balance-1)*3+2,:)=0;
aa(t+(Balance-1)*3+3,:)=0;
aa(:,t+(Balance-1)*3+1)=0;
aa(:,t+(Balance-1)*3+2)=0;
aa(:,t+(Balance-1)*3+3)=0;
%aa(t+Balance,t+Balance)=1;
aa=aa+sparse(t+(Balance-1)*3+1,t+(Balance-1)*3+1,ones(length(Balance),1),ContrlCount+2*Busnum*3,ContrlCount+2*Busnum*3);
aa=aa+sparse(t+(Balance-1)*3+2,t+(Balance-1)*3+2,ones(length(Balance),1),ContrlCount+2*Busnum*3,ContrlCount+2*Busnum*3);
aa=aa+sparse(t+(Balance-1)*3+3,t+(Balance-1)*3+3,ones(length(Balance),1),ContrlCount+2*Busnum*3,ContrlCount+2*Busnum*3);
%%
dxdy=aa\yy;
%% KLU
%spy(aa)
%dxdy = klu(aa,'\',full(yy));
%%
dX=dxdy(1:ContrlCount);
dY=dxdy(ContrlCount+1:end);
XX=[
dX;
dY;
];
end

18
Unbalance.m Normal file
View File

@@ -0,0 +1,18 @@
function [dP, dQ, YdotSinVmf1, YdotCosVmf1, diag_Vmf1_YdotSin, diag_Vmf1_YdotCos]=Unbalance(Balance,busNum, ...
PQi,PG,QG,QGi,PD,QD,Vmf1,Vaf1,Y,Yangle,r,c,Vf2,If2,Vf0,If0)
t1=sparse(r,c,Vaf1(r)-Vaf1(c) -Yangle,busNum,busNum) ;
YdotSin=Y.* ( spfun(@sin,t1) );
YdotCos=Y.* ( spfun (@cos, t1 ) );
orderedPDPG=sparse(PQi,1,PD-PG,busNum,1); %节点功率按顺序排列
orderedQDQG=sparse(PQi,1,QD-QG,busNum,1); %节点功率按顺序排列
diag_Vmf1_YdotCos=diag(Vmf1)*YdotCos;
diag_Vmf1_YdotSin=diag(Vmf1)*YdotSin;
YdotCosVmf1=YdotCos*Vmf1;
YdotSinVmf1=YdotSin*Vmf1;
diag_Vmf1_YdotCosVmf1=diag_Vmf1_YdotCos*Vmf1;
diag_Vmf1_YdotSinVmf1=diag_Vmf1_YdotSin*Vmf1;
dP=diag_Vmf1_YdotCosVmf1+orderedPDPG+real(Vf2.*conj(If2)+Vf0.*conj(If0));
dQ=diag_Vmf1_YdotSinVmf1+orderedQDQG+imag(Vf2.*conj(If2)+Vf0.*conj(If0));
dP(Balance)=0;
dQ(QGi)=0;
end

33
admmatrix.m Normal file
View File

@@ -0,0 +1,33 @@
function [GB,Y,r,c,Angle] = admmatrix(Busnum,Linei,Linej,Liner,Linex,Lineb,Transfori...
,Transforj,Transforr,Transforx,Transfork0,Branchi,Branchg,Branchb)
%**************************************************************************
% 程序功能 : 子函数——形成节点导纳矩阵Y
% 编 者：
% 编制时间：2010.12
%**************************************************************************
%% 支路导纳计算
G = -sparse(Linei,Linej,Liner./(Liner.^2+Linex.^2),Busnum,Busnum) - sparse(Linej,Linei,Liner./(Liner.^2+Linex.^2),Busnum,Busnum);
G = G - sparse(1:Busnum,1:Busnum,sum(G,2)'); % 计算各线路支路电导
B = sparse(Linei,Linej,Linex./(Liner.^2+Linex.^2),Busnum,Busnum)+sparse(Linej,Linei,Linex./(Liner.^2+Linex.^2),Busnum,Busnum);
B = B - sparse(1:Busnum,1:Busnum,sum(B,2)')+sparse(Linei,Linei,Lineb,Busnum,Busnum)+sparse(Linej,Linej,Lineb,Busnum,Busnum);
%% 变压器支路计算
if Transfori>0
mr = Transforr./(Transforr.^2+Transforx.^2); % 计算变压器支路电导
mx = -Transforx./(Transforr.^2+Transforx.^2); % 计算变压器支路电纳
G = G-sparse(Transfori,Transforj,mr./Transfork0,Busnum,Busnum)-sparse(Transforj,Transfori,mr./Transfork0,Busnum,Busnum)...
+sparse(Transfori,Transfori,mr./Transfork0./Transfork0,Busnum,Busnum)+sparse(Transforj,Transforj,mr,Busnum,Busnum);
B = B-sparse(Transfori,Transforj,mx./Transfork0,Busnum,Busnum)-sparse(Transforj,Transfori,mx./Transfork0,Busnum,Busnum)...
+sparse(Transfori,Transfori,mx./Transfork0./Transfork0,Busnum,Busnum)+sparse(Transforj,Transforj,mx,Busnum,Busnum);
end
%% 接地支路计算
if Branchi>0 % 判断有无接地支路
B = B+sparse(Branchi,Branchi,Branchb,Busnum,Busnum);
G = G+sparse(Branchi,Branchi,Branchg,Busnum,Busnum);
end
%% 化作极坐标形式
GB = G+B.*1i; %将电导，电纳合并，写成复数形式
Y = abs(GB); %求节点导纳幅值
%Y=full(Y);
[r,c] = find(Y);
Angle = angle(GB(GB~=0)); %求节点导纳角度
%Angle=angle(GB);

25
checkSSatisfied.m Normal file
View File

@@ -0,0 +1,25 @@
function [ output_args ] = checkSSatisfied(Balance,phaseABCY,VoltpABC,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP,phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ )
%%利用Fortiscue方法计算得到了相电压和相角后，反推负荷功率，检查是否和给定的功率一致。
busNum=size(phaseABCY,1)/3;
VoltpA=VoltpABC(1,:);
VoltpB=VoltpABC(2,:);
VoltpC=VoltpABC(3,:);
Vp3=sparse(busNum*3,1);
Vp3(1:3:end)=VoltpA;
Vp3(2:3:end)=VoltpB;
Vp3(3:3:end)=VoltpC;
Ip3=phaseABCY*Vp3;
Sp3=Vp3.*conj(Ip3);
%VoltpABC.*conj(IpABC);
pLoadABC=sparse(length(phaseASpotLoadP)*3,1);
pLoadABC(1:3:end)=phaseASpotLoadP+1j*phaseASpotLoadQ;
pLoadABC(2:3:end)=phaseBSpotLoadP+1j*phaseBSpotLoadQ;
pLoadABC(3:3:end)=phaseCSpotLoadP+1j*phaseCSpotLoadQ;
ck=Sp3+pLoadABC;
ck(3*(Balance-1)+1:3*(Balance-1)+3)=0;
if any(real(ck)>1e-5) ||any(imag(ck)>1e-5)
fprintf('反推回的功率不匹配。潮流方程约束不满足。\n');
end
output_args=ck;
end

80
dataRead.m Normal file
View File

@@ -0,0 +1,80 @@
function [ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ...
phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY ...
cap] = dataRead(lineZ,dataFile )
data=dlmread(dataFile);
baseHighU=data(1,2);
fsY0=sparse(0);
fsY1=sparse(0);
fsY2=sparse(0);
zeroEntries=find(data(:,1)==0);
lines=data(zeroEntries(1)+1:zeroEntries(2)-1,:);
[setIJ,nodeNum]=numberNode(lines);
if size(setIJ,2)>size(setIJ,1)
setIJ=setIJ';
end
Balance=data(1,1);
Balance=nodeNum(setIJ==Balance);
%%three-phase
phaseABCY=sparse(3*length(nodeNum),3*length(nodeNum));
%把线路型号都给读出来
linePar=lines(:,1);
linePar=unique(linePar);
for I=1:length(linePar)
[fs30,fs31,fs32,retphaseABCY]=lineWithConfig(setIJ,nodeNum,lineZ,lines,linePar(I),baseHighU);
%phaseABCY 三相的导纳矩阵3n x 3n 维
fsY0=fsY0+fs30;
fsY1=fsY1+fs31;
fsY2=fsY2+fs32;
phaseABCY=phaseABCY+retphaseABCY;
end
%% 把电容加到矩阵中
cap=data(zeroEntries(2)+1:zeroEntries(3)-1,:);
if ~isempty(cap)
capNode=nodeNum(ismember(setIJ,cap(:,1)));
end
%电容引起的电流
capG=1j*cap(:,2)/1000;
% fs3Y=3/3*diag([Zl+2*Zm,Zl-Zm,Zl-Zm]);%这里不是1/3，IEEE上的公式错了。
fsY0=fsY0+sparse(capNode,capNode,capG,length(fsY0),length(fsY0));
fsY1=fsY1+sparse(capNode,capNode,capG,length(fsY1),length(fsY1));
fsY2=fsY2+sparse(capNode,capNode,capG,length(fsY2),length(fsY2));
%% spot load
busNum=length(nodeNum);
spotloads=data(zeroEntries(3)+1:zeroEntries(4)-1,:);
spotloads(:,2:end)=spotloads(:,2:end)/1000;
phaseASpotLoadP=sparse(busNum,1);
phaseBSpotLoadP=sparse(busNum,1);
phaseCSpotLoadP=sparse(busNum,1);
phaseASpotLoadQ=sparse(length(nodeNum),1);
phaseBSpotLoadQ=sparse(length(nodeNum),1);
phaseCSpotLoadQ=sparse(length(nodeNum),1);
phaseASpotLoadP( ismember(setIJ,spotloads(:,1) ) )=spotloads(:,2);
phaseBSpotLoadP( ismember(setIJ,spotloads(:,1) ) )=spotloads(:,4);
phaseCSpotLoadP( ismember(setIJ,spotloads(:,1) ) )=spotloads(:,6);
phaseASpotLoadQ( ismember(setIJ,spotloads(:,1) ) )=spotloads(:,3);
phaseBSpotLoadQ( ismember(setIJ,spotloads(:,1) ) )=spotloads(:,5);
phaseCSpotLoadQ( ismember(setIJ,spotloads(:,1) ) )=spotloads(:,7);
%% 补偿电容
cap=data(zeroEntries(2)+1:zeroEntries(3)-1,:);
if ~isempty(cap)
capNode=nodeNum(ismember(setIJ,cap(:,1)));
%数据文件中的补偿电容给的是补偿容量，要将其转换为导纳。
capB=cap(:,2:4)/1000;
cap=struct();
cap.capNode=capNode;
cap.capB=capB;
%把补偿电容加到矩阵中。
offSet=(capNode-1)*3+1;
phaseABCY=phaseABCY+sparse(offSet,offSet,1j*capB(:,1),busNum*3,busNum*3);
offSet=(capNode-1)*3+2;
phaseABCY=phaseABCY+sparse(offSet,offSet,1j*capB(:,2),busNum*3,busNum*3);
offSet=(capNode-1)*3+3;
phaseABCY=phaseABCY+sparse(offSet,offSet,1j*capB(:,3),busNum*3,busNum*3);
else
cap=struct();
cap.capNode=[];
cap.capB=sparse(1,3);
end
end

BIN
feeder104/Cap Data.xls Normal file
View File

Binary file not shown.

BIN
feeder104/Distributed Load Data.xls Normal file
View File

Binary file not shown.

BIN
feeder104/IEEE 13 Node Test Feeder.doc Normal file
View File

Binary file not shown.

BIN
feeder104/Line Data.xls Normal file
View File

Binary file not shown.

BIN
feeder104/Regulator Data.xls Normal file
View File

Binary file not shown.

BIN
feeder104/Spot Load Data.xls Normal file
View File

Binary file not shown.

BIN
feeder104/Transformer data.xls Normal file
View File

Binary file not shown.

BIN
feeder104/UG config.xls Normal file
View File

Binary file not shown.

BIN
feeder104/config.xls Normal file
View File

Binary file not shown.

13
feeder104/data - 副本 (2) - 副本.txt Normal file
View File

@@ -0,0 +1,13 @@
650
0
602 632 633 500
601 650 632 2000
601 632 671 2000
601 671 680 1000
0
0
633 160 120 120 110 90 90
671 385 385 385 220 220 220
0

16
feeder104/data - 副本.txt Normal file
View File

@@ -0,0 +1,16 @@
650
0
600 632 645 500
602 632 633 500
601 633 634 0
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 0
602 684 611 300
601 692 675 500
0

BIN
feeder104/data.xls Normal file
View File

Binary file not shown.

190
feeder104/data1.txt Normal file
View File

@@ -0,0 +1,190 @@
650 4.16
0
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
600 2632 2645 50
602 2632 2633 50
601 2633 2634 60
600 2645 2646 30
601 650 2632 200
602 2684 2652 80
601 2632 2671 200
600 2671 2684 30
600 2671 2692 30
602 2684 2611 30
601 2692 2675 50
600 3632 3645 50
602 3632 3633 50
601 3633 3634 60
600 3645 3646 30
601 650 3632 200
602 3684 3652 80
601 3632 3671 200
600 3671 3684 30
600 3671 3692 30
602 3684 3611 30
601 3692 3675 50
600 4632 4645 50
602 4632 4633 50
601 4633 4634 60
600 4645 4646 30
601 650 4632 200
602 4684 4652 80
601 4632 4671 200
600 4671 4684 30
600 4671 4692 30
602 4684 4611 30
601 4692 4675 50
600 5632 5645 50
602 5632 5633 50
601 5633 5634 60
600 5645 5646 30
601 650 5632 200
602 5684 5652 80
601 5632 5671 200
600 5671 5684 30
600 5671 5692 30
602 5684 5611 30
601 5692 5675 50
600 6632 6645 50
602 6632 6633 50
601 6633 6634 60
600 6645 6646 30
601 650 6632 200
602 6684 6652 80
601 6632 6671 200
600 6671 6684 30
600 6671 6692 30
602 6684 6611 30
601 6692 6675 50
600 7632 7645 50
602 7632 7633 50
601 7633 7634 60
600 7645 7646 30
601 650 7632 200
602 7684 7652 80
601 7632 7671 200
600 7671 7684 30
600 7671 7692 30
602 7684 7611 30
601 7692 7675 50
600 8632 8645 50
602 8632 8633 50
601 8633 8634 60
600 8645 8646 30
601 650 8632 200
602 8684 8652 80
601 8632 8671 200
600 8671 8684 30
600 8671 8692 30
602 8684 8611 30
601 8692 8675 50
0
692 0 0 0
0
634 120 90 120 90 160 110
645 40 30 50 40 60 50
646 100 50 70 40 60 40
652 20 10 60 30 40 40
671 385 220 385 220 385 220
675 68 60 290 212 285 190
692 260 200 40 40 70 40
611 80 20 60 40 30 20
2634 120 90 120 90 160 110
2645 40 30 50 40 60 50
2646 100 50 70 40 60 40
2652 20 10 60 30 40 40
2671 385 220 385 220 385 220
2675 68 60 290 212 285 190
2692 260 200 40 40 70 40
2611 80 20 60 40 30 20
3634 120 90 120 90 160 110
3645 40 30 50 40 60 50
3646 100 50 70 40 60 40
3652 20 10 60 30 40 40
3671 385 220 385 220 385 220
3675 68 60 290 212 285 190
3692 260 200 40 40 70 40
3611 80 20 60 40 30 20
4634 120 90 120 90 160 110
4645 40 30 50 40 60 50
4646 100 50 70 40 60 40
4652 20 10 60 30 40 40
4671 385 220 385 220 385 220
4675 68 60 290 212 285 190
4692 260 200 40 40 70 40
4611 80 20 60 40 30 20
5634 120 90 120 90 160 110
5645 40 30 50 40 60 50
5646 100 50 70 40 60 40
5652 20 10 60 30 40 40
5671 385 220 385 220 385 220
5675 68 60 290 212 285 190
5692 260 200 40 40 70 40
5611 80 20 60 40 30 20
6634 120 90 120 90 160 110
6645 40 30 50 40 60 50
6646 100 50 70 40 60 40
6652 20 10 60 30 40 40
6671 385 220 385 220 385 220
6675 68 60 290 212 285 190
6692 260 200 40 40 70 40
6611 80 20 60 40 30 20
7634 120 90 120 90 160 110
7645 40 30 50 40 60 50
7646 100 50 70 40 60 40
7652 20 10 60 30 40 40
7671 385 220 385 220 385 220
7675 68 60 290 212 285 190
7692 260 200 40 40 70 40
7611 80 20 60 40 30 20
8634 120 90 120 90 160 110
8645 40 30 50 40 60 50
8646 100 50 70 40 60 40
8652 20 10 60 30 40 40
8671 385 220 385 220 385 220
8675 68 60 290 212 285 190
8692 260 200 40 40 70 40
8611 80 20 60 40 30 20
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
634 160 110 120 90 120 90
645 123 102 170 125 133 110
646 17 15 23 13 21 16
652 128 86 30 10 80 70
675 85 10 68 60 20 12
611 52 42 85 60 70 40
680 12 4 10 6 7 4

14
feeder104/data1三相不平衡.txt Normal file
View File

@@ -0,0 +1,14 @@
650 4.16
0
601 650 632 2000
602 632 631 1800
602 631 630 1700
0
632 0 0 0
0
630 90 70 80 40 30 10
631 40 10 60 30 100 80
632 12 12 120 120 11 11
0

13
feeder104/data2 - 副本.txt Normal file
View File

@@ -0,0 +1,13 @@
650
0
602 632 633 500
601 650 632 2000
601 632 671 2000
601 671 680 1000
0
0
633 160 120 120 110 90 90
671 385 385 385 220 220 220
0

12
feeder104/data2.txt Normal file
View File

@@ -0,0 +1,12 @@
650 4.16
0
602 632 633 500
601 650 632 2000
601 632 671 2000
601 671 680 1000
0
0
671 385 385 385 220 220 220
0

BIN
feeder104/feeder13.zip Normal file
View File

Binary file not shown.

BIN
feeder104/ie13.TCW Normal file
View File

Binary file not shown.

27
feeder104/lineParameter - 副本.txt Normal file
View File

@@ -0,0 +1,27 @@
601
0.3465 0.1560 0.1580
0.1560 0.3375 0.1535
0.1580 0.1535 0.3414
0
1.0179 0.5017 0.4236
0.5017 1.0478 0.3849
0.4236 0.3849 1.0348
0
6.2998 -1.9958 -1.2595
-1.9958 5.9597 -0.7417
-1.2595 -0.7417 5.6386
0
602
0.7526 0.1580 0.1560
0.1580 0.7475 0.1535
0.1560 0.1535 0.7436
0
1.1814 0.4236 0.5017
0.4236 1.1983 0.3849
0.5017 0.3849 1.2112
0
5.6990 -1.0817 -1.6905
-1.0817 5.1795 -0.6588
-1.6905 -0.6588 5.4246
0

40
feeder104/lineParameter.txt Normal file
View File

@@ -0,0 +1,40 @@
600
0.3465 0.1560 0.1580
0.1560 0.3375 0.1535
0.1580 0.1535 0.3414
0
1.0179 0.5017 0.4236
0.5017 1.0478 0.3849
0.4236 0.3849 1.0348
0
6.2998 -1.9958 -1.2595
-1.9958 5.9597 -0.7417
-1.2595 -0.7417 5.6386
0
601
0.3465 0.1560 0.1580
0.1560 0.3375 0.1535
0.1580 0.1535 0.3414
0
1.0179 0.5017 0.4236
0.5017 1.0478 0.3849
0.4236 0.3849 1.0348
0
6.2998 -1.9958 -1.2595
-1.9958 5.9597 -0.7417
-1.2595 -0.7417 5.6386
0
602
0.7526 0.1580 0.1560
0.1580 0.7475 0.1535
0.1560 0.1535 0.7436
0
1.1814 0.4236 0.5017
0.4236 1.1983 0.3849
0.5017 0.3849 1.2112
0
5.6990 -1.0817 -1.6905
-1.0817 5.1795 -0.6588
-1.6905 -0.6588 5.4246
0

3
feeder104/note.txt Normal file
View File

@@ -0,0 +1,3 @@
line data
XFM-1 -1
Switch -2

204
feeder123/data.txt Normal file
View File

@@ -0,0 +1,204 @@
1 4.16
0
6 1 2 175
6 1 3 250
1 1 7 300
6 3 4 200
6 3 5 325
6 5 6 250
1 7 8 200
6 8 12 225
6 8 9 225
1 8 13 300
6 9 14 425
6 13 34 150
2 13 18 825
6 14 11 250
6 14 10 250
6 15 16 375
6 15 17 350
6 18 19 250
2 18 21 300
6 19 20 325
6 21 22 525
2 21 23 250
6 23 24 550
2 23 25 275
6 25 26 350
2 25 28 200
6 26 27 275
6 26 31 225
6 27 33 500
2 28 29 300
2 29 30 350
6 31 32 300
6 34 15 100
6 35 36 650
1 35 40 250
6 36 37 300
6 36 38 250
6 38 39 325
6 40 41 325
1 40 42 250
6 42 43 500
1 42 44 200
6 44 45 200
1 44 47 250
6 45 46 300
4 47 48 150
4 47 49 250
4 49 50 250
4 50 51 250
1 52 53 200
1 53 54 125
1 54 55 275
3 54 57 350
1 55 56 275
6 57 58 250
3 57 60 750
6 58 59 250
5 60 61 550
6 60 62 250
6 62 63 175
6 63 64 350
6 64 65 425
6 65 66 325
6 67 68 200
3 67 72 275
3 67 97 250
6 68 69 275
6 69 70 325
6 70 71 275
6 72 73 275
3 72 76 200
6 73 74 350
6 74 75 400
6 76 77 400
3 76 86 700
6 77 78 100
6 78 79 225
6 78 80 475
6 80 81 475
6 81 82 250
6 81 84 675
6 82 83 250
6 84 85 475
6 86 87 450
6 87 88 175
6 87 89 275
6 89 90 225
6 89 91 225
6 91 92 300
6 91 93 225
6 93 94 275
6 93 95 300
6 95 96 200
3 97 98 275
3 98 99 550
3 99 100 300
6 101 102 225
3 101 105 275
6 102 103 325
6 103 104 700
6 105 106 225
3 105 108 325
6 106 107 575
6 108 109 450
6 109 110 300
6 110 111 575
6 110 112 125
6 112 113 525
6 113 114 325
6 13 52 100
6 18 35 100
6 60 67 100
6 97 101 100
0
0
1 40 20 0 0 0 0
2 0 0 20 10 0 0
4 0 0 0 0 40 20
5 0 0 0 0 20 10
6 0 0 0 0 40 20
7 20 10 0 0 0 0
9 40 20 0 0 0 0
10 20 10 0 0 0 0
11 40 20 0 0 0 0
12 0 0 20 10 0 0
16 0 0 0 0 40 20
17 0 0 0 0 20 10
19 40 20 0 0 0 0
20 40 20 0 0 0 0
22 0 0 40 20 0 0
24 0 0 0 0 40 20
28 40 20 0 0 0 0
29 40 20 0 0 0 0
30 0 0 0 0 40 20
31 0 0 0 0 20 10
32 0 0 0 0 20 10
33 40 20 0 0 0 0
34 0 0 0 0 40 20
35 40 20 0 0 0 0
37 40 20 0 0 0 0
38 0 0 20 10 0 0
39 0 0 20 10 0 0
41 0 0 0 0 20 10
42 20 10 0 0 0 0
43 0 0 40 20 0 0
45 20 10 0 0 0 0
46 20 10 0 0 0 0
47 35 25 35 25 35 25
48 70 50 70 50 70 50
49 35 25 70 50 35 20
50 0 0 0 0 40 20
51 20 10 0 0 0 0
52 40 20 0 0 0 0
53 40 20 0 0 0 0
55 20 10 0 0 0 0
56 0 0 20 10 0 0
58 0 0 20 10 0 0
59 0 0 20 10 0 0
60 20 10 0 0 0 0
62 0 0 0 0 40 20
63 40 20 0 0 0 0
64 0 0 75 35 0 0
65 35 25 35 25 70 50
66 0 0 0 0 75 35
68 20 10 0 0 0 0
69 40 20 0 0 0 0
70 20 10 0 0 0 0
71 40 20 0 0 0 0
73 0 0 0 0 40 20
74 0 0 0 0 40 20
75 0 0 0 0 40 20
76 105 80 70 50 70 50
77 0 0 40 20 0 0
79 40 20 0 0 0 0
80 0 0 40 20 0 0
82 40 20 0 0 0 0
83 0 0 0 0 20 10
84 0 0 0 0 20 10
85 0 0 0 0 40 20
86 0 0 20 10 0 0
87 0 0 40 20 0 0
88 40 20 0 0 0 0
90 0 0 40 20 0 0
92 0 0 0 0 40 20
94 40 20 0 0 0 0
95 0 0 20 10 0 0
96 0 0 20 10 0 0
98 40 20 0 0 0 0
99 0 0 40 20 0 0
100 0 0 0 0 40 20
102 0 0 0 0 20 10
103 0 0 0 0 40 20
104 0 0 0 0 40 20
106 0 0 40 20 0 0
107 0 0 40 20 0 0
109 40 20 0 0 0 0
111 20 10 0 0 0 0
112 20 10 0 0 0 0
113 40 20 0 0 0 0
114 20 10 0 0 0 0
0
0

78
feeder123/lineParameter.txt Normal file
View File

@@ -0,0 +1,78 @@
1
0.4576 0.1560 0.1535
0.1560 0.4666 0.1580
0.1535 0.1580 0.4615
0
1.0780 0.5017 0.3849
0.5017 1.0482 0.4236
0.3849 0.4236 1.0651
0
5.6765 -1.8319 -0.6982
-1.8319 5.9809 -1.1645
-0.6982 -1.1645 5.3971
0
2
0.4666 0.1580 0.1560
0.1580 0.4615 0.1535
0.1560 0.1535 0.4576
0
1.0482 0.4236 0.5017
0.4236 1.0651 0.3849
0.5017 0.3849 1.0780
0
5.9809 -1.1645 -1.8319
-1.1645 5.3971 -0.6982
-1.8319 -0.6982 5.6765
0
3
0.4615 0.1535 0.1580
0.1535 0.4576 0.1560
0.1580 0.1560 0.4666
0
1.0651 0.3849 0.4236
0.3849 1.0780 0.5017
0.4236 0.5017 1.0482
0
5.3971 -0.6982 -1.1645
-0.6982 5.6765 -1.8319
-1.1645 -1.8319 5.9809
0
4
0.4615 0.1580 0.1535
0.1580 0.4666 0.1560
0.1535 0.1560 0.4576
0
1.0651 0.4236 0.3849
0.4236 1.0482 0.5017
0.3849 0.5017 1.0780
0
5.3971 -1.1645 -0.6982
-1.1645 5.9809 -1.8319
-0.6982 -1.8319 5.6765
0
5
0.4666 0.1560 0.1580
0.1560 0.4576 0.1535
0.1580 0.1535 0.4615
0
1.0482 0.5017 0.4236
0.5017 1.0780 0.3849
0.4236 0.3849 1.0651
0
5.9809 -1.8319 -1.1645
-1.8319 5.6765 -0.6982
-1.1645 -0.6982 5.3971
0
6
0.4576 0.1535 0.1560
0.1535 0.4615 0.1580
0.1560 0.1580 0.4666
0
1.0780 0.3849 0.5017
0.3849 1.0651 0.4236
0.5017 0.4236 1.0482
0
5.6765 -0.6982 -1.8319
-0.6982 5.3971 -1.1645
-1.8319 -1.1645 5.9809
0

BIN
feeder13/Cap Data.xls Normal file
View File

Binary file not shown.

BIN
feeder13/Distributed Load Data.xls Normal file
View File

Binary file not shown.

BIN
feeder13/IEEE 13 Node Test Feeder.doc Normal file
View File

Binary file not shown.

BIN
feeder13/Line Data.xls Normal file
View File

Binary file not shown.

BIN
feeder13/Regulator Data.xls Normal file
View File

Binary file not shown.

BIN
feeder13/Spot Load Data.xls Normal file
View File

Binary file not shown.

BIN
feeder13/Transformer data.xls Normal file
View File

Binary file not shown.

BIN
feeder13/UG config.xls Normal file
View File

Binary file not shown.

BIN
feeder13/config.xls Normal file
View File

Binary file not shown.

13
feeder13/data - 副本 (2) - 副本.txt Normal file
View File

@@ -0,0 +1,13 @@
650
0
602 632 633 500
601 650 632 2000
601 632 671 2000
601 671 680 1000
0
0
633 160 120 120 110 90 90
671 385 385 385 220 220 220
0

16
feeder13/data - 副本.txt Normal file
View File

@@ -0,0 +1,16 @@
650
0
600 632 645 500
602 632 633 500
601 633 634 0
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 0
602 684 611 300
601 692 675 500
0

BIN
feeder13/data.xls Normal file
View File

Binary file not shown.

57
feeder13/data1.txt Normal file
View File

@@ -0,0 +1,57 @@
650 4.16
0
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
692 0 0 0
0
634 120 90 120 90 160 110
645 40 30 50 40 60 50
646 100 50 70 40 60 40
652 20 10 60 30 40 40
671 385 220 385 220 385 220
675 68 60 290 212 285 190
692 260 200 40 40 70 40
611 80 20 60 40 30 20
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
634 160 110 120 90 120 90
645 123 102 170 125 133 110
646 17 15 23 13 21 16
652 128 86 30 10 80 70
675 85 10 68 60 20 12
611 52 42 85 60 70 40
680 12 4 10 6 7 4

14
feeder13/data1三相不平衡.txt Normal file
View File

@@ -0,0 +1,14 @@
650 4.16
0
601 650 632 2000
602 632 631 1800
602 631 630 1700
0
632 0 0 0
0
630 90 70 80 40 30 10
631 40 10 60 30 100 80
632 12 12 120 120 11 11
0

13
feeder13/data2 - 副本.txt Normal file
View File

@@ -0,0 +1,13 @@
650
0
602 632 633 500
601 650 632 2000
601 632 671 2000
601 671 680 1000
0
0
633 160 120 120 110 90 90
671 385 385 385 220 220 220
0

12
feeder13/data2.txt Normal file
View File

@@ -0,0 +1,12 @@
650 4.16
0
602 632 633 500
601 650 632 2000
601 632 671 2000
601 671 680 1000
0
0
671 385 385 385 220 220 220
0

BIN
feeder13/feeder13.zip Normal file
View File

Binary file not shown.

BIN
feeder13/ie13.TCW Normal file
View File

Binary file not shown.

27
feeder13/lineParameter - 副本.txt Normal file
View File

@@ -0,0 +1,27 @@
601
0.3465 0.1560 0.1580
0.1560 0.3375 0.1535
0.1580 0.1535 0.3414
0
1.0179 0.5017 0.4236
0.5017 1.0478 0.3849
0.4236 0.3849 1.0348
0
6.2998 -1.9958 -1.2595
-1.9958 5.9597 -0.7417
-1.2595 -0.7417 5.6386
0
602
0.7526 0.1580 0.1560
0.1580 0.7475 0.1535
0.1560 0.1535 0.7436
0
1.1814 0.4236 0.5017
0.4236 1.1983 0.3849
0.5017 0.3849 1.2112
0
5.6990 -1.0817 -1.6905
-1.0817 5.1795 -0.6588
-1.6905 -0.6588 5.4246
0

40
feeder13/lineParameter.txt Normal file
View File

@@ -0,0 +1,40 @@
600
0.3465 0.1560 0.1580
0.1560 0.3375 0.1535
0.1580 0.1535 0.3414
0
1.0179 0.5017 0.4236
0.5017 1.0478 0.3849
0.4236 0.3849 1.0348
0
6.2998 -1.9958 -1.2595
-1.9958 5.9597 -0.7417
-1.2595 -0.7417 5.6386
0
601
0.3465 0.1560 0.1580
0.1560 0.3375 0.1535
0.1580 0.1535 0.3414
0
1.0179 0.5017 0.4236
0.5017 1.0478 0.3849
0.4236 0.3849 1.0348
0
6.2998 -1.9958 -1.2595
-1.9958 5.9597 -0.7417
-1.2595 -0.7417 5.6386
0
602
0.7526 0.1580 0.1560
0.1580 0.7475 0.1535
0.1560 0.1535 0.7436
0
1.1814 0.4236 0.5017
0.4236 1.1983 0.3849
0.5017 0.3849 1.2112
0
5.6990 -1.0817 -1.6905
-1.0817 5.1795 -0.6588
-1.6905 -0.6588 5.4246
0

3
feeder13/note.txt Normal file
View File

@@ -0,0 +1,3 @@
line data
XFM-1 -1
Switch -2

BIN
feeder26/Cap Data.xls Normal file
View File

Binary file not shown.

BIN
feeder26/Distributed Load Data.xls Normal file
View File

Binary file not shown.

BIN
feeder26/IEEE 13 Node Test Feeder.doc Normal file
View File

Binary file not shown.

BIN
feeder26/Line Data.xls Normal file
View File

Binary file not shown.

BIN
feeder26/Regulator Data.xls Normal file
View File

Binary file not shown.

BIN
feeder26/Spot Load Data.xls Normal file
View File

Binary file not shown.

BIN
feeder26/Transformer data.xls Normal file
View File

Binary file not shown.

BIN
feeder26/UG config.xls Normal file
View File

Binary file not shown.

BIN
feeder26/config.xls Normal file
View File

Binary file not shown.

13
feeder26/data - 副本 (2) - 副本.txt Normal file
View File

@@ -0,0 +1,13 @@
650
0
602 632 633 500
601 650 632 2000
601 632 671 2000
601 671 680 1000
0
0
633 160 120 120 110 90 90
671 385 385 385 220 220 220
0

16
feeder26/data - 副本.txt Normal file
View File

@@ -0,0 +1,16 @@
650
0
600 632 645 500
602 632 633 500
601 633 634 0
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 0
602 684 611 300
601 692 675 500
0

BIN
feeder26/data.xls Normal file
View File

Binary file not shown.

76
feeder26/data1.txt Normal file
View File

@@ -0,0 +1,76 @@
650 4.16
0
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
600 2632 2645 50
602 2632 2633 50
601 2633 2634 60
600 2645 2646 30
601 650 2632 200
602 2684 2652 80
601 2632 2671 200
600 2671 2684 30
600 2671 2692 30
602 2684 2611 30
601 2692 2675 50
0
692 0 0 0
0
634 120 90 120 90 160 110
645 40 30 50 40 60 50
646 100 50 70 40 60 40
652 20 10 60 30 40 40
671 385 220 385 220 385 220
675 68 60 290 212 285 190
692 260 200 40 40 70 40
611 80 20 60 40 30 20
2634 120 90 120 90 160 110
2645 40 30 50 40 60 50
2646 100 50 70 40 60 40
2652 20 10 60 30 40 40
2671 385 220 385 220 385 220
2675 68 60 290 212 285 190
2692 260 200 40 40 70 40
2611 80 20 60 40 30 20
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
634 160 110 120 90 120 90
645 123 102 170 125 133 110
646 17 15 23 13 21 16
652 128 86 30 10 80 70
675 85 10 68 60 20 12
611 52 42 85 60 70 40
680 12 4 10 6 7 4

14
feeder26/data1三相不平衡.txt Normal file
View File

@@ -0,0 +1,14 @@
650 4.16
0
601 650 632 2000
602 632 631 1800
602 631 630 1700
0
632 0 0 0
0
630 90 70 80 40 30 10
631 40 10 60 30 100 80
632 12 12 120 120 11 11
0

13
feeder26/data2 - 副本.txt Normal file
View File

@@ -0,0 +1,13 @@
650
0
602 632 633 500
601 650 632 2000
601 632 671 2000
601 671 680 1000
0
0
633 160 120 120 110 90 90
671 385 385 385 220 220 220
0

12
feeder26/data2.txt Normal file
View File

@@ -0,0 +1,12 @@
650 4.16
0
602 632 633 500
601 650 632 2000
601 632 671 2000
601 671 680 1000
0
0
671 385 385 385 220 220 220
0

BIN
feeder26/feeder13.zip Normal file
View File

Binary file not shown.

BIN
feeder26/ie13.TCW Normal file
View File

Binary file not shown.

27
feeder26/lineParameter - 副本.txt Normal file
View File

@@ -0,0 +1,27 @@
601
0.3465 0.1560 0.1580
0.1560 0.3375 0.1535
0.1580 0.1535 0.3414
0
1.0179 0.5017 0.4236
0.5017 1.0478 0.3849
0.4236 0.3849 1.0348
0
6.2998 -1.9958 -1.2595
-1.9958 5.9597 -0.7417
-1.2595 -0.7417 5.6386
0
602
0.7526 0.1580 0.1560
0.1580 0.7475 0.1535
0.1560 0.1535 0.7436
0
1.1814 0.4236 0.5017
0.4236 1.1983 0.3849
0.5017 0.3849 1.2112
0
5.6990 -1.0817 -1.6905
-1.0817 5.1795 -0.6588
-1.6905 -0.6588 5.4246
0

40
feeder26/lineParameter.txt Normal file
View File

@@ -0,0 +1,40 @@
600
0.3465 0.1560 0.1580
0.1560 0.3375 0.1535
0.1580 0.1535 0.3414
0
1.0179 0.5017 0.4236
0.5017 1.0478 0.3849
0.4236 0.3849 1.0348
0
6.2998 -1.9958 -1.2595
-1.9958 5.9597 -0.7417
-1.2595 -0.7417 5.6386
0
601
0.3465 0.1560 0.1580
0.1560 0.3375 0.1535
0.1580 0.1535 0.3414
0
1.0179 0.5017 0.4236
0.5017 1.0478 0.3849
0.4236 0.3849 1.0348
0
6.2998 -1.9958 -1.2595
-1.9958 5.9597 -0.7417
-1.2595 -0.7417 5.6386
0
602
0.7526 0.1580 0.1560
0.1580 0.7475 0.1535
0.1560 0.1535 0.7436
0
1.1814 0.4236 0.5017
0.4236 1.1983 0.3849
0.5017 0.3849 1.2112
0
5.6990 -1.0817 -1.6905
-1.0817 5.1795 -0.6588
-1.6905 -0.6588 5.4246
0

3
feeder26/note.txt Normal file
View File

@@ -0,0 +1,3 @@
line data
XFM-1 -1
Switch -2

BIN
feeder39/Cap Data.xls Normal file
View File

Binary file not shown.

BIN
feeder39/Distributed Load Data.xls Normal file
View File

Binary file not shown.

BIN
feeder39/IEEE 13 Node Test Feeder.doc Normal file
View File

Binary file not shown.

BIN
feeder39/Line Data.xls Normal file
View File

Binary file not shown.

BIN
feeder39/Regulator Data.xls Normal file
View File

Binary file not shown.

BIN
feeder39/Spot Load Data.xls Normal file
View File

Binary file not shown.

BIN
feeder39/Transformer data.xls Normal file
View File

Binary file not shown.

BIN
feeder39/UG config.xls Normal file
View File

Binary file not shown.

Some files were not shown because too many files have changed in this diff Show More