From 942120eb70f8884eb3565476d03d3c26cfe8b694 Mon Sep 17 00:00:00 2001
From: facat <dugg@21cn.com>
Date: Tue, 19 Mar 2013 19:41:52 +0800
Subject: [PATCH] =?UTF-8?q?=E4=BF=AE=E6=94=B9=E4=BA=86=E6=96=B0=E7=9A=84?=
 =?UTF-8?q?=E8=B4=9F=E8=8D=B7=E6=96=87=E4=BB=B6?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: facat <dugg@21cn.com>
---
 0PD0.mat                          | Bin 671 -> 671 bytes
 0QD0.mat                          | Bin 672 -> 668 bytes
 12PD0.mat                         | Bin 668 -> 667 bytes
 12QD0.mat                         | Bin 674 -> 669 bytes
 16PD0.mat                         | Bin 670 -> 668 bytes
 16QD0.mat                         | Bin 662 -> 672 bytes
 4PD0.mat                          | Bin 671 -> 670 bytes
 4QD0.mat                          | Bin 670 -> 666 bytes
 8PD0.mat                          | Bin 671 -> 664 bytes
 8QD0.mat                          | Bin 670 -> 668 bytes
 @ForThesis/MaxBranchDeviation.asv |  43 +++++++++++++++++++++++++-----
 @ForThesis/MaxBranchDeviation.m   |   2 +-
 OPF.m                             |   6 ++---
 13 files changed, 41 insertions(+), 10 deletions(-)

diff --git a/0PD0.mat b/0PD0.mat
index c78db611a77fabf43531296084c3723ceeef1bf6..00bfc6c4a94ba360086f1310bd3445a5223a1444 100644
GIT binary patch
delta 576
zcmbQwI-hldkwkcEih^%qk%EDlf{~Gxk%5(giGq=Vq4C5(?THEO6Ke$OqdhJ&im-k7
zKV^;-@9Wf_q&1U#W36uA3A*)7Q_mzlGh~;^?5z$L^<1WSK0RC%a;oQA*Ax!367$k+
zQ==nor?`p<a~{>v%9v>DQY^+Qxw>yU-ye%Pp9?JH+Y@K}vS{;MFwbM!QUjl7HSx3j
z1+O(#EQ^`Ecf)??N1YQoJfu9;CF@;$%nZ8CH+8h6Ozw5@N>OP)(b<%vRGDHNC^&&B
z$Ya|CrA<>aPMQU|l}(tNV!l#gs*C%#30gtE-zG_I>VM&>yE)>2v7xuU+S%vJ8%i&(
zUf4Nv_J${|!Cqa=H}1cB@_41Z>P4s5tFJg_u=?Nfe*696&mD;sf4?e(-TfW*`NJIE
zdT0A3J2x}Qn@lu6o@#!_LiU}|%bk2K^UOboW$n!flTG1W!khT)#M)&UE99T|8m?dU
zQcv$|*J}OsUlZOiuHN=|gRJPq`4d^Z8H5sl-EdvW9d5CVO)y<xSE=ir3$Gde9s6}-
z{_;0FYTe$Rkb2jskp8aW{-x%pr`q`ABmMULjhlYI{?a|(;_k?Q*{`I2@o!vqm(RRU
z)L!c3`_0Oh-|x1#J#XA<=O=jY-Jz_<+}na~b#7if*>3$==F*^cgB_<kCri#`+<&6&
zzW9_dpWhpmQjE%Fzb5+SX=I!59FPsz|Hne|z;e^a%Rb#y2s#|9{OkK(iS=7f_OR-i
os2!^)nkQVdZ=tDp<U<$f@H2JS-X55*uYSwTTlEKfFw-6n08SAQrT_o{

delta 576
zcmbQwI-hldk%Vu4o`P>;k%EDVf`NsViG`Jsp@NZtq4C5(?THEO6Ke$Ob3KF$1lT^*
z&kS=pJ%K~BGvVx<Y~Sp)s%0%*XQrklSw+uQR$Lx=Q~K-DY0|;$zHGPKr)ZR3ei(J`
z`q3?Gd7X?Fi=5uHNytfYhG5`fuBw`s4~`dpf4HOlu!MH}$q1fbQXPG%IoF-y=dC}~
z<r?@VV_8PweeS>m%`P`MJrA^}2-I&_D(!TGSIMI><OXxKBhPOKrXY{z5(Ntu$qB4C
z6__n$8f{LrGikIu4q!RDph3uj(aVBUE5KPZde#4(H$L@ei+@s1IM<ris&>I^RpIx;
zfA=V`wh1qDI>>lRr0BKqzTKJKXH3pqsrxl^e%Y^IH_JAipUv*B<kmIYw`!>v&jjcC
z<d?HImtAQ5!8^C>m03`H%vz82ua+Ap`G$8}ddRcotDfGx`{mcvy2{sm0<E*%>b}Ij
zT=@1-b=r4hYu7(Y@p9LF^-`xUeOGHN%FX9(-_N%4|K$rk8!kBg;rjXc_9wBTr(Tcu
zY(9Qs_rF%(U(XnqTKTTB+}nQ4W!|KywR8FUc+KiR?7P)b{rvXjwQs*Ra(<j%+__bD
zN&mIho?G6HpHq*WSj8FJ5PV$rYna)k=+(tuFLM{Y`jGbbzNh%q$lx=E-{(#V6S@E4
z$OGlRB_(t1j(hH%^kB=Ma`U^j=I5Rm=}YP?%l*vCXjPj0cyWwBe`wrolg0mEC?@v)
n^*@%({!^_q(p21x{cg|OEB}m2Caz0*@w6oQ#(#b-&%@jRff*F|

diff --git a/0QD0.mat b/0QD0.mat
index 19bc6502363257634c8a5801340efb43b3116f6f..733cc7343811124660c410dacc864a2b1b3cb03b 100644
GIT binary patch
delta 581
zcmZ3$I)`<FkwkcEih^%qk%EDlf{~Gxk%5(grGk-xq4C5(?THEO6Kgo?MVJ^ED(0A;
zkM_7+D8l|>zFu<kq?xfC7cV+b5A~aSHp;MX^I93fji*YU?Quxl*yEsFEV}9Fw{pGd
zXYS^hw5(X5q@9@5aZF-w&$d)SqlqrvhvL?>7|pKTTk_!W!#VATe=ajVc5bQBHmSM%
znNG)w#DAZ<!=ve4;j{Ww{H^bd{M+gemfc`Xy1}_S!904OVx!H8cBTnVif<StMGnSI
zILsKp&?v;$B*oXv)^{M#f-TB|F|~l})Iz3`0M6wLxxNJOng*Pn_PSng+o#0|hW1Qd
z_O*N|;<Ha!T^HK%<(PBW()h2lEKi1SFFT|A?PBe3;a$SK6O2ATdGcjL|B_9!sxo(X
z*6+F5e@*vJY0#Pl-@V(`_6h}mJ}X@*<60H;vr>rNw&Yg3@IwaHMH?R-`;mL;<(7w&
zWb|9LuGRQU<e2Pxdi9UP#F-WPyQ|_{_LdZG@V5K;M$GnIfBM|L&llDz7cYIkYWK?9
zrEi_uuK1jvysC9QvrKEcu{leu*ee@jJF`=pFZsNv|Fy*-d7)bWukiSK?yJ|Pw*Ptl
z_KU;*uK1v|rBfbWyTHD)_uFiqKMVKEWG;96)6Auu`tyFh@wM%fKW1(Coc`WdsB4no
zbw%S`E;a8;L-CVpQqjt?ep$QUbeR5qsc^FArhs_mimWf?=ccaBE4t(1pKbAOf&8mG
uPOICL_r80(CgFsP6R+LZwCPQY7o_jqP}HP6G5XiVPaz)v8F+5*%>e+AJr=_N

delta 585
zcmbQkx`1_pk%Vu4o`P>;k%EDVf`NsViG`Jsg@TcRq4C5(?THEO6Kgo?C72i(D(0A;
zkM@|%D8l~XytjDt5rZ>r5nYkdQRl9P=^azlSv&3JD$BXK>sItMl|=JDY}q1t|JseT
z)iceNn0k#Smnw>$Tr*)(%Tu+=g#kJK-CL)-s6^E3%l;^|c>b{9pR`2VPLZ<)t1n(Z
zq@{QGP1GCh&!rtvLWec@IH!MGuTamn?J(biW`;5+wiKh&=4H-e9E#E|8`!@!NvLe#
zFM7amqTzwUlS7Fo+7~OdGik7}&S6p&Y3A};AaW&uE!}}RN`qZ?bIku@Lv8zGy3gDX
zY;vD|`LmVJk|>X9Cz?uA=FR%D+{EhSgI8v3to)w8zb{EvyW+Qb{+E4aeXWm<=(pL|
z_so8g6Z*L$yL+}&y4St4%QuMBe7iS8XWPW8$vIo}A1(U0eQVXc?TtyM`|nqjlszi*
zt=+P075kF#OJ0-vcX~#>wYhjR*7x7+*Vn$>I3rU3C%U&Kb!GmWKa=ZLEu6OK@8;s`
zvpg36GL#iHns|PW?WC)hT_lwPpIpCeEt0Nb;&W$3y#z-{#r*8KW~Zb}(zaD>5edi<
zJ|ed9i}9Ry+$vET;dVD~<(E!W{<r&s_VlH_;??HS%OWSOGrzxgR=odN5znW0woUxm
z(RTTE;?gY#ncioAw>{;#s#|_lmU42EbLxwI-BuORhrfkc_o_`j|MCBv#hok<PoMb1
zyy|+brrf71mo`jLkpA|0?UipkgV(+<3TWE&=}^-huRX_?Pu+fO(M$Pk_gyCeOMM)1

diff --git a/12PD0.mat b/12PD0.mat
index 85525f4c728e5cd84eb981c8e487c878db363530..65f1f295e266fe12413c38b827d6d4681a589fd3 100644
GIT binary patch
delta 580
zcmbQkI-7NZkwkcEih^%qk%EDlf{~Gxk&%_Lg@TcRq4C5(?THEO6Kgo?g_#%_D(0A;
zkM_7+D8l|>{^UQerCz;w)a){0`csp+e#WamraoICm~9n2{r|E{-J3m}G{PRpRe9}N
z+vDY)FvDuX*0sG4mu>AnCa1&k#z}4467EAQG-luE?VHY5V=?D*fyMohV>1E|L<oIN
z*`zmjqwMM0T(gNSqMd7Yd2QSDeS%B<ByW`z>yyesll-}Eyy|rda!FY-V~T*KrfO=6
zx~95l3V)!0fa*>6kQDpC$;~d#Tb8t4nPe5@FV!(mMzOa=-MzEXWMbN-Lt(G$^|pOW
zK2=o5C2_jDX743?hVljJKi+;9xRs;(Jn6&12dC%$xGmDM&o+2j<O~Ouy>YcKi#Wdj
ztPecdncVy%LGtvyT*=mR*_T$j=B3Pj<YD#JC2f`Etw5_Lp5+I=t@&fPYkpnjgr~9x
z`!$cZ-oNV@GPifBjK{YG#Wfn+YPUE(y}jq<2GfQY>vGL|jPqi)&sbh>$Kzbtaryfc
zRzvkH(a+)S1~b3KNmV>?C|Tx_+;#h+q29)McV4^HuYFK<y`ouGa_*<KQLlXqp4(o(
z6}NlBVJU5G<+7ZM+`p$SkvsGKtKO>TcQ?eF><`#-zwzDz?_i5N1yge>9@RfRk`QvD
zezn{o>x1i6|DBPrNPc~+%hKMZbkZ8B#p0a*9?uP)-}K^Ce9g@oEtY-t3PvlxJdBel
tTzI0DX>G-}@Hb4-<=u}<io}=Quc@<>E_wLwm;EXyJ@)?$T+4Rd0{{?#8Iu42

delta 581
zcmbQuI)`<Fk%Vu4o`P>;k%EDVf`NsVsiBpLp@NZtq4C5(?THEO6Kgo?MVJ^ED(0A;
zul5LL6k+@D`b@zYohGS=YaHKPjmZt$ZI|Z0*)Y}Z<dvM+rYy<_7vx+_5N#3CSu%B|
zD6i&xZ56ZS)6+Ia&enEXq_AR|A=jw`2OPAFB5UoeK0GcsC*S@#)rhBf%5ekji<cF1
zc@M7;`}OHf!xYz*Xa3o~I{MdONmIQ@iUK=p4x{Q518?(fhs*++*@_aa-!#l!z$*4>
zp+Hx|0v0bvel`tu<u|Rh0i0?p+N3NPwiXDSNnq5za!}EN$1H8-zvoHQ>YwO8aX;|>
z>zgZGEA3|Oxcu_;>W+)6^VjbdwY(eh#>04a*WxXedz?OnNh^N;Jek$-u|Q>T-mKWi
z^?OB1UrPTgnXA^DcX&xe3cvq;pSYb5_qhHFn3;Ib*nxM`=QVfUEq})N@NdI?*VPtx
zXJnf=39ouRv43)&P@(ba&&lPP%Qo-QnNZ_Z+4gNsb}pk;>TQebS`|(UtJWoddG4$g
z6~1$2{L@9=+gfIC{qdsx*Tv>lPYmY1o;mf^4Fk>7^-H<Z^X&iT`o1`FE$!Cbln|D$
z$|=7;S-+@!zF1s#Y40v8`?E)N^L_5*&U*U#@weFKDL-PKSgg+F-mkx<?uY5yr4ii<
z_bNUbZM0c_rBKBoFmqd2Wm*qc*6bgLwmMGg=kat2l{vbZW!_iG^Dia`Ed2S3_u@w3
uTHRQoz6i<l_hO1=`l~iIH$4__x?)!mIxX{oht2%JPZo2R)=xaN%>e)<U>8mR

diff --git a/12QD0.mat b/12QD0.mat
index fe326ed76cbbd9f4cfa44c5e3d934dc69bf1ebca..371b4f27ac6a0f7cfda6a05918b39619cfee1d45 100644
GIT binary patch
delta 582
zcmZ3)I+t~VkwkcEih^%qk%EDlf{~Gxk&%^&v4W9-q4C5(?THEO6Kgo?MVS~FD(0A;
zkM@|%D8l|hzjUGTR*j-FmfWXq>ugUwclK6%=I5ZZ5r?1o%;lA_+i`R{hZFxq*R{8r
zi+A<s2%B$@Oh|Qp(--ifS%0I~%IQ@a%%V}YE?YUQ<E%bBE;x7i;h)QrZJGWFGtOR_
z%$cu!_)YX%>0<MnThu0A%wCfHo3nnBN0F!8Np;pG6FfWwH`gWGYpMk)DHjzrOqtNU
zWQs<|q$!gRvn-h;qZk$Bn|0D~k+OQnJP8k5OV^$yjWaxyEq!Ns7-uS`T{^Yq`u~_S
z$NDu_|I|!4llRgvr*J{Uf1ikbB2K#XN9{zuKYrY;9ri9pVAa19d|UtNr5xI}Y?krj
z`Wv5DYfZ@C+x2da(#1JO49k~nV3Oi{sJQWv&EsSJ6TYQ<xvely!q{e2PvMt)iNbC3
zs|0PoT*{xv{rJ||oAZ`bGaJ<0|7UsnSxf1^o$(*+I3^$ATY0*9;!HnBzfg{Q*S~o?
z>G52%%JiI7^ymfC`>Jm_@}H(@M|&-c>B#TiA659JUU|l*82hrDtcOcpJ!)^@jkj;T
z{ltzh{#<!l*@r4cq4+i4_f8%Z6!tIKz`Vgk=6vnvA53!ZCjI~T=V^fPqRv|>wFajB
zFIWBaTjBIx?B=(`BQ4iH=9Ds+y=Csqztw)J<I}(Q7c{2TX>Afa<upUSMmy=}b_2<p
vD^-uzNKae0<?c(_>bomGhiyn*na9D^Q<Jjm!qpYK=ZbZ%{?9mb+WO}J3d$Kh

delta 587
zcmbQsx`=gxk%Vu4o`P>;k%EDVf`NsVsiBpLxq^{_q4C5(?THEO6Kgo?rI;8PD(0A;
zkM_9CD8lxke%910VN=A)7ORxZycsj|&7|d1#nfh6nG|ns{U+e4ZD{2CVcNlf!0RQ^
zR?m~P&pr#$aVt@bxOO_u>}_JeW;vsnbrV;I-F*7_Qlor(u^iv?WMiIUy_Upj3)eSn
zFy+}?y-xpJc;d-D0czG$t@q0}2Gn;nMm#v+*U=m+a9-z*M}zHwR+k3{<8H8fZjhhU
z$nRvqz>;LalUX2gNs+nlFy{$2&L>WCb6RRFSXmWVB}EQpPB<)?w(#HcB(Zv)`KS9E
z_HODrCb_c9P4{8)d0po%FN1#etF50nXHBo#?+os_Tc4d?x8vQU`{F6Nb$2|UaUQg)
z|2tnQNXzA--ecKE;r-d%XVzZ%zk0<!j|&ca#o4qY4wgk7_4D}F_1iNuHZD@;jbiyV
zk?4Ou&g;F8-QFdCXHh=m#%s@Vmhr#sGiN$&V0&}I^Z&JxzVEU;Q*ZMb_o{5OdKI_I
z!{EcqmlHNk@wt|~t8)t9FP*Z#rZ;A3{J1C<yF04hJASsW&ADq<yDr9FyS#mE*P^t=
zJo(!*I>qPgwbZ<Hac@yvXwAA=OTHdmf5|3`i^EcN=0x5lmpOue_RL%MGw66wJ@3kz
z!;<~$<}N=JR2_e5y4a_uH<;Vm)?fY=v8i<L>qRryZ0eY*-?L<g#K~g$J-&i$f$qOT
zbXzNHlDlU3P4ZSUy8qn5^vi{#v9DzVCVu?=_vHNEDk~wS&d)jX{eH4*EI9lc09z>*
A!2kdN

diff --git a/16PD0.mat b/16PD0.mat
index 37a72207cc6e1a9fc5b323e5805d6eeccc0973f1..e4e2f654ec0466a9a0ca963b9b49e888304dc7c7 100644
GIT binary patch
delta 581
zcmbQoI)`<FkwkcEih^%qk%EDlf{~Gxk+GGbxq^{_q4C5(?THEO6Kgo?MVJ^ED(0A;
zkM)?$D8lyOeDBQWy7jkQGo5EF-Mn_|7NhTKQx_IrPW<)MY<7tvYwD5=<C8}^HcG5b
z6?e{!6p7A#x?SecHn(+>+`9wxm0l<rO8mXa#J=uAcb~m%g~goWLW}olZHebTN~le)
z5);fxopEUDKD*hT$F%H%i(~ff;;$F5n>fK{)9g2jKbJIk%-kwn-Pv>Eq}?Kg-<^|W
z6kM0gkWo@yGNob4#N;O|mY!N2vobtPFHVdKa(~6??KhcikvdzFO79aM!-;BXYyUk@
zx>hf;?#c27_T9&%*2>RM$ej4Qw#16-W{;A|nqCJ5$<<fwdlsBJ6Q+08@4(N8QNHKD
z)teh|6yMt>vODoqQOx#@>lL1IK6Bm2kSD*fYln56%ID>F@%JQJKD$3$_fEEO#-`^b
zwucvYe6Kry`N@$DYn{&oo3d<J@VdbImQ8|i_qScYAFr6SEV{UxX-1Bml}&+S*G$gY
z(`uz$w>_<Vm2=wI=GWpgU-x`DoMWi=C`mr1c-4>IdI@&d`rQ?uCh>DG(my%nAgAvK
zrT5PQ7<PUC{^tD;{STMlZuC3bcc3P#?)kD!f;MJVZ|Al-F8Qr+?uW>MqS<D$+TRZ^
zJI~{5Ua6Q@8>{)6|8ev#Ta^$ug^9;c)vs)RQEWMD@5aK8N#)L`!xve^f7<?C%bM#d
uTS?n)i#wClbXc8a_!f4(3fVFH@U-<yvyL*`$7mJ*ebB)2pFyW!?_mHQF%wz<

delta 583
zcmbQkI*)aNk%Vu4o`P>;k%EDVf`NsVsgaeTp@NZtq4C5(?THEO6Kgo?#h4fvD(0A;
zkM_7+D8l|>zVx+evD2hPT4qgHy7N|%rON9RH8t;Bb9U}ovTnjkX8l%!7YjQ!m_-QK
zUdbsjG`-op^XL4sNXw#2QOW{}N!^QP-E@*l76?iH^7T?9e|xbU-}Bw`46mI^&Jo?k
zFv+#2Q~dX-oQ5T;HOp+2B5xHRXs?&K!Dy{0{wYEH$>yihyA%069Js0%GNf$a=W^t8
z-5}4$;K(_}NpgxatHuJMzGki~1x#5E%<2opUM-N5dLUxrBy6(z>96_SYyW7(oat|P
zd#f~dkGq*rl~eoGB7WBG`)B`?D-T{6p8sp2p<mec1#44pc-y>mn5griMQu}Q@y@OF
z+fP_nuH$uW*Z5_;)bZx`rhR6s?!97Ro-$Wzb@Y!XD)m2qAOF?9xJ1wQqHXY*i){Dr
zJh-OyJZ0tnq{;rb{STg=Zem(weAchuq*m}|ruUKS%R6ssK7SC~`Mke7D%PC&yo>Rb
ze^=(;m~+*C|JfMM<?@O*TMU;NNp|M{TJZZ$=3@5xy)zXLf0Z^-^A+FxSxm@}UDY;G
zJT!FGxn=g-ga1VRo4dtFdS}+^<tc{7nvYZO^Gr;X;`Prxz2n8g>c-Y5H~rY>PMZ^`
zb}vF>jcM81XL}cI+*Y~zpMbm7!C8Bb7_u)d=}=zj9v^%q+E;U#<E$-w>)*b}T-JN{
w)wGSQee1;3@1Bz@Dqf<oV*cI<PZ!+VzkglL7K@ilkLk_bufMcDW!=8L069ky4FCWD

diff --git a/16QD0.mat b/16QD0.mat
index d03a88e6ffb12a86d47dcf3564d753359821997c..fe41c6ba5ef90f7ab54a7dd1b6925e08dbdbd43b 100644
GIT binary patch
delta 585
zcmbQnx`1_pkwkcEih^%qk%EDlf{~Gxk+GGLp@NZtq4C5(?THEO6Kgo?C72i(D(0A;
z&-R$iD8lwYf9}G-Vn@%UXu*X`ucgggH*Joa*H%@{yh!8l2Pd*3nuNYCR!a5Vm2+>~
z%$lp)lJdN7^i*>Crq`<|AL{gA;k~x#j81gOl9I~3B@Z4Joa;Yaxs&JcoUWu9OJ8gh
zIGj8~x45=&AA|QZSL4>5i+=9qUsSIs8I&NL(-z{<T$1+iSeQT)hf#ubPP@+y#>ff0
z){ab;4hpS%C$Rs1z;HrQp3hNHq)~{^*-)l&@q#8H3$~>T+8!4)MXvb2bH-8m$Im_~
zE9gHl{=K|pip7j9@rb!keEDOw<hD&dH#h#{)<d3}n_kK!9&!GrRLvg$Ol|6;m(}~~
z_f0#dpC_<G#jNUwM*08s3EZhFh3i;O`F@!s{qN%S(&_g?erDfW_wAB)bW5O1#3xVp
z>6!DlOx!*(H?$^aA^)S&H)p0y<as-<dGGz0B^PJCJ37ZCcvivom^*Km7)@5L3!T(m
zclzAUQ?Gcw=$R<YIAa{RI!tK8y{rYtrBk!ZEN8u_KmF!;*1X+~f0o{yHSwPmS6bn@
zuH!mef^;o|74nxPIPGvhvGKxByOW!LD{)_4+rM_#scV-UmVW&5_H|m9&i-#z6)QUB
zddjY>UdmP_C&9Jb=<OBx`sUSM8z!!+U6y@beOils>(yUz?2+2lJ;`(V9)0Grn%Y^u
zVz-N>+|-<9zuwm?sJpAp$_%)$M*sSZjH7Rk`0{5jd$~1t`u{y0i|gGU98?7W%)1#w

delta 575
zcmZ3$I*oOLk%Vu4o`P>;k%EDVf`NsVsgaeTxq^{_q4C5(?THEO6Kgo?`Is0OD(0A;
zkIlH0D8lxk{8ar}y`&b=#oHG3zRZf=c7BdD|H7?0PG@u1Ub38$utkDbYKr5=HiOAr
zr>C|}>q_EQ*FGw-qm#q<<StHqu2riiMZHtvSfqFDraOC`%(<NwbJDq!b2^SkXnr+1
zG}ZXphFj*(w3npt&d(H?w6(V9gJ-?2prGMO<=-djg3e#dH<`rj;>0&e?UGNFhjym&
zTNkgCCEX1&YDGmYCBmwfj($R-7X@@*s@6U+d8wX$BC@FOMaMiI0cOdR*I(_$Z~q9~
zFl#=Ood@5ve)ji=mQ5||U9$ItnCSX%I{uF@-*1_|+nBd>Pny;trw6kxzq|3i`o|}O
zdS%_3xh_lhJ%0P>@^KILoj?1X#Q#|wnP#*3(9g=%`)d#IIG%2OwE2uqleJW-siyf6
zCtsV7uUUECzcrtvB!A`AgvWRF)ca2EO`e>_f2FY?JM~~!>z04+52wiPxs(-n&y@MG
z-%aJ{d;JSrqh6;@mA<g=r^1Kmgt|+{T85Peew(eTf5f`(@y><DdH?bby<7HsWn8!T
zoZ@|}@{cubj9!_Syj=QNmER^?F^2_tkv4IOu1d<WH>Z6)D7Cg=@0EC`zW>E1HRd>7
zx?&cn)~~c&dS4cI>xDvNna+n&e66yFc{!3YkF4lbsR?LXWGilKcP4D_J&|cUCtTGw
ndwx_;=Ui&Z(JV3RDbqR^@BDY$V&A6ncY8Lo{AakiZ)-mQI{pr1

diff --git a/4PD0.mat b/4PD0.mat
index 7181fd56ee38e6430fd58b2132cd9d31ed008477..3c6017a176ad079d2506aa1f5a3696c2822d0eb9 100644
GIT binary patch
delta 583
zcmbQwI*)aNkwkcEih^%qk%EDlf{~Gxk%5(|fr62Nq4C5(?THEO6Kgo?#h4fvD(0A;
zkM^0wD8l~1f5Vk6McIaDH*V8d+U=AgUi{WxX!5lUv+k<!#>GsJe7#BR)wIuzSzBYY
zFEspMnkeYl@W^?V<GS8!sb|@*b!~m7w5io?tHCn6l07Bc?w>nXX<_|emFMuB2L@-H
zUd>E#)H^J4{>#ru=SyjYGi!}LzS&pbU*M^zspPb2+8X88C)v{8xp4{#8s79!NfFjm
zoqbYl)6^*wi#9bkczCj|R0&Nnzv(QJqJPsr>x5%bhm4xLhi>4c#Y_6GOip%jy>cez
z@A;!={?9O-JfCUZoH)Vf9#(>Vn(N<sl`r1rlIOL=vS8Jj$9*F9hZ2hOkMEsqr?E%y
zUVYhaX{qqKRe|N=hfdo3yv(D1>)coC`*+hW-n8nuE*?3Dcg2*A-ey-?*_9ZViSPfr
zW^-Dgx!?UOUF%-mT@hxz?A&o(2|wPy6`AG=DW0>!pI1(paLn#YNaO?o+c{G{#h%!8
z|68Z+W49M?choPr(DCh4;~Ozi?rG*f?%hmiIvn+vyI#mRYG0phK9_VstU=KMnM;bh
zH7jm@yUCQU^~=~I#bxiG-q`n>6W)a0PYPcq%6c){I=`S)p#Jxz%e(G0l~t?SXC^$H
zU1g@rEwypYgZ?P(+6I}2$+i0rFn)A*6#T!Q_ZH7hi}k0Q{{&{QGk@UfUesiK>rU(Q
vtk5lmhkCe0%h$=ZT${16<m%)vY{K2*|8kf9krA7f_q5#gPFMZ8rAK1{GJ6)P

delta 584
zcmbQoI-hldk%Vu4o`P>;k%EDVf`NsVsezTDse+M#q4C5(?THEO6Kgo?#hDlwD(0A;
z&-S>SD8lyP`jHJzGiR(*n|AT&%{@MsZ|2%uS+jxLICt~e%wrKDaVw_W5ohns@tW*8
zty4$l|3csGIYzTDFH*>zk|x3$D3rviC?d9S-K7^5FAww|wv^|4&Y9jbH}2StqE`!5
zT8>9ts)?WVRVBmG;Pgo${WsMf1~==S9(7JonXuNyMd;?E-jqpBM^vURnV73GH${N+
zq>A+kJtb9z>W-c_%IqFymM(n)mX=<hI+`>*?JWI5Tzt7Y+C(NP6*Y-$KKg5Zx9;B^
zNp|+EbA(d2inzo!ukri1)$C#?V{OKi2kHDmzqsO6%x@K_UHTk*{ekm)*7-+F-W<!Q
z|L5NL;ce-b`%^Q|mzgSlmW+SOx1aUuvwxRL#k580cRyA)yShL}z;;IQjb@AB8&6jS
zn-}I}JnOvu^vHVlb6+I)%eHOZJ=-cJ+|XdJ+T)ki>{sTceO_*8eeCEg!BiJkzph1B
z`=-Y2h^folb&T`!iuivW_xTR<d&bXN%lZAx(qH%LE%=IeeLDAd{_BFLsi~<t79}49
zF8Eqb$eyVqe&^wX%UbWEoun?VlG$Fe`>jR9rD=C3MDr&f)jC&uyKB$FB@ew{NmuOn
zTUg(+M=bx)nTdwh(F^wesLr``qiw-KUYnQmwF*>Y&u%=!yzs?NwSZ?WX=PvIUcUC*
yV_df)?{8A`|Im%{<=ZA6w!C>{`xHUuXRFSCzP|k0;g~1AH(zXeBp*CuuL1yXx*S6Q

diff --git a/4QD0.mat b/4QD0.mat
index 9edccdcc0e7dd526b75642c17551fd1d85282042..5d5ef3a38490893e5e53e3038e320799722faeed 100644
GIT binary patch
delta 579
zcmbQoI*WCJkwkcEih^%qk%EDlf{~Gxk%5(|se+M#q4C5(?THEO6Kgo?g_sx^D(0A;
zkM^0wD8l~XyhKah)ZU2Mwx#+^&bpUp{9ep<gK1LKqM(JMeXe#X-Z5()edBb|y~|Xk
zk-fF8+cU6j=fzuVEVoUIb>rVywNPB$_Q9<jzcqc+&waan(~>{^`Tkjn5}&!+9tJ)a
z&NDphA@=Lj9cM1FLW#SKFCW#-D5yVSwrTDg1?$d^H%duqFWoe6`Yu`0Gi5T@ruH|g
zQ&Z$usu=}2J7js<1`3E;dhiMITKbec(Xe!z)iKXxve+dro+lQTj($SQnaR^$)$67G
zJbtM4k4V6u<je2p$4)=+w(Gh3qU`~>)0b|a^QyUOmf){r%o5MPuF$zuZ6Q2c|5qAk
z{esGfr`~t>_M8*jYx7}={Ok0Ti#}bh=R^cg##<=nd8l=-|M9`<*FG&fD~7YXKQ~3b
zF3mAvnas=I#4+`27Sm*%toWJfmTSv6wy56y{8X{SHZv~bkWx|cg5GLfhK*)Al9f^|
zzy3W=FW(n;v_13Pl&i0rww&5?P-NTT2YtRz_wBAfFnRi>Y!N;iagVyI^_w2%ADwW1
z_cYho-vP;q?ZWwMUY&H0G51lOyUnm$`cKH8Mf>lzb)N89eEq50y37yV=I;-L1wMLs
zB+Yu0Yp3$OzqKn?8@9G;|NeMw>fR4VkK2Q0e4Hc9e$f2sUyU8Tmw)x2Drpz*3=BJH
rYBjO9H?n8GY=o5g@++;+=6qQ9%9ppyy5R8|-{_ymo9j1T*tQn{n~N3z

delta 583
zcmbQmI*)aNk%Vu4o`P>;k%EDVf`NsVsezS&g@TcRq4C5(?THEO6Kgo?#h4fvD(0A;
zkM<A_kYN8XIjOkia$eHS7m5Xz#^(B~eIF@Lo1q)HxpZdKx^_bW;bkXTWM{v1%Dbv|
zw?r)C!eJebv)*4GL~*N0rp^rdthk9?KY3~1!|sdsz1i#j{Iuh1|NPwW*cQW{T-RgH
zF+#~Rvv*p)PCUIUBCO15%f5QXjr9}Qm@GJ6Dzxxe95K6R!4s*#WSHX2_ryug;{k_<
zqI@6Y!NLj7j3Q2g8Ve-67Vw;M<b3|XCZ?f#K@-=C7OfpD(-+Rlo%4V1jc4*Zv;PJb
z6!hmmEL`8XZfS&5O{O`oYC-9&`Fuuv9(R1#w9m`0pDH%x>x5Mg&KBLc<GVe4n)3X5
z@q+l@1y6r_YxlB$63$+q@#5x}ln;0Pj7?kGMV4+Wo3v70WxtWtiQl&t{+L`nS14$C
z%F)`lPabOPGCA8`eV6~?%&IRTpF{QJJa#$H@h|LJETheO=b7i;XL?inRpMlM%Puk|
zE_thV-sSQ{>B`S_;;B4Rziac)ru>!_n$Z?}<y(D5j{7+u<xXqu#UDe}vZmQ@ee-Sm
zr5*mo>-CPh`EHJTe`eC$H~-}CyKpX5zG}GKy*0+vRBF=H-QO}-+J#(y7##ibaPy=$
z`k%`#iO$WvHn(KL=V?>+zu52pK62Ua^Us&;U7h~xxl&qJ+D(o9ODEe*x~D!zN;!A-
w5AE(%bD!|kpE)NNwC&P|R}qPS7foJv>G`B!?Wh+a;`y2f3NPAUmECg_09zpxVgLXD

diff --git a/8PD0.mat b/8PD0.mat
index 486594f1308d1392b608082245614c46dc617d27..5d83ebee8c0018426b60d1cd9b57d3f89e5cfaae 100644
GIT binary patch
delta 577
zcmbQwI)inBkwkcEih^%qk%EDlf{~Gxk)f5Txq^{_q4C5(?THEO6Kgo?1(+BZD(0A;
zkM<BQ5MleU|72U6#YBZ_4T)p7_MDj+y)kL7XWz9`LC5c;<>|4XW=iaE;jcc`q8Quj
z)GBy+j^EszZ4tBAEO+Ge;$AJtrN*0-<0*Af&g=J{k_X2h&S^jVlS?vL_P502Mb8<R
zcWgZMcAxqEra7`Z*DUngw3**YkE8y91cxI_Y5<>T&i6KnRSry73uLk!xQr|q)gAek
zJ2J9tEf5PS;5&Dibw&H>3mj1w9F1QbWN^!9)V{!_X5cIxaU<{V`J><J&fGuO-!QTF
z*REwU=3g`Hx`c`qU;W#AH9xA>^JrCa#5wNZ9IbboX2%5|o$K5w89QT*_oN@UUl^>d
z|F|>MG}iG~P>`9Xz;BO@3U@N{KV6vlvFH7&%l;eHZ+-gyN5eKW&O0LN&lJsF^R6GR
zcJ;itZ1tB9B4_o@+$^ptKY4O%>YS7r|6U!--<@=M%baU}9!+gIVI92E@MPM{cH7Uv
z)2jTU9L=mVJ_a$y+x~s?TVm_J-LLIR?|VD<T|IQZK4*%Q46EDcOAKdYH`(UseZBl<
z`MMj|CI$wr`}0sGbn=v~rP6JGFF%>M@=EbF#a%xh94NZ|&-2yph5N++8aPh8wIo-;
zz$r1kFze}A?x+{rzj;MY*thokw)-+iY%AmMe6VnOw<UkGn(U&F?p6K^3>_{lp1X_h
rDGR6U*Z9xI$yXKxNa%bkt?Dc{P!(8T?IAaB`^@VuU-%Dd@7o0czOxvn

delta 584
zcmbQiI-hldk%Vu4o`P>;k%EDVf`NsVsezS=rGk-xq4C5(?THEO6Kgo?#hDlwD(0A;
zkM_9CD8lxkp7*q%hrgrwk_FF}y$#x(BISNrbG6Uft)DhU_FiKXGh|sFk$pB|s~(SV
zv~P-8w59R(Wo~@SCyTiSUUKGW4QX4_v{7n??TzP4neF+W+sK@oe#D@sPWA8_C;yg~
ztvs8P7wVr^i`dw4ROBg__;<&IgY{((9JmS^*;h_DrgQf}Gee6))7%Z*M-&(}PbkP5
zCFl#JIP;xQlou&r;94Ni^^ik`RaB<=aRAfO0Jf(K*my$@aHg&N_dID@eTx1Q_XAn~
z=e)mcXKp%Y`m$B8yAQAato!G{HM2)CzeV%b$SPD8XY}d5TOQOOeDY+C$1=mpv&U!E
z@1AmJt9`_=_p0%4oc7cobDe&sI!|gbW6*a$ANCX8qL1_+KWo}9^7ZEGEw7_8O3Ho}
z?^@dbYRaaJO;`QGSr*?rp{DLB@_O}~JbPZNAG^Lw)>S3TX?tC8)kyneJzto2;*aOk
z*1WkXmh`<QG<5$b;n&C59v3_cp7gHnWsXGF@|Vu_^)`oA@!k67Wc2a(l<tUME1y4Z
zOZ#LYC2Q<D_sg!JpG)-LGhX#p+pS(~E3&b6gQ#1~6xFP=(mTF(Oxj}A|8w5-im=ln
z<sS?=ck$=TSxTza-b%ea_co(c=WL_>`>Wl=eZQZVI#sQ4E@K_PJWrt4n$KsBZBq!%
w+V|0f%Um*KQIh=YoqSer*F`Tr`*P9FTdQYP=2@Fg<@uWtbozxnU)jM908rEpVgLXD

diff --git a/8QD0.mat b/8QD0.mat
index f8fc07681076cbe2db8feb4fb861668599a018c3..c73bbab8b6ae7efe5a5ecf8f62ab71e00cabe133 100644
GIT binary patch
delta 581
zcmbQoI)`<FkwkcEih^%qk%EDlf{~Gxk&%^wp@NZtq4C5(?THEO6Kgo?MVJ^ED(0A;
zkM)>bD8lw(zjS1mD|hpu)lO$pSKplD{pIe=C22Qi-M%wx_OvaZR<~$$yj7a0_bh8u
zxmm97rK0E?U2d+PA|kbK+)g+ctUa+r?}}(f=9&X*j`#EZv6xd_XmQ@^Sj66jjV@nI
zH)Po+Pdd8K=(a*)=Y-7eCEmZ;6<X?fZm`W%XnQ+>)%wP#-nkpNmrUT@?Z{@S&>j=u
zAoqlsouzeycvqtQn|7msX0|O0*qJ3B3d~p_bS9xuYen1P0+y^BMSss94Xg8sI@8~v
z^m)FXV{3C<o%z?-_t-MmZrrrf@s9KCP0FWk@iST_&(*)6`{m;HN$v9_(>w0fuKW@4
zrGDAeEz3eYpIkS6w64#^V&e28{v}OfTW70nd-Yau;`|z>)mqCneO4B>@6mkA)f<&}
zXI&Y~i}n($j5l8|*_cl^T2XlM?#zbWCdMugS81HTx%bKK1FK)FIIilI?(>`7Y$AQH
zBIwyJkK@VJEt@;#q;~&|sMhXXo_3Sl|FPw~wtDF)_3T^qHeHQRf4uASmGHFRw^x1-
z+7u{xHGRj*3l__S@4Qzy?(pVkKqdQ1w-i%9Pwk)2cNEM!=&jp4zyImAt&<-`p4mTh
zzr&$5ai_0;$-AXgZ<00t&7yF1wPXHIf4}&<bmhC2M*=H<HK#9L#xPgU*R$M3KiL11
v)1Q|H{<*A^d4s=K$VgdT_!bp-$!~s`IiH2K?UMbEFGkd>U9@MIwzC)jm~R_;

delta 583
zcmbQkI*)aNk%Vu4o`P>;k%EDVf`NsVsezTLxq^{_q4C5(?THEO6Kgo?#h4fvD(0A;
z&-D;45McXIe>*1c>ZHIav&6)1b>F;}mZ-ZfC|c?CwP`o!>@_%4az^9QnS(t_ZijM;
zrO)}6^tCvwoU`~4Q}UEViN;o*MV?Q838-=|6unwi^YX#*!tW1vxF3<QKAqUpKaYdc
zvnO}%Tkc|8omq!dUooA$W5*J6qNQHv2FL6K<DAAEvwCIrZH+8THpr$Zw8z|Fv{vLV
zJ1nPggZ*a!hnGfEQh+1p3P)!5g(7{9k}~YFQ<~+p9<a=53C`K~?|BlpJy+Ui<%F{D
zDJGl0as7=9(id3!C7^4+!JP8LE`quDUx~%(%T_A!UD?8N(stP{VPkfUf}7%@>q_eb
zl8PRwdw0#Ri<7-%w!iQ)tG>U*gMtYub6?i&Tr;<8S^V@9m#Z>rFSop$ySe=rf9ZSC
zPj9|t*H2T`&|caXo|U|GYoC^Tjm<+z_THD1q}UFbN$FOk#9zDSD%jm^J@1+;XH)ml
z_}gz)1U=*DRli@d?O^-q!;AHP>8`)9W0ke%=6W~j%UQNl5{r54qoaI<^?xmHo3rZR
z{&kYQiVv0@Keke{a8J*^J-zKE&VHBbZ-1Tnc$Z&YvCg_vdtJNOF1{?;IG4{~>R51o
zys~TWv>$2qS8czt*x&e*V!M9Ys^!t0!ko80?XT7}E1FPmymw~L^m+9ud*h-OUf+B!
tK-V^4-LCZ6a+ynRw8^im4U_M*KGb_s()Y@$6}zl!5)A(```_BS900my3p@Y-

diff --git a/@ForThesis/MaxBranchDeviation.asv b/@ForThesis/MaxBranchDeviation.asv
index aeaa522..ac15e4c 100644
--- a/@ForThesis/MaxBranchDeviation.asv
+++ b/@ForThesis/MaxBranchDeviation.asv
@@ -1,10 +1,41 @@
-function [ output_args ] = MaxBranchDeviation( this, Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Branchg,Branchb,Transfork0,Volt0,UAngel0,Volt,UAngel)
+function [ output_args ] = MaxBranchDeviation( ~, Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0,Volt0,UAngel0,Volt,UAngel,FileName,PD0,QD0)
 %% 最大支路功率偏差
 % 支路功率包括线路和变压器
-%% 
-[dispLineloss0 dispTransloss0]=Lineloss(Linei,Linej,Liner,Linex,Lineb2,Transi,Transj,Transr,Transx,Branchg,Branchb,k0,Volt0,Angle0);
-[dispLineloss dispTransloss]=Lineloss(Linei,Linej,Liner,Linex,Lineb2,Transi,Transj,Transr,Transx,Branchg,Branchb,k0,Volt,Angle);
-real(dispLineloss0 - dispLineloss)./real(dispLineloss0)
-
+%% 先用加了误差的负荷功率计算潮流值
+[Busnum,Balance,PQstandard,Precision,~,~,~,~,~,kmax,~ ,...
+    ~,~,~,~,~,~,~,Pointpoweri,PG,QG,PD,QD,PVi,PVu,~,~,~,~,~,~,~]= openfile2(FileName);
+PD=PD0;
+QD=QD0;
+%% 形成节点导纳矩阵
+[~,Y,r,c,Angle] = admmatrix(Busnum,Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,...
+    Transforx,Transfork0,Branchi,Branchg,Branchb);
+[P0,Q0,U,Uangle] = Initial(PG,PD,PQstandard,Pointpoweri,QG,QD,Busnum);                         %求功率不平衡量
+%disp('迭代次数i  最大不平衡量');
+%% 循环体计算
+for i = 0:kmax
+    [Jacob,PQ,U,Uangle] = jacobian(Busnum,Balance,PVi,PVu,U,Uangle,Y,Angle,P0,Q0,r,c);  %形成雅克比矩阵
+   % disp('第一次雅克比');
+    m = max(abs(PQ));
+    m=full(m);
+    %fprintf('       %u         %.8f \n',i,m);
+    if m > Precision                                                                                                                                 %判断不平衡量是否满足精度要求
+        [Uangle,U] = solvefun(Busnum,Jacob,PQ,Uangle,U);                                                           %求解修正方程，更新电压变量
+    else
+        disp(['收敛，迭代次数为',num2str(i),'次']);
+        break                                                                                                                                             %若满足精度要求，则计算收敛
+    end
+end
+
+Volt0=U;
+UAngel0=Uangle;
+
+[dispLineloss0 dispTransloss0]=Lineloss(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0,Volt0,UAngel0);
+[dispLineloss dispTransloss]=Lineloss(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0,Volt,UAngel);
+t1=(dispLineloss0(:,3) - dispLineloss(:,3))./dispLineloss0(:,3);
+t2=(dispTransloss0(:,3) - dispTransloss(:,3))./dispTransloss0(:,3);
+t11=dispLineloss0(:,3)>1e-5;% 太小的值不计算
+t22=dispTransloss0(:,3)>1e-5;% 太小的值不计算
+t3=abs([t1(t11);t2(t22)]);
+output_args=max(t3(t3~=Inf));
 end
 
diff --git a/@ForThesis/MaxBranchDeviation.m b/@ForThesis/MaxBranchDeviation.m
index fadccac..ac15e4c 100644
--- a/@ForThesis/MaxBranchDeviation.m
+++ b/@ForThesis/MaxBranchDeviation.m
@@ -3,7 +3,7 @@ function [ output_args ] = MaxBranchDeviation( ~, Linei,Linej,Liner,Linex,Lineb,
 % 支路功率包括线路和变压器
 %% 先用加了误差的负荷功率计算潮流值
 [Busnum,Balance,PQstandard,Precision,~,~,~,~,~,kmax,~ ,...
-    ~,~,~,~,~,~,~,Pointpoweri,PG,QG,~,~,PVi,PVu,~,~,~,~,~,~,~]= openfile2(FileName);
+    ~,~,~,~,~,~,~,Pointpoweri,PG,QG,PD,QD,PVi,PVu,~,~,~,~,~,~,~]= openfile2(FileName);
 PD=PD0;
 QD=QD0;
 %% 形成节点导纳矩阵
diff --git a/OPF.m b/OPF.m
index b59f5ce..df3f0d6 100644
--- a/OPF.m
+++ b/OPF.m
@@ -7,7 +7,7 @@ clear
 thesis=ForThesis(1,45);
 for II=1:1
     [kmax,Precision,UAngel,Volt,Busnum,PVi,PVu,Balance,Y,Angle,P0,Q0,r,c,GB,Linei,Linej,Transfori,Transforj,GenU,GenL,GenC,PG,QG,PD,QD,CenterA,PGi,PVQU,PVQL,Liner,Linex,Lineb,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0]= ...
-        pf('E:\算例\柳金Ⅰ926_21671693_2012-09-06\newFIle20.txt');
+        pf('E:\算例\柳金Ⅰ926_21671693_2012-09-06\newFIle4 - 计算用.txt');
     %pf('D:\Project\青秀降损项目\最小化潮流\最小潮流算例\原始\津头站津视922（3-1）_0.5_120%.txt');
     %pf('D:\Project\最小化潮流\最小潮流算例\仙海919.txt');
     %pf('c:/file31.txt'); 
@@ -55,8 +55,8 @@ for II=1:1
     %% 加误差
     PD0(Loadi)=PD0(Loadi).*(1+normrnd(0,0.05,length(Loadi),1));
     QD0(Loadi)=QD0(Loadi).*(1+normrnd(0,0.05,length(Loadi),1));
-    load('200PD0.mat');
-    load('200QD0.mat');
+    load('20PD0.mat');
+    load('20QD0.mat');
 %     PG0(PGi)=PG0(PGi).*(1+normrnd(0,0.01,length(PGi),1));
 %     QG0(PVi)=QG0(PVi).*(1+normrnd(0,0.01,length(PVi),1));
     %% 读变压器容量