Classical Mechanics Examples (Lagrange Multipliers)

Size: px
Start display at page:

Download "Classical Mechanics Examples (Lagrange Multipliers)"

Transcription

1 Classical Mechanics Examples (Lagrange Multipliers) Dipan Kumar Ghosh Physics Department, Inian Institute of Technology Bombay Powai, Mumbai September 3, Introuction We have seen that the effect of holonomic constraints is to reuce the number of egrees of freeom. If there are N particles, we expect 3N egrees of freeom if the motion is unconstraine. If there are k holonomic constraints, the number of egrees of freeom becomes 3N k an we nee this many inepenent generalise coorinates to escribe the system. For non-holonomic constraints, however, the number of generalise coorinates require is more than the number of egrees of freeom. If there are m non-holonomic constraints, in aition to the k holonomic constraints, the number of generalise coronates is still 3N k, i.e. m more than the number of egrees of freeom. These aitional quantities are eliminate from the equations of motion by metho of Lagrange Multipliers. 1.1 A mathematical Example: Fin extremal value of f(x, y) = xy subject to the constraint h(x, y) = x 8 + y 1 = 0 This is an equation to al ellipse with major axis an minor axis. The level curves for f(x, y) = xy are hyperbola. Suppose the level curve for f(x, y) intersects the constraint curve h(x, y) at P. If we move along the constraint curve to its right, the value of xy increases whereas if we move to the left its value ecreases. Thus the point P cannot be a point of minimum or a maximum. Clearly, when the point is an extremum, the level curve of f an the constraint curve touch each other an they have a common tangent. Hence at the common point their normals are parallel (or antiparallel) to each other. Thus the set of points S which satisfies h(x, y) = 0 an f + λ h = 0 for some λ contains the extremal value of f, subject to the given constraint.

2 c D. K. Ghosh, IIT Bombay We efine Setting graient of F equal to zero implies f x + λ h ˆx + x which gives us F (x, y, λ) = f(x, y) + λh(x, y) ( f y + λ h y f x + λ h x = 0 f y + λ h y = 0 ) ŷ + hλ = 0 For the present problem, we have, x F (x, y, h) = xy λ 8 + y 1 In this case from (1) we get h = 0 (1) y λ 4 x = 0 x λy = 0 x 8 + y 1 = 0

3 c D. K. Ghosh, IIT Bombay 3 which gives λ = ±. Substituting the possible values of λ into the above equations, we fin two maximum values of f = xy = at (1, ), ( 1, ) an two minimum values xy = at (1, ), ( 1, ). Euler Lagrange Equation with constraint The metho of Lagrange multiplier is use in classical mechanics in hanling situations where the number of ynamical variables happens to be more than the number of egrees of freeom. We ha seen that the principle of least action is expresse in the form δ L = 0 from which we ha obtaine the Euler Lagrange equation as j [ L ] δq j = 0 Using the fact that the generalise coorinates are inepenent, we erive the Euler Lagrange equation for each pair pf generalise coorinate an velocity L = 0 () In the presence of non-holonomic constraints, the generalise coorinates are not inepenent as their number is greater than the number of egrees of freeom an consequently () is not vali. In such situations (an also in case care is not taken to reuce the number of generalise coorinates using the holonomic constraints), the metho of unetermine multiplier is useful. (The metho has limitations, for instance, it cannot be use in cases where the constraints are state as inequalities) Suppose there are m number of non-holonomic constraints involving the generalise coorinates in ifferential form n a r,j q j + b r,j = 0 (3) j=1 a r,j an b r,j may epen on j an t. Here r is an inex which runs from 1 to m an (3) is actually m equations, one for each value of r. We can get correct equations motion if the varie paths are virtual isplacements from actual motion in which case the constraint is n j=1 a r,jq = 0 because vital isplacements take place over constant time. We can then rewrite the principle of least action as n j j=1 [ L + ] λ r a r,j δq j = 0 Note that the aitional term n j=1 m λ ra r,j δq j is actually zero an hence we can put it insie the integral. δq j s are not inepenent an satisfy (3). Since we have m

4 c D. K. Ghosh, IIT Bombay 4 unetermine multipliers λ r, we can choose them such that the first m terms, i.e. j = 1 to m is each zero. Suppose we choose λ j such that for j = 1,,... m, the equation to be satisfie is L + λ r a r,j = 0 (4) Note that the last term in the above equation is no longer zero as the sum over j is missing an we have simply reistribute the term which was zero in various ways. With the λ r etermine by (4), we are left with n j=m+1 [ L + ] λ r a r,j δq j = 0 However, now our q j are inepenent an we have, as a consequence, for j = m + 1,..., n L + (4) an (5) allows us to write a single equation for j = 1,,... n, L + λ r a r,j = 0 (5) λ r a r,j = 0 (6) (4) etermines the m values of λ an (5) gives us n equations of motion. Define Q j = λ r a r,j (7) Equation (5) can now be written as for j = 1,,... n. L = Q j (8) The right han sie of the above is seen to be the generalise force corresponing to the constraint conitions..1 Example : A sliing mass on a paraboloi of revolution Consier a particle of mass m moving frictionlessly on a paraboloi of revolution given by x + y = az.

5 c D. K. Ghosh, IIT Bombay 5 z φ ^z ^φ m ρ^ x The Lagrangian of the system is L = T V = 1 m[ ρ + ρ ϕ + ż ] + mgz (9) where ρ = x + y is the particle from the z axis. The system is holonomic with two egrees of freeom since the constraint ρ = az can be use to eliminate one of the coorinates. Let us write the constraint as ρδρ aδz = 0 (10) Comparing (10) with (3) we ientify a ρ = ρ, a ϕ = 0 an a z = a. As there is only one constraint equation, there is only one Lagrange multiplier λ an the Euler Lagrange equation is L = Q j = λa j The Euler Lagrange equations for the three coorinates are L ρ ρ = λ ρ L ϕ ϕ = 0 L ż z = λ a These equations simplify to

6 c D. K. Ghosh, IIT Bombay 6 m( ρ ρ ϕ ) = λρ m (ρ ϕ) = 0 = mρ ω 0 = L m z = mg λa We can solve these equations to etermine ρ(t), z(t), ϕ(t) an λ. Consier the case where the mass goes aroun in a horizontal circe at a height h from the bottom. z = h. This also implies ρ = constant = ha. We then have (from the last equation) λ = mg a. ϕ = constant = ω 0 = g (from the first equation, using ρ = 0) a λ = mg a = m ω 0 The force which is responsible for the circular motion is the component N sin ψ of the normal force. tan ψ = z ρ = ρ The component of the normal force towars the centre of a the circle is ρ N sin ψ = N 4ρ + a mρω 0 The vertical component of N must be equal to mg a N cos ψ = N 4ρ + a = mg Thus which is what we expecte. λ = mg a = N 4ρ + a = 1 mω 0. A hoop rolling own the top of a cyliner without slipping Consier the case of a hoop rolling own the top of a cyliner without slipping. The constraint is holonomic till the mass slies off. Let r be the istance from O to C. θ is the polar coorinate of O an ϕ is the angle of rotation of the hoop about its own axis. R is the raius of the cyliner. O θ R C V=0

7 c D. K. Ghosh, IIT Bombay 7 The Lagrangian is L = 1 m(ṙ + r θ ) + 1 ma ϕ mgr cos θ Constraints are : (i) r = R + a (ii) (R + a) θ = a ϕ The first is a holonomic constraint : f 1 = r (R + a) an the secon is a non-holonomic (though integrable) constraint f = (R + a) θ + a ϕ = 0. With the former we associate λ (associate with normal reaction) while with the latter we associate µ (associate with the tangential force require for rolling motion). The Euler Lagrange equations (for r, θ an ϕ) are L The equations give = λ f 1 + µ f m r mr θ + mg cos θ = λ mr θ + mrṙ θ mgr sin θ = µ(r + a) ma ϕ = µa Hence µ = ma ϕ = m(r + a) θ Substitute this in the secon equation an use ṙ = 0. θ = g sin θ (R + a) We multiply the above equation by θ an integrate θ = Use this in the raial equation (put ṙ = 0) g (1 cos θ) R + a λ = mr θ + mg cos θ g = m(r + a) (1 cos θ) + mg cos θ R + a = mg( cos θ 1) For θ = 0, λ = mg, the normal force. For θ > 60, λ becomes negative an contact with the cyliner is lost..3 Example : A Rolling hoop on an incline Consier a hoop rolling own without slipping on an incline. The problem is one of holonomic constraint an has been solve earlier. However, we will use this problem to illustrate the metho of Lagrange multiplier. In this case there is just one egree of freeom an there is just one generalise coorinate, i.e. the angle θ by which the hoop rolls. The istance move along the incline is x = Rθ.

8 c D. K. Ghosh, IIT Bombay 8 x R θ α Let us, however, use both x an θ are our generalise coorinates. We then have one holonomic constraint h(x, θ) = x Rθ (11) An observation may be mae here that (11) is not in the form (3). A holonomic constraint is expresse in the form f(q 1, q,..., q n, t) = 0 We can write it in a ifferential form n f δq k + f δt = 0 (1) q k t Comparing this with (3) we have k=1 a r,k = f q k b r,t = f t Thus in the present case, where we have only one constraint (11) (we rop the inex r) The Lagrangian is then given by b = 0 a x = (x Rθ) = 1 x a θ = (x Rθ) = R θ L = 1 mẋ + 1 mr θ mg(l x) sin α where l is the length of the incline. The Euler-Lagrange equations for the coorinate x an θ are given by L ẋ x λ 1 = 0 i.e. mẍ mg sin α λ = 0 (14) an L θ θ + λ R = 0 (13)

9 c D. K. Ghosh, IIT Bombay 9 i.e. Equations (14), (15) along with the equation to the constraint mr θ + λr = 0 (15) ẋ R θ = 0 (16) provie the necessary equations for solution of the problem. From (15) an (16) we get λ = mẍ. Substituting these in (14) we get ẍ = g sin α mg sin α λ = The force of constraint corresponing to the motion along x irection is given by λ h x = mg sin α mg sin α 1 =.4 Rolling with Spinning introucing a non-holonomic constraint Let us now allow the hoop to spin about the vertical axis, in aition to rolling. Let us take the irection own the incline as the x-irection an the irection perpenicular to it along the plane as the y- irection. (17) (18) R θ ϕ The kinetic energy of the hoop is given by α T = 1 mv + 1 mr θ + 1 mr = mr θ mr ϕ Here we have use the fact that the hoop spins about its iameter an the moment of inertia about the iameter is mr /. The Lagrangian is L = mr θ mr ϕ mg(l x) sin α ϕ

10 c D. K. Ghosh, IIT Bombay 10 Here y is cyclic an ẏ oes not enter into the Hamiltonian. Hence we ignore the equation of motion for y. Since the hoop rolls without sliing, the velocity of the rim is v = R θ. When the hoop has spanne by an angle ϕ, we have ẋ = R θ cos ϕ (19) ẏ = R θ sin ϕ (0) These are non-integrable constraints because these cannot be use to connect x with θ an ϕ. The problem has two egrees of freeom. Let us freeze the motion an make virtual isplacement δθ of the angle θ(a isplacement in tilt angle is of secon orer). We then have δx = R cos ϕδθ As before we can get a x = 1 an a θ = R cos ϕ. We can now write the Lagrange equations of motion an get the following equations: mg sin α λ = 0 (1) mr θ + λr cos ϕ = 0 () mr ϕ = 0 (3) These equations etermine the motion completely in terms of θ an ϕ. Note that x has been alienate from the equations of motion.

5th International Conference on Advanced Design and Manufacturing Engineering (ICADME 2015)

5th International Conference on Advanced Design and Manufacturing Engineering (ICADME 2015) 5th International Conference on Avance Design an Manufacturing Engineering (ICADME 25) Research on motion characteristics an application of multi egree of freeom mechanism base on R-W metho Xiao-guang

More information

Impact of changing the position of the tool point on the moving platform on the dynamic performance of a 3RRR planar parallel manipulator

Impact of changing the position of the tool point on the moving platform on the dynamic performance of a 3RRR planar parallel manipulator IOSR Journal of Mechanical an Civil Engineering (IOSR-JMCE) e-issn: 78-84,p-ISSN: 0-4X, Volume, Issue 4 Ver. I (Jul. - Aug. 05), PP 7-8 www.iosrjournals.org Impact of changing the position of the tool

More information

Linear First-Order PDEs

Linear First-Order PDEs MODULE 2: FIRST-ORDER PARTIAL DIFFERENTIAL EQUATIONS 9 Lecture 2 Linear First-Orer PDEs The most general first-orer linear PDE has the form a(x, y)z x + b(x, y)z y + c(x, y)z = (x, y), (1) where a, b,

More information

(a) Find the equation of the plane that passes through the points P, Q, and R.

(a) Find the equation of the plane that passes through the points P, Q, and R. Math 040 Miterm Exam 1 Spring 014 S o l u t i o n s 1 For given points P (, 0, 1), Q(, 1, 0), R(3, 1, 0) an S(,, 0) (a) Fin the equation of the plane that passes through the points P, Q, an R P Q = 0,

More information

2.7 Implicit Differentiation

2.7 Implicit Differentiation 2.7 Implicit Differentiation [ y] = = = [ x] 1 Sometimes we may be intereste in fining the erivative of an equation that is not solve or able to be solve for a particular epenent variable explicitly. In

More information

Math 209 (Fall 2007) Calculus III. Solution #5. 1. Find the minimum and maximum values of the following functions f under the given constraints:

Math 209 (Fall 2007) Calculus III. Solution #5. 1. Find the minimum and maximum values of the following functions f under the given constraints: Math 9 (Fall 7) Calculus III Solution #5. Find the minimum and maximum values of the following functions f under the given constraints: (a) f(x, y) 4x + 6y, x + y ; (b) f(x, y) x y, x + y 6. Solution:

More information

Module13:Interference-I Lecture 13: Interference-I

Module13:Interference-I Lecture 13: Interference-I Moule3:Interference-I Lecture 3: Interference-I Consier a situation where we superpose two waves. Naively, we woul expect the intensity (energy ensity or flux) of the resultant to be the sum of the iniviual

More information

Bends, Jogs, And Wiggles for Railroad Tracks and Vehicle Guide Ways

Bends, Jogs, And Wiggles for Railroad Tracks and Vehicle Guide Ways Ben, Jogs, An Wiggles for Railroa Tracks an Vehicle Guie Ways Louis T. Klauer Jr., PhD, PE. Work Soft 833 Galer Dr. Newtown Square, PA 19073 lklauer@wsof.com Preprint, June 4, 00 Copyright 00 by Louis

More information

Double Integrals, Iterated Integrals, Cross-sections

Double Integrals, Iterated Integrals, Cross-sections Chapter 14 Multiple Integrals 1 ouble Integrals, Iterated Integrals, Cross-sections 2 ouble Integrals over more general regions, efinition, Evaluation of ouble Integrals, Properties of ouble Integrals

More information

Exercises of PIV. incomplete draft, version 0.0. October 2009

Exercises of PIV. incomplete draft, version 0.0. October 2009 Exercises of PIV incomplete raft, version 0.0 October 2009 1 Images Images are signals efine in 2D or 3D omains. They can be vector value (e.g., color images), real (monocromatic images), complex or binary

More information

4.2 Implicit Differentiation

4.2 Implicit Differentiation 6 Chapter 4 More Derivatives 4. Implicit Differentiation What ou will learn about... Implicitl Define Functions Lenses, Tangents, an Normal Lines Derivatives of Higher Orer Rational Powers of Differentiable

More information

Section 19. Thin Prisms

Section 19. Thin Prisms Section 9 Thin Prisms 9- OPTI-50 Optical Design an Instrumentation I opyright 08 John E. Greivenkamp Thin Prism Deviation Thin prisms introuce small angular beam eviations an are useful as alignment evices.

More information

Section 20. Thin Prisms

Section 20. Thin Prisms OPTI-0/0 Geometrical an Instrumental Optics opyright 08 John E. Greivenkamp 0- Section 0 Thin Prisms Thin Prism Deviation Thin prisms introuce small angular beam eviations an are useful as alignment evices.

More information

Parts Assembly by Throwing Manipulation with a One-Joint Arm

Parts Assembly by Throwing Manipulation with a One-Joint Arm 21 IEEE/RSJ International Conference on Intelligent Robots an Systems, Taipei, Taiwan, October, 21. Parts Assembly by Throwing Manipulation with a One-Joint Arm Hieyuki Miyashita, Tasuku Yamawaki an Masahito

More information

10. WAVE OPTICS ONE MARK QUESTIONS

10. WAVE OPTICS ONE MARK QUESTIONS 1 10. WAVE OPTICS ONE MARK QUESTIONS 1. Define wavefront.. What is the shape of wavefront obtaine from a point source at a (i) small istance (ii) large istance? 3. Uner what conitions a cylinrical wavefront

More information

Applications of Triple Integrals

Applications of Triple Integrals Chapter 14 Multiple Integrals 1 Double Integrals, Iterated Integrals, Cross-sections 2 Double Integrals over more general regions, Definition, Evaluation of Double Integrals, Properties of Double Integrals

More information

Constrained Optimization and Lagrange Multipliers

Constrained Optimization and Lagrange Multipliers Constrained Optimization and Lagrange Multipliers MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Constrained Optimization In the previous section we found the local or absolute

More information

FINDING OPTICAL DISPERSION OF A PRISM WITH APPLICATION OF MINIMUM DEVIATION ANGLE MEASUREMENT METHOD

FINDING OPTICAL DISPERSION OF A PRISM WITH APPLICATION OF MINIMUM DEVIATION ANGLE MEASUREMENT METHOD Warsaw University of Technology Faculty of Physics Physics Laboratory I P Joanna Konwerska-Hrabowska 6 FINDING OPTICAL DISPERSION OF A PRISM WITH APPLICATION OF MINIMUM DEVIATION ANGLE MEASUREMENT METHOD.

More information

Research Article Inviscid Uniform Shear Flow past a Smooth Concave Body

Research Article Inviscid Uniform Shear Flow past a Smooth Concave Body International Engineering Mathematics Volume 04, Article ID 46593, 7 pages http://x.oi.org/0.55/04/46593 Research Article Invisci Uniform Shear Flow past a Smooth Concave Boy Abullah Mura Department of

More information

Calculation on diffraction aperture of cube corner retroreflector

Calculation on diffraction aperture of cube corner retroreflector November 10, 008 / Vol., No. 11 / CHINESE OPTICS LETTERS 8 Calculation on iffraction aperture of cube corner retroreflector Song Li (Ó Ø, Bei Tang (», an Hui Zhou ( ï School of Electronic Information,

More information

Figure 1: 2D arm. Figure 2: 2D arm with labelled angles

Figure 1: 2D arm. Figure 2: 2D arm with labelled angles 2D Kinematics Consier a robotic arm. We can sen it commans like, move that joint so it bens at an angle θ. Once we ve set each joint, that s all well an goo. More interesting, though, is the question of

More information

Math 131. Implicit Differentiation Larson Section 2.5

Math 131. Implicit Differentiation Larson Section 2.5 Math 131. Implicit Differentiation Larson Section.5 So far we have ealt with ifferentiating explicitly efine functions, that is, we are given the expression efining the function, such as f(x) = 5 x. However,

More information

TEST-05(Solution) TOPIC: OPTICS COMPLETE

TEST-05(Solution) TOPIC: OPTICS COMPLETE Q. boy is walking uner an incline mirror at a constant velocity V m/s along the x-axis as shown in figure. If the mirror is incline at an angle with the horizontal then what is the velocity of the image?

More information

To graph the point (r, θ), simply go out r units along the initial ray, then rotate through the angle θ. The point (1, 5π 6. ) is graphed below:

To graph the point (r, θ), simply go out r units along the initial ray, then rotate through the angle θ. The point (1, 5π 6. ) is graphed below: Polar Coordinates Any point in the plane can be described by the Cartesian coordinates (x, y), where x and y are measured along the corresponding axes. However, this is not the only way to represent points

More information

Optimal path planning in a constant wind with a bounded turning rate

Optimal path planning in a constant wind with a bounded turning rate Optimal path planning in a constant win with a boune turning rate Timothy G. McGee, Stephen Spry an J. Karl Herick Center for Collaborative Control of Unmanne Vehicles, University of California, Berkeley,

More information

Calculus III. Math 233 Spring In-term exam April 11th. Suggested solutions

Calculus III. Math 233 Spring In-term exam April 11th. Suggested solutions Calculus III Math Spring 7 In-term exam April th. Suggested solutions This exam contains sixteen problems numbered through 6. Problems 5 are multiple choice problems, which each count 5% of your total

More information

Chapter 35 Homework (due 12/03/13)!!

Chapter 35 Homework (due 12/03/13)!! Chapter 35 Homework (ue 12/03/13) 35.6 35.20 35.21 35.29 35.51 35.55 35.66 35.68 page 1 Problem 35.6 An unerwater scuba iver sees the sun an in the parent angle of 45 above the horizontal. What is the

More information

0607 CAMBRIDGE INTERNATIONAL MATHEMATICS

0607 CAMBRIDGE INTERNATIONAL MATHEMATICS CAMBRIDGE INTERNATIONAL EXAMINATIONS International General Certificate of Seconary Eucation MARK SCHEME for the May/June 03 series 0607 CAMBRIDGE INTERNATIONAL MATHEMATICS 0607/4 Paper 4 (Extene), maximum

More information

0607 CAMBRIDGE INTERNATIONAL MATHEMATICS

0607 CAMBRIDGE INTERNATIONAL MATHEMATICS PAPA CAMBRIDGE CAMBRIDGE INTERNATIONAL EXAMINATIONS International General Certificate of Seconary Eucation MARK SCHEME for the May/June 0 series CAMBRIDGE INTERNATIONAL MATHEMATICS /4 4 (Extene), maximum

More information

Transient analysis of wave propagation in 3D soil by using the scaled boundary finite element method

Transient analysis of wave propagation in 3D soil by using the scaled boundary finite element method Southern Cross University epublications@scu 23r Australasian Conference on the Mechanics of Structures an Materials 214 Transient analysis of wave propagation in 3D soil by using the scale bounary finite

More information

New Geometric Interpretation and Analytic Solution for Quadrilateral Reconstruction

New Geometric Interpretation and Analytic Solution for Quadrilateral Reconstruction New Geometric Interpretation an Analytic Solution for uarilateral Reconstruction Joo-Haeng Lee Convergence Technology Research Lab ETRI Daejeon, 305 777, KOREA Abstract A new geometric framework, calle

More information

Lagrange multipliers October 2013

Lagrange multipliers October 2013 Lagrange multipliers 14.8 14 October 2013 Example: Optimization with constraint. Example: Find the extreme values of f (x, y) = x + 2y on the ellipse 3x 2 + 4y 2 = 3. 3/2 1 1 3/2 Example: Optimization

More information

CONSTRUCTION AND ANALYSIS OF INVERSIONS IN S 2 AND H 2. Arunima Ray. Final Paper, MATH 399. Spring 2008 ABSTRACT

CONSTRUCTION AND ANALYSIS OF INVERSIONS IN S 2 AND H 2. Arunima Ray. Final Paper, MATH 399. Spring 2008 ABSTRACT CONSTUCTION AN ANALYSIS OF INVESIONS IN S AN H Arunima ay Final Paper, MATH 399 Spring 008 ASTACT The construction use to otain inversions in two-imensional Eucliean space was moifie an applie to otain

More information

On a Method of Finding Homoclinic and Heteroclinic Orbits in. Multidimensional Dynamical Systems

On a Method of Finding Homoclinic and Heteroclinic Orbits in. Multidimensional Dynamical Systems Applie Mathematics & Information Sciences 4(3) (), 383 394 An International Journal c Diie W Publishing Corporation, U. S. A. On a Metho of Fining Homoclinic an Heteroclinic Orbits in Multiimensional Dynamical

More information

Differentiation Using Product and Quotient Rule 1

Differentiation Using Product and Quotient Rule 1 Differentiation Using Prouct an Quotient Rule 1 1.. ( + 1)( + + 1) + 1 + +. 4. (7 + 15) ( 7 + 15) 5. 6. ( + 7) (5 + 14) 7. 9. + 4 ( 1) (10 + ) 8 + 49 6 4 + 7 10. 8. 1 4 1 11. 1. 1 ( + 1) ( 1) 1. 04 + 59

More information

Physics 235 Chapter 6. Chapter 6 Some Methods in the Calculus of Variations

Physics 235 Chapter 6. Chapter 6 Some Methods in the Calculus of Variations Chapter 6 Some Methods in the Calculus of Variations In this Chapter we focus on an important method of solving certain problems in Classical Mechanics. In many problems we need to determine how a system

More information

Animação e Visualização Tridimensional. Collision Detection Corpo docente de AVT / CG&M / DEI / IST / UTL

Animação e Visualização Tridimensional. Collision Detection Corpo docente de AVT / CG&M / DEI / IST / UTL Animação e Visualização Triimensional Collision Detection Collision Hanling Collision Detection Collision Determination Collision Response Collision Hanling Collision Detection Collision Determination

More information

Implicit and Explicit Functions

Implicit and Explicit Functions 60_005.q //0 :5 PM Page SECTION.5 Implicit Differentiation Section.5 EXPLORATION Graphing an Implicit Equation How coul ou use a graphing utilit to sketch the graph of the equation? Here are two possible

More information

More on the Ray Matrix Formalism

More on the Ray Matrix Formalism More on the Ray Matrix Formalism Tuesay, 10/17/2006 Physics 158 Peter Beyersorf Document info 12. 1 Class Outline Paraxial Ray Matrices General Imaging Systems 12. 2 Compoun Optical Systems Compoun optical

More information

Lagrange multipliers 14.8

Lagrange multipliers 14.8 Lagrange multipliers 14.8 14 October 2013 Example: Optimization with constraint. Example: Find the extreme values of f (x, y) = x + 2y on the ellipse 3x 2 + 4y 2 = 3. 3/2 Maximum? 1 1 Minimum? 3/2 Idea:

More information

Interference and diffraction are the important phenomena that distinguish. Interference and Diffraction

Interference and diffraction are the important phenomena that distinguish. Interference and Diffraction C H A P T E R 33 Interference an Diffraction 33- Phase Difference an Coherence 33-2 Interference in Thin Films 33-3 Two-Slit Interference Pattern 33-4 Diffraction Pattern of a Single Slit * 33-5 Using

More information

Multiple Integrals. max x i 0

Multiple Integrals. max x i 0 Multiple Integrals 1 Double Integrals Definite integrals appear when one solves Area problem. Find the area A of the region bounded above by the curve y = f(x), below by the x-axis, and on the sides by

More information

Tracking and Regulation Control of a Mobile Robot System With Kinematic Disturbances: A Variable Structure-Like Approach

Tracking and Regulation Control of a Mobile Robot System With Kinematic Disturbances: A Variable Structure-Like Approach W. E. Dixon e-mail: wixon@ces.clemson.eu D. M. Dawson e-mail: awson@ces.clemson.eu E. Zergeroglu e-mail: ezerger@ces.clemson.eu Department of Electrical & Computer Engineering, Clemson University, Clemson,

More information

Attitude Dynamics and Control of a Dual-Body Spacecraft using Variable-Speed Control Moment Gyros

Attitude Dynamics and Control of a Dual-Body Spacecraft using Variable-Speed Control Moment Gyros Attitue Dynamics an Control of a Dual-Boy Spacecraft using Variable-Spee Control Moment Gyros Abstract Marcello Romano, Brij N. Agrawal Department of Mechanical an Astronautical Engineering, US Naval Postgrauate

More information

Grad operator, triple and line integrals. Notice: this material must not be used as a substitute for attending the lectures

Grad operator, triple and line integrals. Notice: this material must not be used as a substitute for attending the lectures Grad operator, triple and line integrals Notice: this material must not be used as a substitute for attending the lectures 1 .1 The grad operator Let f(x 1, x,..., x n ) be a function of the n variables

More information

SPECIAL TECHNIQUES-II

SPECIAL TECHNIQUES-II SPECIAL TECHNIQUES-II Lecture 19: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Method of Images for a spherical conductor Example :A dipole near aconducting sphere The

More information

Try It. Implicit and Explicit Functions. Video. Exploration A. Differentiating with Respect to x

Try It. Implicit and Explicit Functions. Video. Exploration A. Differentiating with Respect to x SECTION 5 Implicit Differentiation Section 5 Implicit Differentiation Distinguish between functions written in implicit form an eplicit form Use implicit ifferentiation to fin the erivative of a function

More information

Physics INTERFERENCE OF LIGHT

Physics INTERFERENCE OF LIGHT Physics INTERFERENCE OF LIGHT Q.1 State the principle of superposition of waves an explain the concept of interference of light. Ans. Principle of superposition of waves : When two or more waves, traveling

More information

Unit 2: Locomotion Kinematics of Wheeled Robots: Part 3

Unit 2: Locomotion Kinematics of Wheeled Robots: Part 3 Unit 2: Locomotion Kinematics of Wheeled Robots: Part 3 Computer Science 4766/6778 Department of Computer Science Memorial University of Newfoundland January 28, 2014 COMP 4766/6778 (MUN) Kinematics of

More information

(c) 0 (d) (a) 27 (b) (e) x 2 3x2

(c) 0 (d) (a) 27 (b) (e) x 2 3x2 1. Sarah the architect is designing a modern building. The base of the building is the region in the xy-plane bounded by x =, y =, and y = 3 x. The building itself has a height bounded between z = and

More information

Advanced method of NC programming for 5-axis machining

Advanced method of NC programming for 5-axis machining Available online at www.scienceirect.com Proceia CIRP (0 ) 0 07 5 th CIRP Conference on High Performance Cutting 0 Avance metho of NC programming for 5-axis machining Sergej N. Grigoriev a, A.A. Kutin

More information

Local Path Planning with Proximity Sensing for Robot Arm Manipulators. 1. Introduction

Local Path Planning with Proximity Sensing for Robot Arm Manipulators. 1. Introduction Local Path Planning with Proximity Sensing for Robot Arm Manipulators Ewar Cheung an Vlaimir Lumelsky Yale University, Center for Systems Science Department of Electrical Engineering New Haven, Connecticut

More information

Algebraic transformations of Gauss hypergeometric functions

Algebraic transformations of Gauss hypergeometric functions Algebraic transformations of Gauss hypergeometric functions Raimunas Viūnas Faculty of Mathematics, Kobe University Abstract This article gives a classification scheme of algebraic transformations of Gauss

More information

Particle Swarm Optimization Based on Smoothing Approach for Solving a Class of Bi-Level Multiobjective Programming Problem

Particle Swarm Optimization Based on Smoothing Approach for Solving a Class of Bi-Level Multiobjective Programming Problem BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 17, No 3 Sofia 017 Print ISSN: 1311-970; Online ISSN: 1314-4081 DOI: 10.1515/cait-017-0030 Particle Swarm Optimization Base

More information

Lagrange Multipliers. Lagrange Multipliers. Lagrange Multipliers. Lagrange Multipliers. Lagrange Multipliers. Lagrange Multipliers

Lagrange Multipliers. Lagrange Multipliers. Lagrange Multipliers. Lagrange Multipliers. Lagrange Multipliers. Lagrange Multipliers In this section we present Lagrange s method for maximizing or minimizing a general function f(x, y, z) subject to a constraint (or side condition) of the form g(x, y, z) = k. Figure 1 shows this curve

More information

MATH203 Calculus. Dr. Bandar Al-Mohsin. School of Mathematics, KSU

MATH203 Calculus. Dr. Bandar Al-Mohsin. School of Mathematics, KSU School of Mathematics, KSU Theorem The rectangular coordinates (x, y, z) and the cylindrical coordinates (r, θ, z) of a point P are related as follows: x = r cos θ, y = r sin θ, tan θ = y x, r 2 = x 2

More information

To graph the point (r, θ), simply go out r units along the initial ray, then rotate through the angle θ. The point (1, 5π 6

To graph the point (r, θ), simply go out r units along the initial ray, then rotate through the angle θ. The point (1, 5π 6 Polar Coordinates Any point in the plane can be described by the Cartesian coordinates (x, y), where x and y are measured along the corresponding axes. However, this is not the only way to represent points

More information

Refinement of scene depth from stereo camera ego-motion parameters

Refinement of scene depth from stereo camera ego-motion parameters Refinement of scene epth from stereo camera ego-motion parameters Piotr Skulimowski, Pawel Strumillo An algorithm for refinement of isparity (epth) map from stereoscopic sequences is propose. The metho

More information

Kinematics of Wheeled Robots

Kinematics of Wheeled Robots CSE 390/MEAM 40 Kinematics of Wheeled Robots Professor Vijay Kumar Department of Mechanical Engineering and Applied Mechanics University of Pennsylvania September 16, 006 1 Introduction In this chapter,

More information

Computer Graphics Chapter 7 Three-Dimensional Viewing Viewing

Computer Graphics Chapter 7 Three-Dimensional Viewing Viewing Computer Graphics Chapter 7 Three-Dimensional Viewing Outline Overview of Three-Dimensional Viewing Concepts The Three-Dimensional Viewing Pipeline Three-Dimensional Viewing-Coorinate Parameters Transformation

More information

Winter 2012 Math 255 Section 006. Problem Set 7

Winter 2012 Math 255 Section 006. Problem Set 7 Problem Set 7 1 a) Carry out the partials with respect to t and x, substitute and check b) Use separation of varibles, i.e. write as dx/x 2 = dt, integrate both sides and observe that the solution also

More information

3. The three points (2, 4, 1), (1, 2, 2) and (5, 2, 2) determine a plane. Which of the following points is in that plane?

3. The three points (2, 4, 1), (1, 2, 2) and (5, 2, 2) determine a plane. Which of the following points is in that plane? Math 4 Practice Problems for Midterm. A unit vector that is perpendicular to both V =, 3, and W = 4,, is (a) V W (b) V W (c) 5 6 V W (d) 3 6 V W (e) 7 6 V W. In three dimensions, the graph of the equation

More information

Goals: Course Unit: Describing Moving Objects Different Ways of Representing Functions Vector-valued Functions, or Parametric Curves

Goals: Course Unit: Describing Moving Objects Different Ways of Representing Functions Vector-valued Functions, or Parametric Curves Block #1: Vector-Valued Functions Goals: Course Unit: Describing Moving Objects Different Ways of Representing Functions Vector-valued Functions, or Parametric Curves 1 The Calculus of Moving Objects Problem.

More information

6 Gradient Descent. 6.1 Functions

6 Gradient Descent. 6.1 Functions 6 Graient Descent In this topic we will iscuss optimizing over general functions f. Typically the function is efine f : R! R; that is its omain is multi-imensional (in this case -imensional) an output

More information

MATH2111 Higher Several Variable Calculus Lagrange Multipliers

MATH2111 Higher Several Variable Calculus Lagrange Multipliers MATH2111 Higher Several Variable Calculus Lagrange Multipliers Dr. Jonathan Kress School of Mathematics and Statistics University of New South Wales Semester 1, 2016 [updated: February 29, 2016] JM Kress

More information

A new fuzzy visual servoing with application to robot manipulator

A new fuzzy visual servoing with application to robot manipulator 2005 American Control Conference June 8-10, 2005. Portlan, OR, USA FrA09.4 A new fuzzy visual servoing with application to robot manipulator Marco A. Moreno-Armenariz, Wen Yu Abstract Many stereo vision

More information

Math 265 Exam 3 Solutions

Math 265 Exam 3 Solutions C Roettger, Fall 16 Math 265 Exam 3 Solutions Problem 1 Let D be the region inside the circle r 5 sin θ but outside the cardioid r 2 + sin θ. Find the area of D. Note that r and θ denote polar coordinates.

More information

Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.

Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed. Preface Here are my online notes for my Calculus I course that I teach here at Lamar University. Despite the fact that these are my class notes, they shoul be accessible to anyone wanting to learn Calculus

More information

Vision-based Multi-Robot Simultaneous Localization and Mapping

Vision-based Multi-Robot Simultaneous Localization and Mapping Vision-base Multi-Robot Simultaneous Localization an Mapping Hassan Hajjiab an Robert Laganière VIVA Research Lab School of Information Technology an Engineering University of Ottawa, Ottawa, Canaa, K1N

More information

Homogeneous coordinates, lines, screws and twists

Homogeneous coordinates, lines, screws and twists Homogeneous coordinates, lines, screws and twists In lecture 1 of module 2, a brief mention was made of homogeneous coordinates, lines in R 3, screws and twists to describe the general motion of a rigid

More information

Q. No. 1 Newton postulated his corpuscular theory of light on the basis of

Q. No. 1 Newton postulated his corpuscular theory of light on the basis of Q. No. 1 Newton postulate his corpuscular theory of light on the basis of Newton s rings Option Rectilinear propagation of light Colour through thin films Dispersion of white light into colours. Correct

More information

Holy Halved Heaquarters Riddler

Holy Halved Heaquarters Riddler Holy Halve Heaquarters Riler Anonymous Philosopher June 206 Laser Larry threatens to imminently zap Riler Heaquarters (which is of regular pentagonal shape with no courtyar or other funny business) with

More information

SECTION 2. Objectives. Identify appropriate coordinate systems for solving problems with vectors.

SECTION 2. Objectives. Identify appropriate coordinate systems for solving problems with vectors. SECTION 2 Plan an Prepare Preview Vocabular Scientific Meaning The wor simultaneous is use for phenomena that occur together at the same time. Ask stuents to list some simultaneous phenomena, such as the

More information

A Versatile Model-Based Visibility Measure for Geometric Primitives

A Versatile Model-Based Visibility Measure for Geometric Primitives A Versatile Moel-Base Visibility Measure for Geometric Primitives Marc M. Ellenrieer 1,LarsKrüger 1, Dirk Stößel 2, an Marc Hanheie 2 1 DaimlerChrysler AG, Research & Technology, 89013 Ulm, Germany 2 Faculty

More information

BIJECTIONS FOR PLANAR MAPS WITH BOUNDARIES

BIJECTIONS FOR PLANAR MAPS WITH BOUNDARIES BIJECTIONS FOR PLANAR MAPS WITH BOUNDARIES OLIVIER BERNARDI AND ÉRIC FUSY Abstract. We present bijections for planar maps with bounaries. In particular, we obtain bijections for triangulations an quarangulations

More information

State Indexed Policy Search by Dynamic Programming. Abstract. 1. Introduction. 2. System parameterization. Charles DuHadway

State Indexed Policy Search by Dynamic Programming. Abstract. 1. Introduction. 2. System parameterization. Charles DuHadway State Inexe Policy Search by Dynamic Programming Charles DuHaway Yi Gu 5435537 503372 December 4, 2007 Abstract We consier the reinforcement learning problem of simultaneous trajectory-following an obstacle

More information

To find the maximum and minimum values of f(x, y, z) subject to the constraints

To find the maximum and minimum values of f(x, y, z) subject to the constraints Midterm 3 review Math 265 Fall 2007 14.8. Lagrange Multipliers. Case 1: One constraint. To find the maximum and minimum values of f(x, y, z) subject to the constraint g(x, y, z) = k: Step 1: Find all values

More information

CMPUT 412 Motion Control Wheeled robots. Csaba Szepesvári University of Alberta

CMPUT 412 Motion Control Wheeled robots. Csaba Szepesvári University of Alberta CMPUT 412 Motion Control Wheeled robots Csaba Szepesvári University of Alberta 1 Motion Control (wheeled robots) Requirements Kinematic/dynamic model of the robot Model of the interaction between the wheel

More information

Fundamentals of the stiffness method. Fundamentals of the stiffness method. Fundamentals of the stiffness method

Fundamentals of the stiffness method. Fundamentals of the stiffness method. Fundamentals of the stiffness method CHAPER 6 russ Analsis using Stiffness Metho Objectives เข าใจว ธ ของ stiffness metho ใช ว ธ stiffness metho ก บ russ, BM & rame จะพ ดถ งในบท หน า unamentals of the stiffness metho he stiffness metho: Is

More information

COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates

COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates Department of Computer Science and Software Engineering The Lecture outline Introduction Vectors and matrices Translation

More information

The Reconstruction of Graphs. Dhananjay P. Mehendale Sir Parashurambhau College, Tilak Road, Pune , India. Abstract

The Reconstruction of Graphs. Dhananjay P. Mehendale Sir Parashurambhau College, Tilak Road, Pune , India. Abstract The Reconstruction of Graphs Dhananay P. Mehenale Sir Parashurambhau College, Tila Roa, Pune-4030, Inia. Abstract In this paper we iscuss reconstruction problems for graphs. We evelop some new ieas lie

More information

Classifying Facial Expression with Radial Basis Function Networks, using Gradient Descent and K-means

Classifying Facial Expression with Radial Basis Function Networks, using Gradient Descent and K-means Classifying Facial Expression with Raial Basis Function Networks, using Graient Descent an K-means Neil Allrin Department of Computer Science University of California, San Diego La Jolla, CA 9237 nallrin@cs.ucs.eu

More information

14.5 Directional Derivatives and the Gradient Vector

14.5 Directional Derivatives and the Gradient Vector 14.5 Directional Derivatives and the Gradient Vector 1. Directional Derivatives. Recall z = f (x, y) and the partial derivatives f x and f y are defined as f (x 0 + h, y 0 ) f (x 0, y 0 ) f x (x 0, y 0

More information

AN INVESTIGATION OF FOCUSING AND ANGULAR TECHNIQUES FOR VOLUMETRIC IMAGES BY USING THE 2D CIRCULAR ULTRASONIC PHASED ARRAY

AN INVESTIGATION OF FOCUSING AND ANGULAR TECHNIQUES FOR VOLUMETRIC IMAGES BY USING THE 2D CIRCULAR ULTRASONIC PHASED ARRAY AN INVESTIGATION OF FOCUSING AND ANGULAR TECHNIQUES FOR VOLUMETRIC IMAGES BY USING THE D CIRCULAR ULTRASONIC PHASED ARRAY S. Monal Lonon South Bank University; Engineering an Design 103 Borough Roa, Lonon

More information

PLANCESS RANK ACCELERATOR

PLANCESS RANK ACCELERATOR PLANCESS RANK ACCELERATOR PHYSICS FOR JEE MAIN & ADVANCED Wave Optics 4000+questions with topic wise exercises 000+ problems of IIT-JEE & AIEEE exams of last 5 years 4 Levels of Exercises categorize into

More information

SURVIVABLE IP OVER WDM: GUARANTEEEING MINIMUM NETWORK BANDWIDTH

SURVIVABLE IP OVER WDM: GUARANTEEEING MINIMUM NETWORK BANDWIDTH SURVIVABLE IP OVER WDM: GUARANTEEEING MINIMUM NETWORK BANDWIDTH Galen H Sasaki Dept Elec Engg, U Hawaii 2540 Dole Street Honolul HI 96822 USA Ching-Fong Su Fuitsu Laboratories of America 595 Lawrence Expressway

More information

Dual Arm Robot Research Report

Dual Arm Robot Research Report Dual Arm Robot Research Report Analytical Inverse Kinematics Solution for Moularize Dual-Arm Robot With offset at shouler an wrist Motivation an Abstract Generally, an inustrial manipulator such as PUMA

More information

f xx (x, y) = 6 + 6x f xy (x, y) = 0 f yy (x, y) = y In general, the quantity that we re interested in is

f xx (x, y) = 6 + 6x f xy (x, y) = 0 f yy (x, y) = y In general, the quantity that we re interested in is 1. Let f(x, y) = 5 + 3x 2 + 3y 2 + 2y 3 + x 3. (a) Final all critical points of f. (b) Use the second derivatives test to classify the critical points you found in (a) as a local maximum, local minimum,

More information

A Classification of 3R Orthogonal Manipulators by the Topology of their Workspace

A Classification of 3R Orthogonal Manipulators by the Topology of their Workspace A Classification of R Orthogonal Manipulators by the Topology of their Workspace Maher aili, Philippe Wenger an Damien Chablat Institut e Recherche en Communications et Cybernétique e Nantes, UMR C.N.R.S.

More information

Applications. Human and animal motion Robotics control Hair Plants Molecular motion

Applications. Human and animal motion Robotics control Hair Plants Molecular motion Multibody dynamics Applications Human and animal motion Robotics control Hair Plants Molecular motion Generalized coordinates Virtual work and generalized forces Lagrangian dynamics for mass points

More information

Robotics (Kinematics) Winter 1393 Bonab University

Robotics (Kinematics) Winter 1393 Bonab University Robotics () Winter 1393 Bonab University : most basic study of how mechanical systems behave Introduction Need to understand the mechanical behavior for: Design Control Both: Manipulators, Mobile Robots

More information

Kinematic Analysis of a Family of 3R Manipulators

Kinematic Analysis of a Family of 3R Manipulators Kinematic Analysis of a Family of R Manipulators Maher Baili, Philippe Wenger an Damien Chablat Institut e Recherche en Communications et Cybernétique e Nantes, UMR C.N.R.S. 6597 1, rue e la Noë, BP 92101,

More information

Chapter 11. Parametric Equations And Polar Coordinates

Chapter 11. Parametric Equations And Polar Coordinates Instructor: Prof. Dr. Ayman H. Sakka Chapter 11 Parametric Equations And Polar Coordinates In this chapter we study new ways to define curves in the plane, give geometric definitions of parabolas, ellipses,

More information

Regularized Laplacian Zero Crossings as Optimal Edge Integrators

Regularized Laplacian Zero Crossings as Optimal Edge Integrators Regularize aplacian Zero rossings as Optimal Ege Integrators R. KIMME A.M. BRUKSTEIN Department of omputer Science Technion Israel Institute of Technology Technion ity, Haifa 32, Israel Abstract We view

More information

Unit #5 - Implicit Differentiation, Related Rates Section 3.7

Unit #5 - Implicit Differentiation, Related Rates Section 3.7 Unit #5 - Implicit Differentiation, Relate Rates Section 3.7 Some material from Calculus, Single an MultiVariable b Hughes-Hallett, Gleason, McCallum et. al. Copright 005 b John Wile & Sons, Inc. This

More information

X y. f(x,y,d) f(x,y,d) Peak. Motion stereo space. parameter space. (x,y,d) Motion stereo space. Parameter space. Motion stereo space.

X y. f(x,y,d) f(x,y,d) Peak. Motion stereo space. parameter space. (x,y,d) Motion stereo space. Parameter space. Motion stereo space. 3D Shape Measurement of Unerwater Objects Using Motion Stereo Hieo SAITO Hirofumi KAWAMURA Masato NAKAJIMA Department of Electrical Engineering, Keio Universit 3-14-1Hioshi Kouhoku-ku Yokohama 223, Japan

More information

1. Suppose that the equation F (x, y, z) = 0 implicitly defines each of the three variables x, y, and z as functions of the other two:

1. Suppose that the equation F (x, y, z) = 0 implicitly defines each of the three variables x, y, and z as functions of the other two: Final Solutions. Suppose that the equation F (x, y, z) implicitly defines each of the three variables x, y, and z as functions of the other two: z f(x, y), y g(x, z), x h(y, z). If F is differentiable

More information

Graphics and Interaction Transformation geometry and homogeneous coordinates

Graphics and Interaction Transformation geometry and homogeneous coordinates 433-324 Graphics and Interaction Transformation geometry and homogeneous coordinates Department of Computer Science and Software Engineering The Lecture outline Introduction Vectors and matrices Translation

More information

Loop Scheduling and Partitions for Hiding Memory Latencies

Loop Scheduling and Partitions for Hiding Memory Latencies Loop Scheuling an Partitions for Hiing Memory Latencies Fei Chen Ewin Hsing-Mean Sha Dept. of Computer Science an Engineering University of Notre Dame Notre Dame, IN 46556 Email: fchen,esha @cse.n.eu Tel:

More information

LinkA LinkB LinkC (X,Y) Θ

LinkA LinkB LinkC (X,Y) Θ Following Performance Evaluation of a Mechanicall Couple Platoon Tetsua SHINOHARA, Masafumi OKADA an Ken-ichi YOSHIMOTO Department ofmechano-informatics, Universit oftoko 7-3-, Hongo, Bunko-ku, Toko, 3-8656

More information