Let us start with a two dimensional case. We consider a vector ( x,

Similar documents
Two Coupled Oscillators / Normal Modes

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

IB Physics Kinematics Worksheet

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

Math 334 Fall 2011 Homework 11 Solutions

KINEMATICS IN ONE DIMENSION

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

Chapter 6. Systems of First Order Linear Differential Equations

15. Vector Valued Functions

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

Math 315: Linear Algebra Solutions to Assignment 6

10. State Space Methods

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Some Basic Information about M-S-D Systems

Let ( α, β be the eigenvector associated with the eigenvalue λ i

BEng (Hons) Telecommunications. Examinations for / Semester 2

Announcements: Warm-up Exercise:

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

15. Bicycle Wheel. Graph of height y (cm) above the axle against time t (s) over a 6-second interval. 15 bike wheel

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

Final Spring 2007

Roller-Coaster Coordinate System

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

SOLUTIONS TO ECE 3084

DIFFERENTIAL GEOMETRY HW 5

Chapter 12: Velocity, acceleration, and forces

Chapter 2. First Order Scalar Equations

Solutions from Chapter 9.1 and 9.2

INSTANTANEOUS VELOCITY

4.5 Constant Acceleration

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Linear Response Theory: The connection between QFT and experiments

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

Homework sheet Exercises done during the lecture of March 12, 2014

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

Lecture Notes 2. The Hilbert Space Approach to Time Series

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

18 Biological models with discrete time

HW6: MRI Imaging Pulse Sequences (7 Problems for 100 pts)

Chapter 7: Solving Trig Equations

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013

Introduction to Physical Oceanography Homework 5 - Solutions

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

ON THE DEGREES OF RATIONAL KNOTS

Chapter 3 Kinematics in Two Dimensions

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Trajectory planning in Cartesian space

System of Linear Differential Equations

Linear Dynamic Models

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

From Particles to Rigid Bodies

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

The equation to any straight line can be expressed in the form:

Matlab and Python programming: how to get started

BU Macro BU Macro Fall 2008, Lecture 4

MATH 2050 Assignment 9 Winter Do not need to hand in. 1. Find the determinant by reducing to triangular form for the following matrices.

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2002 November 14 Exam III Physics 191

Traveling Waves. Chapter Introduction

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

3, so θ = arccos

k 1 k 2 x (1) x 2 = k 1 x 1 = k 2 k 1 +k 2 x (2) x k series x (3) k 2 x 2 = k 1 k 2 = k 1+k 2 = 1 k k 2 k series

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

Morning Time: 1 hour 30 minutes Additional materials (enclosed):

Math 116 Practice for Exam 2

CHAPTER 12 DIRECT CURRENT CIRCUITS

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Lie Derivatives operator vector field flow push back Lie derivative of

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims

Exercises: Similarity Transformation

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?

The Arcsine Distribution

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

Problemas das Aulas Práticas

After the completion of this section the student. Theory of Linear Systems of ODEs. Autonomous Systems. Review Questions and Exercises

The average rate of change between two points on a function is d t

6.003 Homework #9 Solutions

Differential Geometry: Revisiting Curvatures

Non-uniform circular motion *

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

AP Calculus BC Chapter 10 Part 1 AP Exam Problems

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

1. VELOCITY AND ACCELERATION

Fishing limits and the Logistic Equation. 1

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Elements of Computer Graphics

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

1 Solutions to selected problems

The Paradox of Twins Described in a Three-dimensional Space-time Frame

Transcription:

Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our prime ineres is in hree dimensions, where proeins are embedded. We expec from roaion marices o have he following properies: is a square real marix such ha = I (which means ha he roaion marix does no change he lengh of he vecor. See below) and de ( ) = 1 (which means ha we do no allow mirror images). Le us sar wih a wo dimensional case. We consider a vecor ( x, y ) ha we wish o roae in he plane by an angle (in he direcion opposie o he clock). The norm of he vecor is kep fixed. The diagram below describes he operaion of ineres I is possible o hink abou he roaion of he vecor ( x, y ) as a roaion of wo vecors ( x,0) and oo!). Consider firs he vecor componen x ( x,0) and he effec of he roaion on i. Roaing by he vecor ( x,0) modifies he old ino a new vecor ( xcos, xsin ) Similarly for he y componen of he vecor ( 0, y ) we find ha couner clockwise roaion by an angle generaes he vecor ysin, ycos( ). The complee new 0, y which we will finally add up o a single operaion (i is more convenien vecor afer he roaion will be xcos ysin, xsin + ycos ( ). The above operaion which we did by using geomery and some inuiion can be pu in a more elegan form using a roaion marix cos = sin sin cos

y y cos sin x xcos sin sin cos = y xsin + cos Some properies of cos sin cos sin sin cos sin cos cos sin cos sin sin cos = = + 1 0 = = = I sin cos cos sin sin cos + 0 1 ( ) cos sin cos de = = cos + sin = 1 sin Roaion back by should bring us o he place in which we sared (i.e., = I =. Indeed we have cos sin cos sin = = = sin cos sin cos Le = φ1+ φ2, we inuiively expec ha he roaion by a sum can be described as sequence of roaions, i.e. ( ) = ( φ ) ( φ ) A proof: 1 2 cos( φ1) sin ( φ1) cos( φ2) sin ( φ2) ( φ1) ( φ2) = = sin ( φ1) cos( φ1) sin ( φ2) cos( φ2) cos( φ1) cos( φ2) sin ( φ1) sin ( φ2) cos( φ1) sin ( φ2) sin ( φ1) cos( φ2) sin ( φ1) cos( φ2) + cos( φ1) sin ( φ2) sin ( φ1) sin ( φ2) + cos( φ1) cos( φ2) cos( φ1+ φ2) sin ( φ1+ φ2) cos sin sin ( φ + φ ) cos( φ + φ ) sin cos = = = 1 2 1 2 Consider a vecor V. A roaion marix does no change he norm of he vecor. Hence =. This is rivial o show once we remember ha V V V V V V herefore have V V = V V = V V = V V =, we

The disance beween any wo poins is no changing by roaion. Consider wo poins in he plane: ( x1, y 1) and ( x2, y 2). The disance beween he wo poins, d 12, is defined as x1 x2 he norm of he vecor difference d12 = = ( x1 x2 ) + ( y1 y2 ) bu we jus y1 y2 showed ha he lengh of a vecor does no change upon roaion. The disance is he lengh of he vecor difference and canno change by roaion. We define a rigid body ransformaion as a change in he coordinaes ha does no affec he disance beween any wo poins of he rigid body (we also require ha no mirror imaging akes place). Noe ha no all he marices saisfying = I are proper roaions. For example for 1 0 = 0 1 we sill have 1 0 1 0 1 0 = = = I. However for he 0 1 0 1 0 1 second condiion menioned above ( de ( ) = 1) we have -- de ( ) = 1 ( 1) = 1 which is no a proper roaion. Anoher way of hinking abou he condiion on he deerminan above is o consider small roaions. As we have seen a roaion by an angle can be divided o a sequence of small roaions. We can make a large number of such divisions and approach infiniesimal roaions. We wrie = lim ( ) N N i = 1,.., N or 1 N ( ) = ( N) ( N) i= 1,..., N N 1 Such decomposiion should always be possible for a real roaion (one has o be careful hough when applying such limis. The marix of an infiniesimal roaion, as wrien, does no saisfy he condiion we se for a roaion). Wihou proof we commen ha marices wih deerminan 1 canno be decomposed o a produc of infiniesimal roaion. This can be used as ye anoher definiion of roaion marices (we mus be able o decompose hem o a produc of smaller roaions). 1 0 The marix we jus considered, reflecs he Y axis hrough he origin and 0 1 generaes a mirror image. I is impossible o find a sequence of small roaions ha akes a vecor coninuously from he wo reflecion saes of he Y axis. Roaions in hree dimensions Roaions in hree dimensions are a simple exension of he roaions in wo dimensions we jus discussed. The difference is ha in hree dimensions we mus define a roaion axis. The wo dimensional roaions we considered so far can be hough as a special case

of roaions in hree dimensions in which he axis of roaion is he Z axis. Consider he following marix: cos sin 0 = sin cos 0 0 0 1 when applied ono a vecor in hree dimensions v ( x, y, z) cos sin 0 x xcos ysin v = sin cos( ) 0 y = xsin ( ) + y cos = he resul is: 0 0 1z z which is exacly he same resul we have for a roaion in wo dimensions in he xy plane. Here Z is he roaion axis which (in general) remains unchanged when we applied he roaion on a vecor along is direcion. We can find he roaion axis of an arbirary roaion marix in hree dimension by solving he following linear problem e = 1 e where e is a vecor along he axis of roaion. Since he axis of roaion is no changing when muliplies i (oherwise i is NOT he axis of roaion), i mus be an eigenvecor of wih an eigenvalue of 1. Deermining e is a linear problem ha can be solved (for example) using he Gaussian eliminaion we discussed earlier. Noe ha if we muliply e by any scalar we sill obain a valid soluion. We herefore choose e o be real (If here is a soluion i can be made real since he marix and he eigenvalue are real) and normalized i o one: ee= 1. The axis of roaion is herefore deermined wih only wo parameers wo angles (he hird is deermined by he normalizaion). The wo parameers can be chosen as he relaive posiion of he axis of roaion wih respec o he Z axis -- cos( Θ ) and he angle of he axis of roaion wih respec o X axis. The las is obained by projecing he axis of roaion firs o he XY plane and hen measuring he relaive posiion wih respec o he X axis -- cos( Ψ ). The skech below demonsraes he process and he way we define a roaion in hree dimensions. As before is a roaion in a plane. The plane is now deermined as he se of vecors perpendicular o he axis of roaion. I.e. he se of vecors perpendicular o e Θ Z e Ψ Y X

A roaion in hree dimensions is herefore compleely deermined by hree parameers: ( ΘΨ,, ). Of course he choice of he 3 parameer is unique and here are oher ways of uniquely deermining a roaion in space, a common example are he Euler angles ha can be found in any exbook on classical mechanics (e.g. Goldsein). As we noed before he deerminan of a roaion marix is 1. Wha does i say on he range of eigenvalues ha a roaion marix may have? We spend considerable ime on roaion marices because we have a very concree quesion abou hese marices. Given wo proein srucures (modeled as rigid bodies) wha kind of roaion we can apply on one of hem so ha he wo srucures will overlap as bes as possible. Hence we are back o proeins and biology