Passive Bilateral Teleoperation with Constant Time Delays
|
|
- Quentin Patrick
- 5 years ago
- Views:
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 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 informationAnalytical 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 informationHigh-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 informationNew 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 informationGlobally 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 informationLarge-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 informationHybrid 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 informationWHY 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 informationAn 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 informationChapter 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 informationDigital 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 informationYukawa 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 informationCompetitive 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 informationThe 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 informationTransfer 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 informationCHAPTER 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 informationNeuro-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 informationBounds 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 informationImproved 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 informationAppendix 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 informationTHE 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 informationModule 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 informationChapter 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 informationDynamic 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 informationAdaptive 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 informationAn 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 informationA 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 informationThe 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 informationECE559VV 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 information3.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 informationIMECE 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 informationKinematics 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 informationDistributed 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 informationOn 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 informationIntroduction 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 informationNMT 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 informationVisualization 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 informationHigh 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 informationNon-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 informationIterative 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 informationNumerical 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 informationAPPENDIX 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 informationChapter 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 informationENTROPIC 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 informationOutput 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 informationprinceton 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 informationAdvanced 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 informationWeek3, 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 informationExplicit 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 informationJournal 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 informationProblem 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 informationPhysics 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 informationLecture 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 informationNote 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 informationAdaptive 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 information829. 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 informationELASTIC 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 informationPhysics 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 informationAdaptive 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 informationRobust 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 informationYong 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 informationAppendix 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 informationDesign 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 informationSome 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 informationFirst 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
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 informationTemperature. 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 informationOn 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 informationA 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 informationHongyi 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 informationTracking 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 informationPerron 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 informationMechanics 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 informationPhysics 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 informationComplex 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 informationComposite 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 informationGENERIC 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 informationThe 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 informationNON-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 informationA 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 informationCSci 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 informationBounds 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 informationChapter 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 informationParameter 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 informationTornado 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 informationOutline. 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 informationPhysics 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 informationControl 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 informationPHZ 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 information12. 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 informationUncertainty 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 informationCOEFFICIENT 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 informationWork 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 informationResearch 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 informationDO 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 informationStructure 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 informationIndeterminate 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 informationfind (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 informationIntroduction. - 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 informationClassical 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