VISUAL SERVOING with applications in medical robotics

Size: px
Start display at page:

Download "VISUAL SERVOING with applications in medical robotics"

Transcription

1 3 rd SUMMER EUROPEAN UNIVERSITY Surgial robotis Montpellier, September 2007 VISUAL SERVOING with appliations in medial robotis Florent Nageotte, Mihel de Mathelin Louis Pasteur Uniersity Strasbourg LSIIT Control, Vision & Robotis Researh Group UMR CNRS

2 Oeriew Part I : Fundamentals of isual seroing Bakground and definitions Seroing arhitetures and lassifiation Position-based isual seroing Image-based isual seroing Visual seroing without feature extration Part II : Medial robotis appliations Laparosopi surgery Internal organ motion traking 2

3 Oeriew Part I : Fundamentals of isual seroing Bakground and definitions Seroing arhitetures and lassifiation Position-based isual seroing Image-based isual seroing Visual seroing without feature extration Part II : Medial robotis appliations Laparosopi surgery Internal organ motion traking 3

4 Visual seroing priniple Control the end-effetor of a robot using a ision sensor (2D sensor). E.g : position the effetor of a robot w.r.t. an objet Similar to registration but the eloity of the robot is generally ontinuously updated Minimize a task funtion dependent of the error between the urrent pose of the robot and the referene pose Creates a irtual link between robot and objet Image proessing, omputer ision and ontrol issues 4

5 I.1 Bakground and definitions A. Coordinates and pose Coordinates of point P with respet to oordinate frame i : i P Position and orientation of frame i with respet to frame j : Pose j p i Tx T y T z α β γ translation etor rotation angles j j T R i i origin of frame i w.r. to frame j rotation matrix 5

6 6 B. Coordinate transformations Coordinates of point i P with respet to oordinate frame j : Coordinates of etor i V with respet to frame j : Homogeneous transformation from frame i to frame j : I.1 Bakground I.1 Bakground and and definitions definitions P R T P i i j i j j + V R V i i j j V H V P H P T R H i i j j i i j j i j i j i j i k k j i j H H H

7 7 C. Veloity of a rigid objet Veloity srew of frame i with respet to frame j in frame j oordinates: Veloity of point P rigidly attahed to frame i with respet to frame j expressed in frame j oordinates : I.1 Bakground I.1 Bakground and and definitions definitions ( ) ( ) ( ) ( ) i j j i j i i j i j j i j j j i j j j V T P R V P P + + Ω + Ω ) ( ( ) ( ) ( ) i j j i j j z y x z y x i j j V r Ω eloity rotational eloity translatio nal ω ω ω

8 8 D. Camera projetion model Perspetie projetion : oordinates of point P with respet to the amera frame I.1 Bakground I.1 Bakground and and definitions definitions image plane P Y X Z λ o x y u optial enter amera frame + + u C z y k z x k u u z y z x y x z y x P λ λ λ 0 0 p λ foal length P pixels

9 9 Intrinsi amera parameters obtained by alibration Mathematial model of a perfet optial system / physial phenomenon (distorsions, et.) Numerial sensor : aliasing, dynamial behaiour (limited bandwidth) Other models apply for other types of «isual» sensors, e.g., C-arm, CTsan, ultrasound probe, I.1 Bakground I.1 Bakground and and definitions definitions D. Camera projetion model u z y x K z y x k u k u z λ λ

10 Oeriew Part I : Fundamentals of isual seroing Bakground and definitions Seroing arhitetures and lassifiation Position-based isual seroing Image-based isual seroing Visual seroing without feature extration Part II : Medial robotis appliations Laparosopi surgery Internal organ motion traking 10

11 I.2 Classifiation I.21 Camera position A. Eye-in-hand onfiguration H t b H e e H e H must be known (alibration) or the referene position must be learned (showing) 11

12 I.2 Classifiation I.21 Camera position B. External amera onfiguration H e must be measured or b H must be known b H e b H H e H t 12

13 I.2 Classifiation I.22 Control leel A. Indiret isual seroing Cartesian ontroller p or r Trajetory generator q Joint ontrol Power amplifiers Control law - joint angles q Image proessing T Robot+amera image Suitable for slow isual seroing (T -1 < 50 Hz) Control law is easier to design 13

14 I.2 Classifiation I.22 Control leel B. Diret isual seroing - Control law τ or q Low leel joint ontrollers q Image proessing T Robot+amera image Suitable for fast isual seroing (T -1 > 50 Hz) Control law design is more omplex (robot dynamis must be taken into aount) 14

15 I.2 Classifiation I.23 Feedbak ariables A. Position-based isual seroing (3D isual seroing) p d - Control Power amplifiers p Robot+amera 15 Pose estimation Features extration A model of the objet must be known or multiple images should be used Calibration errors may indue large pose estimation errors Control law design is easier Possible loss of target for large errors f T image

16 I.2 Classifiation I.23 Feedbak ariables B. Image-based isual seroing (2D isual seroing) f d - Control Power amplifiers 16 Features extration Smaller omputational burden Eliminates errors due to alibration More omplex ontrol law Workspae limits an be hit for large errors f T Robot+amera image

17 I.2 Classifiation I.23 Feedbak ariables C. Hybrid isual seroing (e.g. 2D1/2 isual seroing) f pd p pd - - Hybrid ontrol law Power amplifiers p p f p Partial pose estimation Features extration Smaller errors due to alibration Simplified model of the target Better properties of motion (in the image and in 3D) f T Robot+amera image 17

18 I.2 Classifiation I.24 Bandwidth of the isual sero-loop A. Slow isual seroing Sampling frequeny < 50Hz Indiret isual seroing Robot transfer funtion model without dynamis Proportional ontrol law (P) B. Fast isual seroing Sampling frequeny > 50 Hz Diret isual seroing Dynamial model of the robot must be taken into aount More adaned ontrol laws : PID, preditie, robust, non-linear, 18

19 Oeriew Part I : Fundamentals of isual seroing Bakground and definitions Seroing arhitetures and lassifiation Position-based isual seroing Image-based isual seroing Visual seroing without feature extration Part II : Medial robotis appliations Laparosopi surgery Internal organ motion traking 19

20 I.3 Position-based isual seroing I.31 Indiret isual seroing A. Control law p d e p - Control law p or r Cartesian ontroller pˆ p + δp Pose estimation fˆ joint angles q Features extration T Robot+amera image 20 e p generally expressed as translation etor and Look-then-moe strategy (T ery large, asynhronous) : Pseudo-ontinuous strategy : p T J r Control law : p θu p ep r T k J 1 p ep

21 I.4 Image-based isual seroing I.41 Indiret isual seroing B. Obtaining referene pd Referene position of the end-effetor effetor is gien w.r.t an objet Computation of the orresponding position of the amera w.r.t the objet (eye-in-hand) requires the knowledge of e H Learning of the image in the final position and pose estimation of the amera 21

22 22 Case of a eye-in-hand system u T e p θ I.3 I.3 Position Position-based based isual isual seroing seroing I.31 Indiret isual I.31 Indiret isual seroing seroing C. Interation matrix Jp Ω Ω + Ω Ω e e e e e e e e e e p V R R EC R R R R R R EC V R R C e 0 ) ( ) (

23 D. Pose estimation Camera alibration : Tsaï (IEEE Trans. Rob. Aut., 1987), Zhang ( ICCV 99) Pose estimation : Analytial : 3 or 4 points : Tsaï (o-planar target) Numerial Iteratie : I.3 Position-based isual seroing I.31 Indiret isual seroing - DeMenthon (IEEE Trans. PAMI 1992, Int. J. Comp. Vision 1995), - VVS (Marhand, Eurographis 2002), - Minimization of reprojetion error (gradient, Leenberg-Marquardt, et.) 23

24 I.3 Position-based isual seroing I.31 Indiret isual seroing E. Stability and robustness Stability is not an issue : Look-then-moe strategy is always stable and low speed ision loop: r r p k J T p J T p 1 e p Exponential onergene J p depends on amera position w.r. to end-effetor (eye-in-hand onfiguration) Stability if J J ˆ 1 P P > 0 Measurement error is an issue: an be ery large! δp Main soures of unertainty : Improement: use learning Camera intrinsi parameters of p by showing when possible Camera position w.r. to end-effetor if eye-in-hand onfiguration Camera position w.r. to the robot base and robot kinemati hain if external amera, exept if end-effetor pose is estimated by ision Feature detetion error : δf f ˆ f 24

25 I.3 Position-based isual seroing I.32 Diret isual seroing A. Control law p d e p or Control - law q τ Low-leel joint ontrollers pˆ p + δp Robot+amera Pose estimation Features extration T image p p R q Control law : 1 T 1 q k J R ( q) J p e p T T p J r J J ( q) 25

26 I.3 Position-based isual seroing I.32 Diret isual seroing B. Stability and robustness (1) Stability may be an issue: Low speed ision loop : q q p kj T p J R ( q) J 1 R ( q) J T p 1 e p Exponential onergene High speed ision loop : Linearized approah : 1 p J s T p J R ( q) F( s, q) q q ( s) F( s, q) q ( s) Joint-leel eloity feedbak loops hae a linearizing and deoupling effet 1 R T 1 p q J ( q) J p Control law: with using a LPV disrete-time model of the ision loop p omputed 26 This approah works in pratie with 6DOF ision loop!

27 I.3 Position-based isual seroing I.32 Diret isual seroing B. Stability and robustness (2) LPV disrete-time model p d ontrol algorithm pˆ + p δp p J 1 1 ( ) T R q J p aquisition and image proessing Tz d 1 z 1 p q J T p F(z, q) robot identified model J R (q) q GPC of a 6DOF robot ision loop : J. Gangloff & M. de Mathelin (Adaned Robotis, ol 17, no 10, dé. 2003) 27

28 I.3 Position-based isual seroing I.32 Diret isual seroing B. Stability and robustness (3) Non linear approah : rigid link robot manipulator model τ M ( q) q + C( q, q ) q + g( q) + fr ( q, q ) Inertia Coriolis, entripetal graity frition -PD ontrol sheme(arimoto): τ g( q) K q K e with - Passiity-based sheme (Slotine & Li): τ M ( q) ζ + C( q, q ) ζ + g( q ) K e ζ q q K p p e q Λ d e q Problems: - asymptoti stability is proen with no frition - qd and q d are generally unknown - joint flexibilities and baklash q e q q q d High omplexity for 6DOF! 28

29 Oeriew Part I : Fundamentals of isual seroing Bakground and definitions Seroing arhitetures and lassifiation Position-based isual seroing Image-based isual seroing Visual seroing without feature extration Part II : Medial robotis appliations Laparosopi surgery Internal organ motion traking 29

30 I.4 Image-based isual seroing I.41 Indiret isual seroing A. Control law f d e f - Control law p or r Cartesian ontroller q fˆ f + δf Features extration pixel noise T Robot+amera image Look-then-moe strategy : selet or to derease a ost funtion of Pseudo-ontinuous strategy : p f J I r q Control law : r k + e f J I e f 30

31 B. Obtaining referene d I.4 Image-based isual seroing I.41 Indiret isual seroing f Computation of the perspetie projetion of the objet in the referene position requires a model of the objet requires the knowledge of e H requires an aurate model of amera Learning of the image in the final position 31

32 I.4 Image-based isual seroing I.41 Indiret isual seroing C. Image jaobian or interation matrix Image Jaobian J I : p x n matrix (p feature etor dimension, n # dof of effetor). Depends on kind of features used. Exemple : ase of a point P seen by a amera attahed to the effetor of a robot u P X C Z C Z C Z O P X C Y C X O Y O 32

33 33 ( ) z y x z y x o r ω ω ω ( ) ( ) o o V P P + Ω z y x P z y u x y x z u x z y k u u k z z k z u u k z z z y x P ) ) ( ) ( ( ) ( ) ( ω ω λ ω ω λ ω λ ω image oordinates in pixel: (u,) I.4 I.4 Image Image-based based isual isual seroing seroing I.41 Indiret isual I.41 Indiret isual seroing seroing ) ( ) ( ) )( ( ) ( ) ( ) ( 0 ) ( ) ( ) ( ) )( ( ) ( o u u u u u u r k u u k k u u k k z z k k k k u u k k u u z u u z k u + + λ λ λ λ λ λ λ λ Ω Ω e e I e e e e e o V J u f V R EC R R r 0 ) (

34 I.4 Image-based isual seroing I.41 Indiret isual seroing J I depends on intrinsi parameters and on pose! Also depends on pose of amera w.r.t effetor J I must be estimated Ĵ I J I must be full rank for ontrolling all degrees of freedom (at least 3 points required) End-effetor with n dofs only n features omponents an be ontrolled in the image. Expressions of image jaobian for lassial geometrial objets : straight lines, spheres, ylinders et. in F. Chaumette, thesis 90 and Espiau, TRA 92 and for any planar objets using moments in Tahri O and Chaumette, F, IEEE TR 2005, 21(6) 34

35 I.4 Image-based isual seroing I.41 Indiret isual seroing D. Stability and robustness (1) Low speed : f k J I J r + I e f r Exponential onergene Behaes as a gradient desent optimization : may be stuk in loal minima! 2 nd order optimization Malis, E., Improing isionbased ontrol using effiient seond-order minimization tehniques, 2004 The norm of e f dereases, but some omponents may inrease during motion (possible loss of features) If more features than dofs, least square of error minimization no ontrol in the artesian spae better behaiour for small errors 35

36 I.4 Image-based isual seroing I.41 Indiret isual seroing D. Stability and robustness (2) Main soure of unertainty : Ji f + k J Jˆ e Exponential onergene if ˆ + J I J I > 0 Easily guaranteed : oarse alibration and oarse pose estimation are generally suffiient. No proof of global stability in funtion of k Error in the final position independent of errors on amera / effetor position independent of erros on intrinsi parameters of the amera if referene image has been learned funtion of pixel noise Pixel noise attenuation : pik more features smaller Problem of reahing the workspae limits Hybrid isual seroing Malis, E., 21/2D isual seroing TRA 99 I δf I f 36

37 I.4 Image-based isual seroing I.42 Diret isual seroing A. Control law f d e f Control - law q or τ Low-leel joint ontrollers fˆ f + δf Features extration pixel noise T Robot+amera image f J r J J ( q ) q I I R Control law : q k J 1 + R ( q) J I e f ˆ 37

38 I.4 Image-based isual seroing I.42 Diret isual seroing B. Stability and robustness (1) Stability may be an issue: Low speed ision loop : q q f High speed ision loop : Linearized approah : 1 f J s I I J R R 1 R k J J ( q) J ( q) Jˆ ( q) F( s, q) q + I e f q ( s) Exponential onergene if F( s, q) q ( s) ˆ > 0 + J I J I 1 q J R ( q) r Control law: with using a LPV disrete-time model of the ision loop r omputed This approah works in pratie with 6DOF ision loops! 38

39 I.4 Image-based isual seroing I.42 Diret isual seroing B. Stability and robustness (2) LPV disrete-time model image-based identifiation e p + Ĵ I ontrol algorithm e f + - f d fˆ r δf + f J 1 R ( q) q aquisition and image proessing Tz d 1 z 1 f J I F(z, q) robot identified model J R (q) q GPC of a 6DOF robot ision loop : J. Gangloff & M. de Mathelin (Adaned Robotis, ol 17, no 10, dé. 2003) 39

40 I.4 Image-based isual seroing I.42 Diret isual seroing B. Stability and robustness (3) High speed ision loop : Non linear approah : rigid link robot manipulator model τ M ( q) q + C( q, q ) q + g( q) + fr ( q, q ) PD ontrol sheme : (R. Kelly, IEEE Trans. Rob. Aut., ol 12, 1996) T R τ g( q) K q J ( q) K Jˆ p + I e f Stability is proed only for a 2 DOF robot with no frition > Preious restritions apply 40

41 Oeriew Part I : Fundamentals of isual seroing Bakground and definitions Seroing arhitetures and lassifiation Position-based isual seroing Image-based isual seroing Visual seroing without feature extration Part II : Medial robotis appliations Laparosopi surgery Internal organ motion traking 41

42 Case of planar objets Images linked by a homography u p Gp 1 u G 1 I.5 Visual seroing without feature extration g G g g with det(g) 1 Computation of G using SSD minimization Example : ESM (S. Benhimane and E. Malis 2004) Computation of eulidean 1 homography m K GKm Hm g g g g g g

43 Control law I SSD + r ESM I.5 Visual seroing without feature extration H Control law r Cartesian ontroller I I T Robot+amera image r λν 0 0 e λ ω e ν ω with e e e ν ω e ν ( H I) m and [ ] T and e ω H H m 1 K p Camera Camera oordinates of a seleted point in the referene image 43

44 Stability and robustness I.5 Visual seroing without feature extration S. Benhimane and E. Malis, Homography-based 2D isual seroing, IEEE ICRA 2006 No need of interation matrix omputation Proof of loal stability Pratially : good robustness w.r.t intrinsi parameters Large domain of stability Reently extended to non-planar objets 44

45 Bibliography Hager, G. and S. Huthinson. Speial issue on ision-based ontrol of robot manipulators. IEEE Trans. Rob. Autom., ol 12, no 5, Corke, P. Visual ontrol of robots. Researh studies Press Ltd., Somerset, UK, B. Espiau and F. Chaumette, A new approah to isual seroing in robotis, IEEE Transations on robotis and automation, ol 8, no 3, 1992, S. Benhimane and E. Malis, Homography-based 2D isual seroing, IEEE ICRA

46 Oeriew Part I : Fundamentals of isual seroing Bakground and definitions Seroing arhitetures and lassifiation Position-based isual seroing Image-based isual seroing Visual seroing without feature extration Part II : Medial robotis appliations Laparosopi surgery Internal organ motion traking 46

47 Visual seroing in medial robotis Main appliations In laparosopi surgery In ophtalmology For physiologial motion ompensation In US imaging diagnosis and interention M.A.Vitrani and G.Morel, ICRA 2005 : guidane of surgial instruments A. Krupa and F. Chaumette, IROS 2005 : ontrol of US probe 47

48 Laparosopi surgery Small inisions, long-shaft instruments Endosopi ision system Adantages : Shorter reoery time Lower infetion risk Smaller osts Diffiulties : Tiring gestures Indiret isual feedbak, lak of depth information Inerted motions limited to 4 DOF Poor hapti feedbak 48

49 Laparosopi surgery Classial roboti ontrol of the endosope AESOP (Computer Motion/Intuitie Surgial) EndoAssist (Armstrong-Healthare) Voie ontrolled Controlled with head motions 49

50 Laparosopi surgery Vision-based ontrol of the endosope R. Taylor (95), A. Casals (95), G. Wei and G. Hirzinger (97), S. Voros and Ph. Cinquin (2006), et. More ision guided systems than real isual seroing Centering the image on marked instruments Traking of instruments Centering the image on speifi physiologial markers 50

51 Laparosopi surgery Visual seroing based ontrol of the instrument (Krupa 2003) Goal : surgial instrument guidane Autonomous reoery of the instruments Bring the instrument in the field of iew Centering of the instrument in the image Automati positionning Put the instrument at a desired loation w.r.t to tissues Proide depth informations Motiations : - inrease safety - improe ergonomy - speed up surgial proedures 51

52 Laparosopi surgery Experimental setup Surgial robot amera optis 52 instrument and Laser pointing troart 1 Laser spots organ troart 2 External amera onfiguration Laser projetion to add information to the sene

53 Laparosopi surgery Experimental setup laser pointing deie 4 laser spots To be used with standard instruments (4 mm) and standard troarts (10 mm) 3 optial markers (detetion of the instrument) 53

54 Laparosopi surgery Kinematis in laparosopi surgery - onstraints: 4 DOF W ( ω ω ω ) op x y z z - AESOP: 6 DOF (2 passie joints) Τ ω x ω y q2 q θ z ω z d 1 q q 4 q q 6 kinemati model : W J ( q, d ) q op op 1 q J q ˆ W 1 op(, d1) op 54

55 Laparosopi surgery Visual features detetion Goal: to obtain the image oordinates of the enters of laser spots and optial markers refleting lights diffused laser pattern 55

56 Robust detetion Laparosopi surgery Blinking laser spots and optial markers synhronized with frame aquisition Een frames : laser on markers off Odd frames : markers on laser off Seletion by high-pass filtering laser marker ligne N : Not a feature 56

57 Robust detetion Laparosopi surgery 57

58 Laparosopi surgery Control of the position of the spot 2D diret isual seroing sheme (Image-based) z s p ω + - Z ˆ B { B} k s p 1 J ω ω z J op ( q, d1) 1 q Robot ontroller AESOP 58 Xˆ B Yˆ B sp O Q - ontrolled ariable : S ( ) T p up p - ontrol input : ω ( ω ) T x ωy - ontrol law : p k ( p p) Features detetion S S S J ω ω image amera

59 Estimation of interation matrix Open loop estimation of Hypothesis : standstill amera T T Laparosopi surgery Stability if ω x y x st ω 0 ω 0 Jˆ J ω u ω11 ω x Jˆ Jˆ Jˆ ω11 ω12 ω Jˆ ˆ ω21 Jω22 J ω st ω12 ω 22 y J Jˆ 0 ω 1 ω > Jˆ p T u p ω T y Jˆ Jˆ ω 21 ω x y p T p ω T 59

60 Laparosopi surgery 60

61 Laparosopi surgery Autonomous omplex gestures (F. Nageotte 2005) Automati stithing in laparosopi surgery Trajetory planning in the amera frame Trajetory following using the endosopi amera Aoid registration proedure between sensors 61

62 Laparosopi surgery Autonomous omplex gestures (F. Nageotte 2005) 2D isual seroing : Marked instrument Referene images : projetion of the needle-holder in the endosopi image fine amera alibration required Controlled ariables : marker points + ylinder (10 to 16) Control input : 4 DOFs of end-effetor 62

63 Laparosopi surgery Autonomous omplex gestures (F. Nageotte 2005) Complex interation matrix with arying size. Hypothesis : small traking error using known referene pose required matrix omputed no pose estimation Predition of feature appearane and disappearane for update 63

64 Laparosopi surgery Automati suturing 64

65 Oeriew Part I : Fundamentals of isual seroing Bakground and definitions Seroing arhitetures and lassifiation Position-based isual seroing Image-based isual seroing Visual seroing without feature extration Part II : Medial robotis appliations Laparosopi surgery Internal organ motion traking 65

66 Physiologial motions ompensation Objeties Filter Stabilized image Visual feedbak Controller 66

67 Physiologial motions ompensation Breathing motion ompensation (R. Ginhoux 2003) Preditie ontrol law : take into aount periodi disturbane GPC ontrol law with periodi output perturbation model Dynami model of the robot is required : linearized model obtained by pseudo-random exitation Controlled ariable : distane to the organ Control input : first axis of the robot (partiular onfiguration) Identifiation of the interation gain 67

68 Physiologial motions ompensation In io experiment 68

69 Physiologial motions ompensation Beating heart motion ompensation (R. Ginhoux 2003) Fast robot Fast amera (500hz) In io experiments 2 types of motion: - slow (~0.25 Hz) - fast (~1.5 Hz) 69

70 Physiologial motions ompensation Beating heart motion ompensation Heart motion Respiratory omponent Adaptie filtering of 12 harmonis of the heart beat rate 70

71 Physiologial motions ompensation Beating heart motion ompensation Predition of the perturbation referene orretion GPC ontrol law 3 ontrolled ariables : distane to the organ in the image z, spot position on the organ surfae (x,y) Control input : 3 axes of the robot (q1, q2, q3) Linearized dynami model of the axes of the robot onsidered as dynamially deoupled Hypothesis : interation matrix partly deoupled x 0 y 0 z L 31 L L L13 q 1 L 23 q 2 0 q 3 Estimated around working position in open loop 71

72 In io experiment Physiologial motions ompensation 72

73 Physiologial motions ompensation Breathing motion ompensation with a flexible endosope Flexible endosopes : gastrosopy, olosopy, transluminal operations Hard to manually ontrol, orientation and translation oordination ery omplex Baklashes, dead zones Autonomous following of the organs 73

74 Physiologial motions ompensation Breathing motion ompensation with a flexible endosope 74

75 Physiologial motions ompensation Breathing motion ompensation with a flexible endosope Controlled ariables : position of an interest area in the image (u, ) Control input : q 1 and q 2 u L11 L12 q 1 Interation matrix estimated beforehand q Repetitie ontrol law rejetion of breathing disturbane and model errors after some disturbane periods L 21 L

76 Physiologial motions ompensation Breathing motion ompensation with a flexible endosope 76

77 Bibliography A. Krupa. Visual seroing in laparosopi surgery. Ph.D. thesis, Uniersity Louis Pasteur, Strasbourg I, june 2003 R. Ginhoux. Compensation des mouements physiologiques en hirurgie robotisée par ommande préditie. Ph.D. Thesis, Strasbourg Uniersity, deember 2003 F. Nageotte. Contributions à la suture assistée par ordinateur en hirurgie mini-inasie. Ph.D. Thesis, Strasbourg Uniersity, noember 2005 A. Krupa, J. Gangloff, C. Doignon, M. de Mathelin, G. Morel, J. Leroy, L. Soler, J. Maresaux. Autonomous 3D positioning of surgial instruments in robotized laparosopi surgery using isual seroing. IEEE Transations on Robotis, ol 19, no 5, pages , 2003 R. Ginhoux, J. Gangloff, M. de Mathelin, L. Soler, M. Arenas Sanhez and J. Maresaux. Atie Filtering of Physiologial Motion in Robotized Surgery Using Preditie Control. IEEE Transations on Robotis, ol 21, no 1, pages 67-79, 2005 L. Ott, Ph. Zanne, F. Nageotte and M. De Mathelin, Problemati of transluminal operations, Surgetia 2007 F. Nageotte, Ph. Zanne, M. De Mathelin and Ch. Doignon, Visual seroing based traking for endosopi path following, IROS

78 Perspeties New sensors (isual seroing under CT san?) Seroing without features extration on deformable objets High speed isual seroing with high resolution 78

Force/Vision Based Active Damping Control of Contact Transition in Dynamic Environments

Force/Vision Based Active Damping Control of Contact Transition in Dynamic Environments Fore/Vision Based Ative Damping Control of Contat Transition in Dynami Environments Tomas Olsson, Rolf Johansson, Anders Robertsson Department of Automati Control, Lund Inst. of Teh., Lund University,

More information

Speed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion

Speed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion Speed-feedbak Diret-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion Y. Yamamoto, T. Nakamura 2, Y. Takada, T. Koseki, Y. Aoyama 3, and Y. Iwaji 3

More information

Control Theory association of mathematics and engineering

Control Theory association of mathematics and engineering Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology

More information

Control of industrial robots. Control of the interaction

Control of industrial robots. Control of the interaction Control o industrial robots Control o the interation Pro. Paolo Roo (paolo.roo@polimi.it) Politenio di Milano Dipartimento di Elettronia, Inormazione e Bioingegneria Introdution So ar we have assumed that

More information

Visual Servoing Based on Multirate Sampling Control

Visual Servoing Based on Multirate Sampling Control (), () Visual ervoing Based on Multirate ampling Control tability Improvement by Intersample Observer and Periodi Disturbane Rejetion Hiroshi Fujimoto (Nagaoka University of Tehnology), Yoihi Hori (The

More information

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install

More information

We are IntechOpen, the first native scientific publisher of Open Access books. International authors and editors. Our authors are among the TOP 1%

We are IntechOpen, the first native scientific publisher of Open Access books. International authors and editors. Our authors are among the TOP 1% We are IntehOpen, the first natie sientifi publisher of Open Aess books 3,35 18, 1.7 M Open aess books aailable International authors and editors Downloads Our authors are among the 151 Countries deliered

More information

NUMERICAL SIMULATION OF ATOMIZATION WITH ADAPTIVE JET REFINEMENT

NUMERICAL SIMULATION OF ATOMIZATION WITH ADAPTIVE JET REFINEMENT Paper ID ILASS8--7 ILASS 28 Sep. 8-, 28, Como Lake, Italy A44 NUMERICAL SIMULATION OF ATOMIZATION WITH ADAPTIVE JET REFINEMENT Anne Bagué, Daniel Fuster, Stéphane Popinet + & Stéphane Zaleski Université

More information

The Lorenz Transform

The Lorenz Transform The Lorenz Transform Flameno Chuk Keyser Part I The Einstein/Bergmann deriation of the Lorentz Transform I follow the deriation of the Lorentz Transform, following Peter S Bergmann in Introdution to the

More information

An Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems

An Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems An Integrated Arhiteture of Adaptive Neural Network Control for Dynami Systems Robert L. Tokar 2 Brian D.MVey2 'Center for Nonlinear Studies, 2Applied Theoretial Physis Division Los Alamos National Laboratory,

More information

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry 9 Geophysis and Radio-Astronomy: VLBI VeryLongBaseInterferometry VLBI is an interferometry tehnique used in radio astronomy, in whih two or more signals, oming from the same astronomial objet, are reeived

More information

FREQUENCY DOMAIN FEEDFORWARD COMPENSATION. F.J. Pérez Castelo and R. Ferreiro Garcia

FREQUENCY DOMAIN FEEDFORWARD COMPENSATION. F.J. Pérez Castelo and R. Ferreiro Garcia FREQUENCY DOMAIN FEEDFORWARD COMPENSATION F.J. Pérez Castelo and R. Ferreiro Garia Dept. Ingeniería Industrial. Universidad de La Coruña javierp@ud.es, Phone: 98 7.Fax: -98-7 ferreiro@ud.es, Phone: 98

More information

Grasp Planning: How to Choose a Suitable Task Wrench Space

Grasp Planning: How to Choose a Suitable Task Wrench Space Grasp Planning: How to Choose a Suitable Task Wrenh Spae Ch. Borst, M. Fisher and G. Hirzinger German Aerospae Center - DLR Institute for Robotis and Mehatronis 8223 Wessling, Germany Email: [Christoph.Borst,

More information

Control of industrial robots

Control of industrial robots Control o industrial robots Control with eternal sensors Pro. Paolo Roo (paolo.roo@polimi.it) Politenio di Milano Dipartimento di Elettronia, Inormazione e Bioingegneria Introdution So ar we have assumed

More information

A Spatiotemporal Approach to Passive Sound Source Localization

A Spatiotemporal Approach to Passive Sound Source Localization A Spatiotemporal Approah Passive Sound Soure Loalization Pasi Pertilä, Mikko Parviainen, Teemu Korhonen and Ari Visa Institute of Signal Proessing Tampere University of Tehnology, P.O.Box 553, FIN-330,

More information

Radiation Transport in Castro

Radiation Transport in Castro Radiation Transport in Castro May 9 Louis Howell P. O. Box 88 Liermore CA 9455 This work performed under the auspies of the U.S. Department of Energy by under Contrat DEAC57NA744 LLNLPRES4597 The radiation

More information

Finite-time stabilization of chaotic gyros based on a homogeneous supertwisting-like algorithm

Finite-time stabilization of chaotic gyros based on a homogeneous supertwisting-like algorithm OP Conferene Series: Materials Siene Engineering PAPER OPEN ACCESS Finite-time stabilization of haoti gyros based on a homogeneous supertwisting-like algorithm To ite this artile: Pitha Khamsuwan et al

More information

Chapter 35. Special Theory of Relativity (1905)

Chapter 35. Special Theory of Relativity (1905) Chapter 35 Speial Theory of Relatiity (1905) 1. Postulates of the Speial Theory of Relatiity: A. The laws of physis are the same in all oordinate systems either at rest or moing at onstant eloity with

More information

The Thomas Precession Factor in Spin-Orbit Interaction

The Thomas Precession Factor in Spin-Orbit Interaction p. The Thomas Preession Fator in Spin-Orbit Interation Herbert Kroemer * Department of Eletrial and Computer Engineering, Uniersity of California, Santa Barbara, CA 9306 The origin of the Thomas fator

More information

the following action R of T on T n+1 : for each θ T, R θ : T n+1 T n+1 is defined by stated, we assume that all the curves in this paper are defined

the following action R of T on T n+1 : for each θ T, R θ : T n+1 T n+1 is defined by stated, we assume that all the curves in this paper are defined How should a snake turn on ie: A ase study of the asymptoti isoholonomi problem Jianghai Hu, Slobodan N. Simić, and Shankar Sastry Department of Eletrial Engineering and Computer Sienes University of California

More information

The gravitational phenomena without the curved spacetime

The gravitational phenomena without the curved spacetime The gravitational phenomena without the urved spaetime Mirosław J. Kubiak Abstrat: In this paper was presented a desription of the gravitational phenomena in the new medium, different than the urved spaetime,

More information

Maximum Entropy and Exponential Families

Maximum Entropy and Exponential Families Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It

More information

Task space trajectory tracking control of an underwater vehicle-manipulator system under ocean currents

Task space trajectory tracking control of an underwater vehicle-manipulator system under ocean currents Indian Journal of Geo-Marine Sienes Vol. 4(6), Otober 013, pp. 675-683 ask spae trajetory traking ontrol of an underwater ehile-anipulator syste under oean urrents M Santhakuar Disipline of Mehanial Engineering,

More information

A model for measurement of the states in a coupled-dot qubit

A model for measurement of the states in a coupled-dot qubit A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:

More information

Robust Flight Control Design for a Turn Coordination System with Parameter Uncertainties

Robust Flight Control Design for a Turn Coordination System with Parameter Uncertainties Amerian Journal of Applied Sienes 4 (7): 496-501, 007 ISSN 1546-939 007 Siene Publiations Robust Flight ontrol Design for a urn oordination System with Parameter Unertainties 1 Ari Legowo and Hiroshi Okubo

More information

Controller Design Based on Transient Response Criteria. Chapter 12 1

Controller Design Based on Transient Response Criteria. Chapter 12 1 Controller Design Based on Transient Response Criteria Chapter 12 1 Desirable Controller Features 0. Stable 1. Quik responding 2. Adequate disturbane rejetion 3. Insensitive to model, measurement errors

More information

Kinematics and inverse dynamics analysis for a general 3-PRS spatial parallel mechanism Yangmin Li and Qingsong Xu

Kinematics and inverse dynamics analysis for a general 3-PRS spatial parallel mechanism Yangmin Li and Qingsong Xu Robotia (2005) volume 23, pp. 219 229. 2005 Cambridge University Press doi:10.1017/s0263574704000797 Printed in the United Kingdom Kinematis and inverse dynamis analysis for a general 3-PRS spatial parallel

More information

Advances in Radio Science

Advances in Radio Science Advanes in adio Siene 2003) 1: 99 104 Copernius GmbH 2003 Advanes in adio Siene A hybrid method ombining the FDTD and a time domain boundary-integral equation marhing-on-in-time algorithm A Beker and V

More information

An Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances

An Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances An aptive Optimization Approah to Ative Canellation of Repeated Transient Vibration Disturbanes David L. Bowen RH Lyon Corp / Aenteh, 33 Moulton St., Cambridge, MA 138, U.S.A., owen@lyonorp.om J. Gregory

More information

A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS

A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS Vietnam Journal of Mehanis, VAST, Vol. 4, No. (), pp. A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS Le Thanh Tung Hanoi University of Siene and Tehnology, Vietnam Abstrat. Conventional ship autopilots are

More information

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite. Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need

More information

20 Doppler shift and Doppler radars

20 Doppler shift and Doppler radars 20 Doppler shift and Doppler radars Doppler radars make a use of the Doppler shift phenomenon to detet the motion of EM wave refletors of interest e.g., a polie Doppler radar aims to identify the speed

More information

QCLAS Sensor for Purity Monitoring in Medical Gas Supply Lines

QCLAS Sensor for Purity Monitoring in Medical Gas Supply Lines DOI.56/sensoren6/P3. QLAS Sensor for Purity Monitoring in Medial Gas Supply Lines Henrik Zimmermann, Mathias Wiese, Alessandro Ragnoni neoplas ontrol GmbH, Walther-Rathenau-Str. 49a, 7489 Greifswald, Germany

More information

An iterative least-square method suitable for solving large sparse matrices

An iterative least-square method suitable for solving large sparse matrices An iteratie least-square method suitable for soling large sparse matries By I. M. Khabaza The purpose of this paper is to report on the results of numerial experiments with an iteratie least-square method

More information

Contact State Estimation using Multiple Model Estimation and Hidden Markov Models

Contact State Estimation using Multiple Model Estimation and Hidden Markov Models Thomas J. Debus Pierre E. Dupont Department of Aerospae and Mehanial Engineering Boston University Boston, MA 02445, USA Robert D. Howe Division of Engineering and Applied Sienes Harvard University Cambridge,

More information

Optimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach

Optimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach Amerian Journal of heoretial and Applied tatistis 6; 5(-): -8 Published online January 7, 6 (http://www.sienepublishinggroup.om/j/ajtas) doi:.648/j.ajtas.s.65.4 IN: 36-8999 (Print); IN: 36-96 (Online)

More information

Instytut Fizyki Doświadczalnej Wydział Matematyki, Fizyki i Informatyki UNIWERSYTET GDAŃSKI

Instytut Fizyki Doświadczalnej Wydział Matematyki, Fizyki i Informatyki UNIWERSYTET GDAŃSKI Instytut Fizyki Doświadzalnej Wydział Matematyki, Fizyki i Informatyki UNIWERSYTET GDAŃSKI I. Bakground theory. 1. Laser doppler anemometry LDA (Doppler model): a) LDA priniples; b) single and multi-hannel

More information

Journal of Theoretics Vol.5-2 Guest Commentary Relativistic Thermodynamics for the Introductory Physics Course

Journal of Theoretics Vol.5-2 Guest Commentary Relativistic Thermodynamics for the Introductory Physics Course Journal of heoretis Vol.5- Guest Commentary Relatiisti hermodynamis for the Introdutory Physis Course B.Rothenstein bernhard_rothenstein@yahoo.om I.Zaharie Physis Department, "Politehnia" Uniersity imisoara,

More information

Planning with Uncertainty in Position: an Optimal Planner

Planning with Uncertainty in Position: an Optimal Planner Planning with Unertainty in Position: an Optimal Planner Juan Pablo Gonzalez Anthony (Tony) Stentz CMU-RI -TR-04-63 The Robotis Institute Carnegie Mellon University Pittsburgh, Pennsylvania 15213 Otober

More information

CSC2515 Winter 2015 Introduc3on to Machine Learning. Lecture 5: Clustering, mixture models, and EM

CSC2515 Winter 2015 Introduc3on to Machine Learning. Lecture 5: Clustering, mixture models, and EM CSC2515 Winter 2015 Introdu3on to Mahine Learning Leture 5: Clustering, mixture models, and EM All leture slides will be available as.pdf on the ourse website: http://www.s.toronto.edu/~urtasun/ourses/csc2515/

More information

Aircraft CAS Design with Input Saturation Using Dynamic Model Inversion

Aircraft CAS Design with Input Saturation Using Dynamic Model Inversion International Journal of Control, Automation, and Systems Vol., No. 3, September 003 35 Airraft CAS Design with Input Saturation sing Dynami Model Inversion Sangsoo Lim and Byoung Soo Kim Abstrat: This

More information

19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007 INVESTIGATION ON THE DYNAMIC BEHAVIOR OF A COMPOUND PLANETARY GEAR DRIVE

19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007 INVESTIGATION ON THE DYNAMIC BEHAVIOR OF A COMPOUND PLANETARY GEAR DRIVE 9 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, -7 SEPTEMBER 007 INVESTIGATION ON THE DYNAMIC BEHAVIOR OF A COMPOUND PLANETARY GEAR DRIVE PACS: 43.40.-r Dhouib, Seddik ; Hbaieb, Riadh ; Chaari, Fakher

More information

arxiv: v1 [physics.gen-ph] 5 Jan 2018

arxiv: v1 [physics.gen-ph] 5 Jan 2018 The Real Quaternion Relativity Viktor Ariel arxiv:1801.03393v1 [physis.gen-ph] 5 Jan 2018 In this work, we use real quaternions and the basi onept of the final speed of light in an attempt to enhane the

More information

AC : A GRAPHICAL USER INTERFACE (GUI) FOR A UNIFIED APPROACH FOR CONTINUOUS-TIME COMPENSATOR DESIGN

AC : A GRAPHICAL USER INTERFACE (GUI) FOR A UNIFIED APPROACH FOR CONTINUOUS-TIME COMPENSATOR DESIGN AC 28-1986: A GRAPHICAL USER INTERFACE (GUI) FOR A UNIFIED APPROACH FOR CONTINUOUS-TIME COMPENSATOR DESIGN Minh Cao, Wihita State University Minh Cao ompleted his Bahelor s of Siene degree at Wihita State

More information

EE 321 Project Spring 2018

EE 321 Project Spring 2018 EE 21 Projet Spring 2018 This ourse projet is intended to be an individual effort projet. The student is required to omplete the work individually, without help from anyone else. (The student may, however,

More information

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six

More information

The Special Theory of Relativity

The Special Theory of Relativity The Speial Theory of Relatiity Galilean Newtonian Relatiity Galileo Galilei Isaa Newton Definition of an inertial referene frame: One in whih Newton s first law is alid. onstant if F0 Earth is rotating

More information

A New Version of Flusser Moment Set for Pattern Feature Extraction

A New Version of Flusser Moment Set for Pattern Feature Extraction A New Version of Flusser Moment Set for Pattern Feature Extration Constantin-Iulian VIZITIU, Doru MUNTEANU, Cristian MOLDER Communiations and Eletroni Systems Department Military Tehnial Aademy George

More information

Journal of Physical Mathematics

Journal of Physical Mathematics Journal of Physial Mathematis Researh Artile Artile Journal of Physial Mathematis Makanae, J Phys Math 207, 8: DOI: 0.472/2090-0902.00025 OMICS Open International Aess Verifying Einstein s Time by Using

More information

Einstein s theory of special relativity

Einstein s theory of special relativity Einstein s theory of speial relatiity Announements: First homework assignment is online. You will need to read about time dilation (1.8) to answer problem #3 and for the definition of γ for problem #4.

More information

Numerical simulation of a one-dimensional shock tube problem at supercritical fluid conditions

Numerical simulation of a one-dimensional shock tube problem at supercritical fluid conditions International Journal of Physial Sienes Vol. 3 (1), pp. 314-30, Deember, 008 Available online at http://www.aademijournals.org/ijps ISSN 199-1950 008 Aademi Journals Full ength esearh Paper Numerial simulation

More information

23.1 Tuning controllers, in the large view Quoting from Section 16.7:

23.1 Tuning controllers, in the large view Quoting from Section 16.7: Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output

More information

Transient wave propagation analysis of a pantograph- catenary system

Transient wave propagation analysis of a pantograph- catenary system Journal of Physis: Conferene Series PAPER OPEN ACCESS Transient wave propagation analysis of a pantograph- atenary system To ite this artile: Kyohei Nagao and Arata Masuda 216 J. Phys.: Conf. Ser. 744

More information

τ = 10 seconds . In a non-relativistic N 1 = N The muon survival is given by the law of radioactive decay N(t)=N exp /.

τ = 10 seconds . In a non-relativistic N 1 = N The muon survival is given by the law of radioactive decay N(t)=N exp /. Muons on the moon Time ilation using ot prouts Time ilation using Lorentz boosts Cheking the etor formula Relatiisti aition of eloities Why you an t eee the spee of light by suessie boosts Doppler shifts

More information

STABLE DESIGN OF H OUTPUT FEEDBACK CONTROL SYSTEMS WITH TIME DELAY

STABLE DESIGN OF H OUTPUT FEEDBACK CONTROL SYSTEMS WITH TIME DELAY 3 th World Conferene on Earthquake Engineering Vanouver, B.C., Canada August -6, 4 Paper No. 44 SABLE DESIGN OF H OUPU FEEDBACK CONROL SYSEMS WIH IME DELAY Chi-Chang LIN, Chang-Ching CHANG, Huang-Lin CHEN

More information

LATTICE BOLTZMANN METHOD FOR MICRO CHANNEL AND MICRO ORIFICE FLOWS TAIHO YEOM. Bachelor of Science in Mechanical Engineering.

LATTICE BOLTZMANN METHOD FOR MICRO CHANNEL AND MICRO ORIFICE FLOWS TAIHO YEOM. Bachelor of Science in Mechanical Engineering. LATTICE BOLTZMANN METHOD FOR MICRO CHANNEL AND MICRO ORIFICE FLOWS By TAIHO YEOM Bahelor of Siene in Mehanial Engineering Ajou University Suwon, South Korea 2005 Submitted to the Faulty of the Graduate

More information

Stabilization of the Precision Positioning Stage Working in the Vacuum Environment by Using the Disturbance Observer

Stabilization of the Precision Positioning Stage Working in the Vacuum Environment by Using the Disturbance Observer Proeedings of the 4th IIAE International Conferene on Industrial Appliation Engineering 216 Stabilization of the Preision Positioning Stage Working in the Vauum Environment by Using the Disturbane Observer

More information

IDENTIFICATION AND CONTROL OF ACOUSTIC RADIATION MODES

IDENTIFICATION AND CONTROL OF ACOUSTIC RADIATION MODES IDENTIFICATION AND CONTROL OF ACOUSTIC RADIATION MODES Arthur P. Berkhoff University of Twente, Faulty of Eletrial Engineering, P.O. Box 217, 7 AE Enshede, The Netherlands email: a.p.berkhoff@el.utwente.nl

More information

Sensitivity Analysis in Markov Networks

Sensitivity Analysis in Markov Networks Sensitivity Analysis in Markov Networks Hei Chan and Adnan Darwihe Computer Siene Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwihe}@s.ula.edu Abstrat This paper explores

More information

Distance and Orientation Measurement of a Flat Surface by a Single Underwater Acoustic Transducer

Distance and Orientation Measurement of a Flat Surface by a Single Underwater Acoustic Transducer Distane and Orientation Measurement of a Flat Surfae by a Single Underwater Aousti Transduer Vinent Creuze To ite this version: Vinent Creuze. Distane and Orientation Measurement of a Flat Surfae by a

More information

Variation Based Online Travel Time Prediction Using Clustered Neural Networks

Variation Based Online Travel Time Prediction Using Clustered Neural Networks Variation Based Online Travel Time Predition Using lustered Neural Networks Jie Yu, Gang-Len hang, H.W. Ho and Yue Liu Abstrat-This paper proposes a variation-based online travel time predition approah

More information

Relativistic Addition of Velocities *

Relativistic Addition of Velocities * OpenStax-CNX module: m42540 1 Relativisti Addition of Veloities * OpenStax This work is produed by OpenStax-CNX and liensed under the Creative Commons Attribution Liense 3.0 Abstrat Calulate relativisti

More information

Simplified Modeling, Analysis and Simulation of Permanent Magnet Brushless Direct Current Motors for Sensorless Operation

Simplified Modeling, Analysis and Simulation of Permanent Magnet Brushless Direct Current Motors for Sensorless Operation Amerian Journal of Applied Sienes 9 (7): 1046-1054, 2012 ISSN 1546-9239 2012 Siene Publiations Simplified Modeling, Analysis and Simulation of Permanent Magnet Brushless Diret Current Motors for Sensorless

More information

Finite Element Analysis on Pulsed Laser Forming of Sheet Metal

Finite Element Analysis on Pulsed Laser Forming of Sheet Metal 5 th International & 26 th All India Manufaturing Tehnology, Design and Researh Conferene (AIMTDR 214) Deember 12 th 14 th, 214, IIT Guwahati, Assam, India Finite Element Analysis on Pulsed Laser Forming

More information

Navigation Performance Enhancement using Online Mosaicking

Navigation Performance Enhancement using Online Mosaicking Navigation Performane Enhanement using Online osaiking Vadim Indelman Supervisors: Prof. Pini Gurfil, AE, ehnion Prof. Ehud Rivlin, CS, ehnion Advisor: Dr. Hetor Rotstein, Rafael Nov. Introdution Navigation:

More information

The Geometric Interpretation of Some Mathematical Expressions Containing the Riemann ζ-function

The Geometric Interpretation of Some Mathematical Expressions Containing the Riemann ζ-function Mathematis Letters 6; (6): 4-46 http://www.sienepublishinggroup.om/j/ml doi:.648/j.ml.66. The Geometri Interpretation of Some Mathematial Epressions Containing the Riemann ζ-funtion Yuriy N. Zayko Department

More information

Critical Reflections on the Hafele and Keating Experiment

Critical Reflections on the Hafele and Keating Experiment Critial Refletions on the Hafele and Keating Experiment W.Nawrot In 1971 Hafele and Keating performed their famous experiment whih onfirmed the time dilation predited by SRT by use of marosopi loks. As

More information

Physical Laws, Absolutes, Relative Absolutes and Relativistic Time Phenomena

Physical Laws, Absolutes, Relative Absolutes and Relativistic Time Phenomena Page 1 of 10 Physial Laws, Absolutes, Relative Absolutes and Relativisti Time Phenomena Antonio Ruggeri modexp@iafria.om Sine in the field of knowledge we deal with absolutes, there are absolute laws that

More information

A Differential Equation for Specific Catchment Area

A Differential Equation for Specific Catchment Area Proeedings of Geomorphometry 2009. Zurih, Sitzerland, 3 ugust - 2 September, 2009 Differential Equation for Speifi Cathment rea J. C. Gallant, M. F. Huthinson 2 CSIRO Land and Water, GPO Box 666, Canberra

More information

ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW. P. М. Меdnis

ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW. P. М. Меdnis ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW P. М. Меdnis Novosibirs State Pedagogial University, Chair of the General and Theoretial Physis, Russia, 636, Novosibirs,Viljujsy, 8 e-mail: pmednis@inbox.ru

More information

Exploring the feasibility of on-site earthquake early warning using close-in records of the 2007 Noto Hanto earthquake

Exploring the feasibility of on-site earthquake early warning using close-in records of the 2007 Noto Hanto earthquake Exploring the feasibility of on-site earthquake early warning using lose-in reords of the 2007 Noto Hanto earthquake Yih-Min Wu 1 and Hiroo Kanamori 2 1. Department of Geosienes, National Taiwan University,

More information

Determination of the Aerodynamic Characteristics of Flying Vehicles Using Method Large Eddy Simulation with Software ANSYS

Determination of the Aerodynamic Characteristics of Flying Vehicles Using Method Large Eddy Simulation with Software ANSYS Automation, Control and Intelligent Systems 15; 3(6): 118-13 Published online Deember, 15 (http://www.sienepublishinggroup.om//ais) doi: 1.11648/.ais.1536.14 ISSN: 38-5583 (Print); ISSN: 38-5591 (Online)

More information

Stellar Aberration, Relative Motion, and the Lorentz Factor

Stellar Aberration, Relative Motion, and the Lorentz Factor ong Beah 010 PROCEEDINGS of the NP 1 Stellar berration, Relatie Motion, and the orentz Fator Joseph. Rybzyk 139 Stetson Drie, Chalfont, P 18914-3751 e-mail: jarybzyk@erizon.net Presented are the results

More information

Complexity of Regularization RBF Networks

Complexity of Regularization RBF Networks Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw

More information

Special Theory of Time- Asymmetric Relativity 1 2

Special Theory of Time- Asymmetric Relativity 1 2 Part I Speial Theory of Time- Asymmetri Relatiity 1 The expanding-unierse osmology is founded on the assumption that Einstein s Relatiity is appliable to the entire unierse. This osmology settles diffiulties

More information

Speed Regulation of a Small BLDC Motor using Genetic-Based Proportional Control

Speed Regulation of a Small BLDC Motor using Genetic-Based Proportional Control World Aademy of Siene, Engineering and Tehnology 47 8 Speed Regulation of a Small BLDC Motor using Geneti-Based Proportional Control S. Poonsawat, and T. Kulworawanihpong Abstrat This paper presents the

More information

Identifying Phase Synchronization Clusters in International Business Cycles

Identifying Phase Synchronization Clusters in International Business Cycles Identifying Phase Synhronization Clusters in International Business Cyles Andreas Groth in ooperation with Mihael Ghil, Stephane Hallegatte, and Patrie Dumas Geosienes Department and Environmental Researh

More information

v = fy c u = fx c z c The Pinhole Camera Model Camera Projection Models

v = fy c u = fx c z c The Pinhole Camera Model Camera Projection Models The Pinhole Camera Model Camera Projetion Models We will introdue dierent amera projetion models that relate the loation o an image point to the oordinates o the orresponding 3D points. The projetion models

More information

Breakdown of the Slowly Varying Amplitude Approximation: Generation of Backward Traveling Second Harmonic Light

Breakdown of the Slowly Varying Amplitude Approximation: Generation of Backward Traveling Second Harmonic Light Claremont Colleges Sholarship @ Claremont All HMC Faulty Publiations and Researh HMC Faulty Sholarship 1-1-003 Breakdown of the Slowly Varying Amplitude Approximation: Generation of Bakward Traveling Seond

More information

Gyrokinetic calculations of the neoclassical radial electric field in stellarator plasmas

Gyrokinetic calculations of the neoclassical radial electric field in stellarator plasmas PHYSICS OF PLASMAS VOLUME 8, NUMBER 6 JUNE 2001 Gyrokineti alulations of the neolassial radial eletri field in stellarator plasmas J. L. V. Lewandowski Plasma Physis Laboratory, Prineton University, P.O.

More information

Collinear Equilibrium Points in the Relativistic R3BP when the Bigger Primary is a Triaxial Rigid Body Nakone Bello 1,a and Aminu Abubakar Hussain 2,b

Collinear Equilibrium Points in the Relativistic R3BP when the Bigger Primary is a Triaxial Rigid Body Nakone Bello 1,a and Aminu Abubakar Hussain 2,b International Frontier Siene Letters Submitted: 6-- ISSN: 9-8, Vol., pp -6 Aepted: -- doi:.8/www.sipress.om/ifsl.. Online: --8 SiPress Ltd., Switzerland Collinear Equilibrium Points in the Relativisti

More information

Electromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution.

Electromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution. arxiv:physis/99536v1 [physis.lass-ph] 15 May 1999 Eletromagneti radiation of the travelling spin wave propagating in an antiferromagneti plate. Exat solution. A.A.Zhmudsky November 19, 16 Abstrat The exat

More information

The experimental plan of displacement- and frequency-noise free laser interferometer

The experimental plan of displacement- and frequency-noise free laser interferometer 7th Edoardo Amaldi Conferene on Gravitational Waves (Amaldi7) Journal of Physis: Conferene Series 122 (2008) 012022 The experimental plan of displaement- and frequeny-noise free laser interferometer K

More information

ON-LINE ESTIMATION OF THE VENTILATION RATE OF GREENHOUSES

ON-LINE ESTIMATION OF THE VENTILATION RATE OF GREENHOUSES ON-LINE ESTIMATION OF THE VENTILATION RATE OF GREENHOUSES J. Bontsema *, E.J. van Henten *, J.G. Kornet *, J. Budding ** and T. Rieswijk ** * Agrotehnology and Food Innovations B.V., Greenhouse Tehnology

More information

Formal Specification for Transportation Cyber Physical Systems

Formal Specification for Transportation Cyber Physical Systems Formal Speifiation for Transportation Cyber Physial Systems ihen Zhang, Jifeng He and Wensheng Yu Shanghai Key aboratory of Trustworthy Computing East China Normal University Shanghai 200062, China Zhanglihen1962@163.om

More information

REGULATION AND INPUT DISTURBANCE SUPPRESSION FOR PORT-CONTROLLED HAMILTONIAN SYSTEMS. Luca Gentili Arjan van der Schaft

REGULATION AND INPUT DISTURBANCE SUPPRESSION FOR PORT-CONTROLLED HAMILTONIAN SYSTEMS. Luca Gentili Arjan van der Schaft REGULATION AND INPUT DISTURBANE SUPPRESSION FOR PORTONTROLLED HAMILTONIAN SYSTEMS Lua Gentili Arjan van der Shaft ASYDEIS University of Bologna viale Risorgimento 4136 Bologna Italy email: lgentili@deisuniboit

More information

Relativistic Analysis of Doppler Effect and Aberration based on Vectorial Lorentz Transformations

Relativistic Analysis of Doppler Effect and Aberration based on Vectorial Lorentz Transformations Uniersidad Central de Venezuela From the SeletedWorks of Jorge A Frano June, Relatiisti Analysis of Doppler Effet and Aberration based on Vetorial Lorentz Transformations Jorge A Frano, Uniersidad Central

More information

Astr 5465 Mar. 29, 2018 Galactic Dynamics I: Disks

Astr 5465 Mar. 29, 2018 Galactic Dynamics I: Disks Galati Dynamis Overview Astr 5465 Mar. 29, 2018 Subjet is omplex but we will hit the highlights Our goal is to develop an appreiation of the subjet whih we an use to interpret observational data See Binney

More information

Hankel Optimal Model Order Reduction 1

Hankel Optimal Model Order Reduction 1 Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both

More information

Singular Event Detection

Singular Event Detection Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate

More information

MICROSCALE SIMULATIONS OF CONDUCTIVE / RADIATIVE HEAT TRANSFERS IN POROUS MEDIA

MICROSCALE SIMULATIONS OF CONDUCTIVE / RADIATIVE HEAT TRANSFERS IN POROUS MEDIA MICROCALE IMULAION OF CONDUCIVE / RADIAIVE HEA RANFER IN POROU MEDIA J.-F. hovert, V.V. Mourzenko, C. Roudani Institut PPRIME-CNR Context, motivation moldering in porous media 400K (measured) Mirosale

More information

GEOMETRY FOR 3D COMPUTER VISION

GEOMETRY FOR 3D COMPUTER VISION GEOMETRY FOR 3D COMPUTER VISION 1/29 Válav Hlaváč Czeh Tehnial University, Faulty of Eletrial Engineering Department of Cybernetis, Center for Mahine Pereption 121 35 Praha 2, Karlovo nám. 13, Czeh Republi

More information

ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES.

ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES. ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES. All systems with interation of some type have normal modes. One may desribe them as solutions in absene of soures; they are exitations of the system

More information

19 Lecture 19: Cosmic Microwave Background Radiation

19 Lecture 19: Cosmic Microwave Background Radiation PHYS 652: Astrophysis 97 19 Leture 19: Cosmi Mirowave Bakground Radiation Observe the void its emptiness emits a pure light. Chuang-tzu The Big Piture: Today we are disussing the osmi mirowave bakground

More information

Sparsity based Ground Moving Target Imaging via Multi-Channel SAR

Sparsity based Ground Moving Target Imaging via Multi-Channel SAR Sparsity based Ground Moving Target Imaging via Multi-Channel SAR Di Wu, Mehrdad Yaghoobi and Mike Davies Shool of Engineering University of Edinburgh UK, EH9 3JL Email: {D.Wu, m.yaghoobi-vaighan, mike.davies}@ed.a.uk

More information

Determining Aerodynamic Loads Based on Optical Deformation Measurements

Determining Aerodynamic Loads Based on Optical Deformation Measurements Determining Aerodynami Loads Based on Optial Deformation Measurements Tianshu Liu, Western Mihigan University, Kalamazoo, MI D. A. Barrows, A. W. Burner & R. D. Rhew NASA Langley Researh Center, Hampton,

More information

Breakdown of the Special Theory of Relativity as Proven by Synchronization of Clocks

Breakdown of the Special Theory of Relativity as Proven by Synchronization of Clocks Breakdown of the Speial Theory of Relativity as Proven by Synhronization of Cloks Koshun Suto Koshun_suto19@mbr.nifty.om Abstrat In this paper, a hypothetial preferred frame of referene is presumed, and

More information

Entropic Graphs for Manifold Learning

Entropic Graphs for Manifold Learning Entropi Graphs for Manifold Learning Jose A. Costa and Alfred O. Hero III Department of Eletrial Engineering and Computer Siene Uniersity of Mihigan, Ann Arbor, MI 9 Email: josta@umih.edu, hero@ees.umih.edu

More information

MOVING OBJECTS OBSERVATION THEORY IN PLACE OF SPECIAL RELATIVITY

MOVING OBJECTS OBSERVATION THEORY IN PLACE OF SPECIAL RELATIVITY Inquiry, ol. 8, no., Deember 007, pp. 4 49 IIGSS Aademi Publisher MOVING OBJECTS OBSERVATION THEORY IN PLACE OF SPECIAL RELATIVITY LI ZIFENG Petroleum Engineering Institute, Yanshan Uniersity, Qinhuangdao,

More information

TAP 702-6: Binary stars

TAP 702-6: Binary stars TAP 702-6: Binary stars Orbiting binary stars: A type of ariable star. This type of ariable star onsists of two stars orbiting around eah other. When the dier star is in front of the brighter one, the

More information