From 6fdf435f666d89571e1765988eb501ed38085f64 Mon Sep 17 00:00:00 2001 From: facat Date: Thu, 4 Apr 2013 22:05:27 +0800 Subject: [PATCH] =?UTF-8?q?=E7=94=A8matlab=E6=9D=A5=E7=94=BB=E9=82=A3?= =?UTF-8?q?=E5=87=A0=E4=B8=AA=E5=9B=BE=E3=80=82?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: facat --- 20PD0.mat | Bin 672 -> 671 bytes 20QD0.mat | Bin 670 -> 671 bytes FormLz.asv | 24 +++++++++++++++--------- JSMJZM.asv | 12 ++++++++++++ JSMJZM.m | 18 ++++++++++++++++++ JSNZN.asv | 14 ++++++++++++++ JSNZN.m | 18 ++++++++++++++++++ OPF.m | 4 +++- OPF_Init.asv | 4 ++-- OPF_Init.m | 8 ++++---- case1V.mat | Bin 1274 -> 1284 bytes case2V.mat | Bin 1276 -> 1284 bytes case3V.mat | Bin 1250 -> 1284 bytes caseM.mat | Bin 1053 -> 1061 bytes caseR.mat | Bin 1016 -> 1057 bytes plotVolt.m | 20 +++++++++++++++----- 16 files changed, 101 insertions(+), 21 deletions(-) create mode 100644 JSMJZM.asv create mode 100644 JSMJZM.m create mode 100644 JSNZN.asv create mode 100644 JSNZN.m diff --git a/20PD0.mat b/20PD0.mat index abd9e7b5b5adc77140e5bcc1ceb593da6b750fd4..f8e0b19279e388edd7c7386e21bc59c10b94b694 100644 GIT binary patch delta 584 zcmZ3$I-hldkwkD}iGpupk%F;-f}y#Uv7wcLp@NZtq4C5(?THEO6Kgo?#hDlwD(0A; zkM_9CD8lw(|Ju`=Vp1=4%#fIvcJ`%Fbn=$|9PR0MQ+;F4aRtq8_E>m(_VkTAr!6zm znKW_ok{tfEb2hd=oODT#XPeks&K zxo9fSWYf#fpLeSyigjD+oSu32yP#OTr@CZ{@JW@~C&dn>Kg;9rFgvN}v}x*-Nx3Th zIqFJEYC@YPH@s2W-8tutDp!yP-vs9;9GXgAMLjD#3^f&(y7=}@Hd>_k^a+b``2S+V z&-QwAPx>>qo(}oDyl$>b;pu{ZODi}gZkcuA?VaTovnEZSDD*vf!{>x+x3&bj{L=Ez z-cnz;Ljn1eu>J=iA@=VXlYP$GE!BDjh2} z*k9OCC}6%N-pnG~*}$*os9gH(1JZM!RQ#VY;q%UW>rU6SHePz2eeQ<5(u(-NwR=Px z)@T1U_IH)I#2+{>@^07PSuWr67c$GezQ1{0iT1LoGW8QLU1Qam{wZ3;Ok%eFhQ*K6 zBA%2!yjE=cVqMAIkQ!q}b2+|czBJ zNuyUJIVbMdww*T*xy7(B)xG629q-wOOHD%IT7r$*2 zbvDh_IN`UcyFue*-lldD0ntEVUQLy#AfJ*YEiXS9dh&S$zNdoq5q(XFncgHJsD@)%ay_L}aJv;Wh3NB`t9J01PA5BOI7O1xgeH^;QC zzTnuaC#Jgl`j1@Om^xQ}cjE34L%SKz8rS{%v?^;g_p9goWovee>j$~|hAa0AtUnVy zarM_nA#XQpy}mAR=H~i*_xU>}irZ}ui+FS3#Pb`OnJOJ0{om)AnQu6;?daEa4gy6R zkIPuy>yECskBN+U(D3kiz7HLr?2*53=4Y`X~W$Mj~-1}71cH9#^GU-t5*5KFj xPCu2~bj`ExRwm-K1Z&ODT~B3WtSLl?DwjDha^A^@dX7-9eb diff --git a/20QD0.mat b/20QD0.mat index e294606a4eef4e798faec8d15a3d50dde2983a18..0296b5f27a7f0d1a79b06acb1ed2260ec6d0f422 100644 GIT binary patch delta 584 zcmbQoI-hldkwkD}iGpupk%F;-f}y#Uv7wcLp@NZtq4C5(?THEO6Kgo?#hDlwD(0A; zkM_9CD8lw(|I1wmG?aB5b+;I8zml1Hd(%hFv$0!Ly<;-Q|Hd%Oojg&NJ=HuO#B`(rWZbAiQux3)yRHHkB}Ub!wD zrgC_R&ey9qJq}H4HTEjJ_357TA2*LB;URXGD@mNliWpA=4*H|YN}`zHO*0RT{2^ilB$cZ-{hd8w#b$Lch0yf zZ?5xnd4pc$8vXYhB!b?}xnc0=pki|CJLczxHjB4^na&fdF)jMJ{ii(@H>A1x?v*{R ztq%-dd1U2*q=awhj%?xHe{1RHh>Sx|&+6xWNjH78>)*3&AK#j7xV}_mS4-TUzMo=O zuYP0@yLXwV@@Dah#kY3+ESG*@@N4(sOIv63b28F_zbyY1nXJd2t$$i%qj7!Z zKe_xbO51PBeAoP^v?IBR@tP0!FRvn&$1~I3JXpfU;#F_F=GIAFW#+FIYM=K1owW4m zqa{y56+UT2ESy=r%u&{-?ffc-a&5Ne{8QG+CKJNcILvQqCuFa9m3lo(ApK_M<>>Q{ z64{^L{tN7!e)wFheSLY7#YfN5J+~S&d@oG54R5t;-KbY{s`ve0#mdhYjQ6Y-mEUWh y;^=F-|FOTL&ZGJ!WkF!R~0a^E&_qpBqX5 delta 583 zcmbQwI*)aNkwkcEih^%qk%EDlf}y#Ufq|8=nSzmlq4C5(?THEO6Kgo?#h4fvD(0A; zkM)>bD8lxkKJswRwHvLC4;QK2o^$qWf9lq+8Iq|f+O@QO z>rAKaE!K=EN%U-smQ=jjdN3#=Ft9hqF>&G~$&KAx9=o&G$(-A1F(+Lz`OOKZV=Px5 zpG*==ex3g7Pt7cy8KTF(w4S>8x3*r${n805B~_zIeos`I&)yA~;Il|YYm)0GPrj3? zS(_SjRAzThYY>@`n__*_v1Ccd8>OQz{$3~IHuXz*Y8p=byQF7EXX}*7%uAZiUi}}F z^Y6H-@lVYMQl%O{vwCICpH54toy%|g@8YMc$IL>GUbp4C{%d0^zutMfTaDAx+urf& z*MC*|>inPk4@1 zFBiRYV)tq>DNFf#+rn-C1>RtZjnDhEHqUxqjj?)Hc5T7!YN;S0>4g&{ZLjg^yG%U5 z+jyb-dHj=_a5vkS*yPN2?0z4Xz48y1V`uu)+BLC0lby|6*4{jI2TQBnag{UM6Sc*K zz5NP)O#Pp=YwZs%wd1kO8{hj@h)Vr-;pW5 zJ}%_6cH46Bl0)=0.8*PDL(PDL>0); -PDL(PDL<0)=1.2*PDL(PDL<0); -PDL(PDL==0)=-.2; -t1=([PL',QL',PDL',VoltL])'; +PDL(PDL>0)=0.700*PDL(PDL>0); +PDL(PDL<0)=1.300*PDL(PDL<0); +PDL(PDL==0)=-0.400; +%PDL=-10*ones(length(Loadi),1); +QDL=QD0(Loadi); +QDL(QDL>0)=0.700*QDL(QDL>0); +QDL(QDL<0)=1.300*QDL(QDL<0); +QDL(QDL==0)=-0.400; +% QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF; +t1=([PDL',QDL',VoltL])'; t2=Mat_G-Init_L'-t1; Lz=t2; diff --git a/JSMJZM.asv b/JSMJZM.asv new file mode 100644 index 0000000..86be91d --- /dev/null +++ b/JSMJZM.asv @@ -0,0 +1,12 @@ +JSM=[0.0575 0.051 0.0029]; +JZM=[0.003 0.02629 0.0582]; +plot(1:3,JSM,'k'); +axis([0 4 0 0.07]) +text(1,JSM(1)+0.001,'0.0575') +text(2,JSM(2)+0.001,'0.051') +text(3,JSM(3)+0.001,'0.0029') +hold on +plot(1:3,JZM,'k--'); +text(1,JZM(1)+0.001,'0.003') +text(2,JZM(2)+0.001,'0.02629') +text(3,JZM(3)+0.001,'0.0029') \ No newline at end of file diff --git a/JSMJZM.m b/JSMJZM.m new file mode 100644 index 0000000..936502b --- /dev/null +++ b/JSMJZM.m @@ -0,0 +1,18 @@ +JSM=[0.0575 0.051 0.0029]; +JZM=[0.0003 0.02629 0.0582]; +for I=0.01:0.01:0.06 + line([0.5 3.5],[I I],'Color',[220 220 220]/255) +end +hold on +plot(1:3,JSM,'k'); +ylabel('统计误差') +xlabel('情况') +axis([0.5 3.5 0 0.065]) +text(1,JSM(1)+0.001,'0.0575') +text(2,JSM(2)+0.001,'0.0510') +text(3,JSM(3)+0.001,'0.0029') +hold on +plot(1:3,JZM,'k--'); +text(1-.3,JZM(1)+0.001,'0.0003') +text(2,JZM(2)-0.001,'0.0263') +text(3,JZM(3)-0.001,'0.0029') diff --git a/JSNZN.asv b/JSNZN.asv new file mode 100644 index 0000000..ea4d42c --- /dev/null +++ b/JSNZN.asv @@ -0,0 +1,14 @@ +JSN=[0.1213 0.1092 0.0587]; +JZN=[0.0971 0.1314 0.1408]; +plot(1:3,JSM,'k'); +ylabel('值') +xlabel('情况') +axis([0.5 3.5 0 0.065]) +text(1,JSM(1)+0.001,'0.1213') +text(2,JSM(2)+0.001,'0.0510') +text(3,JSM(3)+0.001,'0.0029') +hold on +plot(1:3,JZM,'k--'); +text(1-.3,JZM(1)+0.001,'0.0003') +text(2,JZM(2)-0.001,'0.0263') +text(3,JZM(3)-0.001,'0.0029') \ No newline at end of file diff --git a/JSNZN.m b/JSNZN.m new file mode 100644 index 0000000..8d1451f --- /dev/null +++ b/JSNZN.m @@ -0,0 +1,18 @@ +JSN=[0.1213 0.1092 0.0587]*100; +JZN=[0.0971 0.1314 0.1408]*100; +for I=2:2:14 + line([0.5 3.5],[I I],'Color',[220 220 220]/255) +end +hold on +plot(1:3,JSN,'k'); +ylabel('最大相对偏差/%') +xlabel('情况') +axis([0.5 3.5 0 16]) +text(1,JSN(1)+0.3,'12.13') +text(2,JSN(2)+0.1,'10.92') +text(3,JSN(3)+0.6,'5.87') +hold on +plot(1:3,JZN,'k--'); +text(1,JZN(1)-0.5,'9.71') +text(2,JZN(2)-0.5,'13.14') +text(3,JZN(3)-0.5,'14.08') diff --git a/OPF.m b/OPF.m index 2daf42d..5a67bfc 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/newFIle20 - 计算用2 - (最终用这个做20点的计算).txt'); % pf('E:\算例\柳金Ⅰ926_21671693_2012-09-06\newFIle16.txt'); %pf('D:\Project\青秀降损项目\最小化潮流\最小潮流算例\原始\津头站津视922（3-1）_0.5_120%.txt'); %pf('D:\Project\最小化潮流\最小潮流算例\仙海919.txt'); @@ -56,6 +56,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)); +% save('20PD0.mat','PD0'); +% save('20QD0.mat','QD0'); load('20PD0.mat'); load('20QD0.mat'); % PG0(PGi)=PG0(PGi).*(1+normrnd(0,0.01,length(PGi),1)); diff --git a/OPF_Init.asv b/OPF_Init.asv index 453a5f4..724c726 100644 --- a/OPF_Init.asv +++ b/OPF_Init.asv @@ -21,9 +21,9 @@ wQG=0; %randPDind=randInt(1:10); randPDind=0; wPD=1/.05^2*ones(size(Loadi,1),1); -wPD(1:2:end)=0; +% wPD(1:2:end)=0; wQD=1/.05^2*zeros(size(Loadi,1),1); -wQD(1:2:end)=0; +% wQD(1:2:end)=0; %wD(randPDind)=0;%一些负荷不约束 %wD(7)=0; % wD(11)=0; diff --git a/OPF_Init.m b/OPF_Init.m index 4ccea57..6ffb2e4 100644 --- a/OPF_Init.m +++ b/OPF_Init.m @@ -20,10 +20,10 @@ wQG=0; %randInt=randperm(size(Loadi,1)); %randPDind=randInt(1:10); randPDind=0; -wPD=1/.05^2*ones(size(Loadi,1),1); -% wPD(1:2:end)=0; -wQD=1/.05^2*ones(size(Loadi,1),1); -% wQD(1:2:end)=0; +wPD=1/.05^2*zeros(size(Loadi,1),1); +wPD(1:2:end)=0; +wQD=1/.05^2*zeros(size(Loadi,1),1); +wQD(1:2:end)=0; %wD(randPDind)=0;%一些负荷不约束 %wD(7)=0; % wD(11)=0; diff --git a/case1V.mat b/case1V.mat index 886f9f76ed65fb034856a4aca1c068205c2949fa..d303cf982fdbcbf1aa8325a1780f1baedfb50ea4 100644 GIT binary patch delta 1202 zcmeyx*}^r!NFum2Pr)~_NWs`p!N}0c)Y!_{RKdu=(0F2?_QVACi8UPcH7pDa6?5KR z?hd{ZD#HFj{=9>L$%N%Lw_G1hRQed8{Jn)m=fq-%xTYiljz)!!Jwm*>8YN6Vr^*x_ zD@w98a~L%oHFj}OIN+h;%El!6e3Q?y!wSFh-u3)vkNe9EZ}&gAjjv-q;g2mqiOL39?d(S|C_1rJuGIKd8l*cy5C)u2YKB# z^4)s)vs_|kf?>9j57=UKllU~^M7>fN)nzRInZ$@%rO+Y#coeP1eEJ5th?s(4MV zyo=lIuT$RPg7f-4RjhHUTO;Z%j2~|9_@1%D`g63kscTSa?&^sX<>>`Rz`e|o)6@e1oMwGg3v!&4vgH+Q4nQ5`nO!#-CE^&B}&)9%>|YeVrl#%Ui-|Fo zt)F|Eo-W>6e9q`p|FoHB?9Q&xFyDISm5g$7-QOt*4T0I4>(_Kg3NG&YZL?eD&%W#J zQ$9~QCE3f{VYEr#DC+jRD@Xpby`Q<>t8dq-jB}bR^z{2DKArLTTfi67Uq4*<*Y{sM zzR%gH&Z58K^8L7VNnfHn7CZSeE)LzcWk>OZi<7)_3gXQcXRmtHp)n`>h3%BeWpgjS z54X>M_VD+upR-9WtNuz@$NRT0mfn_bx0`gSw6*fY;+mMvKSf{McRnHSa8Dxn z&WylQ8HR^09~U}*e7W`Q8ISi&n~?V2^~&wyw$(~nk{@+s*2)&n+uwJ?d*SRrN7tpz z8>R=A?*D({XLir4$I747H~jCe&D{LBbcVjX^482Jhb);7&q#8c{_Dv3WfQJneP6%L z-_h3R^}4nt_tj4+-rn_LPj$$OZJeJpbKkcgUo9EAX2tjAof4KOy6#%-OIvbNWL~xG z^P4N~M9cnV`7OpW-OlRNi9I%M*({%zS+0_)$mzJa{Cv>VcS(Be&$DB{Esj60qxY8m ziPq6yU$_1E+p|i{XB*GCNxvR%n07gJ^OuKKR|V_$eu+FaU%zYihj~k<%)B$jec9#` zz0c2oUt8T%B`3RH>~G6H^E$^DWTps|}KGxRFp7{iJnA%_iORqMzr9*Y8-eW`5wQ ZQ%6}8uW9d|bhnI?cTdH6d#+T&IRHLoPp|+0 delta 1191 zcmZqS`o%fHNWwQiPr)~_NWsWd!NAhW$kNKtK*7ks(0F2?_QVACi8UPcMJx;q6?5KR z?hd{ZD#HF@enXPZiUl1i68x!4pKj$kJ)4`Gd)?{TrKOeYKDjx_v2}D>3-Bm4b4@yG z@ni#Yljx_8j1x}kJjV@!l$Nz#y2L5)=u=QP)33aDK7VxMGn!y+pdS;ZgTbEf{fK~&F^h+lRqyQgP;5;v^;0T5_Vn$AT?#08gl@i?o38e8 z%9WD3oIR7=PG0(aL}T5ANt*t*|Ac9(-debDas9)IQ+!Ui+fQE}@M*)t%njO$j%g~r z^joIldGThraSE4D3)mnfMc zq-^zkeS6HaJJAmtW7hmp4(KfXcKb|Ltn1&BjKIY*U+;W*ee7`O6qoAru|WwnyGnFd z9^26#`P|gD{z?DZXDo#s%#*s=u69&Di@!E4A%5elpUcF!&ZgG5+jPsn=&t#Ez}(I0 zNY}6EPW=t*YZtwj(KeiaRI=yPk*Z(a=F@LSK(gKFlk%0 z?a8%vC#F`WOI*?Zn$9JyI^~Ji_S78<>t*fc+wAGDG=6rPHF?G5ow+4OtnY)5Oy6rQ z=fmk9R#TQPp;{h(um7dXT|KeouXw&)HQ4ralHRAuQAtNjgjB!gpS)G~ELcgtM^@{2 zdE0H*@4tQr7Dc|wo7-V`rhb;om37|_RZ8*iYFA9RY}tH1xi0xufka%{u(|=fd9;-|eiwYrSjk*_rP3w|3mwpmb)l_SApNdjEa7 ze`@2(1AD!Hbe`Ye{_$ulPp0uB50~OF&PQv!@3(TE(e|!=^RB68(TB3gSy^#+pFG{L Nr2gJ#@u)0@Hvl)`Qtbc$ diff --git a/case2V.mat b/case2V.mat index 0f7417947b36fb67d98932cb81e00f01e62bec10..e823a8e3e329423d44e11b32ecc2dd34d848dc49 100644 GIT binary patch delta 1202 zcmeyv*}^r!NFum2Pr)~_NWs`p!N}0c)Y!_@RKdu=(0F2?_QVACi8UPcH7pDa6?5KR z?lr$!F2nv|ezoqQ4kx)sMYd^j$y@V&2C1CVS-5MF;zW+w9y;#tt}PL`zE&br?wm<) zYEDkOcjXBIx1?)sx2@!U#aT-~v)J-I_WPWF?Q`!I$2~8$yZ`?CKfC$Af7h+kI&2Wj z#j}|$P%tj~xIy~Tz-gxq^@}at=GAj{^!iO$_g^XY&;QT!|9QOs8Gh5&o0x0#bCaHuu5~UX!zHVSNapc~TsQjWq??t0s4imEh3!g-zVZlKb zpMyMS4l?N^8b=(IPT)05;4Vwx=i9(4bCC01%i0aR_nO`wXqMd|Jg2dNd&eQY8@ze- z5BRDNFn?=OuQ+ITga6(KvmXbn3m9r1D*R})XE~r~!Q>;w1F>)Bdw3-UXbQrg}!kFj(uHU`fEI z(pQh4{y9D+rgQ!^o_ZJ0=TFU!@9mx9aqM@casB-OBRSP?IZ}&K3_sud#oJL`(Z(k( zuz&X2uYCvY&Tm@%`FEmM-M2ooRY$9)yDdC+zoIX8^Yfdh{JeIAuX`tPuqyXU@0YME zl3yQQ-|VK5c{$^sne@*;E!A0>!m6R?%Cej-jz4=DZggY+>y7hgUvlnTuww6|{#w_- zNH67Y>v?xX*Lz)vclK=Gu6*ro$s9W~w?798Z1gH$T>ZJwQ~x8&xvgu{`erKi8-2Td z)OnsqNKThrz@k~EbBcwY&lbNvXZp3-=JRvoyLTx0=*@Hezc+IBxwDF48}-$6Viu{K zxZ0{L6u9Zm2g}KYQzzN4%jQlt*5ov~a@y?V@*|5UOXnCmoLqRbej4NFAg^~%{>yz) zyr+M8hIv|y>&d-y?#!G!b;^?7YFUHbo}FJc3$>T`O^V>&x6(lTd6Hhiu?yEjW_8W4 zmA`kpH(PrRlh?24+{><}QzY&t&E)FR(A0nZ*pfeJ(#Jc~tUG?0D{pW87#M!9OI2ro z|L>@u!3%ot&wf<)yOr}-$Ai@R@YKiU?$dU-Ut7q(e?mbhg3$^4x(w?K{O$<)?Y8NRhti*2aSk|WR09p3L|m|qvDxp1d_f}P91H=bvne&xE~CG{ox+nmzx zc17QB?=p_6sGl>-Rd}`|XW)E3(fCypCTuCyT3o;Dxyq7)grc`^T3aecPQjVOyBS?}b0%(k&4)&BIH7svn0 zB-4D+-5u9cE?|l+Pg3G254GtAV*JPVqZs)kegBJBWCp@?kyG#D|ud(1m0*0da?VY&tG0&>(4WdpO>HC`~Caw?>g;? zGrZRu^e`?`+~6abw%Y96bE`9_>tot4rYCJO4!QMQ`BwY?|Nk!AzqP$Te}8)J1y` zM-xNK1NMIc42=RT2NjwGBG?xSussxDe61+7NI@z|LF|zO&&LgnLJ7(q4qTlZIL(}e zPAEtQDTrQDkbSbCf#*bv&WSeN1uZfsSVf;4{5|3Po_;)$+msRMvFFG5fY+7$9Z{*NCQY|}oegFHcd7Io!^sc`+EwH`hik`3Q z8>tKX4O_npZT7Oi+4-xxFZGSioAkZ0E4DcA`P7@?Vu#TKB4Nmw9ih@z~2UVtIgK{0EEHNs$xZbGs|5YRzuA zl~(gy)^hgLzMA>>PAt^AQ~z9C^-)tyzSk9-#et$lUOx)cYc(y;TxRK0HZR_*F!ynA zQFZpQTN!8MwWC5U6C6E_mDfbJ*nF3?H)g;8|MmI(DV-bQ6MUb2ZV>P~_uW!<=Cr_1 zPrSXqYMk@5y&eBY_nP#ty^p?G*1W%}(7B)U+4*O+U!rI0o?H6srhlo#^-#6?*IR!C zDrtXSb*ZoZg5KXBCZFUVzt=Bwt4n+Q?vR-NoYmfwrixyxy%9C{L#D*O>%BK$x9quf z^Xwm^SFJOzi^#kwF3hQSDpygRyRAHVYDa(6)!eTZ+npDMHOE?9_Y8Fkw|w;V<>To^ z1$DWT*Okpau32gI#QgMb$CIA#V}HKV^U|u{tA1tQ{KfI#?k472dz@@HPTo~^^u_X9 zxsxxsc%bR`A6d>V-`4WkS@D(0@n2>e3+%$a z?fGT&ab}obV@2UtRAnSJX>naZc3y Nzvi>J>R(pITmXcJLbCt> diff --git a/case3V.mat b/case3V.mat index 3c80983ea8ccbdbd1e2c3c931cde55a5f5ebe34a..38606e3b270adeab7a2afd4f8aa1016c571a239f 100644 GIT binary patch delta 1202 zcmaFF*}^r!NFum2Pr)~_NWs`p!N}0c)YQtrM8U|w(0F2?_QVACi8UPcH7pDa6?5KR z?lr$!F2nvIe)kS7x>fUga(=Uo3e=lu8Id+(qBUbla~_IZPz zjEM6+{IBkx&Ycgt7b$d)QOslQ z8P8v=xu4_v*RV!LuorERd-BlM!P>BE0n4ccOsf{Ky>GV4Xi&>I$abNP-+{d(QTfF| zwF~Ub4-TZ=U~1bS+LItH(abrcNpeOrYebW1#No0VjMfL)$_}(RU#I0lUqI_Qe7W zjRGtO1(<3$S`9eb&nq%1J(O^9G^t@U6yQGokmZP@&F5$P&#u|?u=bkQA)g7?bEALB z{d4jR`#=&0$`$N!o*$(b=snSYl?;e=q%tL{Gc zv*u6jI$;>raqajY-sY40X64n)moEOTckJlOz^U9fp5K`?tGe1)wQ{E4#QKC7i9*i{ z4OX6}ujU5rox*Z;sau*-c;ug1^*g3b68OJ(`M)DpneSF}>ztW?CZ|&Hiq}cee$9PS z&#J$agsMCb*!!klOXv8{h{Vg&ggZ>!JhM11dfZy_J?HwXD);x78CfOH_Xg)wH?8=3 zougXr`QBAAc6a|wn3Pf)mSb}5ear;j=~r&o)mggoUG(0Yx?bXWzuXJ)Jx($AEiGnG z^E)|r$ndQ1!!!Ud{m?r`XI+qEZJ<5kp~ zia7n8B~N-j3QpDgbI;`AGaKV8C+p{D6y56o)1E(R{qbpmiy|x^9dos-EKW*`WmT06 zXYUc7z<>X-bk2`?lhm&5xUTTY`qNaa+!f0|%y=5!v}f(#QdhlezkRLt-)OEXR^2vl zceF-k;K>)oKYAxA^?rJ=K+EgYA2B!mq`j#uzr&nNzZDt=oOR4I)|XnVvFOOELl;fN z>nksBT__@cCcY-^SE~7^y2IP<&35vVk`!FezLR@>)a}<|@8Z9&zNLOY@X!mNCqq z88PBp?9KkOzc(mauOjfhS}A4EGnsqGSG-XFsrBac>|18BcQ)QCpZ!Vn)FFd?{$09K z2Q6*c+ufdjSXXrFVC^|2uXAc*%YP=vZMiJXSF$^2VmWoQtp4{$pj;HY*1JSZY>U delta 1167 zcmZqSdc--wNWwQiPr)~_NWsWd!NAhW$lS`rLcz$u(0F2?_QVACi8UPcQ7jA$6?5KR z?hXzKm0|yIS~^it`9Q3w*@TEVP5sK#jvO)f-W>YGCOF~gQg53h4(?$`D=ZFNNfGwA z+2CU2#p%by{8+5D&68_Z3ZoXYT;tHho6(r;saNXI4jb2G*x4sowng>2>m_?*H}wF3bO{mJ8pLUa8T{wLwzV60a}M&oImjrJXl>DGyrV&TN27Jc0lp6`OmQ28_Z(!e zzsDH=p*i*ddp$>k14mN=M`Hj7Gb4whjFJP(L;==|0?Zo)*gsBS=-ePZDS=vu*tYK5V>|0))_u!v--nIWfSxTdv=hqjAsowcy zyxsI-{>AsJ-Z$RrxccUinET0?>n9>~_rhtVYrI$M#W?hx*4ru8dHt8AnfCtWk9J)$JLn)hX};3lGQOJ% zabI70Ej;pJvEp^3FRxjpw|H*2?^Uy4ukIy>#peo@FXfr-exg#ed&xcr)t*Vs_R>#% z-!DpeFfsjno5P#s(Ua=b3fBJ92`gB8r(NhWv+Cz$C$Bfrf2!)&x=smsyZDa%9?i`D z%h!FsELJhTSEZ^xY0{DltGvXXm6t9$P$5zKZs+H}&jhDTP&vKG^UnH>_Z~fZDd2U^ z>igGNM?QA-U*QYh-(W4>9%OU<&EIn+8#R5S{p*0AUo-#s-JEXx zD8EK?C8vDeZ`HO*au1g`nd#SiRUBHuC7HF{_twt4pZ@JnsmcBlo-6mFNe?-z?<_Q`iMXPh|xXJ1^wBUS&0cm6*6eW+64{NvEDiY?N` zA7)vPb>_q84;A1;k%FiMf@5xq^{_q4C5(?THEO6Kgo?=Q1-eRLpS= z;`W)GF7kKb2?yi5-h~FQjI6bf__|!{<=wj}ZO4w*jVm=$Uv#YC+!FEhB44fnhZF0r ziMnHd74&E|DV$m?x{{wwWWuWlhjaSc`M;Q$ee~ZqXHNC|?{kaa*QXcD zn4f%X5xzkBJJZ|yJnZF%0uT>#;;-U(!3qlS=0ze$STRy)NRRE9N$j zcd_I%>7Faus=7z(4Q^R)v~%Z~6FK$O%S%TyJ{0|QxKiqN?Y-Bg=tll{-7T-q8SZ)@ zXl}>;J?ioFbYI1P2L%r{iQet`!}RWETxOcc+oaDT5B|+Id~SYj$H|PTlb!ZS{aicO zwBAj)B(&6d@8hi>9_~B(taq2SaOTx`h5F(x>yH;q4SXd1Z|Cxf{12yfR-D+L|06_* z$?o~ya-k~+KBvtK;0)nt+9JMi{*x8oSDyOJ&ob-rAG!Y=KC;_n|8o9iv55DXm(v)3 zaA)<^73B}TFTbt&yu8R~;Q^1MC+h4K{=b}NBe)^sWBsz8(wjYZ_}+btE$8@qHs*eb zb^pa5kIb&WWA3@-<-2fI{DtEyvcl@r7qjo2*c98^lj8g7LDiY?8+(9;W?)L5wBX6A-JXF6uTTA|&(KKfj=@Z{SU$~!%lx5kBEikfxB z*z5m_E8jl`y6xMfUw`iJm1^TpIm|K57Of4|f=?^!1FD|Jq_+LtXW-@oJB&lJ0S-`iQ)JzkfO+SWIxKm5MI z?vwMC^KT7a{o(!g|6Xzo+r98BVVSx~U(#O*bX># zm)GmZ2l?-A$kYkF(Y?NY#g3AeKp)pbx7Ujr{JQQu_hr}5A6xW-ex?};sb~n@o-84{ z&t4_AqgFWUclh=hAM2CKsyQ}JzbAO>Mt#wbx%;Yi9N@LdDi=9#xBqqdcPVDyZvRb8sb2M_A!jCT8`{zC5HB;={*ca?) zL;u_MmVbZ$#%E`K)y{IpYhE+u3a@@Lo!z>k+qmZ#m*MiRV~+Zge>9B`C6%1WZu&o6 zhF$BqhV1We;R$QYOLFzT)yBKl9;k%Ez_f`O%#iJ_ILnSzmlq4C5(?THEO6Kgo?r!q4zRLt?6 z%pY($RpPk&VnJ0VDW@29>q|cW-o3I4h~ZaQxI!U0@8BjWR_O;W|CSXV@(P$PJB7*W z(1ruNmNw#9>Am%Shc5G|b8EinV_1DzuVRL%^`{N{#LBbxKIfjcIs%uqbJ&|v)w^BHc{XS12g z4E}1Iw3#{W$P~rQ*}JA4G>J=|DzajR&8c&#CeI8)f1UF`eely!bGg!tne+a>Eh>+) zx*QNUH*EItFRvb|AM5F;-orn$-u&3>!u)F<@9mu`-76Q} zGD~&)x2tb~G}{;H3opMPX@1_^J#8Jk_dcJOCXL5$$JSTBUuyI>Wv!*{v=bqv%j%s4 zYox?fe}2hHTNNGoK$t&H;?m3Vb?KAt1?AM2O9wA}yDoTf!lsB{r{5O(xW3daGJI7# zt=09`6d$e4Yoa;-Gp%108+=0A<3~ygU+KS-nrqqCT=(oM`@Ysk*KqQ#3dze0GpbJa zS$*xBY^ys%d)JZ7`dYSMelF>a?B6z~{ahNQ_{&;K<4jL&m0ZcyEa3;UeWVv35&dVOzVVQp>x|&gu*(1SOxL7s_xXL3 z{I+@$+tjKB_i|k%wrkCJdHK)cQ^_6b+wbb+*G+u+k^lb=ZKs$=wTIU8f6nVN{OaoP z{6+G7Gw+$}+TYE8(fYqzE9%ql2kt%|*WJWRGOh;iUccbqj1^a3ILxmOtP59r; r%{s7A;n^N$rzi1dccawW{~h?3{;aMb?|XTw)!F*je|RUFKeYq^_p0+L diff --git a/caseR.mat b/caseR.mat index d20c4f25216b3e492dbd2020d2c1c7240edfcc46..67c6031dfa40011a072807a62b6b2efaf5a363cf 100644 GIT binary patch delta 973 zcmeytzK~;rkwkE5o`P>;k%FiHVh=g@TcRq4C5(?THEO6Kgo?XEHM|RLpS= z;`Eu!F7j_{7`NBroQemXA?fX0&ViB}{T-X#r}k<#1v?5XW;Kaf{c+U_!M5937s5QR zRs{&`Y!lXez_@ggU)ItP%b?d%Gj3(vyWi*@GxgUEpC6O$@BiQX{onU}~+G9=ZGMqH{>OQVtB>PbNv_#^c z8HJMl_Rv1G-j zrmPcr))uPQJu591oml+Ln&Z;S=d5O%e&;<%Hw}>W-Iw!+@v^2!$)Z+>D8JDY)(~O-w&+yIL%I3J$H|N5{{cpOXBQuXo{W_a@6=(jt zt99`n8pc0O-Wc;AG2bWfTTN(LZ;ZnJEC2s%o_tpIs{Q^=nX>-E!uvi{PJ5!-(q;+giv`*W>OHD~Oext+f5g6ry^ajFPU_*FIW~X#`0|Px|``&*B$29SidAXVu?e={#MysZFi0xYOiKiNhvNW z<#f+=nxWrsRKH$D_V?QJCAU}Ty;iH&xvI>wS!``k;QrIGtL8waUB<*n5X6&q3~s=2flBplpwC*{%k z<7+Q+cnb!J*VpXrc*W}HS^46x`k^h+F;|6{55;cFeSK$jioMGHdO7Eh*~w+CEOE-U zdspr1J)-Cl{?hqp^5oQXGr|9twmu7t7d>?21B3hidz1G>K8koG@0a_=#&?g_uDO+8 ztq)CIW$Rz3BKs(yfBn`g^^HZ7miD}AJ>3-i{L-I`wr`f(&g7UbzU=q&3wsOuOMcAF zah-7fO3I9vdTZCF2kT9~cQiS>p836S|MiGf$1eX-3*RGr#%$w+>2mS&S0_IFmav`g z{au|iIU#RebR`^24qbJ{Zf@)H*qzt%V>q_&&-M_$STpx$(fN{F{w-Hm?Dea>{Z{xd zw=Z+~^QpQbUMKP<$HphB&t=@U?V;bT*PR!fmjBQ@)Y7t}{?Up5hSvMjWD~zi7hlU# qIrK>3yL4d5`*Z5vr@MPOKAU~fjdm-mKio5Mp)PZ!9pkC_0!9E)?%Rq0 delta 931 zcmZ3;@q>MWk%Vu4o`P>;k%Ez_f`O%#v4NGbp@NZtq4C5(?THEO6Kgo?3z!)gD(1LO z;tc3XlsInhrPXy%Li0$8&(Q>#)Gqey@aR5QvE41YyUjLP9gB?g&ONU3R70CbcDA5! z>K4moJVk!*X2qOc_GE?m^-ZbM>}TD~e>VI6_xaEN+nxVxcmMI3!kFc}$8R=%yHo$K z_ScseaTVMDF4k2K%k3#YP`=f-zq9`J62)k@3bPXumenl3A{R;AJ!}7cRp1GY;1$6Z z9ullyc;D%Ms^YeV)6ErLz2=p$8*e18cTlkH~i^Y^kUibt@5W2UfnKyeS3q7 z|Gs_9wUu9f{QXw?;mCsvXYH3%CKPo>2ipC3IO)}wsT^*<{drgWwbtkUT=2NGzmELXdvM(R`XwU=_r?cPXLH#7OQ<~@b4A;%zSOd9(fX77$~QT!-DR2?|4Kw| zaT{aW*Y!tx>bCb@s?<$2xs`C~>X|f?7oV1CGci6}`)O~4LR!kojFsgF|MWd`c_F)C zw@5p;TYIA3)Wq#x6a7yV-woc%V6rC1upw>f?{=*x-gjdYt}I$V>&4v3&I{hZnItv; zl;lfs*?yA+|E8a+@9NesJkPW=uxdj8KDpYHd-pb7Ja4{Aeb&sLB9&SrJqPbSzV1a2 z>ov0HJ!ISebMw#LsuM0I++QCkc_FHH=QR7SK+m!!4T-AOh>Pap-{p%H)e?nX-HkrV z`RKXftQ9BwYGqEIh`sdBr`=MNt)Ig{vEJv6G{clnMH{Q zl8^L0o_F%q#dyD^j^~_X?xlR+!uPAR{M%EHj5t59-2X3o)UJ9(-Ej>wyK0@<#~8|P z-}#JtjmfcF+P4k&M?GA6Xn|_atlyTr+aK6ZS3h$m z_z$}O{3SwvEzi3elUj0WZ?N(1ma2<;GWF+dlRxmEyXyMwYVjJ`TyvA`4fV{2*-eu5 zyR&4hqBh#^+tz=mylI|owR(Opcg~TTZ56V&bldc|2Y&p+_${|{WmYUl#?Q)SdzlJL yFK-u)TyTCPYxm*