DARWIN-OP HUMANOID ROBOT KINEMATICS. Robert L. Williams II, Ph.D., Mechanical Engineering, Ohio University, Athens, Ohio, USA

Size: px
Start display at page:

Download "DARWIN-OP HUMANOID ROBOT KINEMATICS. Robert L. Williams II, Ph.D., Mechanical Engineering, Ohio University, Athens, Ohio, USA"

Transcription

1 Proceedngs of the ASME 1 Internatonal Desgn Engneerng Techncal onferences & omputers and Informaton n Engneerng onference IDET/IE 1 August 1-15, 1, hcago, IL, USA DET1-765 DARWIN-OP HUMANOID ROBOT KINEMATIS Robert L. Wllams II, Ph.D., wllar4@oho.edu Mechancal Engneerng, Oho Unversty, Athens, Oho, USA ABSTRAT Ths paper presents knematcs analyss for the DARwIn- OP (Dynamc Anthropomorphc Robot wth Intellgence Open-Platform) robot. Ths s a -dof humanod walkng robot developed by Vrgna Tech, Purdue, and the Unversty of Pennsylvana and marketed by Robots Inc. The robot verson analyzed n ths paper s 455 mm tall and has a mass of.8 kg. Oho Unversty has two of these unts for robot applcatons research and teachng. Presented are a descrpton of DARwIn-OP, the Denavt- Hartenberg Parameters for each seral chan (-dof head pan/tlt, 3-dof arms, and 6-dof legs), specfc length parameters, jont angle lmts, plus forward pose knematcs equatons and partal nverse pose knematcs solutons, wth examples. 1. INTRODUTION DARwIn-OP s a unque open-platform humanod robot sponsored by NSF, developed collaboratvely by Vrgna Tech ( Purdue, and the Unversty of Pennsylvana, and marketed by Robots Inc. ( Thanks to NSF support, Oho Unversty was one of the nsttutons who receved two DARwIn humanod robots for research and teachng. DARwIn s beng used to play robot soccer, but we are nvestgatng other humanod robot applcatons. DARwIn-OP s ntended to be open-platform, wth a communty of users contrbutng to robot developments followng the ntal collaborator s developments. The software and control archtecture s all open-source and free to the world. Multple computer platforms and varous programmng languages are sad to be allowed. Programmng enhancements should be shared wth the DARwIn communty. In addton, the desgn of the humanod robot s open-source, wth complete AD fles and parts lsts avalable for free so one could buld the robot from scratch. One could make mprovements and nnovatons to the exstng robot desgn, whch should also be shared wth the DARwIn communty. Ths paper s organzed as follows. Frst we descrbe the degrees-of-freedom (dof) DARwIn-OP robot concept and hardware. Then we present the Denavt-Hartenberg (DH) parameters for each seral chan of the humanod robot (-dof head pan/tlt, 3-dof arms, and 6-dof legs). Included n ths secton are the specfc length parameters and jont angle lmts for the 455 mm tall DARwIn-OP. Next we present forward pose knematcs equatons for each seral chan, wth examples. Last we present partal nverse pose knematcs solutons for each seral chan of the robot. All analyses n ths paper are from well-known seral robotcs technques (rag, 5). We could not fnd ths knematcs analyss from any source so ths wll be part of our contrbuton to the DARwIn communty.. DARwIn-OP DESRIPTION The -dof DARwIn Humanod Robot s shown n Fgure 1, n the standard pose wth all zero jont angles. In the hardware shown n Fgure 1, each arm has three sngle-dof revolute (R) jonts, each leg has sx sngle-dof R jonts, and the pan/tlt head has two sngle-dof R jonts. Therefore, the overall robot has dof. Each arm has a -dof (offset-u-jont) shoulder jont and a 1-dof elbow jont, for a total of 3-dof per arm. Each leg has a 3-dof (S-jont wth 3 ntersectng R jonts) hp jont, a 1-dof knee jont, and a -dof (U-jont) ankle jont for a total of 6-dof per leg. The last hp jont, the knee jont, and the frst ankle jont are all about parallel Z axes. The -dof head (U-jont) enables pan and tlt for the camera va an azmuth R jont and an elevaton R jont. Each DARwIn-OP robot jont s drven by a Dynamxel MX-8 servomotor controlled by an nternal M-73 Robots servo-controller. The DARwIn-OP sensors nclude an HD dgtal vdeo camera, tr-axs gyroscope and accelerometer, and two mcrophones. The DARwIn-OP on-board computer s a 1.6 GHz Intel Atom Z53 FtP wth 4 GB SSD and the Ubuntu Lnux operatng system. DARwIn s programmed n ++ and 1 opyrght 1 by ASME Downloaded From: on /18/16 Terms of Use:

2 s able to operate both tethered and battery-powered. The LPO 3ELL 11.1v 1mAh battery provdes up to 3 operaton mnutes (Merln Robotcs, 1). Secton 3.3. Refer to the lengths defntons n Fgure 3 for these parameters. The basc torso reference frame s shown n red n Fgure 3. Ths {} dextral artesan reference frame, stands for hest. Frame {} s movng and serves as the reference frame for all fve seral chans. A fxed world {Wo} artesan reference frame (not shown) may be used for DARwIn-OP absolute knematc referencng. Fgure 1. DARwIn-OP Humanod Robot, Zero Pose The numberng conventon for the jont servomotors s shown n Fgure. The jont angle values wll be denoted by, = 1,,,, accordng to the numberng scheme n Fgure. Agan, the pose shown s for all zero jont angles. For notatonal smplcty n artesan coordnate frame defnton, frame numbers {}, {1}, {}, etc. wll be recycled for the fve seral chans. Thus, there wll be fve Z axes, four Z 3 axes, two Z 6 axes, and so on. Ths s to avod the onerous noton of Z RA, Z LL, etc. (representng rght arm and left leg). Accordng to rag (5) notaton for the Denavt- Hartenberg Parameters (1955), jont angle s supposed to rotate about R-jont axs Z. Agan for smplcty n notaton, ths rule wll be volated, n order to number the jont angles as shown n Fgure. 3. DARwIn-OP PARAMETERS The DARwIn-OP Denavt-Hartenberg (DH, Denavt and Hartenberg, 1955) parameters are presented n ths secton, for each of the seral chans (head, arm, and leg). DH parameters have been adopted for standard knematcs analyss n seralchan robots (rag, 5). The physcal length parameters are gven frst, followed by the jont angle lmts, and then the DH parameters for the DARwIn-OP humanod robot. We present DH parameters for the -dof head, the 3-dof rght arm, and the 6-dof rght leg (the 3-dof left arm and 6-dof left leg are symmetrc to ther counterparts, as seen n Fgures 1 and ). 3.1 DARwIn-OP robot lengths Ths secton presents the varous lnk lengths for the DARwIn-OP Humanod Robot. These are specfc hardware desgn values assocated wth the DH Parameters gven n Fgure. DARwIn-OP Servomotor Numberng (Robots, 11) Torso Lengths These values are not DH Parameters but wll be ntroduced later n the Forward Pose Knematcs solutons. The subscrpt h for L h s lower-case and ndcates head. Length Value (mm) L h 5.5 L 8. L TX 5. L TY 1. L TZ Head Lengths These values are not DH Parameters but wll be ntroduced later n the Forward Pose Knematcs solutons. Length Value (mm) H X 33. H Y 34.4 H Z.5 opyrght 1 by ASME Downloaded From: on /18/16 Terms of Use:

3 Length Value (mm) F X 14. F Y 15. (1. front) F Z 66. L X 5. L Z DARwIn-OP Jont Angle Lmts Ths secton presents the knematc jont angle lmts for the DARwIn-OP Humanod Robot. These are specfc hardware desgn values assocated wth the DH Parameters gven n the followng secton. Refer to the jont numberng gven n Fgure. All unts n Table I are degrees. Fgure 3. DARwIn Robot Lengths Defntons Rght and Left Arm Lengths These DH Parameters are shown n Fgure 3. L H s not a DH Parameter but wll be used later n the Forward Pose Knematcs solutons. The subscrpt H for L H s upper-case and ndcates hand. Length Value (mm) L L 6. L L H Rght and Left Leg Lengths These DH Parameters are shown n Fgure 3. L F s not a DH Parameter but wll be used later n the Forward Pose Knematcs solutons. Length Value (mm) L L L F Rght and Left Foot Lengths DARwIn s feet (see Fgure 3) are 14 mm long, 66 mm wde and 15 mm thck (taperng to 1 mm thck at the front edge). The feet are centered wth the ankle jont front-to-back (L X = 5 mm), but not sde-to-sde (L Z = 3 mm from the nner edge of each foot to the ankle jont). Table I. DARwIn-OP Jont Lmts Jont Jont name Axs MIN MAX Rght Arm 1 shoulder ptch Z shoulder roll Z elbow Z 3 16 Left Arm shoulder ptch Z shoulder roll Z elbow Z 3 16 Rght Leg 7 hp yaw Z hp roll Z 6 9 hp ptch Z knee Z ankle ptch Z ankle roll Z Left Leg 8 hp yaw Z hp roll Z 6 1 hp ptch Z knee Z ankle ptch Z ankle roll Z Head 19 head pan Z head tlt Z DARwIn-OP DH Parameters The standardzed Denavt-Hartenberg (DH) Parameters (1955) are used to descrbe the lnks/jonts geometry of a seralchan robot. In ths secton the DH Parameters are presented for each of the DARwIn-OP seral chans,.e. pan/tlt head, rght arm, and rght legs. rag (5) conventon s used for the DH Parameters n ths paper. All angular unts n the DH Parameter tables below are degrees. 3 opyrght 1 by ASME Downloaded From: on /18/16 Terms of Use:

4 The artesan reference frame defntons for three DARwIn-OP seral chans (-dof head, 3-dof rght arm, 6-dof rght leg) are gven n Fgures 4, 5, and 6. It s mportant to note that all of these DH parameters fgures are shown n the zero-angle poses. Therefore t looks lke many movng coordnate frame axes lne up, but ths wll change as jont angles rotate away from zero n all cases. The 3-dof left arm and 6-dof left leg are symmetrc to the rght arm and leg, so these cases are not presented here Two-dof Pan/Tlt Head The artesan reference frame defntons for the two-dof pan/tlt head are shown n Fgure 4. Table II gves the DH parameters for the two-dof pan/tlt head. Table III. Three-dof Rght Arm DH Parameters 1 a 1 d L 1 L L Sx-dof Rght Leg The artesan reference frame defntons for the sx-dof rght leg are shown n Fgure 6. Table IV gves the DH parameters for the sx-dof rght leg. Fgure 4. Two-dof Pan/Tlt Head oordnate Frames Table II. Two-dof Pan/Tlt Head DH Parameters 1 a 1 d Three-dof Rght Arm The artesan reference frame defntons for the three-dof rght arm are shown n Fgure 5. Table III gves the DH parameters for the three-dof rght arm. Fgure 6. Sx-dof Rght Leg oordnate Frames Fgure 5. Three-dof Rght Arm oordnate Frames (Top and Front Vews) Table IV. Sx-dof Rght Leg DH Parameters 1 a 1 d L L opyrght 1 by ASME Downloaded From: on /18/16 Terms of Use:

5 4. DARwIn-OP FORWARD POSE KINEMATIS In general, the Forward Pose Knematcs (FPK) problem for a seral-chan robot s stated: Gven the jont values, calculate the pose (poston and orentaton) of the end frame of nterest. For seral-chan robots, the FPK problem set up and soluton s straght-forward. It s based on substtutng each lne of the Denavt-Hartenberg Parameters table nto the equaton below (rag, 5), gvng the pose of frame {} wth respect to ts nearest neghbor frame { 1} back along the seral chan: c s a 1 1 sc 1 cc 1 s 1 ds 1 T ss 1 cs 1 c 1 dc 1 1 The equaton above represents pose (poston and orentaton) of frame {} wth respect to frame { 1} by usng a 4x4 homogeneous transformaton matrx. The upper left 3x3 matrx s the rotaton matrx gvng the orentaton and the upper rght 3x1 vector s the poston vector. Then transformaton equatons are used to fnd the pose of the overall end-frame of nterest wth respect to the base reference frame, to complete the FPK soluton for each seral chan. 4.1 Two-dof Pan/Tlt Head FPK Expressons The statement of the FPK problem for the two-dof pan/tlt head seral chan of the DARwIn-OP humanod robot s: Gven, calculate 19, T and X. where X s the pontng vector along the X unt vector drecton, expressed n {} coordnates. Substtute each row of the DH parameters n Table II nto 1 the above equaton for T to obtan the two neghborng homogeneous transformaton matrces as a functon of the jont angles. c19 s19 s19 c19 1T 1 1 c s 1 1 T s c 1 Where the followng abbrevatons were used: c cos, s sn, for = 19,. Now substtute these two neghborng homogeneous transformaton matrces nto the followng homogeneous transform equaton to derve the FPK result. c19c c19s s19 1 s19c s19s c19 T 1T19 T s c 1 Note that the fourth column result, the poston vector X gvng the orgn of {} from the orgn of {}, expressed n {} coordnates, s all zero components snce the pan/tlt axes are ntersectng. There are no knematc lengths or length offsets nvolved. To derve the second FPK result X, use the followng transform equaton, where X 1 1 T. c19c s19c X T X s 1 The desred pontng vector s contaned n the frst three elements of X,.e. the 1 n the fourth locaton s merely a placeholder for homogeneous transformaton matrx multplcaton. The basc pan/tlt head FPK result s T. To calculate the pose of each eye frame (RE and LE for rght and left eye, respectvely) wth respect to the DARwIn chest reference frame {}, the followng transform equatons are used. RET T T 19, RET LET T T 19, LET 1 1 H X 1 L h T 1 H Y 1 RET 1 H Z H X 1 H Y LET 1 H Z 1 The overall transform equatons above can be evaluated numercally. onstants L h, H X, H Y, and H Z were gven earler. Note that T, T RE, and T LE are not evaluated by any row n the DH parameter table, snce there s no varable assocated wth these fxed homogeneous transformaton matrces based on constant lengths and orentaton. Instead, they are determned by nspecton, usng the rotaton matrx and 5 opyrght 1 by ASME Downloaded From: on /18/16 Terms of Use:

6 poston vector components of the homogeneous transformaton matrx defnton. In the above expressons for T RE and T LE the orentaton of each eye frame {RE} and {LE} was arbtrarly assgned as dentty,.e. dentcal to the orentaton of {}. Ths can easly be changed as requred. 4. Three-dof Rght Arm FPK Expressons The statement of the FPK problem for the three-dof rght arm seral chan of the DARwIn-OP humanod robot s: Gven 1, 3, 5, calculate T 3 and T H. Where {H} s the rght-arm end-effector (hand) frame and {} s the DARwIn reference frame n the chest. As stated earler, for notatonal smplcty n artesan coordnate frame defnton, frame numbers {}, {1}, {}, etc. wll be recycled for the fve seral chans (head, rght and left arms, rght and left legs). Therefore, n ths paper many repeated frame numbers are context-dependent, whch must be sorted out n programmng the overall humanod robot. Substtute each row of the DH parameters n Table III nto the equaton for 1 T to obtan the three neghborng homogeneous transformaton matrces as a functon of the jont angles. c1 s1 s3 c3 L1 s1 c1 1T 1 1 L 3 1 T c3 s3 1 1 c5 s5 L 1 3T s5 c5 1 Where the followng abbrevatons were used: c cos, s sn, for = 1,3,5. Now substtute these three neghborng homogeneous transformaton matrces nto the followng homogeneous transform equaton to derve the FPK result. 1 3T 1T 1 T 3 3T 5 s1s5csc sc 1 5css cc 1 3 Lc 1 1Lcs 1 3Ls 3 1 cs 1 5ssc cc 1 5sss sc 1 3 Ls 1 1Lss 1 3Lc 3 1 cc 3 5 cs 3 5 s3 Lc 3 1 T H The basc rght-arm FPK result s 3 T. To calculate, the pose of the rght-arm end-effector frame {H} wth respect to the DARwIn chest reference frame {}, the followng transform equaton s used. 3 HT T 3T 1, 3, 5 HT 1 1 LH 1 T 3 1 L 3 1 L HT The overall transform equaton above can be evaluated numercally. onstants L, L 3, and L H were gven earler. Note that T and 3 T H are not evaluated by any row n the DH parameter table, snce there s no varable assocated wth these fxed homogeneous transformaton matrces based on constant lengths and orentaton. Instead, they are determned by nspecton, usng the rotaton matrx and poston vector components of the homogeneous transformaton matrx defnton. For use n nverse pose knematcs, here s the resultng poston vector of the hand tp wth respect to the {} frame: X 3 H 3T 1, 3, 5 XH Lc 1 1Lcs 1 3L3( s1sc 1 5cs 1 3s5) LH ( ss 1 5csc 1 3 5) Ls 1 1 Lss 1 3 L3( c1 cc 1 5 sss 1 3 5) LH ( cs 1 5 ssc 1 3 5) Lc 3Lcs 3 3 5LH cc Sx-dof Rght Leg FPK Expressons The statement of the FPK problem for the sx-dof rght leg seral chan of the DARwIn-OP humanod robot s: Gven 7, 11, 9, 13, 17, 15, calculate T 6 and T F. Where {F} s the rght-leg end-effector (foot) frame and {} s the DARwIn chest reference frame. As stated earler, for notatonal smplcty n artesan coordnate frame defnton, frame numbers {}, {1}, {}, etc. wll be recycled for the fve seral chans (head, rght and left arms, rght and left legs). Therefore, n ths paper many repeated frame numbers are context-dependent, whch must be sorted out n programmng the overall humanod robot. Substtute each row of the DH parameters n Table IV nto 1 the equaton for T to obtan the sx neghborng homogeneous transformaton matrces as a functon of the jont angles. c7 s7 s11 c11 s7 c7 1T T c11 s opyrght 1 by ASME Downloaded From: on /18/16 Terms of Use:

7 c9 s9 c13 s13 L4 1 3T 3 s13 c13 4 s9 c9 T c17 s17 L5 c15 s15 4 s17 c17 5T T s15 c Where the followng abbrevatons were used: c cos, s sn, for = 7,11,9,13,17,15. Now substtute these three neghborng homogeneous transformaton matrces nto the followng homogeneous transform equaton to derve the FPK result T 1T 7 T 11 3T 9 4T 13 5T 17 6T15 Snce there are three parallel Z axes (3,4,5) n the rght leg, we can group the above matrx multplcatons as follows for a sgnfcant smplfcaton T 1T 7 T 11 5T 9, 13, 17 6T15 where: cabc sabc L4c9L5cab 1 5T sabc cabc L4s9 L5s ab 1, abc abc were used. The rght leg FPK result s: r11 r1 r13 px r1 r r3 p y 6T r31 r3 r33 p z 1 where the abbrevatons cab cos9 13, sab sn 9 13 c cos, s sn where: r11 s7sabc c7s11cabc c15 c7c11s15 r c s s s c c s c s r s s c c c 1 7 abc 7 11 abc abc r1 s7sabc c7s11cabc s15 c7c11c15 r c s s s c s s c c r s c c s c 7 abc 7 11 abc abc r s c c s s r c c s s s r c s 13 7 abc 7 11 abc 3 7 abc 7 11 abc abc px L4s9L5sab s7 L4c9L5cab c7s11 p L s L s c L c L c s s p L c L c c y ab ab 7 11 z ab 11 The basc rght-leg FPK result s 6 T. To calculate T F, the pose of the rght-leg end-effector frame {F} wth respect to the DARwIn chest reference frame {}, the followng transform equaton s used. 6 FT T 6T 7, 11, 9, 13, 17, 15 FT 1 LTX 1 L TY T 1 L TZ 1 1 LF 6 1 FT 1 1 The overall transform equaton above can be evaluated numercally. onstants L TX, L TY, L TZ, and L F were gven earler. Note that T and 6 T F are not evaluated by any row n the DH parameter table, snce there s no varable assocated wth these fxed homogeneous transformaton matrces based on constant lengths and orentaton. Instead, they are determned by nspecton, usng the rotaton matrx and poston vector components of the basc homogeneous transformaton matrx defnton. 4.4 FPK Examples Two-dof Pan/Tlt Head, 1, : Gven T Three-dof Rght Arm,, 1,,3 Gven T H : X opyrght 1 by ASME Downloaded From: on /18/16 Terms of Use:

8 4.4.3 Sx-dof Rght Leg,,,,, 1, 3,, 35, 45,6 Gven T F These FPK results were valdated by comparng the symbolc formulae derved n ths paper wth a numercal approach n each case. Fgure 7 shows the MATLAB graphcal results for the zero confguraton (see Fgures 1 and ). Fgure 8 shows the MATLAB graphcal results for the FPK examples confguraton (symmetry s lost snce we are usng dentcal, rather than symmetrc, jont angles for the rght and left arms and legs). Fgure 8. DARwIn-OP MATLAB Smulaton, examples confguraton 5. DARwIn-OP INVERSE POSE KINEMATIS In general, the Inverse Pose Knematcs (IPK) problem for a seral-chan robot s stated: Gven the pose (poston and orentaton) of the end frame of nterest, calculate the jont values. For seral-chan robots, the IPK soluton starts wth the FPK equatons. The soluton of coupled nonlnear algebrac equatons s requred and multple soluton sets generally result. The head IPK soluton s complete but the rght arm and rght leg IPK solutons are partal below. The FPK examples can also serve as IPK examples when the nput s reversed. Fgure 7. DARwIn-OP MATLAB Smulaton, zero confguraton 5.1 Two-dof Pan/Tlt Head IPK Soluton The statement of the IPK problem for the two-dof pan/tlt head seral chan of the DARwIn-OP humanod robot s: Gven X, calculate 19,. X s the pontng vector along the X unt vector drecton, expressed n {} coordnates; t must be entered as a unt vector. Its FPK expresson s recalled below: c19c x X s19c y s z where x, y, z are the gven components of the desred unt vector X. 8 opyrght 1 by ASME Downloaded From: on /18/16 Terms of Use:

9 The IPK soluton for the pan/tlt head s obtaned n the followng order. Frst use a rato of the y to x equatons to elmnate c : 19 atan( yx, ) where atan s the quadrant-specfc nverse tangent functon. Next use a rato of the z to x equatons: atan( zx, / c ) 19 or, equvalently, use a rato of the z to y equatons: atan( zy, / s ) 19 If 19 9 use the y equaton soluton, and f 19,18 use the x equaton soluton. If nether of these artfcal sngularty condtons exst, both solutons wll yeld the same result for. There s a sngle soluton set,.e. no multple solutons n ths case. 5. Three-dof Rght Arm IPK Soluton (partal) The statement of the IPK problem for the three-dof rght arm seral chan of the DARwIn-OP humanod robot s: Gven,,. T H, calculate Where {H} s the rght-arm end-effector hand frame and {} s the DARwIn chest reference frame. The gven pose T H represents 6-dof, whle the unknown jont space n ths case s only 3-dof. Therefore, the entre T H cannot be gven, but only a 3-dof subspace. Let us specfy the poston vector and calculate the unknown jont angles: Gven X H, calculate 1, 3, 5. From FPK, here are the equatons to solve: L1c1 Lcs 1 3L3( s1sc 1 5css 1 3 5) LH ( ss 1 5csc 1 3 5) x L1s1 Lss 1 3 L3( c1 cc 1 5 sss 1 3 5) LH ( cs 1 5 ssc 1 3 5) XH y Lc 3 Lcs 3 3 5LHcc 3 5 z 1 where x, y, z are the gven components of the desred tp poston vector X H. We must frst smplfy the gven X H to X H usng the followng transform equaton: 1 X H T XH Where constant transform matrx T was gven n the FPK secton for the rght arm. The above set of three equatons fully coupled n the three unknowns,, s nonlnear (transcendental). In order to smplfy ths set of equatons, use the followng transform equaton: x 1 3 X H y 1T1 T3 3T5 XH z x T T y T X z H xc s yss zc Ls L Ls L c H 5 xc c ys c zs L c xs yc L L L c L s H 5 Ths soluton has not yet been completed. 5.3 Sx-dof Rght Leg IPK Soluton (partal) The statement of the IPK problem for the sx-dof rght leg seral chan of the DARwIn-OP humanod robot s: Gven,,,,,. T F, calculate Where {F} s the rght-leg end-effector foot frame and {} s the DARwIn chest reference frame. We can frst smplfy the gven T F to T 6 usng the followng transform equaton: T T FT FT Where constant transform matrces T and 6 T F were gven n the FPK secton for the rght leg. Let the followng symbols represent the numercal values for T 6, derved from the orgnally-gven T F. r11 r1 r13 x r1 r r3 y 6T r31 r3 r33 z 1 These numercal values are equated to the FPK soluton as a functon of the sx unknown jont angles: T 1T 7 T 11 3T 9 4T 13 5T 17 6T15 The above represents 16 equatons (4 trval) n the sx unknown jont angles. The three poston vector component equatons are all ndependent, but there are only three ndependent equatons among the 9 equatons of the rotaton matrx. The equatons are coupled and nonlnear (transcendental). 9 opyrght 1 by ASME Downloaded From: on /18/16 Terms of Use:

10 Lookng at the FPK terms gven earler, we need some smplfcaton pror to solvng these equatons. There are three possbltes: 1. Notce that the hp s sphercal (.e. frames {}, {1}, {}, {3} share a common orgn pont);. FPK already took advantage of the three parallel Z axes (last hp Z 9, knee Z 13, and frst ankle Z 17 ); 3. Snce the hp s sphercal and the - dof ankle also rotates about a common pont, the dstance from the hp center to the ankle center s only a functon of one jont angle, 13. The thrd smplfcaton wll be exploted frst. 6 5 P P x y z 3 3 P L L c L s P L L L L c x y z cos 13 1 P L L LL Note ths s equvalent to the Law of osnes, lookng at the plane contanng the hp, knee, ankle, and both leg lengths L 4 and L 5. Ths soluton has not yet been completed. 6. ONLUSION Ths paper has presented knematc analyss for the -dof DARwIn-OP humanod walkng robot, ncludng the Denavt- Hartenberg Parameters for each seral chan (-dof head pan/tlt, 3-dof arm, and 6-dof leg), specfc length parameters, jont angle lmts, plus forward pose knematcs equatons and partal nverse pose knematcs solutons, wth examples. Future work plans nclude fnshng the IPK solutons. REFERENES J.J. rag, 5, Introducton to Robotcs: Mechancs and ontrol, Thrd Edton, Pearson Prentce Hall, Upper Saddle Rver, NJ. J. Denavt and R.S. Hartenberg, 1955, A Knematc Notaton for Lower-Par Mechansms Based on Matrces, Journal of Appled Mechancs: Merln Robotcs, 1, ndex.php/bpedal-robots/darwn-op-1.html Robots, 11, DARwIn-OP Quck Start Operaton and Programmng Gude, 1 opyrght 1 by ASME Downloaded From: on /18/16 Terms of Use:

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

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

More information

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

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 4. Moton Knematcs 4.2 Angular Velocty Knematcs Summary From the last lecture we concluded that: If the jonts

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

MEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut

MEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut MEV442 Introducton to Robotcs Module 2 Dr. Santhakumar Mohan Assstant Professor Mechancal Engneerng Natonal Insttute of Technology Calcut Jacobans: Veloctes and statc forces Introducton Notaton for tme-varyng

More information

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

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

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

Iterative General Dynamic Model for Serial-Link Manipulators

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

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

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

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

More information

Indeterminate pin-jointed frames (trusses)

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

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

Mathematical Preparations

Mathematical Preparations 1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system. Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and

More information

Programming Project 1: Molecular Geometry and Rotational Constants

Programming Project 1: Molecular Geometry and Rotational Constants Programmng Project 1: Molecular Geometry and Rotatonal Constants Center for Computatonal Chemstry Unversty of Georga Athens, Georga 30602 Summer 2012 1 Introducton Ths programmng project s desgned to provde

More information

Robot Modeling and Kinematics Errata R. Manseur. List of errors and typos reported as of 4/1/2007:

Robot Modeling and Kinematics Errata R. Manseur. List of errors and typos reported as of 4/1/2007: Robot Modelng and Knematcs Errata R. Manseur. Lst of errors and typos reported as of 4//7:. Fle CF4.wrl mentoned on page 45 s not on the CD. Fle ZYZ_Sm.wrl mentoned on page 67 s not on the CD. 3. Pg. 43.

More information

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars

More information

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

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

More information

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016 ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

APPENDIX A Some Linear Algebra

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

More information

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

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

More information

12. The Hamilton-Jacobi Equation Michael Fowler

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

More information

Spin-rotation coupling of the angularly accelerated rigid body

Spin-rotation coupling of the angularly accelerated rigid body Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s

More information

PHYS 705: Classical Mechanics. Canonical Transformation II

PHYS 705: Classical Mechanics. Canonical Transformation II 1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

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

More information

Department of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution

Department of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution Department of Statstcs Unversty of Toronto STA35HS / HS Desgn and Analyss of Experments Term Test - Wnter - Soluton February, Last Name: Frst Name: Student Number: Instructons: Tme: hours. Ads: a non-programmable

More information

MEM 255 Introduction to Control Systems Review: Basics of Linear Algebra

MEM 255 Introduction to Control Systems Review: Basics of Linear Algebra MEM 255 Introducton to Control Systems Revew: Bascs of Lnear Algebra Harry G. Kwatny Department of Mechancal Engneerng & Mechancs Drexel Unversty Outlne Vectors Matrces MATLAB Advanced Topcs Vectors A

More information

THE SUMMATION NOTATION Ʃ

THE SUMMATION NOTATION Ʃ Sngle Subscrpt otaton THE SUMMATIO OTATIO Ʃ Most of the calculatons we perform n statstcs are repettve operatons on lsts of numbers. For example, we compute the sum of a set of numbers, or the sum of the

More information

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to ChE Lecture Notes - D. Keer, 5/9/98 Lecture 6,7,8 - Rootndng n systems o equatons (A) Theory (B) Problems (C) MATLAB Applcatons Tet: Supplementary notes rom Instructor 6. Why s t mportant to be able to

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω

More information

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department

More information

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

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

More information

Linear Regression Analysis: Terminology and Notation

Linear Regression Analysis: Terminology and Notation ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented

More information

An Interactive Optimisation Tool for Allocation Problems

An Interactive Optimisation Tool for Allocation Problems An Interactve Optmsaton ool for Allocaton Problems Fredr Bonäs, Joam Westerlund and apo Westerlund Process Desgn Laboratory, Faculty of echnology, Åbo Aadem Unversty, uru 20500, Fnland hs paper presents

More information

Linear Momentum. Center of Mass.

Linear Momentum. Center of Mass. Lecture 6 Chapter 9 Physcs I 03.3.04 Lnear omentum. Center of ass. Course webste: http://faculty.uml.edu/ndry_danylov/teachng/physcsi Lecture Capture: http://echo360.uml.edu/danylov03/physcssprng.html

More information

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

More information

Section 8.1 Exercises

Section 8.1 Exercises Secton 8.1 Non-rght Trangles: Law of Snes and Cosnes 519 Secton 8.1 Exercses Solve for the unknown sdes and angles of the trangles shown. 10 70 50 1.. 18 40 110 45 5 6 3. 10 4. 75 15 5 6 90 70 65 5. 6.

More information

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

Compilers. Spring term. Alfonso Ortega: Enrique Alfonseca: Chapter 4: Syntactic analysis

Compilers. Spring term. Alfonso Ortega: Enrique Alfonseca: Chapter 4: Syntactic analysis Complers Sprng term Alfonso Ortega: alfonso.ortega@uam.es nrque Alfonseca: enrque.alfonseca@uam.es Chapter : Syntactc analyss. Introducton. Bottom-up Analyss Syntax Analyser Concepts It analyses the context-ndependent

More information

In this section is given an overview of the common elasticity models.

In this section is given an overview of the common elasticity models. Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

11. Dynamics in Rotating Frames of Reference

11. Dynamics in Rotating Frames of Reference Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons

More information

on-line NotesBook Supplement

on-line NotesBook Supplement on-lne NotesBook Supplement Robot Mechancs Dr. Robert L. Wllams II Mechancal Engneerng Oho Unversty NotesBook Supplement for EE/ME 49/59 Mechancs and Control of Robotc Manpulators 18 Dr. Bob Productons

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

Lecture Note 3. Eshelby s Inclusion II

Lecture Note 3. Eshelby s Inclusion II ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product The Vector Product Week 11: Chapter 11 Angular Momentum There are nstances where the product of two vectors s another vector Earler we saw where the product of two vectors was a scalar Ths was called the

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

EN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics

EN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics N40: ynamcs and Vbratons Homewor 7: Rgd Body Knematcs School of ngneerng Brown Unversty 1. In the fgure below, bar AB rotates counterclocwse at 4 rad/s. What are the angular veloctes of bars BC and C?.

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

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

More information

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information

The optimal delay of the second test is therefore approximately 210 hours earlier than =2.

The optimal delay of the second test is therefore approximately 210 hours earlier than =2. THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple

More information

Department of Electrical & Electronic Engineeing Imperial College London. E4.20 Digital IC Design. Median Filter Project Specification

Department of Electrical & Electronic Engineeing Imperial College London. E4.20 Digital IC Design. Median Filter Project Specification Desgn Project Specfcaton Medan Flter Department of Electrcal & Electronc Engneeng Imperal College London E4.20 Dgtal IC Desgn Medan Flter Project Specfcaton A medan flter s used to remove nose from a sampled

More information

AP Physics 1 & 2 Summer Assignment

AP Physics 1 & 2 Summer Assignment AP Physcs 1 & 2 Summer Assgnment AP Physcs 1 requres an exceptonal profcency n algebra, trgonometry, and geometry. It was desgned by a select group of college professors and hgh school scence teachers

More information

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

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

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI 2017 2nd Internatonal Conference on Electrcal and Electroncs: echnques and Applcatons (EEA 2017) ISBN: 978-1-60595-416-5 Study on Actve Mcro-vbraton Isolaton System wth Lnear Motor Actuator Gong-yu PAN,

More information

Unit 5: Quadratic Equations & Functions

Unit 5: Quadratic Equations & Functions Date Perod Unt 5: Quadratc Equatons & Functons DAY TOPIC 1 Modelng Data wth Quadratc Functons Factorng Quadratc Epressons 3 Solvng Quadratc Equatons 4 Comple Numbers Smplfcaton, Addton/Subtracton & Multplcaton

More information

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

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

More information

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

COMPLEX NUMBERS AND QUADRATIC EQUATIONS COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not

More information

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN

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

More information

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results. Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson

More information

10. Canonical Transformations Michael Fowler

10. Canonical Transformations Michael Fowler 10. Canoncal Transformatons Mchael Fowler Pont Transformatons It s clear that Lagrange s equatons are correct for any reasonable choce of parameters labelng the system confguraton. Let s call our frst

More information

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

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

More information

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

More information

Physics 111: Mechanics Lecture 11

Physics 111: Mechanics Lecture 11 Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton

More information

CS 468 Lecture 16: Isometry Invariance and Spectral Techniques

CS 468 Lecture 16: Isometry Invariance and Spectral Techniques CS 468 Lecture 16: Isometry Invarance and Spectral Technques Justn Solomon Scrbe: Evan Gawlk Introducton. In geometry processng, t s often desrable to characterze the shape of an object n a manner that

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

Quantum Mechanics for Scientists and Engineers. David Miller

Quantum Mechanics for Scientists and Engineers. David Miller Quantum Mechancs for Scentsts and Engneers Davd Mller Types of lnear operators Types of lnear operators Blnear expanson of operators Blnear expanson of lnear operators We know that we can expand functons

More information

MAE140 - Linear Circuits - Winter 16 Midterm, February 5

MAE140 - Linear Circuits - Winter 16 Midterm, February 5 Instructons ME140 - Lnear Crcuts - Wnter 16 Mdterm, February 5 () Ths exam s open book. You may use whatever wrtten materals you choose, ncludng your class notes and textbook. You may use a hand calculator

More information

ENGI9496 Lecture Notes Multiport Models in Mechanics

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

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

σ τ τ τ σ τ τ τ σ Review Chapter Four States of Stress Part Three Review Review

σ τ τ τ σ τ τ τ σ Review Chapter Four States of Stress Part Three Review Review Chapter Four States of Stress Part Three When makng your choce n lfe, do not neglect to lve. Samuel Johnson Revew When we use matrx notaton to show the stresses on an element The rows represent the axs

More information

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

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

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

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

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

More information

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

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

More information

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1

More information

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

Quantum Mechanics I - Session 4

Quantum Mechanics I - Session 4 Quantum Mechancs I - Sesson 4 Aprl 3, 05 Contents Operators Change of Bass 4 3 Egenvectors and Egenvalues 5 3. Denton....................................... 5 3. Rotaton n D....................................

More information

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

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

More information

Affine and Riemannian Connections

Affine and Riemannian Connections Affne and Remannan Connectons Semnar Remannan Geometry Summer Term 2015 Prof Dr Anna Wenhard and Dr Gye-Seon Lee Jakob Ullmann Notaton: X(M) space of smooth vector felds on M D(M) space of smooth functons

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

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

More information

Homework Notes Week 7

Homework Notes Week 7 Homework Notes Week 7 Math 4 Sprng 4 #4 (a Complete the proof n example 5 that s an nner product (the Frobenus nner product on M n n (F In the example propertes (a and (d have already been verfed so we

More information

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015 Lecture 2. 1/07/15-1/09/15 Unversty of Washngton Department of Chemstry Chemstry 453 Wnter Quarter 2015 We are not talkng about truth. We are talkng about somethng that seems lke truth. The truth we want

More information

Lecture 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

More information

Fundamental loop-current method using virtual voltage sources technique for special cases

Fundamental loop-current method using virtual voltage sources technique for special cases Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information