Optimal Control Problems in Eye Movement. Bijoy K. Ghosh Department of Mathematics and Statistics Texas Tech University, Lubbock, TX
|
|
- Rudolph Cannon
- 5 years ago
- Views:
Transcription
1 Optimal Control Problems in Eye Movement Bijoy K. Ghosh Department of Mathematics and Statistics Texas Tech University, Lubbock, TX
2
3 We study the human oculomotor system as a simple mechanical control system. It is a well known physiological fact that all eye movements obey Listing s Law which states that eye rotations are always restricted to a given fixed plane called the Listing s plane.
4 In fact, Listing s constraint is sufficient to ensure that corresponding to any specific gaze direction, the torsion is predetermined. In other words, orientation of the eye is completely specified by its gaze. Thus objects in space are not viewed with a rotational ambiguity.
5 All these ideas are very interesting and biomechanists are interested in the study of how the brain is able to satisfy the Listing s constraint through muscle based actuation of the eye ball. The subject matter of my talk is the following: Assume that the Listing s Law is satisfied, how can we rotate the eye ball optimally from one gaze direction to
6 It is natural to start with SO (3), the space of 3 by 3 rotation matrices, as the configuration space of the eye ball. We denote by List the subset of which is described as follows: -- SO (3) List = R SO(3) v =[v 1,v 2 ] R 2 \{0} 3R v 1 v 2 0 = v 1 v 2 0
7 It is natural to start with SO (3), the space of 3 by 3 rotation matrices, as the configuration space of the eye ball. We denote by List the subset of which is described as follows: -- SO (3) The space List is the set of all 3 by 3 rotation matrices that has axis of rotation perpendicular to the vector
8 It turns out that elements of List are rotation matrices of the form : q2 0 + q1 2 q2 2 q3 2 2(q 1 q 2 q 0 q 3 ) 2(q 1 q 3 + q 0 q 2 ) 2(q 1 q 2 + q 0 q 3 ) q0 2 + q2 2 q1 2 q3 2 2(q 2 q 3 q 0 q 1 ) 2(q 1 q 3 q 0 q 2 ) 2(q 2 q 3 + q 0 q 1 ) q0 2 + q3 2 q1 2 q2 2 where q 3 =0, and (q 0,q 1,q 2 ) is a unit vector. It is not at all clear, apriori, why this would be the case.
9 It turns out that there is a correspondence between rotation matrices and unit quaternions. A quaternion Q is a 4 tuple of real numbers written as a 0 ~1 + a 1 ~ i + a2 ~ j + a3 ~ k and we have two obvious maps: vec : Q R 3, a 7 (a 1,a 2,a 3 ) called the vector part of the quaternion.
10 It turns out that there is a correspondence between rotation matrices and unit quaternions. A quaternion Q is a 4 tuple of real numbers written as a 0 ~1 + a 1 ~ i + a2 ~ j + a3 ~ k and we have two obvious maps: scal : Q R, a 7 a 0 called the scalar part of the quaternion.
11 The space of unit quaternions will be identified with the unit sphere in R 4, and is denoted by S 3. Each q S 3 can be written as q =cos(α/2) ~1 +sin(α/2) n 1 ~ i +sin(α/2) n2 ~ j +sin(α/2) n3 ~ k where α [0, π] and n =(n 1,n 2,n 3 ) is a unit vector in R 3. We denote by rot the standard map from S 3 into SO 3, which maps q to a rotation around the axis n by a clockwise angle α. There are two explicit ways to describe this map. First it is easy to verify that, rot (q)(v 1,v 2,v 3 )=vec ³ q ³v 1 ~ i + v2 ~ j + v2 ~ j 1 q
12 Alternatively, we can think of rot(q) as follows: rot(q) = q0 2 + q1 2 q2 2 q 2 3 2(q 1 q 2 q 0 q 3 ) 2(q 1 q 3 + q 0 q 2 ) 2(q 1 q 2 + q 0 q 3 ) q0 2 + q2 2 q1 2 q3 2 2(q 2 q 3 q 0 q 1 ) 2(q 1 q 3 q 0 q 2 ) 2(q 2 q 3 + q 0 q 1 ) q0 2 + q3 2 q1 2 q2 2
13 The main point is that counterclockwise rotation by an angle α with respect to an axis n =(n 1,n 2,n 3 ), an unit vector, is given precisely by the unit quaternion cos (α/2) ~1+sin (α/2) n 1 ~ i+sin (α/2) n2 ~ j+sin (α/2) n3 ~ k where Listing s constraint is imposed by restricting n 3 =0
14 Alternative we can say that unit quaternions with the fourth coordinate 0 precisely corresponds to rotation matrices that satisfy the Listing s constraint. Thus we have a map: ρ :[0, π] [0, 2π] S 3 cos (φ/2) ρ (θ, φ) = cos (θ)sin(φ/2) sin (θ)sin(φ/2) 0
15 Riemannian metric on List We assume that eye is a perfect sphere and its inertia tensor is equal to. I 3 3 This is associated with the left invariant Riemannian metric on given by SO 3 hω (e i ), Ω (e j )i I = δ i,j Where Ω (e k )= 0 δ 3,k δ 2,k δ 3,k 0 δ 1,k δ 2,k δ 1,k 0 and {δ i,j } denotes the Kronecker delta function
16 To compute the Riemannian metric on List induced from, we define SO 3 g 11 = h θ, θ i g 12 = h θ, φ i g 22 = h φ, φ i
17 Using the map ρ :[0, π] [0, 2π] S 3 cos (φ/2) ρ (θ, φ) = cos (θ)sin(φ/2) sin (θ)sin(φ/2) 0 it follows that we have the following diagram: List T (θ,φ) List ρ S 3 ρ Tρ(θ,φ) S 3
18 We can compute the following: µ µ ρ θ µ ρ φ = sin (φ/2) 1 sin (θ)sin(φ/2) 2 cos (θ)cos(φ/2) 1 cos (θ)sin(φ/2) 2 sin (θ)cos(φ/2) 0 0
19 The Riemannian metric on List is easily computed by taking inner products and we obtain the following: g 11 =sin 2 (φ/2) g 12 =0 g 22 = 1 4 Finally we obtain g =sin 2 (φ/2) dθ dφ2
20 Using a Riemannian metric one is able to measure distance between points on a manifold. It is a standard calculation in differential geometry to treat the Riemannian Metric as the kinetic energy of the eye in motion and compute the path taken by the eye that would minimize this energy, while disregarding the potential energy. This path would be called the Geodesic
21 Equations of geodesics are given in the form of a differential equation. For our problem it is given by θ + 1 tan (φ/2) θ φ =0 φ sin (φ) φ 2 =0 θ Recall that is the angle for the axis of rotation with respect to the horizontal line on the Listing s plane. φ on the other hand is the angle of rotation with respect to the given axis.
22 If we define state variables as x 1 = θ ; x 3 = φ ; x 2 = θ x 4 = φ then the geodesic is given by: ẋ 1 = x 2 ẋ 2 = x 2 x 4 cot (x 3 /2) ẋ 3 = x 4 ẋ 4 = x 2 2 sin (x 3 ) This dynamical system has a singularity when is an even multiple of φ and equilibrium when is an odd multiple of and is zero. φ π φ π
23 The space of (θ, φ) isactuallythe surface of a doughnut in R 3. The geodesic is defined as a dynamical system on this doughnut. This dynamical system has a singularity when is an even multiple of φ and equilibrium when is an odd multiple of and is zero. φ π φ π We would like to ask for example How to deal with this singularity? Is the equilibrium stable, unstable or saddle?
24 It turns out that to deal with the singularity we can define state variables as x 1 = θ ; x 3 = φ ; and the geodesic is given by: ẋ 1 = x 2 tan (x 3 /2) ẋ 2 = x 2 x cot 2 (x 3 /2) ẋ 3 = x 4 tan (x 3 /2) x 2 = θ/ tan (φ/2) x 4 = φ/ tan (φ/2) ẋ 4 = x 2 2 sin (x 3 )tan(x 3 /2).5 x 4 2 sec 2 (x 3 /2) This dynamical system has a singularity when is an odd multiple of φ and equilibrium when is an even multiple of and is zero. φ π x 4 π
25 Thus there is no way to get out of singularity but one can define the geodesic on open charts that cover the whole space and switch between charts.
26 So far we have not considered any potential energy and also did not consider any external torque as a control input. If we consider the potential energy given by V (θ, φ) =A cos 2 (φ/2) The geodesic equations are modified as follows: θ + θ φ cot (φ/2) = csc 2 (φ/2) τ θ φ φ 2 sin (φ)+ 1 2 sin (φ) 2A =4τ φ
27 We can now wish to control the state from (θ 0, 0, φ 0, 0) to (θ 1, 0, φ 1, 0) in T units of time, while minimizing the control energy R T 0 h (τ θ (t)) 2 +(τ φ (t)) 2i dt We use the Hamiltonian H (z,λ) =λ ż 1 2 τ 2 θ + τ 2 φ and the Hamilton s equations: ż = H λ, λ = H z
28 We obtain a state space system given by: ż 1 ż 2 ż 3 ż 4 λ 1 λ 2 λ 3 λ 4 = z 2 z 2 z 4 cot (z 3 /2) + λ 2 csc 4 (z 3 /2) z 4 z 2 sin (z 3 )+2Acos (z 3 )+16λ 4 0 λ 1 + λ 2 z 4 cot (z 3 /2) 2λ 4 z 2 sin (z 3 ) λ 2 z 2 cot (z 3 /2) + λ 3 where : = 1 2 λ 2z 2 z 4 csc 2 (z 3 /2) + λ 2 2 cot (z 3/2) csc 4 (z 3 /2) λ 4 z λ 4A cos (z 3 )
29 The optimal controls are given by: τ θ = λ 2 sin 2 (z 3 /2) τ φ =4λ 4
30 An example of an optimal and a geodesic path
31 θ and φ as a function of time for the optimal path
32 The trajectories of the corresponding optimal torques
33 x0 = [1.5,.1e-3,3,.0001] V=0,-1,-5 e-5
34 x0 = [1.5,.1e-3,3,.001] V=0, -.1, -1, -10, -100 e-5
35 V=0
36 V=.1
37 V=1
38 V=10
39 V=100
40 V=1000
41
Physics 411 Lecture 7. Tensors. Lecture 7. Physics 411 Classical Mechanics II
Physics 411 Lecture 7 Tensors Lecture 7 Physics 411 Classical Mechanics II September 12th 2007 In Electrodynamics, the implicit law governing the motion of particles is F α = m ẍ α. This is also true,
More informationTutorial General Relativity
Tutorial General Relativity Winter term 016/017 Sheet No. 3 Solutions will be discussed on Nov/9/16 Lecturer: Prof. Dr. C. Greiner Tutor: Hendrik van Hees 1. Tensor gymnastics (a) Let Q ab = Q ba be a
More information1 Vector fields, flows and Lie derivatives
Working title: Notes on Lie derivatives and Killing vector fields Author: T. Harmark Killingvectors.tex 5//2008, 22:55 Vector fields, flows and Lie derivatives Coordinates Consider an n-dimensional manifold
More informationAppendix to Lecture 2
PHYS 652: Astrophysics 1 Appendix to Lecture 2 An Alternative Lagrangian In class we used an alternative Lagrangian L = g γδ ẋ γ ẋ δ, instead of the traditional L = g γδ ẋ γ ẋ δ. Here is the justification
More informationMath review. Math review
Math review 1 Math review 3 1 series approximations 3 Taylor s Theorem 3 Binomial approximation 3 sin(x), for x in radians and x close to zero 4 cos(x), for x in radians and x close to zero 5 2 some geometry
More informationA Lie Group. These notes introduce SU(2) as an example of a compact Lie group. SU(2) = { A A a 2 2 complex matrix,deta = 1,AA = A A = 1l }
A Lie Group The Definition These notes introduce SU2) as an example of a compact Lie group. The definition of SU2) is SU2) { A A 2 complex matrixdeta 1AA A A 1l } In the name SU2) the S stands for special
More informationRotational & Rigid-Body Mechanics. Lectures 3+4
Rotational & Rigid-Body Mechanics Lectures 3+4 Rotational Motion So far: point objects moving through a trajectory. Next: moving actual dimensional objects and rotating them. 2 Circular Motion - Definitions
More informationEuler Characteristic of Two-Dimensional Manifolds
Euler Characteristic of Two-Dimensional Manifolds M. Hafiz Khusyairi August 2008 In this work we will discuss an important notion from topology, namely Euler Characteristic and we will discuss several
More informationRotational Kinetic Energy
Lecture 17, Chapter 10: Rotational Energy and Angular Momentum 1 Rotational Kinetic Energy Consider a rigid body rotating with an angular velocity ω about an axis. Clearly every point in the rigid body
More informationÜbungen zu RT2 SS (4) Show that (any) contraction of a (p, q) - tensor results in a (p 1, q 1) - tensor.
Übungen zu RT2 SS 2010 (1) Show that the tensor field g µν (x) = η µν is invariant under Poincaré transformations, i.e. x µ x µ = L µ νx ν + c µ, where L µ ν is a constant matrix subject to L µ ρl ν ση
More information2.10 Saddles, Nodes, Foci and Centers
2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one
More informationAssignment 2. Goldstein 2.3 Prove that the shortest distance between two points in space is a straight line.
Assignment Goldstein.3 Prove that the shortest distance between two points in space is a straight line. The distance between two points is given by the integral of the infinitesimal arclength: s = = =
More informationExact Solutions of the Einstein Equations
Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.
More informationTHE COMPOUND ANGLE IDENTITIES
TRIGONOMETRY THE COMPOUND ANGLE IDENTITIES Question 1 Prove the validity of each of the following trigonometric identities. a) sin x + cos x 4 4 b) cos x + + 3 sin x + 2cos x 3 3 c) cos 2x + + cos 2x cos
More information5 Trigonometric Functions
5 Trigonometric Functions 5.1 The Unit Circle Definition 5.1 The unit circle is the circle of radius 1 centered at the origin in the xyplane: x + y = 1 Example: The point P Terminal Points (, 6 ) is on
More informationHandout 7: Torque, angular momentum, rotational kinetic energy and rolling motion. Torque and angular momentum
Handout 7: Torque, angular momentum, rotational kinetic energy and rolling motion Torque and angular momentum In Figure, in order to turn a rod about a fixed hinge at one end, a force F is applied at a
More informationChristoffel Symbols. 1 In General Topologies. Joshua Albert. September 28, W. First we say W : λ n = x µ (λ) so that the world
Christoffel Symbols Joshua Albert September 28, 22 In General Topoloies We have a metric tensor nm defined by, Note by some handy theorem that for almost any continuous function F (L), equation 2 still
More informationLecture 41: Highlights
Lecture 41: Highlights The goal of this lecture is to remind you of some of the key points that we ve covered this semester Note that this is not the complete set of topics that may appear on the final
More informationClassical Mechanics Comprehensive Exam Solution
Classical Mechanics Comprehensive Exam Solution January 31, 011, 1:00 pm 5:pm Solve the following six problems. In the following problems, e x, e y, and e z are unit vectors in the x, y, and z directions,
More informationIntroduction and Vectors Lecture 1
1 Introduction Introduction and Vectors Lecture 1 This is a course on classical Electromagnetism. It is the foundation for more advanced courses in modern physics. All physics of the modern era, from quantum
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationMath Section 4.3 Unit Circle Trigonometry
Math 10 - Section 4. Unit Circle Trigonometry An angle is in standard position if its vertex is at the origin and its initial side is along the positive x axis. Positive angles are measured counterclockwise
More information16. Elastic Scattering Michael Fowler
6 Elastic Scattering Michael Fowler Billiard Balls Elastic means no internal energy modes of the scatterer or of the scatteree are excited so total kinetic energy is conserved As a simple first exercise,
More informationPhysics 351, Spring 2015, Final Exam.
Physics 351, Spring 2015, Final Exam. This closed-book exam has (only) 25% weight in your course grade. You can use one sheet of your own hand-written notes. Please show your work on these pages. The back
More informationPart II. Classical Dynamics. Year
Part II Year 28 27 26 25 24 23 22 21 20 2009 2008 2007 2006 2005 28 Paper 1, Section I 8B Derive Hamilton s equations from an action principle. 22 Consider a two-dimensional phase space with the Hamiltonian
More informationThe Higgins-Selkov oscillator
The Higgins-Selkov oscillator May 14, 2014 Here I analyse the long-time behaviour of the Higgins-Selkov oscillator. The system is ẋ = k 0 k 1 xy 2, (1 ẏ = k 1 xy 2 k 2 y. (2 The unknowns x and y, being
More informationMath 223 Final. July 24, 2014
Math 223 Final July 24, 2014 Name Directions: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 Total 1. No books, notes, or evil looks. You may use a calculator to do routine arithmetic computations. You may not use your
More informationPhysics 411 Lecture 8. Parametrized Motion. Lecture 8. Physics 411 Classical Mechanics II
Physics 411 Lecture 8 Parametrized Motion Lecture 8 Physics 411 Classical Mechanics II September 14th 2007 We have our fancy new derivative, but what to do with it? In particular, how can we interpret
More information( 3 ) = (r) cos (390 ) =
MATH 7A Test 4 SAMPLE This test is in two parts. On part one, you may not use a calculator; on part two, a (non-graphing) calculator is necessary. When you complete part one, you turn it in and get part
More informationSymmetries 2 - Rotations in Space
Symmetries 2 - Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system
More informationMassachusetts Institute of Technology Department of Physics. Final Examination December 17, 2004
Massachusetts Institute of Technology Department of Physics Course: 8.09 Classical Mechanics Term: Fall 004 Final Examination December 17, 004 Instructions Do not start until you are told to do so. Solve
More informationSince 1 revolution = 1 = = Since 1 revolution = 1 = =
Fry Texas A&M University Math 150 Chapter 8A Fall 2015! 207 Since 1 revolution = 1 = = Since 1 revolution = 1 = = Convert to revolutions (or back to degrees and/or radians) a) 45! = b) 120! = c) 450! =
More informationUsing this definition, it is possible to define an angle of any (positive or negative) measurement by recognizing how its terminal side is obtained.
Angle in Standard Position With the Cartesian plane, we define an angle in Standard Position if it has its vertex on the origin and one of its sides ( called the initial side ) is always on the positive
More informationCP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017
CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and
More informationFind the indicated derivative. 1) Find y(4) if y = 3 sin x. A) y(4) = 3 cos x B) y(4) = 3 sin x C) y(4) = - 3 cos x D) y(4) = - 3 sin x
Assignment 5 Name Find the indicated derivative. ) Find y(4) if y = sin x. ) A) y(4) = cos x B) y(4) = sin x y(4) = - cos x y(4) = - sin x ) y = (csc x + cot x)(csc x - cot x) ) A) y = 0 B) y = y = - csc
More informationPhysics 312, Winter 2007, Practice Final
Physics 312, Winter 2007, Practice Final Time: Two hours Answer one of Question 1 or Question 2 plus one of Question 3 or Question 4 plus one of Question 5 or Question 6. Each question carries equal weight.
More informationStabilization and Passivity-Based Control
DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive
More informationRobot Control Basics CS 685
Robot Control Basics CS 685 Control basics Use some concepts from control theory to understand and learn how to control robots Control Theory general field studies control and understanding of behavior
More informationPhys 7221 Homework # 8
Phys 71 Homework # 8 Gabriela González November 15, 6 Derivation 5-6: Torque free symmetric top In a torque free, symmetric top, with I x = I y = I, the angular velocity vector ω in body coordinates with
More informationClassical Mechanics III (8.09) Fall 2014 Assignment 3
Classical Mechanics III (8.09) Fall 2014 Assignment 3 Massachusetts Institute of Technology Physics Department Due September 29, 2014 September 22, 2014 6:00pm Announcements This week we continue our discussion
More informationChapter 12. Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx
Chapter 1 Lecture Notes Chapter 1 Oscillatory Motion Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx When the mass is released, the spring will pull
More informationGravitation: Tensor Calculus
An Introduction to General Relativity Center for Relativistic Astrophysics School of Physics Georgia Institute of Technology Notes based on textbook: Spacetime and Geometry by S.M. Carroll Spring 2013
More informationTutorial Exercises: Geometric Connections
Tutorial Exercises: Geometric Connections 1. Geodesics in the Isotropic Mercator Projection When the surface of the globe is projected onto a flat map some aspects of the map are inevitably distorted.
More informationMath 6455 Nov 1, Differential Geometry I Fall 2006, Georgia Tech
Math 6455 Nov 1, 26 1 Differential Geometry I Fall 26, Georgia Tech Lecture Notes 14 Connections Suppose that we have a vector field X on a Riemannian manifold M. How can we measure how much X is changing
More informationChapter 8 continued. Rotational Dynamics
Chapter 8 continued Rotational Dynamics 8.4 Rotational Work and Energy Work to accelerate a mass rotating it by angle φ F W = F(cosθ)x x = rφ = Frφ Fr = τ (torque) = τφ r φ s F to x θ = 0 DEFINITION OF
More informationI. HARTLE CHAPTER 8, PROBLEM 2 (8 POINTS) where here an overdot represents d/dλ, must satisfy the geodesic equation (see 3 on problem set 4)
Physics 445 Solution for homework 5 Fall 2004 Cornell University 41 points) Steve Drasco 1 NOTE From here on, unless otherwise indicated we will use the same conventions as in the last two solutions: four-vectors
More informationLecture II: Rigid-Body Physics
Rigid-Body Motion Previously: Point dimensionless objects moving through a trajectory. Today: Objects with dimensions, moving as one piece. 2 Rigid-Body Kinematics Objects as sets of points. Relative distances
More informationcarroll/notes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general
http://pancake.uchicago.edu/ carroll/notes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity. As with any major theory in physics, GR has been
More informationPhysics 106b/196b Problem Set 9 Due Jan 19, 2007
Physics 06b/96b Problem Set 9 Due Jan 9, 2007 Version 3: January 8, 2007 This problem set focuses on dynamics in rotating coordinate systems (Section 5.2), with some additional early material on dynamics
More informationKinematics. Chapter Multi-Body Systems
Chapter 2 Kinematics This chapter first introduces multi-body systems in conceptual terms. It then describes the concept of a Euclidean frame in the material world, following the concept of a Euclidean
More informationArc Length and Riemannian Metric Geometry
Arc Length and Riemannian Metric Geometry References: 1 W F Reynolds, Hyperbolic geometry on a hyperboloid, Amer Math Monthly 100 (1993) 442 455 2 Wikipedia page Metric tensor The most pertinent parts
More informationFrom now on angles will be drawn with their vertex at the. The angle s initial ray will be along the positive. Think of the angle s
Fry Texas A&M University!! Math 150!! Chapter 8!! Fall 2014! 1 Chapter 8A Angles and Circles From now on angles will be drawn with their vertex at the The angle s initial ray will be along the positive.
More informationIs there a magnification paradox in gravitational lensing?
Is there a magnification paradox in gravitational lensing? Olaf Wucknitz wucknitz@astro.uni-bonn.de Astrophysics seminar/colloquium, Potsdam, 26 November 2007 Is there a magnification paradox in gravitational
More informationProblem 1. Mathematics of rotations
Problem 1. Mathematics of rotations (a) Show by algebraic means (i.e. no pictures) that the relationship between ω and is: φ, ψ, θ Feel free to use computer algebra. ω X = φ sin θ sin ψ + θ cos ψ (1) ω
More informationMath Section 4.3 Unit Circle Trigonometry
Math 10 - Section 4. Unit Circle Trigonometry An angle is in standard position if its vertex is at the origin and its initial side is along the positive x axis. Positive angles are measured counterclockwise
More informationPAPER 311 BLACK HOLES
MATHEMATICAL TRIPOS Part III Friday, 8 June, 018 9:00 am to 1:00 pm PAPER 311 BLACK HOLES Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY
More informationGeneral Relativity and Differential
Lecture Series on... General Relativity and Differential Geometry CHAD A. MIDDLETON Department of Physics Rhodes College November 1, 2005 OUTLINE Distance in 3D Euclidean Space Distance in 4D Minkowski
More informationUNIVERSITY OF DUBLIN
UNIVERSITY OF DUBLIN TRINITY COLLEGE JS & SS Mathematics SS Theoretical Physics SS TSM Mathematics Faculty of Engineering, Mathematics and Science school of mathematics Trinity Term 2015 Module MA3429
More informationChapter 8 continued. Rotational Dynamics
Chapter 8 continued Rotational Dynamics 8.4 Rotational Work and Energy Work to accelerate a mass rotating it by angle φ F W = F(cosθ)x x = s = rφ = Frφ Fr = τ (torque) = τφ r φ s F to s θ = 0 DEFINITION
More informationBackground Trigonmetry (2A) Young Won Lim 5/5/15
Background Trigonmetry (A) Copyright (c) 014 015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1. or
More informationSome aspects of the Geodesic flow
Some aspects of the Geodesic flow Pablo C. Abstract This is a presentation for Yael s course on Symplectic Geometry. We discuss here the context in which the geodesic flow can be understood using techniques
More informationSurfaces of Arbitrary Constant Negative Gaussian Curvature and Related Sine-Gordon Equations
Mathematica Aeterna, Vol.1, 011, no. 01, 1-11 Surfaces of Arbitrary Constant Negative Gaussian Curvature and Related Sine-Gordon Equations Paul Bracken Department of Mathematics, University of Texas, Edinburg,
More informationMath 122 Test 3. April 17, 2018
SI: Math Test 3 April 7, 08 EF: 3 4 5 6 7 8 9 0 Total Name Directions:. No books, notes or April showers. You may use a calculator to do routine arithmetic computations. You may not use your calculator
More information4 Riemannian geometry
Classnotes for Introduction to Differential Geometry. Matthias Kawski. April 18, 2003 83 4 Riemannian geometry 4.1 Introduction A natural first step towards a general concept of curvature is to develop
More informationGeneral Relativity and Cosmology Mock exam
Physikalisches Institut Mock Exam Universität Bonn 29. June 2011 Theoretische Physik SS 2011 General Relativity and Cosmology Mock exam Priv. Doz. Dr. S. Förste Exercise 1: Overview Give short answers
More informationCHAIN RULE: DAY 2 WITH TRIG FUNCTIONS. Section 2.4A Calculus AP/Dual, Revised /30/2018 1:44 AM 2.4A: Chain Rule Day 2 1
CHAIN RULE: DAY WITH TRIG FUNCTIONS Section.4A Calculus AP/Dual, Revised 018 viet.dang@humbleisd.net 7/30/018 1:44 AM.4A: Chain Rule Day 1 THE CHAIN RULE A. d dx f g x = f g x g x B. If f(x) is a differentiable
More informationPhysics 4B Winter 2018 Test 2
Physics 4B Winter 018 Test Name: Feb 15, 018 Please show your work! Answer as many questions as you can, in any order. Calculators are allowed. Books, notes, and internet connected devices are not allowed.
More informationControllability Analysis of A Two Degree of Freedom Nonlinear Attitude Control System
Controllability Analysis of A Two Degree of Freedom Nonlinear Attitude Control System Jinglai Shen, Amit K Sanyal, and N Harris McClamroch Department of Aerospace Engineering The University of Michigan
More informationDepartment of Architecture & Civil Engineering
MODE ANSWER age: 1 4. The students are given approximately 4 hours of lectures devoted to this topic. Thus the emphasis in the answer must be in demonstrating an understanding of the physical principals
More information2016 FAMAT Convention Mu Integration 1 = 80 0 = 80. dx 1 + x 2 = arctan x] k2
6 FAMAT Convention Mu Integration. A. 3 3 7 6 6 3 ] 3 6 6 3. B. For quadratic functions, Simpson s Rule is eact. Thus, 3. D.. B. lim 5 3 + ) 3 + ] 5 8 8 cot θ) dθ csc θ ) dθ cot θ θ + C n k n + k n lim
More informationTTK4150 Nonlinear Control Systems Solution 6 Part 2
TTK4150 Nonlinear Control Systems Solution 6 Part 2 Department of Engineering Cybernetics Norwegian University of Science and Technology Fall 2003 Solution 1 Thesystemisgivenby φ = R (φ) ω and J 1 ω 1
More informationTrajectory-tracking control of a planar 3-RRR parallel manipulator
Trajectory-tracking control of a planar 3-RRR parallel manipulator Chaman Nasa and Sandipan Bandyopadhyay Department of Engineering Design Indian Institute of Technology Madras Chennai, India Abstract
More informationCircular motion. Aug. 22, 2017
Circular motion Aug. 22, 2017 Until now, we have been observers to Newtonian physics through inertial reference frames. From our discussion of Newton s laws, these are frames which obey Newton s first
More informationAppendix B. Solutions to Chapter 2 Problems
Appendix B Solutions to Chapter Problems Problem Problem 4 Problem 3 5 Problem 4 6 Problem 5 6 Problem 6 7 Problem 7 8 Problem 8 8 Problem 9 Problem Problem 3 Problem et f ( ) be a function of one variable
More informationThe Cross Product. MATH 311, Calculus III. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan The Cross Product
The Cross Product MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Introduction Recall: the dot product of two vectors is a scalar. There is another binary operation on vectors
More informationEscape velocity and Schwarzschild s Solution for Black Holes
Escape velocity and Schwarzschild s Solution for Black Holes Amir Ali Tavajoh 1 1 Amir_ali3640@yahoo.com Introduction Scape velocity for every star or planet is different. As a star starts to collapse
More informationPHY6426/Fall 07: CLASSICAL MECHANICS HOMEWORK ASSIGNMENT #1 due by 9:35 a.m. Wed 09/05 Instructor: D. L. Maslov Rm.
PHY646/Fall 07: CLASSICAL MECHANICS HOMEWORK ASSIGNMENT # due by 9:35 a.m. Wed 09/05 Instructor: D. L. Maslov maslov@phys.ufl.edu 39-053 Rm. 4 Please help your instructor by doing your work neatly.. Goldstein,
More informationWilliam P. Thurston. The Geometry and Topology of Three-Manifolds
William P. Thurston The Geometry and Topology of Three-Manifolds Electronic version 1.1 - March 00 http://www.msri.org/publications/books/gt3m/ This is an electronic edition of the 1980 notes distributed
More informationGeometric Modelling Summer 2016
Geometric Modelling Summer 2016 Exercises Benjamin Karer M.Sc. http://gfx.uni-kl.de/~gm Benjamin Karer M.Sc. Geometric Modelling Summer 2016 1 Dierential Geometry Benjamin Karer M.Sc. Geometric Modelling
More informationGauß Curvature in Terms of the First Fundamental Form
Math 4530 Supplement March 6, 004 Gauß Curvature in Terms of the First Fundamental Form Andrejs Treibergs Abstract In these notes, we develop acceleration formulae for a general frame for a surface in
More informationA little taste of symplectic geometry
A little taste of symplectic geometry Timothy Goldberg Thursday, October 4, 2007 Olivetti Club Talk Cornell University 1 2 What is symplectic geometry? Symplectic geometry is the study of the geometry
More informationUniversity Calculus I. Worksheet # 8 Mar b. sin tan e. sin 2 sin 1 5. b. tan. c. sec sin 1 ( x )) cos 1 ( x )) f. csc. c.
MATH 6 WINTER 06 University Calculus I Worksheet # 8 Mar. 06-0 The topic covered by this worksheet is: Derivative of Inverse Functions and the Inverse Trigonometric functions. SamplesolutionstoallproblemswillbeavailableonDL,
More informationMath 122 Test 3. April 15, 2014
SI: Math 1 Test 3 April 15, 014 EF: 1 3 4 5 6 7 8 Total Name Directions: 1. No books, notes or 6 year olds with ear infections. You may use a calculator to do routine arithmetic computations. You may not
More information1 First and second variational formulas for area
1 First and second variational formulas for area In this chapter, we will derive the first and second variational formulas for the area of a submanifold. This will be useful in our later discussion on
More informationClassical Oscilators in General Relativity
Classical Oscilators in General Relativity arxiv:gr-qc/9709020v2 22 Oct 2000 Ion I. Cotăescu and Dumitru N. Vulcanov The West University of Timişoara, V. Pârvan Ave. 4, RO-1900 Timişoara, Romania Abstract
More information1. Geometry of the unit tangent bundle
1 1. Geometry of the unit tangent bundle The main reference for this section is [8]. In the following, we consider (M, g) an n-dimensional smooth manifold endowed with a Riemannian metric g. 1.1. Notations
More informationhas a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity.
http://preposterousuniverse.com/grnotes/ has a lot of good notes on GR and links to other pages. General Relativity Philosophy of general relativity. As with any major theory in physics, GR has been framed
More informationTorque and Simple Harmonic Motion
Torque and Simple Harmonic Motion Recall: Fixed Axis Rotation Angle variable Angular velocity Angular acceleration Mass element Radius of orbit Kinematics!! " d# / dt! " d 2 # / dt 2!m i Moment of inertia
More informationINSTRUCTOR SAMPLE E. Check that your exam contains 25 questions numbered sequentially. Answer Questions 1-25 on your scantron.
MATH 41 FINAL EXAM NAME SECTION NUMBER INSTRUCTOR SAMPLE E On your scantron, write and bubble your PSU ID, Section Number, and Test Version. Failure to correctly code these items may result in a loss of
More informationIntroduction to Elementary Particle Physics I
Physics 56400 Introduction to Elementary Particle Physics I Lecture 2 Fall 2018 Semester Prof. Matthew Jones Cross Sections Reaction rate: R = L σ The cross section is proportional to the probability of
More informationDIFFERENTIATION RULES
3 DIFFERENTIATION RULES DIFFERENTIATION RULES Before starting this section, you might need to review the trigonometric functions. DIFFERENTIATION RULES In particular, it is important to remember that,
More information1.7 Sums of series Example 1.7.1: 2. Real functions of one variable Example 1.7.2: 2.1 General definitions Example 2.1.3: Example 2.1.
7 Sums of series We often want to sum a series of terms, for example when we look at polynomials As we already saw, we abbreviate a sum of the form For example and u + u + + u r by r u i i= a n x n + a
More information7.6 Additional Problems
7in x 10in Felder c07_online.tex V3 - January 21, 2015 10:11 A.M. Page 26 26 Chapter 7 Linear Algebra II (nline) 7.6 Additional Problems 7.137 The picture below shows a rubber sheet with three small stones
More informationTwo simple ideas from calculus applied to Riemannian geometry
Calibrated Geometries and Special Holonomy p. 1/29 Two simple ideas from calculus applied to Riemannian geometry Spiro Karigiannis karigiannis@math.uwaterloo.ca Department of Pure Mathematics, University
More informationTHEODORE VORONOV DIFFERENTIAL GEOMETRY. Spring 2009
[under construction] 8 Parallel transport 8.1 Equation of parallel transport Consider a vector bundle E B. We would like to compare vectors belonging to fibers over different points. Recall that this was
More informationMatFys. Week 5, Dec , 2005
MatFys Week 5, Dec. 1-18, 005 Lectures During this week s lectures we shall finish our treatment of classical mechanics and begin the quantum theory part of the course with an introduction to Hilbert spaces.
More informationDYNAMICAL SYSTEMS
0.42 DYNAMICAL SYSTEMS Week Lecture Notes. What is a dynamical system? Probably the best way to begin this discussion is with arguably a most general and yet least helpful statement: Definition. A dynamical
More informationParameterizing the Pythagorean Triples
Parameterizing the Pythagorean Triples Saratoga Math Club January 8, 204 Introduction To ask for one or two Pythagorean triples is not a hard nor deep question at all. Common examples, like (3, 4, 5),
More informationImplicit Differentiation and Inverse Trigonometric Functions
Implicit Differentiation an Inverse Trigonometric Functions MATH 161 Calculus I J. Robert Buchanan Department of Mathematics Summer 2018 Explicit vs. Implicit Functions 0.5 1 y 0.0 y 2 0.5 3 4 1.0 0.5
More informationPhysics 480/581. Homework No. 10 Solutions: due Friday, 19 October, 2018
Physics 480/58 Homework No. 0 Solutions: due Friday, 9 October, 208. Using the coordinate bases for -forms, and their reciprocal bases for tangent vectors, and the usual form of the Schwarzschild metric,
More information