From 9385ef2269cf891ec43aabdda53ce071be3f57bd Mon Sep 17 00:00:00 2001 From: "dmy@lab" Date: Sat, 16 May 2015 15:49:16 +0800 Subject: [PATCH] =?UTF-8?q?=E5=A2=9E=E5=8A=A0=E4=BF=AE=E6=94=B9=E7=9A=8469?= =?UTF-8?q?=E8=8A=82=E7=82=B9=E7=AE=97=E4=BE=8B?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: dmy@lab --- func_deltF.m | 6 + func_deltdeltF.m | 5 +- modified-feeder69/data.txt | 154 +++++++++++++++++++++++++ modified-feeder69/feeder69.txt | 172 ++++++++++++++++++++++++++++ modified-feeder69/feeder69.xlsx | Bin 0 -> 25087 bytes modified-feeder69/feeder69ieee.txt | 166 +++++++++++++++++++++++++++ modified-feeder69/lineParameter.txt | 40 +++++++ modified-feeder69/~$feeder69.xlsx | Bin 0 -> 165 bytes run.m | 44 ++++--- 9 files changed, 567 insertions(+), 20 deletions(-) create mode 100644 modified-feeder69/data.txt create mode 100644 modified-feeder69/feeder69.txt create mode 100644 modified-feeder69/feeder69.xlsx create mode 100644 modified-feeder69/feeder69ieee.txt create mode 100644 modified-feeder69/lineParameter.txt create mode 100644 modified-feeder69/~$feeder69.xlsx diff --git a/func_deltF.m b/func_deltF.m index 5def6d4..e94101f 100644 --- a/func_deltF.m +++ b/func_deltF.m @@ -6,6 +6,12 @@ function deltF=func_deltF(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,busNum % t5=sum(t4,2); % PBal=t5(Balance); % PPG=([PQ(1),PBal])';%暂时用土办法处理一下 + +% wV1r=1000*ones(length(wV1r),1); +% wV1i=1000*ones(length(wV1i),1); +% wI1r=3300*ones(length(wI1r),1); +% wI1i=3300*ones(length(wI1i),1); + %% deltF=[ %zeros(busNum*2,1); diff --git a/func_deltdeltF.m b/func_deltdeltF.m index 8116328..ef0f2a7 100644 --- a/func_deltdeltF.m +++ b/func_deltdeltF.m @@ -1,5 +1,8 @@ function deltdeltF=func_deltdeltF(busNum,fsY11,Loadi,wV1r,wV1i,wI1r,wI1i,V1measurement,V1r,V1i) - +% wV1r=1000*ones(length(wV1r),1); +% wV1i=1000*ones(length(wV1i),1); +% wI1r=3300*ones(length(wI1r),1); +% wI1i=3300*ones(length(wI1i),1); deltdeltF=[ %zeros(busNum*2,busNum*2+length(Loadi)*2); sparse(Loadi,Loadi,-4*( -2*V1r(Loadi).*V1r(Loadi)+ abs(V1measurement).^2-(V1r(Loadi).^2+ V1i(Loadi).^2) )./[wV1r.*wV1r],busNum,busNum),sparse(Loadi,Loadi+busNum,-4*V1r(Loadi).*( -2*V1i(Loadi))./[wV1r.*wV1r],busNum,busNum+length(Loadi)*2); diff --git a/modified-feeder69/data.txt b/modified-feeder69/data.txt new file mode 100644 index 0000000..72a5a40 --- /dev/null +++ b/modified-feeder69/data.txt @@ -0,0 +1,154 @@ +1000 12.66 +0 +222 1000 1 1 +222 1 2 1 +222 2 3 3 +222 3 4 50 +222 4 5 730 +222 5 6 760 +222 6 7 185 +222 7 8 99 +222 8 9 1638 +222 9 10 374 +222 10 11 1442 +222 11 12 2060 +222 12 13 2088 +222 13 14 2116 +222 14 15 393 +222 15 16 748 +222 16 17 10 +222 17 18 655 +222 18 19 421 +222 19 20 683 +222 20 21 28 +222 21 22 318 +222 22 23 692 +222 23 24 1497 +222 24 25 617 +222 25 26 346 +222 2 27 9 +222 27 28 128 +222 28 29 795 +222 29 30 140 +222 30 31 702 +222 31 32 1678 +222 32 33 3416 +222 33 34 2948 +222 2 277 9 +222 277 288 128 +222 288 65 210 +222 65 66 60 +222 66 67 4 +222 67 68 1456 +222 68 69 620 +222 69 70 82 +222 70 88 18 +222 88 89 217 +222 89 90 2 +222 3 35 7 +222 35 36 172 +222 36 37 580 +222 37 38 164 +222 7 40 185 +222 40 41 663 +222 8 42 384 +222 42 43 406 +222 43 44 567 +222 44 45 562 +222 45 46 3180 +222 46 47 1567 +222 47 48 608 +222 48 49 772 +222 49 50 1015 +222 50 51 194 +222 51 52 290 +222 52 53 1421 +222 53 54 2082 +222 10 55 402 +222 55 56 9 +222 11 57 1478 +222 57 58 9 +0 +692 0 0 0 +0 +5 36.2 34.9 36.1 34.8 36.3 35 +6 51.5 48 51.4 47.9 51.6 48.1 +7 25 18 24.9 17.9 25.1 18.1 +8 48 45.3 47.9 45.2 48.1 45.4 +9 36.9 33.9 36.8 33.8 37 34 +10 48.3 34.7 48.2 34.6 48.4 34.8 +11 48.3 34.7 48.2 34.6 48.4 34.8 +12 9.33 5.17 9.23 5.07 9.43 5.27 +13 9.33 5.17 9.23 5.07 9.43 5.27 +15 15.2 10 15.1 99 15.3 10.1 +16 58 49.7 57.9 49.6 58.1 49.8 +17 20 11.7 19.9 11.6 20.1 11.8 +19 37.7 36.3 37.6 36.2 37.8 36.4 +20 38 2.7 37.9 26.9 38.1 27.1 +21 46.4 45.8 46.3 45.7 46.5 45.9 +23 47.3 44.7 47.2 44.6 47.4 44.8 +25 32.2 30.9 32.1 30.8 32.3 31 +26 32.2 30.9 32.1 30.8 32.3 31 +27 8.67 6.2 8.57 6.1 8.77 6.3 +28 36.2 33.7 36.1 33.6 36.3 33.8 +277 46.7 44.2 46.6 44.1 46.8 44.3 +288 8.67 6.18 8.57 6.08 8.77 6.28 +32 4.67 3.33 4.57 3.23 4.77 3.43 +33 44.5 42.7 44.4 42.6 44.6 42.8 +34 32.9 33.2 32.8 32.1 33 32.3 +36 64.3 56.8 64.2 56.7 64.4 56.9 +37 128 91.5 12.8 91.4 12.8 91.6 +38 128 91.5 12.8 91.4 12.8 91.6 +40 43.5 39.4 43.4 39.3 43.6 39.5 +41 45.5 42.2 45.4 42.1 45.6 42.3 +42 31.5 31.2 31.4 31.1 31.6 31.3 +43 46.8 44.3 46.7 44.2 46.9 44.4 +44 8 5.73 7.9 5.63 8.1 5.83 +48 33.3 24 33.2 23.9 33.4 24.1 +50 415 296 415 296 415 296 +51 78.7 75.7 78.6 75.6 78.8 75.8 +53 75.7 54 75.6 53.9 75.8 54.1 +54 87.7 82 87.6 81.9 87.8 82.1 +55 74 72.3 73.9 72.2 74.1 72.4 +56 36 34.3 35.9 34.2 36.1 34.4 +57 39.3 36.7 39.2 36.6 39.4 36.8 +58 80.3 77.7 80.2 77.6 80.4 77.8 +70 43.3 42.8 43.2 42.7 43.4 42.9 +89 13.1 8.77 13 86.7 13.2 8.87 +90 51.1 46.8 51 46.7 51.2 46.9 +66 35.5 33.2 35.4 33.1 35.6 33.3 +67 35.5 33.2 35.4 33.1 35.6 33.3 +68 7.07 3.67 6.97 3.57 7.17 3.77 +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/modified-feeder69/feeder69.txt b/modified-feeder69/feeder69.txt new file mode 100644 index 0000000..7a2d917 --- /dev/null +++ b/modified-feeder69/feeder69.txt @@ -0,0 +1,172 @@ +69 68 1 1.00E-06 9 +0 +1 1 1 +0 +1 1 2 3.12E-06 1.50E-05 0 0 +2 2 3 3.12E-06 1.50E-05 0 0 +3 3 4 9.36E-06 4.49E-05 0 0 +4 4 5 0.001566052 0.000366868 0 0 +5 5 6 0.002283567 0.002325993 0 0 +6 6 7 0.002377155 0.002422078 0 0 +7 7 8 0.000575259 5.86E-05 0 0 +8 8 9 0.000307595 0.00031321 0 0 +9 9 10 0.005109948 0.003377932 0 0 +10 10 11 0.001167988 0.00077242 0 0 +11 11 12 0.004438605 0.002933697 0 0 +12 12 13 0.00642643 0.004242692 0 0 +13 13 14 0.00651378 0.004305084 0 0 +14 14 15 0.00660113 0.004362486 0 0 +15 15 16 0.001226637 0.000811103 0 0 +16 16 17 0.002354694 0.001544839 0 0 +17 17 18 2.93E-05 2.00E-05 0 0 +18 18 19 0.002043979 0.001351422 0 0 +19 19 20 0.001313987 0.000861017 0 0 +20 20 21 0.002131329 0.001408823 0 0 +21 21 22 8.73E-05 5.74E-05 0 0 +22 22 23 0.000992665 0.000656369 0 0 +23 23 24 0.002160653 0.001428789 0 0 +24 24 25 0.004671953 0.00308843 0 0 +25 25 26 0.001927305 0.001274055 0 0 +26 26 27 0.001080639 0.000713771 0 0 +27 3 28 2.75E-05 0.000134768 0 0 +28 28 29 0.000399312 0.001952886 0 0 +29 29 30 0.002481975 0.001640924 0 0 +30 30 31 0.000437996 0.000289501 0 0 +31 31 32 0.002189978 0.001447507 0 0 +32 32 33 0.005234733 0.003513947 0 0 +33 33 34 0.010656644 0.007045364 0 0 +34 34 35 0.009196659 0.006080776 0 0 +35 3 36 2.75E-05 0.000134768 0 0 +36 36 37 0.000399312 0.001952886 0 0 +37 37 38 0.000656993 0.001534856 0 0 +38 38 39 0.000189673 0.000442987 0 0 +39 39 40 1.12E-05 2.62E-05 0 0 +40 40 41 0.004544048 0.010617961 0 0 +41 41 42 0.001934168 0.004520963 0 0 +42 42 43 0.000255809 0.000596473 0 0 +43 43 44 5.74E-05 0.000144751 0 0 +44 44 45 0.000679455 0.001713299 0 0 +45 45 46 5.62E-06 1.50E-05 0 0 +46 4 47 2.12E-05 0.000104819 0 0 +47 47 48 0.00053096 0.002599273 0 0 +48 48 49 0.001808135 0.008848508 0 0 +49 49 50 0.000512867 0.002509428 0 0 +50 8 51 0.000579003 0.000590233 0 0 +51 51 52 0.002070808 0.001390106 0 0 +52 9 53 0.00108563 0.001105596 0 0 +53 53 54 0.001266568 0.001290277 0 0 +54 54 55 0.001773196 0.00180564 0 0 +55 55 56 0.001755102 0.00178817 0 0 +56 56 57 0.009920412 0.006659779 0 0 +57 57 58 0.004952095 0.003281847 0 0 +58 58 59 0.001897981 0.001255338 0 0 +59 59 60 0.002408976 0.001462481 0 0 +60 60 61 0.003166421 0.003225694 0 0 +61 61 62 0.000607703 0.000618934 0 0 +62 62 63 0.000904692 0.000920914 0 0 +63 63 64 0.004432989 0.004515971 0 0 +64 64 65 0.006495062 0.006616104 0 0 +65 11 66 0.001255338 0.000762437 0 0 +66 66 67 2.93E-05 1.75E-05 0 0 +67 12 68 0.004613304 0.003049747 0 0 +68 68 69 2.93E-05 2.00E-05 0 0 +0 +0 +0 +0 +0 +0 A A B B C C +1 -1 1 0 0 0 0 0 0 +2 -1 2 0 0 0 0 0 0 +3 -1 3 0 0 0 0 0 0 +4 -1 4 0 0 0 0 0 0 +5 -1 5 0 0 0.0362 0.0349 0.0361 0.0348 0.0363 0.035 0 0 0 0 +6 -1 6 0 0 0.0515 0.048 0.0514 0.0479 0.0516 0.0481 0.0026 0.0022 0 0 +7 -1 7 0 0 0.025 0.018 0.0249 0.0179 0.0251 0.0181 0.0404 0.03 0 0 +8 -1 8 0 0 0.048 0.0453 0.0479 0.0452 0.0481 0.0454 0.075 0.054 0 0 +9 -1 9 0 0 0.0369 0.0339 0.0368 0.0338 0.037 0.034 0.03 0.022 0 0 +10 -1 10 0 0 0.0483 0.0347 0.0482 0.0346 0.0484 0.0348 0.028 0.019 0 0 +11 -1 11 0 0 0.0483 0.0347 0.0482 0.0346 0.0484 0.0348 0.145 0.104 0 0 +12 -1 12 0 0 0.00933 0.00517 0.00923 0.00507 0.00943 0.00527 0.145 0.104 0 0 +13 -1 13 0 0 0.00933 0.00517 0.00923 0.00507 0.00943 0.00527 0.008 0.005 0 0 +14 -1 14 0 0 0.008 0.0055 0 0 +15 -1 15 0 0 0.0152 0.01 0.0151 0.099 0.0153 0.0101 0 0 0 0 +16 -1 16 0 0 0.058 0.0497 0.0579 0.0496 0.0581 0.0498 0.0455 0.03 0 0 +17 -1 17 0 0 0.02 0.0117 0.0199 0.0116 0.0201 0.0118 0.06 0.035 0 0 +18 -1 18 0 0 0.06 0.035 0 0 +19 -1 19 0 0 0.0377 0.0363 0.0376 0.0362 0.0378 0.0364 0 0 0 0 +20 -1 20 0 0 0.0464 0.0458 0.0463 0.0457 0.0465 0.0459 0.001 0.0006 0 0 +21 -1 21 0 0 0.0464 0.0458 0.0463 0.0457 0.0465 0.0459 0.114 0.081 0 0 +22 -1 22 0 0 0.005 0.0035 0 0 +23 -1 23 0 0 0.0473 0.0447 0.0472 0.0446 0.0474 0.0448 0 0 0 0 +24 -1 24 0 0 0.028 0.02 0 0 +25 -1 25 0 0 0.0322 0.0309 0.0321 0.0308 0.0323 0.031 0 0 0 0 +26 -1 26 0 0 0.0322 0.0309 0.0321 0.0308 0.0323 0.031 0.014 0.01 0 0 +27 -1 27 0 0 0.014 0.01 0 0 +28 -1 28 0 0 0.026 0.0186 0 0 +29 -1 29 0 0 0.026 0.0186 0 0 +30 -1 30 0 0 0 0 0 0 +31 -1 31 0 0 0 0 0 0 +32 -1 32 0 0 0 0 0 0 +33 -1 33 0 0 0.014 0.01 0 0 +34 -1 34 0 0 0.0195 0.014 0 0 +35 -1 35 0 0 0.006 0.004 0 0 +36 -1 36 0 0 0.026 0.01855 0 0 +37 -1 37 0 0 0.026 0.01855 0 0 +38 -1 38 0 0 0 0 0 0 +39 -1 39 0 0 0.024 0.017 0 0 +40 -1 40 0 0 0.024 0.017 0 0 +41 -1 41 0 0 0.0012 0.001 0 0 +42 -1 42 0 0 0 0 0 0 +43 -1 43 0 0 0.006 0.0043 0 0 +44 -1 44 0 0 0 0 0 0 +45 -1 45 0 0 0.03922 0.0266 0 0 +46 -1 46 0 0 0.03922 0.0263 0 0 +47 -1 47 0 0 0 0 0 0 +48 -1 48 0 0 0.079 0.0564 0 0 +49 -1 49 0 0 0.3847 0.2745 0 0 +50 -1 50 0 0 0.3847 0.2744 0 0 +51 -1 51 0 0 0.0405 0.0283 0 0 +52 -1 52 0 0 0.0036 0.0027 0 0 +53 -1 53 0 0 0.00435 0.0035 0 0 +54 -1 54 0 0 0.0264 0.019 0 0 +55 -1 55 0 0 0.024 0.0172 0 0 +56 -1 56 0 0 0 0 0 0 +57 -1 57 0 0 0 0 0 0 +58 -1 58 0 0 0 0 0 0 +59 -1 59 0 0 0.1 0.072 0 0 +60 -1 60 0 0 0 0 0 0 +61 -1 61 0 0 1.244 0.888 0 0 +62 -1 62 0 0 0.032 0.023 0 0 +63 -1 63 0 0 0 0 0 0 +64 -1 64 0 0 0.227 0.162 0 0 +65 -1 65 0 0 0.059 0.042 0 0 +66 -1 66 0 0 0.018 0.013 0 0 +67 -1 67 0 0 0.018 0.013 0 0 +68 -1 68 0 0 0.028 0.02 0 0 +69 -1 69 0 0 0.028 0.02 0 0 +0 +1 1 1 1 0 0 0 0 +0 1 0 +0 +1 1 1 1 1 +0 -1 1 +0 +1 1 0.999 1.001 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 diff --git a/modified-feeder69/feeder69.xlsx b/modified-feeder69/feeder69.xlsx new file mode 100644 index 0000000000000000000000000000000000000000..f951c2685caf498baec75c1071cd8aafa3c99293 GIT binary patch literal 25087 zcmWIWW@Zs#U}NB5U|>*WC^3s#R>jD`(89{VAk4tQ6z!a!SCX1n5+71okXjt8SCN~e zu!x08V*vvb!zdWhA+RRotRIu1z`of(_-4oPMt?}?NqpE)6(aP3G0S=b@7--Lc{uXl z=Qt+uIB~>mySeT2U!U}4_A%z`A{?EhYgDUJ!$lsH8eBWC{n4=YtxXq`XXWCU4*{LL z3ncjR@_$|SWzn&?m8_bz=!c5xD=Deg)Q$%~Q^bN6chp~waa$^$#rJ97?b-=e6-@rJ zS%yaz{O2*=S`xBnO=b1!Czq=Nd$xU#TIIDuV(!(hh!ruFpWBxRrx~w4$g(#_V%L!w z4w7paR}`lDOq0CK(s4tAt8#*IIY+aaG}Bmwuf-pk5+yUE;)_ zy)jO0mW~};3LCEbRkPV1jaki-<`lE{^J?xNy|qs5VSfEtdK0zP6Z@-s=0{h}adMd_ z#d%Kr!{tS}CF({3Dw{UsUvfUWFa62H;*H4}>9ZJ%8={{{TK!;SVEF%^nE^G+ZVfe1 z{=&e(;KRhgz|FwH6kn8@Q>?ECB1g08h|ZyK&##%}>wnmQr|rGRpY{{quil*?;N)HI z=&@Tdq^`m94*%BP#VXQw_p8eYgh(|ryVvin`*Z5{U7OxlCltRuGFj6tBx0?+t$c0U z^k0Rg>{@C?)7%QzJUk*;cJ=GkZ!_J!DmBXX_G_))v0|b3$@|J~mp{5}6gaOJF2yOl zRHjYLc4lMLG2h60Nk7&xUUdq7yS8n^cJl?VCZC^Oe|_4rm%R2JmOgJT%@jRRoRofZ zUCoIc#jii5LW46bGdCp7E3#6mj8Dlq@Y!Yde&N40x_)y`1Sy#a{V+Q3eq43I*ZQa} zCLd4Ftvq7Rd3g1S4=0k)sdEC*}D?*-82N z*`S&YmS-3>7BDcv18;;f)=WB`b;N+j^?U2T=8O9ZjcNnh3ezqI&B&cy*HAFU!Y7G6 zy7tNY>>#U?2PWigGW0m~U9NNY?fUYyKUCzU(nV7|xti9Lr6y;yzCU{Uta8VJ4DXz| zs;tbcH=bV)-yNPc@$vG6S#}*OlA4(g-^pLfb1R@-vC#IW#mYJFCnUejzPZ>g$+u|o zGX6{FO%?1fi?jzd+E3Z8(6zAt-2d#|tFlaNt^Zz}_>^1B=);9$tS*iF)WSYQt4|bs zlKT7XtZy>sohJW|IQ8-DG>J{7lKbY^vpF9TKE6`)iA7|~&Uj(BAtP#_==fHyh^vKX_vEjpz0Z*wWCizby>!7#SGOGcz#oGcYhzCT*AGA{Au`wu{mxPM>dnIhGr474};=r@dr;BXa?F-+{~tJC212 zG4bw|Jffy|d+PJ5_n}XxSFb&xUJw@*d?7_Hn5ivV@$9w*<=gYq#i#r@(35aVO5+8i z&0d}P_NT31mq}-rvoQAPXUIBDcAfa#G~?OxZH&$4Q%h&D-r1G=M4H{_hm=Bw=+eSH zJtmXe`Z#7*bG5D2I{C2o=~^~nHv#r8R-OCTH+>YI^~}HV@|m}{(vGdo6G&2?%G1s` zMd5fG811%QVl;ErkDmE>`p2H!KfK(gnRB&X zly3X?@$2dD`SY(!Rj-|-kS5x|u)SudozSegyMFKo8Q&}JoA4##hQxnehne4VYxEa> z-oNzC_LTn$y&Ga33pM!M3(aNhJp{8vzUz2?sI#eNFN$t&iaVmumcB^&#lnvtI}=sY z3~Co@x~hDt&~#6pJ6~U9A{u) zn8C=vAOcAm#Tki3sVTuFMVWc&#h}vV$OOkMWO@8yv`s}mrRY%V7U0cn^^KruQ6Y1>7&u?EYZm_8F zxIwxad$`bzpU2s?7cJYU!>4xjQMC5ZJ+Z~YZr$Cj^JwYn}paIbE#^mU7xkm6Ho-X{m=CROSBb{5b7`JTlg z|Ks#M7u)NX{I=e{vG&@3ta-vZ*U0-I0|SE-BLjmn#N(hmQJj&QT2c%xV!<3kP!NC% z*p~i-U55-rT)#8_bFQeKe65woIoMojN9c3rq+hq#Zcm@^UF2`|W|`CjPoEZ?e`+@S zx6MlFG`@*CS9USAK1mkO*zlqHmhOJXjTW!&h;ZuKOh~+}w{I4+pPyIf@d**4+#6d1 zRsKpGt~+vbV;|?r33ktJdUjSCoc!4>ePPY2tlUrW>|6Bxw!YVSQXAfEWg8Uv?@eE= zyiMNX{pUPBop159oRDC*$iwo4>)THaNsh0}zs<|8f z`kTgT$(a-PNDDKmR_cg&U7j@E>7}>4*KG$TU6a42vzIN=e4_uw!}yy=_D;{syh~KG zBhpOf_u0)Y{_?Wy+nuu=R_zuKTaM~(uUYoQ%0h0VVdrw)@2?ibXYAk;*N{^_m2{^_ ze8%qQ+wRxD*SRI6*tdcCYx2IkRs#ylUpqOE(qS} z&Da>Eb+t=SY2)i78cds~s|gz&IyZTKedV{awdK{Bd)*gyMifLEEIPSwOK-0hm$jXS zH%t772DNjFi*&`Mq&6>R+@~4SuB7+C`Lymz`|Ueww{>jWqrrq@HinuV;&E4ElXfYTteS1Kl9@p8SuVdxa`M%$i)@c>S%r1i@0YXL`~T(U_wOPfyk+4n{n6ku z!&>~~#Nt-h3+7K*CtJ#?{ak(Wfe@WU8&KX$>@-fSruDnW6X?Jfe- zTYD8+Zhbg=&vg0W6>G1@)?_FpIw)i*+8?-5@VL@MTw~KrnFSrM`0l$ZZJBr{UPW2< z`eP;~+oMw~_ZmLC6ZSla>F?F1I>sYM<-XbO>hMWXI?;IcnO>5Qd)ZA#fkhc>Sn4>iXh5yY<72ToIl-CXTE6CL4Sd6nt3Mk+$K8{_5TrQ_`0$zcbU2PrE*J&9@^zT5`A6 z&v38Kne#5R4S&9p)< ze%5IhryjJt>XBchcGZ|cO!zd<3FW&}V$93@zv=mI4e!|E=;Uoaxn)|LqSRv@>#vI^ zU+yV2=3nXPzbna7sz~bVvDqD4MJ5X0a!8KV%y=WS^-k2)_Qx9qXKlU1`R30 zb$6yYPmsELok) ze17Y?omZ?pLSE~9oD}(8RE?!K>fQ`TU!%C`i7RhLuJn8>BF-{PJxlt}wpq$+=2hIC zDgSEjnJ7y~o$FJzK3=`sUT1x<&D_CeR{@{R@}|`RhweHzC0v|7Pi6TgM&~EXPF(Bw z@Z^l}rK0?)f3}H-{8hFFutayaBa_Dy=e>Ig=*!l zznAxA&$qWZlFYXLPnAEj?g;;Q>mN6^%52S-Yp)J+Ffi=$U|`^dROiJdl{u;42F#V1 z=>FR-B6aSM>?5XC_|*D83KhDNa$QJhftRk(dcQ|0zO|h>%ilCj2#eXm#3Hn6hfTwd zl~dSPu-F6?{x;3I@sImlkG}C{PqWF_B6oM)>Xkh|^ZlInXDrjxa?0n~f3khpCww$d zqQv#-EFCUC1FqwF25-C8-`OOU8`vRiEpalBOuhaScZlK7q zhblYS9xW*_Ds(?*j9_z2XDUF6sKP2 zyd$}~z)f&zT)IeAqKZm_ltbHjmTm8sgwL`4@!|gd_ls7mZF^sJ{AHx5x%HOC5edM^H5dE-MC9p&grKv-m1%Oz2=-U`_7AFdb{H@OO{oy(LP++ zc`N^2-OZ^S{j5jT}~N`oAo2>AzV^BiF{hK3`$~ zv`&(pqdD#EjAmYmHLutH)vDYv>EHF63a4#E`J6>K%JS-_pJP0ocVp9qiB~e$KX-ZgTn@O%nqDQV!6~fHrDGwm>HYS;^nfDPfW6B~h4$qgzIyqD zU06trv&YM&(o#njKJOFkddVPO6?sG2gn!LJ#&^@5+E$s$iDb7#U0bk~c_Djkae&>W z#?HBpT7k7!oBU6Dte1Wfyz7@kXr|fOL+taIrX-$onm1ubhMjVJ_RM8B!)CoLyM4vf z_4>w%?K9nXAG~2%bS!pbd3%?h*1a1r`{j{kOV-t1~=N}`9$QA8p42x{ut93nH@7MK(QCocPxr;6w&-Oi4*FUg`H79#9 z@1;N5YxbqaUk;fl_r+@4b3f+$oC$}ob_$iZe9Lz)c(J5CgXN=+xOG-S_dq=8T1tWWE-PsT{vvZ1LnuuZGW@z^3xofs<;Oj@*(u(siQ0VmrrI zQ;vhz?wH-xb&5!iIB9ib+cSHKxZoLaqxdpV?Q` zV=&=Wr{FUoOE(EG7G;KgN36@Omx)f%*%ZLH&%=T7&?^VWkT}j&LEdY(Sz9f2@E@yDQNDRdMhlRt9R43ui|f~Zq+@IF~Ntmtn$j9TUoOo&elAc z)jesYX;$l`D-*dUpM1@=XO1o{L$2JJZD%rW>)KUD48P#fC=WKIJM&I{q2I<>P-r8yY$ z8@Wj&3X3{TF}P5bo!H_#Wxqz^g0I&Oam?0u6Tsr2AI+_tvq-5`h%ImX)(uU*vQNd7 z&OA_AuydC06lRUtg;u)CJsVD1tv&W~OOL@JwbxrbvgKwQ9^#mtbLE2NVm6SZ8}D4> zEQ3QFvvnBR)ET?;{p4%Ey9y>s3SDfn?q3-ssp7X;ZhgU~>7JT%8n&85Dc)l8K4>Q` zexd1vSP<8($a*_cC_;v55`}}!(zn^vaIOWzn}Ky=l%X`=p3n%wb478o-gp}-rg7v>BhiE z&mZQj{Ih!j|8tEqMY`ptsmGry01_zaPDgFGXwCl7?k_>JDFE|NFQ0^pw!b4`MGm@3n7cf8SQw^!@CE?~e}L zZR^@6Ui~3#mPW7{bN%Ewue+`|a>?J1;@Xh;X4$PN>!rKyY%qz;eBrzMZicz(;zY@= zUH7m2W^YRkxY2w3@VXO|7A?{L>iD2{?~$*v=Xs0P6}7$C>>wl)x2R=`vEsz(d|&^u zV{7a1DHi`OGG}1$_GDy`fV6dxhAt3;7E5aH7t5Njd-gl*?f(at9)~Th{q@~)-Hq*U zluw3T?Nw2)HlDdBGsiNpH@8yN*Jx+`-@3ajvs4%+-7P&onWahLCdB;_o?|+v6o}AwAZys-duV&x>s^6cU ztiOGH^8bpj`~P46UVqR2&*zWdZ{J>?KHZ+{n)&2YP>B%CwW0TYW*Vg}kx}4uUxb^ie{(t-D>|)F@ zE`7DD{?7gX#aa)UHkGY7`&g%P)!DP{%Xi!7&AoDc{cE`-x&QN5Kl*Y% z%N6?U+WdPqclZ4I_we8O^Z;&?Zy#@;U6S&Zd&{MBKVSXk7k>X*YVFU*R||gp`zgQT zI=j+-#ZPu>^0GVoRxS^24Bj|-iF0X3{okAA$!F@_=5T+t)Y^At-_4-kX?n%!-dWo{ zveG>lrRMv*3JhBOt32B4ZPC~HA(^H2-_rj2zPcDFeX%>Z@J@Gb==Ju6cXjvwdb8W{ zOVWw>@Vkj`W^I1|YsSkvqED?e!gg(+BvpGL(R*|I&YO!g^L<`j42-_`X5|juvwS}% zpDlAendY5!C8y`g*S)sq&Wia~CI|KR-Rey9Kl;<~_P5`EUmSnCqT2e&`wO>yjVs<} zlpfz5^i+Jpl^5JT+0V0T@1Ku%*z&5^)7ENLz}aB;wX5qS%50~F|8l&1wA{9MPOd4> z??m}t&cJ^=A5U33?aGR)?Ye92&a4&F1zV%#_IBlBv(U}?7e6dL{e9QfGdIotK76p8 ziKpme?K#=yIa`;wtzFkR$4n>a@~rMqFTEhOpBy)Zw?~xvb>0biYmt;1_bbr+MNZ<| zS!w3K8W--rWV!LwG3}x|zjD2|MC-d&-@Bsr(<5bL(R`tu>yPbw`MqvMa;2D84O^JK z`IoD8FZilzCfi)i3txGpD<(hc&OzI+i!)A2&#qC~#9*}CC2#fJ8}c>F-y5W^EZhBQ z!?&%zZ`*(JaVZ_Q^bB2l{?WuLVFkwTJMuu)jIiZ`NAA)_;Rgf{x~%( z)LULYX=#;p_oG+727GslH>*zYzR>S8J)!S~^Ho!4@dM#KWjj|0XNNDm-R>NBH`l^l zDCqA#&I?!hB(DlxQhnjR%sayU%feT$%Xl_)eP!g32z$Hk{@dF;A{Oc^y_h!@th?9p z{D5tCgG&kLuZb^(t_s$!Xp4;8|;n} z_mTy>_et}`@}1kz8u|Dr&6V?DEB!$_paQxV-4zUwWe2g zSA97U=KQ}nDfH|)rrl}68<@Xe$+_8M4i1hRhO!y56}C4nzL4yXNL+jK?f2d|uJfIa zrqvdtK+3(lXNc0kY3-YAx7_s$TrJ8xERe0wOiSMET1AY0CVt2Nzm>^zfYr>?Pc z6MJVm`+?k!=lzNgAug_VX+JBOJ$K{er6mg%E|g^wJ;1QPqW^|^V|~37T;u1 zT+O!FWb*=-YLMJGCxy;VyA$Jp5ioiyX6!)fAtKJ@KnmQ}uYWRty zU1h?{n7{uKRapl1+Kj#zx)ZXvr(gNRd}iV^>!hO}Ud;Y`@5T|q=B%=NojcBUzi0?j zHjK@_-0zaFx-9wMs@vDwFG{>R5)kgZ-P)vz%SN5;>J6Pq+2)Ft1^%@<+W*x+aiL_n zsNBFN$$rx3^>YNXgWon*F-as#yzU6CZO~bkZJuaVa9-Y_vN!kF#Gk6ml&gZT9?)34 zu}7GB+Q+;%!|sh8Cufu% zSXcM+(ye(4Y)Ku8a zxqmC4_Ew?&hOtLZ7(^a2zu?S#wo7#G!#iD~ciW0XrXC4@VRq26;Czs5$w-`L*1%)9A}#%6 z^mfkujn(gRJMRm9QNFBNr(D&`T76)d>YP(5b;@Bo*33F=cWhVDX1_JBm{TUkZ09`L zSpAPd=%nfk_eJGGC(kZ5>RZ2hS7v!5=eZ5nT_k5^Wt2_(J+~`3UV;1luH228e3~c4 z_BKv=dD;V1(A+buHQN5(>iuLR8>1suy1id3BtDcLxm4jNv+8bS$n{*O(g*Wys|Ut{ zeB@H0PxQ?4UT` z0Wt1gjMn)-wsZ7trNpI{r#R$JQIw4GI=WJS0`KIzlZk>y4^2q-~d;UZ!2d^8J@@0so!i%c^(4Ii#Tfg>93z)8w|Le~)x= zem0OS`+9F0DDM}nzt8d=l=lzq1?7Ft?;U(Qw%@+KHgUUBTkTf%Eq8xbm>gc_xMtPO z$O+f;nJOR5o2}~@2Py&e9px@FzW?SJ@b=(_AR+Cl>hgwQHR0AOaKswM?dN>lSpAD( zlVg5|UVpv8?v@Il<1zbhr;2~t-1>DxHq+~V(LDDxtN#9U7QXH0Yr$~UdrfbpGnlV>>qS zrs@mjUr!b{t`4nUHmzifMB)xhZ@G=D+|qR`|Eeuvf4B3Z&`WTS3;oRro4c~yyy!Zp=H~Yjvqg0_Hg^b#zk?Dwu8dI^M1!Y-*ZP|1rn1h zG&{B)oj9xNe&yr7g>R3nvU3nS$vb<3%#P=E(>)+%ujL~DHtQp2qt7l|v3J$W-({cf zeDvopnzVG*T-Tz*YAWC)4bE8nLMH|H&U8A;*ZBNkP)*nU)w?unZoI#}wr=`_v$8S^ z4h21~6`mw5qcw@QYBH#pS}1d|kb`aC%44s}KmAhLC$TI?-YoBd_@t%3Y`Gt_2^>)T zu&1C-@=fq@exZY4AF>=TnIkOa%(0sP{~X;l1=qLTvreBYtQ7W@^%Hx7Gt2fz+!iS& zE1Wl&{k0H?O|n=dzu`CEdFP0>-L0a<$-K+wY`GS=E2~WGu&?^3>B^1Tv)(Satl%hf zzRDLb>aGrhkpK0N{&?zyZUZTWgf@FJ6&bg873Y!?1sg= zS^u6W>3lw;Z}XnNd~t05B8_D`+jzMZkHkIZjb`->jaTS?-X?H2^t~d3|i# z(sY&RCW|S)Pu?%tc_GbS)ktp&2fxhWqE*fZ#ja^`I9MItEAv*f{&g{1weg#a4?4o) zCo;_mn18^&*zhu8b7!x)-S+YQtccdOJ>gGwoZs%)wzsyDZDQ&C*>9E~v8&AG+FR_l=78tN z2%}@~m9}HOPC-`%q4V%e4zOdgMOZ^@=M- z=yjKzdNR{G#$ek#w`=>-FMhka>UZul_6_H4^>V+uT(`~pJYVSlD!Jtn{ie3>BKR#P zhP^4AwrkPx$X8u*=_|FX?reOwxt7^)ofNwIj_Jl#Gi|d?Z{=^6YOKA+ zn7LJb*EPo3|ED}v$hvod=LnL)diJ}!!1m;0&OwQqaiZJpTaqGcZ9%oj_`Z~j?+q(58baewBGT^*u^;#rfz z1ZKY}xRD~gU>4I3!QFuoM)z5_YJa=0w3X$m{`9slNwQ(3dzHen*XJDL*!!I4+n$;m zX#&i@ZPb?9x4u|4S@_R#7k6zJ_qm09K8LFoZ&8tZG%+p0c<&sg9ZF#_%uAofR9Le4 zeLE`^cnA^U9p;Q*DbH2bKUQ=&+JwIZIeo7NS52GjQPE^%8^RaWtQ>8p9WKNLOEU7vVDDJ*l^wv*S*e(qH|zIW%@q|ybujhk@C~o$jxQUMgcb)xHvY=R|XtGnu-PvaUH9sn=?f7!{yx6t5iC*7| znyS~HSy%V^=@F0U<^|<(S)SMR#KMvXZZTuci+2Y%I}+4oiE(@Lt*uA3&Fxk5B>YEZ1^@+Z>n*iax3J`Hui1ailLwD3H#}pgw_p95W0iG}xV+fP z%f504{pXj}yk2@Nd#c&fz0pggj<@KESz0~%H0e@7T*P}l&9krEGH=afe_4HQLWI=z z$Qr(LD=vD3>s%}~zx8MNk^O2hX^FheCLiDSHgM}-y_L6X)dSs@MVyK&Vo$aRTv=xz z8}n)VY(CcK;_tMRj(YY6#jfuRWxShH@-A`qLyk9}AB0v1yI+sjn6bRN!U z-E&W*-~OY&?D}=gPgt&gp3`t!>geN;Ih)d+rN#N&J!PQoFI~0s>CW4yMYT80YrbE( z^V*xYs-{(zJD-GbE*ARon&t1_Wyk&b-fo*Yqu|iR(&aW`hdABJZ&`l$AQZbw!F`jR zW8?Jc9v-4{eC3b5zN?*n`}k|@_Y2Yb+dE|ocDY}etrI@u=$D5Rxb?$#RGR2lMY%QY z{_EG&H_zDY`fdP>Tcn+efuh*$;i(n;N~Qi3tN|UZ`^xUEcCs` zG>>_Tk60ETe=GB%kk6ZI{e$`Y`peE$-n;fDS%2jgBle}!%)%tj>X_ftmAn`ran@Vb zD_Z-@#>aE>)UTYly*w+TuQ6A8nMrKPXS2&9HGlKCrr%EVGuUg435;wYIlx_4Tg!HWB%#LtMIg{OL*S+*(SEag!oMo3u zT?^lgWBNYtZWd_WnZL2HX!#GtuT=^xs%qud>?I96)|Yk7KeDvim1pj_v(@BdveEJx zx-v8LXU=2ZS-4(m-w)m*!wb`YO;fAmDKYxAmf`6>zHE~NCQ}@lHysT7nm5}w=h6Ja zop;)bn+&;cST11nmYBGZQ7?Bc+lH5!y=*!~;TZ?Ac3x<3eaFQ+;hHX6o|bw_ESQe+hm3a~Ab>KB4>t&W>yB;5%`)BE528q{kd6jFbcd|}s zcdeb6EtA8Uwco# z%!@Ry1F>sbTlw#v*;B!6oCEeZ8HhHaQ{!wAGtt(q9(*!IRFz7K}uYNJHQ7=Bo zJ5Rgh=b`8#7v`P1d$f1eJnY>aqwG|&bG@^ptKlnOWeJ1Xp6>T=xf%QR8&-YUQ2O(a z;k`#iwoA?KPg%fTuEZW!p%y4q{(2?1zm4&Fv*;5MyB+dQ3j8*k5ntD?e5>%sO=SLEOIPZIZE@6*Bd%mn>wG*3JDgrEGcG<0%O%7$jFbiB+q(SbS}n z=);2njK{mxA6{7*pPKo2`>|AyHQVMnvPP`y5V;&@zT!a8{R+8jdjc4zY5Tov%$u zQcUcznPT%zE7VJ+LfhO`CWPE;{H-N9FLRoSSLx+F;b*wBomYmGY_=(Rv1i^ay~#w%WOIn)Qd)&4P*N4+)m# zZT8C^oa5kFX0bf{#xnijGJHqh)&ZvC8f{LHfyM2|HtD?q#wR8SCHc4L0NM{d{6YX3PTxea(cOQs17w zE&A%VTjs}d+o<)4^SuoYG|Xj|dMNZkW`iKxo8Y6pGSh7Dd3ViCcRMm`di3gdoC{$4EG!#>+!)^w>?f?c6^w=76xHGlcyY|ByKi>~>3TXtMM_ha$8iN}PuX#Poi z?vTH!_f-JTo9}y+60XW@>AsVcQLnIl%cf+G-@n|HZ&;LyUfce{U3F)b_q+}JZ~l3| z^URdOy*DMcvi+;93X3lCy71U+;^QsLSR~?4*%dvK-I(N_9Ju(&XNSBym#b{{rpBGn zl{|BwwdnNkuM74C`bQ_vstEPJ{Uv){`=RJAd%4SJ-HfGwJfAnMab?AH4Z~%ci6W-X z{>2S5_EkE)Jt=l?SqZ;<)veRLTMlhH!&1b(=fz9cyggT5-TZUtvzcm0`N=&N%T&vs z{CT;g+4Z&gFZ;xr;F62~rF-*xjsLBXzLf1{;AVA~S54plsZ7OTRksz>md;&bowJT- z^^`7kmnq6#8A8X*LfzW}qxv?xu1mb?EWdfli@9eHAKa1DQu#H*EI0GT*_K{NsBJVo zH(~L++T+Z#S6+N0#8+jqn`v*Hp^REOOMu0c~SgK_ne!=bI#o} z7j+lzyO?*)VsZDbPrH|~EWZ9$=Ec{eY$1Jxkri*_zE3)}_uQTePJ=&@hCTb*($?!w z3j4C_aqRWiXLc{UzLdAeMkD0&lZg8n(%Exz-C~R5f}1YBF5~;M%Uj~(vJH0@pNKv* zLtOoA*Yy^|_YKFhm3Q7M?VNi1D&zHb4gOs#n4eXxTfbwK&;1WArE{K&@oOc>C3r{v z6mocW_503*8|rJ)y7b!qx~+~>UTfxb$L+Px<1=hpO{bR03qG)F=6Kq(@n!lM{RONH zjhX8joMx1i)Y%oR;5;ia!{dK!xw(8?&7ap_zI{CWHaun@pNbvth01_^>kitzt(YP` z>BXmQMQ7#fznI;Ln(obB9kXkKSl<0GkK)-ZSDV**+FvvKYvO2sxg_v*!RC#NwyAN{ zH9zkaelSI-=CjL3MXR{LJzt8rG8I2XDfR|9_PoC!F}cs*|3t#BN>{lzk8Of;{7iCh zsrP+s%MsApzGzjYxUaOHv-XmcJN;^o75GIf>e?^e`f$oNr>MGR8Qaf_B%JTPns@sI zdv(~YZ4>W5_1|53+fD7p=F;0{7k2$|c)g9+WJCIEu|8Iv2}@t@Kh58!C&QC{-}9LF8d-i!K00IV&3?-KvgL;MUk+Rh-=o^Ru6&R7yd@6Dt9G(7v$C5|m8w1y;Ng$HF@<4n{F20%^X-x~ zeA^x+&CrN!JhiYua9y+Fq|+A^+d_G}oH!31u;gvd5iGDcX19}NTA!R}ys?hPi$fA( zQx;ugu{yZ4t=04Rh6M$JraBgP7XP@z_v`YSZgC~&^)`%Ts}+?(6PcKSj$mCQX6T6@n= zYqvBwXYjIin?ogQO4jEd^Jm)+DzdBUTDWY~@p>;C>%$wtC$m@l0cY2HA(!b1tCStJ z<6eE(b8NSBK%)aB5w8rbe|H^A^)>zqC@%r}_KeJXgYgRbNCu=~+;4T?Si(@zdS8{Y!u8 zxvsF7AaK-nrRM#Mg|jo~Fx@TJOp0uISNtIG<4;weUuUxS>}R<@XS(|ivfQ0t5xndEp~C`+7L5D`VrSl%h(*a<4og0G#vz&UX{w{G zo15A0J#|qrHjTE9X1f#Si#6)(4?j`CZ_uK)T67lwzfAAs47;TrSHH#jmRyt-_i<;K zcBST2V0Vl6xlNC@F4BG$`%&v^sY19r*B*YAsd`?&4A&$Uy8E1+lNBMMW^cXoW#)b% z{Y^*WzvOqmS@Y}0ykqUcPaI}s&$<;?_S@@*^S^l-{qB58mQRcVv>eYY5BQ$bXSVmT z*sJ%sYn*RNmG_n3o&4&3kjv=>|N8WfXGhTSm8@%S-e1aL+Wsq;W9|`Km+ABLGwfqI)#v|v*dp2=<@kH<1kv|154dDBeJtLOKfA8ETJb(AcKOE=l?Z8L_UTGxo);>(@&z(vIocPw zyKk&=e5~1BR8-jXHAA*!4zt>o+~3_tR_RBIExdGP!l&J(0hxVY-gv7oG3+#6U-+7N z*;`H1rAJEj4NaEi8Z$?09=-lg^`M^?&(BA@S7_eR>2+_|boTbt&r3Vki%vK5eZKS2 zs}BoW7u1DdmfEpynt0_39oURX8yF1Rqoox$4Ypt83(Bw z+T}tHUV6^u1`d_mSR$@&4HNJBuvRSUt+UUyOw*mAy`HmVR-Y2K4rmI9J#$GVJN-tY zM*oj*7qX@=cpSIz$duz%Cj_!iuRJ=PXW^0VC8EFmvRb}QQ&fBAI%%6jQj*!+wl_9H zmzhPbraaG#yw|Sri}7;FIoHV3+e5Y{s?OePHpOE4!h zYmqrJt!+({qXTLr4+wY(Xh?DSZ7Tcplf~!u6!)7qG!u_D_g|<`PvJVF>9gQM{WHZ&M?d@F`~;tcKbBm*yF2Q7Z~Kzfr3IQ>t>h1Kp1SG9$*cD*psifXP5NH! zp)Pe}LpcwRpJ#n~*!8zBcpSkWvF>A*nj3?+RvUMM?FUnxo`|0IIr~H3J>6*amP z#-yTm-WoO6ma_UxK7H?i^0GH?ZhrFdnXJonyfx`@ZWiMQzH6mi@9#8;<&|rG+Uu>u zFh#*^Y6QcHqZIYKq~C`N|%#3%+YjWv@T9*-LY4(b}S>Piun?X)rt0&3F{`JZQ?si9%sEYpRYq z?OOF?)2pq|b?w&vvYQeup<~Zhzd{D#fK(ahx;h* zWS9{zbhGw=uUXx*_pNW*R-epUtM9*i)+fscJ^p8YK6bmW^&^fxo7`rOVXpshdtTglgho<3O#~M6ZAA zp0R8^+W5^nGFCajb_4&^MQ>(K+3-lx_sn0v6J6pH48>&L$K6X`FyG;r*M<|l503fW zH*J~P@Kox2M$)@<&tL26&iu8EfAxiNiRy(ZL|O4WGgk>l~aPW!?n@_*M)dVe9PiK9U0PLqq@dQf5!*0 z`OR&9lf$X)`$bTS`oKhH-vTwe=L8Uu9)wXUdK-P5Am`k&P+C zCP^#HB!j2oULE`AZ(OJ=G9!l3Y*(iHC$B7hvE>@(Zzrk#^KdZr3)5~C2;iNnbT|4! z<@8N_*DUUdu5WgJ9+~&CwBUu!Y2AH%eS!yyaw85;>bZWfF}o-(M)>-ZKRbT#oa)M) z?!bJI#q6h0nbhj%MmGY~;$1%cT5`Zh`tQc%OY*yy9O2@w@r#^dKTW8uxogqFDPgVu zf|_RS`Fi-pj6F9+qiQ5%yWUh-KW^m+%+)YDt!vlA=bEUYAYg0lS<4jb!^+vyAAIcFbzNaTqkf)7bILZaQz`Z{6}Wpn zq8^x6#U5XBeQNNPond|fap7JPOLFDB{QgBKa#yqpeRxpHc-86T%|behf6_kTUFhyIKBpsw?2&Z0jX=PdZMams}J*d~;JEq zE$>l=8!!6rs{h?MN6FCs^3}v}FM<2-+Go93yXizQ)1%;KTb{jY=MQQ9yt4RXh3Lme z&%R9xWWDw7R`rTk5%WtMcuVHnTuR+};je2RC->hs6_b~{JlXyux@msM{{`zBnl?S( zGt)G!j`PE2$pEK6Yg-&2Utjx_@y!&jFUCs)-uKE~uf6p6l;)|&r#K%e$6lGTV2Y@e z-#U-f`O8_N567xyByqk?wp0K9L5ny3$+`<4t`;qP*nLM-YF_dFRs0;h=Pea{dC#XB znC=O1@;tnKs(EI`!UId%^1rNlH#siUFYtJ{rDJ8NW}?lC8BAU;LoOJoOT72v;^r6V zS*@j)sOehl`1s^O?wv;`oN~xZtnaq0otPba>$&R7uMD5d&aJ)m(5bpLvMBDnnEQgI z@mKX4zb`eYS<$&W!0FoO{!Z)ZvbvY1&CKC5(Ja-d-hXi;*FO=1-JppN3hoWpe1ShS8IB&p0j+4!MiDc=LS9bcSSe$QwgtV7m=z-kZYl$3$0`vGvIvuyZ;4+A=*Q`_AGhEyYrW z<{tf}Z`vk_XUz!t)cfvI?WI^R&7kQEpF9x_T^^;$$i5|{he7qEwD^LjemAK`3p=hD?*Hp%uKgqUvP~uD)s5k&l=A;FZdxca*Y26S z=OxAg>%&*mf4QAf8 zkF_%5dKlKN-(kbMdEW}=3pMA?uUebS~ISwKc|>=CpkbPPhf51zUN#U z`u>P|aWo0EvZgH#E}dVpDyZv9md}qT3ZZ_X(@HP6&(rW!bGYj%@#jJ7NwHUcrE9-< z@4N8cb;48g*iGA(PET5Rai15f4r_1RsuU&h%jq81r4BfUSX?#X_Ln>=z1R9|rH1!E zjWZ8-zTsWN^zm721Ft-mQ(AebZXi>4}AM1&}+f8 zIsF?gIYN3Ot&gvN;(XZri{CHmWCPVD~hO||-a|31FoKa4nE-`F-?IyvX= zw)u1KhT2lrZT@hj$Byav$+?+!%t=~{JtKd*dOtlX+&g1Yz3-O(0MT3O2@-A& zD-XPweQo7p@AlrfpbyL0SeANhJ>F2D*~~R7G4a!@*YXzL>y-{`oaT}}UQzJ$(XIay zF&USintpc5^E`gPy{LP4Sp5A-egB_4E&Hgc(7VY%&gZz=L)o?W>b+AEZmbjN)1CiS z@Y%T+leMyI%&PST`?@5xr_Q^S+b4NKT>jXN?x%^1gTGDyM`tAu+w{}`N8>pQ&J_eZ zt!Fx+k;paG?9CMKteo!qFHeYP>|ZUSa`bQ3O-=E*GY${WTRnOBbedA>g1<5CYS-$D zO;1#YeqX#%>PgLn>eQukUQdd$QHYm)bs~MmXYtz!-_D2eckeek`n=fWkxZhN+p`T9 zk8=IH^2UF;!O`iB`+}ReCM{%Hcy-3FShmVpzb6S7IJX+rYp`o`N)}cZCw8@c+RdbB z-gkJ5(~&(fwwla;HY{0}a;SZmalSIsf-3>XQsx(LKlMA!@aprLnaPXe51vfVILht! z$|HWIP3ojG=`TK&EBRble&Imr`Uf&*R!Ilv2h~ki;#YC_sp;r?A;9VGS1p%GlTw(z z%;%rbWxLw)r}26lg#?yf)}n4kHtzKn+8KWTV%Q(Q(D_ooa-o@_<740ZDzl4^m^UAL z9y7C|aX|rpyj`{ae+d#n*>idh+~i>orr{XT*in_qK8^&E7Kq?1HA|We)?g65c4r z9Gma(TiSM8`;z0+f}ThmV0y5r=T#T8*MtkgAC5nr5blzf#KQO zJP*z*nZFCS)lJR`DAp=mxSF#~Z0FO8p9}VV+;@q|d0F|SX>EV~j-Nh%X!10Lt&t7? zp8uJzu~0xUUQ0jQ<51j!fGW+neIYfkzFbLM{9WwN%B4+z=D%uVK9wZl7Ry(AYj)9z zkIfe!6dl+v*Uwh_ST>C7S+19Uf1*jn&Q=bu9ya~$2fEiDRN5|bv~z0G#OPaz%o5-vS(Z7%>3Kwh-1e!HlHFNapPg(Fca;eg!-*3a5j<7X>r<$xfN~vtd&3&dr7%ttvW`T)s4!wK}bhytLrn0l`|{ ztz~;wEG$(kn7E~1^`-Z`1t%=o^X@*&&trdb@~6w!RHlH}F>IClO{Q3B-?#V{GC{xn zEKAlWUxE7<%3B#srDM0XFF8FeXo(v0l?mx#%~F#@-g_14`POry_FO$sd;7{h zv8S0+^{0df{mdzQyvJi<+sg%eI4+v(&h@+5X}|Q7o@>@O>vKx{Y4HV4`=^`v_)Yv4 zuwsuki@)oXc@y?TojMxp7PL~eC7x>yU(&~`fk)%k?bzhwuJHfG>+aK656cVmtvKu- zm-lg*`WgmHdEI?|JvTg49*8=`es8g~e{Om-TjA5UWkH;mfUk9O?*EJQcehc=FoEI+?{ZsJ!>&vH{K1U{}ow?}oe20Wxr1>p_F0BQ5 zQx<%e*>YA^=yv^)b1c6X7TrBC??>L0OPUoYIJxA4nBKkl(w+L};VXeX>sG`F?RL5l zl3Tg|g7uzHPlZypaH?N<8MoGX#>PcOo>Ledi=SRmInHS$wP0b3;tu8DnjHnD5rNyE zJPOz?Qfa{=-zT(0e7f*4hA-W_TDuQ-?NYtoc57+LCyQ!lPNPTvn}lD*3Ovx$5?{xs ztu;MH<{9JPFKJ&M-yG|PReFx@#RuqYq3IF^~C2zb1ukC zQZ|i|3A>`Z)9i}g(&_3yrThZ2nN_bm>}FdO4ASy+QnvHkpp+?wb~36fwl)fFa8pr# zklXD%H%3Oo^_bh@qAhGIjy4MIda$DU&%;@{-OM+HyP21skl0(k{rhZxmwkTw?cws>|J@!+aDL}6 zcW^QH*ZBAJUX$pJ05|ns7vqMcK-8g9fY%>On!xx%Z zp6Z;N#8Y1_E0OuXR$z+pa*@0Xb;~v1SJ zg~!%OEXaDZeCf5fHB+k@|D>(hA;b0O`>qE%EXqZj-0H+ERew*KC%(P>t5|5|EZxT? zlIg3hVkfVj^W^n)rvGdqiD%B=UulQ+Oh9&~^`~2y85lIg7#IW~8$gp2bCR7i5;OBa z=_+T&?JVI$5x4iO_Z83VJ?m)Ao*C8kKsdv%+}X;kF(t^z@ypv-QU1jBwaQ;kyZ`<6 zC_bqA`-jIB61DHo9sYOaeD}+px#tD^(U-ZYg(c^T{PL@;}**{=HjPGEcm}|MKlk5^rY)_VBH*?^|!T@?TBS+>LL` ze_LhGnJOb%tyZqJ0@vNQt49Y+!I%I=~LCooik7QDy^SBbIYE%xJ#d^B0W#PdHN(x zOSSV%oYN(v>X#m;->53`{PQfdA4nk{?cKCRW2pLpHdE-Fs{XGZ1TCyP_R zP0G*x*>NWB(?z4|lkTVAsLJbZo?f%%3`FToZ@Z{8{ht{#5lWBT`?OX^KG4NO#9hF( zMfX@L>r_jYQeOf=5?tH|o(8U2I_r!ey8NPi>jSJM;#ZSLZ0vWz<=~HENKYFpb8Tm%&ZWgJBW^>pOHw9rv@2N6Li2@0^MUnk`!{_KES>wfw*B*3^&yQh2u zTZCs;Yn1I+TXEM+Q#F zau*L~iBP&~@87BH5tAu9r^426b|90V++{D@=Y9L)6H*=bq~4t=%4t1yv!(aCb1Ul% z19v_Xj$`IIaj~O{XSUJ8bmf#3k1Oh2Kf7zfv&1I}f7xN(d6<9a+UTo1e~g*+vDTJ_ z^P*Y%K*b#ss8B^UEdE56c{UIw60GxAQeCKCM6NGvAe-6<>e6 zds4$v3mxz5El#_Z?lG|X@q*jO@0rw^2eE+ycR9Yq8tj;`^14NFZg>E}pqi4(bPMz)WZ|mpi z74v`X`g?xU{>$FmgQ||6^ZD?8zP3FmrGihiW@OT3z;z5gB*qXVmP6?SyipC<8j5@< zHbN@{!)9g%2FTIB;L|A4HBMqjJ}4Mj;|(_O;f@f!ppzm{jtxdP2mP>Egh?`-U~`a7 zK-Y|ZoG3!GIX74{#^IvqCRn4L%ZV_-fCp>>$a?S@7dTGoL^rS;?U+l1fkwPw1Hs2# zA{}~(ZUFkxVh98N3xf?n95jZml@INt4usa%GGMJpo&o#G5LD8kn}XVpK$yY6Fgb*Q z0i#!et{uI#h%n$y7}yagZAWz7=rtEY_q-^uZnSy}YdwZ;C~8@cFp`1cU^J4Cpe28R VH!B-R4KD*PLkcqkLu(v}2LRVX@w@;4 literal 0 HcmV?d00001 diff --git a/modified-feeder69/feeder69ieee.txt b/modified-feeder69/feeder69ieee.txt new file mode 100644 index 0000000..dbb9240 --- /dev/null +++ b/modified-feeder69/feeder69ieee.txt @@ -0,0 +1,166 @@ +69 68 1 20 0.1 +1.00E-06 2 +1 1 +0 +1 1 2 3.12E-06 1.50E-05 0 0 +2 2 3 3.12E-06 1.50E-05 0 0 +3 3 4 9.36E-06 4.49E-05 0 0 +4 4 5 0.001566052 0.000366868 0 0 +5 5 6 0.002283567 0.002325993 0 0 +6 6 7 0.002377155 0.002422078 0 0 +7 7 8 0.000575259 5.86E-05 0 0 +8 8 9 0.000307595 0.00031321 0 0 +9 9 10 0.005109948 0.003377932 0 0 +10 10 11 0.001167988 0.00077242 0 0 +11 11 12 0.004438605 0.002933697 0 0 +12 12 13 0.00642643 0.004242692 0 0 +13 13 14 0.00651378 0.004305084 0 0 +14 14 15 0.00660113 0.004362486 0 0 +15 15 16 0.001226637 0.000811103 0 0 +16 16 17 0.002354694 0.001544839 0 0 +17 17 18 2.93E-05 2.00E-05 0 0 +18 18 19 0.002043979 0.001351422 0 0 +19 19 20 0.001313987 0.000861017 0 0 +20 20 21 0.002131329 0.001408823 0 0 +21 21 22 8.73E-05 5.74E-05 0 0 +22 22 23 0.000992665 0.000656369 0 0 +23 23 24 0.002160653 0.001428789 0 0 +24 24 25 0.004671953 0.00308843 0 0 +25 25 26 0.001927305 0.001274055 0 0 +26 26 27 0.001080639 0.000713771 0 0 +27 3 28 2.75E-05 0.000134768 0 0 +28 28 29 0.000399312 0.001952886 0 0 +29 29 30 0.002481975 0.001640924 0 0 +30 30 31 0.000437996 0.000289501 0 0 +31 31 32 0.002189978 0.001447507 0 0 +32 32 33 0.005234733 0.003513947 0 0 +33 33 34 0.010656644 0.007045364 0 0 +34 34 35 0.009196659 0.006080776 0 0 +35 3 36 2.75E-05 0.000134768 0 0 +36 36 37 0.000399312 0.001952886 0 0 +37 37 38 0.000656993 0.001534856 0 0 +38 38 39 0.000189673 0.000442987 0 0 +39 39 40 1.12E-05 2.62E-05 0 0 +40 40 41 0.004544048 0.010617961 0 0 +41 41 42 0.001934168 0.004520963 0 0 +42 42 43 0.000255809 0.000596473 0 0 +43 43 44 5.74E-05 0.000144751 0 0 +44 44 45 0.000679455 0.001713299 0 0 +45 45 46 5.62E-06 1.50E-05 0 0 +46 4 47 2.12E-05 0.000104819 0 0 +47 47 48 0.00053096 0.002599273 0 0 +48 48 49 0.001808135 0.008848508 0 0 +49 49 50 0.000512867 0.002509428 0 0 +50 8 51 0.000579003 0.000590233 0 0 +51 51 52 0.002070808 0.001390106 0 0 +52 9 53 0.00108563 0.001105596 0 0 +53 53 54 0.001266568 0.001290277 0 0 +54 54 55 0.001773196 0.00180564 0 0 +55 55 56 0.001755102 0.00178817 0 0 +56 56 57 0.009920412 0.006659779 0 0 +57 57 58 0.004952095 0.003281847 0 0 +58 58 59 0.001897981 0.001255338 0 0 +59 59 60 0.002408976 0.001462481 0 0 +60 60 61 0.003166421 0.003225694 0 0 +61 61 62 0.000607703 0.000618934 0 0 +62 62 63 0.000904692 0.000920914 0 0 +63 63 64 0.004432989 0.004515971 0 0 +64 64 65 0.006495062 0.006616104 0 0 +65 11 66 0.001255338 0.000762437 0 0 +66 66 67 2.93E-05 1.75E-05 0 0 +67 12 68 0.004613304 0.003049747 0 0 +68 68 69 2.93E-05 2.00E-05 0 0 +0 +0 +0 +1 0 0 0 0 +2 0 0 0 0 +3 0 0 0 0 +4 0 0 0 0 +5 0 0 0 0 +6 0 0 0.0026 0.0022 +7 0 0 0.0404 0.03 +8 0 0 0.075 0.054 +9 0 0 0.03 0.022 +10 0 0 0.028 0.019 +11 0 0 0.145 0.104 +12 0 0 0.145 0.104 +13 0 0 0.008 0.005 +14 0 0 0.008 0.0055 +15 0 0 0 0 +16 0 0 0.0455 0.03 +17 0 0 0.06 0.035 +18 0 0 0.06 0.035 +19 0 0 0 0 +20 0 0 0.001 0.0006 +21 0 0 0.114 0.081 +22 0 0 0.005 0.0035 +23 0 0 0 0 +24 0 0 0.028 0.02 +25 0 0 0 0 +26 0 0 0.014 0.01 +27 0 0 0.014 0.01 +28 0 0 0.026 0.0186 +29 0 0 0.026 0.0186 +30 0 0 0 0 +31 0 0 0 0 +32 0 0 0 0 +33 0 0 0.014 0.01 +34 0 0 0.0195 0.014 +35 0 0 0.006 0.004 +36 0 0 0.026 0.01855 +37 0 0 0.026 0.01855 +38 0 0 0 0 +39 0 0 0.024 0.017 +40 0 0 0.024 0.017 +41 0 0 0.0012 0.001 +42 0 0 0 0 +43 0 0 0.006 0.0043 +44 0 0 0 0 +45 0 0 0.03922 0.0266 +46 0 0 0.03922 0.0263 +47 0 0 0 0 +48 0 0 0.079 0.0564 +49 0 0 0.3847 0.2745 +50 0 0 0.3847 0.2744 +51 0 0 0.0405 0.0283 +52 0 0 0.0036 0.0027 +53 0 0 0.00435 0.0035 +54 0 0 0.0264 0.019 +55 0 0 0.024 0.0172 +56 0 0 0 0 +57 0 0 0 0 +58 0 0 0 0 +59 0 0 0.1 0.072 +60 0 0 0 0 +61 0 0 1.244 0.888 +62 0 0 0.032 0.023 +63 0 0 0 0 +64 0 0 0.227 0.162 +65 0 0 0.059 0.042 +66 0 0 0.018 0.013 +67 0 0 0.018 0.013 +68 0 0 0.028 0.02 +69 0 0 0.028 0.02 +0 +1 1 -100 100 0 0 0 0 +0 +1 10 10 10 0 10 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 diff --git a/modified-feeder69/lineParameter.txt b/modified-feeder69/lineParameter.txt new file mode 100644 index 0000000..b8e6825 --- /dev/null +++ b/modified-feeder69/lineParameter.txt @@ -0,0 +1,40 @@ +222 +0.3465 0.1560 0.1560 +0.1560 0.3465 0.1560 +0.1560 0.1560 0.3465 +0 +1.0179 0.5017 0.5017 +0.5017 1.0179 0.5017 +0.5017 0.5017 1.0179 +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/modified-feeder69/~$feeder69.xlsx b/modified-feeder69/~$feeder69.xlsx new file mode 100644 index 0000000000000000000000000000000000000000..1e59d2c39c31c328d4a871218087efa2b06d9e42 GIT binary patch literal 165 ZcmZQh$*oi%9WXPbFyt~+GAPgm0ssv=4_g2L literal 0 HcmV?d00001 diff --git a/run.m b/run.m index fd5d8dc..2770951 100644 --- a/run.m +++ b/run.m @@ -4,10 +4,10 @@ function [JMeasurement,AME_Volt,AME_VAngle,AME_PD,AME_QD,AME_mVolt,AME_mPD,AME_m close all clc clear -lineZ=readLineZ('feeder13\lineParameter.txt'); +lineZ=readLineZ('E:\算例\modified-feeder69\lineParameter.txt'); [ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ... phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY ... - cap]=dataRead(lineZ,'feeder13\data1.txt'); + cap]=dataRead(lineZ,'E:\算例\modified-feeder69\data.txt'); a=exp(1j*2*pi/3); Tp2f=1/3*[1 1 1; 1 a a^2; @@ -160,9 +160,9 @@ mphaseCSpotLoadQ=phaseCSpotLoadQ.*(1+normrnd(0,sigma,length(phaseCSpotLoadQ),1)) % mVoltpB=VoltpB.*(1+normrnd(0,sigma,length(VoltpB),1)); % mVoltpC=VoltpC.*(1+normrnd(0,sigma,length(VoltpC),1)); % -mVoltpA=abs(rVoltpA).*(1+normrnd(0,sigma,length(VoltpA),1)); -mVoltpB=abs(rVoltpB).*(1+normrnd(0,sigma,length(VoltpB),1)).*exp(1j*-120/180*pi); -mVoltpC=abs(rVoltpC).*(1+normrnd(0,sigma,length(VoltpC),1)).*exp(1j*+120/180*pi); +mVoltpA=abs(rVoltpA).*(1+normrnd(0,sigma/3,length(VoltpA),1)); +mVoltpB=abs(rVoltpB).*(1+normrnd(0,sigma/3,length(VoltpB),1)).*exp(1j*-120/180*pi); +mVoltpC=abs(rVoltpC).*(1+normrnd(0,sigma/3,length(VoltpC),1)).*exp(1j*+120/180*pi); % % mVoltpA=sparse(ones(busNum,1)); % mVoltpB=sparse(ones(busNum,1)).*exp(1j*-120/180*pi); @@ -224,13 +224,15 @@ IBiW=calIWi(real(mVoltpB),imag(mVoltpB),mphaseBSpotLoadP*sigma,mphaseBSpotLoadQ* ICrW=calIWr(real(mVoltpC),imag(mVoltpC),mphaseCSpotLoadP*sigma,mphaseCSpotLoadQ*sigma ); ICiW=calIWi(real(mVoltpC),imag(mVoltpC),mphaseCSpotLoadP*sigma,mphaseCSpotLoadQ*sigma ); -wI1r=sqrt((IArW+ICrW+IBrW)/9 ); +% wI1r=sqrt((IArW+ICrW+IBrW)/9 ); +wI1r=sqrt(1/9*(IArW+0.5^2.*IBrW+0.8668^2.*IBiW+0.5^2.*ICrW+0.866^2.*ICiW)); wI1r=wI1r(Loadi); -wI1i=sqrt((IAiW+ICiW+IBiW)/9 ); +% wI1i=sqrt((IAiW+ICiW+IBiW)/9 ); +wI1i=sqrt(1/9*(IAiW+0.5^2.*IBiW+0.8668^2.*IBrW+0.5^2.*ICiW+0.866^2.*ICrW)); wI1i=wI1i(Loadi); -wV1r=sqrt(2*( (abs(mVoltpA).*sigma).^2+(abs(mVoltpB).*sigma).^2+(abs(mVoltpC).*sigma).^2 )/9 ); -wV1i=sqrt(2*( (abs(mVoltpA).*sigma).^2+(abs(mVoltpB).*sigma).^2+(abs(mVoltpC).*sigma).^2 )/9 ); +wV1r=sqrt(2*( (abs(mVoltpA).*sigma).^2+(abs(mVoltpB).*sigma).^2+(abs(mVoltpC).*sigma).^2 )/9.*abs(mfV1) ); +wV1i=sqrt(2*( (abs(mVoltpA).*sigma).^2+(abs(mVoltpB).*sigma).^2+(abs(mVoltpC).*sigma).^2 )/9.*abs(mfV1) ); wV1r=wV1r(Loadi); wV1i=wV1i(Loadi); % [ V1r,V1i,I1r,I1i ]=IPMLoop(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,BalI1r,BalI1i,busNum,Loadi,fsY11,Balance,1 ); @@ -260,13 +262,15 @@ I2measurement=mIf2(Loadi); % wI2r=abs(real(I2measurement)).*sigma; % wI2i=abs(imag(I2measurement)).*sigma; -wI2r=sqrt((IArW+ICrW+IBrW)/9 ); +% wI2r=sqrt((IArW+ICrW+IBrW)/9 ); +wI2r=sqrt(1/9*(IArW+0.5^2.*IBrW+0.8668^2.*IBiW+0.5^2.*ICrW+0.866^2.*ICiW)); wI2r=wI2r(Loadi); -wI2i=sqrt((IAiW+ICiW+IBiW)/9 ); +% wI2i=sqrt((IAiW+ICiW+IBiW)/9 ); +wI2i=sqrt(1/9*(IAiW+0.5^2.*IBiW+0.8668^2.*IBrW+0.5^2.*ICiW+0.866^2.*ICrW)); wI2i=wI2i(Loadi); -wV2r=sqrt(2*( (abs(mVoltpA).*sigma).^2+(abs(mVoltpB).*sigma).^2+(abs(mVoltpC).*sigma).^2 )/9 ); -wV2i=sqrt(2*( (abs(mVoltpA).*sigma).^2+(abs(mVoltpB).*sigma).^2+(abs(mVoltpC).*sigma).^2 )/9 ); +wV2r=sqrt(2*( (abs(mVoltpA).*sigma).^2+(abs(mVoltpB).*sigma).^2+(abs(mVoltpC).*sigma).^2 )/9.*abs(mfV2) ); +wV2i=sqrt(2*( (abs(mVoltpA).*sigma).^2+(abs(mVoltpB).*sigma).^2+(abs(mVoltpC).*sigma).^2 )/9.*abs(mfV2) ); wV2r=wV2r(Loadi); wV2i=wV2i(Loadi); @@ -298,13 +302,15 @@ I0measurement=mIf0(Loadi); % wI0r=abs(real(I0measurement)).*sigma; % wI0i=abs(imag(I0measurement)).*sigma; -wI0r=sqrt((IArW+ICrW+IBrW)/9 ); +% wI0r=sqrt((IArW+ICrW+IBrW)/9 ); +wI0r=sqrt(1/9*(IArW+IBrW+ICrW)); wI0r=wI0r(Loadi); -wI0i=sqrt((IAiW+ICiW+IBiW)/9 ); +% wI0i=sqrt((IAiW+ICiW+IBiW)/9 ); +wI0i=sqrt(1/9*(IAiW+IBiW+ICiW)); wI0i=wI0i(Loadi); -wV0r=sqrt(2*( (abs(mVoltpA).*sigma).^2+(abs(mVoltpB).*sigma).^2+(abs(mVoltpC).*sigma).^2 )/9 ); -wV0i=sqrt(2*( (abs(mVoltpA).*sigma).^2+(abs(mVoltpB).*sigma).^2+(abs(mVoltpC).*sigma).^2 )/9 ); +wV0r=sqrt(2*( (abs(mVoltpA).*sigma).^2+(abs(mVoltpB).*sigma).^2+(abs(mVoltpC).*sigma).^2 )/9.*abs(mfV0) ); +wV0i=sqrt(2*( (abs(mVoltpA).*sigma).^2+(abs(mVoltpB).*sigma).^2+(abs(mVoltpC).*sigma).^2 )/9.*abs(mfV0) ); wV0r=wV0r(Loadi); wV0i=wV0i(Loadi); % matlabpool local 3 @@ -314,8 +320,8 @@ noLoadi=[1,8,10,12]; % noLoadi=[1:3:max(Loadi)]; % noLoadi=[1,19]; % noLoadi=[4,3,2,5,6,7,11,14,22]; -% noLoadi=Loadi; -noLoadi=[1]; +noLoadi=Loadi; +% noLoadi=[1]; wI0r(ismember(Loadi,noLoadi))=1e100000; wI0i(ismember(Loadi,noLoadi))=1e100000; wI1r(ismember(Loadi,noLoadi))=1e100000;