From dd123299310a22758c49662bef9c5be10599033d Mon Sep 17 00:00:00 2001
From: Oliver <oliver-mx@gmx.de>
Date: Sun, 11 Aug 2024 03:11:02 +0200
Subject: [PATCH] Paretofront - case1 - epsilon constraint method
---
DATA_case1.mat | Bin 2768 -> 9930 bytes
DATA_case2.mat | Bin 826 -> 1672 bytes
Skript.m | 49 +++++++++++++++++++++++++++++++------------------
fun_1.m | 12 ++++++------
4 files changed, 37 insertions(+), 24 deletions(-)
diff --git a/DATA_case1.mat b/DATA_case1.mat
index 004ec10e0d4293447c5ef96ff1b845873931e076..f65feaed2a210fae64321e106a7291e55f5217d1 100644
GIT binary patch
delta 9551
zcmca0ddhc#i9~Q|o`Pd(x`Lshf`O5hk)f4=k%Ezdk;%kB<%tQb^{uK55Wwi^%g?~T
z(7?^WP%%gM8mCR@Z8`P==Vg;l1k5s?qIhSC&qe;FPu}dk?&&JX>1xtd@LI}ykxTQ!
zDNc$rz1mckOfU)9sp)lb=^de^Ci}MDuYUgd=KRh3=bNto8{a&qqr+q2nvM>YzthEj
z-naih`N#gFFD^UC2mY#mW)kk=(zdVt;-9sQ?-=`nc%HGSOl7STKCtyuLFj7{wyMnM
zb*E1n?NBx>*z4TxSUvZ?s;VfP`TUn+%cP7V_LhX@G;1*jda_T-I+u`fYS;dx<%S!+
z$5uVQlf<?ytoP&UtCDk0N=t4jN%-(7s<U#T9P{+aQop9o<M}jy%9A>?#r03C?)N>c
zku+Ii8Fe$#U`5W^n3a+z9vqB$J4Hg8Eqh-1Dzi4j3%7pU*m1bevHqoS|0@eNaU0&2
zeNJqgT|1NSOqR^JQ*^DnB<aIl-%BzxdG@@DKVSXg*n?jir7zlZNks6R_m;hL!ePTn
z6W2|X4nB#wnw`DfNTcRuwAw|7*TFo&^?Do&K0EkSi^>^tls@ZO`u$SfOkJ(0|C6@u
zz8O<lFT%XLRsT-WFG-(=O|i!>)?GO|(@sc!;a3yu^=B-+K7U=aOI+pPqR8mmcmAm~
zM{9}Qb>_OVc*<G!^M*GT+<vepU8?)z-0e#)6-%D@Gv~F^;?qebRS#KzC$TA~roDW7
zOsB`Y{_34xZ3&q@JJtU@N)nTN?)+xil}mqqJqmR>p`^IfV8zafjvps$Zu=Y?=_8?^
z6#sET<I_zdb2skO*>d%Jsa)+Tg>^~?Lv=b_=U+N`^7n3uI}#O7)^F%{eZ29&7w6iN
zIjj66t2%YcvhJAs|2}cx^KO^U2Rz;WY%8;iU+y^O@B8j{{jB4UPSsqTs=8>4?6jG6
z3+J?Y&AZkkCMo86+{!5`fAa@><19yZ+2<F7a~+N4Hd-uySnZP})P1u2KPTtAq>a0k
zt+V$@bw@6GJ-u{Kx?%XTTX+21kN&DMy8VfrJ8eyQ`6|V3xu-XTOO0fn?9~(bI(y@c
zYi~>L<}k6JH(%HpWu{ZlrS|XRp0MWTeYbzaTfdE*6LW1!X8MVTB}Z4R2x1po|0(X?
z>!>4+8~;xJmmw~e?OR)YNP<s0%hh~W<tC4r>q;AnpM7vO+VpU%eD}odbw-QqWZ1mJ
z??uZTIHa(2>2j+H%}+lWY+bZhr|PLr-TaP67t;3Jt$kM~wpsnn-@UpM51*=k&bzj?
zJ<?^n-oM0&=N2w}b~4_0QbbPnyoE|r&pQTXF7q&%dz@qc@%u})EYIwIZ&4X$*8M={
z=xlvy%b2~>m1`GGIJ&DrqSZfFM&;Cf+f9o;?^?fL|M&ZHd}ZvfY~Ekpu%o!pVsUYY
z^Pfw73eWhs|ISq07w^guadyj!`sj@MdiVOsx~WPxbm!VAevs_0ZWHo+^q6BuoVxJR
zyj%Z|OO*zPa_r0WIDYqi`ouL`cT~I!Y0|ygRDN)Nq=-9L{>^NGmYibiW&gX6_RLd~
z&5w1u?!3pbUxj`4&%>SUJe=<&e0!{{`<%brOxt#9%7NcoKizj%D+)MWv*I#$Q}CyS
zUV6qa>n-$;+fRshOHjH!%{aLI@Uo@vjdLxUO1{0nm;c3@x8tcg%R04!=mgn+7v>x;
zn0Wk4_iGNjxn*^e_s(#hKk-yg|2&rF{P*H%7v4|%vUQ*L^?rrm%)a0!Sy6W1GmfU!
zH`Q)T-Dj64qh+(*Vm|-Er{8@qZ~ppB;HR<o>sX6M?aa3Ni+p_N%XhwzIws&DUlqAP
z@X*7|8Sl2*x^t+x%lXU^IeGA^)0M@FoexfV{`LJ-xMlCBE1kmUFMP9}yDO!?IDj*%
z=<IWwgNjQzRvRw2ym30p`bwPA34ZVB)kTemfA8VUUOr7gBlp&d^^qR?S6@!3_fT#8
zmMD}Eb51eWKeynvZ&UrL`on7aE1zzteDH6Ar|`jmYpZIWtMSbWoIlrZn}7jp*W@z0
zM*-rg$2_(OKj9MLS1m2DDEV=k`@UY&yN#)n_4YqZFnvDf(cw6zub1D4o7O!j2)y@X
z-A@L`jXPS8++HtpZI7Oo|9__sk36bR`Aeza{P~zSwNAJ~_Rb;!WycBOruE<6=y@vk
ziFmvXmGmv}o-obA%a3ETYiP4?5sU8aJH~OE!M7_}?1EJ+YBjY@Q&$N{eQTN`EzDV_
zC3nj(HkjkG@b>RNFIaG$eK*;Arm{=b4&S*uqz=w{TWI^ZrzOHJWa12?ZiPjQ-xqD!
zDNx0p!o-^qBf4AMUv#tJqaQ+Rd%xDRv6MSA=lq?{9VjnQedb`tfs-FzfAd=?aKG&8
zl`S(b6spM03E1k^!7=$I^IS2Oe08y~7dyEFm;dfR?v~wjPrc`BTE60(ZG4_$w}fIQ
zt?j<MJC@^R&%)h@zp=Fb{&R=9jC<m=KVcu{9b1^ZY1a3pwoVy&CcYV>LJ_IQ%vsqv
z4eRI3-nx9!v0RTmkzQdP0%z@npLgr2UO2O(?6!JH#{tj2_a>RUKHbpHox7!`<leGN
z&rROFF4Hi!{>CM^L|}8J)Z6>A>y_o>FZ?+s?lk$@Y~dR%)yMZNv@?AEam_o{E%hn4
z587p&-o_o+R}%ZTW^Uci`BmFa<-gXkWok|hTz^0Kczv;+TGB0+PZRg4GGDP;UFs8l
zq4bqbK&i}<x_z0y=RCfgDZ{<?Rr`wCb=`k=8FxNN`gZJ>|CQ!bYZIMMv6y$=_kDQV
z_`pfi#iATX9)<4z@HSOqSJcW=mp=OJW83#nTKDq2m6vLl<W9|1ivJ%_zU8unVN>fh
zcH8~I?A1QMYpPqE>%aP)FS&QT);M+Bg&)}`*VP{ST{h$9oAY;Nwy<2fWWORcW}3;)
zZN(h5VcAD+f6)1D^4zX{)^`5BkAGj^SUxc}_PFCME2XE#r7_=nE~M>>OZ)Y=OK#b}
zk1@V`=e%1SruueW{EGK?QmvJZEbf?`ocPb{ui>_vj@C1sKUKdJ-V<%7QU5D;(d+Fs
zXDn_CzS3H}o<HTVWfk8|*3a+iWEbu;X;`;?&!$|%D?zsI&+je%ee}E0#@PESlEtI>
z^kUtwSAYH)&Uq#ARM)?D&h|rhU;TeuAF`+_((dkjmsT$Miq}Q=y{`ILxUB8{{LAxi
z`UM&0_*g#v$HH5KSFV?ezolT+>zo%^-}F86)ua<GPwzjm+!LWa_2_eT=EUT!?u9o*
zLv|_6t(vrF_l-k|VyR9aZ`b_#V!Z$BlV1@vA31YQE|N17o46`Jy|I4RFZ}>@&9n2P
zDyJXablq+Kk3af-%lE%pk@WD9z?YfZV>0*5HnxqADc}3^WcA&%`^)E9dB10!Y9+LP
z#o5(f_0E@G?v!6+c4vzJ-9HCKl9eMnWLB-zfAIEqgz>@ouCC|Uo+`bPooaRdedPv^
z^U;;7XIL3~n62rXx4nJum%r-+cFT(F$XWWlc&q0|Z_ce}3+KOJtD1K^m+#!;+MwH)
zUi|YBdpMu#&$=M@8mXO9TIZei*`<4F%GcYMt*`&n-MHg-eY4)Z_opA{nP05FR`;{D
z_URShGxF-!wi_S#7~Xq*h3V%t&B?DNYxbqo&X7O9BjE8m-|!cIpEmEjTbjk5zq0w_
zOP(@O(RnMMoUV?l)?fVH^lyLls~=w+YbEE;xXUB;S9;g#6-R@;xy14k7S!F{|DJ7Y
zg!BI}b;jo4+46QPpY4BMpZWV>`{$Q`%Um93@qWFndiPo8mwm58*Gk2n?swPTVCO%n
zaZAnRPqVLN{|j1G|7|_Pb@e|l&sA?}zbbLuov9}@e#zP6=RZU>hfY7u%2AOS80~&n
zbz4sV%}mB;zFs$Kx~sfu!ip+(FtEx`mihT5GUsE!ly(MJRZ01upZcG+UH)C4%ktst
z<$uSjOU|y{=4Z+{>Hid-j0M-qgS|G03LZ!;t2BMjx?R-kM*6o6e4d)mbF==xxpVpS
zw~95|JP+@t%{DeSIVN{8&P`fsW%Tm$IX9o>^oQwpm;d?GR)1bd{PL#u_Y)4xKEMC<
z-E!55jZasE|9PzaewMnJZTRQukKIydrV7{F@81?KGi^<LAOAc}?f0Kw-&d+RI#<wI
z!o=eGy{F%;Z+3_`{wn*yKQX^=EvR{JzyNBVKTb|ausCqyzyY=;mJ$jk92=Z|`Z&}u
zbx2A^NGS0<W_I|)*~7CqVFDw=4I$=kkiH-ZkUrm=oF<|1R~QPMmG&~8oU-InOZ};k
ziCN0ULMJ7(R2Hr9Xi@q+r6u6lsVOOPrwYy_r^y*rb56XlV8Z$B`||&uug>0EdO!O0
z_H(uA^K0MRR@)27%Fg}EXSJ*MwY%7^|JMKBcfb6<|Nr}qb$@;Rk2|Dm__=S-{JZ(i
zF865L3G)mYLo;T!P5ZLTVfX3vduAonKj{3LUVmc8Dy{GJi*?pcy*tyo_LULi-QFEX
zrqmx}p7!PZ`+sjFxMnX|(tOV<;rWu1+|Z*<4ToIz?cb;TlGn_~d`7bC;fHmX*zF3<
z->f>jTf%fLhsMd1@4u;~)qC*IT+p}0j=6E8ytuy2k#jEsbo=+~G6ojKg`Qc@axn3>
za`5+gcI?xp*H2w_Ebc~b&c@rZ$^H(r&qU|^$!kCKZ@EnA%B+tQPL>+nzQw=$(Ee*-
zbNXJdv)Ji$)ru!-%c}kt6>5hbJq~_TJwIW2wNC7>S;mZi<@o&18*sHR)ago}anV~r
zI<ELqRpZ>3OZ(PTm@`RDQktHA!tD?@UuN#VyXzh2sM<d+kty@8-+Qq<#g=o$-Yr^h
z<~oV&Dcdx?!pE=Cq4m|%cKwaN56zYOd8NI<YUu>7)O;QjQO&t)Ui-IB`tv3C?8+Nk
zm$w|<%qBEpR$9vQ*Sh)(EKH8q7pGax(f8WF=dEm6eb-fmhlY$>Ti!gJlhDb$bivQJ
zzyBB?_&3k{Q(MiOP3EhuPPYoz3;Y!}*Pd7Z@qz2Pz~rqbT55W?2j7w3aC_ykTL!YO
z4pp5sdj#u+9F~ah%zS@>>Fbo@`I1^2x?jZol@a%LU^m*|F`-`NP;}*4sdOi^9Q&Ku
zyO#4=t8w3td91=DB{BQS=U>xUc*-_%+tqJs)XG!U`6_6=^vRlLCry?mYcJ07G_U_E
z<TGi;^1W(B58mARedO}o8?`nTYd2UnH(b$`N|jWLKX{r`zE#FcZ%_TR-4D~`%bwkp
zI`cS^<;vXW>z!q)dD-8do3vx~XN6xvm%r{ve6!5&Zugo~>I=60*p%|FG3A43nb_77
zr`T_c8va{%{)X*`evPvg{tiV4Z;O<DThgE|yS(1}OE*(Ne&y5?rq06l7Yo<!C}O&#
z7!xU&=gqoaU)!JiFPG!YU%xK5vcCD|7IR-R)Ysv*-R2wKrpp-XPky`^)UGtc)c&a4
z#+F!)(<073%#4*W^&h{alv>qz-CaJnaHj-MPSgIa{|@(wd|CT6m76ir;QCR`>3bU$
z`|f|8C-q*jUN_dqZJOs>`|`f*6Ds>BH!Nw@&OIlZdSK;cXX&R)?Hk@r*Ppa;_QuU!
zRst6;c{^0{urB>&nXl+KbB@)<-7I(HWW96JkNbVN7Wr|fBa_y)i}HJ>Oym??YWlVO
z#8QLng3t92Hp^a@<<)#w$aG7{qUea$PuUy0iq!r8>9+cQ>r|@eXWekxVo#BC<)ntD
zZ?n(KUq02e#Q*Q7F5%P*{}c-{0ywPCPqhtQx0xj*VP=5-Z$s{PJ>{Q%Y@6O3*pRw;
z6-U|nSgoD+^aU<>y_l8~aGq_F%e?iCnoN$p!n0T(-H=U;4f^|icEcahw+~8U%`$2W
ztL!gzKKNC1==$4BH(vhRwWxl<HZO<peTMbnVP_h}<b$=7w{m4{cG%4Cmss8~-9qrN
zgNtv=ncuHOEe%3`JMUb%hRLNz&*O|nl6s29*4&9RJDh^!&d;-*%AK|@=sBNPXwySY
z`|zj6H)3D>y0ySpWs&UOH`BM+KXJOy(Y-!BN_EN&*1S`fkFA^kYx(wu<_-_(`qP*6
zJ)R#7yuXdVmnn2YX3^)jdCYIln6Y`jd8DY6mb`HNub-6xUxNO!2Um#N1l!bT-x0Ci
zs50-<^4pym^}7}Sx^~t`Xl_oMaw~-Q(0ie}e_Ide?a7}0_1RsQErEYeO**-|hC}9p
z{Qd`Vyz3XOl9-#j;r7e3=il$2-w`u!UiMy{`X~l9<7)Xx-LS)}CY3yYq5G!dtWY**
zdq>AB=C#i<Sqo~SHMY+U<y~WxqL*?1@H?})Q__~G2yXpfX1?uCG{?L7q1V^(UUT@J
zwJda_b=lLFJ4bAdT(;Dl|IvH(v{OXXouoCz`@~9SpSFMZ`PexXyTt$P!7fj7o>g28
z+`utQ?(Nh1No!ZPD$faxKE%9HS~q^y)lij1`S<-6Rc~}H_^n<m^Uy{1%*UJJTy2Tw
zLcanZPnQy0TFWb3`c%8kzA)1@H?s6V!j~o17G2M*ekrQ{d?X|I^o`R>{&}K@w=wVj
zQuF;v^Znn8_hlE}_*vo4w)%yTVE8kmyUORjFIak`%139J+?@J0|G1Lc-q^o(s%KS|
zm8ABrPkHhAiEYV=t~*&urUzW!F1Xllek1T}ua1$vaNyJVr!D_2;Xknc_xtvPi`hKp
z@jW#<cko@+jiQ}fRR!n9>qT*M{b9&hW4iWxIrEOXi>I};7VWf4`Ip6cTWF%f`sFbi
z-kd78*sE9EHfuV6=lNvyQ}r8-{W22&)^&HJ-QChOziPU}uD;ËHGXO5-Yu0Gf|
zzck^&USYx5M=kq2+JX{1cq=wKZs+V-9WU|GJXt*`;hW-t-VT%7r=KRtO*v#M`b@94
zO|9bIspF6DHr|!izWc$uzhlL`cVEulm>T49<Kcx@#hWA-39q{=`|fC-@BWJ}^~dIS
zaP3Lc7M^oD;q%3o#yM}&EEd0gw8paM?x)_LdGC(9Y+3MOyXL=e$%>Hoe{bx0uFN4E
zl^kq$IcERMY2wokcfL3+@Ol3d>!%Ck-ldoe{1x2AKUG9HZNq8v6Fn0Bii<9%|Fw>*
z*JnSwuC|}8Q@7ybzn-s)S>Nocuu-}t)Zvj?fBQr4y2TuIPw&=Ud3KR0a(eumc~OV*
zCO%phFw@6n5|@?pu0D>C-%~w*+DAWFkh$s4Pt!LR^Om%TrV0dl`z&Tt6>qGbu+(Qk
zsrHlWr+$bUuKfQzX|wp8;^K}I4_&hx4>C0+`#06ik<VcbjNGQh$W+!N^GjRO)um+1
z=cyODH*(c?aMkTwuyH?cmyqVSlQK7Ke@(EO8tUfK?la}`_S>NijmoZ>(<7yO3jal@
zeoQjASY5Jh$E*$>^>zPxtD_m6ABTzWv(IPVnC+UfTlU7+k1Ic1Gh6q(d(&mbgEP5S
zsKx#Xzi^k$=>FTMzrQ7zXZ?+ydCARX+qAN2ork76OqVbCs8_#Jb{F@TZ*xVC`dQiR
zt6SFIA!GY;Z+T>R(`*0d#_d;bC9X4ZQx0RYc9(5;KGW>7WZxG{d)=i%YFpRDu9){W
zVeX&gInHdm_Xb~*sy^1KVYJ!HJ7tkj$gdO@XRlk0y948o%Vj*<c>C$nQzC+@J_hI0
zI6@m`*}sc!jDP#!=hi&$`fefLIeYh~-k+ho=&jtkrR--`FL-KH9{lzppMY10p6HbW
zep7?e#J(2){h9dd=u6)u(H6EX^9(n=YxttOPVoEr<mF30doIeG*>P^)r>17*&<75M
z=CeIsHLPg97dn^iM(>19nfjNWC*p1_R197qwXGs&)lRz$jC-G(+luULk+oH>zq?9;
zn{z44(v_VvLzr*p-P@rbZT&E|vTXV9&{;byUY)ry!J|Wc;zG^em8Nf$?BfjVu5(Vh
zp3ZEYv*Gt~DL>EWDhm%<hh|Osdv$_^dB~j?*E#n@-cI<Q`^Mv@mixYs69qnp{rr<O
zPin)P7rUyLY!c4-ef-h;yNC7eSQZ|hI9Fv${jCk*{$cysgqHYuwCk5UoxZ;N)7f`7
z=I);y_$5J})3PW2XQ)sh^JVF8F~78~vd+I@o4&0&`FW1@;(J9CT5g%LC;69NZIIt`
zUT*(Q)-;`R_4TYt$M+abQ`T^Gy45T0uiO_U@pa~#k3KJ%HeThHFP@+9{oTjUFLxjJ
z^r%1dXZN1^GTj+&9h@P(C61ZXUoCu^ek1ps)Qr5z$_wX<ynM3i_*RD<=R2lV?tYk2
zlC&^OrtGkwzA^vG2?27#7vokg6?3Vcq~%h&tI^&y$>f1(*&`?YqfX~1U2vT$EdH9K
z^uWP#=Lr23(nfQa?F(C;oW8Q|&oh5np`fYP_=2s2*^Sfp{HT9_wWN90l}j46W@QWQ
zJIsB0Jv3f?np3^GUbdolg-+Yw*KNV7=O^)tCvATzwL0r^s!)+@`r56pUnlH;Huv3>
zWZQs$!R>dKh*%$Hme28zQC8%hu04}k`gLRTlkhM1_OWwlotaw~dSfE<GO6Oc2?2{=
zT0dD5!vDCcbW&`l?F89<^<RE?GnW;Ap7Y&B#^c3_jSrqn?|pU9&sk-C{Tg<kb$6rn
zau-MPACLQTskLQN?PZC#TCW}V#k{ULXvP(yC)am<VS>0&+&jG&j!uijp0De>V9MaQ
zJV5ge6YGNfH*P-^%JRD2MRx8}ney@E>hl~@37@kc-k5z*BttWD(qDtVyO(a&Z#=eF
zQDsVMsn@>qYZ-qnn;JYdTXRCiR0+44ht6ep{df>-!RdQm#9v~HNkj3XrzYQ}SU(-#
z68!wxhSzB!XQtQgo40)XJ%ew5a~n1-4|}=zc-2Gq-`5UH*C%f8c(Kk@b&|mQ#-r|S
zM!g?==Uw@CX`lFyxyffD!Zz$S-lBQQWNM2^y_nkk@7AVGBCmVqXGZDAX#9D&!Ka?>
z`BbswyG74_RrP3{J1$~!=p*lW;j7215@&xE-}+qVM&&iRgswB{3;pdFUozQCKfLiI
z=3xFlJBf!oHox#cq<7=Z)qTpUDqGgJJ<OPvB7Ehk<mx3azc^ma&D%G3&W)A(E%k1M
zdOGzS_hBxpznjN&%KGeb_vT*;n-^=%N|i3#JugGHEKfnvI-d1r&l)uYo*8p9R=;EB
z{@k&con_-|rRQ$%GnEy$^*y})s`&Ne+8JBRW8OAjx*)o<efNgyTT(uHFBKO)4_f|f
zvg@pa2@g8%*q^KmIIHyJQ$b>R_P(>{vnRAH@|&vjHPXCU?51M<r-oenn%dgu6(_{r
zoU;w+;a;}x{NvCi-cD7<C9Ns5R8+q^_Z=+0n>g#|yGyg~zy30{W!=|hN9O#x#B)57
zm)(!=+u><}_fN;JG098OwJd0C=~8+B+b-_@b{%h<2fKEB6>cfd{KRgMIKAjin8;Ut
z&Z%j;R_=(bW7*FBG=28@+NS!q+j&&>8z-b+*!M=&)8B1cQc@fDdON1(BkUg~@B2Tn
zN_z9TNiAMw%C3mVkta^P{5Eyh>&Nnr3#KeBS^Gy=MCgex-=!z<o~J7Inn_E}^_g!n
z`Tkud&h@8r{Qv!AnZ9mDxXq%I&l|MY8;kEe>9IxsVcn6bViqFO3r?%qF~sq^$k+S$
zU%J!Y@%2-s^@$aGZ!NDe+sDeXeGl*3jZzH#9-l&I8ooK@awPWf_lXmqD4U)(cGf+x
zbb65dECU9s^S+U9>@U9IuUoaIMP6-j{@ecbDIc~wy;?2(Q@x7)!QFFg%O>c(`4;>5
zMp!{l&p!2+MZ$X94#fw5oR#-UA!U~F_qN(2@9NKAF3+3n;P%h#V%FinwFf?h7ngce
z3M9R&f3YvTJN9CuOV>6n$;cdA@uJpCT$ZbTJ=%Vk{b1Aqg<n3FXDfcy@8{2)7})T0
z{i(`pW;=ElWt7;}3#?G*D>D{JR9}Aa<c}j0dLBJJXXeSVz46QQ`682VI#%qT-6t{I
zRrq;l|Ju5=`jjd88O0f1YkBL|#dLkXXUx)+d*3DU@`mUaId>mrJG=4w@4qCnDX-b>
zWcj|O$DBW;-kzSg+GW=L<3~6DxYpRyVYB4P>AKtx^ImKHzj0ad!=wTZzwTu_`D-`V
z{X5^=!zcSprq?&K`Q?^<d%1g^1>C0RTzC1rrzT;R^~*QjZjQg|UHkrS&7JVFFfMmi
zorJ|xuE?lAnUCE4-U(MZxw*f6llR$g%R&2@Hm}8ACQkU7Tw{EZ*;-P5?}Gh?DJ6af
zlIJr9Tgzy!(C*kM&Gx!$*X3vDObR{C>h^r=>akcA`D7mFR_9Ay-h0zN8Axc(<4xGn
zI`^bp+->FeCpsn^6u(rwOSax)Pujg;=EEO?d*4slEn2qnRu%iBu9F$w-%O0RS6YWe
znDjl(mukGnzH@J^$?w$*cjqj5d^JU+Lgo6WL%-^ZFTPB;({h1>#asPX&XwbRPnGxb
zEu7jTBJ`~Ku+Z&;)*m0(h3ym+*ys57;IhtG-+61*!~a*=aogD3@4Mv~XVzAKy#A=5
z$J@Y=_>d!dAEuYSkeb?~@{V2Fif3!{F{#yQ6AyEInD_PGV#V&*_wOxsHNH>Ta$r+o
zWc6BEpI<lfIc+2xi`zZlr)GRR&pjpo%qQ)XFH;krvS0qjcf<F;j=gFvix<DVcKEG+
z*@qu&TaBe}Tv+ww`eI&n>jRHI=B|5pOR(cY{m)bPR2??oU;9#fLt?6Pr017ZA2nHn
zz>;I1Z@%ED>59{O`TWSc@@ThGHEC(DiFXY2pWkAXGV{MZHQY_1w(Y8Nx9DzhzoXN|
zcH9W~U9|0H;|#tDyEpvWUhrC@;NP@`|9ykqRqZbMm)<_tcE(!SL+q&R1g(eved2ON
zEc9%i?fO|SY~9<nz(`rqtmRMXBgx;^3h5;Wp8uX1m>|t@V~<7e+-vI^*!(7b{g~9*
zU0+}B;F5eS_vGru3nGPfn)fZ^yE<*J6sJ(mi8F@-kM}MY_Uqt}>#Vlml4I4>V~xG>
zH|XU#OZCHT6?|1yyX4-osI$$RIq5~<1>b%1_eS*3y(d3$etl$7&I`HSnGdS(`d!E<
z{wjMyw)Ne=rud?zH>S$Xe7AGesW!2G!|1MiZ(9~d|6Be%+HIlX%{Bi$&B~a$Sn?Bf
z&8~!9aH+RlB>W=W?>w(x-=Pn0o~Xsn+ShvU(3*}3jM701F5jIp?K<b>L_w#UJjnr`
z|4TRZE{ngESsD5K(9^Fw>VM?U`*G(cr&-3<Ba`ZRd|0+DlW5-X`FqNjj!j9<Ut||;
z7kh5K@S?$y$J?(+oLISkTc!JDd8>ofNiU!1nr(TU^qJf6o8E<=Z&FL6i+TDaUQP<h
z{1tSLX`6Z3?hT5|4ovW>-l%t@?&p*B*RqoxGRo@2Qu^*bY}hctxc|n5f46V**FT^0
znn}c#(WT$;f~J^I*~O(V4}ZO$Cc?f^F)jP#o1<T?9P^$(np|s`EHJm|?u-RZb8pR`
zVENm_Ol0@MvM8G++*R}ILY|$LZI(KAX+^Nf-;-?pf71Inb38tl9k@6<|N3hd<K~R4
z%Mvo;FCQz)9PoR3%<u8D60?-$bE@?gSN`2oUsK|`?~0*?^^1F1=ezHA=^Ea#oX=}G
zOG9SXzgw(7eWN6^wUa*tOsIHpOYM!?sa0QX1f_44h278j|19r@WSY>5eWfxoU%2J`
zFMiqF%3Z&`CO64>|EZ>zYg%_7tUkB!!>xmb7py0K$p6S)Hd$#Wlj@YO7U#d8UbeEd
z`SIc1skgl9A3u1v@#N99?aTV#tzD<!mg+qHOj7msnD=ewAG@Dl_c^8*cFjv~h1Blp
zdW%0v-oF_z(SEDziS^PZiOYAaK5P-f_tf`XN^<%PMFGPNyX5vyxR|whO5Tih1y?nA
zWZEw}yqi@t>wh)h)Hf^m@1E&9(|k5(Qr7Oq&5MKe%rj-8R%|t{w`bbmt#S3_+b5~+
zIv@J;-xQml$ede|6Kligna#wQSH7|PV}L^X+jRGB>#i>@H~QMJZ|m>(>#h}c<V4J|
ztY)>YUbgn`qf2izG&^ldBUiBhn*N^uOl@5M{bu2;w)73ZwR7HQtvuByAs>?GZ*iBo
zkA3fp&(n@8R<*eQ?-VH$UskWt?!$9)qRo7T>l<!X{SEo%*VI3k>*Mcq3$-5`v$^gD
z)Ldzkz507a%%pEI>*t4f&UXDC^l<gn-k*o1H$?BAWOn!UvA%EHOr+v`-Y^FGGk^W^
zIAq7An2XY;Yd_7ueQ&m)@wXOPQ(a-TRa#XGidGpbKa~-UnHpc$m3>t2!>@guYwJ_W
zo){MUtjXWOp7Qdva_N~B6TZA~u|9QFPx0KICpDbbQG4zeOv%r!-BVk3>bdO7539aE
zwBGCSX5Q!A4`<$_xc_T2w>-}saNvW4_h;>REgrrsx$T9A-Zg#uR_n9*&5PV0Ym#1S
za?SYAZNz;npz*=xZgU6z-T%%#x@VKdmy+^Pp+4toO0CUyZ|UWiRvw7lT={sFO4*|@
znXbpXEd`!FIcw?{&K1$lVLL~ESyOs)N`jW;?iFi%1CzOAI(`T3?KxM>m1E9vJ1W4x
z_3Y=UD!Iw0{_ku$Zr!<F*J47_Tjg^XzKa;hbIlNa9~sY|^8ZztWLf#^I^(53ixu{t
zdSCAKTVq1~!p#S>j!TH@<>#5ERyE1iKl3wEZD;m8Hrpm9a#6#b5Qk}3ctlN>Fc)vt
zOI#lGQR&MO33jj2b8>YjZgbmgs|*daUty3VWnWb*U$#oOC_QP%wgtu|I{NR{U7z|d
ze6_Ibv=txLr1t%OwB=HLpYiU+8v@I@i%woKuU>zhds%z_kNV^5tdEKMeq5yTCgSkZ
zJ0@|tj(<<4CjOVHK6o^?cZ-F8YurKokVxJ)J}HY+$|qR{E?ja`KFCGiWk%TiS~Hi{
z)lNQ&Xa9En-+AT3{^GdLb2^THfBf#np12=hN;<B;l@y&JvPj$H|D1p4|MI)8T$~A-
ziZF&xVN6I!)JjNTz9hpVAjJ}~utuO3I)%aen0ZFS8}TzdixWh8ngkgbKFTq>Y`)D=
z!&q5-zxiuwtdPvftLHBm6faBpV#Pk;pD53c$J09*jq4`9Jv>cMQX=58s_xg%0tNa$
za~LD-m23Z2H#O9Kso^;ttg!dO7e4=IDaXzE{%=UQ*#AbKQNv#Ou$KAtRXZe>>~ChQ
e-znC(WIpF9cfIR29z6UpPo+2*7^LzHmIDBRf9NIv
delta 2344
zcmX@*dqH%9iG*8Grh;Q>x`KkGf}yFEfq|8wse+M#k;%kB<%tQb_1&Bd5Wwi^%g?~T
zu!4z!p<<5i*=mpALJ{^4*52ELPKRbmrQS@tq@w##FtYCQrL6Dw?(a<r40D>Kn#rtd
zT>apLlb6t6XJ5{(4^56<I5h7<_feCQnf*~SH_v&t#i&UrFy{63#|IwHIdA!=E?w@h
z!C60^$4;wzeun>V`CMOrSyb2~*lD+IwqzRL+pQ%t-mlJ{!oP7wXZaD?%Cm>|>M7Wi
zb{$QQ4fk-@@BVqQO<sPgOrG@dW6E}$3ykV^3Pqe=sN-a0uGHwvw>sv>s`-bnZ+YB#
za;D6l?P~d5k!KIr?OhiieEeX}-H(N?TK&~p>o5K)<9o|m*%!@gQ`GyaKKau8FLSqL
zPfVQDUQ>AV$6W@wtG%4dK1npKs>^=!Ay1C)YxBx?HBCz%Ov*mL)h0|de#`S&9M7b_
z%>DktK6#z`8OOye)0T2Y`s?<y_14v|>2FYbyKB~?IJvBsyH9a(XqdI{T~Vf3@G*Pp
zx)Wjv`jJ6DlJyUsopW9Ip(~ROTi3q&<x7M;N|!!6?az{7v%U8I1)mAwKYm{Rs>tZ*
z+kRx3T@%-1@y^vDF88h;-6qw^aCO1)+TaBouKvqX3l=ssZh!Y?W`*Fc#bxG?^dE$*
zP_$HvJ@VuCgD>q9nZ#QB+V@=SnDI<9evv>z?UdMhCGVE>sPNKbN{qWd@3l*P+7+QU
z!=%1Eg~^ERQdHZ@nFjMeg>WDIxvKW|@|BBs_NveR!V-Aw#HNXxmFli~mi>}f37Gfo
zwvos8&fQLKMWHOWSp5GNzB(&xwBAp;|A<8O?9N;I5kD1oPHH!>P0m?5@BhD>0p{nE
zQ_u7=8JFg7_xC@o(tfzGf3EJ8qW9bKcY9f?OtqQXUmvuae@V;twT~_&&0V`WaQDp<
z(y6-LOQtYCx?lZ^`>v^;j+IR4YxNcSM<<#-Q!t&jcaG@O1>eHvo{n36qwoI76^iws
zM61OBO0<uY6B05QlFSsIHC$3S!1dIaaU+|6vC#$t7l~u+4L=xpuCOzF*~pX$(qYF6
z(&2lhp4;F`wg|(A{pqK->aWiWn{{#3iN|eIY|bU^4xD^z?Kban;iW-$RK7g8^!A-f
zzn5zCOS4lNl_86wHTK1>b&oq;rL#$A(JT$Fl~1yEJ$^m={h6QV&g7pvGt=^OAV(9&
zRh5pJ-V-|B|Ga(i$NqoQrk`oyTIlsdN10>oO!I8Et0sB%^~WB*R+zo=WkbQ^42SZN
z<(pkNKkNxx@3HIF<ga^AAJ2U5EPvMKaOs5@wb|7sKV|OhkFa?YH-m51$6t+a<YN}?
z4_jNdh)=52*nWoaM7E-Z1xGL3$q9(e-@7S6yX@;;*;&!+UN4v{w)O7i?|BUGy`L`#
zkF={fQT<$CO8?!=t-WE*^{!Sg=Db*)Wb;_8rqVEUyUCMyi*J{AZLy!1C@$#MdwufD
z1r_T`9YekJ-fY#*s!xhwTF(0?$NkubYU{20uNHc+A8q?E&0zjL=HT=`4Th<XdKc%O
zy!3fi>F=#qfAFao301B6>RX`QV7$E~UoTp*?TyU2RVvnR;=-T!UYwco<(pXj{Ugzz
z7y`2X*B$-2e1Yt-aE|SX&-awPSs~xub7SS=_}ZKO2WHy+@%glg>4>YCRL}km+;-84
zX@}L7XNEpu%nme6**bAXpWju7m8U0fUYBg1^6IzcKMBdy)i0(7#LhpU`aSpU?8#XU
z&e_`=bPnlNhR&&8Frh~%wa1?=GnFB0^1S*xm#jRNi|<=m+jQ69x!=JL%tl>NXIg|s
z1m_s-^GTVwWWkmiVgFYL&#kz-YR)YuH$Odo(;Ca938&rXir+LnbYD-JciGkJ=XTX-
zYA*lGEUC&8bH%E9r_7#Exy`%(&R*bP%aFG_@wuu|ap%?L3w+)iF0DTuoKeIZyVX&l
zIsTjsXMN|#V|wl4p))i69lyoiw$SI%R*20pe>n9rkKCUV5A2gSoZgygYJS_vv9u||
z?bI|c1M$7Qldf++ux#D+noOy(;`8V4FO>Ckysght922}sVaCpf6AG4yavk3Hrl+Pu
z_HgvpH|C10s_a3Yi<V_+be-enzcBf5U&{aU%)B?=p86V8A0ZZap@wbl$;H<lB0i;+
zEvU8T5uBUJyZ1(K`J0>KxBVQ|ow_zl=52j2A+UDd-et^o9C2?sRjxg>-TP?SbPHd{
zOD9AP+c--v1igB`<-3Er#qoFZKD<5r_matx@`~;i9ISR?(IHozKP9!Us=rix!}NUe
z(p25DNxxQ2dir|hiKJikFP1zAX6DIH)A#(poO_l0%4#;TWOK=dXW}GWUCnwFx4kme
zymO{7YE{^phEFT66hv{Av9EVEQ0(tYSpRHw-E_Ab$LFfsU%kTHes^{9nQw{Lw-o&L
zp6aMorGK+>*3HWt%zcHYHg?FJs42g1lKhjuSfu#%qYRTP<!0K`rLVG>>(!gO$^X2=
zygieNz0<60may3V(>a+d(vPI29PnT(R7<_mBkc9S^~<*>r$5Zw7kxnM{P%^5_E#tR
zm2KoH(Y~Ey^ZZNGC$Hq$fm&s|V|u6dsb^|53I!HlUoNTQDb3UKxccFom2MLR*>@Mr
zyOFl1qm55ug@3lcS@YJpV)Jd6^YfY3uMyP<*;u^0pgj8p3-8)*|9(asTPh|X+_mx9
z<Jt#nICdX8b?E;EtNWKYmCLP?rry{V5WcNx(v-5$Py6^6RX@*+DLv>{+)!;}HMg&4
zua~va?Wn}S7oqj?b@n}m#V<ouv(xTs<WGo6U(|hUb<2hi1=4d;!e>3VdBc8t<yXzw
zJR9qidzM}iDSN3DzFMQYaiPqPpxc|KzyIwxtJLE2yr=6o=N+u^p5XIx@lnT`>+jf$
zE^I%w<Ls^r{F>KVJkL)!zUDXolxb^@fC^M?coF=A@koxt8OEiK4Wg&mFG@L_IdR~S
zK&FKJ0hW~c2@-Y;4DYR(V?kQ{85kHQo2ZIx{-~0`SlPzBT<ZO~O+_l<_RfdT%$#)b
zwT#2BaDxYZ@oH>mez<PqkKUD%;NUxb%a!xa2lg5ivZeg(`SSHWC)2NspAB{V8^1VR
zwzio$S@yogpTkYA_M7&zEcnyIv-(cV>PJZv{_(N>c^btssZLEperL>Q1w)&Hso4q)
K47-jP-Ua}z3vyNf
diff --git a/DATA_case2.mat b/DATA_case2.mat
index e66e8aa3eabf7b54bb1562ff923c4f07abfb300f..e34ea4de08f603b4aa3e5a1111416bbfb74c624e 100644
GIT binary patch
delta 1125
zcmdnR*1<c$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-ronojh^9O?@mAEC@9ck`ow|*!b8uc$wA&*mYcE
zNOEISYx7WMWAjtE&A6iDclOJl>zmH7tEet7zI}IJ!u~?fdnK>ew12C$-KPF_{Vr<<
zpPwJ?ulRjEdSvMw_P$M1J{Qcj&d5&`JXo*vZc**E(9f$G_dF{o_1aTpey`<&+>IHE
zrS%^gUfBJfygul4vGDmD+JDdR-kz4n6E1EY{r}j~lQHZ%H+e&+7t3zhzqMlJwMVnK
zzumL?&R>&e)V!KILN{dk_olxNPpnS9pS7_4?dhN|-`54~&-ov^CuP^XZq=B%d-U>e
z1f_msPQ119^tB?HjpuKvluUhkKJ~ZFVyDx=AHxbtkJbD2t=_B^!S7?jcVFjj&6T=a
zv%kf@?~B`-zH867-@csDJ9bU~Klj`9Vq>`(TlOs7o&V49x5RS2ypze--_Ea;{rUUt
zO6{0f`PEJr7pnHKA4%%ju`k9Z^Z&0;bI;uRzx|DG-~SMM`|G^x)}<X>XX-S~hy6}g
z{%x6L`|SEVvy-Pay{xzW{g!=`l=y6=LzUM0@7IQ>&*VGNHMOI5>fQwTnxyiHuiMza
zt)4Z#^plUvirE(m)AyA!eQxyGfBLrTw)ida*>g2>U;a2z*EhNEUZyMmj$JcW?cJ0p
zKZCzy$^U!l!WB2;&ZWtRpZ<CL+u7^u^JeZ8S$Iw*>G-05D|OhXiAcWTs!FXr`tSMo
z6>ZuJ{oi_jyYryyiNlJ?e-0`{?wI%WfaCse9n;TVe=~bdk)7e&`ZWEe8u|C@<R0Fy
z4)@Z(lkvc6!JhM>6C&TW{T2E;f#uwZ=iiK_zLo#jVR`tfVNX+_0au$214CFoZ^~pT
z#^;RolbM)~)qIrh=Kl9<<!tLotM50y{+0fC_d<5ZU*QH1`r_5t&iruQ#vi>aCBeaW
z`j#u_oe%6aC}d0d+w<k?drqcb7e5>7_&0uWx@>JTbF%b(jeqW!&NsQ*Z`#kY;7<?F
p>N_#3A0<ur$H(^PX%x$(IyDXXoiU#k3~dUgW-Bl-yxnFL2>=0CIi~;s
delta 249
zcmeC+-NiP+M8d5oQ^Bz`T|vQ8!O+Oc$jr*XP{GK+$Yf%m^27vI>%WW)5Wwi^%g?~T
zpw7U+P%-Cmaza7|Lz0=ovxZ9w2e_UZGj3!PFlLGrbQ7p!VE)R;uw9CA=j4NoHj|TB
zWQC3hWFBEqk(BT=;E`uwh*x1+J2{T!Iiu-hH`ZfCh3xBI&+mQTCby*RpD53c$J09*
zjq4`9Jv>cMQX=58s_xg%0tNa$a~LD-m23Z2H#O9Kso^;ttg!dO7e4=IPmY;){h#<{
zbHc^`H~Nek_R5E~%&)K7A+cnCGh_WuvBo9yIZwIkUAOVz;g5ML#lgTJ(QmX50HNVn
Aga7~l
diff --git a/Skript.m b/Skript.m
index b820493..15b4176 100644
--- a/Skript.m
+++ b/Skript.m
@@ -39,7 +39,7 @@ load("DATA_case1.mat")
% load("DATA_case1.mat")
% save DATA_case1.mat X1 Y1 X1_P Y1_P time1 time1_P
%
-% X1 - set of operating conditions for epsilon constraint method (2x40)
+% X1 - set of operating conditions for epsilon constraint method (2x200)
% Y1 - system output using X2 (5x40)
%
% X1_P - set of operating conditions from paretosearch (2x200)
@@ -52,23 +52,25 @@ startTime=datetime("now");
A= [-1 1]; b= -.1; Aeq=[]; beq=[]; lb = [30;30]; ub = [70;70];
option_mesh = 1e4; option_BVP = 1e-6; option_data = 1;
foptions = optimoptions('fmincon','Display','off','Algorithm','interior-point', 'StepTolerance', 1e-12, 'OptimalityTolerance',1e-4, 'MaxFunEvals',100);
-%foptions = optimoptions('fmincon','Display','iter-detailed','Algorithm','interior-point', 'StepTolerance', 1e-12, 'OptimalityTolerance',1e-3, 'MaxFunEvals',100);
rng default
%
-FWmin=linspace(0.0065, 0.3861, 40);
-X0=[linspace(60,69.99,40);linspace(55,50,40)];
+FWmin=linspace(0.0065, 0.3861, 200);
+X0=[linspace(60,69.99,200);linspace(55,50,200)];
%
-parfor i=1:40
+parfor i=1:200
epsilon=1;
minFW=-FWmin(i);
x0=X0(:,i);
- while epsilon > .9
+ while epsilon > .1
[x_neu, minSEC] = fmincon(@(x)fun_1(x,option_data,'SEC',option_mesh,option_BVP),x0 ,A,b,Aeq,beq,lb,ub,@(x)nonlcon(x,'FW' ,option_data,option_mesh,option_BVP,-minFW ),foptions);
[x0, minFW] = fmincon(@(x)fun_1(x,option_data,'FW' ,option_mesh,option_BVP),x_neu,A,b,Aeq,beq,lb,ub,@(x)nonlcon(x,'SEC',option_data,option_mesh,option_BVP,-minSEC),foptions);
- epsilon = max(norm(x0-x_neu)),
+ y_neu = fun_1(x_neu,option_data,'Pareto',option_mesh,option_BVP),
+ y_0 = fun_1(x0,option_data,'Pareto',option_mesh,option_BVP),
+ epsilon = min(norm(x0-x_neu), norm(y_neu-y_0));
end
- X1(:,i)=x0; %x_neu;
+ X1(:,i)=x0;
Y1(:,i)=fun_1([X1(:,i)],option_data,'sol',option_mesh,option_BVP);
+ i,
end
endTime=datetime("now");
time1 = endTime - startTime;
@@ -79,7 +81,7 @@ scatter(Y1(1,:),Y1(2,:),'MarkerEdgeColor',[0 0 0],'MarkerFaceColor','r'); hold o
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);
-% Optimisation - case2 (epsilon constraint method)
+%% Optimisation - case2 (epsilon constraint method)
%
% Create Paretofront for SWRO+ERD
%
@@ -99,24 +101,27 @@ load("DATA_case2.mat")
%
startTime=datetime("now");
A= [-1 1]; b= -.1; Aeq=[]; beq=[]; lb = [30;30]; ub = [70;70];
-option_mesh = 1e4; option_BVP = 1e-6; option_data = 1;
+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, 40);
-X0=[linspace(37,70,40);linspace(36,50,40)];
+FWmin=linspace(0.0015, 0.3769, 16);
+X0=[linspace(37,70,16);linspace(36,50,16)];
%
-parfor i=1:40
+parfor i=1:16
epsilon=1;
minFW=-FWmin(i);
x0=X0(:,i);
- while epsilon > .9
+ while epsilon > .1
[x_neu, minSEC] = fmincon(@(x)fun_1(x,option_data,'SEC',option_mesh,option_BVP),x0 ,A,b,Aeq,beq,lb,ub,@(x)nonlcon(x,'FW' ,option_data,option_mesh,option_BVP,-minFW ),foptions);
[x0, minFW] = fmincon(@(x)fun_1(x,option_data,'FW' ,option_mesh,option_BVP),x_neu,A,b,Aeq,beq,lb,ub,@(x)nonlcon(x,'SEC',option_data,option_mesh,option_BVP,-minSEC),foptions);
- epsilon = max(norm(x0-x_neu)),
+ y_neu = fun_1(x_neu,option_data,'Pareto',option_mesh,option_BVP),
+ y_0 = fun_1(x0,option_data,'Pareto',option_mesh,option_BVP),
+ epsilon = min(norm(x0-x_neu), norm(y_neu-y_0));
end
- X2(:,i)=x0; %x_neu;
- Y2(:,i)=fun_1([X1(:,i)],option_data,'sol',option_mesh,option_BVP);
+ X2(:,i)=x0;
+ Y2(:,i)=fun_1([X2(:,i)],option_data,'sol',option_mesh,option_BVP);
+ i,
end
endTime=datetime("now");
time2 = endTime - startTime;
@@ -126,8 +131,16 @@ scatter(Y2(1,:),Y2(2,:),'MarkerEdgeColor',[0 0 0],'MarkerFaceColor','b'); hold o
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);
+ A(i)=Y0(1);
+ B(i)=Y0(2);
+end
+scatter(A,B,'MarkerEdgeColor',[0 0 0],'MarkerFaceColor','g'); hold on
diff --git a/fun_1.m b/fun_1.m
index 8fdab10..7179491 100644
--- a/fun_1.m
+++ b/fun_1.m
@@ -157,20 +157,20 @@ SWRO_Recovery =(J_f(end)./J_d(1))*100; % SWRO freshwater recovery [%]
PRO_Recovery = (1- Q_f(end)./Q_f(1))*100; % PRO recovery rate [%]
% warning flags:
if FW < 0
- fprintf(2,' \nERROR: Waring: Freshwater production is negative! \n');
+ %fprintf(2,' \nERROR: Waring: Freshwater production is negative! \n');
end
if SWRO_Recovery > 100
- fprintf(2,' \nERROR: Waring: SWRO recovery is greater then 100 %% \n');
+ %fprintf(2,' \nERROR: Waring: SWRO recovery is greater then 100 %% \n');
end
if SWRO_Recovery < 0
- fprintf(2,' \nERROR: Waring: SWRO recovery is negative! \n');
+ %fprintf(2,' \nERROR: Waring: SWRO recovery is negative! \n');
end
if version(6)>1
if PRO_Recovery > 100
- fprintf(2,' \nERROR: Waring: PRO recovery is greater then 100 %% \n');
+ %fprintf(2,' \nERROR: Waring: PRO recovery is greater then 100 %% \n');
end
if PRO_Recovery < 0
- fprintf(2,' \nERROR: Waring: PRO recovery is negative! \n');
+ %fprintf(2,' \nERROR: Waring: PRO recovery is negative! \n');
end
end
%% Output of System
@@ -199,7 +199,7 @@ W_p4 = 1/LP_eff * (p_f(1)-pE)*(Q_f(1)*Z)./local_ro_f(1); W_p4 = W_p4*W_r;
if pERD < 0; pERD=pMax; end
W_p3 = 1/HP_eff * (P_d(1)-pERD)*(J_ERD*swro_Z)/rho_ERD; W_p3=W_p3*swro_W_r; W_p4=0; W_t=0;
if W_p3 < 0
- fprintf(2,' \nERROR: Waring: Pump 3 generates generates power! \n');
+ %fprintf(2,' \nERROR: Waring: Pump 3 generates generates power! \n');
sol.stats.maxerr =1000;
end
end
--
GitLab