commit 627ba1eb7218b019e805d73a56ea201c49675427 Author: dmy@lab Date: Fri May 8 22:02:19 2015 +0800 提交第一版 Signed-off-by: dmy@lab diff --git a/.gitignore b/.gitignore new file mode 100644 index 0000000..21229f7 --- /dev/null +++ b/.gitignore @@ -0,0 +1 @@ +*.asv \ No newline at end of file diff --git a/AssignXX.m b/AssignXX.m new file mode 100644 index 0000000..9353ff3 --- /dev/null +++ b/AssignXX.m @@ -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 \ No newline at end of file diff --git a/CalPGQG.m b/CalPGQG.m new file mode 100644 index 0000000..5e69e08 --- /dev/null +++ b/CalPGQG.m @@ -0,0 +1,25 @@ +function [ output_args ] = CalPGQG(Balance,phaseABCY,VoltpABC,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP,phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ ) +%PGQG +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 + diff --git a/FormAlphaD.m b/FormAlphaD.m new file mode 100644 index 0000000..38816b9 --- /dev/null +++ b/FormAlphaD.m @@ -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 \ No newline at end of file diff --git a/FormAlphaP.m b/FormAlphaP.m new file mode 100644 index 0000000..a5757cc --- /dev/null +++ b/FormAlphaP.m @@ -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 \ No newline at end of file diff --git a/FormG.m b/FormG.m new file mode 100644 index 0000000..d24743a --- /dev/null +++ b/FormG.m @@ -0,0 +1,8 @@ +function Mat_G=FormG(Volt,PD,QD,Loadi) + +Mat_G=[ + sparse(PD); + sparse(QD); +% Volt; + ]; +end \ No newline at end of file diff --git a/FormH.m b/FormH.m new file mode 100644 index 0000000..b0a55a8 --- /dev/null +++ b/FormH.m @@ -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 \ No newline at end of file diff --git a/FormLw.m b/FormLw.m new file mode 100644 index 0000000..3a258f8 --- /dev/null +++ b/FormLw.m @@ -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 \ No newline at end of file diff --git a/FormLx.m b/FormLx.m new file mode 100644 index 0000000..8167945 --- /dev/null +++ b/FormLx.m @@ -0,0 +1,4 @@ +function Lx=FormLx(deltF,deltH,Init_Y) +t1=deltF-deltH*Init_Y'; +Lx=t1; +end \ No newline at end of file diff --git a/FormLxComa.m b/FormLxComa.m new file mode 100644 index 0000000..c38db16 --- /dev/null +++ b/FormLxComa.m @@ -0,0 +1,4 @@ +function LxComa=FormLxComa(Lx) +%t1=deltG*(Init_Z'+Init_W');%% +LxComa=Lx; +end \ No newline at end of file diff --git a/FormLz.m b/FormLz.m new file mode 100644 index 0000000..913e331 --- /dev/null +++ b/FormLz.m @@ -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 \ No newline at end of file diff --git a/FormYY.m b/FormYY.m new file mode 100644 index 0000000..cbad106 --- /dev/null +++ b/FormYY.m @@ -0,0 +1,10 @@ +function YY=FormYY(Lul,Lz,Ly,Luu,Lw,Lx) +YY=[ + Lx; + -Ly; + -Lz; + -Lw; + -Lul; + -Luu; +]; +end \ No newline at end of file diff --git a/IEEE30.dat b/IEEE30.dat new file mode 100644 index 0000000..419d29d --- /dev/null +++ b/IEEE30.dat @@ -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 diff --git a/IEEE300.dat b/IEEE300.dat new file mode 100644 index 0000000..00a73d9 --- /dev/null +++ b/IEEE300.dat @@ -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 diff --git a/IEEE4.dat b/IEEE4.dat new file mode 100644 index 0000000..65c80a5 --- /dev/null +++ b/IEEE4.dat @@ -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 diff --git a/Jacobi.m b/Jacobi.m new file mode 100644 index 0000000..73f31fd --- /dev/null +++ b/Jacobi.m @@ -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 +%% ƽڵDz +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 \ No newline at end of file diff --git a/Modification.m b/Modification.m new file mode 100644 index 0000000..0a8b7f4 --- /dev/null +++ b/Modification.m @@ -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 \ No newline at end of file diff --git a/Modify.m b/Modify.m new file mode 100644 index 0000000..7524f42 --- /dev/null +++ b/Modify.m @@ -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 \ No newline at end of file diff --git a/OPF.m b/OPF.m new file mode 100644 index 0000000..4e17a37 --- /dev/null +++ b/OPF.m @@ -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); +%׼ +%ƽڵ01 +fsY2(Balance,:)=0; +fsY2(:,Balance)=0; +fsY2=fsY2+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum); +%ƽڵ01 +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 diff --git a/OPF_Init.m b/OPF_Init.m new file mode 100644 index 0000000..6e5beda --- /dev/null +++ b/OPF_Init.m @@ -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 \ No newline at end of file diff --git a/README.md b/README.md new file mode 100644 index 0000000..2c4fb7d --- /dev/null +++ b/README.md @@ -0,0 +1 @@ +把原来内点法的程序给改了，改成只有等式约束的普通WLS方法。 \ No newline at end of file diff --git a/Solv.m b/Solv.m new file mode 100644 index 0000000..f658297 --- /dev/null +++ b/Solv.m @@ -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 \ No newline at end of file diff --git a/SolveIt.m b/SolveIt.m new file mode 100644 index 0000000..25c291b --- /dev/null +++ b/SolveIt.m @@ -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 \ No newline at end of file diff --git a/Unbalance.m b/Unbalance.m new file mode 100644 index 0000000..e8bfcf6 --- /dev/null +++ b/Unbalance.m @@ -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 \ No newline at end of file diff --git a/admmatrix.m b/admmatrix.m new file mode 100644 index 0000000..16b536f --- /dev/null +++ b/admmatrix.m @@ -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); \ No newline at end of file diff --git a/checkSSatisfied.m b/checkSSatisfied.m new file mode 100644 index 0000000..07b83ab --- /dev/null +++ b/checkSSatisfied.m @@ -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 + diff --git a/dataRead.m b/dataRead.m new file mode 100644 index 0000000..4135a1e --- /dev/null +++ b/dataRead.m @@ -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/3IEEEϵĹʽˡ +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))); +%ļеIJݸDzҪתΪɡ +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 + diff --git a/feeder104/Cap Data.xls b/feeder104/Cap Data.xls new file mode 100644 index 0000000..7d9834d Binary files /dev/null and b/feeder104/Cap Data.xls differ diff --git a/feeder104/Distributed Load Data.xls b/feeder104/Distributed Load Data.xls new file mode 100644 index 0000000..21ec9b4 Binary files /dev/null and b/feeder104/Distributed Load Data.xls differ diff --git a/feeder104/IEEE 13 Node Test Feeder.doc b/feeder104/IEEE 13 Node Test Feeder.doc new file mode 100644 index 0000000..2be1b2a Binary files /dev/null and b/feeder104/IEEE 13 Node Test Feeder.doc differ diff --git a/feeder104/Line Data.xls b/feeder104/Line Data.xls new file mode 100644 index 0000000..2299c58 Binary files /dev/null and b/feeder104/Line Data.xls differ diff --git a/feeder104/Regulator Data.xls b/feeder104/Regulator Data.xls new file mode 100644 index 0000000..b919d6b Binary files /dev/null and b/feeder104/Regulator Data.xls differ diff --git a/feeder104/Spot Load Data.xls b/feeder104/Spot Load Data.xls new file mode 100644 index 0000000..103e81e Binary files /dev/null and b/feeder104/Spot Load Data.xls differ diff --git a/feeder104/Transformer data.xls b/feeder104/Transformer data.xls new file mode 100644 index 0000000..2bf79b5 Binary files /dev/null and b/feeder104/Transformer data.xls differ diff --git a/feeder104/UG config.xls b/feeder104/UG config.xls new file mode 100644 index 0000000..bca247a Binary files /dev/null and b/feeder104/UG config.xls differ diff --git a/feeder104/config.xls b/feeder104/config.xls new file mode 100644 index 0000000..e27b846 Binary files /dev/null and b/feeder104/config.xls differ diff --git a/feeder104/data - 副本 (2) - 副本.txt b/feeder104/data - 副本 (2) - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder104/data - 副本 (2) - 副本.txt @@ -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 + + diff --git a/feeder104/data - 副本.txt b/feeder104/data - 副本.txt new file mode 100644 index 0000000..fed6b5d --- /dev/null +++ b/feeder104/data - 副本.txt @@ -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 + diff --git a/feeder104/data.xls b/feeder104/data.xls new file mode 100644 index 0000000..8260943 Binary files /dev/null and b/feeder104/data.xls differ diff --git a/feeder104/data1.txt b/feeder104/data1.txt new file mode 100644 index 0000000..03957c8 --- /dev/null +++ b/feeder104/data1.txt @@ -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 + diff --git a/feeder104/data1三相不平衡.txt b/feeder104/data1三相不平衡.txt new file mode 100644 index 0000000..843b67a --- /dev/null +++ b/feeder104/data1三相不平衡.txt @@ -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 + + diff --git a/feeder104/data2 - 副本.txt b/feeder104/data2 - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder104/data2 - 副本.txt @@ -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 + + diff --git a/feeder104/data2.txt b/feeder104/data2.txt new file mode 100644 index 0000000..ef1d29b --- /dev/null +++ b/feeder104/data2.txt @@ -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 + + diff --git a/feeder104/feeder13.zip b/feeder104/feeder13.zip new file mode 100644 index 0000000..f8f5cca Binary files /dev/null and b/feeder104/feeder13.zip differ diff --git a/feeder104/ie13.TCW b/feeder104/ie13.TCW new file mode 100644 index 0000000..152a9d6 Binary files /dev/null and b/feeder104/ie13.TCW differ diff --git a/feeder104/lineParameter - 副本.txt b/feeder104/lineParameter - 副本.txt new file mode 100644 index 0000000..b24a66b --- /dev/null +++ b/feeder104/lineParameter - 副本.txt @@ -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 + diff --git a/feeder104/lineParameter.txt b/feeder104/lineParameter.txt new file mode 100644 index 0000000..8358036 --- /dev/null +++ b/feeder104/lineParameter.txt @@ -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 + diff --git a/feeder104/note.txt b/feeder104/note.txt new file mode 100644 index 0000000..80a3e26 --- /dev/null +++ b/feeder104/note.txt @@ -0,0 +1,3 @@ +line data +XFM-1 -1 +Switch -2 \ No newline at end of file diff --git a/feeder123/data.txt b/feeder123/data.txt new file mode 100644 index 0000000..0f9fa11 --- /dev/null +++ b/feeder123/data.txt @@ -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 diff --git a/feeder123/lineParameter.txt b/feeder123/lineParameter.txt new file mode 100644 index 0000000..bc6b5d9 --- /dev/null +++ b/feeder123/lineParameter.txt @@ -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 \ No newline at end of file diff --git a/feeder13/Cap Data.xls b/feeder13/Cap Data.xls new file mode 100644 index 0000000..7d9834d Binary files /dev/null and b/feeder13/Cap Data.xls differ diff --git a/feeder13/Distributed Load Data.xls b/feeder13/Distributed Load Data.xls new file mode 100644 index 0000000..21ec9b4 Binary files /dev/null and b/feeder13/Distributed Load Data.xls differ diff --git a/feeder13/IEEE 13 Node Test Feeder.doc b/feeder13/IEEE 13 Node Test Feeder.doc new file mode 100644 index 0000000..2be1b2a Binary files /dev/null and b/feeder13/IEEE 13 Node Test Feeder.doc differ diff --git a/feeder13/Line Data.xls b/feeder13/Line Data.xls new file mode 100644 index 0000000..2299c58 Binary files /dev/null and b/feeder13/Line Data.xls differ diff --git a/feeder13/Regulator Data.xls b/feeder13/Regulator Data.xls new file mode 100644 index 0000000..b919d6b Binary files /dev/null and b/feeder13/Regulator Data.xls differ diff --git a/feeder13/Spot Load Data.xls b/feeder13/Spot Load Data.xls new file mode 100644 index 0000000..103e81e Binary files /dev/null and b/feeder13/Spot Load Data.xls differ diff --git a/feeder13/Transformer data.xls b/feeder13/Transformer data.xls new file mode 100644 index 0000000..2bf79b5 Binary files /dev/null and b/feeder13/Transformer data.xls differ diff --git a/feeder13/UG config.xls b/feeder13/UG config.xls new file mode 100644 index 0000000..bca247a Binary files /dev/null and b/feeder13/UG config.xls differ diff --git a/feeder13/config.xls b/feeder13/config.xls new file mode 100644 index 0000000..e27b846 Binary files /dev/null and b/feeder13/config.xls differ diff --git a/feeder13/data - 副本 (2) - 副本.txt b/feeder13/data - 副本 (2) - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder13/data - 副本 (2) - 副本.txt @@ -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 + + diff --git a/feeder13/data - 副本.txt b/feeder13/data - 副本.txt new file mode 100644 index 0000000..fed6b5d --- /dev/null +++ b/feeder13/data - 副本.txt @@ -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 + diff --git a/feeder13/data.xls b/feeder13/data.xls new file mode 100644 index 0000000..8260943 Binary files /dev/null and b/feeder13/data.xls differ diff --git a/feeder13/data1.txt b/feeder13/data1.txt new file mode 100644 index 0000000..27e0b6e --- /dev/null +++ b/feeder13/data1.txt @@ -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 + diff --git a/feeder13/data1三相不平衡.txt b/feeder13/data1三相不平衡.txt new file mode 100644 index 0000000..843b67a --- /dev/null +++ b/feeder13/data1三相不平衡.txt @@ -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 + + diff --git a/feeder13/data2 - 副本.txt b/feeder13/data2 - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder13/data2 - 副本.txt @@ -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 + + diff --git a/feeder13/data2.txt b/feeder13/data2.txt new file mode 100644 index 0000000..ef1d29b --- /dev/null +++ b/feeder13/data2.txt @@ -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 + + diff --git a/feeder13/feeder13.zip b/feeder13/feeder13.zip new file mode 100644 index 0000000..f8f5cca Binary files /dev/null and b/feeder13/feeder13.zip differ diff --git a/feeder13/ie13.TCW b/feeder13/ie13.TCW new file mode 100644 index 0000000..152a9d6 Binary files /dev/null and b/feeder13/ie13.TCW differ diff --git a/feeder13/lineParameter - 副本.txt b/feeder13/lineParameter - 副本.txt new file mode 100644 index 0000000..b24a66b --- /dev/null +++ b/feeder13/lineParameter - 副本.txt @@ -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 + diff --git a/feeder13/lineParameter.txt b/feeder13/lineParameter.txt new file mode 100644 index 0000000..8358036 --- /dev/null +++ b/feeder13/lineParameter.txt @@ -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 + diff --git a/feeder13/note.txt b/feeder13/note.txt new file mode 100644 index 0000000..80a3e26 --- /dev/null +++ b/feeder13/note.txt @@ -0,0 +1,3 @@ +line data +XFM-1 -1 +Switch -2 \ No newline at end of file diff --git a/feeder26/Cap Data.xls b/feeder26/Cap Data.xls new file mode 100644 index 0000000..7d9834d Binary files /dev/null and b/feeder26/Cap Data.xls differ diff --git a/feeder26/Distributed Load Data.xls b/feeder26/Distributed Load Data.xls new file mode 100644 index 0000000..21ec9b4 Binary files /dev/null and b/feeder26/Distributed Load Data.xls differ diff --git a/feeder26/IEEE 13 Node Test Feeder.doc b/feeder26/IEEE 13 Node Test Feeder.doc new file mode 100644 index 0000000..2be1b2a Binary files /dev/null and b/feeder26/IEEE 13 Node Test Feeder.doc differ diff --git a/feeder26/Line Data.xls b/feeder26/Line Data.xls new file mode 100644 index 0000000..2299c58 Binary files /dev/null and b/feeder26/Line Data.xls differ diff --git a/feeder26/Regulator Data.xls b/feeder26/Regulator Data.xls new file mode 100644 index 0000000..b919d6b Binary files /dev/null and b/feeder26/Regulator Data.xls differ diff --git a/feeder26/Spot Load Data.xls b/feeder26/Spot Load Data.xls new file mode 100644 index 0000000..103e81e Binary files /dev/null and b/feeder26/Spot Load Data.xls differ diff --git a/feeder26/Transformer data.xls b/feeder26/Transformer data.xls new file mode 100644 index 0000000..2bf79b5 Binary files /dev/null and b/feeder26/Transformer data.xls differ diff --git a/feeder26/UG config.xls b/feeder26/UG config.xls new file mode 100644 index 0000000..bca247a Binary files /dev/null and b/feeder26/UG config.xls differ diff --git a/feeder26/config.xls b/feeder26/config.xls new file mode 100644 index 0000000..e27b846 Binary files /dev/null and b/feeder26/config.xls differ diff --git a/feeder26/data - 副本 (2) - 副本.txt b/feeder26/data - 副本 (2) - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder26/data - 副本 (2) - 副本.txt @@ -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 + + diff --git a/feeder26/data - 副本.txt b/feeder26/data - 副本.txt new file mode 100644 index 0000000..fed6b5d --- /dev/null +++ b/feeder26/data - 副本.txt @@ -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 + diff --git a/feeder26/data.xls b/feeder26/data.xls new file mode 100644 index 0000000..8260943 Binary files /dev/null and b/feeder26/data.xls differ diff --git a/feeder26/data1.txt b/feeder26/data1.txt new file mode 100644 index 0000000..ad27e57 --- /dev/null +++ b/feeder26/data1.txt @@ -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 + diff --git a/feeder26/data1三相不平衡.txt b/feeder26/data1三相不平衡.txt new file mode 100644 index 0000000..843b67a --- /dev/null +++ b/feeder26/data1三相不平衡.txt @@ -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 + + diff --git a/feeder26/data2 - 副本.txt b/feeder26/data2 - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder26/data2 - 副本.txt @@ -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 + + diff --git a/feeder26/data2.txt b/feeder26/data2.txt new file mode 100644 index 0000000..ef1d29b --- /dev/null +++ b/feeder26/data2.txt @@ -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 + + diff --git a/feeder26/feeder13.zip b/feeder26/feeder13.zip new file mode 100644 index 0000000..f8f5cca Binary files /dev/null and b/feeder26/feeder13.zip differ diff --git a/feeder26/ie13.TCW b/feeder26/ie13.TCW new file mode 100644 index 0000000..152a9d6 Binary files /dev/null and b/feeder26/ie13.TCW differ diff --git a/feeder26/lineParameter - 副本.txt b/feeder26/lineParameter - 副本.txt new file mode 100644 index 0000000..b24a66b --- /dev/null +++ b/feeder26/lineParameter - 副本.txt @@ -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 + diff --git a/feeder26/lineParameter.txt b/feeder26/lineParameter.txt new file mode 100644 index 0000000..8358036 --- /dev/null +++ b/feeder26/lineParameter.txt @@ -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 + diff --git a/feeder26/note.txt b/feeder26/note.txt new file mode 100644 index 0000000..80a3e26 --- /dev/null +++ b/feeder26/note.txt @@ -0,0 +1,3 @@ +line data +XFM-1 -1 +Switch -2 \ No newline at end of file diff --git a/feeder39/Cap Data.xls b/feeder39/Cap Data.xls new file mode 100644 index 0000000..7d9834d Binary files /dev/null and b/feeder39/Cap Data.xls differ diff --git a/feeder39/Distributed Load Data.xls b/feeder39/Distributed Load Data.xls new file mode 100644 index 0000000..21ec9b4 Binary files /dev/null and b/feeder39/Distributed Load Data.xls differ diff --git a/feeder39/IEEE 13 Node Test Feeder.doc b/feeder39/IEEE 13 Node Test Feeder.doc new file mode 100644 index 0000000..2be1b2a Binary files /dev/null and b/feeder39/IEEE 13 Node Test Feeder.doc differ diff --git a/feeder39/Line Data.xls b/feeder39/Line Data.xls new file mode 100644 index 0000000..2299c58 Binary files /dev/null and b/feeder39/Line Data.xls differ diff --git a/feeder39/Regulator Data.xls b/feeder39/Regulator Data.xls new file mode 100644 index 0000000..b919d6b Binary files /dev/null and b/feeder39/Regulator Data.xls differ diff --git a/feeder39/Spot Load Data.xls b/feeder39/Spot Load Data.xls new file mode 100644 index 0000000..103e81e Binary files /dev/null and b/feeder39/Spot Load Data.xls differ diff --git a/feeder39/Transformer data.xls b/feeder39/Transformer data.xls new file mode 100644 index 0000000..2bf79b5 Binary files /dev/null and b/feeder39/Transformer data.xls differ diff --git a/feeder39/UG config.xls b/feeder39/UG config.xls new file mode 100644 index 0000000..bca247a Binary files /dev/null and b/feeder39/UG config.xls differ diff --git a/feeder39/config.xls b/feeder39/config.xls new file mode 100644 index 0000000..e27b846 Binary files /dev/null and b/feeder39/config.xls differ diff --git a/feeder39/data - 副本 (2) - 副本.txt b/feeder39/data - 副本 (2) - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder39/data - 副本 (2) - 副本.txt @@ -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 + + diff --git a/feeder39/data - 副本.txt b/feeder39/data - 副本.txt new file mode 100644 index 0000000..fed6b5d --- /dev/null +++ b/feeder39/data - 副本.txt @@ -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 + diff --git a/feeder39/data.xls b/feeder39/data.xls new file mode 100644 index 0000000..8260943 Binary files /dev/null and b/feeder39/data.xls differ diff --git a/feeder39/data1.txt b/feeder39/data1.txt new file mode 100644 index 0000000..809b3c7 --- /dev/null +++ b/feeder39/data1.txt @@ -0,0 +1,95 @@ +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 +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 +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 + diff --git a/feeder39/data1三相不平衡.txt b/feeder39/data1三相不平衡.txt new file mode 100644 index 0000000..843b67a --- /dev/null +++ b/feeder39/data1三相不平衡.txt @@ -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 + + diff --git a/feeder39/data2 - 副本.txt b/feeder39/data2 - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder39/data2 - 副本.txt @@ -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 + + diff --git a/feeder39/data2.txt b/feeder39/data2.txt new file mode 100644 index 0000000..ef1d29b --- /dev/null +++ b/feeder39/data2.txt @@ -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 + + diff --git a/feeder39/feeder13.zip b/feeder39/feeder13.zip new file mode 100644 index 0000000..f8f5cca Binary files /dev/null and b/feeder39/feeder13.zip differ diff --git a/feeder39/ie13.TCW b/feeder39/ie13.TCW new file mode 100644 index 0000000..152a9d6 Binary files /dev/null and b/feeder39/ie13.TCW differ diff --git a/feeder39/lineParameter - 副本.txt b/feeder39/lineParameter - 副本.txt new file mode 100644 index 0000000..b24a66b --- /dev/null +++ b/feeder39/lineParameter - 副本.txt @@ -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 + diff --git a/feeder39/lineParameter.txt b/feeder39/lineParameter.txt new file mode 100644 index 0000000..8358036 --- /dev/null +++ b/feeder39/lineParameter.txt @@ -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 + diff --git a/feeder39/note.txt b/feeder39/note.txt new file mode 100644 index 0000000..80a3e26 --- /dev/null +++ b/feeder39/note.txt @@ -0,0 +1,3 @@ +line data +XFM-1 -1 +Switch -2 \ No newline at end of file diff --git a/feeder52/Cap Data.xls b/feeder52/Cap Data.xls new file mode 100644 index 0000000..7d9834d Binary files /dev/null and b/feeder52/Cap Data.xls differ diff --git a/feeder52/Distributed Load Data.xls b/feeder52/Distributed Load Data.xls new file mode 100644 index 0000000..21ec9b4 Binary files /dev/null and b/feeder52/Distributed Load Data.xls differ diff --git a/feeder52/IEEE 13 Node Test Feeder.doc b/feeder52/IEEE 13 Node Test Feeder.doc new file mode 100644 index 0000000..2be1b2a Binary files /dev/null and b/feeder52/IEEE 13 Node Test Feeder.doc differ diff --git a/feeder52/Line Data.xls b/feeder52/Line Data.xls new file mode 100644 index 0000000..2299c58 Binary files /dev/null and b/feeder52/Line Data.xls differ diff --git a/feeder52/Regulator Data.xls b/feeder52/Regulator Data.xls new file mode 100644 index 0000000..b919d6b Binary files /dev/null and b/feeder52/Regulator Data.xls differ diff --git a/feeder52/Spot Load Data.xls b/feeder52/Spot Load Data.xls new file mode 100644 index 0000000..103e81e Binary files /dev/null and b/feeder52/Spot Load Data.xls differ diff --git a/feeder52/Transformer data.xls b/feeder52/Transformer data.xls new file mode 100644 index 0000000..2bf79b5 Binary files /dev/null and b/feeder52/Transformer data.xls differ diff --git a/feeder52/UG config.xls b/feeder52/UG config.xls new file mode 100644 index 0000000..bca247a Binary files /dev/null and b/feeder52/UG config.xls differ diff --git a/feeder52/config.xls b/feeder52/config.xls new file mode 100644 index 0000000..e27b846 Binary files /dev/null and b/feeder52/config.xls differ diff --git a/feeder52/data - 副本 (2) - 副本.txt b/feeder52/data - 副本 (2) - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder52/data - 副本 (2) - 副本.txt @@ -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 + + diff --git a/feeder52/data - 副本.txt b/feeder52/data - 副本.txt new file mode 100644 index 0000000..fed6b5d --- /dev/null +++ b/feeder52/data - 副本.txt @@ -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 + diff --git a/feeder52/data.xls b/feeder52/data.xls new file mode 100644 index 0000000..8260943 Binary files /dev/null and b/feeder52/data.xls differ diff --git a/feeder52/data1.txt b/feeder52/data1.txt new file mode 100644 index 0000000..52691ec --- /dev/null +++ b/feeder52/data1.txt @@ -0,0 +1,114 @@ +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 +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 +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 + diff --git a/feeder52/data1三相不平衡.txt b/feeder52/data1三相不平衡.txt new file mode 100644 index 0000000..843b67a --- /dev/null +++ b/feeder52/data1三相不平衡.txt @@ -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 + + diff --git a/feeder52/data2 - 副本.txt b/feeder52/data2 - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder52/data2 - 副本.txt @@ -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 + + diff --git a/feeder52/data2.txt b/feeder52/data2.txt new file mode 100644 index 0000000..ef1d29b --- /dev/null +++ b/feeder52/data2.txt @@ -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 + + diff --git a/feeder52/feeder13.zip b/feeder52/feeder13.zip new file mode 100644 index 0000000..f8f5cca Binary files /dev/null and b/feeder52/feeder13.zip differ diff --git a/feeder52/ie13.TCW b/feeder52/ie13.TCW new file mode 100644 index 0000000..152a9d6 Binary files /dev/null and b/feeder52/ie13.TCW differ diff --git a/feeder52/lineParameter - 副本.txt b/feeder52/lineParameter - 副本.txt new file mode 100644 index 0000000..b24a66b --- /dev/null +++ b/feeder52/lineParameter - 副本.txt @@ -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 + diff --git a/feeder52/lineParameter.txt b/feeder52/lineParameter.txt new file mode 100644 index 0000000..8358036 --- /dev/null +++ b/feeder52/lineParameter.txt @@ -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 + diff --git a/feeder52/note.txt b/feeder52/note.txt new file mode 100644 index 0000000..80a3e26 --- /dev/null +++ b/feeder52/note.txt @@ -0,0 +1,3 @@ +line data +XFM-1 -1 +Switch -2 \ No newline at end of file diff --git a/feeder65/Cap Data.xls b/feeder65/Cap Data.xls new file mode 100644 index 0000000..7d9834d Binary files /dev/null and b/feeder65/Cap Data.xls differ diff --git a/feeder65/Distributed Load Data.xls b/feeder65/Distributed Load Data.xls new file mode 100644 index 0000000..21ec9b4 Binary files /dev/null and b/feeder65/Distributed Load Data.xls differ diff --git a/feeder65/IEEE 13 Node Test Feeder.doc b/feeder65/IEEE 13 Node Test Feeder.doc new file mode 100644 index 0000000..2be1b2a Binary files /dev/null and b/feeder65/IEEE 13 Node Test Feeder.doc differ diff --git a/feeder65/Line Data.xls b/feeder65/Line Data.xls new file mode 100644 index 0000000..2299c58 Binary files /dev/null and b/feeder65/Line Data.xls differ diff --git a/feeder65/Regulator Data.xls b/feeder65/Regulator Data.xls new file mode 100644 index 0000000..b919d6b Binary files /dev/null and b/feeder65/Regulator Data.xls differ diff --git a/feeder65/Spot Load Data.xls b/feeder65/Spot Load Data.xls new file mode 100644 index 0000000..103e81e Binary files /dev/null and b/feeder65/Spot Load Data.xls differ diff --git a/feeder65/Transformer data.xls b/feeder65/Transformer data.xls new file mode 100644 index 0000000..2bf79b5 Binary files /dev/null and b/feeder65/Transformer data.xls differ diff --git a/feeder65/UG config.xls b/feeder65/UG config.xls new file mode 100644 index 0000000..bca247a Binary files /dev/null and b/feeder65/UG config.xls differ diff --git a/feeder65/config.xls b/feeder65/config.xls new file mode 100644 index 0000000..e27b846 Binary files /dev/null and b/feeder65/config.xls differ diff --git a/feeder65/data - 副本 (2) - 副本.txt b/feeder65/data - 副本 (2) - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder65/data - 副本 (2) - 副本.txt @@ -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 + + diff --git a/feeder65/data - 副本.txt b/feeder65/data - 副本.txt new file mode 100644 index 0000000..fed6b5d --- /dev/null +++ b/feeder65/data - 副本.txt @@ -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 + diff --git a/feeder65/data.xls b/feeder65/data.xls new file mode 100644 index 0000000..8260943 Binary files /dev/null and b/feeder65/data.xls differ diff --git a/feeder65/data1.txt b/feeder65/data1.txt new file mode 100644 index 0000000..5700d93 --- /dev/null +++ b/feeder65/data1.txt @@ -0,0 +1,133 @@ +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 +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 +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 + diff --git a/feeder65/data1三相不平衡.txt b/feeder65/data1三相不平衡.txt new file mode 100644 index 0000000..843b67a --- /dev/null +++ b/feeder65/data1三相不平衡.txt @@ -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 + + diff --git a/feeder65/data2 - 副本.txt b/feeder65/data2 - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder65/data2 - 副本.txt @@ -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 + + diff --git a/feeder65/data2.txt b/feeder65/data2.txt new file mode 100644 index 0000000..ef1d29b --- /dev/null +++ b/feeder65/data2.txt @@ -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 + + diff --git a/feeder65/feeder13.zip b/feeder65/feeder13.zip new file mode 100644 index 0000000..f8f5cca Binary files /dev/null and b/feeder65/feeder13.zip differ diff --git a/feeder65/ie13.TCW b/feeder65/ie13.TCW new file mode 100644 index 0000000..152a9d6 Binary files /dev/null and b/feeder65/ie13.TCW differ diff --git a/feeder65/lineParameter - 副本.txt b/feeder65/lineParameter - 副本.txt new file mode 100644 index 0000000..b24a66b --- /dev/null +++ b/feeder65/lineParameter - 副本.txt @@ -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 + diff --git a/feeder65/lineParameter.txt b/feeder65/lineParameter.txt new file mode 100644 index 0000000..8358036 --- /dev/null +++ b/feeder65/lineParameter.txt @@ -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 + diff --git a/feeder65/note.txt b/feeder65/note.txt new file mode 100644 index 0000000..80a3e26 --- /dev/null +++ b/feeder65/note.txt @@ -0,0 +1,3 @@ +line data +XFM-1 -1 +Switch -2 \ No newline at end of file diff --git a/feeder78/Cap Data.xls b/feeder78/Cap Data.xls new file mode 100644 index 0000000..7d9834d Binary files /dev/null and b/feeder78/Cap Data.xls differ diff --git a/feeder78/Distributed Load Data.xls b/feeder78/Distributed Load Data.xls new file mode 100644 index 0000000..21ec9b4 Binary files /dev/null and b/feeder78/Distributed Load Data.xls differ diff --git a/feeder78/IEEE 13 Node Test Feeder.doc b/feeder78/IEEE 13 Node Test Feeder.doc new file mode 100644 index 0000000..2be1b2a Binary files /dev/null and b/feeder78/IEEE 13 Node Test Feeder.doc differ diff --git a/feeder78/Line Data.xls b/feeder78/Line Data.xls new file mode 100644 index 0000000..2299c58 Binary files /dev/null and b/feeder78/Line Data.xls differ diff --git a/feeder78/Regulator Data.xls b/feeder78/Regulator Data.xls new file mode 100644 index 0000000..b919d6b Binary files /dev/null and b/feeder78/Regulator Data.xls differ diff --git a/feeder78/Spot Load Data.xls b/feeder78/Spot Load Data.xls new file mode 100644 index 0000000..103e81e Binary files /dev/null and b/feeder78/Spot Load Data.xls differ diff --git a/feeder78/Transformer data.xls b/feeder78/Transformer data.xls new file mode 100644 index 0000000..2bf79b5 Binary files /dev/null and b/feeder78/Transformer data.xls differ diff --git a/feeder78/UG config.xls b/feeder78/UG config.xls new file mode 100644 index 0000000..bca247a Binary files /dev/null and b/feeder78/UG config.xls differ diff --git a/feeder78/config.xls b/feeder78/config.xls new file mode 100644 index 0000000..e27b846 Binary files /dev/null and b/feeder78/config.xls differ diff --git a/feeder78/data - 副本 (2) - 副本.txt b/feeder78/data - 副本 (2) - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder78/data - 副本 (2) - 副本.txt @@ -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 + + diff --git a/feeder78/data - 副本.txt b/feeder78/data - 副本.txt new file mode 100644 index 0000000..fed6b5d --- /dev/null +++ b/feeder78/data - 副本.txt @@ -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 + diff --git a/feeder78/data.xls b/feeder78/data.xls new file mode 100644 index 0000000..8260943 Binary files /dev/null and b/feeder78/data.xls differ diff --git a/feeder78/data1.txt b/feeder78/data1.txt new file mode 100644 index 0000000..386d506 --- /dev/null +++ b/feeder78/data1.txt @@ -0,0 +1,152 @@ +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 +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 +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 + diff --git a/feeder78/data1三相不平衡.txt b/feeder78/data1三相不平衡.txt new file mode 100644 index 0000000..843b67a --- /dev/null +++ b/feeder78/data1三相不平衡.txt @@ -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 + + diff --git a/feeder78/data2 - 副本.txt b/feeder78/data2 - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder78/data2 - 副本.txt @@ -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 + + diff --git a/feeder78/data2.txt b/feeder78/data2.txt new file mode 100644 index 0000000..ef1d29b --- /dev/null +++ b/feeder78/data2.txt @@ -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 + + diff --git a/feeder78/feeder13.zip b/feeder78/feeder13.zip new file mode 100644 index 0000000..f8f5cca Binary files /dev/null and b/feeder78/feeder13.zip differ diff --git a/feeder78/ie13.TCW b/feeder78/ie13.TCW new file mode 100644 index 0000000..152a9d6 Binary files /dev/null and b/feeder78/ie13.TCW differ diff --git a/feeder78/lineParameter - 副本.txt b/feeder78/lineParameter - 副本.txt new file mode 100644 index 0000000..b24a66b --- /dev/null +++ b/feeder78/lineParameter - 副本.txt @@ -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 + diff --git a/feeder78/lineParameter.txt b/feeder78/lineParameter.txt new file mode 100644 index 0000000..8358036 --- /dev/null +++ b/feeder78/lineParameter.txt @@ -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 + diff --git a/feeder78/note.txt b/feeder78/note.txt new file mode 100644 index 0000000..80a3e26 --- /dev/null +++ b/feeder78/note.txt @@ -0,0 +1,3 @@ +line data +XFM-1 -1 +Switch -2 \ No newline at end of file diff --git a/feeder91/Cap Data.xls b/feeder91/Cap Data.xls new file mode 100644 index 0000000..7d9834d Binary files /dev/null and b/feeder91/Cap Data.xls differ diff --git a/feeder91/Distributed Load Data.xls b/feeder91/Distributed Load Data.xls new file mode 100644 index 0000000..21ec9b4 Binary files /dev/null and b/feeder91/Distributed Load Data.xls differ diff --git a/feeder91/IEEE 13 Node Test Feeder.doc b/feeder91/IEEE 13 Node Test Feeder.doc new file mode 100644 index 0000000..2be1b2a Binary files /dev/null and b/feeder91/IEEE 13 Node Test Feeder.doc differ diff --git a/feeder91/Line Data.xls b/feeder91/Line Data.xls new file mode 100644 index 0000000..2299c58 Binary files /dev/null and b/feeder91/Line Data.xls differ diff --git a/feeder91/Regulator Data.xls b/feeder91/Regulator Data.xls new file mode 100644 index 0000000..b919d6b Binary files /dev/null and b/feeder91/Regulator Data.xls differ diff --git a/feeder91/Spot Load Data.xls b/feeder91/Spot Load Data.xls new file mode 100644 index 0000000..103e81e Binary files /dev/null and b/feeder91/Spot Load Data.xls differ diff --git a/feeder91/Transformer data.xls b/feeder91/Transformer data.xls new file mode 100644 index 0000000..2bf79b5 Binary files /dev/null and b/feeder91/Transformer data.xls differ diff --git a/feeder91/UG config.xls b/feeder91/UG config.xls new file mode 100644 index 0000000..bca247a Binary files /dev/null and b/feeder91/UG config.xls differ diff --git a/feeder91/config.xls b/feeder91/config.xls new file mode 100644 index 0000000..e27b846 Binary files /dev/null and b/feeder91/config.xls differ diff --git a/feeder91/data - 副本 (2) - 副本.txt b/feeder91/data - 副本 (2) - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder91/data - 副本 (2) - 副本.txt @@ -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 + + diff --git a/feeder91/data - 副本.txt b/feeder91/data - 副本.txt new file mode 100644 index 0000000..fed6b5d --- /dev/null +++ b/feeder91/data - 副本.txt @@ -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 + diff --git a/feeder91/data.xls b/feeder91/data.xls new file mode 100644 index 0000000..8260943 Binary files /dev/null and b/feeder91/data.xls differ diff --git a/feeder91/data1.txt b/feeder91/data1.txt new file mode 100644 index 0000000..f790669 --- /dev/null +++ b/feeder91/data1.txt @@ -0,0 +1,171 @@ +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 +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 +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 + diff --git a/feeder91/data1三相不平衡.txt b/feeder91/data1三相不平衡.txt new file mode 100644 index 0000000..843b67a --- /dev/null +++ b/feeder91/data1三相不平衡.txt @@ -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 + + diff --git a/feeder91/data2 - 副本.txt b/feeder91/data2 - 副本.txt new file mode 100644 index 0000000..2cf4bae --- /dev/null +++ b/feeder91/data2 - 副本.txt @@ -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 + + diff --git a/feeder91/data2.txt b/feeder91/data2.txt new file mode 100644 index 0000000..ef1d29b --- /dev/null +++ b/feeder91/data2.txt @@ -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 + + diff --git a/feeder91/feeder13.zip b/feeder91/feeder13.zip new file mode 100644 index 0000000..f8f5cca Binary files /dev/null and b/feeder91/feeder13.zip differ diff --git a/feeder91/ie13.TCW b/feeder91/ie13.TCW new file mode 100644 index 0000000..152a9d6 Binary files /dev/null and b/feeder91/ie13.TCW differ diff --git a/feeder91/lineParameter - 副本.txt b/feeder91/lineParameter - 副本.txt new file mode 100644 index 0000000..b24a66b --- /dev/null +++ b/feeder91/lineParameter - 副本.txt @@ -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 + diff --git a/feeder91/lineParameter.txt b/feeder91/lineParameter.txt new file mode 100644 index 0000000..8358036 --- /dev/null +++ b/feeder91/lineParameter.txt @@ -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 + diff --git a/feeder91/note.txt b/feeder91/note.txt new file mode 100644 index 0000000..80a3e26 --- /dev/null +++ b/feeder91/note.txt @@ -0,0 +1,3 @@ +line data +XFM-1 -1 +Switch -2 \ No newline at end of file diff --git a/func_ddf.m b/func_ddf.m new file mode 100644 index 0000000..594161d --- /dev/null +++ b/func_ddf.m @@ -0,0 +1,176 @@ +function ddf=func_ddh(AngleIJMat,GB,Volt,Init_Y,Busnum) +%% deltaPi/deltaThytai_deltaThytaj ǶԽԪ +t1=-Volt'*Volt; +t2=real(GB).*cos(AngleIJMat)+imag(GB).*sin(AngleIJMat); +dPdTidTj=t1.*t2; %%(˶ԽԪص) +dPidTidTj_2=dPdTidTj-diag(diag(dPdTidTj));%ȥ˶ԽԪص +t3=repmat(Init_Y',1,size(Init_Y,2));%yϵ +dPidTidTj_2=dPidTidTj_2.*t3(1:2:size(Init_Y,2),1:2:size(Init_Y,2)); +t2=real(GB).*cos(AngleIJMat)-imag(GB).*sin(AngleIJMat); +t3=t1.*t2; +dPjdTidTj=t3-diag(diag(t3));%ȥ˶ԽԪص +t3=repmat(Init_Y,size(Init_Y,2),1); +dPjdTidTj=dPjdTidTj.*t3(1:2:size(Init_Y,2),1:2:size(Init_Y,2));%yϵ +dPdTidTj=dPidTidTj_2+dPjdTidTj;%շǶԽԪ @@ +%% deltaP/deltaThyta_deltaThyta ԽԪ +t1=sum(-dPidTidTj_2,2); +t2=diag(t1'.*Init_Y(1:2:size(Init_Y,2)));%yϵ +t3=sum(-dPidTidTj_2,1); +t4=diag(t3.*Init_Y(1:2:size(Init_Y,2)));%yϵ +dPdTidTi=t2+t4;%%նԽԪ @@ +%% deltaP/deltaThytai_dVi ԽԪ +t1=ones(Busnum,1)*Volt.*(real(GB).*sin(AngleIJMat)-imag(GB).*cos(AngleIJMat)); +t2=diag(t1); +t3=t1-diag(t2);%ȥ˶ԽԪص +t4=sum(t3,2); +t4=t4'.*Init_Y(1:2:size(Init_Y,2));%yϵ +dPidTidVi=diag(t4); +dPidTjdVj=-t3; +t6=sum(dPidTjdVj,1);%yϵ +t6=t6.*Init_Y(1:2:size(Init_Y,2));%yϵ +dPdTidVi=dPidTidVi+diag(t6);%%նԽԪ @@ +%% deltaP/deltaThytai_dVj ǶԽԪ +t1=ones(Busnum,1)*Volt.*(real(GB).*sin(AngleIJMat)-imag(GB).*cos(AngleIJMat)); +dPidTidVj=t1-diag(diag(t1));%%ȥԽԪص +t2=repmat(Init_Y',1,size(Init_Y,2)); +dPidTidVj=dPidTidVj.*t2(1:2:size(Init_Y,2),1:2:size(Init_Y,2));%yϵ +t1=ones(Busnum,1)*Volt.*real(GB).*(sin(AngleIJMat)+imag(GB).*cos(AngleIJMat)); +dPjdTidVj=t1-diag(diag(t1));%ȥԽԪص +t2=t2'; +dPjdTidVj=dPjdTidVj.*t2(1:2:size(Init_Y,2),1:2:size(Init_Y,2));%yϵ +dPdTidVj=dPidTidVj+dPjdTidVj;%շǶԽԪ @@ +%% deltaP/dVi_deltaThytaj ǶԽԪ +t1=-ones(Busnum,1)*Volt.*real(GB).*(sin(AngleIJMat)-imag(GB).*cos(AngleIJMat)); +dPidVidTj=t1-diag(diag(t1)); +t2=repmat(Init_Y',1,size(Init_Y,2)); +dPidVidTj=dPidVidTj.*t2(1:2:size(Init_Y,2),1:2:size(Init_Y,2));%yϵ +t1=-ones(Busnum,1)*Volt.*(real(GB).*sin(AngleIJMat)+imag(GB).*cos(AngleIJMat)); +dPjdVidTj=t1-diag(diag(t1)); +t2=t2'; +dPjdVidTj.*t2(1:2:size(Init_Y,2),1:2:size(Init_Y,2));%yϵ +dPdVidTj=dPidVidTj+dPjdVidTj;%շǶԽԪ +%% deltaPi/dVi_deltaThyta ԽԪ +dPdVidTi=dPdTidVi;%նԽԪ @@ +%% deltaP/dVi_dVj ǶԽԪ +t1=-(real(GB).*(cos(AngleIJMat)+imag(GB).*sin(AngleIJMat))); +dPidVidVj=t1-diag(diag(t1));%ȥԽԪص +t2=repmat(Init_Y',1,size(Init_Y,2)); +dPidVidVj=dPidVidVj.*t2(1:2:size(Init_Y,2),1:2:size(Init_Y,2));%yϵ +t1=-(real(GB).*(cos(AngleIJMat)-imag(GB).*sin(AngleIJMat))); +dPjdVidVj=t1-diag(diag(t1)); +t2=t2'; +dPjdVidVj=dPjdVidVj.*t2(1:2:size(Init_Y,2),1:2:size(Init_Y,2));%yϵ +dPdVidVj=dPidVidVj+dPjdVidVj;%շǶԽԪ @@ +%% deltaP/dVi_dVi ԽԪ +t1=diag(real(GB)); +t2=t1'.*Init_Y(1:2:size(Init_Y,2));%yϵ +dPidVidVi=-2*diag(t2); +dPidVjdVj=0; +dPdVidVi=dPidVidVi+dPidVjdVj;%նԽԪ @@ +%% APi +APi=zeros(2*Busnum,2*Busnum); +APi(1:2:2*Busnum,1:2:2*Busnum)=dPdTidTj;%%ǶԽ TT +APi(1:2:2*Busnum,2:2:2*Busnum)=dPdTidVj;%%ǶԽ TV +APi(2:2:2*Busnum,1:2:2*Busnum)=dPdVidTj;%%ǶԽ VT +APi(2:2:2*Busnum,2:2:2*Busnum)=dPdVidVj;%%ǶԽ VV +APi(1:2:2*Busnum,1:2:2*Busnum)=dPdTidTi;%%Խ +APi(1:2:2*Busnum,2:2:2*Busnum)=dPdTidVi;%%Խ +APi(2:2:2*Busnum,1:2:2*Busnum)=dPdVidTi;%%Խ +APi(2:2:2*Busnum,2:2:2*Busnum)=dPdVidVi;%%Խ +%% deltaQ/deltaThyta_deltaThyta ǶԽԪ +t1=-Volt'*Volt; +t2=real(GB).*(sin(AngleIJMat)-imag(GB).*cos(AngleIJMat)); +dQidTidTj=t1.*t2;%ȥԽԪˣҪ +t3=repmat(Init_Y',1,size(Init_Y,2)); +dQidTidTj=dQidTidTj.*t3(2:2:size(Init_Y,2),2:2:size(Init_Y,2));%yϵ +t2=real(GB).*(sin(AngleIJMat)+imag(GB).*cos(AngleIJMat)); +dQjdTidTj=t1.*t2; +t3=t3'; +dQjdTidTj=dQjdTidTj.*t3(2:2:size(Init_Y,2),2:2:size(Init_Y,2));%yϵ +dQdTidTj=dQidTidTj+dQjdTidTj;%շǶԽԪ +%% deltaQ/deltaThyta_deltaThyta ԽԪ +t1=dQidTidTj-diag(diag(dQidTidTj)); +t2=sum(t1,2); +t3=t2'.*Init_Y(2:2:size(Init_Y,2));%yϵ +dQidTidTi=diag(t3); +t1=-Volt'*Volt; +t2=real(GB).*(sin(AngleIJMat)-imag(GB).*cos(AngleIJMat)); +t3=t1.*t2; +t4=t3-diag(diag(t3)); +t5=sum(t4,1); +t6=t5.*Init_Y(2:2:size(Init_Y,2));%yϵ +dQjdTidTi=diag(t6); +dQdTidTi=dQjdTidTi+dQidTidTi;%նԽԪ +%% deltaQ/deltaThyta_deltaV ǶԽԪ +t1=-Volt; +t2=real(GB).*(cos(AngleIJMat)+imag(GB).*sin(AngleIJMat)); +t3=ones(Busnum,1)*t1.*t2; +t4=repmat(Init_Y',1,size(Init_Y,2)); +dQidTidVj=t3.*t4(2:2:size(Init_Y,2),2:2:size(Init_Y,2));%yϵ +t2=real(GB).*(cos(AngleIJMat)-imag(GB).*sin(AngleIJMat)); +t3=Volt'*ones(1,Busnum).*t2; +t4=t4'; +dQjdTidVj=t2.*t3.*t4(2:2:size(Init_Y,2),2:2:size(Init_Y,2));%yϵ +dQdTidVj=dQidTidVj+dQjdTidVj;%շǶԽԪ +%% deltaQ/deltaThyta_deltaV ԽԪ +t1=sum(dQidTidVj,2)-diag(dQidTidVj); +t2=t1'.*Init_Y(2:2:size(Init_Y,2));%yϵ +dQidTidVi=diag(t2); +t1=-Volt'*ones(1,Busnum).*(real(GB).*cos(AngleIJMat)+imag(GB).*sin(AngleIJMat)); +t2=t1-diag(diag(t1)); +t3=sum(t2,1); +t4=t3.*Init_Y(2:2:size(Init_Y,2));%yϵ +dQjdTidVi=diag(t4); +dQdTidVi=dQidTidVi+dQjdTidVi; +%% deltaQ/deltaV_deltaV ǶԽԪ +t1=real(GB).*(sin(AngleIJMat)-imag(GB).*cos(AngleIJMat)); +t2=repmat(Init_Y',1,size(Init_Y,2)); +t3=t1.*t2(2:2:size(Init_Y,2),2:2:size(Init_Y,2));%yϵ +dQidVidVj=t3; +t1=-real(GB).*sin(AngleIJMat)-imag(GB).*cos(AngleIJMat); +t2=t2'; +t3=t1.*t2(2:2:size(Init_Y,2),2:2:size(Init_Y,2));%yϵ +dQjdVidVj=t3; +dQdVidVj=dQidVidVj+dQjdVidVj; +%% deltaQ/deltaV_deltaV ԽԪ +t1=2*(real(GB).*(sin(AngleIJMat)-imag(GB).*cos(AngleIJMat))); +t2=diag(t1); +t3=t2'.*Init_Y(2:2:size(Init_Y,2));%yϵ +dQidVidVi=diag(t3); +t1=-real(GB).*sin(AngleIJMat)-imag(GB).*cos(AngleIJMat); +t2=diag(t1); +t3=t1-diag(t2); +t4=sum(t3,1); +t5=t4.*Init_Y(2:2:size(Init_Y,2));%yϵ +dQjdVidVi=diag(t5); +dQdVidVi=dQidVidVi+dQjdVidVi; +%% deltaQ/deltaV_deltaThyta ǶԽԪ +t1=real(GB).*sin(AngleIJMat)+imag(GB).*cos(AngleIJMat); +t2=-ones(Busnum,1)*Volt.*(t1); +t3=repmat(Init_Y',1,size(Init_Y,2)); +t4=t2'.*t3(2:2:size(Init_Y,2),2:2:size(Init_Y,2));%yϵ +dQidVidTj=t4; +t1=real(GB).*sin(AngleIJMat)-imag(GB).*cos(AngleIJMat); +t2=-ones(Busnum,1)*Volt.*(t1); +t3=t3'; +t4=t2.*t3(2:2:size(Init_Y,2),2:2:size(Init_Y,2));%yϵ +dQjdVidTj=t4; +dQdVidTj=dQidVidTj+dQjdVidTj; +%% deltaQ/deltaV_deltaThyta ԽԪ +dQdVidTi=dQdTidVi; +%% AQi +AQi=zeros(2*Busnum,2*Busnum); +AQi(1:2:2*Busnum,1:2:2*Busnum)=dQdTidTj;%%ǶԽ TT +AQi(1:2:2*Busnum,2:2:2*Busnum)=dQdTidVj;%%ǶԽ TV +AQi(2:2:2*Busnum,1:2:2*Busnum)=dQdVidTj;%%ǶԽ VT +AQi(2:2:2*Busnum,2:2:2*Busnum)=dQdVidVj;%%ǶԽ VV +AQi(1:2:2*Busnum,1:2:2*Busnum)=dQdTidTi;%%Խ +AQi(1:2:2*Busnum,2:2:2*Busnum)=dQdTidVi;%%Խ +AQi(2:2:2*Busnum,1:2:2*Busnum)=dQdVidTi;%%Խ +AQi(2:2:2*Busnum,2:2:2*Busnum)=dQdVidVi;%%Խ +%% ddf +t=[zeros(2*14-2*5,2*14); + zeros(2*5,2*14-2*5),AQi+APi; +]; +ddf=t; +end \ No newline at end of file diff --git a/func_ddg.m b/func_ddg.m new file mode 100644 index 0000000..4f33610 --- /dev/null +++ b/func_ddg.m @@ -0,0 +1,4 @@ +function ddg=func_ddg(PGi,PVi,Busnum,RestraintCount,Loadi,PD,QD) + +ddg=0; +end \ No newline at end of file diff --git a/func_ddh.m b/func_ddh.m new file mode 100644 index 0000000..3c07a62 --- /dev/null +++ b/func_ddh.m @@ -0,0 +1,56 @@ +function ddh=func_ddh(Volt,Init_Y,Busnum,Y,UAngel,r,c,YAngle,Loadi,ContrlCount) +%ѭд +%ContrlCount=size(PVi,1)+size(PGi,1)+size(Loadi,1)+Busnum*2; + +%AngleIJ=AngleIJMat-angle(GB); +mat_AngleIJ=sparse(r,c,UAngel(r)-UAngel(c)-YAngle,Busnum*3,Busnum*3); +mat_INV_AngleIJ=mat_AngleIJ'; +yP=Init_Y(1:size(Init_Y,2)/2);%ʱ 20111227 +yQ=Init_Y(size(Init_Y,2)/2+1:size(Init_Y,2));%ʱ 20111227 +t1=-sparse(1:Busnum*3,1:Busnum*3,Y.*cos(mat_INV_AngleIJ)*sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3)*yP',Busnum*3,Busnum*3); +t2=sparse(1:Busnum*3,1:Busnum*3,sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3)*yP',Busnum*3,Busnum*3)*Y.*cos(mat_AngleIJ); +t3=(t1+t2)*sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3); +t4=-(sparse(1:Busnum*3,1:Busnum*3,Y.*cos(mat_AngleIJ)*Volt,Busnum*3,Busnum*3) -sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3)*Y.*cos(mat_INV_AngleIJ) )*sparse(1:Busnum*3,1:Busnum*3,sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3)*yP',Busnum*3,Busnum*3); +ddPdTdT=t3+t4;%ok1 +t1=(-sparse(1:Busnum*3,1:Busnum*3,Y.*sin(mat_AngleIJ)*Volt,Busnum*3,Busnum*3)+sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3)*Y.*sin(mat_INV_AngleIJ) )*sparse(1:Busnum*3,1:Busnum*3,yP,Busnum*3,Busnum*3); +t2= -sparse(1:Busnum*3,1:Busnum*3, sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3)*yP' ,Busnum*3,Busnum*3)*Y.*sin(mat_AngleIJ)+sparse(1:Busnum*3,1:Busnum*3,Y.*sin(mat_INV_AngleIJ)*sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3)*yP',Busnum*3,Busnum*3); +ddPdVdT=t1+t2;%ok1 +t1=diag( Y.*sin(mat_INV_AngleIJ)*diag(Volt)*yP'); +t2=diag(yP)*Y.*sin(mat_AngleIJ)*diag(Volt); +t3=-diag(yP)*diag(Y.*sin(mat_AngleIJ)*Volt); +t4=-Y.*sin(mat_INV_AngleIJ)*diag( diag(Volt)*yP' ); +ddPdTdV=t1+t2+t3+t4;%ҵIJһ +t1=Y.*cos(mat_INV_AngleIJ)*diag(yP); +t2=diag(yP)*Y.*cos(mat_AngleIJ); +ddPdVdV=t1+t2; +t1=-diag(Y.*sin(mat_AngleIJ)*Volt); +t2=diag(Volt)*Y.*sin(mat_INV_AngleIJ); +t3=(t1+t2)*diag( diag(Volt)*yQ' ); +t4=-diag( diag(Volt)*yQ' )*Y.*sin(mat_AngleIJ); + +t5=diag(Y.*sin(mat_INV_AngleIJ)*diag(Volt)*yQ'); +t6=-(t4+t5)*diag(Volt); +ddQdTdT=t3+t6;%ok1 +t1=(diag(Y.*cos(mat_AngleIJ)*Volt)-diag(Volt)*Y.*cos(mat_INV_AngleIJ) )*diag(yQ); +t2=+diag( diag(Volt)*yQ' )*Y.*cos(mat_AngleIJ)-diag(Y.*cos(mat_INV_AngleIJ)*diag(Volt)*yQ'); +ddQdVdT=t1+t2; +t1=Y.*cos(mat_INV_AngleIJ)*diag(diag(Volt)*yQ'); +t2=diag(yQ)*diag(Y.*cos(mat_AngleIJ)*Volt); +t3=-diag(Y.*cos(mat_INV_AngleIJ)*diag(Volt)*yQ'); +t4=-diag(yQ)*Y.*cos(mat_AngleIJ)*diag(Volt); +ddQdTdV=t1+t2+t3+t4; +t1=Y.*sin(mat_INV_AngleIJ)*diag(yQ); +t2=diag(yQ)*Y.*sin(mat_AngleIJ); +ddQdVdV=t1+t2; +% t=[ddPdVdV+ddQdVdV,ddPdTdV+ddQdTdV; +% ddPdVdT+ddQdVdT,ddPdTdT+ddQdTdT; +% ]; +t=[ddPdVdV+ddQdVdV,ddPdVdT+ddQdVdT; + ddPdTdV+ddQdTdV,ddPdTdT+ddQdTdT; +]; +sizeLoadi=size(Loadi,1)*2; +ddh=[ + sparse(sizeLoadi,ContrlCount); + sparse(2*Busnum*3,sizeLoadi),-t; + ]; +end \ No newline at end of file diff --git a/func_deltF.m b/func_deltF.m new file mode 100644 index 0000000..840fc36 --- /dev/null +++ b/func_deltF.m @@ -0,0 +1,11 @@ +function deltF=func_deltF(wVolt,wPD,wQD,PD0,PD,QD,QD0,Volt,mVolt,Busnum,Loadi) +t3=2*wPD.*(PD-PD0); +t4=2*wQD.*(QD-QD0); +t5=2*wVolt.*(Volt-mVolt); +deltF=[sparse(t3); + sparse(t4); + sparse(t5); + sparse(1*Busnum*3,1); +]; + +end \ No newline at end of file diff --git a/func_deltG.m b/func_deltG.m new file mode 100644 index 0000000..23430f3 --- /dev/null +++ b/func_deltG.m @@ -0,0 +1,32 @@ +function deltG=func_deltG(Busnum,Loadi) +%% + +sizeLoadi=size(Loadi,1); +%% + + +%% +dg3_dPD=sparse(1:size(Loadi,1),1:size(Loadi,1),ones(size(Loadi,1),1),size(Loadi,1),size(Loadi,1)); +dg4_dPD=sparse(size(Loadi,1),length(Loadi)); +dg5_dPD=sparse(size(Loadi,1),Busnum*3); +%% +dg3_dQD=sparse(length(Loadi),length(Loadi)); +dg4_dQD=sparse(1:size(Loadi,1),1:size(Loadi,1),ones(size(Loadi,1),1),size(Loadi,1),size(Loadi,1)); +dg5_dQD=sparse(size(Loadi,1),Busnum*3); + +%% +dg3_dx=sparse(2*Busnum*3,sizeLoadi); +dg4_dx=sparse(2*Busnum*3,length(Loadi)); +dg5_dx=[sparse(1:Busnum*3,1:Busnum*3,ones(Busnum*3,1),Busnum*3,Busnum*3); + sparse(Busnum*3,Busnum*3); + ]; +%% + +deltG=[dg3_dPD,dg4_dPD,dg5_dPD; + dg3_dQD,dg4_dQD,dg5_dQD; + dg3_dx,dg4_dx,dg5_dx; +]; +deltG=[dg3_dPD,dg4_dPD; + dg3_dQD,dg4_dQD; + dg3_dx,dg4_dx;]; +end \ No newline at end of file diff --git a/func_deltH.m b/func_deltH.m new file mode 100644 index 0000000..6d5e4d2 --- /dev/null +++ b/func_deltH.m @@ -0,0 +1,7 @@ +function deltH=func_deltH(Busnum,Volt,Y,UAngel,r,c,Angle,Loadi) +dH_dPD=[sparse(1:size(Loadi,1),Loadi,-ones(size(Loadi,1),1),size(Loadi,1),Busnum*3) sparse(size(Loadi,1),Busnum*3)]; +dH_dQD=[sparse(size(Loadi,1),Busnum*3) sparse(1:size(Loadi,1),Loadi,-ones(size(Loadi,1),1),size(Loadi,1),Busnum*3)]; +dH_dx = jacobian_M(Busnum,Volt,Y,Angle,UAngel,r,c); %γſ˱Ⱦ +deltH=[dH_dPD;dH_dQD;dH_dx']; + +end \ No newline at end of file diff --git a/func_deltdeltF.m b/func_deltdeltF.m new file mode 100644 index 0000000..004023e --- /dev/null +++ b/func_deltdeltF.m @@ -0,0 +1,14 @@ +function deltdeltF=func_deltdeltF(wVolt,wPD,wQD,ContrlCount,Loadi) + +%ContrlCount=size(PVi,1)+size(PGi,1)+size(wD,1)+Busnum*2; %P,Q,Volt thetaЩƱ +C=[wPD' wQD']; +CVolt=wVolt; +sizeC=length(C); +diagC=sparse(1:sizeC,1:sizeC,C,sizeC,sizeC); +deltdeltF=[ + diagC*2,sparse(sizeC,ContrlCount-sizeC); + %sparse(ContrlCount-sizeC,ContrlCount); + diag(CVolt)*2,sparse(length(CVolt),ContrlCount-length(CVolt)); + sparse(length(CVolt),ContrlCount); + ]; +end \ No newline at end of file diff --git a/jacobian.m b/jacobian.m new file mode 100644 index 0000000..596b585 --- /dev/null +++ b/jacobian.m @@ -0,0 +1,40 @@ +function [Jacob,PQ,U,Uangle]=jacobian(Busnum,Balance,PVi,PVu,U,Uangle,Y,Angle,P0,Q0,r,c) +%************************************************************************** +% : ӺγſɱȾJacobian +% ߣ +% ʱ䣺2010.12 +%************************************************************************** +%% ֱſ˱ȾH,L,N,Jй޹P,Q +AngleIJ = Uangle(r) - Uangle(c)- Angle'; +U(PVi) = PVu; +U(Balance)=1; +temp1= -sparse(1:Busnum,1:Busnum,U,Busnum,Busnum)*Y*sparse(1:Busnum,1:Busnum,U,Busnum,Busnum); % ſ˱Ⱦõм +temp2 = sum(temp1.*sparse(r,c,sin(AngleIJ)),2); +temp3 = sum(temp1.*sparse(r,c,cos(AngleIJ)),2); +temp4=sparse(1:Busnum,1:Busnum,temp2,Busnum,Busnum); +temp5=sparse(1:Busnum,1:Busnum,temp3,Busnum,Busnum); +H = temp1.*sparse(r,c,sin(AngleIJ))-temp4; +L = temp1.*sparse(r,c,sin(AngleIJ))+temp4; +N = temp1.*sparse(r,c,cos(AngleIJ))+temp5; +J = -temp1.*sparse(r,c,cos(AngleIJ))+temp5; + +Q = Q0+temp2'; %йP +P = P0+temp3'; %޹Q +%% ƽڵpvڵ +H(:,Balance) = 0; +H(Balance,:) = 0; +%H(Balance,Balance) = 100; % ƽڵӦĶԽԪһ +H=H+sparse(Balance,Balance,ones(1,length(Balance)),Busnum,Busnum); +L(:,PVi) = 0; +L(PVi,:) = 0; +L = L+sparse(PVi,PVi,ones(1,length(PVi)),Busnum,Busnum); % PVڵӦĶԽԪΪ1 +J(:,Balance) = 0; +J(PVi,:) = 0; +N(:,PVi) = 0; +N(Balance,:) = 0; +Q(PVi) = 0; % pvڵ޹ƽ +P(Balance) = 0; % ƽڵйʲƽ +%% ϳPQſɱȾ +PQ = cat(2,P,Q); % γɹʲƽ +Jacob = cat(1,cat(2,H,N),cat(2,J,L)); % γJacobian +end \ No newline at end of file diff --git a/jacobian_M.m b/jacobian_M.m new file mode 100644 index 0000000..fb999de --- /dev/null +++ b/jacobian_M.m @@ -0,0 +1,20 @@ +function [Jacob]=jacobian_M3(Busnum,Volt,Y,Angle,UAngel,r,c) +%************************************************************************** +% : ӺγſɱȾJacobian +% ߣ +% ʱ䣺2010.12 +%************************************************************************** +%%ͼ6¥д +%% ֱſ˱ȾH,L,N,Jй޹P,Q +AngleIJ=UAngel(r)-UAngel(c)-Angle; +mat_AngleIJ=sparse(r,c,AngleIJ,Busnum*3,Busnum*3); +mat_IvAngleIJ=mat_AngleIJ'; +H=sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3)*Y.*sin(mat_IvAngleIJ)*sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3)-sparse(1:Busnum*3,1:Busnum*3,Y.*sin(mat_AngleIJ)*Volt,Busnum*3,Busnum*3)*sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3); +N=-sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3)*Y.*cos(mat_IvAngleIJ)*sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3)+sparse(1:Busnum*3,1:Busnum*3,Y.*cos(mat_AngleIJ)*Volt,Busnum*3,Busnum*3)*sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3); +J=sparse(1:Busnum*3,1:Busnum*3,Y.*cos(mat_AngleIJ)*Volt,Busnum*3,Busnum*3)+Y.*cos(mat_IvAngleIJ)*sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3); +L=sparse(1:Busnum*3,1:Busnum*3,Y.*sin(mat_AngleIJ)*Volt,Busnum*3,Busnum*3)+Y.*sin(mat_IvAngleIJ)*sparse(1:Busnum*3,1:Busnum*3,Volt,Busnum*3,Busnum*3); +t1=[J,L; + H,N; +]'; +Jacob=-t1; +end \ No newline at end of file diff --git a/lineWithConfig.m b/lineWithConfig.m new file mode 100644 index 0000000..16c2a87 --- /dev/null +++ b/lineWithConfig.m @@ -0,0 +1,106 @@ +function [ fs30,fs31,fs32,phaseABCY ] = lineWithConfig(setIJ,nodeNum,lineZ,lines,config,baseHighU ) +phase3Entry=find(lines(:,1)==config); +phase3Line=lines(phase3Entry,:); +phase3Line(:,4)=phase3Line(:,4)/1000; +entry=find(lineZ(:,1)==config); +if isempty(entry) + fs30=0; + fs31=0; + fs32=0; + phaseABCY=0; + fprintf('warning: line type %d not found.\n',config) + return; +end +phase3R=lineZ(entry+1:entry+3,:); +phase3X=lineZ(entry+5:entry+7,:); +phase3B2=lineZ(entry+9:entry+11,:); +phase3Z=phase3R+1j*phase3X; +phase3Y=inv(phase3Z); +%ֵ +baseS=1; +baseY=baseS/(baseHighU^2); +phase3Y=phase3Y/baseY; +diagEle=diag(phase3Y);%ԽԪ +diagEle=sum(diagEle)/3; +offDiagEle=phase3Y-diag(diag(phase3Y)); +offDiagEle=sum(sum(offDiagEle))/6; +%öԽǺͷǶԽ +phase3Y(1,1)=diagEle; +phase3Y(1,2)=offDiagEle; +phase3Y(2,1)=offDiagEle; +phase3Y(3,1)=offDiagEle; +phase3Y(1,3)=offDiagEle; +phase3Y(1,3)=phase3Y(1,2); +phase3Y(3,1)=phase3Y(2,1); +phase3Y(3,2)=phase3Y(3,1); +phase3Y(2,3)=phase3Y(1,3); +phase3Y(2,2)=phase3Y(1,1); +phase3Y(3,3)=phase3Y(1,1); +%Fortescue任 +%ֱд +a=exp(1j*2*pi/3); +Tp2f=1/3*[1 1 1; + 1 a a^2; + 1 a^2 a]; +Tf2p=inv(Tp2f); +% +% Tf2p=[1 1 1; +% 1 a^2 a; +% 1 a a^2]; +% fs3Y=Tp2f*phase3Y*Tf2p; +% Zl=phase3Y(1,1); +% Zm=phase3Y(2,1); +Zl=diagEle; +Zm=offDiagEle; +fs3Y=3/3*diag([Zl+2*Zm,Zl-Zm,Zl-Zm]);%ﲻ1/3IEEEϵĹʽˡ +fs3Y(abs(fs3Y)<1e-5)=0; +fs3Y=sparse(fs3Y); +lineNodeI=zeros(length(phase3Line(:,2) ) ,1); +lineNodeJ=zeros(length(phase3Line(:,2) ) ,1); +for I=1:length(lineNodeI) + lineNodeI(I)=nodeNum(setIJ==phase3Line(I,2) ); + lineNodeJ(I)=nodeNum(setIJ==phase3Line(I,3) ); +end +busNum=length(nodeNum); +%ർɾ +phaseABCY=sparse(busNum*3,busNum*3); +for I=1:length(phase3Line(:,4)) + [phase3Yr,phase3Yc,phase3Yv]=find(phase3Y./phase3Line(I,4)); + offsetI=phase3Yr+3*(lineNodeI(I)-1); + offsetJ=phase3Yc+3*(lineNodeJ(I)-1); + phaseABCY=phaseABCY+sparse(offsetI,offsetJ,-phase3Yv,busNum*3,busNum*3); + phaseABCY=phaseABCY+sparse(offsetJ,offsetI,-phase3Yv,busNum*3,busNum*3); + %1 + diagY=sum(phaseABCY(:,1:3:end),2); + offsetI=1:busNum*3; + offsetJ=zeros(busNum*3,1); + for J=0:busNum*3-1 + offsetJ(J+1)=floor(J/3)*3+1; + end + phaseABCY=phaseABCY-sparse(offsetI,offsetJ,diagY,busNum*3,busNum*3); + %2 + diagY=sum(phaseABCY(:,2:3:end),2); + offsetJ=zeros(busNum*3,1); + for J=0:busNum*3-1 + offsetJ(J+1)=floor(J/3)*3+2; + end + phaseABCY=phaseABCY-sparse(offsetI,offsetJ,diagY,busNum*3,busNum*3); + %3 + diagY=sum(phaseABCY(:,3:3:end),2); + offsetJ=zeros(busNum*3,1); + for J=0:busNum*3-1 + offsetJ(J+1)=floor(J/3)*3+3; + end + phaseABCY=phaseABCY-sparse(offsetI,offsetJ,diagY,busNum*3,busNum*3); +end +fs30=sparse(lineNodeI,lineNodeJ,-fs3Y(1,1)./phase3Line(:,4),busNum,busNum); +fs30=fs30+sparse(lineNodeJ,lineNodeI,-fs3Y(1,1)./phase3Line(:,4),busNum,busNum); +fs30=fs30-sparse(1:busNum,1:busNum,sum(fs30,2)); +fs31=sparse(lineNodeI,lineNodeJ,-fs3Y(2,2)./phase3Line(:,4),busNum,busNum); +fs31=fs31+sparse(lineNodeJ,lineNodeI,-fs3Y(2,2)./phase3Line(:,4),busNum,busNum); +fs31=fs31-sparse(1:busNum,1:busNum,sum(fs31,2)); +fs32=sparse(lineNodeI,lineNodeJ,-fs3Y(3,3)./phase3Line(:,4),busNum,busNum); +fs32=fs32+sparse(lineNodeJ,lineNodeI,-fs3Y(3,3)./phase3Line(:,4),busNum,busNum); +fs32=fs32-sparse(1:busNum,1:busNum,sum(fs32,2)); +end + diff --git a/numberNode.m b/numberNode.m new file mode 100644 index 0000000..2cc497c --- /dev/null +++ b/numberNode.m @@ -0,0 +1,8 @@ +function [ setIJ,nodeNum ] = numberNode( lines ) +%% ڵ +setI=lines(:,2); +setJ=lines(:,3); +setIJ=union(setI,setJ); +nodeNum=(1:length(setIJ))'; +end + diff --git a/readLineZ.m b/readLineZ.m new file mode 100644 index 0000000..5f71d9c --- /dev/null +++ b/readLineZ.m @@ -0,0 +1,4 @@ +function [ data ] = readLineZ( lineParaFile ) +data=dlmread(lineParaFile); +end + diff --git a/test.m b/test.m new file mode 100644 index 0000000..e73dac0 --- /dev/null +++ b/test.m @@ -0,0 +1,42 @@ +close all +clear +clc +AME_VoltSum=0; +AME_VAngleSum=0; +JMeasurementSum=0; +AME_PDSum=0; +AME_QDSum=0; +AME_mVoltSum=0; +AME_mPDSum=0; +AME_mQDSum=0; +totalTimeSum=1000; +N=200; +loopN=1; +while 1 + [JMeasurement,AME_Volt,AME_VAngle,AME_PD,AME_QD,AME_mVolt,AME_mPD,AME_mQD,isConverged,totalTime]=OPF(); + if isConverged==0 + %I=I-1; + continue + end + loopN=loopN+1; + if loopN>N + break; + end +% [JMeasurement]=run(); + AME_VoltSum=AME_VoltSum+AME_Volt; + AME_VAngleSum=AME_VAngleSum+AME_VAngle; + AME_PDSum=AME_PDSum+AME_PD; + AME_QDSum=AME_QDSum+AME_QD; + AME_mVoltSum=AME_mVoltSum+AME_mVolt; + AME_mPDSum=AME_mPDSum+AME_mPD; + AME_mQDSum=AME_mQDSum+AME_mQD; + totalTimeSum=min([totalTimeSum,totalTime]); +end +JMeasurementSum=JMeasurementSum/N; +AME_VoltSum=AME_VoltSum/N; +AME_VAngleSum=AME_VAngleSum/N; +AME_PDSum=AME_PDSum/N; +AME_QDSum=AME_QDSum/N; +AME_mVoltSum=AME_mVoltSum/N; +AME_mPDSum=AME_mPDSum/N; +AME_mQDSum=AME_mQDSum/N;