Interpolated Rigid-Body Motions and Robotics
|
|
- Kerrie Strickland
- 5 years ago
- Views:
Transcription
1 Interpolate Rigi-Boy Motions an Robotics J.M. Selig Faculty of Business, Computing an Info. Management. Lonon South Bank University, Lonon SE AA, U.K. Yaunquing Wu Dept. Mechanical Engineering. Shanghai Jiaotong University, Shanghai, P.R. China. Abstract This work looks at several problems concerne with interpolating rigi-boy motions an their application in robotics. In particular a new type of motion is introuce. Two recently propose interpolation methos are shown to prouce the same results. We also iscuss how it might be possible to control a robot in such a way as to follow one of these interpolate motions. Inex Terms Rigi-boy motion, path planning, interpolate motion. I. INTRODUCTION This work has several objectives, the first of which is to correct a small mistake in [9]. This concerns the relations which gives the acceleration of points on the eneffector of a serial robot. The correction is given in the next section together with a brief introuction to the geometrical methos use in the paper. The secon objective is to introuce a new type of special rigi boy motion which we call the Bishop s motion. This is moele on the Frenet-Serret motion given in [3]. These motions coul be use for robot path planning. They shoul be particularly useful when human operators specify the esire path of the robot s en-effector since they are base on curves in 3 imensions an are hence easy to visualize. The Frenet-Serret motion is reviewe in section III an in the following section the Bishop s motion is introuce. Next we look at some connections with robotics. Specifically, we stuy how it is possible to guie a serial robot arm along a particular esire motion. Two answers are sketche, a numerical metho for the inverse kinematics an a non-linear control metho. Finally, we look at two recently propose interpolation schemes, [] an [4]. Although these methos approach the problem ifferently, the first focuses on the Lie groups while the secon concentrates on the paths of points in space, we are able to show that the two methos prouce ientical curves. II. SCREW THEORY Screw theory concerns the group of rigi boy motions SE3 an its Lie algebra se3. Elements of the groups, that is rigi boy motions, can be written as 4 4 matrices, R t g with R a 3 3 rotation matrix an t a translation vector. In this representation of the group it action on points can be written as follows. We enote a point in space by a vector p x, y, z T, then we exten this vector by aing an extra, p x y. z After the rigi transformation the position of the new point will be given by, p R t p Rp + t. Elements of the Lie algebra of the group correspon to velocities or small motions, they are calle twists or screws. In general, if we think of a path in the group as a continuously parameterise sequence of group elements gt, then for any parameter value we have a Lie algebra element given by, S t gt gt. If the group elements gt are given as 4 4 matrices then the Lie algebra element will have the form, Ω v S, where Ω is a 3 3 anti-symmetric matrix corresponing to the angular velocity of the motion. The vector v is a characteristic linear velocity of the motion. The action of the group on its Lie algebra elements is given by the conjugation, S gsg. This correspons to a left translation in the group. To see this suppose gt is a path in the group as above, now lefttranslate this by a constant group element h to get hgt. The Lie algebra elements of this path are now, h t gt gt h hsh. Clearly a right translation has no effect on the screw.
2 It is often convenient to write screws as 6-imensional Ω v vectors, corresponing to a screw S we can write a vector, ω s, v where the vector ω correspons to the anti-symmetric matrix Ω in the following way, for any vector u we have Ωu ω u. This 6-imensional representation of screws is useful in many cases. In particular it allows us to efine elements of the ual Lie algebra. Elements of this ual vector space are calle wrenches an written, M W, F where F is a force an M is a moment. These coul also be linear an angular momenta, if we were consiering ynamics. If a rigi boy is moving with instantaneous velocity given by the screw s an is acte on by a total force an moment given by the wrench W then the instantaneous power exerte on the boy is given by, power W T s M ω + F v. In robotics we can associate a screw to each joint in of the robot. The possible motions about the joint are given by the exponential, gθ e θs, where θ is the joint parameter, an angle if the joint is revolute. The superscript in the above refers to the position of the joint screw in the home position. The exponential of a matrix is given by the stanar formula, e S I + S + 2! S2 + 3! S3 + In general the position of a point p attache to the eneffector of the robot will be given by, pt e θs e θ 2S 2 e θ 3S 3 e θ 4S 4 e θ 5S 5 e θ 6S 6 p. Although the home positions of the joint screws Si are constants the current positions of these joint screws change, the kinematics of these are given by, S S S 2 e θs S 2 e θs S 3 e θ S e θ 2 S 2 S 3 e θ 2S 2 e θ S. S 6 e θ S e θ 2 S 2 e θ 5 S 5 S 6 e θ 5S 5 e θ S So the time erivatives of these screw are, t S since S is constant. t S 2 θ e θ S S S2 S2S e θ S θ [S, S 2 ] the last equality follows because S commutes with e θs. t S 3 θ e θ S S e θ 2S 2 S 3 e θ 2S 2 e θ S an so forth until, θ e θ S e θ 2 S 2 S 3 e θ 2S 2 S e θ S + θ 2 e θs e θ 2S 2 S 2 S 3e θ2s 2 e θ S θ 2 e θs e θ 2S 2 S 3 S 2e θ2s 2 e θ S θ [S, S 3 ] + θ 2 [S 2, S 3 ] t S 6 θ [S, S 6 ] + θ 2 [S 2, S 6 ] + θ 3 [S 3, S 6 ] + θ 4 [S 4, S 6 ] + θ 5 [S 5, S 6 ]. In [9, Chap 3.] a formula is given for the acceleration of a point attache to the en-effector of a serial robot. This formula contains an error, so we present a correcte version here. Differentiating the position of a point on the en-effector gives, ṗt θ Se θ S e θ 6 S 6 p + θ 6 e θs S 6 e θ6s 6 + p an this simplifies in the usual way to, ṗt θ S + + θ pt 6 S 6 This gives the velocity of the point. If we ifferentiate again to get the acceleration we have, { pt 6 θ i i S i + θ 6 i<j 6 i θj [S i, S j ] + θ 2 } pt i i S i We can simplify the above a little by expaning the square term, remembering that S i an S j o not generally commute. We get, pt { 6 i θi S i + θ i 2S2 i + 2 i<j 6 θ i θj S i S j } pt The terms θ i S i + θ i 2S2 i are the accelerations woul expect from motion about a single joint. A more comprehensive account of this view of screw theory can be foun in []. III. FRENET-SERRET MOTION In [3] Bottema an Roth stuy a number of special motions, one of which is the Frenet-Serret motion. Such a motion is etermine by a space-curve pt. Now in a Frenet-Serret motion a point in the moving boy moves along the curve an the coorinate frame in the moving boy remains aligne with the tangent t, normal n, an
3 binormal b, of the curve. Using the 4-imensional representation of SE3 the motion can be specifie as, Rt pt gt, where pt is the curve an the rotation matrix has the unit vectors t, n an b as columns, Rt t n b. 2 The famous Frenet-Serret relations are, ṫ vκn, 3 ṅ vκt + vτb, 4 ḃ vτn, 5 where v, κ an τ are respectively the spee, curvature an torsion of the curve. Our work here will be simplifie by introucing the Darboux vector ω vτt + vκb which has the properties that, ṫ ω t, ṅ ω n, ḃ ω b, 6 see [7, II.4] for example. This means that we can write, Ṙ ΩR, 7 where Ω is the 3 3 anti-symmetric matrix corresponing to ω. Hence we have that, Q t g g Ω vt ω p, 8 remember that ṗ vt. Using the Frenet-Serret relations 3 5 above the erivative of the velocity is, Q Ω vt ω p. 9 Here, ω vτ + v τt + vκ + v κb. In the 6-vector representation this is, ω q. vt ω p IV. BISHOP S MOTION In [2] Bishop gave an alternative metho to associate a moving frame to points on a curve in 3 imensions. In the same way that the Frenet-Serret frame etermines a special rigi boy motion etermine by a curve, the Bishop frame can also be use to efine a special motion. A point in the rigi boy follows a curve an an orthonormal frame in the boy stays aligne with the Bishop frame. Such a motion will be calle a Bishop s motion. There are some applications of the Bishop frame in Computer graphics to thicken curves. The Bishop frame is use because it oesn t twist about the curve, see Fig. b. In the Frenet-Serret frame the torsion cannot be efine at points where the curvature vanishes, this is not a problem for Bishop s frame. Another feature of the Bishop frame is that it is not unique, however once the initial frame has been selecte the frames for the rest of the curve are uniquely etermine. This suggests that the Bishop s motion efine above may be useful for robot path planning. The frame equations for Bishop s frame are: ṫ vk n + vk 2 n 2, ṅ vk t, 2 ṅ 2 vk 2 t. 3 As usual v is the spee of the curve an ṗ vt, where t is the unit tangent vector. The vectors n an n 2 are unit normal vectors an together with the tangent vector t, they form an orthonormal frame. So for example n n 2 t an so forth. The parameters k an k 2 are curvature-like functions. In terms of the Frenet-Serret frame the normal vectors n an n 2 are given by, n cos φn sin φb, 4 n 2 sin φn + cos φb. 5 where the angle φ is given by the inefinite integral, φ vτ t. 6 So, as mentione above, a curve oes not uniquely etermine a Bishop frame, there is a single rotational freeom in efining the Bishop frame, given by the constant of integration. But if we choose the unit normal vectors n an n 2 at t then the Bishop frame for the rest of the curve is unique. The path in the group etermine by a Bishop s motion will be the same as for the Frenet-Serret motion, see above, but now the rotation matrix will be given by, Rt t n n 2. 7 To compute the velocity of a Bishop s motion we nee an analogue of the Darboux vector. This is given by the vector, It is simple to verify that, a vk 2 n + vk n 2. 8 ṫ a t, ṅ a n, an ṅ 2 a n 2. Note that it is possible to show that a vκb. The velocity is, Q t g g A vt a p, 9 where, as usual, capital A represents the 3 3 antisymmetric matrix corresponing to the vector a. Proceeing as in the previous section we can compute the erivative of the velocity screw, as a 6-imensional vector this is ȧ q. 2 vt ȧ p
4 a b Fig.. a a Frenet-Serret motion an b a Bishop s motion. Both these motions were base on the same cubic spline curve which is approximate by the black tangent lines in either case. The initial frame for the Bishop s motion was chosen to coincie with the initial Frenet-Serret frame. So the frame at the bottom left of each iagram are the same. Notice how the blue binormal line in a rotates through 8 egrees as we move along the curve. However, in b the blue normal only turns through 9 egrees. V. ROBOT KINEMATICS AND CONTROL In this section we relate the esire motions to a couple of common problems in robotics. In general we want to rive the robot so that its en-effector follows the esire path. In the first case we look at computing the inverse kinematics along such a path. Let us write z 6, for the velocity screw of the en-effector of a robot. In terms of the joint variables this velocity screw is, 6 θ i s i z 6. 2 i This equation shows that the columns of the robot s Jacobian are given by the robot s joint screws s i. Now we can think of the inverse Jacobian of the manipulator as compose of six wrenches W, W 2,..., W 6 such that, { Wi T, if i j s j, if i j. 22 For many commercially available esigns of robot these quantities can be compute symbolically, see [, 6.7]. Setting the velocity of the robot s en-effector to the esire velocity gives, θ s + θ 2 s θ 6 s 6 q. 23 Pairing this equation with the wrenches W i, the rows of the inverse Jacobian of the robot, yiels, θ i W T i q, i, 2, These ifferentially equations shoul be straightforwar to set up an solve numerically using the Runge-Kutta metho for example. Notice that we nee to know the starting configuration of the robot but we only nee the velocity of the esire motion. Of course, like most methos in robotics, this metho will fail near singularities. As an alternative approach to guiing the robot along a esire path we can esign a close-loop controller for the robot. The scheme here is base on a non-linear feeback control scheme introuce in [8]. The iea of this is to control the en-effector of the robot along a specifie path without having to perform any inverse kinematic computations. The key is to arrange for the feeback to be such that the close loop ynamics of the system are ifferential equations in the robot s joint variables that are satisfie by the esire path. This iea is similar to the Passive Velocity Fiel Control metho introuce by Li an Horowitz[5]. However the metho presente here is somewhat simpler, since we o not introuce any extra egrees of freeom. Equation2 above can be ifferentiate to give, 6 θ i s i + i j<k 6 θ j θk [s j, s k ] ż Away from singularities we can pair this equation with the wrenches W i, we get, θ i + Γ ijk θj θk W T i ż6, i, 2,..., 6 26 where summing over repeate inices is assume. The quantities Γ ijk are given by, { Γ ijk 2 WT i [s j, s k ] j k 2 WT i [s 27 j, s k ] j k. If we replace the velocity of the en-effector with the velocity of the esire motion we obtain an equation for the esire motion in the joint space variables of the robot, θ i + Γ ijk θj θk W T i q, i, 2,..., Notice, the time erivative of a velocity screw is not exactly an acceleration. If the close-loop ynamics of the robot take this form then the en-effector will follow the esire motion, subject to consierations of stability an errors of course. Depening on whether we want to follow a Frenet-Serret motion of a Bishop s motion we must use q as given by either or 2. Now we turn to the ynamics of the robot. These can be represente by the equations, A ij θj + B ijk θj θj τ i, i, 2,..., 6 29
5 again summation over repeate inices is assume, see [, Chap. 3] for example. Here A ij is the generalise mass matrix an B ijk represents the coupling an Coriolis terms. The joint torques, which we assume that the control system can apply to robot s motors, are represente by τ i. For simplicity gravity will be ignore. Now suppose that our control system measures the joint variables an their rates, the motor torques can be set to, τ i B ijk A il Γ ljk θ j θk + A il W T l q, i,..., 6 3 This is a non-linear feeback control law. The close-loop equations of the system are foun by substituting for τ i in the robot ynamics equation 29 above. The result will be the equation for the esire motion 28 above multiplie by the positive-efinite mass matrix A ij. Notice that the evelopment of the control law above is base on the task space equation, ż 6 q. 3 For a stable controller it is probably better to base the metho on the equation, ż 6 q + λz 6 q, 32 where λ is a positive constant. Clearly, z 6 q is a solution to this equation. Moreover, if we assume that z 6 q + e, where e is an error vector, then the ynamics of the error vector obey, ė λe. 33 So the errors in the velocity ecrease exponentially. The joint space version of equation 32 is, θ i + Γ ijk θj θk W T i q + λw T i q θ i, i,..., 6 34 The corresponing close loop control law will be, τ i B ijk A il Γ ljk θ j θk + A il Wl T q + λa il Wl T q θ l, i,..., 6 35 VI. PROJECTION BASED INTERPOLATION Recently, two new methos of proucing rigi boy motions have been propose, see [] an [4]. Both these methos rely on projecting affine motions to rigi-boy motions. Belta an Kumar s metho interpolates the motions irectly to prouce a curve in the group of affine motions GA + 3 an then projects these motions to the nearest rigi-boy motion. On the other han, the metho escribe by, Hofer, Pottmann an Ravani is base on the paths of points in space. In this metho several points on the rigi boy are selecte an interpolation curves are efine for each of the points. Now at some particular instant the positions of the interpolate points may not be a rigi transformation of the original points. To prouce a rigi boy motion a least squares problem is solve to fin the positions of the points subject to the constraint that they are rigily relate to the original points. We show here that these methos are essentially the same. We begin with a number of rigi boy motions we wish to interpolate or approximate, assume these are given by, Ri t i, i,..., n. Now choose a number, say k, of points a j, at least 4 non-coplanar points accoring to [4]. The knot-points for the interpolation are then, b j i R i a j + t i, i,..., n. 36 So we get a set of interpolate curve, n p j t f i tb j i, j,..., k 37 i where f i t are the interpolating functions. In Belta an Kumar s metho we woul interpolate in the space of matrices to get, Mt t Xt Then clearly we have, p j t a j Xt n i f n itr i i f itt i. 38, j,..., k. 39 In general the matrix Xt will lie in the group GA + 3, this is because the 3 3 block Mt n i f itr i will generally be an element of GL + 3. However, it may happen that for some particular values of t, etmt, some conitions for this are stuie in []. The corresponing problem for points, that is values of t where the set of points p j t may be coplanar, oes not seem to have been anticipate in [4]. In Hofer et al s metho we fin the rigi motion at some time t by minimising the quantity k j pj t Ra j t 2. That is we seek a rotation matrix R an a translation vector t such that this expression is minimal. A solution to this problem can be foun in []. The first step is to choose the origin of coorinates to be at the centroi of the a j points, this is so that k j aj. In these coorinates the minimal translation is given by, t k k p j. 4 j To simplify notation we have roppe the explicit epenence on t, but we are assuming some efinite value of t. Now using the fact that p j Ma j + an the fact that the sum of the a j s is zero we have, t. 4 The minimal rotation can be foun from the polar ecomposition of the matrix, P k p j a j T. 42 j
6 a b Fig. 2. a a cubic projecte interpolate motion base on the en-points of the Bishop s motion given in Fig. b an b the Bishop s motion from Fig. b for comparison. This simplifies to, P M k a j a j T. 43 j The solution is the rotation R such that P RK where K is symmetric. Details on computing this an its relation to the singular value ecomposition of P can be foun in almost any numerical methos text. Belta an Kumar s metho is not really a single metho but several epening on a positive efinite symmetric matrix W. The metho requires us to fin the polar ecomposition of the matrix prouct MW. So if we ientify W with the matrix, k W a j a j T, 44 j the curves in SE3 given by the two approaches will be ientical. On the other han, Belta an Kumar give particular attention to the case where W I the ientity matrix. Now if we choose the points a j carefully then we can make the matrix k j aj a j T I. A suitable choice here woul be the vertices of a regular tetraheron, a 2/ 6 /2, a 2 / 6 / 2 3 /2, 3 a 3 / 6 / 2 /2, an a /2 In fact, from Schur s lemma, any collection of points a j which are symmetrical with respect to some finite subgroup of rotation crystallographic point groups, will give a matrix W that is a scalar multiple of the ientity. An example of such a motion is shown in Fig. 2a. VII. CONCLUSIONS In sections III an IV above we introuce two new types of rigi motion paths. The first of these, the Frenet- Serret motion, has been known in the mechanisms community for many years but oes not seem to have been applie to robots. The Bishop s motion oes not seem to be familiar to the mechanisms community. Neither of these motions is really an interpolate motion since once the path of a point is chosen the rotational motion is etermine. In the Bishop s motion we also have the freeom to choose an initial rotation. Nevertheless, these simple motions may be useful in practical circumstances. The purpose of section V was really to show that for applications to robotics we nee to know the velocity of the rigi motion curve. An sometimes also the erivative of the velocity. The control metho presente was rather simple an not to be taken too seriously. However, the iea of esigning the control system of a robot in such a way that it follows a esire trajectory using a knowlege of the geometry of the esire path rather than just employing inverse kinematics is a worthwhile goal. This means that interpolation techniques which result in curves whose velocities are ifficult to calculate are less esirable for applications in robotics. At the moment it is ifficult to see how to compute the velocity of a curve given by the two schemes stuie in section VI. REFERENCES [] C. Belta an V. Kumar, An SVD-projection metho for interpolation on SE3, IEEE Trans. Robotics an Automation, vol. 8, pp , 22. [2] R.L. Bishop, There is more than one way to frame a curve, Am. Math. Monthly, vol. 82 pp , 975. [3] O. Bottema an B. Roth, Theoretical Kinematics, Dover Publications, New York, 99. [4] M. Hofer, H. Pottmann, an B. Ravani. From curve esign algorithms to the esign of rigi boy motions. The Visual Computer vol. 2, pp , 24. [5] P.Y. Li an R. Horowitz, Passive velocity fiel control of mechanical manipulators, IEEE Trans. Robotics an Automation, vol. 5, pp , 999. [6] D. Marsh, Applie Geometry for Computer Graphics an CAD 2n. e., Springer Verlag, Lonon, 25. [7] B. O Neill, Elementary Differential Geometry, Acaemic Press, New York, 966. [8] J.M. Selig an A.I. Ovseevitch, Manipulating robots along helical trajectories, Robotica vol. 4 pp , 996. [9] J.M. Selig. e. Geometric founations of robotics Worl Scientific, Singapore, 2. [] J.M. Selig. Three Problems in Robotics, Proc. Inst. Mech. Eng. part C: J. Mechanical Engineering Science, vol 26, pp.73 8, 22. [] J.M. Selig. Geometric Funamentals of Robotics. Springer Verlag, New York, 25.
Interpolated Rigid-Body Motions and Robotics
Interpolated Rigid-Body Motions and Robotics J.M. Selig London South Bank University and Yuanqing Wu Shanghai Jiaotong University. IROS Beijing 2006 p.1/22 Introduction Interpolation of rigid motions important
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationFree rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012
Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationTutorial Test 5 2D welding robot
Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.
More informationNumerical Integrator. Graphics
1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationDarboux s theorem and symplectic geometry
Darboux s theorem an symplectic geometry Liang, Feng May 9, 2014 Abstract Symplectic geometry is a very important branch of ifferential geometry, it is a special case of poisson geometry, an coul also
More informationPhysics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1
Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More information3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes
Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationIntegration Review. May 11, 2013
Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In
More informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationTotal Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*
51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationand from it produce the action integral whose variation we set to zero:
Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationCalculus in the AP Physics C Course The Derivative
Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.
More informationG j dq i + G j. q i. = a jt. and
Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationCHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold
CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the
More informationFrom Local to Global Control
Proceeings of the 47th IEEE Conference on Decision an Control Cancun, Mexico, Dec. 9-, 8 ThB. From Local to Global Control Stephen P. Banks, M. Tomás-Roríguez. Automatic Control Engineering Department,
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationMath 1B, lecture 8: Integration by parts
Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores
More informationDifferentiability, Computing Derivatives, Trig Review. Goals:
Secants vs. Derivatives - Unit #3 : Goals: Differentiability, Computing Derivatives, Trig Review Determine when a function is ifferentiable at a point Relate the erivative graph to the the graph of an
More informationDifferentiability, Computing Derivatives, Trig Review
Unit #3 : Differentiability, Computing Derivatives, Trig Review Goals: Determine when a function is ifferentiable at a point Relate the erivative graph to the the graph of an original function Compute
More informationEuler Equations: derivation, basic invariants and formulae
Euler Equations: erivation, basic invariants an formulae Mat 529, Lesson 1. 1 Derivation The incompressible Euler equations are couple with t u + u u + p = 0, (1) u = 0. (2) The unknown variable is the
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationTRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS
TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email
More informationShort Intro to Coordinate Transformation
Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationComputing Derivatives
Chapter 2 Computing Derivatives 2.1 Elementary erivative rules Motivating Questions In this section, we strive to unerstan the ieas generate by the following important questions: What are alternate notations
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationBEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica
More informationMany problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx
Math 53 Notes on turm-liouville equations Many problems in physics, engineering, an chemistry fall in a general class of equations of the form w(x)p(x) u ] + (q(x) λ) u = w(x) on an interval a, b], plus
More informationMathematical Review Problems
Fall 6 Louis Scuiero Mathematical Review Problems I. Polynomial Equations an Graphs (Barrante--Chap. ). First egree equation an graph y f() x mx b where m is the slope of the line an b is the line's intercept
More informationNotes on Lie Groups, Lie algebras, and the Exponentiation Map Mitchell Faulk
Notes on Lie Groups, Lie algebras, an the Exponentiation Map Mitchell Faulk 1. Preliminaries. In these notes, we concern ourselves with special objects calle matrix Lie groups an their corresponing Lie
More informationVIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK
AIAA Guiance, Navigation, an Control Conference an Exhibit 5-8 August, Monterey, California AIAA -9 VIRTUAL STRUCTURE BASED SPACECRAT ORMATION CONTROL WITH ORMATION EEDBACK Wei Ren Ranal W. Bear Department
More informationModule FP2. Further Pure 2. Cambridge University Press Further Pure 2 and 3 Hugh Neill and Douglas Quadling Excerpt More information
5548993 - Further Pure an 3 Moule FP Further Pure 5548993 - Further Pure an 3 Differentiating inverse trigonometric functions Throughout the course you have graually been increasing the number of functions
More informationSome vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10
Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the
More information23 Implicit differentiation
23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For
More information1 Heisenberg Representation
1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.
More informationAn Explicit Formulation of Singularity-Free Dynamic Equations of Mechanical Systems in Lagrangian Form Part Two: Multibody Systems
Moeng, Ientification an Control, Vol 33, No 2, 2012, pp 61 68 An Expcit Formulation of Singularity-Free Dynamic Equations of Mechanical Systems in Lagrangian Form Part Two: Multiboy Systems Pål Johan From
More informationAccelerate Implementation of Forwaring Control Laws using Composition Methos Yves Moreau an Roolphe Sepulchre June 1997 Abstract We use a metho of int
Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SISTA/TR 1997-11 Accelerate Implementation of Forwaring Control Laws using Composition Methos 1 Yves Moreau, Roolphe Sepulchre, Joos Vanewalle
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationMake graph of g by adding c to the y-values. on the graph of f by c. multiplying the y-values. even-degree polynomial. graph goes up on both sides
Reference 1: Transformations of Graphs an En Behavior of Polynomial Graphs Transformations of graphs aitive constant constant on the outsie g(x) = + c Make graph of g by aing c to the y-values on the graph
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationDIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10
DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10 5. Levi-Civita connection From now on we are intereste in connections on the tangent bunle T X of a Riemanninam manifol (X, g). Out main result will be a construction
More informationFinal Exam Study Guide and Practice Problems Solutions
Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making
More informationProblem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs
Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationPhysics 251 Results for Matrix Exponentials Spring 2017
Physics 25 Results for Matrix Exponentials Spring 27. Properties of the Matrix Exponential Let A be a real or complex n n matrix. The exponential of A is efine via its Taylor series, e A A n = I + n!,
More informationTOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH
English NUMERICAL MATHEMATICS Vol14, No1 Series A Journal of Chinese Universities Feb 2005 TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH He Ming( Λ) Michael K Ng(Ξ ) Abstract We
More informationFurther Differentiation and Applications
Avance Higher Notes (Unit ) Prerequisites: Inverse function property; prouct, quotient an chain rules; inflexion points. Maths Applications: Concavity; ifferentiability. Real-Worl Applications: Particle
More information4. Important theorems in quantum mechanics
TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this
More informationCalculus Class Notes for the Combined Calculus and Physics Course Semester I
Calculus Class Notes for the Combine Calculus an Physics Course Semester I Kelly Black December 14, 2001 Support provie by the National Science Founation - NSF-DUE-9752485 1 Section 0 2 Contents 1 Average
More informationu t v t v t c a u t b a v t u t v t b a
Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying
More informationIn the usual geometric derivation of Bragg s Law one assumes that crystalline
Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this
More informationGyroscopic matrices of the right beams and the discs
Titre : Matrice gyroscopique es poutres roites et es i[...] Date : 15/07/2014 Page : 1/16 Gyroscopic matrices of the right beams an the iscs Summary: This ocument presents the formulation of the matrices
More informationOptimal Variable-Structure Control Tracking of Spacecraft Maneuvers
Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign
More informationImplicit Lyapunov control of closed quantum systems
Joint 48th IEEE Conference on Decision an Control an 28th Chinese Control Conference Shanghai, P.R. China, December 16-18, 29 ThAIn1.4 Implicit Lyapunov control of close quantum systems Shouwei Zhao, Hai
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationPhysics 2212 GJ Quiz #4 Solutions Fall 2015
Physics 2212 GJ Quiz #4 Solutions Fall 215 I. (17 points) The magnetic fiel at point P ue to a current through the wire is 5. µt into the page. The curve portion of the wire is a semicircle of raius 2.
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More informationKinematics of Rotations: A Summary
A Kinematics of Rotations: A Summary The purpose of this appenix is to outline proofs of some results in the realm of kinematics of rotations that were invoke in the preceing chapters. Further etails are
More informationMath 1271 Solutions for Fall 2005 Final Exam
Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly
More informationx f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.
Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?
More informationPermanent vs. Determinant
Permanent vs. Determinant Frank Ban Introuction A major problem in theoretical computer science is the Permanent vs. Determinant problem. It asks: given an n by n matrix of ineterminates A = (a i,j ) an
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationNewton Euler equations in general coordinates
Newton Euler equations in general coorinates y ertol ongar an Frank Kirchner Robotics Innoation Center, DFKI GmbH, remen, Germany Abstract For the computation of rigi boy ynamics, the Newton Euler equations
More informationTopic 2.3: The Geometry of Derivatives of Vector Functions
BSU Math 275 Notes Topic 2.3: The Geometry of Derivatives of Vector Functions Textbook Sections: 13.2 From the Toolbox (what you nee from previous classes): Be able to compute erivatives scalar-value functions
More informationExperimental Robustness Study of a Second-Order Sliding Mode Controller
Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans
More informationYear 11 Matrices Semester 2. Yuk
Year 11 Matrices Semester 2 Chapter 5A input/output Yuk 1 Chapter 5B Gaussian Elimination an Systems of Linear Equations This is an extension of solving simultaneous equations. What oes a System of Linear
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More informationAN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A
AN INTRODUCTION TO AIRCRAFT WIN FLUTTER Revision A By Tom Irvine Email: tomirvine@aol.com January 8, 000 Introuction Certain aircraft wings have experience violent oscillations uring high spee flight.
More informationA Spectral Method for the Biharmonic Equation
A Spectral Metho for the Biharmonic Equation Kenall Atkinson, Davi Chien, an Olaf Hansen Abstract Let Ω be an open, simply connecte, an boune region in Ê,, with a smooth bounary Ω that is homeomorphic
More informationAdaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements
Aaptive Gain-Scheule H Control of Linear Parameter-Varying Systems with ime-delaye Elements Yoshihiko Miyasato he Institute of Statistical Mathematics 4-6-7 Minami-Azabu, Minato-ku, okyo 6-8569, Japan
More informationON THE RIEMANN EXTENSION OF THE SCHWARZSCHILD METRICS
ON THE RIEANN EXTENSION OF THE SCHWARZSCHILD ETRICS Valerii Dryuma arxiv:gr-qc/040415v1 30 Apr 004 Institute of athematics an Informatics, AS R, 5 Acaemiei Street, 08 Chisinau, olova, e-mail: valery@ryuma.com;
More informationSolutions to Practice Problems Tuesday, October 28, 2008
Solutions to Practice Problems Tuesay, October 28, 2008 1. The graph of the function f is shown below. Figure 1: The graph of f(x) What is x 1 + f(x)? What is x 1 f(x)? An oes x 1 f(x) exist? If so, what
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationIntegration by Parts
Integration by Parts 6-3-207 If u an v are functions of, the Prouct Rule says that (uv) = uv +vu Integrate both sies: (uv) = uv = uv + u v + uv = uv vu, vu v u, I ve written u an v as shorthan for u an
More informationUnit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule
Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin
More informationThe Pitch, the Angle of Pitch, and the Distribution Parameter of a Closed Ruled Surface
Appl. Math. Inf. Sci. 9 No. 4 809-85 05 809 Applie Mathematics & Information Sciences An International Journal http://x.oi.org/0.785/amis/09048 The Pitch the Angle of Pitch the Distribution Parameter of
More informationA. Incorrect! The letter t does not appear in the expression of the given integral
AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question
More informationTRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM
265 Asian Journal of Control, Vol. 4, No. 3, pp. 265-273, September 22 TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM Jurachart Jongusuk an Tsutomu Mita
More information