}0vBvq; AߑLAS{zoCmØw
+Nn`Vt䟣;bgO&l
~/ 5@>rѱУ AU
3\-tq
7 ^R{֘1q/Dg{`K*=jF;^3ǎMCmwR4
7 WDX%"G GPTc`EڹЙ2Ұ0 !Mr XG,@Wh /t@Lj+tlB* c4\k3xǺc_n71LO X,
K6؝wz\Ml$ y"e WBXG'=9)mkgw?5H7V xi
+HgXyf??0C 3Dڈ0TbQezpϤ29PXqwag?A>rK25NPxzu >c_ޓ,n,[#Y& "J">Uw#Yz|ꊷriؖOl͚ƑQdZПHnZl*̣-h'YV ĊKDҁr%=adK ":ՏmL4>=0Ll=MW\._vLk3SYL S6IJ>8nKR F |MGԒ/G*-_ C-@6@Vv2 O@@\:`G*Dl?`ݥ-K: 6@
BiL%p-znYgC ƸV~'!(fQEcз\=?Yc pccGWyo''1{(Zv3zGO3@;=*EѿMܹJ6G+j`L o~1W,Z4gz6P`7VQKڷt $J녅dFk-(Q1@N^K]l1QK8 @yI F)Mtg1%&CQ8(;1>F^ a~} UCҟo"d:t F hcϰ^ Q= ?$BVp
+ evcB+{W}Ln
+ZeEZ$X@z`'YZ@F c_
+E՛YT/
+K_R ^2i9 ֶ{3z uc6pAk䋟HTeNg0 J҃ɬDzS`Ӱ58XZ\k4O Xr;JvGk`
w SuEx\JgYBD=0]pP}A5dvpX۩'0zx]~Tt,l&sX4ٺ,,1OjUIEjFçUNl(nͥם+ %qr97DpR* U[Ѡ08L]NqT3*Ӥ|[ao\`>C NO&REm$v
|J67
+. '4j0[u "ՊoTZAl.
/r-RU VkpP 4ŷN2a VϵSeXqQ:aX
W[v\`l~r4=ҳN)2*?W\K8=N L*ȵ`4g1D|[6M8NMnخߓ *FF<$ǜ c\][*6YrlH1bq)(dZy%CC6 6BeƎ!1 pN/jcQiF_6K6kuie0֠}rAz$O%t,
++)Y Q|-I{0Mw7/{qS*3
X "~M톋'dBF0_AoE.|A ;hřI`"k6:1m,EkFB@FK nhJ]mNmAROUÑ}a5 LV'o4: ic쪀0BS=+V5(-]) EᎾxئ/S{
oP*xq#r&1 fߥ`A 守SE?u =îWg4 w7`-O; fmE#U@ `q
i}8=ۅyix O<(aꒁ`F& g۵ܬ; 잲4b1(= MXQ_[0OC=Y`'od'>"lFz`ty[ ȊtcJfHQf ^<7U7p c75C2+r
lb}" U0YDz{C` o>^ԀUY
+6ئ+Ϙmi*br+YS`~&@Q0N?l%BǝD[)43DN@!Dluӯ
+|e
+H<% x{ '}
M"Dz H)j)A"֟:v$uIh ;^4[qaeܺaz9M(HPa%_>af|;P8ΈG X vV:fDٙ?pq0lyW*45u^}πʀ,QvV%┭x
+ZwY폦@2aJ;sc.蹁aF{FK .tq QU՝xh/jB ǃڢ܊w9 섎i/&%bJ^tΫ^!x#tq
+0ʁ)[PtN[eY$cE E\EE5܈Q,*ce2Ӈ\v"8ZN<$f&%X3.S
+ !k;T4JD
+!pi.vu0ror#ڌcpbb^!CR+b#r1PbdTOM/i#`V
+;`+ARq-RM?Cdr;rqh7maɅ95#
+khʘ SJii׆\` ף"6'uϠ/`/Oyk[d XCmz^=*BڋyT?-<[D!ix9Bw,Cr7j85 dWh$MNo ^p5F:qYT4؟G <# GDE,uc0/
+~c#b4xt&]LѾDf°?VLؑlsa1"P~Sȝ✊dac:^0BbOsҥ)sY`f,Nw1(6\7qt1>L*Ceii."ꂱ. WȂX
m>! ϭ!7ۭY]Ai-I [xK ̫s"m " wFBuظ,٨h\ B "VpE`|d
+kZ\h OFS92zm ld eb?No"\
+5JJk9P w$X N.CLh b^Y1gמl8IrY<{T#t7 hz^ȦLNQ#+ʹ_}Ј_y8Yn|8F& {},[
+3Zi}L#΅u1-j
+ o 1° yg.BxL"L˵ W ؈6zW7tiQ2t 9M IJ hk?=wa@`> &yx=pt0+CNo J0` .L[E /lc;% )^%*(Y cpYy߇8h܆0uQ <; ;E&S0u(A.q;5tl9;oi{JZȰ.
+w\`7 ą nbû,0HP7hۥ:ESK15h;pvhuL_n#"J`,;`A iCW
X5>Eau} ;S3mйߢcc!:bAL ƴ#om;
F4*0?_ciT 䳮k̻g2k 3Pf3c+]>9!0<{~#1SFad>D3fB F r żSaŸrf\1p@谛roΧ~ ysqMh7.&?a[Oτ.`G
+\^;fw # N4;MEս̄{ڊ:F}K1ض 5XXD&^LlD+Nrs뛆f[E; é)tH]LCʢE]tvAekoȵv|wMzW [B
=k0oQ9*Dn>N&A| >eWE{4kfbQ
+z 4U)CNq(PnM,uRDce*Ri_၂_^̖6/ ^%d dUgm*.0?'b=^Dm/h,?e-p!.ʸ_Q@ l'μ lmȾhw척^m
+v v`rqTe
+i]ڍT&BuArqPЖ_:9A
,'
5ɶ?E3 t/o0Ͱƈљ)>`L1agn[NN8=I@=oMڌJ
+e radqǶva@_kͨ&&n9gCZXڵ#ov$ [8vZga.YR a *%S]8j`'ӈ!V0`4KC2 ?x
dFY Z
%> ! 7!FmUKSF+>ئ[@V#CCBaf6RU&W-V#Bɀ۪@d1\
Mb# {^c@&v6l)"sbzbcèXJ;{w@MVi. 4O1l9p̞F}s=f)HxavG\y܉V(]M}uH#*tM (3"9 375HJ`/F%k˧5l1s4uEW̗_%4Y%}0"D"|nY0j- fp.,Y6ʑ lsl'N
LzN`'&m-xyKrZ;Ҁ<":*^-^`d5y@"G(4 X?( 17%$ Yfe*ʏ,:\mM#ri h9uec'0]G(U6ak`YN .'MGy BI l#rb)e*8 0j@atnYo5kZ-Gרh Llc> )(?LD68-G|i`FD-pf: 6D+IDZZèACt Ac!#].:!leRE}D
LjaÞ4 R-;2&}`+24n͆8Ve Q. mm } 줈ieo0S
+wk\(՛PXZo!ަ=d0;ϏX͞[ :~P%
EYW`|$(!>ܗ-5)'u74y h H[=[rU]Njz{Vdt!s%1U}6ﰂw$D]vK
bC.Yjȉ, &gcO 1M?4a"ޠPjH(-#R ˦\E"F4=CI l#R5.6> CFJ;T7f
+KH죮L+#db%>{ S5WkkI"5 rE1] 60 ZJ"h2@\ρC(]'k]ҷVsɌA& %oaTSj7%7P JHh hﰓŗsiΡRNYwq5LFvd, (k@M&+m1D tٰĉ:GZ2Kt¹ rIRt}JhdѽIvkQA3syLںyHTjOuO/?3<+S
ϹƋU#q
+)dx~='9M\siŗ-,{(Xinx'R`;tyxP-r } VɈƤ[ '1b_P|$1V&ڻUO\ ^L"]ubj$Y=Husͷ2]& 6$RC#"/Lm z ID
+݊ݣh ]0'?
A@ :i})[%Ub-j¬ՏSs$疐GC{8v l)q@#o8bT֎CW4$y3
Pva^B(mK~ޖsB"ɲzTk{-g D&{8vUQH%%
(Ujx.G6t] fj`jD\Y9hf*`=$o"1W^
+'Я'1OZ m/
8A-[dBǹӔkB..X &aL*\zblNi8͇̈W4AoA# n:`-
:z¾ʌ`*ɝ#J%nUnXLuVum)#_ѝ19±7pw d'&]/Ww˒{r4CY|:CH? 6<2ZNJ.j%* o4T9a!@ KxG_=h;pDF^+LɐJtR'#aHb)S?'w10EHlHFDGxX h0D~'
+95Fla3O hG$͵+,|ei,@1VucfU;h>e+?K&ޫ""{C_I5- EweqϩF ^=]黨eU? _ r:MCP0++QFCi@
!;ѝ+qLWЖVph <9# # ӛQ1Dx$^gZ Sj"&nH\Ipg{u;y-hu,6~i]zJ@7Wn@ T" \ 5$_ūr3J/Q5whvL֍&uuv%mxV/͝xZH^?VZ6kf}7Z0 +
lwS?9Ek [&&ɭFTZ,Vaڂ͗;v<
Kĝ^W SLMIUT_92qUU ^9S $=c
+ ?EiM-*o/"Ln6瘇\N>4 И>+5LhкVXi||S_SΊy3e/P>eh ` \!¯fMM=帙 X@oб;6Xuo2># U{Zg)UCg28Eo"A% \W ;E:g>#%hy:;h "\Zꥒ8 | , $R
+9S`VeϹOT l*q'Az@8W6eG,F_$y0vpf|+"
:>7`,PJz9doei JՕNK jjERw Q )p7h=kϟ hȑQ7`_W%oN*XCE,~ W6vo[ hw(lCup<Esj=?/DB1OC+M cd )˼QjDZlBa[tes a Rȗw
JqzW飲ٚ&!h$!-t
+|b{zO"DI8|L9I MͯDi5~,!1Ԗ h /@ ̮z7)sTSccj0w^ ]vCH/,5 cDDr4$vwl; ~Z{tL NJKiMA2xX~*JyV,X'QgDjޡ e#`.01t:X~r=ll M{KDI "{
e Ϭy'GI5Y z+A7`U2 !\K(( >_~30:EO^` O N+id@q
L wwd.8Mv`&Mȅd}n:_\) vY E+ {+´ J\/eBe((9"b/
L 40ڟ+^Bz|u*e Ë&< G,8
+HdZؓM@ĥؐ/NSY%Rw 9L(Ou^' 4NWV
+ZxPD-'aj8VҢu'wj0kQkDPmFkQ :F;{)`o!'Sz݃мǰ;Kifg* ,L";LZJr<HԄA, \Q?ό9R̂&+ f'MP
+aˊ"
+Z[uL {(@sRtM{acW,ۋ߾Ӳ| tB0F"̀xI|!8Yz55i$!ǫATAI rA*El@㘽!J0C x!۴g6Y.RLSSaJ=*iNM*n
+!O8V @lKIi `>0)w=0tav+ EOG[wՙӳضǷSݠ')% b20( k3|
#˟jCA~K#0 af%G
+Koou ,Va9eA
+ lprUsDa*[ Z"F91&0 \j|WZtq4.ncfTR<ޓ6>F ķd⎇nM:O2,?wH6ҧ Z>5? ^ _ R^ɴPq^[I@Tә1ӊist &W ~ٓ l4T] J'_і5f:K_trP
+:'w+TnLULrmdt]);<%' l 5:wip;`W( 3ςS|9O)2vw
f 1 lpyA~ uʃ 7l=˯t,Rw0B[ux Pu{C"9{zRdB=+AjFS*5bHoTa(
+,o
+{ؐϹoOXg;\^~YGAxc0 $Y9! mbbq98pI`=ؕzQ&{@踬 w( 9Hfj&YpkۑvcGDxɕ
+3*rv<-<^HFNJce܌UɾfWkru*[W RMf)L dsAO1
+.4 p<~>l7s0{RE wZy&RcX1O2}3ܻRۖYw_ޚލGĢ{tJ cxXĹS
+lf&K8{ty,Rmŕ+PIIf~
+(;v~L:0HRAsL?1=UfPKYi8ᄬ)HN^nhzpuv#~P~lb>^:
"B+Fi4R?B. iU*!Br =twѵO ˇ<7%0 rUCdK?^o"c_@ŭ,i~5Eq0}J ;Nͮ87O_eO$ ܙ:֩mw層*jiH&3##`&|z&7KqƔGhգJ=m$n.R s4e 2:TƐ\^dGz^9&= y}]%}K)?J )6~ J|X R?q!$W0:C{d=iaW)TD@uh˖!;uWչ[)
#b]j!A?n\,`7Ò:GyKyd.h)Wc/o,9=Z4qbr1J&@#
A #-~ -H
%E!*S ut'!qy Db9@[ٹ/ eLɕ =fĄ\-]gI:Bk~v? %ٮ}~\hm~练fR$l7U0'xD^ԯbCK*< CwAv'Pf_I@uiW wc"^~??u|H%VV*C 棼;ڑ3Cb .ảiȪi`ҤW'"xh~x!yD-eŧh[2H K=0 c6v7/Pz9rӻt~[#m-}*f0MkkL#JԎ5jKҽ9J~i@d$ N8GQϵ }|y.7!HI
+jζ' =GX:Bd /zOp+~F6V`WBX=9ŲBJ)U̥K9CDN
+;)=`hu\+bj<`om؉;d݂mҕt55#d mRP):B~K4@Ul]}V%pIs i!F g^J{NL??)ٻJ?1[24Tz.JHTΟC0J^LgVzXʉvFroAbw8IenXE6Jm]Qiy ZC贮G#++O 1颪ܹg0%4)K8 ,[N+ۤ%YZ;dr +G0ߏF /jMG Χ 1 *SӔvEZ]v.NSY w}e{ pZF+UĪXݣ+ Xϛf%Ԕd Lp% a 2.y݊tZY O?4vϞQ"1 X|xuqng{,aBf6<.#1&oxh G%2 ȭ-O&Տ.y{,bH;za~!K F?VU ws:ؖd !DL\}v-
+ϹEu^&yYSA_n9Nb؏;"LZzܣQTn͓]M B}e+5Ց]sK>6> n?o>vK04M\GGP
C8QP!cm-EJ~\uJ*5-/?@-˭>)^#;{fSCKbK2Kg4eH˪|qE=jvlͶhO@pD=Xl҇>y0-ZGt{4-OLX36M|=+]dzF]~$7F?":)B@D Z걔1f%AX֘`?ź\LHcb18Ʊ Y; tF(=2G>;.Re c!F5P3ݫT Eܿspa)8_ i6~al2 gM$tޥ>T6Vx- o *k5^GX/eϲÓnc,[ K;j-@WyBe6h!BbO`(JKp+soay} dER۱~G:C-, ?D0, ^WED P& -a6o?aG]=,aM Quu*yGYh: y52 #栒ǺTR.mR*vcJ䋹 mE/GqPR喣.J1\}$PK&zu0='L(KP cjM [GǠ^ѧKHխ_ϵޜ$+RBq9O~06 ǷHFŗ毇,fnhn.v݇AkzC{cl|EGխUiMVtbzN2r"k &_~o%>oF&gOۗzӔ~8'> o*Oۯ~w!KxdфnA[+` d;$b(*fO]o
cٔLC{ :ovasme3f(WZRWOX};(۬<<0xaY&/L/ RǓ،#uAQ9|֠[<"E-IP<8gD?gIc 9܅XJq {L:I:t a%HrEpۭ,$&VP5bw R"p!3'őI3WOp.RhgZ;B8,bx7yJ}h[sv7
++fBlt1x@ަznz
?>]7-]hJi nm@@&=M K2'6n]oBr). ",bH8H+G+KN^NIL ݆ҭl '})(% y*]2X}ބW'
Nn [nxFi,VkN; :4<ޞIכ# Z["߇}+q[i#^֛P ebN5;)z]fYxaw#o|X=ę1=}_(_R}>m/*]F{Z|j5 -7XfYyn3:
+ˮC7a-6 6* ~bj&V^YԻ_M
+Mcڬ+icpa%o
+\Ca7 g젗3a(XD A `
+$_=!z-.oڒFnV6)d?O[/!Id az -a쎡::fL_$* v`(0VP t춫5r? `JdGXDjn), bnFDa0Rr䊿"|}в%l KaD$v*woER.\|m0eH&-$
+XA5|(@AXۡC9gl'a#
``:?xsHa lRn) UX^;1fa%F8}-rpf,y'"oo$4ѼԓcƸ/$;ƈe{ w"6/ yt}dA˲$4PDHw*Es<`-mR1A]T
8 qH&ۤZ|. bv@6CVbt*6MȆ#
r` MV_sC6A{6J`xbwkd,L
^~C(EB 羁)ƌ܂p 7Dhc6>1Ms~@q4[x{gc\^Jn^(]n1+2cx;Z a#Dp[1"`-NuA+noxMU#&N?ji\2s}EXG#`ehYrY4MY5`|oPb;-v8)T*oNB[ubM$V h"k2e˅&],lpVrԍ,d=F%tntsfbP-gn*4yW+QU02b!S(f
+)4BwvT|0|/Na:.͒Kmj`2q*Ju:Fb.uz#ޚfI_>_`h,W+)0:U 3#ljuѴHXJ]
+A#&:Ӿ)N^*-qf qih>q{hTDOjB+e\Fi0B @h!Bi Q堹(t,:8@Hܿ`u0J
+|4!ڈ`2-ԗ) 1G6-VNbׁ;!3@:| l@GLqD /+UP8I2ʐm)zQrӾmlɠ5 wc'%MHgY>Uv8VJ& 2ʚD.Ued@"^|e%2= h+JK&X@j't%8l tMU
+b
]xd CUǨد+G㗍+yĥ&PND$B"
+ kB iHˏ}0N|qwW $i7ellHy&Aڨq!#{*ÎJh h 0ts§ּ*o'%~Śsװn9 M a&ԨE2C1P }-
+-uqqZ\A Z<ݎj@d 3d&dhxrePmv@+sA'
+r5^|N3T# XJ 39٨ IN/ đVDK5|t|E q \h3 H AP2$l)5 Q.@Y.hE0kḮMj!]NV҂Gd ]Zh 9g+!@ƢŒЦYy˶jg[caOB H
I hV|+6
+թ~+vt/_H[@|c1ڈHDOۚ -|Kzc|:ݙ7AIbk:QJH^ax,b[͢ѕȝs`=һ \
+ZOڼÐ 'hh 4Q^=WqgN Hܡ}rDOd%CΥ;ZS)ԝPqa(fniBծ!\]3[-IK˥r^ ZNtoR{l5y^8VJVWq2_.,N #W݉vV f{|ä(Wqf F; q"#6KչUqPTU"*VUA댅,/&(>ijDQhSJMIHD,`
F`WXec.T%vSKx2Uн"-G&IBiAP0Oz${D3d( (3Rd4/8 B]h6eHf \L^#GF B8L;{5͚,A.tH,K,0mFsdžn P[jJ+uy-I̙Qp8b!E=zxJ
+T;߂yWra'tw
+H UR_JǪ6c/ !ob4r
+p!aD[䎈|/u8ŶKf.V=>#&D;d{ *KL ٽX>} 5ǎ 5- #*V[͢ ?d lLK#yOտȀ}?-4S`~ubN$/u5hVXm!v^V4'bY ۻIU
+}6c{8 DrCR!_%lݷRE&m5 []4GHpc35' ݤ* VՋoOIT`9(l Aw 4ĒmagN7a=i1~5
Q/WOJp %RTGR(S-qZ'4]#Cw z#2ժ LM< /UvF'P /y1I8l0Mhҫ[CD4[&+AWdOэCnnn,ych^>}"Hr`-&1`DL
*cy G˳)*WL%. +FOL\b1>E4 *M&j(Br[sV|PӿS薻Z(w%"kIl "}{ |D
+ `6/4`.OW q,$$S5/_T'^N5Ҁ3A[T)Bi.i x<`" 9؍iz
+
9G
c*1 Hl VE3JM.p eA`Ke˒8+ e%1D}5ZrQ*G N d)G yTp;!оЄ( Ts Y)rxoyrK!0!(ns,lT]nO,__iPy x. .SqǩEke@]E?tF ߠ] 5aO+ kۥԼ1n+c"Pf [BaD~u98wK1cN&u*`dL YP3H~Nb]QfE0l#h/ "D*u]2{bPt
+ TBϊ"/VV0v rdrqQ LxOw%P$1Ŗ>Zw< w7nCw/ƽ}HeV땏"rzU")ӿhQmf؊cN 6K 3v~nRJ_Th.3 b]{XRLLrǷlƨ`ɯ4۔W&z۲_:8OI- %!H ZID~y މ'+h_X"p5ɢ T+ZшHn;,&:\SEA(΄벩pW[EkłA>E] &6]LS|Vg +}<\#~yp y}t0ĽMFb!W? 05vQĩ:GW
d6m"[-J`aH/VpƋ`"':'D 8h_tkTp>OIJ|ƍI/Jzb.i#MNz>,X(c&
+y/)d{;\/ lɃWX}H.8#w .$c t~Qm9s}*
+#g c2oh㻅E!l~r9vg[$)MZNLn1>,}ɥJK^ vlP$Oke
+Se1f.)JA
U`$dD~JØ ל56,5I{!8 # W=jҋKKD FM &zb/xS72b!Pak 4es=^LZݤJnfƷu s:/ӱ ur;DT)0~0R ഘ$ѧؐ3-\NMCyxV?lUfrCWꤑ[ $"Lnض(MllL9}4I.4_CNQ&EjDԤc9?" '#HJLGYC%,KW K66%< nKRL!-KT""g a?EurdaFGrbElNj[}&>`&8 8e";FiIaE2;ř"K?}v%Hl|Iy 5}'6 ($qjbt(+.2m_d%? t gJ~DuIGb0阘_JY#e.hLZ:fR=b^UǴ4k$rHB4PZoզ-1O`u<$nOڟ' tr #NQN^ >]x3Մ9"{;U 3a]6x n1|IIJl~v & JȊsjkNٜ do˃ l0#m:ȹYcco 'xCJX^*4#F IΗݭcL۶TDRMGfL K/5 _x9 HtctUoSȲz^g6ѦuKy2s{,::̰ m9_Dʪ nNO̞0hl>eWn`δq C^OZvf\&XƄ=ii09EҴ1pwN%F1[ȏ5*,_4(E BC05MeW lVZf{.gCplk V1i"Lmqw!!p IHzڴj07&V;r-JW^>gOgV$i+Tƶ6`)IcXLq|>xxKam5M#p8c"}pӧu, y禫*tW{v,\fe
jKK"s
ʅ *)249#v/$X kŖ.kL5yvO?F+ֶ$ُE,V2@iM bTQ hwfi0:b@9;̫ 8GaF6zbT
LA.K)_I΅[%$M۪" i/zj3squ_hK\vB`0+P +<}Hxda,"#sRy5߱(NWk8S%жTuG͏"= BNӤg2\8g,@4ƕ;i k31XM<>mAHre 7q{#niP:u+
],gjtDD'Vw4KHDO-i{0ch覨"78Nn4<o U8Q
yF^Y
+>H l#k*aL%ZNZ/^WRǝ_ W.sJMEOeTݱ'HsSjT%+
È]?vVR W ,Ȟ; u|ITjcʉ|HP& kzba`:Tv䷠wq[y85DtM{d'B>_
HvS_uX]BggL 2E`|'ArK TD5Oĩ wX3}SpgPu&s2ZTIE+}:c :H> }Y9w&Jo'$. xUQct ]~s&D^[?Q%YyrڪUIDFkF2NƐ4JȺT'*[ˁq0E _"M%&a.?hM[GO & . j\IW? Ak=U8Ϝt+dVS&t?DHjq^`#w="
+ڣaRK?D|NL! [Z )9 (]\뎭(Il(e xĪwX=H"m7
+
ZÏkA mV wbߓv 5tؓ, jW(ʁpXL]|^b7 Ib R)⛤:]GQxsCYl91aH"'W瀤1l81WM2VFG4;B qV')!ٌpCŘ]GOLx!{ȉ2BXW f!5> Yc4 q c[O/! #p%|ZȐ*h!s=#=cU\hv+3I MOȊ6c0Zu^ $SzEbqT߇!dO"I>w[w57Xa
+enGq+IҌ ͅl
++KX($B!mظ jL]) T W$?&fNWc&5/I%$\ҳ$t\w18`6g
+| $`;b Ѻ=vbj
+X|:3Ɉa&BwbOku%fjԧT
+6[cO- 5X]KMFF86R7VIuO H$O#B~uT:1 (F5ٰм*ot5v
?2>%!jZ4[b ۞TfU F' e58xcu_bbo^?O\"fdPΤd$q< ʅ(+=s1X!n.u69K!8:LԮa H!8 E'}Ҙ{uyBl
+A?Ko`®ĔK)F|.ļ#!q0P'"VK*WW.lrh-b
mC1Vgn%}3Nz<X6&TMe'yBib> cqnA'6S(UM[BOq;)?5;Z"JW (pkRXr mnӁCf 6)iRDQYIY5Q{9NK\4蚝,d.f=dh(V5J\z@PCġߛ2AƼĞʚ ̯oH kziW=WRSg]S. {8cf9HXE\X G,#&[Yv'͘JQ;`:5e' 'B){I-5Eg8)<ߕf9+Yq+tNKvuB<9q+'iÒ ^z*)w2t&
$:e֮uHcijF}*Y'6tӗ ֨PM(.ӁRVlY*zjn3N(>.:5'ZbuIvD /]%ތE>:?[d*Ր}-h{1'[G* 0鑣Pn8=
+Os~x`#R!ԤVSە5s:h?)fC1CnH7YޫwlՂ:hnƮȾ ̲:sA(֔[%tr^[ŧ֛꙽-ٯC]$_/4֜F@ ;0k}4Vq٥nKBy&F ɅG 4r t R3Gl50t$*`æ^=3Mkl=_ 5aׇ+.h,&ئsL(͒?2yffIϣ ?9"ݝuL1&L-X%vӛ18
ԻQVfׇEC{eG{'g> Rd0WisTԦʗɼ D{:OX?mcJyX {\@0Xiչսs?]_#WU_# qy4}`yCG6#"'}[ZS~;/4Jak1)}ܟr?.LBHjZi@cRE&^
+aGey"jhETdH$#sq.*V#R[lp ; & X㲎L./6ME`>ڬr( _
蓼κ>4PKo0s):0K14ϑj#C>7qpVcHP!tMXE"]juL=8MUUe"h=xd<5 hiRv]jEeRx'/(k|)l#
+ ˕&zMB{y[{U(Պ1=-C=DRl
3P[H:#@8 e6 M`,'m /ex0eO[RF-o5>M^:&nRCeY (1"( cf7~V5[#@ ٌH_劮6i픓P#XewùXo 5:O5 P>H
apzB+~uH &؎oKN$jyxF*0j4K)ȫϔgQ K |Ttj06{e5 7)ݵxql._].i)/d`hOsD- bl(FMCîS4!J5mjGXE]2Vby&8*GZ`Vɣ<%EDLEXe"
Vie ^!잶,t
0i;vmJ,UO-m[Q-I'vR [EջOcdb ۧ j &U')%1zr}v?ri`i" i
+R$MW#,UnjqgNFm K+э:Θ3)a,")Z7]H d6#g٨EO<25Ob.5jX@#]b\JgKn&F&JtUd(8Cnn-R@lh H
Oq+j6Eh)~um`4]r0?E|"1w> GIXc :1,|V"Q\YNJKMty}Gb%u uqC?PbCJd*ss[Xϲ`GS|zC={{= ,툽Jnj"A,/̶#:_5GZuXQ2HYү.>SSw\5+`nuZbT{00 5
.猁 TzpO2/rJ{KϪ(7gJ|2ؤ,-W29pAw[@qܱ qL0cCH5Hڑ!q6`7s#w4:m`:1xg 3 B e?gh@yX]kG2KW'5iCO"yQ9fȓ!TH얮$ˆ(YLmZ dZ#MXFZ%}[T̊ ,G.O 'r '*R0ִFBȡUm8:xwb(fIҲ-uL.r͡Ӫ + 2F:Gڐ;m`T_ H¿W Ҙ xQ;;%~W~{z$N.efqQ #pD7wF <:=}lw1}j0H6UIYk+)O& >W~AU`P ]S]xu]rq(O9fo1 'ye S[G4SuO4j#x_4,XImx:4&>R%G" -4Dz%^
+aWB-!x 4`!8K;eʔG"%==`* FB>q&I!١_ j芭R>3H+KJCnbkO&<[Z:Ln- 17R- H{ϻs1'ى=̄ 嫯gxf#NT^G uV` utbO^}旳~KI)SzM/ڑb3PH%{"zO :*鵾?~yZ Y>"KD=y`1@JSUˑ~ۓn+^y@[џOg3l M `M.ozL K7O5p[uIV"zJqޠ@a\u:ԧڋRj
!T1HNju :"Ybqeߦ$jVY8cJ41XzH Kn&-NČ[LMwC=u=SH{(UG3fEXގed*/rQAiOm^p} LXR7cx\vOC [7vqđC 'ˈ[9t !}H^=%H2}d7Cc 98NUpG^r81 ܂^e QĀ$ԷU2;Ӟ nLzV=nEi_SPQ.V_ .|5ōe|VfyE鑔.W'V=Jҙ|QRU\̤]R`gbI*.IF%jr(c,坆:Kr[ITUOy :]wZ_T6u
ɳt֥ZOy-iC{EJ}"LmuWz*1xK~oC4x|߬pq>f;#6B!vY7#0g.vW9A8 k^(3f4 KQ߷x7~*o?iUĶe|SwUX҆ut0"A}AIelΓ`Q9d}F*_fD ]ņ}o`8a3IL0CrG5]!C!ؑ:
v/ ˽^~$)
!i2&
+y
l^~,ioCd ﱣaʾLZq2U!/^n稚ɑff)@\bS
RwCZ[o)g*;15z%owe J 9Ȳ\Zoт$vF3n0G+[/f?X*[ƫB:^/7nfS$&T ulD3Cx٣Ǎ:?A*x!.L8z?Kj/q_˽ %ȦL4( yd7 w껕īhjtFDYSX 6sMc"ˁ!B'g,L(9oW/#2B?3I]5p) WY2 @AMQy#e,P,.O0'&ª2
+'lZ+_}ۚ cp} ;QSPm!Tu֩e} Bx |+lDn@;&^a_av,^a>tw#x72+u`PB02!nwJTADMRjS$_exzQbuiM2$ʎxvCLJ~ˍ >:3q REig79!yѽ$boCʡd:j^50yo2
+p`TfbhnG4EX&
_e\] ]A{2@nW;=[=۸Xv/{nPo0O xE{0%Cx#umJ I9lˊa!M iiߌ Uo<GcZa
+4,)ll]`+mGaEe, {In胣ocǛQ Z풅WbHAFAaKȇ4x;1IJʓ~߈5gCJLKaSRD|C]3,礫M"A˒[6USa-O .>MW9l c+c&?gxTR r0{Dױ\Stc1~hH~?c\hG ;m.;s?y,_0Ъ4V#5zHy"R[M O=yqa.D1i3 3ߢj-iVp_Ո/1uQGQJODm cb45hrH,zZuϱ۟֝?~v'W܉'DOO褍H }=L
+endstream
+endobj
+185 0 obj
+110989
+endobj
+186 0 obj
+[184 0 R]
+endobj
+187 0 obj
+<<
+ /Resources 188 0 R
+ /Type /Page
+ /MediaBox [0 0 226 173]
+ /CropBox [0 0 226 173]
+ /BleedBox [0 0 226 173]
+ /TrimBox [0 0 226 173]
+ /Parent 189 0 R
+ /Contents 186 0 R
+>>
+endobj
+190 0 obj
+<<
+ /Type /FontDescriptor
+ /FontName /EAAAAA+mwa_cmr10
+ /FontBBox [-43 -250 1008 750]
+ /Flags 33
+ /CapHeight 683
+ /Ascent 750
+ /Descent -250
+ /ItalicAngle 0
+ /StemV 0
+ /MissingWidth 500
+ /FontFile2 191 0 R
+ /CIDSet 192 0 R
+>>
+endobj
+191 0 obj
+<<
+ /Length1 6388
+ /Length 193 0 R
+ /Filter /FlateDecode
+>>
+stream
+xY pי_-KaI>tXm uX!˲-.9|>[ 9JҐ 4I3lɲ)mnI#L٤nt6#IB,%H{] 67i -[N#oX5{Rk
+\=;V@胉خ& (o %+'f/"3s[q>M--{R_3__|% ur~naA PrQ7>-gFdPeXv pkW^ϱy2ok 4(<,Ll- X A D !@-X,;]aϱ$+cfdtAfEUyr?w|Sxpw
+\x
+Wx (H̡L#((KBUD:mF]*,P*䂌XxuXNvmzt-e/URgHGџ4y=nP5^OCF2`]&#fL2kK鰧,cS}۬ DqX&23mVsXoNofhy_ :bm7>7CI14z=㱁
+:3tGjrg|gdh{
+؎ Sv0XIF!W+DSILit3 k PXFYSX:'^A
b(5 bm\S1{b]Xo[&_u[^yihkm,Uu57 < k v#ͣNbZ Ty ˻ (`A)\B[R Ƃ}E8c7L
+2
+M! iܞ?jaU[%[If' }3μi}JeLNb95G<5(`KU ,!5g^ZuW L=u2@ o-ӗp}X --qJ@ b6+n.lA=XjrJ8A
+RNCw۟t‵͛L^Mi/)B4jm?1jsǞ.}eM}ϴ1gZ bāRДPQJ -]Y!sÁ
9 8tm'5 =kѩk=5߽R_x+}Mǚ5E]#wG z=ѨgGL3S&{{]x*? ."Qpzp6ɉMo+( As')#S&^ #*I}C=d-hDڎ&-/3ځ GG$|ph퉬lqe+_?VB0V uu*NI]*RF S؉P
+2lau,qHIֈwu/}:vINm@֖!4QBjs# f3L4A'鯢VA ./TZAVeSI
+PMZIA@krf,ڱxni"hV!~]!!֜ɐWg-mp[ oPH)xl,zc87o?;h?F^aʾ_L}Ӌ;
+KI91E橾h%]6rQdYM! ,fI}}>w]1P!"VM5IzDI!>t rޗʾ\N.~2gvI_Sҵ3ѵ{NST\HK$e+䃢bCH45~Dl{t3h2[({Y &3`+( P v%B>::{Y OB/W/ۡ)Zk+ 7LF)S`,F,4vZ`
+B9R*E,9>"X9@`sJ+S25G
+nWMm6ks]omV從?=D_ouD\MM3ǷeU qױcM
ؾZAcv~>/ǵI]'hI/^T@8Uq45&PU94#(Fq5kUS=hFQ͠;) KI+"0h͋PJ Ds'VW[GjZ=F8Me
^.?/^PG̗)]wnb{h皡|oNS1 *Cxz`1Z z;ۼ |p:Kп,;oRpDKMb|s)*=$g,)χ;Q)O1ڜm]-Ӟ=5#2.
.s/o(9~uyspԧ%D%aAl29-SU`ɔNCX%H$yĜ3MgS|k5]ڡء
+>:=X{ODl|4a~a B$|
HEP*LLc,N#sĀn4=݆3aDÑ j z~sW}ފiG-~/desi{Ѡpߗvw*oz㾖O:6z N ϲPȊ^PW\l Ei9Kn{}[sM>=,U } ML*uIZe"ůȉy1@;PpP%C#W 5~8= "nNjK5Ik}m,Vupk 6Zcu}/'M.?9 HF)VsbUR
+z#.%`xj)E!XfgPEP}$lpww*d|iLY4Hp2=C/# ziΊ&]3 WlIJ 6 mBCȄ/8.* Nmb?^8/Eh3p| uMYG: e)Kw9@y̳'O>Oi
+u8lGF@4K h%OSPHZ4[4<͂ :`!
`v^ Pu=fe1HAz-4 ǖ_`'` Hm"HYmAƱ =w G~i͵ngR+(LJ=뤝ƱNj} s}'[:-A%ZY֧f |KR `I/;0>櫫Rk֧,]s{cnff|t eo 7En1ף01: Id.' }Ljߚ{1aB =@;m u} HymSTTkR^vx]/s c(
+endstream
+endobj
+193 0 obj
+4720
+endobj
+192 0 obj
+<< /Length 194 0 R /Filter /FlateDecode >>
+stream
+xk`edX* 0[
+endstream
+endobj
+194 0 obj
+20
+endobj
+195 0 obj
+<<
+ /Type /Font
+ /Subtype /Type0
+ /BaseFont /EAAAAA+mwa_cmr10
+ /Encoding /Identity-H
+ /ToUnicode 196 0 R
+ /DescendantFonts [197 0 R]
+>>
+endobj
+197 0 obj
+<< /Type /Font
+/BaseFont /EAAAAA+mwa_cmr10
+/CIDToGIDMap /Identity
+/Subtype /CIDFontType2
+/CIDSystemInfo << /Registry (Adobe) /Ordering (UCS) /Supplement 0 >>
+/FontDescriptor 190 0 R
+/DW 0
+/W [ 0 [365 500 500 500 500 651 391 443 526 555 555 443 526 276 750 443 276 916 500 500 276 388 555 776 276 500 ] ]
+>>
+endobj
+196 0 obj
+<< /Length 198 0 R /Filter /FlateDecode >>
+stream
+x]j0E
+-Eg!PR
+^A~,SC- ^+Tܙh&VO&A48Vk Oȿe/.r2NW X
+MP2}q ?]߬"]mbGZot}U" P)'K>
+endobj
+200 0 obj
+<<
+ /Length1 1968
+ /Length 202 0 R
+ /Filter /FlateDecode
+>>
+stream
+xU]hWdj~VlTdc.JI5%[ddv3l\@%VRP|"E&yhb
+Eh+BK>(g&Z[sιw= w#ʡOfH#9iodp?*M7ߚ.,L|+?Js9NFHޞ#EV ɛrŹvd䪂5I\8Yێ~di([fwHpgJVyn-4 c/͚%KF = =+6q*>yx '2Q a'؆HU_um8"oBYG a1>2.RO]Fq#1M
ֆ낵5U~( LB'emꈪ,wȼhs"LJUNUY} vMJI9#k4Jˆ]J&GF#Ԗi u.r}аqp4bK#94*w[;lv#̙bw)*9*_Ӳh$m^7㚪 J6/kˌ.2lAvVb !y- hZݙ-Q0YۜفΰHU!+W4Jb9$%i5f
KOUȤtJ7UUU:}πs.H.?ڒ-2GGTIsaL >nM # EB3z)$|TXnnt1_
;β\;r70ŭF$$KS&; =վvn~;u Io{Q9 FG\ ѽ(
]|?+ObTSa 0Oz Q˾i&tr0Ca }a BՋY1 ho8Uz`PMɣ@x1zV LT\]?i'leW2i4)FYM=.w$Pf<'+$ib$J1PVia6? R>
+stream
+xk
+endstream
+endobj
+203 0 obj
+9
+endobj
+204 0 obj
+<<
+ /Type /Font
+ /Subtype /Type0
+ /BaseFont /EAAAAB+mwb_cmsy10
+ /Encoding /Identity-H
+ /ToUnicode 205 0 R
+ /DescendantFonts [206 0 R]
+>>
+endobj
+206 0 obj
+<< /Type /Font
+/BaseFont /EAAAAB+mwb_cmsy10
+/CIDToGIDMap /Identity
+/Subtype /CIDFontType2
+/CIDSystemInfo << /Registry (Adobe) /Ordering (UCS) /Supplement 0 >>
+/FontDescriptor 199 0 R
+/DW 0
+/W [ 0 [750 776 ] ]
+>>
+endobj
+205 0 obj
+<< /Length 207 0 R /Filter /FlateDecode >>
+stream
+x]Pj0+ l C ԒX}q z0h[lE^X gBpjV &
+:G:g}
+endstream
+endobj
+207 0 obj
+234
+endobj
+208 0 obj
+<<
+ /Type /FontDescriptor
+ /FontName /EAAAAC+mwb_cmmi10
+ /FontBBox [16 -215 882 715]
+ /Flags 33
+ /CapHeight 0
+ /Ascent 715
+ /Descent -215
+ /ItalicAngle 0
+ /StemV 0
+ /MissingWidth 500
+ /FontFile2 209 0 R
+ /CIDSet 210 0 R
+>>
+endobj
+209 0 obj
+<<
+ /Length1 2240
+ /Length 211 0 R
+ /Filter /FlateDecode
+>>
+stream
+xU{lSU{iF.clh
-SBHڭvV6lQ45rY
! >PH$HLԄ5wߧ9~ ADh@-U +?᱑zG'
] ?@ȖGLWȦ|cc=)s|l<1;6WTF[/?&iϏ3'.a/Y Sx G(!܋ G?"؊[
kЅqcIǶscg _)~IWJ_KKҀX/s!6q_hZ ;ֲF\%_ sxiSxԽxIK#0i}n L";,Kʷc )&R~7Z\[S]r.pm\,{XQݪw
2cVꌲ?TЦznlЀ`b&I=nˆVt9~"uS; f'i؛u;o5ΨҮЮ\ܞf87k%x/UUToӲhL3݆S4-2a%<Ѐ!)ZQnSFH-&4zXFLRh^M4\g510{Y>!3b0]
OP&~^=z75m*z:!/2 v4<
6ɃRiF;M(]$zrܮ$8~[U#eq^zօ*ĮjZR&5 *ka+Qȗ
+U
MpDYSjxuuV-jnfkVe{`AMiZRZk|l=s*M6f1;?yg#BTrW zp
+3_ :X5 L:?urϞcC#1G }yjc@_dE. ^dGcO;BB/=a_Z2 "$X|xmwfߢu?R-kcwfM`v-zeෂyZe3t:yO܆ńP]XhT[9"`g;X$%GY6_L*%kΊ21C
X@%X"cE,9Fs#1oHaILޮ$o+~IϥLUN
+BfG4E(*H sF:[β'(Gqfo2TK-FLȓFL[L Q)G4;v`<5W[e05=J˛}r>K&>yc:-Gi99%o?雳}ܓ p6N M^yK"H%`?eH#pP&67B4-bdZ 9O *>YZj >?n`ZxV1 *.n25ZYahy#*EjE%{~j
+endstream
+endobj
+211 0 obj
+1548
+endobj
+210 0 obj
+<< /Length 212 0 R /Filter /FlateDecode >>
+stream
+xk`p
+endstream
+endobj
+212 0 obj
+11
+endobj
+213 0 obj
+<<
+ /Type /Font
+ /Subtype /Type0
+ /BaseFont /EAAAAC+mwb_cmmi10
+ /Encoding /Identity-H
+ /ToUnicode 214 0 R
+ /DescendantFonts [215 0 R]
+>>
+endobj
+215 0 obj
+<< /Type /Font
+/BaseFont /EAAAAC+mwb_cmmi10
+/CIDToGIDMap /Identity
+/Subtype /CIDFontType2
+/CIDSystemInfo << /Registry (Adobe) /Ordering (UCS) /Supplement 0 >>
+/FontDescriptor 208 0 R
+/DW 0
+/W [ 0 [750 639 ] ]
+>>
+endobj
+214 0 obj
+<< /Length 216 0 R /Filter /FlateDecode >>
+stream
+x]Pj0+ l ťCP' K+GPKb- ƅ.13ð:g#ycc&B
+aѺA[wķd(D qpXz ;@| iӫF_h/!.B=t_~)'>N'95 Ԍ- t#2UTMv77`S6FUBeY+쒷qDT QJX9uxl5U|~}@
+endstream
+endobj
+216 0 obj
+235
+endobj
+217 0 obj
+<<
+ /Type /FontDescriptor
+ /FontName /EAAAAD+mwa_cmmi10
+ /FontBBox [-34 -250 1047 750]
+ /Flags 33
+ /CapHeight 683
+ /Ascent 750
+ /Descent -250
+ /ItalicAngle 0
+ /StemV 0
+ /MissingWidth 500
+ /FontFile2 218 0 R
+ /CIDSet 219 0 R
+>>
+endobj
+218 0 obj
+<<
+ /Length1 3164
+ /Length 220 0 R
+ /Filter /FlateDecode
+>>
+stream
+xV{lS}8ľvl6)8N@Ʌ8qC CU TZZt!&h(VmRQ:m4EcM+ :&ev=ν==@ Gey@{gɚC@F'r2t,{pX@
+\Τ呱M=Wu"$?d~k@!>Psr܃rU6moeks ar.'Σ|:=,Xa:>1؞wM܆8u hy>βg0gy1,<߇ x6BH A'iۚj:w[`b~~^d> fY Gf m}B ,姚r#&z'K
+ ?a ށ 98 o||qc p+z> ;ЇwCXɚ6l2V2NXx a~K.zZ (I])A$C^)T]~]
+EB7JBG{Aa.YY緥Y;9Q(.Y-)Cak ֔wĵ=_o] ڙlq/lo?ǰ`T':!'%%(KZ"I. PTq%. ENb*gpPa] .i\\bhBW&)z/ HĒE\
+nLX|[J^ Jn DO^ E³xN >
+Pc"lI [pc[)H4Kc/ijtIJ(4W ;?I2p5 4Fv`C<
+ aԵHm:eSOquՖr&Gd# @6Ok>
+(ۚprx 96rk8
<'j L\'?l%\Ha L;gh-TpzJ Ǝ0MYFkVnI`;,s gЄ`+xI[QD
+ϯ(uUgCD5Q)GLJ&\AKX@}G*z){d%KHpŌ3]35
+
+_躯
+"m|4N[^Pm5 }`r ,B_\E^Å+tut
+Ō4h M1;emxTpT?n
7^["}
nmn?pGGiC έ#)U)Mu0UZ?KI xv@ktyYyL xk?~'n8uIg@3t'==ßusWvxj6iiΑ 9 `}a}xkC"|@kik-s7i9VtK5>Ta 'bZgU8I4ez=ZFt=s?jxe\%yb
7;`kx=/2i5QHqb=FhUF5 :_5TMQ
aj0z/oh[$\\ I6Z}>iuVL&a%'Ep(XW`£C
-_gٳ h[
UG_*vZ@KP%M| c]+ 4I 9(Fn/6
+A%^n<~ VԓH S-%L#/af[¸q&at(])j%LO:KF~ćJZ s_a 0ci`
mE,q" ɢjau"X6Uاq|V ,UKszZJcі:zW\)#"jv889>$@{{WLz'tf:mJHO'g#@P{Нߋ|6==3(9EЏ V$gݹx ^ݲ8m<- 1U"^cX,| (>۾)tw ZM3:Ls/x,J1X ѱWy< a
+endstream
+endobj
+220 0 obj
+2211
+endobj
+219 0 obj
+<< /Length 221 0 R /Filter /FlateDecode >>
+stream
+xk```Pb
+endstream
+endobj
+221 0 obj
+15
+endobj
+222 0 obj
+<<
+ /Type /Font
+ /Subtype /Type0
+ /BaseFont /EAAAAD+mwa_cmmi10
+ /Encoding /Identity-H
+ /ToUnicode 223 0 R
+ /DescendantFonts [224 0 R]
+>>
+endobj
+224 0 obj
+<< /Type /Font
+/BaseFont /EAAAAD+mwa_cmmi10
+/CIDToGIDMap /Identity
+/Subtype /CIDFontType2
+/CIDSystemInfo << /Registry (Adobe) /Ordering (UCS) /Supplement 0 >>
+/FontDescriptor 217 0 R
+/DW 0
+/W [ 0 [365 831 680 612 583 276 ] ]
+>>
+endobj
+223 0 obj
+<< /Length 225 0 R /Filter /FlateDecode >>
+stream
+x]Qn +| VRHTrCMD!_0iV%0cƃmX=uZy`Έ =JKY@pR:+JJ .fn3m8wz4P%\` x&̀ q7')=hZ3j9PK}tΝ $sH3>7P.R
H\,踞0>k:e
`BE@y^?&T&$8
+şhyMRt =sQ+u{S,Zȱ
+endstream
+endobj
+225 0 obj
+262
+endobj
+189 0 obj
+<< /Type /Pages
+/Count 1
+/Kids [187 0 R ] >>
+endobj
+226 0 obj
+<<
+ /Type /Catalog
+ /Pages 189 0 R
+ /Lang (x-unknown)
+>>
+endobj
+188 0 obj
+<<
+ /Font <<
+ /F405 195 0 R
+ /F407 204 0 R
+ /F408 213 0 R
+ /F409 222 0 R
+>>
+ /ProcSet [/PDF /ImageB /ImageC /Text]
+ /ExtGState <<
+ /GS1 2 0 R
+ /GS2 3 0 R
+ /GS3 4 0 R
+ /GS4 5 0 R
+ /GS5 6 0 R
+ /GS6 7 0 R
+ /GS7 8 0 R
+ /GS8 9 0 R
+ /GS9 10 0 R
+ /GS10 11 0 R
+ /GS11 12 0 R
+ /GS12 13 0 R
+ /GS13 14 0 R
+ /GS14 15 0 R
+ /GS15 16 0 R
+ /GS16 17 0 R
+ /GS17 18 0 R
+ /GS18 19 0 R
+ /GS19 20 0 R
+ /GS20 21 0 R
+ /GS21 22 0 R
+ /GS22 23 0 R
+ /GS23 24 0 R
+ /GS24 25 0 R
+ /GS25 26 0 R
+ /GS26 27 0 R
+ /GS27 28 0 R
+ /GS28 29 0 R
+ /GS29 30 0 R
+ /GS30 31 0 R
+ /GS31 32 0 R
+ /GS32 33 0 R
+ /GS33 34 0 R
+ /GS34 35 0 R
+ /GS35 36 0 R
+ /GS36 37 0 R
+ /GS37 38 0 R
+ /GS38 39 0 R
+ /GS39 40 0 R
+ /GS40 41 0 R
+ /GS41 42 0 R
+ /GS42 43 0 R
+ /GS43 44 0 R
+ /GS44 45 0 R
+ /GS45 46 0 R
+ /GS46 47 0 R
+ /GS47 48 0 R
+ /GS48 49 0 R
+ /GS49 50 0 R
+ /GS50 51 0 R
+ /GS51 52 0 R
+ /GS52 53 0 R
+ /GS53 54 0 R
+ /GS54 55 0 R
+ /GS55 56 0 R
+ /GS56 57 0 R
+ /GS57 58 0 R
+ /GS58 59 0 R
+ /GS59 60 0 R
+ /GS60 61 0 R
+ /GS61 62 0 R
+ /GS62 63 0 R
+ /GS63 64 0 R
+ /GS64 65 0 R
+ /GS65 66 0 R
+ /GS66 67 0 R
+ /GS67 68 0 R
+ /GS68 69 0 R
+ /GS69 70 0 R
+ /GS70 71 0 R
+ /GS71 72 0 R
+ /GS72 73 0 R
+ /GS73 74 0 R
+ /GS74 75 0 R
+ /GS75 76 0 R
+ /GS76 77 0 R
+ /GS77 78 0 R
+ /GS78 79 0 R
+ /GS79 80 0 R
+ /GS80 81 0 R
+ /GS81 82 0 R
+ /GS82 83 0 R
+ /GS83 84 0 R
+ /GS84 85 0 R
+ /GS85 86 0 R
+ /GS86 87 0 R
+ /GS87 88 0 R
+ /GS88 89 0 R
+ /GS89 90 0 R
+ /GS90 91 0 R
+ /GS91 92 0 R
+ /GS92 93 0 R
+ /GS93 94 0 R
+ /GS94 95 0 R
+ /GS95 96 0 R
+ /GS96 97 0 R
+ /GS97 98 0 R
+ /GS98 99 0 R
+ /GS99 100 0 R
+ /GS100 101 0 R
+ /GS101 102 0 R
+ /GS102 103 0 R
+ /GS103 104 0 R
+ /GS104 105 0 R
+ /GS105 106 0 R
+ /GS106 107 0 R
+ /GS107 108 0 R
+ /GS108 109 0 R
+ /GS109 110 0 R
+ /GS110 111 0 R
+ /GS111 112 0 R
+ /GS112 113 0 R
+ /GS113 114 0 R
+ /GS114 115 0 R
+ /GS115 116 0 R
+ /GS116 117 0 R
+ /GS117 118 0 R
+ /GS118 119 0 R
+ /GS119 120 0 R
+ /GS120 121 0 R
+ /GS121 122 0 R
+ /GS122 123 0 R
+ /GS123 124 0 R
+ /GS124 125 0 R
+ /GS125 126 0 R
+ /GS126 127 0 R
+ /GS127 128 0 R
+ /GS128 129 0 R
+ /GS129 130 0 R
+ /GS130 131 0 R
+ /GS131 132 0 R
+ /GS132 133 0 R
+ /GS133 134 0 R
+ /GS134 135 0 R
+ /GS135 136 0 R
+ /GS136 137 0 R
+ /GS137 138 0 R
+ /GS138 139 0 R
+ /GS139 140 0 R
+ /GS140 141 0 R
+ /GS141 142 0 R
+ /GS142 143 0 R
+ /GS143 144 0 R
+ /GS144 145 0 R
+ /GS145 146 0 R
+ /GS146 147 0 R
+ /GS147 148 0 R
+ /GS148 149 0 R
+ /GS149 150 0 R
+ /GS150 151 0 R
+ /GS151 152 0 R
+ /GS152 153 0 R
+ /GS153 154 0 R
+ /GS154 155 0 R
+ /GS155 156 0 R
+ /GS156 157 0 R
+ /GS157 158 0 R
+ /GS158 159 0 R
+ /GS159 160 0 R
+ /GS160 161 0 R
+ /GS161 162 0 R
+ /GS162 163 0 R
+ /GS163 164 0 R
+ /GS164 165 0 R
+ /GS165 166 0 R
+ /GS166 167 0 R
+ /GS167 168 0 R
+ /GS168 169 0 R
+ /GS169 170 0 R
+ /GS170 171 0 R
+ /GS171 172 0 R
+ /GS172 173 0 R
+ /GS173 174 0 R
+ /GS174 175 0 R
+ /GS175 176 0 R
+ /GS176 177 0 R
+ /GS177 178 0 R
+ /GS178 179 0 R
+ /GS179 180 0 R
+ /GS180 181 0 R
+ /GS181 182 0 R
+ /GS182 183 0 R
+>>
+>>
+endobj
+xref
+0 227
+0000000000 65535 f
+0000000015 00000 n
+0000000145 00000 n
+0000000197 00000 n
+0000000249 00000 n
+0000000301 00000 n
+0000000353 00000 n
+0000000405 00000 n
+0000000457 00000 n
+0000000509 00000 n
+0000000561 00000 n
+0000000614 00000 n
+0000000667 00000 n
+0000000720 00000 n
+0000000773 00000 n
+0000000826 00000 n
+0000000879 00000 n
+0000000932 00000 n
+0000000985 00000 n
+0000001038 00000 n
+0000001091 00000 n
+0000001144 00000 n
+0000001197 00000 n
+0000001250 00000 n
+0000001303 00000 n
+0000001356 00000 n
+0000001409 00000 n
+0000001463 00000 n
+0000001517 00000 n
+0000001571 00000 n
+0000001625 00000 n
+0000001678 00000 n
+0000001731 00000 n
+0000001784 00000 n
+0000001837 00000 n
+0000001890 00000 n
+0000001943 00000 n
+0000001996 00000 n
+0000002049 00000 n
+0000002102 00000 n
+0000002155 00000 n
+0000002208 00000 n
+0000002261 00000 n
+0000002314 00000 n
+0000002367 00000 n
+0000002420 00000 n
+0000002473 00000 n
+0000002526 00000 n
+0000002579 00000 n
+0000002632 00000 n
+0000002685 00000 n
+0000002738 00000 n
+0000002791 00000 n
+0000002844 00000 n
+0000002897 00000 n
+0000002950 00000 n
+0000003003 00000 n
+0000003056 00000 n
+0000003109 00000 n
+0000003163 00000 n
+0000003217 00000 n
+0000003271 00000 n
+0000003325 00000 n
+0000003378 00000 n
+0000003431 00000 n
+0000003484 00000 n
+0000003537 00000 n
+0000003590 00000 n
+0000003643 00000 n
+0000003696 00000 n
+0000003749 00000 n
+0000003802 00000 n
+0000003855 00000 n
+0000003908 00000 n
+0000003961 00000 n
+0000004014 00000 n
+0000004067 00000 n
+0000004120 00000 n
+0000004173 00000 n
+0000004226 00000 n
+0000004279 00000 n
+0000004332 00000 n
+0000004385 00000 n
+0000004438 00000 n
+0000004491 00000 n
+0000004544 00000 n
+0000004597 00000 n
+0000004650 00000 n
+0000004703 00000 n
+0000004756 00000 n
+0000004809 00000 n
+0000004862 00000 n
+0000004915 00000 n
+0000004968 00000 n
+0000005021 00000 n
+0000005074 00000 n
+0000005127 00000 n
+0000005180 00000 n
+0000005233 00000 n
+0000005286 00000 n
+0000005339 00000 n
+0000005393 00000 n
+0000005447 00000 n
+0000005501 00000 n
+0000005555 00000 n
+0000005609 00000 n
+0000005663 00000 n
+0000005717 00000 n
+0000005771 00000 n
+0000005825 00000 n
+0000005879 00000 n
+0000005933 00000 n
+0000005987 00000 n
+0000006041 00000 n
+0000006095 00000 n
+0000006149 00000 n
+0000006203 00000 n
+0000006257 00000 n
+0000006311 00000 n
+0000006365 00000 n
+0000006419 00000 n
+0000006473 00000 n
+0000006527 00000 n
+0000006581 00000 n
+0000006635 00000 n
+0000006689 00000 n
+0000006743 00000 n
+0000006797 00000 n
+0000006851 00000 n
+0000006905 00000 n
+0000006959 00000 n
+0000007013 00000 n
+0000007067 00000 n
+0000007121 00000 n
+0000007175 00000 n
+0000007229 00000 n
+0000007283 00000 n
+0000007337 00000 n
+0000007391 00000 n
+0000007445 00000 n
+0000007499 00000 n
+0000007553 00000 n
+0000007607 00000 n
+0000007661 00000 n
+0000007715 00000 n
+0000007769 00000 n
+0000007823 00000 n
+0000007877 00000 n
+0000007931 00000 n
+0000007985 00000 n
+0000008039 00000 n
+0000008093 00000 n
+0000008147 00000 n
+0000008201 00000 n
+0000008255 00000 n
+0000008309 00000 n
+0000008363 00000 n
+0000008417 00000 n
+0000008471 00000 n
+0000008525 00000 n
+0000008579 00000 n
+0000008633 00000 n
+0000008687 00000 n
+0000008741 00000 n
+0000008795 00000 n
+0000008849 00000 n
+0000008903 00000 n
+0000008957 00000 n
+0000009011 00000 n
+0000009065 00000 n
+0000009119 00000 n
+0000009173 00000 n
+0000009227 00000 n
+0000009281 00000 n
+0000009335 00000 n
+0000009389 00000 n
+0000009443 00000 n
+0000009497 00000 n
+0000009551 00000 n
+0000009605 00000 n
+0000009659 00000 n
+0000009713 00000 n
+0000009767 00000 n
+0000009821 00000 n
+0000009875 00000 n
+0000120942 00000 n
+0000120966 00000 n
+0000120993 00000 n
+0000136192 00000 n
+0000136053 00000 n
+0000121191 00000 n
+0000121446 00000 n
+0000126286 00000 n
+0000126264 00000 n
+0000126384 00000 n
+0000126404 00000 n
+0000126889 00000 n
+0000126563 00000 n
+0000127316 00000 n
+0000127337 00000 n
+0000127589 00000 n
+0000129001 00000 n
+0000128979 00000 n
+0000129088 00000 n
+0000129107 00000 n
+0000129498 00000 n
+0000129267 00000 n
+0000129810 00000 n
+0000129831 00000 n
+0000130083 00000 n
+0000131751 00000 n
+0000131729 00000 n
+0000131840 00000 n
+0000131860 00000 n
+0000132251 00000 n
+0000132020 00000 n
+0000132564 00000 n
+0000132585 00000 n
+0000132841 00000 n
+0000135172 00000 n
+0000135150 00000 n
+0000135265 00000 n
+0000135285 00000 n
+0000135692 00000 n
+0000135445 00000 n
+0000136032 00000 n
+0000136115 00000 n
+trailer
+<<
+ /Root 226 0 R
+ /Info 1 0 R
+ /ID [ ]
+ /Size 227
+>>
+startxref
+139232
+%%EOF
diff --git a/figs/test_id31_cf_control_alpha.png b/figs/test_id31_cf_control_alpha.png
new file mode 100644
index 0000000..4ee6495
Binary files /dev/null and b/figs/test_id31_cf_control_alpha.png differ
diff --git a/figs/test_id31_cf_control_dy_dz_diff.pdf b/figs/test_id31_cf_control_dy_dz_diff.pdf
new file mode 100644
index 0000000..19a8a7f
Binary files /dev/null and b/figs/test_id31_cf_control_dy_dz_diff.pdf differ
diff --git a/figs/test_id31_cf_control_dy_dz_diff.png b/figs/test_id31_cf_control_dy_dz_diff.png
new file mode 100644
index 0000000..ff276bf
Binary files /dev/null and b/figs/test_id31_cf_control_dy_dz_diff.png differ
diff --git a/figs/test_id31_high_bandwidth_S.pdf b/figs/test_id31_high_bandwidth_S.pdf
new file mode 100644
index 0000000..1862f4e
Binary files /dev/null and b/figs/test_id31_high_bandwidth_S.pdf differ
diff --git a/figs/test_id31_high_bandwidth_S.png b/figs/test_id31_high_bandwidth_S.png
new file mode 100644
index 0000000..f873368
Binary files /dev/null and b/figs/test_id31_high_bandwidth_S.png differ
diff --git a/figs/test_id31_high_bandwidth_T.pdf b/figs/test_id31_high_bandwidth_T.pdf
new file mode 100644
index 0000000..d8bf519
Binary files /dev/null and b/figs/test_id31_high_bandwidth_T.pdf differ
diff --git a/figs/test_id31_high_bandwidth_T.png b/figs/test_id31_high_bandwidth_T.png
new file mode 100644
index 0000000..57092df
Binary files /dev/null and b/figs/test_id31_high_bandwidth_T.png differ
diff --git a/matlab/mat/2023-08-21_14-09_m0_cf_inv_tustin.mat b/matlab/mat/2023-08-21_14-09_m0_cf_inv_tustin.mat
new file mode 100644
index 0000000..f15d1f1
Binary files /dev/null and b/matlab/mat/2023-08-21_14-09_m0_cf_inv_tustin.mat differ
diff --git a/matlab/mat/2023-08-21_14-19_m0_cf_inv_tustin_Dz_60Hz_a_20.mat b/matlab/mat/2023-08-21_14-19_m0_cf_inv_tustin_Dz_60Hz_a_20.mat
new file mode 100644
index 0000000..312d749
Binary files /dev/null and b/matlab/mat/2023-08-21_14-19_m0_cf_inv_tustin_Dz_60Hz_a_20.mat differ
diff --git a/matlab/mat/2023-08-21_14-20_m0_cf_inv_tustin_Dy_60Hz_a_20.mat b/matlab/mat/2023-08-21_14-20_m0_cf_inv_tustin_Dy_60Hz_a_20.mat
new file mode 100644
index 0000000..1a90e71
Binary files /dev/null and b/matlab/mat/2023-08-21_14-20_m0_cf_inv_tustin_Dy_60Hz_a_20.mat differ
diff --git a/matlab/mat/2023-08-21_17-56_m2_cf_10Hz_a2.mat b/matlab/mat/2023-08-21_17-56_m2_cf_10Hz_a2.mat
new file mode 100644
index 0000000..6e5cc0f
Binary files /dev/null and b/matlab/mat/2023-08-21_17-56_m2_cf_10Hz_a2.mat differ
diff --git a/matlab/mat/2023-08-21_17-58_m2_cf_10Hz_a05.mat b/matlab/mat/2023-08-21_17-58_m2_cf_10Hz_a05.mat
new file mode 100644
index 0000000..a9c48e6
Binary files /dev/null and b/matlab/mat/2023-08-21_17-58_m2_cf_10Hz_a05.mat differ
diff --git a/matlab/test_id31_5_experiments.m b/matlab/test_id31_5_experiments.m
index 8441bfc..a25a1c4 100644
--- a/matlab/test_id31_5_experiments.m
+++ b/matlab/test_id31_5_experiments.m
@@ -934,6 +934,129 @@ data_dt_1000ums.Dy_rms_cl = rms(detrend(data_dt_1000ums.Dy_int(i_dt_1000ums)-dat
data_dt_1000ums.Dz_rms_cl = rms(detrend(data_dt_1000ums.Dz_int(i_dt_1000ums), 0));
data_dt_1000ums.Ry_rms_cl = rms(detrend(data_dt_1000ums.Ry_int(i_dt_1000ums), 0));
+%% Complementary Filter Control - Experimental Validation
+% Parameters to analyze results in the frequency domain
+Ts = 1e-4; % Sampling Time [s]
+Nfft = floor(2.0/Ts);
+win = hanning(Nfft);
+Noverlap = floor(Nfft/2);
+
+%% Closed-Loop performance identification (S and T transfer functions)
+% Different parameters w0 in different directions
+data = load("2023-08-21_14-09_m0_cf_inv_tustin.mat");
+
+[S_dy, f] = tfestimate(data.Dy.id_cl, data.Dy.e_dy, win, Noverlap, Nfft, 1/Ts);
+[T_dy, ~] = tfestimate(data.Dy.id_cl, data.Dy.Dy_int, win, Noverlap, Nfft, 1/Ts);
+
+[S_dz, ~] = tfestimate(data.Dz.id_cl, data.Dz.e_dz, win, Noverlap, Nfft, 1/Ts);
+[T_dz, ~] = tfestimate(data.Dz.id_cl, data.Dz.Dz_int, win, Noverlap, Nfft, 1/Ts);
+
+[S_ry, ~] = tfestimate(data.Ry.id_cl, data.Ry.e_ry, win, Noverlap, Nfft, 1/Ts);
+[T_ry, ~] = tfestimate(data.Ry.id_cl, data.Ry.Ry_int, win, Noverlap, Nfft, 1/Ts);
+
+%% Obtained Closed-Loop transfer functions - Different Complementary filters for Dy and Dz
+figure;
+hold on;
+plot(f, abs(S_dy), '-', 'color', colors(1,:),'DisplayName', '$S_{D_y}$');
+plot(f, abs(S_dz), '-', 'color', colors(2,:),'DisplayName', '$S_{D_z}$');
+plot(f, abs(S_ry), '-', 'color', colors(3,:),'DisplayName', '$S_{R_y}$');
+hold off;
+set(gca, 'XScale', 'log'); set(gca, 'YScale', 'log');
+xlabel('Frequency [Hz]'); ylabel('Magnitude');
+ylim([1e-3, 5]);
+xlim([1, 1e3]);
+leg = legend('location', 'southeast', 'FontSize', 8, 'NumColumns', 1);
+leg.ItemTokenSize(1) = 18;
+
+%% Effect of a change of complementary filter parameter alpha
+% w0 = 10, alpha = 2
+data_a2 = load("2023-08-21_17-56_m2_cf_10Hz_a2.mat");
+
+[S_a2, ~] = tfestimate(data_a2.Dy.id_cl, data_a2.Dy.e_dy, win, Noverlap, Nfft, 1/Ts);
+[T_a2, ~] = tfestimate(data_a2.Dy.id_cl, data_a2.Dy.Dy_int, win, Noverlap, Nfft, 1/Ts);
+
+% Analytical formulas
+w0 = 2*pi*10;
+alpha = 2;
+Hh_a2 = (s/w0)^2*((s/w0)+1+alpha)/(((s/w0)+1)*((s/w0)^2 + alpha*(s/w0) + 1));
+Hl_a2 = ((1+alpha)*(s/w0)+1)/(((s/w0)+1)*((s/w0)^2 + alpha*(s/w0) + 1));
+
+% w0 = 10, alpha = 0.5
+data_a05 = load("2023-08-21_17-58_m2_cf_10Hz_a05.mat");
+
+[S_a05, ~] = tfestimate(data_a05.Dz.id_cl, data_a05.Dz.e_dz, win, Noverlap, Nfft, 1/Ts);
+[T_a05, ~] = tfestimate(data_a05.Dz.id_cl, data_a05.Dz.Dz_int, win, Noverlap, Nfft, 1/Ts);
+
+% Analytical formulas
+w0 = 2*pi*10;
+alpha = 0.5;
+Hh_a05 = (s/w0)^2*((s/w0)+1+alpha)/(((s/w0)+1)*((s/w0)^2 + alpha*(s/w0) + 1));
+Hl_a05 = ((1+alpha)*(s/w0)+1)/(((s/w0)+1)*((s/w0)^2 + alpha*(s/w0) + 1));
+
+%% Obtained Closed-Loop transfer functions - Effect of changing alpha
+figure;
+hold on;
+plot(freqs, abs(squeeze(freqresp(Hh_a2, freqs, 'Hz'))), '--', 'color', colors(1,:),'DisplayName', '$\alpha = 2$, $H_H$, $H_L$');
+plot(f, abs(S_a2), '-', 'color', [colors(1,:), 0.5], 'linewidth', 2.5,'DisplayName', '$\alpha = 2$, $S$, $T$');
+plot(freqs, abs(squeeze(freqresp(Hh_a05, freqs, 'Hz'))), '--', 'color', colors(2,:),'DisplayName', '$\alpha = 0.5$, $H_H$, $H_L$');
+plot(f, abs(S_a05), '-', 'color', [colors(2,:), 0.5], 'linewidth', 2.5,'DisplayName', '$\alpha = 0.5$, $S$, $T$');
+plot(freqs, abs(squeeze(freqresp(Hl_a2, freqs, 'Hz'))), '--', 'color', colors(1,:), 'HandleVisibility', 'off');
+plot(f, abs(T_a2), '-', 'color', [colors(1,:), 0.5], 'linewidth', 2.5, 'HandleVisibility', 'off');
+plot(freqs, abs(squeeze(freqresp(Hl_a05, freqs, 'Hz'))), '--', 'color', colors(2,:), 'HandleVisibility', 'off');
+plot(f, abs(T_a05), '-', 'color', [colors(2,:), 0.5], 'linewidth', 2.5, 'HandleVisibility', 'off');
+hold off;
+set(gca, 'XScale', 'log'); set(gca, 'YScale', 'log');
+xlabel('Frequency [Hz]'); ylabel('Magnitude');
+xlim([1, 1e3]);
+ylim([1e-3, 5]);
+leg = legend('location', 'south', 'FontSize', 8, 'NumColumns', 2);
+leg.ItemTokenSize(1) = 18;
+
+%% Trying high bandwidth controllers - 60Hz
+% Z direction
+data = load("2023-08-21_14-19_m0_cf_inv_tustin_Dz_60Hz_a_20.mat");
+[S_dz, ~] = tfestimate(data.Dz.id_cl, data.Dz.e_dz, win, Noverlap, Nfft, 1/Ts);
+[T_dz, ~] = tfestimate(data.Dz.id_cl, data.Dz.Dz_int, win, Noverlap, Nfft, 1/Ts);
+
+% Y direction
+data = load("2023-08-21_14-20_m0_cf_inv_tustin_Dy_60Hz_a_20.mat");
+[S_dy, ~] = tfestimate(data.Dy.id_cl, data.Dy.e_dy, win, Noverlap, Nfft, 1/Ts);
+[T_dy, ~] = tfestimate(data.Dy.id_cl, data.Dy.Dy_int, win, Noverlap, Nfft, 1/Ts);
+
+% Analytical formulas
+w0 = 2*pi*60;
+alpha = 2;
+Hh_60Hz = (s/w0)^2*((s/w0)+1+alpha)/(((s/w0)+1)*((s/w0)^2 + alpha*(s/w0) + 1));
+Hl_60Hz = ((1+alpha)*(s/w0)+1)/(((s/w0)+1)*((s/w0)^2 + alpha*(s/w0) + 1));
+
+%% Bode plot of the Weighting filters and Obtained complementary filters
+figure;
+hold on;
+plot(f, abs(S_dy), '-', 'color', colors(1,:), 'DisplayName', '$S$, $D_y$');
+plot(f, abs(S_dz), '-', 'color', colors(2,:), 'DisplayName', '$S$, $D_z$');
+plot(freqs, abs(squeeze(freqresp(Hh_60Hz, freqs, 'Hz'))), 'k--','DisplayName', '$H_H$');
+hold off;
+set(gca, 'XScale', 'log'); set(gca, 'YScale', 'log');
+xlabel('Frequency [Hz]'); ylabel('Magnitude');
+xlim([1, 1e3]);
+ylim([3e-4, 5]);
+leg = legend('location', 'southeast', 'FontSize', 8, 'NumColumns', 1);
+leg.ItemTokenSize(1) = 18;
+
+%% Bode plot of the Weighting filters and Obtained complementary filters
+figure;
+hold on;
+plot(f, abs(T_dy), '-', 'color', colors(1,:), 'DisplayName', '$T$, $D_y$');
+plot(f, abs(T_dz), '-', 'color', colors(2,:), 'DisplayName', '$T$, $D_z$');
+plot(freqs, abs(squeeze(freqresp(Hl_60Hz, freqs, 'Hz'))), 'k--','DisplayName', '$H_L$');
+hold off;
+set(gca, 'XScale', 'log'); set(gca, 'YScale', 'log');
+xlabel('Frequency [Hz]'); ylabel('Magnitude');
+xlim([1, 1e3]);
+ylim([3e-4, 5]);
+leg = legend('location', 'southwest', 'FontSize', 8, 'NumColumns', 1);
+leg.ItemTokenSize(1) = 18;
+
%% Summary of results
% 1e9*data_tomo_m0_Wz6.Dy_rms_ol, 1e9*data_tomo_m0_Wz6.Dz_rms_ol, 1e6*data_tomo_m0_Wz6.Ry_rms_ol; % Tomo - OL - 6deg/s - 0kg
% 1e9*data_tomo_m0_Wz6.Dy_rms_cl, 1e9*data_tomo_m0_Wz6.Dz_rms_cl, 1e6*data_tomo_m0_Wz6.Ry_rms_cl; % Tomo - CL - 6deg/s - 0kg
diff --git a/test-bench-id31.org b/test-bench-id31.org
index 2f27ced..1eaa40b 100644
--- a/test-bench-id31.org
+++ b/test-bench-id31.org
@@ -3529,6 +3529,7 @@ Several scientific experiments were replicated, such as:
- Diffraction Tomography:continuous $R_z$ rotation using the Spindle and lateral $D_y$ scans performed at the same time using the translation stage (Section ref:ssec:test_id31_scans_diffraction_tomo)
Unless explicitly stated, all closed-loop experiments were performed using the robust (i.e. conservative) high authority controller designed in Section ref:ssec:test_id31_iff_hac_controller.
+Higher performance controllers using complementary filters are investigated in Section ref:ssec:test_id31_cf_control.
For each experiment, the obtained performances are compared to the specifications for the most demanding case in which nano-focusing optics are used to focus the beam down to $200\,nm\times 100\,nm$.
In this case, the goal is to keep the sample's point of interest in the beam, and therefore the $D_y$ and $D_z$ positioning errors should be less than $200\,nm$ and $100\,nm$ peak-to-peak, respectively.
@@ -3942,7 +3943,7 @@ data_tomo_m0_Wz180.Ry_rms_ol = rms(data_tomo_m0_Wz180.Ry_int(1:i_m0) - (y0 + R*s
A comparative analysis was conducted using three tomography scans at $180\,\text{deg/s}$ to evaluate the effectiveness of the HAC-LAC strategy in reducing positioning errors.
The scans were performed under three conditions: open-loop, with decentralized IFF control, and with the complete HAC-LAC strategy.
-For these specific measurements, an enhanced high authority controller was optimized for low payload masses to meet the performance requirements.
+For this specific measurement, an enhanced high authority controller (discussed in Section ref:ssec:test_id31_cf_control) was optimized for low payload masses to meet the performance requirements.
Figure ref:fig:test_id31_hac_cas_cl presents the cumulative amplitude spectra of the position errors for all three cases.
The results reveal two distinct control contributions: the decentralized IFF effectively attenuates vibrations near the nano-hexapod suspension modes (an achievement not possible with HAC alone), while the high authority controller suppresses low-frequency vibrations primarily arising from Spindle guiding errors.
@@ -5041,6 +5042,236 @@ data_dt_1000ums.Dz_rms_cl = rms(detrend(data_dt_1000ums.Dz_int(i_dt_1000ums), 0)
data_dt_1000ums.Ry_rms_cl = rms(detrend(data_dt_1000ums.Ry_int(i_dt_1000ums), 0));
#+end_src
+** Feedback control using Complementary Filters
+<>
+
+# TODO - Add link to section
+A control architecture utilizing complementary filters to shape the closed-loop transfer functions was proposed during the detail design phase.
+Experimental validation of this architecture using the NASS is presented herein.
+
+Given that performance requirements are specified in the Cartesian frame, decoupling of the plant within this frame was achieved using Jacobian matrices.
+Consequently, the control space comprises the directions $D_x$, $D_y$, $D_z$, $R_x$, and $R_y$.
+Control performance in each of these directions can be tuned independently.
+A schematic of the proposed control architecture is illustrated in Figure\nbsp{}ref:fig:test_id31_cf_control.
+
+#+name: fig:test_id31_cf_control
+#+caption: Control architecture in the Cartesian frame. Only the controller corresponding to the $D_z$ direction is shown. $H_L$ and $H_H$ are complementary filters.
+[[file:figs/test_id31_cf_control.png]]
+
+
+# TODO - Add link to 2DoF model
+Implementation of this control architecture necessitates a plant model, which must subsequently be inverted.
+This plant model was derived from the multi-body model incorporating the previously detailed 2-DoF APA model, such that the model order stays relatively low.
+Proposed analytical formulas for complementary filters having $40\,\text{dB/dec}$ were used during this experimental validation.
+# TODO - Add link to the analytical formulas
+
+An initial experimental validation was conducted under no-payload conditions, with control applied solely to the $D_y$, $D_z$, and $R_y$ directions.
+Increased control bandwidth was achieved for the $D_z$ and $R_y$ directions through appropriate tuning of the parameter $\omega_0$.
+The experimentally measured closed-loop sensitivity transfer functions corresponding to these three controlled directions are presented in Figure ref:fig:test_id31_cf_control_dy_dz_diff.
+
+#+begin_src matlab
+%% Complementary Filter Control - Experimental Validation
+% Parameters to analyze results in the frequency domain
+Ts = 1e-4; % Sampling Time [s]
+Nfft = floor(2.0/Ts);
+win = hanning(Nfft);
+Noverlap = floor(Nfft/2);
+
+%% Closed-Loop performance identification (S and T transfer functions)
+% Different parameters w0 in different directions
+data = load("2023-08-21_14-09_m0_cf_inv_tustin.mat");
+
+[S_dy, f] = tfestimate(data.Dy.id_cl, data.Dy.e_dy, win, Noverlap, Nfft, 1/Ts);
+[T_dy, ~] = tfestimate(data.Dy.id_cl, data.Dy.Dy_int, win, Noverlap, Nfft, 1/Ts);
+
+[S_dz, ~] = tfestimate(data.Dz.id_cl, data.Dz.e_dz, win, Noverlap, Nfft, 1/Ts);
+[T_dz, ~] = tfestimate(data.Dz.id_cl, data.Dz.Dz_int, win, Noverlap, Nfft, 1/Ts);
+
+[S_ry, ~] = tfestimate(data.Ry.id_cl, data.Ry.e_ry, win, Noverlap, Nfft, 1/Ts);
+[T_ry, ~] = tfestimate(data.Ry.id_cl, data.Ry.Ry_int, win, Noverlap, Nfft, 1/Ts);
+#+end_src
+
+#+begin_src matlab :exports none :results none
+%% Obtained Closed-Loop transfer functions - Different Complementary filters for Dy and Dz
+figure;
+hold on;
+plot(f, abs(S_dy), '-', 'color', colors(1,:),'DisplayName', '$S_{D_y}$');
+plot(f, abs(S_dz), '-', 'color', colors(2,:),'DisplayName', '$S_{D_z}$');
+plot(f, abs(S_ry), '-', 'color', colors(3,:),'DisplayName', '$S_{R_y}$');
+hold off;
+set(gca, 'XScale', 'log'); set(gca, 'YScale', 'log');
+xlabel('Frequency [Hz]'); ylabel('Magnitude');
+ylim([1e-3, 5]);
+xlim([1, 1e3]);
+leg = legend('location', 'southeast', 'FontSize', 8, 'NumColumns', 1);
+leg.ItemTokenSize(1) = 18;
+#+end_src
+
+#+begin_src matlab :tangle no :exports results :results file replace
+exportFig('figs/test_id31_cf_control_dy_dz_diff.pdf', 'width', 'half', 'height', 'normal');
+#+end_src
+
+Another test was conducted with a $26\,\text{kg}$ payload.
+For this configuration, complementary filters were implemented with $\omega_0 = 2\pi \cdot 10\,\text{rad/s}$, and parameter $\alpha$ was varied.
+The resulting experimentally obtained closed-loop transfer functions are compared against the theoretical complementary filter responses in Figure\nbsp{}ref:fig:test_id31_cf_control_alpha.
+As illustrated in the figure, a close correspondence between the measured closed-loop responses and the target complementary filter magnitude was observed.
+It also shows that the parameter $\alpha$ provides a mechanism for managing the trade-off between low-frequency disturbance rejection performance and the potential amplification of disturbances within the crossover frequency region.
+
+#+begin_src matlab
+%% Effect of a change of complementary filter parameter alpha
+% w0 = 10, alpha = 2
+data_a2 = load("2023-08-21_17-56_m2_cf_10Hz_a2.mat");
+
+[S_a2, ~] = tfestimate(data_a2.Dy.id_cl, data_a2.Dy.e_dy, win, Noverlap, Nfft, 1/Ts);
+[T_a2, ~] = tfestimate(data_a2.Dy.id_cl, data_a2.Dy.Dy_int, win, Noverlap, Nfft, 1/Ts);
+
+% Analytical formulas
+w0 = 2*pi*10;
+alpha = 2;
+Hh_a2 = (s/w0)^2*((s/w0)+1+alpha)/(((s/w0)+1)*((s/w0)^2 + alpha*(s/w0) + 1));
+Hl_a2 = ((1+alpha)*(s/w0)+1)/(((s/w0)+1)*((s/w0)^2 + alpha*(s/w0) + 1));
+
+% w0 = 10, alpha = 0.5
+data_a05 = load("2023-08-21_17-58_m2_cf_10Hz_a05.mat");
+
+[S_a05, ~] = tfestimate(data_a05.Dz.id_cl, data_a05.Dz.e_dz, win, Noverlap, Nfft, 1/Ts);
+[T_a05, ~] = tfestimate(data_a05.Dz.id_cl, data_a05.Dz.Dz_int, win, Noverlap, Nfft, 1/Ts);
+
+% Analytical formulas
+w0 = 2*pi*10;
+alpha = 0.5;
+Hh_a05 = (s/w0)^2*((s/w0)+1+alpha)/(((s/w0)+1)*((s/w0)^2 + alpha*(s/w0) + 1));
+Hl_a05 = ((1+alpha)*(s/w0)+1)/(((s/w0)+1)*((s/w0)^2 + alpha*(s/w0) + 1));
+#+end_src
+
+#+begin_src matlab :exports none :results none
+%% Obtained Closed-Loop transfer functions - Effect of changing alpha
+figure;
+hold on;
+plot(freqs, abs(squeeze(freqresp(Hh_a2, freqs, 'Hz'))), '--', 'color', colors(1,:),'DisplayName', '$\alpha = 2$, $H_H$, $H_L$');
+plot(f, abs(S_a2), '-', 'color', [colors(1,:), 0.5], 'linewidth', 2.5,'DisplayName', '$\alpha = 2$, $S$, $T$');
+plot(freqs, abs(squeeze(freqresp(Hh_a05, freqs, 'Hz'))), '--', 'color', colors(2,:),'DisplayName', '$\alpha = 0.5$, $H_H$, $H_L$');
+plot(f, abs(S_a05), '-', 'color', [colors(2,:), 0.5], 'linewidth', 2.5,'DisplayName', '$\alpha = 0.5$, $S$, $T$');
+plot(freqs, abs(squeeze(freqresp(Hl_a2, freqs, 'Hz'))), '--', 'color', colors(1,:), 'HandleVisibility', 'off');
+plot(f, abs(T_a2), '-', 'color', [colors(1,:), 0.5], 'linewidth', 2.5, 'HandleVisibility', 'off');
+plot(freqs, abs(squeeze(freqresp(Hl_a05, freqs, 'Hz'))), '--', 'color', colors(2,:), 'HandleVisibility', 'off');
+plot(f, abs(T_a05), '-', 'color', [colors(2,:), 0.5], 'linewidth', 2.5, 'HandleVisibility', 'off');
+hold off;
+set(gca, 'XScale', 'log'); set(gca, 'YScale', 'log');
+xlabel('Frequency [Hz]'); ylabel('Magnitude');
+xlim([1, 1e3]);
+ylim([1e-3, 5]);
+leg = legend('location', 'south', 'FontSize', 8, 'NumColumns', 2);
+leg.ItemTokenSize(1) = 18;
+#+end_src
+
+#+begin_src matlab :tangle no :exports results :results file replace
+exportFig('figs/test_id31_cf_control_alpha.pdf', 'width', 'half', 'height', 'normal');
+#+end_src
+
+#+name: fig:test_id31_cf_control_results
+#+caption: Measured closed-loop transfer functions. Different bandwidth can be specified for different directions using $\omega_0$ (\subref{fig:test_id31_cf_control_dy_dz_diff}). The shape can be adjusted using parameter $\alpha$ (\subref{fig:test_id31_cf_control_alpha}).
+#+attr_latex: :options [htbp]
+#+begin_figure
+#+attr_latex: :caption \subcaption{\label{fig:test_id31_cf_control_dy_dz_diff}Chose of bandwidth using $\omega_0$, $m = 0\,\text{kg}$}
+#+attr_latex: :options {0.49\textwidth}
+#+begin_subfigure
+#+attr_latex: :scale 0.9
+[[file:figs/test_id31_cf_control_dy_dz_diff.png]]
+#+end_subfigure
+#+attr_latex: :caption \subcaption{\label{fig:test_id31_cf_control_alpha}Effect of a change of $\alpha$, $m = 26\,\text{kg}$}
+#+attr_latex: :options {0.49\textwidth}
+#+begin_subfigure
+#+attr_latex: :scale 0.9
+[[file:figs/test_id31_cf_control_alpha.png]]
+#+end_subfigure
+#+end_figure
+
+Finally, $\omega_0$ was gradually increased to estimate the maximum bandwidth (i.e. the best low frequency disturbance rejection) that can be achieved with this architecture.
+No payload was used for this test, and the parameter $\omega_0$ was increased for the controllers in the $D_y$ and $D_z$ directions.
+A value $\omega_0 = 2\pi \cdot 60 \,\text{rad/s}$ could be achieved.
+Measured closed-loop transfer functions are shown in Figure ref:fig:test_id31_high_bandwidth, indicating a reduction of disturbances in the considered direction of $1000$ at $1\,\text{Hz}$.
+For higher values of $\omega_0$, the system became unstable in the vertical direction, probably because of the resonance at $250\,\text{Hz}$ that is not well captured with the multi-body model (Figure ref:fig:test_id31_hac_plant_effect_mass).
+
+#+begin_src matlab
+%% Trying high bandwidth controllers - 60Hz
+% Z direction
+data = load("2023-08-21_14-19_m0_cf_inv_tustin_Dz_60Hz_a_20.mat");
+[S_dz, ~] = tfestimate(data.Dz.id_cl, data.Dz.e_dz, win, Noverlap, Nfft, 1/Ts);
+[T_dz, ~] = tfestimate(data.Dz.id_cl, data.Dz.Dz_int, win, Noverlap, Nfft, 1/Ts);
+
+% Y direction
+data = load("2023-08-21_14-20_m0_cf_inv_tustin_Dy_60Hz_a_20.mat");
+[S_dy, ~] = tfestimate(data.Dy.id_cl, data.Dy.e_dy, win, Noverlap, Nfft, 1/Ts);
+[T_dy, ~] = tfestimate(data.Dy.id_cl, data.Dy.Dy_int, win, Noverlap, Nfft, 1/Ts);
+
+% Analytical formulas
+w0 = 2*pi*60;
+alpha = 2;
+Hh_60Hz = (s/w0)^2*((s/w0)+1+alpha)/(((s/w0)+1)*((s/w0)^2 + alpha*(s/w0) + 1));
+Hl_60Hz = ((1+alpha)*(s/w0)+1)/(((s/w0)+1)*((s/w0)^2 + alpha*(s/w0) + 1));
+#+end_src
+
+#+begin_src matlab :exports none :results none
+%% Bode plot of the Weighting filters and Obtained complementary filters
+figure;
+hold on;
+plot(f, abs(S_dy), '-', 'color', colors(1,:), 'DisplayName', '$S$, $D_y$');
+plot(f, abs(S_dz), '-', 'color', colors(2,:), 'DisplayName', '$S$, $D_z$');
+plot(freqs, abs(squeeze(freqresp(Hh_60Hz, freqs, 'Hz'))), 'k--','DisplayName', '$H_H$');
+hold off;
+set(gca, 'XScale', 'log'); set(gca, 'YScale', 'log');
+xlabel('Frequency [Hz]'); ylabel('Magnitude');
+xlim([1, 1e3]);
+ylim([3e-4, 5]);
+leg = legend('location', 'southeast', 'FontSize', 8, 'NumColumns', 1);
+leg.ItemTokenSize(1) = 18;
+#+end_src
+
+#+begin_src matlab :tangle no :exports results :results file replace
+yticks([1e-3, 1e-2, 1e-1, 1]);
+exportFig('figs/test_id31_high_bandwidth_S.pdf', 'width', 'half', 'height', 'normal');
+#+end_src
+
+#+begin_src matlab :exports none :results none
+%% Bode plot of the Weighting filters and Obtained complementary filters
+figure;
+hold on;
+plot(f, abs(T_dy), '-', 'color', colors(1,:), 'DisplayName', '$T$, $D_y$');
+plot(f, abs(T_dz), '-', 'color', colors(2,:), 'DisplayName', '$T$, $D_z$');
+plot(freqs, abs(squeeze(freqresp(Hl_60Hz, freqs, 'Hz'))), 'k--','DisplayName', '$H_L$');
+hold off;
+set(gca, 'XScale', 'log'); set(gca, 'YScale', 'log');
+xlabel('Frequency [Hz]'); ylabel('Magnitude');
+xlim([1, 1e3]);
+ylim([3e-4, 5]);
+leg = legend('location', 'southwest', 'FontSize', 8, 'NumColumns', 1);
+leg.ItemTokenSize(1) = 18;
+#+end_src
+
+#+begin_src matlab :tangle no :exports results :results file replace
+yticks([1e-3, 1e-2, 1e-1, 1]);
+exportFig('figs/test_id31_high_bandwidth_T.pdf', 'width', 'half', 'height', 'normal');
+#+end_src
+
+#+name: fig:test_id31_high_bandwidth
+#+caption: Measured Closed-Loop Sensitivity (\subref{fig:test_id31_high_bandwidth_S}) and Complementary Sensitivity (\subref{fig:test_id31_high_bandwidth_T}) transfer functions for the highest test bandwidth $\omega_0 = 2\pi\cdot 60\,\text{rad/s}$.
+#+attr_latex: :options [htbp]
+#+begin_figure
+#+attr_latex: :caption \subcaption{\label{fig:test_id31_high_bandwidth_S}Sensitivity}
+#+attr_latex: :options {0.49\textwidth}
+#+begin_subfigure
+#+attr_latex: :scale 0.9
+[[file:figs/test_id31_high_bandwidth_S.png]]
+#+end_subfigure
+#+attr_latex: :caption \subcaption{\label{fig:test_id31_high_bandwidth_T}Complementary Sensitivity}
+#+attr_latex: :options {0.49\textwidth}
+#+begin_subfigure
+#+attr_latex: :scale 0.9
+[[file:figs/test_id31_high_bandwidth_T.png]]
+#+end_subfigure
+#+end_figure
+
** Conclusion
:PROPERTIES:
:UNNUMBERED: t
diff --git a/test-bench-id31.pdf b/test-bench-id31.pdf
index 604e487..eeb3ec0 100644
Binary files a/test-bench-id31.pdf and b/test-bench-id31.pdf differ
diff --git a/test-bench-id31.tex b/test-bench-id31.tex
index aca3f8a..af261fd 100644
--- a/test-bench-id31.tex
+++ b/test-bench-id31.tex
@@ -1,4 +1,4 @@
-% Created 2025-04-07 Mon 14:38
+% Created 2025-04-20 Sun 22:21
% Intended LaTeX compiler: pdflatex
\documentclass[a4paper, 10pt, DIV=12, parskip=full, bibliography=totoc]{scrreprt}
@@ -215,6 +215,7 @@ Results are summarized in Table \ref{tab:test_id31_metrology_acceptance}.
The obtained lateral acceptance for pure displacements in any direction is estimated to be around \(+/-0.5\,mm\), which is enough for the current application as it is well above the micro-station errors to be actively corrected by the NASS.
\begin{table}[htbp]
+\caption{\label{tab:test_id31_metrology_acceptance}Estimated measurement range for each interferometer, and for three different directions.}
\centering
\begin{tabularx}{0.45\linewidth}{Xccc}
\toprule
@@ -227,8 +228,6 @@ The obtained lateral acceptance for pure displacements in any direction is estim
\(d_5\) (z) & \(1.33\, mm\) & \(1.06\,mm\) & \(>2\,mm\)\\
\bottomrule
\end{tabularx}
-\caption{\label{tab:test_id31_metrology_acceptance}Estimated measurement range for each interferometer, and for three different directions.}
-
\end{table}
\section{Estimated measurement errors}
\label{ssec:test_id31_metrology_errors}
@@ -739,6 +738,7 @@ Several scientific experiments were replicated, such as:
\end{itemize}
Unless explicitly stated, all closed-loop experiments were performed using the robust (i.e. conservative) high authority controller designed in Section \ref{ssec:test_id31_iff_hac_controller}.
+Higher performance controllers using complementary filters are investigated in Section \ref{ssec:test_id31_cf_control}.
For each experiment, the obtained performances are compared to the specifications for the most demanding case in which nano-focusing optics are used to focus the beam down to \(200\,nm\times 100\,nm\).
In this case, the goal is to keep the sample's point of interest in the beam, and therefore the \(D_y\) and \(D_z\) positioning errors should be less than \(200\,nm\) and \(100\,nm\) peak-to-peak, respectively.
@@ -748,6 +748,7 @@ In terms of RMS errors, this corresponds to \(30\,nm\) in \(D_y\), \(15\,nm\) in
Results obtained for all experiments are summarized and compared to the specifications in Section \ref{ssec:test_id31_scans_conclusion}.
\begin{table}[htbp]
+\caption{\label{tab:test_id31_experiments_specifications}Specifications for the Nano-Active-Stabilization-System}
\centering
\begin{tabularx}{0.45\linewidth}{Xccc}
\toprule
@@ -757,8 +758,6 @@ peak 2 peak & 200nm & 100nm & \(1.7\,\mu\text{rad}\)\\
RMS & 30nm & 15nm & \(250\,\text{nrad}\)\\
\bottomrule
\end{tabularx}
-\caption{\label{tab:test_id31_experiments_specifications}Specifications for the Nano-Active-Stabilization-System}
-
\end{table}
\section{Tomography Scans}
\label{ssec:test_id31_scans_tomography}
@@ -824,7 +823,7 @@ Nevertheless, even with this robust (i.e. conservative) HAC implementation, the
A comparative analysis was conducted using three tomography scans at \(180\,\text{deg/s}\) to evaluate the effectiveness of the HAC-LAC strategy in reducing positioning errors.
The scans were performed under three conditions: open-loop, with decentralized IFF control, and with the complete HAC-LAC strategy.
-For these specific measurements, an enhanced high authority controller was optimized for low payload masses to meet the performance requirements.
+For this specific measurement, an enhanced high authority controller (discussed in Section \ref{ssec:test_id31_cf_control}) was optimized for low payload masses to meet the performance requirements.
Figure \ref{fig:test_id31_hac_cas_cl} presents the cumulative amplitude spectra of the position errors for all three cases.
The results reveal two distinct control contributions: the decentralized IFF effectively attenuates vibrations near the nano-hexapod suspension modes (an achievement not possible with HAC alone), while the high authority controller suppresses low-frequency vibrations primarily arising from Spindle guiding errors.
@@ -1083,6 +1082,74 @@ Alternatively, a feedforward controller could improve the lateral positioning ac
\end{subfigure}
\caption{\label{fig:test_id31_diffraction_tomo}Diffraction tomography scans (combined \(R_z\) and \(D_y\) motions) at several \(D_y\) velocities (\(R_z\) rotational velocity is \(6\,\text{deg/s}\)).}
\end{figure}
+\section{Feedback control using Complementary Filters}
+\label{ssec:test_id31_cf_control}
+
+A control architecture utilizing complementary filters to shape the closed-loop transfer functions was proposed during the detail design phase.
+Experimental validation of this architecture using the NASS is presented herein.
+
+Given that performance requirements are specified in the Cartesian frame, decoupling of the plant within this frame was achieved using Jacobian matrices.
+Consequently, the control space comprises the directions \(D_x\), \(D_y\), \(D_z\), \(R_x\), and \(R_y\).
+Control performance in each of these directions can be tuned independently.
+A schematic of the proposed control architecture is illustrated in Figure~\ref{fig:test_id31_cf_control}.
+
+\begin{figure}[htbp]
+\centering
+\includegraphics[scale=1]{figs/test_id31_cf_control.png}
+\caption{\label{fig:test_id31_cf_control}Control architecture in the Cartesian frame. Only the controller corresponding to the \(D_z\) direction is shown. \(H_L\) and \(H_H\) are complementary filters.}
+\end{figure}
+
+
+Implementation of this control architecture necessitates a plant model, which must subsequently be inverted.
+This plant model was derived from the multi-body model incorporating the previously detailed 2-DoF APA model, such that the model order stays relatively low.
+Proposed analytical formulas for complementary filters having \(40\,\text{dB/dec}\) were used during this experimental validation.
+An initial experimental validation was conducted under no-payload conditions, with control applied solely to the \(D_y\), \(D_z\), and \(R_y\) directions.
+Increased control bandwidth was achieved for the \(D_z\) and \(R_y\) directions through appropriate tuning of the parameter \(\omega_0\).
+The experimentally measured closed-loop sensitivity transfer functions corresponding to these three controlled directions are presented in Figure \ref{fig:test_id31_cf_control_dy_dz_diff}.
+
+Another test was conducted with a \(26\,\text{kg}\) payload.
+For this configuration, complementary filters were implemented with \(\omega_0 = 2\pi \cdot 10\,\text{rad/s}\), and parameter \(\alpha\) was varied.
+The resulting experimentally obtained closed-loop transfer functions are compared against the theoretical complementary filter responses in Figure~\ref{fig:test_id31_cf_control_alpha}.
+As illustrated in the figure, a close correspondence between the measured closed-loop responses and the target complementary filter magnitude was observed.
+It also shows that the parameter \(\alpha\) provides a mechanism for managing the trade-off between low-frequency disturbance rejection performance and the potential amplification of disturbances within the crossover frequency region.
+
+\begin{figure}[htbp]
+\begin{subfigure}{0.49\textwidth}
+\begin{center}
+\includegraphics[scale=1,scale=0.9]{figs/test_id31_cf_control_dy_dz_diff.png}
+\end{center}
+\subcaption{\label{fig:test_id31_cf_control_dy_dz_diff}Chose of bandwidth using $\omega_0$, $m = 0\,\text{kg}$}
+\end{subfigure}
+\begin{subfigure}{0.49\textwidth}
+\begin{center}
+\includegraphics[scale=1,scale=0.9]{figs/test_id31_cf_control_alpha.png}
+\end{center}
+\subcaption{\label{fig:test_id31_cf_control_alpha}Effect of a change of $\alpha$, $m = 26\,\text{kg}$}
+\end{subfigure}
+\caption{\label{fig:test_id31_cf_control_results}Measured closed-loop transfer functions. Different bandwidth can be specified for different directions using \(\omega_0\) (\subref{fig:test_id31_cf_control_dy_dz_diff}). The shape can be adjusted using parameter \(\alpha\) (\subref{fig:test_id31_cf_control_alpha}).}
+\end{figure}
+
+Finally, \(\omega_0\) was gradually increased to estimate the maximum bandwidth (i.e. the best low frequency disturbance rejection) that can be achieved with this architecture.
+No payload was used for this test, and the parameter \(\omega_0\) was increased for the controllers in the \(D_y\) and \(D_z\) directions.
+A value \(\omega_0 = 2\pi \cdot 60 \,\text{rad/s}\) could be achieved.
+Measured closed-loop transfer functions are shown in Figure \ref{fig:test_id31_high_bandwidth}, indicating a reduction of disturbances in the considered direction of \(1000\) at \(1\,\text{Hz}\).
+For higher values of \(\omega_0\), the system became unstable in the vertical direction, probably because of the resonance at \(250\,\text{Hz}\) that is not well captured with the multi-body model (Figure \ref{fig:test_id31_hac_plant_effect_mass}).
+
+\begin{figure}[htbp]
+\begin{subfigure}{0.49\textwidth}
+\begin{center}
+\includegraphics[scale=1,scale=0.9]{figs/test_id31_high_bandwidth_S.png}
+\end{center}
+\subcaption{\label{fig:test_id31_high_bandwidth_S}Sensitivity}
+\end{subfigure}
+\begin{subfigure}{0.49\textwidth}
+\begin{center}
+\includegraphics[scale=1,scale=0.9]{figs/test_id31_high_bandwidth_T.png}
+\end{center}
+\subcaption{\label{fig:test_id31_high_bandwidth_T}Complementary Sensitivity}
+\end{subfigure}
+\caption{\label{fig:test_id31_high_bandwidth}Measured Closed-Loop Sensitivity (\subref{fig:test_id31_high_bandwidth_S}) and Complementary Sensitivity (\subref{fig:test_id31_high_bandwidth_T}) transfer functions for the highest test bandwidth \(\omega_0 = 2\pi\cdot 60\,\text{rad/s}\).}
+\end{figure}
\section*{Conclusion}
\label{ssec:test_id31_scans_conclusion}
@@ -1104,6 +1171,7 @@ Overall, the experimental results validate the effectiveness of the developed co
The identified limitations, primarily related to high-speed lateral scanning and heavy payload handling, provide clear directions for future improvements.
\begin{table}[htbp]
+\caption{\label{tab:test_id31_experiments_results_summary}Summary of the experimental results performed using the NASS on ID31. Open-loop errors are indicated on the left of the arrows. Closed-loop errors that are outside the specifications are indicated by bold number.}
\centering
\begin{tabularx}{\linewidth}{Xccc}
\toprule
@@ -1132,8 +1200,6 @@ Diffraction tomography (\(6\,\text{deg/s}\), \(1\,mm/s\)) & \(\bm{53}\) & \(10\)
\textbf{Specifications} & \(30\) & \(15\) & \(250\)\\
\bottomrule
\end{tabularx}
-\caption{\label{tab:test_id31_experiments_results_summary}Summary of the experimental results performed using the NASS on ID31. Open-loop errors are indicated on the left of the arrows. Closed-loop errors that are outside the specifications are indicated by bold number.}
-
\end{table}
\chapter*{Conclusion}
\label{ssec:test_id31_conclusion}