From 44135f1696eaa8fa11ee1a437bd5208691b2cf15 Mon Sep 17 00:00:00 2001 From: Oliver <oliver-mx@gmx.de> Date: Tue, 13 Aug 2024 09:16:39 +0200 Subject: [PATCH] backup 13.08 --- DATA_case2.mat | Bin 1672 -> 9964 bytes Skript.m | 26 ++++++++++---------------- 2 files changed, 10 insertions(+), 16 deletions(-) diff --git a/DATA_case2.mat b/DATA_case2.mat index e34ea4de08f603b4aa3e5a1111416bbfb74c624e..455e211040c4dcac5fc0190a101cbe34ae13c03f 100644 GIT binary patch delta 9483 zcmeC+ed9a9M8Y>ePr<P?UBS>u!O+6W#L~*ZSi#7^$Yf%m^27wz`emvN5Wwi^%g?~T zu!@_3p<<5iHExs8oAT@himFu_!7j?GN+O<W7q2Vtdj6(Z-nFttNouN4>mKgM{TiZP zOItj+CZAkkIAw{AQq_X#*E_r>s)_u}fBk=d_T19*>2LR*+xvXZ{r9`S^B5X#WQsI2 z)cj=|`lkNA@cfhix7S<P%3VknvfX@HTzK-b=MPKl)-S%d{#38Si>-y(mhK<Elx}Rj z#>d2;diJ+;F2_CpW0$<F@>B{hYfSCCeP+{zJ>Ig%9`CFU=s)u3;Ut&X133}T&I-&} z8k@1TMD<7DwzrdZGd6aAobfN^m9(U})!8rs9p18A-_MsdEfUMqD&btkSHJz_vS4Wm ziFLJ$Jr~3)6x&R&I{D#SldTlX`uE!FHY~T>p#IuwdyB$<uYw&arQ7mjpNfk%$!zW} zoq3yMliJ%5`5=`Ym(ETX3V7YRe$AOF>J^sz3Jp%SDkyuU%?+6K;QIwj&uLd~wfcM2 zO;w3V3NE-c?bp?g;HURavV@<$YQsLap7UGf+v~m(=??|PHmEj6c~$Md*C{AA?bz<k zS66OkQhVLbzv)J3%#JN*Iu(ozwJmq_yY75mzpLd+_TQyS;d@gR&v5Vgl%{=9=j)3D z$8&Y!3XR-tS6rXDCVx_PYU19DTYoH8ZrpR*Q*T9Q>(1t{Lh-MXKOX+Ib*s*scBWNq zvK;5?J+}Buu1~p`x3={3^j%jC)7vX0yZ^K)SINKOyjA=?RO<i9MV||9AIiP5`_H>M zzgA8-`XgF<f0cC;Z+-L$yM2ci-P^Y3w|JKG)NjVsdK+gbmi{ceAjVV{xPSRSNio0Y z&+}_H&U}!T>HTV>i`YWFeTHpUniw6ui+Q?F%U|C6D2KCtk8k&z{jzC^+i%|7H}#t5 zzh41Q^V=hALiU=i@_HMZT)eH!zI*Gm_v?1#D&6?fUj9$eU+nJe%i(kGd2A@JD9k%- zCia{;Z|AiRj(;ty=5LhkKCCs_pPw_+R;&L`kI;sU8+GeqC+^yIJgIh(X!qZLo-fsy zm@Xe%tefu>5u$0l{zrXH(*c3jE1%hPtS(hqZ@aL|WZVCb-)E)j)O^|3A1TnZv8;dA zdhSS@tgjNrDV%F=slB~6>$p<k{d{B7jJx-Kr#qa^x2nB2Yib>r4U_Nrjz{uBI(_ke zYoxcU9r*a%)bGAT!OrjN?{8Fk@%FdG-0z)@zr2n*p8v{{Ig@w(+}#p8w%1$VIJZgV z$5+<voVq+cHm|}{X9<5%Pm8u&lU!h^-C$~>JR|b8jBwrbgHh-CB|pA?@jKYuBKiY+ zb8}OBShr7RZC%(dvq`5DiYFBLUx{Lx`PT7~k+=kp*1QjzvCbPj4*!^ZiIeS@VSdWd z`h{N$OfLNj6lq>BB=+spYR<6P>F+)2c@qNWen~#IDfG)vtCl@?cg=R}T2#K|mlu=! z+kJod)Y(26|4C8#%E4yy|4OZTf`QN^<yg<(Ml<FH6&~q1r~ga)Z1C*9{=1w0$If>= zbYcI(#asN(%ULgo+%}_e@7L~ybGP`<J{qr)b<)`5Oy!5K8>MRFV?}c|t2CR;<NqBy zb-|qam-3r_A39uF`jad2*7o~0W-RHYpAtSl-SD$P=-}sL)o-u$+N_DPD|oPtX>IuJ zIKBOI*sr<jzqk>@G<p6!5kAq60U6(~T}`w+%dGwL&YfzTGnMb7vK;3Be<4@=afO=P zrt@7tpI`hJ_4(LQ-EH0MhL4I;8};8FU%(O>{Ofnt>3uWm?RphwtlB*P<&*{Q-~Z#f z+}F+0x$w{Ie_#F^>mSnQ>E_(x(ByPgqC)bBi|q>K4gbDv?peji%=)wSp>C|n{*qNT z^P`M5Z2AAe`Ioh$*D3SGbIo|t?r^>oI$7xPj(20!j^=~$jb_LGNy+HMXFnJDY1FWO zZ^kZ7_T$n1%GOE}ZzSvUPu<8m{4np*x$64)DF>d%)la-}RPpNnZ`Uq)U%C3zd65BM zb-s+*V=l=vvycB?H_>IieEjm7g>xKh)+{f47Rp!ov25abDT#>E9W3*99DaC0W+lgL zGq&xIRR0~WK9}h5GUl0K)cSef>&heUAAh)cL4aTRoP!&+1#h`Ns#$ZgN?*0!zS(ta z@R5w8Hg;$9y3$)UH%n}bc^|X%u=|Fj^Y+<WGxyJooA)^F;PII=%KPIbHWy_7dbMnS zYyXz!H+#?U>TJt1F^t<Wz5m0!hO<Y%>4^C4JiSu>;o*<rFYmGOi_I#p78hIop>V%- z)^#s?vFB~sHCpRGCW<b4Z~oU@tXq}6%e4Oe2gfg_+hi>cDjnK+&2aV1{`Y@ZYvmY4 zzPYx~dd<D#?DZybXWo>H{rbzlC|Zp(;qqIly_@Ge3)~dE@4gM=q4(u~K5m-lI>ovF zTU{=r%D*1V2mfcao_hLcVQm`InkoG{fy*jmb?eU_N!!dIUAz4Mq<?c=vl0${coWNX z;!fE#8`Ju!2fiKo7QFf}i-1~Yn9JG4uP0`%U3goQG0nX7t#_dMi=MCZ=FinWu;cZ_ zX_i~wyDeI7VXq^sU{SjMVbUz&gRO=0+vjl<X#8Gs%4^naUpL-~dvlpYmM`t)%zRf> zzF#7{M!%ut&6U};S3W9xC#UmqPI!7RUdpxDN@I)m%e6BE{?*^zU8Zq4Mn6l>^oF+T zk<ORK%ccf>{vw=nCbe}1$Lo3fTh>Wyn&&UsKU2A-|NoZA(2}`tCj6YU@nyQAlG3xU zZ@w+FG<p=TzcZvmVEu%9Uk{q3`=9;$!85+&$#b!5TpwAqHs5{RevPYe*~?#(--OK4 zUM(g4bxrGyiyKsT>bEEKU9QhKlsHkqC`*VlI%EEcr)!=(x8M?-D!V5<z;i=PX?Jeh zZUMG~rF^BmEX%cmgzmo=R#Hs(mVf#-Czt8v@aZx_O8KX{uWhkwa_rtLFRIk~qae}Y z-93%i_NR4SM&1kbp0t|zKNe8QdjEW`nA7!DTQve35C7}V&vnx0T)1)bFNfxOt9EJC zSUdJlEfpnMzr)L!`J{!OTh7_1>B2L4UGa%I&Ucz0w?;m7nzeo1b5<w5y;{@u30`%I znBDSc&8s7c1xFvB`K0dj$9nk!_bUE&w)oBdx8u1BLjzx}j$?Q``MvX1{;HM_yB1aH zD9>tcc6m^6FI}-@%^RWAOnJG|YKxU+A@v;-OU(*TEck3a?d;PfE(SXl#Et6aHdjta z&VQvNxXtWI{q;Yyziln+34dB6xNDhP5Bv3^iYJpFzgVW`bo=<3%~rZHvrc;7Y~H0b zLwe2LtWVRL^Fvphxa!feeN+5ttNjHT7hl_+nAI&{d;M#e#M|3@9#>VT3k%K*m)oi_ z^JIH%TYcEy9cwviC#_a3*L|?~i9=E0w#gid>J`n?ul)E?xvoED^Ky>7-MjYfHGiah z^1FN5f_@I8n)y#kitS2R98H6pTu)f(eD{ez;9ghWHK}X?hkNL+6LaH#{Lt21T52pC zf6;i(_aaOAfAaRHLR&ivj$G3(F3^-;!}5Jo@2v$(Vtk$J8~M*)-5hm4aCL6@#CnaR z$G1qIma+>vZR;f%{%g8(;^B{b&R=n3UjMACyvV}r_SV0uUws_cK6=H*_chq}{Qg-x z_FWS%uG3_nYMsBHFSy$J+=XqYo?e)eCHZ&mFC+ckx>hp}UrKL#E9{rMrZlau;KC`* zOR`Jl>gL{;e|~hK+LHVA>4&1im^V!;llh~6B6jwN7Y#3i{iF1K4sY(`eE;q4gU|1c z`09$cGONA$6r6TPs*+FjX`Wc}75*>z#m2`t?9y|;%>1xm(eIN6yY$w4K4E43Q0sGb zqL%5FeP6|LEj{<HE4UbXm#N8ux%Gz4tIV@L8-j%n{`uW|^}71I|A9qs7CG0K#8>`z zi(xFCwx#ER+1uw|jF#=1@JpaN?UINv(^BRn_0KC!=XQObr!D;E&$ev~T7K<$d1k^= zll}A7zR}v}%G=wwui?4QG;4Vt?WOy|>fc>AF0o&+Ss`p+kJZOxf$_$>=6`;@?bX^P z|Ca5_d~U5ScHUlWm5;@Ny~XcxV|Uf-d<zcU<$gC_?BX?-@AhA_o|mZKD}U~}_?OVG ze^>qn7|h)fzyIXMb(RN9>P?o3m>bp`n0>!&C(9N)dr!Ronf+UL9Cz5-)Xw_*@1v|0 zZ|6)a*cdDQA%5DFj!lyfGTPJ!N`PAOzBf5dt~%ak_^>-dds5J(l-!L|CQZ7e>~eBa z%32r2$W7XTEnZGqN~tH~0zEbzsb6#Rd+q)00h1=Y`&xPa_uTS3#qVZTpZ_`c`MEQ5 zXFj(#exo9EGLOCU$d;f@6aGd2dw%rC|NH-gr>)&{Qibc}1BT8cyFNwM9o-$3u2<jn zexcyprt5X@Z@gxX-5KvH`q=3~=waFP!^~Qr{`YR$u;+W7<J5f$y&c>;<qhVSvCs1E zzCLTaePb!(eeTNC4aR&4U+wn`y`P%zv1j%y{x_yJ`?oK%U^#ZnN7Fs*=<mr-#2#*c zBqp|x<%;!hhUXLReST}KSfHn<p72&{`Rjitjq0Q3xXqdPIK(FCFT<02(;{50+8nd0 zZ*#WS-r(JMZqg@_%p->vO}KA=-?n|b%iASd;Wi)m)L7O(J?i)0a|`!-O{tA{wG|J0 zs<1rrDk*w&PxXLe|Bf}Q@7$Pq@_*r!(4~TjjloA#WfcSNnQ2BR&*j*(Y^A=ES+e@7 z36?jPIx1iLQSW)JHHNd*w3N5nuXs=E!vgWAHGXfZ&iL-L%)ZF+Rdcg#cwZai>}z{N z7PoGHaBYR+pE=CCe+E8epT9M;<$lSBWeav1Y&h|}a=Oc&<{sXe^Do^D{~BMTH+v>y zZ@>?2finS8d`!j59=T*0NgVV&Xc_iqUzS~(&Tb#&ZFfb^b7@r8^KSlbbGBlZox`Nv zH-1O3z4^A>O@3FZo3fqVu7wYt3v<fe{(Gi{Hz;w{ex;0@$e;D=GV6}%+?waMXZNd3 z)z6Qa2Y8j=TD5K#-|=}We^$*toadSDRdc4*ZHE@m2AysD8eb`MAG*0`e!~AqRqNf7 zZw|@td>(4O{z&!4oDE+eOM9HJ58aS7&y&5i^TX@5{TmhsUz^}^C-X?_3;oDdiH09) zHYPp#azuT?gKblu&OCHZ;Y^$PrEa%7-)B4%`_RR@En|t6;S|R6yCj4<1K%v^tJ^lM z%0*elJl#(qtWcq#`eBvl_w<VDjB8%S38B^NeI%aSN6wTNR-5>C*2fKVmFFE^WG7so zV(-D_Uy^wIx#Sw*Pf1V9`Iot^I-%AOe0lM~^;6>={@pMZeH!-W-Ha%^Iw?oxm?KV= zZ(W5F-&;)$PhA?RpuWNJpQCNr?_XP0o<GhwVy|w=ZR9WX#PmvOcu1Oazy)vj-(GL# zRQ^21uXdT^u^6``<FV}l2ELUuBW5x0-pQnEQ}4Cm_nkv~|K3!WD(uU=Sa$mS_JCE3 z6ALyeh`ZSRDz_B8A-ZvH)OO41ZXe2C^{bTaYu>Rg%zF9yc*jGL8y$3|cUQLWWBTtT z_3=`?=kvmY)dnka?H<gX#-8)&dBT77H~((2@JWQN`SMwOQbeGa5li+*jnhk?9Ja~t zsE9Qx3lyJRzc0gd&6(}d^XBm%3RY-Q^#AcJn#qLu*(_7x4VQ(Yv@@<dD(h^qZFiQh zYe=7cl=u00#SJ?aRjBB4?e4a`%oNFNq7?Zu%4qs$_XPQ}6(%d)7wAk2{-(`iz5J8! z%-PGsu4%tsk-)k&q@Y8evt)~2!{`4!YI~bE++HDAtz6sZ=22g55|qUi{P2E%(}adj zri!fBr~f|?EX#i5_HtT@)Wr+i1*U3GW?Q#(sqW_MT#7$g^tUg*zUt@0*@qS#_fdbN zs{GHR&Zn7&<IYXv7!H$W+kZxRH_q-bez)QA0^wKot4zFIFEu_tTH*cus?vt4t^3zr z7Thh%rPb&f{gNft;u~+cOTE>V`FGuG4Fg#IWTjntzGkEFhN=BFdM@fwywBJx%xjsZ z_kVY-mlkinv-Izj@5^qeiqHExY2pIm%~nlC-Hr1boaZI2GLGO;vyy#tQGDa$zl(k} z|2CLAW1eMf-6T!sinhJ)zs#KCkg$-W{!`_f3H`o5&k2`2^7t0f_Hr&yK?eW+dhL7R zOj%KZ;&Cey+2bZBe?NbLvv%5r4UIc_R%H2EU#Z<GvaCAvugH`Qv#mlUL>4-$Z@U-H zFTeQ*PZvM$V!N%=oia;zKAV?(W9ui@-Wk%DJ2wQnG5-E{zF<RA+TF<Je1$D-!5{XC zls$dr-xt4ak?`!E*>C4brW{(RcwF23*1pF2sjRW8@gln~7KZIQ>T6YKxIRXEnX>Sz z%6k{pdrv!TsZ(d)&rtUJ|54SyUls|^dZl?-@K*eR)2B{u+qIs@Zr}6Avu`!s4LF{# zD^g)@#q7N$ACrEH-<a_GVD<f_MLV|X{mu}zE?ae9OMiczr}~*q66aQ*NIRJF;d}G- zS*IMgect9!FRxsdvbBe4TZPrcduQhAe6{1gvAz2bpG`31;i!-2Qvdt?J3rGiQ*N#2 zVqxq4PubG%&$QjWIc-<pRnZT-=l=b&JN`!Q>c_#CRA0Dod@ivmcQW%S;=2=5y-!HR z&C9Z`H_80{o=GcLzhK#X_P6&eCd&njfBeY{Emr;U^xuv2JH2=7_l2;(bC)e~QoW;J ztf0j1HBWOE-^*~vW!~&XIS2aQ&+M7`MPlxY(no>MZ_nZ8oAdG3^QhGe%oWqjeoB`; z++w-)Ri)KMTk~&gpN8tS<S*7Tjal3Fxt>=^w`gPc50%eJ9f>myZtQ+KUGjDKhJcXV zoZOXcm$F}flg+%*cvbMWyLA2ZPT^?2JJ0rg7rj&dR(7jJJo9q?qp3%auFukcZd+cy zFtthQVtv$vDeMp1lg`fHa=l?i!@0nV({JcXzG!<U7oj$z>eZ$-v7a7<y!-4WJvYjr z#=>yDKtpWgo_TF8kBl}xsgC^kcD+r<<m3*6{I#w2^45Q3Y&SZ~*mp9$O;!4Rvi{7H z>7{lZHx0t##a1`Zy;XJ9vG2x2scf5KpXr?M0`$(kdRS9n_2sVC)lGacE9dlyJ-cB! z{nW9@-8#a`Q@!p_>gBB1;HJ3#L#(XFO;L%I*U9E@Ls!H<+wEgw%eicJnc#z0kEhl> zsFsbXdvfM?NNMQ*+<%4v|FULsd|GK3B<e0uKehVOI<GZyIg{__{;w6Wep##ayTVRS z*yy3vw`<o0Ln4hbvLdfGKaY9;cdpP4PR%7tCB$`w<M;6&{iUk(<NEr%;7|MH{=GGf z@VT|`RBUa*@uj6R15VG6dLZFaV0F3hUfmu!rz}>J-WwaeZJ!srWM<7=a5Us~1&2t( zwf8DU-5K@FcQmBGW-ikFd;X)gd6(EkPmlk9n$GVDWc(oe-uTPTdkl}Wd9^iyZ}9%! zTevk@S@D$R!O}Z|tOfJmPWyd%JMV_wb+fx-)ytxVc3Q5_QC%o=-q7u4AJY%9tAg*X zt}#?~-+m%LJz4zLF{una2dDjqrv3k?B7P`pX7IC9(ozpk_0}()^StuLSN}f!>t2&C zG_N~X8uHa?!<`+AANB+)Z+I&@;o16w?^5>fnbIu7`PPeXMa~!PmQ2aFo#leP6HXM_ z>TQ(XE%vr1R4r6!*7uq{BGF+SM_;I_?_DieWy8Db+`j{SKX;Y?dOb<tcjInuKOa|x ziSwr|J#(97SJU(*a(c<(<@L)nuRhb_)O0EMH7`hy$xx{K*Jf|u59^PtG0n04XP>5> z_2YJ~=L?7WzKdb=l}_yN=@H+X*Ys)Sg+%?jH%DS>@AK@OIN>ha!^f*uaD1E8Up4<l z29vhh)&8TCZrr^5PXFi6Ndc4gbVeO460`{S>{)yF8qeZ~$IAaa-T2vj+Mi%?pL(Z1 z>TFI;PrDX|UHmLJKhb!`me$WJvlGSFKL3<)R$6FR+OMEnch*Y8-7mCN{<BHa=GJ^k zGwI~xTW6O9s<pQG%+;TADn-0`*PgeR`nILr*s$@}72e&UuW#MFK8fw7ib?VB0KUUD zhAuKyw#RlYWj%6zN#6Qz2iEQ1|8`w6C+E(OQuW&YGCGY%yt|r3?uv2Tv|cvj4~upG zgqghMuRSwr^*%>TbLO0~<t<kxe_rFA&6QC~fj3@qPr75f+cP8L#^n#LS)3}2cTE;= zWc_w7`pojB3GOO3v&-ew6`O@#Rd>G>`tc%ZZ$<ZJ(H&wr78h52EpdLj>yvk}%9TR3 z>@c4_B6U0Jbu73uw>A~LdzbT^>&+$C8Q)G$JlVmq&hKCKUXh?)eUU3RnXG(#`<mk4 z-gsLxdv~_7w4~_uHHK^abXw<%^U0O2$!uluof(~fd_(ZeyEzvR`??(}vN^iG-TUDo z^XnR6+srtoO#ZngPtMvk^8bX4L?f&3%^&;cUNB#gWx31tk59eym%VcKb2n?n&w6A3 zAvS(ek4(J#olv8nY>lNigRKADdT?W(`!!DK@SQb@yw8$TN>ZNWS%zORzu>pcLOSXi zn?j1%fo)pv^VbCao3qBrP58|JAhlil8uw-_m*p<DmRYhjw&PBF?Dq|`?z}HfNzs^5 z@ln29a?2OidcTky_MHvM_57a4HyYnR8usfzM~{k(^ir+T(+BPv{8oM&XFbJtuEo!X zUshGrEt8X(mh#0-+0xc%2an(Mi)vQ_t#|x9IycPVK%U}>Q2iq_d#)ULm-{_aqj{s) z^xYr7lx?t=d3`3iP4>KR@aaf@qpU9;KZ5q{lW|(67m_9Wx<M~CueF}<h9P_X*}OeZ z{oT^f)SkLwa?0`jy5K)=-ttxKT)Fw_2irGy{G@DG^QkKHUJuLue*CFHmD1w~w=niM zFTR<v{15%T=d;rt0dKe4nVLsuK7D%NfoZ|J4>jyO#j$3SEB3znxik0qZFTt*8d>+3 zi9b6dsxYfU!82l)@C+`$6if9)aj6r6rZeh&-LgHtW<K{1J5craZp^8^xd$F(dP<)- z;1}YvB6>|r&y-1L=6<gXcewVn;m_sFxgXBVx_@k!fc4^+kMq~QI-xN|X8)Op07may z>$&^?zhLpa*PP0<^~PKgzMYcSJ>4cQch3KQQCq=r?bk^yF5+)^FL$hceMB#B$Mg(I z#l<gE=hxJitPwgiWs=VJh*HMMPh5WOG<;K)w7yp1uafd?YyRnR_tv$1ef~$XUahpz zAhhK3)yOwBhx5}vT{@AGtF|&+d$rI#qs^TzJtd8mH;??$6T0`_o@?@2!L)@6E3Xuu zTeQ>b!}YYCdo`}`uE`2D7rUM~J?MHokF@3s@6PPJ&Ja<S7xkab%V(AxuG?Ah+ur8J z-e+I0PJg<1fqSFS@BSO%j`!~E__}-1F2^_iKYsm)y|LHm;Ys^_?M+GN*4^FuELqfL z?vB!pJ6222_@VzrOYX?K=uO#k4mola-TL{;?Rbi)39G^Tiw|EmZu|Lsf6|hT+&`+@ z+@&=#N`77c_kGdT54T<}id$VD$FI_pcK^qYQ<cxpbSQ*!c?pTHn7VCE+x)v*H$Ax^ zcjKVi|Hp^5t_OvreUEEz3i01m=j*kCt@3J--?d|+;qTvWxOMl%y~64PGnzDB&UD@0 zw#xDKW%<3`GJ6foFBaEr&5ript(mv?lK!bzPp`9NAGo^2L^Xb|__sGgUoTj;%#E*~ z89C>uFW05K?Y<d%W0(~g{lpH|?G5-ndCtp9k=@HI8?PnlbFEsLROi&Yrt#^ck16jq z#TA_A{3*AjDYjAVb&!JsYiu{Kf6tqhjr&&iEB{|n8*m{1vi6$Rxn5=an{Rn@UCP>h zM`-U>u{=AUxxzv(nm?qw**rOYBlVTkrr$>cC)97sI-QrkH)_VN2^pXHemyvJvU+0u z`egOY|HhU-m02Yl?i`-=AnHVXLjIh6SJ=DmF4eL-ccZnZ?(fy<&W;bO=d~R-x#HAq zYP6-~J@Y+7HMf1$H)>B_3Vy1wbOHbRQ#q-TuO2K~kXVpkA%CTyWS`~R{jctP+}N}+ z$CkA)SCm8GMtvyzntQ$1uWbJP;5$Eae$k6npY1pQzIrTH#KiXJ=BnNzSr_Z$KDC)| zoBe;De|PIhLi?*;J^!Ad3%g`RUiUsLb$BZL$-?G7`<+`_p-m@msQ%5fpYXA}$?BN; zPP3m^8inq^fAGck?}Jh|jfwi)yVn%3J=eO`6|JK3|JWJ})q`AN^(*TXLSrj+Lbjax zsuueqec_@Rzc-y^W-U#)@cjGhNn!Rsro`L{3h$LYSiHU@O<9zcr$jUC%hiLA<pX8k z=l*8=Wl~dpdxLgib!GLu{;pYEr%D%1d!6v<!?MTBjDI&&+@0;6FJRqwIQn~l%gPn8 z^Zwa2te4xi`j&B*`gdoC*UfM17hS%wwDkHkBk$mgfA0FMS{-A1W6I3x8zEnu<>PM~ zCNWyCy|Bgm*r$*sbzje^{`LM^@TPz<|3-!VpOd~j6MQXOZyP@PJ@MrkFY6{fugo1+ zosX|v_)u;cW5N`J&h_bs-+lhQxZ7%1NXfff#pi9;zjE|eKK^+935J}dUe~_m9P`tw zKT>hM$56EJbzebo>UC%N{U>{NJo4Jp^>y=lmW{>%AzLLjFIv=cb}G}G+^-*=oHLu$ zd6GfGDppbH=*G<vWqRui7B1RX6I6A2epP{oSef}ZiH1APXPu7oC6~|o_lm7wT+}iz zN%;NEI}Ghh&2ws`EE(!D`CYUatn-|o9|@_~x>Ro(x!j;+17p*1Ik)Lg`5YdtyuI=E z`8BV^+xaf0FAOnxCcAH2%x=bgJv#fca^xE({pwK9w7OAycgD=80?w|tTe`cAWuz8Z zY-+h@ZO`!cw<hcT-VN4D51uZY=j6(&W#)c(zL|sgk_o?dvL(9x{rtE2`o`N;Pp0*I zh>5CBudv%az5W(M|NCzh-#vJp=Jl!R&sQ$vTK{#WcWa2rhh1Oigyu49>^qf|e_g62 zWBb<jXLmN<e!r!D+p1YzWlwj$oHIXe<^4&|E@$0#%y>K@q`0bKu7YrQc!jiR>1v~$ zGO5`u@77<eFJfX-e4KMQd8)#jZ(?}`_pLo#<v;Nzzwa<*nxmgozkHuO%k|?ui<Sty z32OC>y>mHaNrKJNkma{Hyv=S_Nj+_TVB_jPzvJk+mcQ<2?l6f~w(H;iyZ0W4aKHY} zsD8et^@`~WpIg8AlaZD4GlcEPyeS7wpWo*2W9sXx>Txdke=Katx6G>*FFB@fOHFrh zeyv{WE-bZBZG(C7{s#g9MXL2yHun?0ci&X~^)QT0C4aZ3)~+RjG3TpZD@rsiE3B)E zktmy^Xrz2eiRZ-LO_pZWI|Sz>w|Ojdb5=+W`_TSWvh3a8!-q{g*^W3=XMJs3$ayBj z@`=mc_5*&u<IV}^zgbst-*TIKALG)#`4_u|vY6_wUw9lJ&)0aT?TnDV`kN^x^`E#F z`x`5m|2r#Y#^q};Td<_#>uG^0oK<h{Cbi9-yz;`6+E_LfomHYsb1$*X(TiO8e0E#H zuL&W~8$%zP9G!P~<zp>oFYnLs7B6nHe7hF-+9tj2zygEKsxh3q8Sl@yIQ<*XiCJ}Z zv$R)nrtsG7=Sz}jd;h-kp=!{L&M@h~Y+H$Xf&Y@Xjdgzd9!Na3_jl_JgOKR@ZC<ga z=bgoSOZG5JUg%uq+P86?nuF5o*X!p?wQ;^T>z3;a{W^I<nW~b6fNo#^(v+2)PyRo7 zo@Xk})*G>me_hF*b*#5ICNT#~1n3-|l&0fv%#(ZmwC5gko~LV~lFIidAAd4q*Nz#- z8@6rtR;jO)&twZby8g9=a=YU5N51bG!`|Ex>sfmJ`nn$G%va^cf73ix@@HQ=DcM&1 zb>F4+VNY8pOKa<G`);^^!Hqw5>g=-y!p|P$UpC^gnpT%{QGSE;{kdheInx*|@5PnH z)Li75ve(aV&gR1x*1R;SP~+LX+CFYasFC5KKi2JuZ=a`qs1JH^|L8e8rgK*7N_WMr zem9eC|IecQWEc5$?R<YTa@zcQO=}&i-5QRnJvy>OwQQ9#pVYERJttaqKRBFTpSGbS zao44FCQQeI?kV-3TYE~LhpA<;k@BXuH!VfaTx7a7U(ZiWy0>A^{FmQ1gnsXs_&r_f ze@vQ7)YOG{wtwH+u&lSXe%^%F9CPY7$~->Ne)sUo=(XqadOR*GE%`jjjr)$%^raT_ z`W+giwij`-?$!;L_f2bGyJfr6-(Q@w#S(-Pe+kU7Vfr1j?c6i|$IsLGubw^hf1kl3 zr?xo0#Q$OwreD(U*Og?{e;IW&`TmB>A$+|1ub=BUa%Sz{7w@AEd_V1U{M<ajdXDml z&0EY0uht$7O+GImrEH|VFUThJ9@Fcc`KLSoNj9EfpTE28PW;yQf_yeJb<ebHQ8s0g zIv#aE>&e8;%C|MngwL#-ZGLQ>%9(9<!zOdi+_@}T+<1G__IW#=#V-up@UwjNMF+9n zxy<?ZwC|=mnU&|ytJu6*H08y)={Eg6OsWO-KTgT_UHw|T*=owUmL*$eh|AB<X`SF= ze}4X}#RWw#KOC5-xM%*qTSu-6*qv<owY{=_gIO-?eQweB$9vo#zK*$T^`j;B$vJmv z{<q0)Q!m-Q%euE#;f|pHnsN!&dHY@~%9L-AE)0C}rYPO*$phutsYy|d8&59Ue6v`j zCx8B{J6m4W&pBrI^27BFdTC4Ee~-AlWrOIcsKaqSKl`2XdL8dQWUaaP#MI_+WZ5*< zt#|iH8r{3iWt_{iUwvL*)6Tax^YdiNZ0C#Ye;1SB_%60gb?TR}4L)6#OGFH~Uf3R9 z(j-#$y0$Lgt|Fh;s&4PL$rIP}HVM@KZThFN;%!~<me9zL_3zmJGYI&diJTnA@|@9b zvK#C1st%FtNY$HF?nPg%3u6E0)R=!!@c7TgcIV^iBMp-KPkc-Cp3cLw;^m}izbZL) zh^w7zSTld(zJIok4)#@cZBI2P#AQ^qs}~nNmOlF5<5%Q{7w+GtGfbI3F>&hb?9dpV iOZN{p?2qYoxOAT7sq*w}Ih8i{bw#`^470BpOa}n0NZWn@ delta 1151 zcmaFk+rc}*L?XB}Pr<P?UBS>$!NADM)Y8htOu@*&$Yf%m^27wzdPi0U2w?Q|<!4}E zuw-OlsF?FOIU(W3fg=YFur0BaP%z=x;Plf$lZ{s)jb)QkkA#Fp;~TaG3;c_ZPqzQe zEF*Trx3oC?Q|X5}pM2jOVdiK5@4J3RmDOCCx_d!d{{J4nTk!ws`B^tP`}nHfu3mop z-1mZnKQo_Hoj!5$;iP)kKA*(4<<5Un+QiM?o#XqTr1t%P^vjnY9u&@eHBs~&pY-|9 zZci;LEqyM<Z~gSEMx6cKqt=6me%Hv(d!tnwRDAt(qSoia=`FYG`0O<0UhX=*M^^2F zNVHDntM>&ncYcbC-_h@Uq}Qr_-Vug&_nd>5XZiiMh^byKs_{qZfaAg^tp4ARn8mk0 zc{hDNe=56GPVuZf@wB(kFC^_sx&8i^Lwoe2jX9_G?&I@5WWCdz@n@F!w^FyAzbrUr z>h)OiKH2}x{Qu6+o$btrH6#KU7-ronojloy$EH4(2^NML3CRfzN^E>=9K1|x0_-}j zF(kRMskM11v$6Rp+-6+S@jLtF&-G1b*i}@Q7vH|SFJXV7=e?5GYudln+HO;SyMC9o zgU`>8_E-GA9zC*j4tw9GDW41GT4&@Z3LdOidbg<dTIlE1jC-CHlzQzcGQZdILGH#3 z#nSo@4KM6|PhKDNx>)%94eh^YcyCY3;|UkHj{bja>B$&&otwO&(~D)d?B80k^4g<W z+~4k5edn)9GiqMV9ibaC{d?12hbLAi-_Kgu{`Pdxm+$KW_UHT$-IKCwUbkw@+&y~v zH-b{XF(=+ydHPzB%*OM#R7$2kJ)ip9X0g-h;E!PirN`?1`c`k&is1J#;k&PMx8_RS zt=ZpV-}l9BP2aWW+izda=pDPJ|DXHqda<$Gj4gYX?#};b_*-JRUf#*%>u={*%KrTQ zcBOVqto&-HiwjkI*pDRj?ARA$lllMGr@3ct{onq^x9@+5z5R9Gb?ed&t}}I-=EHs` zEC04kvVC^_o!QCLnqJo1{(j58NlJXS(xFOg{r79b(`WLX=$hJ5J9Te@d`(jM#Mf=? z-&W6>Ui!($WyS0Zh3WfBnLaoA>_2_mbzA(F`0Tlwxi5d5sOy{DcQ4bGf5)yFtM+b6 zl%K(0vgH50bm5Adap%(H!%zP_{_X7b^?5UQiYz>*l5~90zm+=d(?lfSa8=dU9{u<H z`-(R0h5m27zukGz^~7Pt<Ua=$B6rODdcbl2w~py&ufLf+r^wFmZGD>lQjPrkb#f2y zSBHD)-^qAjwP4Tr&<T<6+WrcCoxpPL#Pe^)Qs2sd?65q1)v%{2(15E=hk+q1pEm`R zS+yBp3GD~tksOCJj7uFGL{Ck=t8$jne)4M7<JHpL-2Z;9oNYa6_5H@zztSJ?UdZnF zE8O5gU%VRInIEp(_@j5FBslm^-*V-=^MSnvg={H*d%k>q&&l-b;%7r0|Hdy)m#uAP zPL{r}@z4Fz`6gHUP5W6E{ORFYeJ5u1qofJ{_}KnDjbfQpr=}slGv>2`p-sWmYy}2} Ix7&;&0R$35K>z>% diff --git a/Skript.m b/Skript.m index 15b4176..282b6fc 100644 --- a/Skript.m +++ b/Skript.m @@ -75,7 +75,7 @@ end endTime=datetime("now"); time1 = endTime - startTime; save DATA_case1.mat X1 Y1 X1_P Y1_P time1 time1_P -%scatter +%%scatter f=figure(1); f.Position = [1200 500 800 500]; scatter(Y1(1,:),Y1(2,:),'MarkerEdgeColor',[0 0 0],'MarkerFaceColor','r'); hold on grid on; xlim([-5.501 -.5]); ylim([0 0.48]);view(2); @@ -105,10 +105,11 @@ option_mesh = 1e4; option_BVP = 1e-6; option_data = 2; foptions = optimoptions('fmincon','Display','off','Algorithm','interior-point', 'StepTolerance', 1e-12, 'OptimalityTolerance',1e-4, 'MaxFunEvals',100); rng default % -FWmin=linspace(0.0015, 0.3769, 16); -X0=[linspace(37,70,16);linspace(36,50,16)]; +FWmin=linspace(0.0015, 0.3769, 200); +X0=X1; +%X0=[linspace(37,70,16);linspace(36,50,16)]; % -parfor i=1:16 +parfor i=1:200 epsilon=1; minFW=-FWmin(i); x0=X0(:,i); @@ -132,11 +133,11 @@ grid on; xlim([-5.501 -.5]); ylim([0 0.48]);view(2); legend('Single SWRO unit','SWRO with ERD','Location', 'Northeast');ylabel('FW [m^3/h]','FontSize',16);xlabel('SEC_{net} [kWh/m^3]','FontSize',16); %% test - -X0=[linspace(37,70,16);linspace(36,50,16)]; -A=zeros(1,16);B=A; -parfor i=1:16 - Y0=fun_1([X0(:,i)],option_data,'sol',option_mesh,option_BVP); +X0=X1; +%X0=[linspace(37,70,16);linspace(36,50,16)]; +A=zeros(1,200);B=A; +parfor i=1:200 + Y0=fun_1([X0(:,i)],2,'sol',1e4,1e-4); A(i)=Y0(1); B(i)=Y0(2); end @@ -154,13 +155,6 @@ scatter(A,B,'MarkerEdgeColor',[0 0 0],'MarkerFaceColor','g'); hold on - - - - - - - -- GitLab