Passive Bilateral Teleoperation with Constant Time Delays

Size: px
Start display at page:

Download "Passive Bilateral Teleoperation with Constant Time Delays"

Transcription

1 Proceengs of the 6 IEEE Internatonal Conference on Robotcs an Automaton Orlano, Flora - May 6 Passve Blateral eleoperaton wth Constant me Delays Dongjun Lee an Mark W Spong Coornate Scence Laboratory Unversty of Illnos at Urbana-Champagn 138 W Man St Urbana, IL 6181 USA E-mal: -lee@controlcsluuceu, mspong@uuceu Abstract We propose a novel control framework for blateral teleoperaton of a par of mult-egree-of-freeom (DOF) nonlnear robotc systems uner constant communcaton elays he propose framework utlzes the smple proportonal-ervatve (PD) control, e the master an slave robots are rectly connecte va sprng an amper over the elaye communcaton channels Usng the controller passvty concept, the Lyapunov- Krasovsk technque, an Parseval s entty, we can passfy the combnaton of the elaye communcaton an control blocks altogether robustly, as long as the elays are fnte constants an an upper-boun for the roun-trp elay s known Havng explct poston feeback through the elaye P-acton, the propose framework enforces master-slave poston coornaton whch s often compromse n the velocty-base schemes (eg conventonal scatterng-base teleoperaton) he propose control framework proves humans wth extene physologcal proprocepton so that s/he can affect an sense the remote slave envronments manly relyng on her/hs musculoskeletal systems Experments are performe to valate the propose control framework I INRODUCION Energetcally, as llustrate n Fg 1, a close-loop teleoperator s a two-port system wth the master an slave ports beng couple wth the human operators an slave envronments, respectvely herefore, the foremost an prmary goal of the control (an communcaton) esgn for the teleoperaton shoul be to ensure nteracton safety an/or couple stablty [1] when mechancally couple wth a broa class of slave envronments an humans o ensure such nteracton safety an stablty, energetc passvty (e mechancal power as the supply rate []) of the close-loop teleoperator has been wely utlze as the control objectve [3], [4], [5] In [6], energetc passvty of the elaye teleoperaton s acheve by passfyng the communcaton block wth (possbly unknown) fnte constant tme-elays hs passfcaton was mae possble by applyng scatterng theory In [7], ths scatterng-base result s further extene an the noton of the wave varables was ntrouce Snce these two semnal works, scatterng-base (or wave-base) teleoperaton has been vrtually the only way to enforce energetc passvty of the tme-elaye blateral teleoperaton (eg [3], [8], [9]) In ths paper, we propose a novel control framework for blateral teleoperaton of a par of mult-dof nonlnear robotc systems wth fnte constant communcaton elays he propose framework s base on the smple proportonalervatve (PD) control, e rectly connectng the master an slave robots va sprng an amper over the elaye Human Operator F 1 1 τ F Master 1 Master Slave Comm+ Slave Robot Delay Comm+ Control Robot Control q 1 Fg 1 Sensng q 1 Delay τ Control an Communcaton Close-loop eleoperator q Sensng Close-loop teleoperator as a two-port system q Slave Envron communcaton channels hen, utlzng the controller passvty concept [5], [1], the Lyapunov-Krasovsk technque for elaye systems [11], an Parseval s entty [1], we can enforce energetc passvty of the close-loop teleoperator as long as the communcaton elays are fnte constants an an upper-boun of the roun-trp elay (e sum of the forwar an backwar elays) can be known, even f the elays are asymmetrc (e τ 1 τ n Fg 1) or ther exact estmates are not avalable Relyng on the controller passvty concept, ths passvty s ensure even n the presence of moel parametrc uncertanty (e robust passvty [5], [1]) In [13], we propose a smlar PD-base teleoperaton scheme, whch also enforces passvty However, t requres that the communcaton elays be exactly known an symmetrc (e τ 1 = τ ) In ths paper, we completely remove these two requrements, whch are ffcult to satsfy n practce Compare to the velocty-base schemes (eg conventonal scatterng-base teleoperaton), the man avantage of the propose framework s the explct poston feeback through the elaye P-control acton (e sprng term wth elaye set poston) he lack of such explct poston feeback has been known as the major rawback of the conventonal scatterngbase teleoperaton, n whch, roughly speakng, the velocty nformaton s extracte from the communcate scatterng varables, an then, ntegrate to recover the set poston nformaton herefore, f ths ntegraton becomes naccurate (eg slave robot makes a har contact wth a rg wall or communcaton s blacke out shortly), the master an slave postons may start rftng away from each other (see [14] for example) Havng the explct poston feeback, our propose framework prevents such a poston rftng hs explct poston feeback also enables us to guarantee (e theoretcally prove) asymptotc master-slave poston coornaton In contrast to the scatterng-teleoperaton where the elaye communcaton block s passfe so that the close /6/$ 6 IEEE 9

2 loop teleoperator becomes an nterconnecton of passve submoules (e passfe communcaton, passve control, an passve master/slave robots), the propose framework passfes the combnaton of the communcaton an control blocks altogether Whle the scatterng-base teleoperaton can be use wthout any knowlege on the (fnte constant) elays, our propose framework requres to know an upper-boun of the roun-trp elay However, we thnk that ths s a ml restrcton, snce, n many applcatons, such a roun-trp elay s relatvely easy to measure/estmate [15] Although t acheves at least a level of eal transparency [3], the man goal of the propose control framework s to prove humans wth extene physologcal proprocepton (EPP) [16], e the close-loop teleoperator as a tool, by whch the human operator can affect an sense the remote slave envronments manly relyng on her/hs musculoskeletal systems In ths sense, the propose framework s along the lne of such research works as common passve mechancal tool [4] an vrtual tool for wave-base teleoperaton [9] he rest of ths paper s organze as follows he control problem s formulate n secton II, an the control law s esgne an ts propertes are etale n secton III Expermental results are presente n secton IV an secton V contans some conclung remarks hs paper s a short verson of our recent journal paper [17], to whch we refer reaers for more etals II PROBLEM FORMULAION A Moelng of eleoperators uner Constant me Delay Let us conser a teleoperator consstng of a par of n- egree-of-freeom (DOF) nonlnear robotc systems: M 1 (q 1 ) q 1 (t)+c 1 (q 1, q 1 ) q 1 = 1 (t)+f 1 (t), (1) M (q ) q (t)+c (q, q ) q = (t)+f (t), () where q (t),f (t), (t) R n are the confguratons, human/envronmental force, an controls, an M (q ), C (q, q ) R n n are symmetrc an postve-efnte nerta matrces an Corols matrces, respectvely, st Ṁ (q ) C (q, q ) are skew-symmetrc ( =1, ) Here, we assume that the gravty effects are ether nclue n F 1 (t),f (t) or pre-compensate by the local controls In ths work, we assume the communcaton structure as shown n Fg 1, where the forwar an backwar communcatons are elaye by fnte constant tme-elays τ 1 an τ, respectvely hen, the controls 1 (t), (t) n (1)-() can be efne as functons of the current local nformaton an the elaye remote nformaton (rectly receve from the communcaton lne), e 1 (t) := 1 (q 1 (t), q 1 (t),q (t τ ), q (t τ )) R n, (3) (t) := (q (t), q (t),q 1 (t τ 1 ), q 1 (t τ 1 )) R n (4) B Control Objectves We woul lke to esgn 1 (t), (t) n (3)-(4) to acheve master-slave poston coornaton: f(f 1 (t),f (t)) =, q E (t) :=q 1 (t) q (t), t ; (5) an statc force reflecton: wth ( q 1 (t), q (t), q 1 (t), q 1 (t)), F 1 (t) F (t) (6) For safe nteracton an couple stablty, we woul also lke to enforce the followng energetc passvty of the close-loop teleoperator (1)-(): a fnte constant Rst [ F 1 (θ) q 1 (θ)+f (θ) q (θ) ] θ, t, (7) e maxmum extractable energy from the two-port closeloop teleoperator s boune (see Fg 1) Let us also efne controller passvty [4], [5], [1]: a fnte constant c Rst [ 1 (θ) q 1 (θ)+ (θ) q (θ) ] θ c, t, (8) e the two-port controller n Fg 1 generates only boune amount of energy Lemma 1 [5], [1] For the mechancal teleoperator (1)-(), controller passvty (8) mples energetc passvty (7) Proof: Let us efne the total knetc energy κ f (t) := 1 q 1 (t)m 1 (q 1 ) q 1 (t)+ 1 q (t)m (q ) q (t), (9) then, usng (1)-() wth ts skew-symmetrc property, we have t κ f (t) =[ 1 (t)+f 1 (t)] q 1 (t)+[ (t)+f (t)] q (t) (1) hus, by ntegratng (1) wth the controller passvty conton (8) an the fact that κ f (t), wehave, t, [F 1 (θ) q 1 (θ)+f (θ) q (θ)]θ κ f () c =: Lemma 1 s smple but powerful n the sense that t enables us to analyze energetc passvty (7) of the close-loop teleoperator by examnng only the controller structure whch s often much smpler than that of the close-loop ynamcs Furthermore, by enforcng controller passvty (8), energetc passvty (7) wll be guarantee robustly (robust passvty [5], [1]), because controller passvty (8) oes not epen on the possbly uncertan open-loop ynamcs (1)-() III CONROL DESIGN o acheve the master-slave coornaton (5), blateral force reflecton (6), an energetc passvty (7), we esgn the master an slave controls 1 (t), (t) n (3)-(4) to be 1 (t) := K v ( q 1 (t) q (t τ )) (K + P ɛ ) q 1 (t) K p (q 1 (t) q (t τ )), (11) (t) := K v ( q (t) q 1 (t τ 1 )) (K + P ɛ ) q (t) K p (q (t) q 1 (t τ 1 )), (1) where τ 1,τ are the forwar an backwar fnte constant elays, K v,k p R n n are the symmetrc an postveefnte proportonal (P) an ervatve (D) control gans, K R n n s the sspaton to passfy the elaye P-acton (e wth K p ) n (11)-(1) (to be esgne below), an P ɛ R n n s an atonal ampng ensurng master-slave coornaton 93

3 (5) Inherent evce vscous ampng can substtute ths P ɛ Note that the control laws (11)-(1) contan the explct poston feeback (va the P-acton), whch s generally absent n the conventonal scatterng-base teleoperaton o enforce energetc passvty (7), we esgn sspaton K an P-gan K p n (11)-(1) to satsfy the followng conton: K [ ] sn w(τ 1+τ ) K p K 1 w K p, w R, (13) where, for square matrces A, B, A B mples that A B s postve-semefnte As to be shown below, wth ths conton (13), the elaye P-control acton (e wth K p )s sn wτ passfe by the sspaton K From the fact that τ w, w R, one possble soluton for the conton (13) s K = τ rt K p, (14) where τ rt s an upper-boun of the roun-trp elay τ rt := τ 1 + τ st τ rt τ rt heorem 1 Conser the mechancal teleoperator (1)-() wth the controls (11)-(1) uner the conton (13) 1) (Robust Passvty) he close-loop teleoperator s energetcally passve (e satsfes (7)) regarless of parametrc uncertanty n the open-loop ynamcs (1)-(); ) (Couple Stablty) Suppose that the human operator an slave envronment n Fg 1 are energetcally passve: fnte constants 1, Rst t, F (θ) q (θ) }{{} θ, = {1, }, (15) power nflow to human/envron e the maxmum extractable energy from them are boune hen, q 1 (t), q (t) L hus, f the human an slave envronment are L -stable nput-output mpeance maps, we wll also have F 1 (t),f (t) L ; 3) (Poston Coornaton) Suppose that the human operator an slave envronment are passve n the sense of (15) Suppose further that M jk M (q ), jk (q ) an M jk m (q ) are ml all boune wrt q, where M jk (q ) an q m are the jkth an the m-th components of M (q ) an q, respectvely hen, q E (t) =q 1 (t) q (t) s boune t Moreover, f (F 1 (t),f (t)) = t, (q E (t), q E (t)) ; 4) (Statc Force Reflecton) If ( q 1 (t), q (t), q 1 (t), q (t)), then, F 1 (t) F (t) K p (q 1 q ) Proof: 1) Let us enote the mechancal power generate by the controls (11)-(1) by s c (t) :=1 (t) q 1 (t)+ (t) q (t) = s (t)+s p (t) P (t), (16) where s (t) an s p (t) are the supply rates assocate to the elaye D-acton, an elaye P-acton (+ sspaton K ), respectvely efne by s (t) := q 1 (t)k v q 1 (t)+ q 1 (t)k v q (t τ ) q (t)k v q (t)+ q (t)k v q 1 (t τ 1 ), (17) s p (t) := q 1 (t)k q 1 (t) q 1 (t)k p (q 1 (t) q (t τ )) q (t)k q (t) q (t)k p (q (t) q 1 (t τ 1 )), (18) an P (t) s the followng quaratc form: ( ) [ ) q1 (t) Pɛ ]( q1 (t) P (t) := (19) q (t) P ɛ q (t) We want the total controller supply rate s c (t) n (16) to satsfy the controller passvty (8) Let us frst conser the elaye D-acton supply rate s (t) n (17) hen, usng the fact that, for (, k) ={(1, ), (, 1)}, q (t)k v q k (t τ k ) q (t)k v q (t)+ q k (t τ k )K v q k (t τ k ), () we can show that s (t) =1 1 [ q (t)k v q (t) q (t τ )K v q (t τ ) ] = t V v(t), (1) where V v (t) s a Lyapunov-Krasovsk functonal for elaye systems [11] efne by 1 V v (t) := q (t + θ)k v q (t + θ)θ () =1 τ hen, by ntegratng the nequalty (1), we have s (θ)θ V v (t)+v v (), (3) e energy generaton by the elaye D-acton s always boune by the energy store n V v (t) (e V v (t) s the storage functon for the supply rate s (t)) Now, let us conser the supply rate s p (t) n (18) hen, we can rewrte s p (t) n (18) st: wth N := {(1, ), (, 1)}, s p (t) = q 1 (t)k q 1 (t) q (t)k q (t) q 1 (t)k p (q 1 (t) q (t)) q 1 (t)k p (q (t) q (t τ )) q (t)k p (q (t) q 1 (t)) q (t)k p (q 1 (t) q 1 (t τ 1 )) = q 1 (t)k q 1 (t) q (t)k q (t) t V p(t) q (t)k p (q k (t) q k (t τ k )) (4) (,k) N where V p (t) s the P-acton sprng energy efne as V p (t) := 1 q E(t)K p q E (t), (5) wth q E (t) =q 1 (t) q (t) gven n (5) Let us efne the truncate sgnal q t(θ) of q (t) st { q (θ) t q (θ) f θ [,t] :=, (6) otherwse 94

4 =1, hen, the energy generaton by the supply rate s p (t) n (4) can be wrtten as: s p (θ)θ = V p (t)+v p () [ + θ q 1 t (θ)k p [ + θ q t (θ)k p =1 q t (ξ)ξ θ τ q t (θ)k q (θ)θ t q t (ξ)ξ }{{} q t 1(ξ)ξ :=g t(θ) θ τ1 q t 1(ξ)ξ } {{ } :=g1 t(θ) ] ] θ θ (7) Let us enote the Fourer transform of q t ( = 1, ) by V t(w) := q t (θ)e θ θ = q (θ)e θ θ, where j = 1 hen, the Fourer transform of g t (θ) n (7) (e G t (w) := gt (θ)e θ θ) s gven by G t (w) = 1 e τ V t (w) Let us enote the complex conjugate transpose of a complex vector C n by (e = ) hen usng the Parseval s entty [1], (7) can be rewrtten as s p (θ)θ = V p (t)+v p () (8) =1 (,k) N 1 π 1 π V t (w)k V t (w)w V t = V p (t)+v p () 1 π where N = {(1, ), (, 1)} an [ H(w) = K p K e τ e τ 1 1 e τ k (w)k p Vk t (w)w ( ) V t ( ) 1 (w) V t V t H(w) 1 (w) (w) V t (w) K p e τ 1 e τ K w, ] C n n For more etals on these ervatons, see [17] Snce H(w) s Hermtan (e H (w) = H (w) =H(w)) wth postveefnte block agonal matrces K, followng [18, pp473], H(w) (e H(w) s postve-semefnte) f an only f K eτ1 e τ e τ e τ1 K p K 1 K p, (9) whch s nothng but the conton (13) hus, wth the conton (13) st H(w), we can show from (8) that s p (θ)θ V p (t)+v p (), (3) e energy generaton by the supply rate s p (t) s always boune by the P-acton sprng energy V p (t) n (5) hus, by summng up (3) an (3) wth (16) an the fact that V v (t) an V p (t) t, we can prove controller passvty (8) st for all t, [ 1 (θ) q 1 (θ)+ (θ) q (θ)]θ V v (t)+v v () V p (t)+v p () P (θ)θ V v () + V p () =: c, (31) where the term V v () wll be zero f we start from zero veloctes (e ( q 1 (t), q (t)) = t (, ]), whle the term V p () woul be small, f the ntal coornaton error q E () = q 1 () q () s small Fnally, from Lemma 1, energetc passvty (7) of the close-loop teleoperator follows ) By ntegratng (1) wth (31) an (15), we have, t, κ f (t)+v v (t)+v p (t) (3) κ f () + V v () + V p () P (θ)θ + 1 +, where P (t) n (19) Here, snce V v (),V p (), 1, an κ f () are all boune, κ f (t) s boune hus, q 1 (t), q (t) are also boune t (e q 1 (t), q (t) L ) herefore, f the human operator an the slave envronment are L -stable mpeance maps, F 1 (t),f (t) L 3) Bouneness of q E (t) =q 1 (t) q (t) s a rect consequence of (3) wth the efnton of V p (t) n (5) Frst step of the poston-coornaton proof s to show that ( q 1 (t), q (t)) Suppose that F 1 (t),f (t) =, t hen, from (3) wth 1 = =an the bouneness of P ɛ,m 1 (q 1 ),M (q ),wehave: t, κ f (t) κ f () + c κ f () + c γ P (θ)θ κ f (θ)θ, (33) where c = V v () + V p () as gven n (31) an γ>s a constant scalar Here, snce κ f (t), the term κ f (θ)θ s monotoncally ncreasng an upper boune, thus, t converges to a lmt herefore, followng Barbalat s lemma, f κ f (t) s unformly contnuous, κ f (t) wll also converge to (e ( q 1 (t), q (t)) ) o show ths, let us conser t κ f (t) hen, from (1) wth F 1 (t) =F (t) =,wehave t κ f (t) = 1 (t) q 1 (t)+ (t) q (t), where q 1 (t), q (t) are boune from tem of ths theorem Also, 1 (t), (t) are boune, snce, n ther efntons (11)-(1), for (,k) ={(1, ), (, 1)}, 1) q (t) q k (t τ k ) s boune (wth boune q (t)); an ) q (t) q k (t τ k )=q E (t)+ τ k q k (t + θ)θ s also boune (wth boune q (t), q E (t) an τ k ) hus, t κ f (t) s boune, an κ f (t) s unformly contnuous herefore, κ f (t) an ( q 1 (t), q (t)) he next step of the proof s to show that ( q 1 (t), q (t)) Let us conser the ynamcs (1)-() wth F 1 (t) =F (t) =, where, as shown n the above paragraph, the controls 1 (t), (t) n (11)-(1) are boune Also, from the bouneness assumpton of Mjk (q ), the Corols terms C m (q, q ) q ( = 1, ) n (1)-() are boune hus, the acceleratons 95

5 q 1 (t), q (t) are also boune t Now, let us conser the acceleraton q (t) n (1)-() (wth F (t) =): q = M 1 (q )C (q, q ) q + M 1 (q ) (t), =1,, (34) where the tme-ervatves of the terms n the RHS are all boune, ue to the bouneness of q (t), q (t), q E (t), t M 1 (q )= M 1 (q ) t M (q )M 1 (q ) (from the bouneness of Mjk (q ) ), an m t [C (q, q ) q ] (from the bouneness of M jk (q ), q, an q ) hs mples that the RHS ml of (34) s unformly contnuous hus, q 1 (t), q (t) are also unformly contnuous herefore, followng Barbalat s lemma, ( q 1 (t), q (t)) as ( q 1 (t), q (t)) Now, let us conser the ynamcs (1)-() wth the controls 1 (t), (t) n (11)-(1) an F 1 (t) = F (t) = hen, snce ( q 1 (t), q (t), q 1 (t), q (t)), we have K p (q (t) q k (t τ k )), (, k) = {(1, ), (, 1)} hs conton can be rewrtten as K p (q E (t) + τ k q k (t + θ)θ), where the secon term n the parenthess goes to zero, because q k (t) an τ k s fnte herefore, snce K p s postveefnte, q E (t) (e q 1 (t) q (t)) 4) Suppose that ( q 1 (t), q (t), q 1 (t), q (t)) hen, from the ynamcs (1)-() an ther controls (11)-(1), we have: F 1 (t) K p (q 1 q ), F (t) K p (q q 1 ), (35) wth q (t τ) q (t) q as q (t), q (t) In heorem 1, the negatve sgns n the passvty conton for the human an slave envronment (15) come from the fact that the power nflows to those systems are gven by F (t) q (t), e the prouct of the reacton force F (t) an the nteracton velocty q (t) he bouneness assumpton on M jk (q ), Mjk (q ) an M jk m (q ) n heorem 1 s guarantee, ml f the robot s confguraton space s compact an ts nerta matrx s smooth Such compact confguraton space an smooth nerta are possesse by many practcal robotc systems (eg revolute jont robots) he conton (14) (or (13)) enables us to passfy the elaye P-acton n (11)-(1) by the sspaton K hs elaye P-acton contans an explct poston feeback nformaton, the lack of whch s recognze as the man cause of the master-slave poston rft n the conventonal scatterng-base teleoperaton In contrast, the elaye D-acton n (11)-(1) s tself passve wth the Lyapunov-Krasovsk functon V v (t) n () as ts storage functon (see (3)) Snce the conton (14) can be acheve as long as the elays τ 1,τ are fnte constants an ther roun-trp elay (e τ 1 + τ ) s upperboune, energetc passvty (7) can also be ensure wth such fnte constant elays, even f they are not exactly known or asymmetrc (e τ 1 τ ) In [13], we propose a smlar PD-base scheme, whch also enforces passvty However, t requres the elays to be symmetrc an exactly known he conton (14) mposes the followng mplcatons on the system performance: 1) wth the same statc force reflecton performance (e same K p, see tem 4 of heorem 1), the moton aglty (e less K ) woul be compromse as the elays becomes longer, snce, n the conton (14), the requre sspaton K s proportonal to the roun-trp elay τ rt ; an ) wth the elays fxe, there s a trae-off between the statc force reflecton performance (e K p ) an moton aglty (e K ), snce, uner the conton (14), a large K p (e sharp force reflecton) requres a large K (e poor moton aglty), or a small K (e agle free-moton) permts only a small K p (e poor force reflecton) he key step n the proof of heorem 1 s the use of the Parseval s entty n (8) whch we assume to be true A suffcent conton for the Parseval s entty to hol s that q 1 (t), q (t) L [1] As the followng Lemma shows, ths suffcent conton s guarantee n many practcal stuatons where the human an slave envronment are passve n the sense of (15), the nerta matrces of the master an slave robots (1)-() have boune partal-ervatve wrt q (e M jk (q ) are boune), an the master an slave veloctes m an the coornaton error are ntally boune Due to the space lmtaton, we omt the proof here an refer reaers to [17] he Parseval s entty was also use n [19] to ensure the energetc passvty of haptc-nterfaces uner zero-orer-hol Lemma Suppose that the human an slave envronment are passve n the sense of (15) an efne L -stable mpeance maps (e f q 1 (t), q (t) L, F 1 (t),f (t) L ) Suppose further that Mjk (q ) are boune hen, f q m 1 (), q (),q E () are boune, q 1 (t), q (t) L, thus, the Parseval s entty n (8) hols an the tems 1-3 of heorem 1 are ensure IV EXPERIMEN For the experment, we use a par of rect-rve planar -DOF seral-lnks revolute-jont robots More etals on the harware/software constructon of ths system can be foun n [17] along wth more experment an smulaton results he control laws (11)-(1) are erve for the Cartesanspace ynamcs so that, followng tem 4 of heorem 1, the Cartesan (statc) force reflecton can be acheve We nstall an alumnum wall n the slave envronment, an 5[ms] samplng rate s obtane We set the elays st (τ 1,τ ) = (1, 18)[sec] (e τ rt = 3[sec]) he P-gan K p an sspaton K n (11)-(1) are esgne accorng to the conton (14) Atonal ampng P ɛ s omtte n the control mplementaton, as we leave the evce frcton uncompensate We conser the followng scenaro: 1) ntally, a human stablzes the teleoperator; ) then, wthout seeng, moves the slave close to the alumnum wall an keeps pushng the master untl s/he perceves the wall; 3) makes a har contact; an 4) fnally, retracts the slave from the wall Expermental results are gven n Fg As shown by the force profle n Fg, the human operator can perceve the alumnum wall through the force reflecton Also, when the contact s remove (e free-moton), the master an slave postons become coornate wth each other hese freemoton an contact behavor are all stable, as we enforce passvty (7) through the conton (14) After the contact, 96

6 the human operates the master so that the slave returns to ts startng poston hs leas nto bumps n the human force n Fg after aroun 6[sec] In Fg, there are some errors n both the force reflecton an poston coornaton (eg aroun 8[sec]) hese errors are ue to the substantal (bearng) Coulomb frcton of the robots, whch we foun can go up to 1[N] For nstance, suppose n (1)-() that F 1 = F h μ h an F = F e μ e, where F h,f e are the human/contact force, an μ h,μ e are the frctons Durng the contact, the frctons always oppose the human/contact forces hen, from tem 4 of heorem 1, F h F e (μ e + μ h ), e (μ e + μ h ) causes force reflecton error (eg aroun 4[sec]) hese errors were not observe when we performe a smulaton wthout such frctons [17] Although we not observe any substantal problems urng the experment, the requre sspaton K (from (14)) may cause sluggsh system behavor an lmt the framework s usablty, especally when the task requres agle operaton but the elays are large How to further mnmze ths requre K whle enforcng passvty s a topc for future work V CONCLUSIONS he man avantage of the propose framework s n that we can explot the benefts of explct poston feeback whle stll retanng passvty Also, ue to ts smple PD-base structure whch are pervasve n many nustral control systems, we thnk that the propose framework woul make a goo mpact n the real practce of teleoperaton applcatons We beleve that the propose framework s promsng for the Internet teleoperaton, where ts explct poston feeback woul o an mportant role to recover the poston coornaton n the presence of packet-loss an tme-varyng elays ACKNOWLEDGMEN Research partally supporte by the Offce of Naval Research uner grants N an N , the Natonal Scence Founaton uner grants IIS , CCR -9 an ECS-1-41, an the College of Engneerng at the Unversty of Illnos REFERENCES [1] J E Colgate Couple stablty of multport systems - theory an experments ransactons of the ASME, Journal of Dynamc Systems, Measurement an Control, 116(3):419 48, 1994 [] J C Wllems Dsspatve ynamcal systems part1: general theory Arch Ratonal Mech Anal, 45():31 351, 197 [3] D A Lawrence Stablty an transparency n blateral teleoperaton IEEE ransactons on Robotcs an Automaton, 9(5):64 637, 1993 [4] D J Lee an P Y L Passve blateral feeforwar control of lnear ynamcally smlar teleoperate manpulators IEEE ransactons on Robotcs an Automaton, 19(3): , 3 [5] D J Lee an P Y L Passve blateral control an tool ynamcs renerng for nonlnear mechancal teleoperators IEEE ransactons on Robotcs, 1(5): , 5 [6] R J Anerson an M W Spong Blateral control of tele-operators wth tme elay IEEE ransactons on Automatc Control, 34(5):494 51, 1989 [7] G Nemeyer an J J E Slotne Stable aaptve teleoperaton IEEE Journal of Oceanc Engneerng, 16(1):15 16, 1991 [8] N Chopra, M W Spong, S Hrche, an M Buss Blateral teleoperaton over the nternet: the tme varyng elay problem In Proceengs of Amercan Control Conference, pages , 3 Fg x Force [N] y Force [N] q1 q [m] Master an Slave x Forcng: τ rt =3sec, Half Kp Master an Slave y Forcng: τ rt =3sec, Half Kp F1x Fx 1 F1y Fy tme [sec] Master Slave Poston Error: τ rt =3sec, Half Kp tme [sec] Expermental results wth the roun-trp elay τ rt =3[sec] [9] G Nemeyer an J J E Slotne elemanpulaton wth tme elays Internatonal Journal of Robotcs Research, 3(9):873 89, 4 [1] Dongjun Lee Passve Decomposton an Control of Interactve Mechancal Systems uner Coornaton Requrements Doctoral Dssertaton, Unversty of Mnnesota, 4 [11] K Gu an S Nculescu Survey on recent results n the stablty an control of tme-elay systems ASME Journal of Dynamc Systems, Measurements, an Control, 15: , 3 [1] R R Golberg Fourer ransforms he Cambrge Unversty Press, New York, NY, 1961 [13] D J Lee an M W Spong Passve blateral control of teleoperators uner constant tme-elay In Proceengs of the IFAC Worl Congress, 5 [14] S Hrche an M Buss Packet loss effects n passve telepresence systems In Proceengs of the IEEE Conference on Decson an Control, pages , 4 [15] O Gurewtz an M S Estmatng one-way elays from cyclc-path elay measurements In Proceengs of IEEE INFOCOM, pages , 1 [16] D S Chlress Control strategy for upper-lmb prostheses Proceengs of the th Annual Internatonal Conf of the IEEE Engneerng n Mecne an Bology Socety, (5):73 75, 1998 [17] D J Lee an M W Spong Passve blateral teleoperaton wth constant tme-elay IEEE ransactons on Robotcs, 6 o appear Preprnt avalable at lee/leespongdelayro4pf [18] R A Horn an C R Johnson Matrx analyss Cambrge Unversty Press, Cambrge, UK, 1985 [19] J E Colgate an G Schenkel Passvty of a class of sample-ata systems: applcaton to haptc nterfaces Journal of Robotc Systems, 14(1):37 47,

ENGI9496 Lecture Notes Multiport Models in Mechanics

ENGI9496 Lecture Notes Multiport Models in Mechanics ENGI9496 Moellng an Smulaton of Dynamc Systems Mechancs an Mechansms ENGI9496 Lecture Notes Multport Moels n Mechancs (New text Secton 4..3; Secton 9.1 generalzes to 3D moton) Defntons Generalze coornates

More information

Analytical classical dynamics

Analytical classical dynamics Analytcal classcal ynamcs by Youun Hu Insttute of plasma physcs, Chnese Acaemy of Scences Emal: yhu@pp.cas.cn Abstract These notes were ntally wrtten when I rea tzpatrck s book[] an were later revse to

More information

High-Order Hamilton s Principle and the Hamilton s Principle of High-Order Lagrangian Function

High-Order Hamilton s Principle and the Hamilton s Principle of High-Order Lagrangian Function Commun. Theor. Phys. Bejng, Chna 49 008 pp. 97 30 c Chnese Physcal Socety Vol. 49, No., February 15, 008 Hgh-Orer Hamlton s Prncple an the Hamlton s Prncple of Hgh-Orer Lagrangan Functon ZHAO Hong-Xa an

More information

New Liu Estimators for the Poisson Regression Model: Method and Application

New Liu Estimators for the Poisson Regression Model: Method and Application New Lu Estmators for the Posson Regresson Moel: Metho an Applcaton By Krstofer Månsson B. M. Golam Kbra, Pär Sölaner an Ghaz Shukur,3 Department of Economcs, Fnance an Statstcs, Jönköpng Unversty Jönköpng,

More information

Globally Stable Output Feedback Synchronization of Teleoperation with Time-Varying Delays

Globally Stable Output Feedback Synchronization of Teleoperation with Time-Varying Delays Globally Stable Output Feedback Synchronzaton of Teleoperaton wth Tme-Varyng Delays Yuan Yang Department of Mechancal Engneerng Unversty of Vctora Vctora, Brtsh Columba 25-853-322 Emal: yangyuan@uvc.ca

More information

Large-Scale Data-Dependent Kernel Approximation Appendix

Large-Scale Data-Dependent Kernel Approximation Appendix Large-Scale Data-Depenent Kernel Approxmaton Appenx Ths appenx presents the atonal etal an proofs assocate wth the man paper [1]. 1 Introucton Let k : R p R p R be a postve efnte translaton nvarant functon

More information

Hybrid Virtual-Proxy Based Control Framework for Passive Bilateral Teleoperation over the Internet

Hybrid Virtual-Proxy Based Control Framework for Passive Bilateral Teleoperation over the Internet Hybrd Vrtual-Proxy Based Control Framewor for Passve Blateral Teleoperaton over the Internet Ke Huang and Dongjun Lee Abstract Let us consder the control problem of blateral teleoperaton over the Internet.

More information

WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION?

WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION? ISAHP 001, Berne, Swtzerlan, August -4, 001 WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION? Masaak SHINOHARA, Chkako MIYAKE an Kekch Ohsawa Department of Mathematcal Informaton Engneerng College

More information

An Approach to Formation Maneuvers of Multiple Nonholonomic Agents using Passivity Techniques

An Approach to Formation Maneuvers of Multiple Nonholonomic Agents using Passivity Techniques An Approach to Formaton Maneuvers of Multple onholonomc Agents usng Passvty echnques Wu Fan, Geng Zhyong he State Key Laboratory for urbulence an Complex Systems, Department of Mechancs an Aerospace Engneerng,

More information

Chapter 2 Transformations and Expectations. , and define f

Chapter 2 Transformations and Expectations. , and define f Revew for the prevous lecture Defnton: support set of a ranom varable, the monotone functon; Theorem: How to obtan a cf, pf (or pmf) of functons of a ranom varable; Eamples: several eamples Chapter Transformatons

More information

Digital Signal Processing

Digital Signal Processing Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over

More information

Yukawa Potential and the Propagator Term

Yukawa Potential and the Propagator Term PHY304 Partcle Physcs 4 Dr C N Booth Yukawa Potental an the Propagator Term Conser the electrostatc potental about a charge pont partcle Ths s gven by φ = 0, e whch has the soluton φ = Ths escrbes the

More information

Competitive Experimentation and Private Information

Competitive Experimentation and Private Information Compettve Expermentaton an Prvate Informaton Guseppe Moscarn an Francesco Squntan Omtte Analyss not Submtte for Publcaton Dervatons for te Gamma-Exponental Moel Dervaton of expecte azar rates. By Bayes

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Neuro-Adaptive Design - I:

Neuro-Adaptive Design - I: Lecture 36 Neuro-Adaptve Desgn - I: A Robustfyng ool for Dynamc Inverson Desgn Dr. Radhakant Padh Asst. Professor Dept. of Aerospace Engneerng Indan Insttute of Scence - Bangalore Motvaton Perfect system

More information

Bounds for Spectral Radius of Various Matrices Associated With Graphs

Bounds for Spectral Radius of Various Matrices Associated With Graphs 45 5 Vol.45, No.5 016 9 AVANCES IN MATHEMATICS (CHINA) Sep., 016 o: 10.11845/sxjz.015015b Bouns for Spectral Raus of Varous Matrces Assocate Wth Graphs CUI Shuyu 1, TIAN Guxan, (1. Xngzh College, Zhejang

More information

Improved delay-dependent stability criteria for discrete-time stochastic neural networks with time-varying delays

Improved delay-dependent stability criteria for discrete-time stochastic neural networks with time-varying delays Avalable onlne at www.scencedrect.com Proceda Engneerng 5 ( 4456 446 Improved delay-dependent stablty crtera for dscrete-tme stochastc neural networs wth tme-varyng delays Meng-zhuo Luo a Shou-mng Zhong

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

THE GUARANTEED COST CONTROL FOR UNCERTAIN LARGE SCALE INTERCONNECTED SYSTEMS

THE GUARANTEED COST CONTROL FOR UNCERTAIN LARGE SCALE INTERCONNECTED SYSTEMS Copyrght 22 IFAC 5th rennal World Congress, Barcelona, Span HE GUARANEED COS CONROL FOR UNCERAIN LARGE SCALE INERCONNECED SYSEMS Hroak Mukadan Yasuyuk akato Yoshyuk anaka Koch Mzukam Faculty of Informaton

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Chapter 7: Conservation of Energy

Chapter 7: Conservation of Energy Lecture 7: Conservaton o nergy Chapter 7: Conservaton o nergy Introucton I the quantty o a subject oes not change wth tme, t means that the quantty s conserve. The quantty o that subject remans constant

More information

Dynamic Modeling of a Synchronous Generator Using T-S Fuzzy Approach

Dynamic Modeling of a Synchronous Generator Using T-S Fuzzy Approach e-issn : 0975-0 Hee-Jn Lee / Internatonal Journal of Engneerng an echnology (IJE) Dynamc oelng of a Synchronous Generator Usng -S Fuzzy Approach Hee-Jn Lee Department of Electronc Engneerng Kumoh Natonal

More information

Adaptive Consensus Control of Multi-Agent Systems with Large Uncertainty and Time Delays *

Adaptive Consensus Control of Multi-Agent Systems with Large Uncertainty and Time Delays * Journal of Robotcs, etworkng and Artfcal Lfe, Vol., o. (September 04), 5-9 Adaptve Consensus Control of Mult-Agent Systems wth Large Uncertanty and me Delays * L Lu School of Mechancal Engneerng Unversty

More information

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,

More information

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force. The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

ECE559VV Project Report

ECE559VV Project Report ECE559VV Project Report (Supplementary Notes Loc Xuan Bu I. MAX SUM-RATE SCHEDULING: THE UPLINK CASE We have seen (n the presentaton that, for downlnk (broadcast channels, the strategy maxmzng the sum-rate

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

IMECE PASSIVE TELEOPERATION OF A MULTI DEGREE OF FREEDOM HYDRAULIC BACKHOE USING DYNAMIC PASSIVE VALVE

IMECE PASSIVE TELEOPERATION OF A MULTI DEGREE OF FREEDOM HYDRAULIC BACKHOE USING DYNAMIC PASSIVE VALVE Proceedngs of IMECE 03 003 ASME Internatonal Mechancal Engneerng Congress Washngton, D.C., November 15 1, 003 IMECE003-43719 PASSIVE TELEOPERATION OF A MULTI DEGREE OF FREEDOM HYDRAULIC BACKHOE USING DYNAMIC

More information

Kinematics of Fluid Motion

Kinematics of Fluid Motion Knematcs of Flu Moton R. Shankar Subramanan Department of Chemcal an Bomolecular Engneerng Clarkson Unversty Knematcs s the stuy of moton wthout ealng wth the forces that affect moton. The scusson here

More information

Distributed Exponential Formation Control of Multiple Wheeled Mobile Robots

Distributed Exponential Formation Control of Multiple Wheeled Mobile Robots Proceedngs of the Internatonal Conference of Control, Dynamc Systems, and Robotcs Ottawa, Ontaro, Canada, May 15-16 214 Paper No. 46 Dstrbuted Exponental Formaton Control of Multple Wheeled Moble Robots

More information

On a one-parameter family of Riordan arrays and the weight distribution of MDS codes

On a one-parameter family of Riordan arrays and the weight distribution of MDS codes On a one-parameter famly of Roran arrays an the weght strbuton of MDS coes Paul Barry School of Scence Waterfor Insttute of Technology Irelan pbarry@wte Patrck Ftzpatrck Department of Mathematcs Unversty

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 NMT EE 589 & UNM ME 48/58 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 48/58 7. Robot Dynamcs 7.5 The Equatons of Moton Gven that we wsh to fnd the path q(t (n jont space) whch mnmzes the energy

More information

Visualization of 2D Data By Rational Quadratic Functions

Visualization of 2D Data By Rational Quadratic Functions 7659 Englan UK Journal of Informaton an Computng cence Vol. No. 007 pp. 7-6 Vsualzaton of D Data By Ratonal Quaratc Functons Malk Zawwar Hussan + Nausheen Ayub Msbah Irsha Department of Mathematcs Unversty

More information

High resolution entropy stable scheme for shallow water equations

High resolution entropy stable scheme for shallow water equations Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal

More information

Non-negative Matrices and Distributed Control

Non-negative Matrices and Distributed Control Non-negatve Matrces an Dstrbute Control Yln Mo July 2, 2015 We moel a network compose of m agents as a graph G = {V, E}. V = {1, 2,..., m} s the set of vertces representng the agents. E V V s the set of

More information

Iterative General Dynamic Model for Serial-Link Manipulators

Iterative General Dynamic Model for Serial-Link Manipulators EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 1. Introducton Iteratve General Dynamc Model for Seral-Lnk Manpulators In ths set of notes, we are gong to develop a method for computng a general

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

ENTROPIC QUESTIONING

ENTROPIC QUESTIONING ENTROPIC QUESTIONING NACHUM. Introucton Goal. Pck the queston that contrbutes most to fnng a sutable prouct. Iea. Use an nformaton-theoretc measure. Bascs. Entropy (a non-negatve real number) measures

More information

Output Feedback Stabilization of Networked Control Systems with Random Delays

Output Feedback Stabilization of Networked Control Systems with Random Delays Preprnts of the 8th IFAC Worl Congress Mlano (Italy) August 28 - September 2 2 Output Feebac Stablzaton of Networe Control Systems wth Ranom Delays Shou-Wan Gao Gong-You ang College of Informaton Scence

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

Advanced Mechanical Elements

Advanced Mechanical Elements May 3, 08 Advanced Mechancal Elements (Lecture 7) Knematc analyss and moton control of underactuated mechansms wth elastc elements - Moton control of underactuated mechansms constraned by elastc elements

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Explicit bounds for the return probability of simple random walk

Explicit bounds for the return probability of simple random walk Explct bouns for the return probablty of smple ranom walk The runnng hea shoul be the same as the ttle.) Karen Ball Jacob Sterbenz Contact nformaton: Karen Ball IMA Unversty of Mnnesota 4 Ln Hall, 7 Church

More information

Journal of System Design and Dynamics

Journal of System Design and Dynamics Vbraton Control of Safe Robot Arm wth MR-Based Passve Complant Jont * Seung-kook YUN**,, Seong-Sk YOON**,, Sungchul KANG**, and Munsang KIM** **Intellgent Robotcs Research Center, Korea Insttute of Scence

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

Physics 3A: Linear Momentum. Physics 3A: Linear Momentum. Physics 3A: Linear Momentum. Physics 3A: Linear Momentum

Physics 3A: Linear Momentum. Physics 3A: Linear Momentum. Physics 3A: Linear Momentum. Physics 3A: Linear Momentum Recall that there was ore to oton than just spee A ore coplete escrpton of oton s the concept of lnear oentu: p v (8.) Beng a prouct of a scalar () an a vector (v), oentu s a vector: p v p y v y p z v

More information

Lecture 14: Bandits with Budget Constraints

Lecture 14: Bandits with Budget Constraints IEOR 8100-001: Learnng and Optmzaton for Sequental Decson Makng 03/07/16 Lecture 14: andts wth udget Constrants Instructor: Shpra Agrawal Scrbed by: Zhpeng Lu 1 Problem defnton In the regular Mult-armed

More information

Note 10. Modeling and Simulation of Dynamic Systems

Note 10. Modeling and Simulation of Dynamic Systems Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada

More information

Adaptive sliding mode reliable excitation control design for power systems

Adaptive sliding mode reliable excitation control design for power systems Acta Technca 6, No. 3B/17, 593 6 c 17 Insttute of Thermomechancs CAS, v.v.. Adaptve sldng mode relable exctaton control desgn for power systems Xuetng Lu 1, 3, Yanchao Yan Abstract. In ths paper, the problem

More information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. An adaptive method for inertia force identification in cantilever under moving mass 89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Adaptive Synchronization of Rössler System Based on. Dynamic Surface ControlPF

Adaptive Synchronization of Rössler System Based on. Dynamic Surface ControlPF FP Mar. 6 Journal of Electronc Scence an Technology of Chna Vol. No. Aaptve Synchronzaton of Rössler System Base on * Dynamc Surface ControlPF LIU Tong-shuan, XU Hao, GUAN Xn-png Insttute of Electrcal

More information

Robust Dynamic Programming for Discounted Infinite-Horizon Markov Decision Processes with Uncertain Stationary Transition Matrice

Robust Dynamic Programming for Discounted Infinite-Horizon Markov Decision Processes with Uncertain Stationary Transition Matrice roceengs of the 2007 IEEE Symposum on Approxmate Dynamc rogrammng an Renforcement Learnng (ADRL 2007) Robust Dynamc rogrammng for Dscounte Infnte-Horzon Markov Decson rocesses wth Uncertan Statonary Transton

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

Design and Analysis of Landing Gear Mechanic Structure for the Mine Rescue Carrier Robot

Design and Analysis of Landing Gear Mechanic Structure for the Mine Rescue Carrier Robot Sensors & Transducers 214 by IFSA Publshng, S. L. http://www.sensorsportal.com Desgn and Analyss of Landng Gear Mechanc Structure for the Mne Rescue Carrer Robot We Juan, Wu Ja-Long X an Unversty of Scence

More information

Some modelling aspects for the Matlab implementation of MMA

Some modelling aspects for the Matlab implementation of MMA Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton

More information

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force. Secton 1. Dynamcs (Newton s Laws of Moton) Two approaches: 1) Gven all the forces actng on a body, predct the subsequent (changes n) moton. 2) Gven the (changes n) moton of a body, nfer what forces act

More information

( ) = : a torque vector composed of shoulder torque and elbow torque, corresponding to

( ) = : a torque vector composed of shoulder torque and elbow torque, corresponding to Supplementary Materal for Hwan EJ, Donchn O, Smth MA, Shamehr R (3 A Gan-Fel Encon of Lmb Poston an Velocty n the Internal Moel of Arm Dynamcs. PLOS Boloy, :9-. Learnn of ynamcs usn bass elements he nternal

More information

Temperature. Chapter Heat Engine

Temperature. Chapter Heat Engine Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the

More information

On Liu Estimators for the Logit Regression Model

On Liu Estimators for the Logit Regression Model CESIS Electronc Workng Paper Seres Paper No. 59 On Lu Estmators for the Logt Regresson Moel Krstofer Månsson B. M. Golam Kbra October 011 The Royal Insttute of technology Centre of Excellence for Scence

More information

A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON

A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON PIOTR NAYAR AND TOMASZ TKOCZ Abstract We prove a menson-free tal comparson between the Euclean norms of sums of nepenent ranom vectors

More information

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k) ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of

More information

Tracking with Kalman Filter

Tracking with Kalman Filter Trackng wth Kalman Flter Scott T. Acton Vrgna Image and Vdeo Analyss (VIVA), Charles L. Brown Department of Electrcal and Computer Engneerng Department of Bomedcal Engneerng Unversty of Vrgna, Charlottesvlle,

More information

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Perron Vectors of an Irreducible Nonnegative Interval Matrix Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecture 3 Contnuous Systems an Fels (Chapter 3) Where Are We Now? We ve fnshe all the essentals Fnal wll cover Lectures through Last two lectures: Classcal Fel Theory Start wth wave equatons

More information

Physics 141. Lecture 14. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 14, Page 1

Physics 141. Lecture 14. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 14, Page 1 Physcs 141. Lecture 14. Frank L. H. Wolfs Department of Physcs and Astronomy, Unversty of Rochester, Lecture 14, Page 1 Physcs 141. Lecture 14. Course Informaton: Lab report # 3. Exam # 2. Mult-Partcle

More information

Complex Numbers, Signals, and Circuits

Complex Numbers, Signals, and Circuits Complex Numbers, Sgnals, and Crcuts 3 August, 009 Complex Numbers: a Revew Suppose we have a complex number z = x jy. To convert to polar form, we need to know the magntude of z and the phase of z. z =

More information

Composite Nonlinear Feedback-based Bounded Formation Control of Multi-quadrotor Systems

Composite Nonlinear Feedback-based Bounded Formation Control of Multi-quadrotor Systems Composte Nonlnear Feeback-base Boune Formaton Control of Mult-quarotor Systems Zhcheng Hou, Isabelle Fanton To cte ths verson: Zhcheng Hou, Isabelle Fanton. Composte Nonlnear Feeback-base Boune Formaton

More information

GENERIC CONTINUOUS SPECTRUM FOR MULTI-DIMENSIONAL QUASIPERIODIC SCHRÖDINGER OPERATORS WITH ROUGH POTENTIALS

GENERIC CONTINUOUS SPECTRUM FOR MULTI-DIMENSIONAL QUASIPERIODIC SCHRÖDINGER OPERATORS WITH ROUGH POTENTIALS GENERIC CONTINUOUS SPECTRUM FOR MULTI-DIMENSIONAL QUASIPERIODIC SCHRÖDINGER OPERATORS WITH ROUGH POTENTIALS YANG FAN AND RUI HAN Abstract. We stuy the mult-mensonal operator (H xu) n = m n = um + f(t n

More information

The Noether theorem. Elisabet Edvardsson. Analytical mechanics - FYGB08 January, 2016

The Noether theorem. Elisabet Edvardsson. Analytical mechanics - FYGB08 January, 2016 The Noether theorem Elsabet Evarsson Analytcal mechancs - FYGB08 January, 2016 1 1 Introucton The Noether theorem concerns the connecton between a certan kn of symmetres an conservaton laws n physcs. It

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON

A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON PIOTR NAYAR AND TOMASZ TKOCZ Abstract We prove a menson-free tal comparson between the Euclean norms of sums of nepenent ranom vectors

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

Bounds for the Laplacian spectral radius of graphs Huiqing Liu a ; Mei Lu b a

Bounds for the Laplacian spectral radius of graphs Huiqing Liu a ; Mei Lu b a Ths artcle was ownloae by: [Tsnghua Unversty] On: 16 December 2009 Access etals: Access Detals: [subscrpton number 912295224] Publsher Taylor & Francs Informa Lt Regstere n Englan an Wales Regstere Number:

More information

Chapter Newton s Method

Chapter Newton s Method Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve

More information

Parameter Estimation for Dynamic System using Unscented Kalman filter

Parameter Estimation for Dynamic System using Unscented Kalman filter Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,

More information

Tornado and Luby Transform Codes. Ashish Khisti Presentation October 22, 2003

Tornado and Luby Transform Codes. Ashish Khisti Presentation October 22, 2003 Tornado and Luby Transform Codes Ashsh Khst 6.454 Presentaton October 22, 2003 Background: Erasure Channel Elas[956] studed the Erasure Channel β x x β β x 2 m x 2 k? Capacty of Noseless Erasure Channel

More information

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1]

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1] DYNAMIC SHORTEST PATH SEARCH AND SYNCHRONIZED TASK SWITCHING Jay Wagenpfel, Adran Trachte 2 Outlne Shortest Communcaton Path Searchng Bellmann Ford algorthm Algorthm for dynamc case Modfcatons to our algorthm

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

Control of Uncertain Bilinear Systems using Linear Controllers: Stability Region Estimation and Controller Design

Control of Uncertain Bilinear Systems using Linear Controllers: Stability Region Estimation and Controller Design Control of Uncertan Blnear Systems usng Lnear Controllers: Stablty Regon Estmaton Controller Desgn Shoudong Huang Department of Engneerng Australan Natonal Unversty Canberra, ACT 2, Australa shoudong.huang@anu.edu.au

More information

PHZ 6607 Lecture Notes

PHZ 6607 Lecture Notes NOTE PHZ 6607 Lecture Notes 1. Lecture 2 1.1. Defntons Books: ( Tensor Analyss on Manfols ( The mathematcal theory of black holes ( Carroll (v Schutz Vector: ( In an N-Dmensonal space, a vector s efne

More information

12. The Hamilton-Jacobi Equation Michael Fowler

12. The Hamilton-Jacobi Equation Michael Fowler 1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and

More information

Uncertainty and auto-correlation in. Measurement

Uncertainty and auto-correlation in. Measurement Uncertanty and auto-correlaton n arxv:1707.03276v2 [physcs.data-an] 30 Dec 2017 Measurement Markus Schebl Federal Offce of Metrology and Surveyng (BEV), 1160 Venna, Austra E-mal: markus.schebl@bev.gv.at

More information

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN Int. J. Chem. Sc.: (4), 04, 645654 ISSN 097768X www.sadgurupublcatons.com COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN R. GOVINDARASU a, R. PARTHIBAN a and P. K. BHABA b* a Department

More information

Work is the change in energy of a system (neglecting heat transfer). To examine what could

Work is the change in energy of a system (neglecting heat transfer). To examine what could Work Work s the change n energy o a system (neglectng heat transer). To eamne what could cause work, let s look at the dmensons o energy: L ML E M L F L so T T dmensonally energy s equal to a orce tmes

More information

Research on Time-history Input Methodology of Seismic Analysis

Research on Time-history Input Methodology of Seismic Analysis Transactons, SRT 19, Toronto, August 2007 Research on Tme-hstory Input ethoology of Sesmc Analyss Jang Nabn, ao Qng an Zhang Yxong State Key Laboratory of Reactor System Desgn Technology, Nuclear Power

More information

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. EE 539 Homeworks Sprng 08 Updated: Tuesday, Aprl 7, 08 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. For full credt, show all work. Some problems requre hand calculatons.

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

find (x): given element x, return the canonical element of the set containing x;

find (x): given element x, return the canonical element of the set containing x; COS 43 Sprng, 009 Dsjont Set Unon Problem: Mantan a collecton of dsjont sets. Two operatons: fnd the set contanng a gven element; unte two sets nto one (destructvely). Approach: Canoncal element method:

More information

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method Stablty Analyss A. Khak Sedgh Control Systems Group Faculty of Electrcal and Computer Engneerng K. N. Toos Unversty of Technology February 2009 1 Introducton Stablty s the most promnent characterstc of

More information

Classical Mechanics Symmetry and Conservation Laws

Classical Mechanics Symmetry and Conservation Laws Classcal Mechancs Symmetry an Conservaton Laws Dpan Kumar Ghosh UM-DAE Centre for Excellence n Basc Scences Kalna, Mumba 400085 September 7, 2016 1 Concept of Symmetry If the property of a system oes not

More information