Small oscillations and normal modes
|
|
- Earl Perkins
- 5 years ago
- Views:
Transcription
1 Chapter 4 Small oscillations and normal modes 4.1 Linear oscillations Discuss a generalization of the harmonic oscillator problem: oscillations of a system of several degrees of freedom near the position of equilibrium remember for s = 1 L = 1 Mq q Vq, T > } {{} T q 0 minimum of the potential energy, x = q q 0 displacement expand Vq and Mq Vq = Vq 0 + dv }{{} dq q q d V q=q0 dq q q q=q0 const }{{} =0 }{{} k>0 Vq = const+ 1 kx +Ox 3 if k = 0 non-linear oscillation, higher derivatives important Mq = Mq 0 +Ox = m+ox restrict to orders Ox,Oẋ in general also Oxẋ in L for linear oscillations L = 1 mẋ 1 kx procedure is called: linearization of the original Lagrangian 4.1 results: linear differential equation with constant coefficients mẍ+kx = 0 6
2 ansatz for the solution x = Acosωt+ϕ algebraic equation mω +k = 0 ω = k m m > 0,k > 0: ω > 0 A amplitude, ω frequency, ωt+ϕ phase, ϕ initial phase now consider s degrees of freedom L = T Vq 1,...q s, T = 1 M ij q q i q j since q i q j is symmetric relative to i j, the coefficients M ij can be chosen symmetric M ij q = M ji q denote by q i0,i = 1,...,s the point of the minimum for the potential energy V expand V relative to x i = q i q i0 Vq = const+ 1 k ij x i x j +Ox 3 k, k ij = V q i q j = k ji qi0 i,j for x i 0 the potential energy increases with respect to its minimum at x i = 0 quadratic form k ij x i x j 0 expand M ij q i,j i,j M ij q = m ij +Ox k, m ij = M ij q k0 = m ji since T > 0 we have another quadratic form m ij ẋ i ẋ j 0 again restrict ourselves to Ox i x j,oẋ i ẋ j i,j L = 1 m ij ẋ i ẋ j k ij x i x j i,j 63
3 Introduce a matrix notation column vector x = x 1. x s and the transposed vector T denotes transposition as a raw vector mass matrix ˆm and matrix of elasticity ˆk m m 1s ˆm =... m s1... m ss introduce the scalar product x T = x 1,..., x s ˆk = k k 1s... k s1... k ss x,y x T y = s x i y i i=1 form of the Lagrangian with ẋ, ˆmẋ 0 and x, ˆkx 0 ˆk and ˆm are symmetric matrices: equation of motion L = 1 ẋ, ˆmẋ 1 x, ˆkx ˆm T = ˆm, ˆkT = ˆk ansatz for the solution d L dt ẋ = L x ˆmẍ+ˆkx = 0 x = Acosωt+ϕ non-trivial solutions for ω ˆm+ˆkA = 0 eigenvalue equation det ω ˆm+ˆk = 0 characteristic equation 64
4 assume ω 1,...,ω s eigenvalues solutions of the characteristic equation put the eigenvalues ω α into the eigenvalue equation, find the eigenvectors A α if A α is a solution, then aa α is also a solution introduce the normal mode vector x α t ω α ˆm+ˆkA α = 0 4. x α t = A α Q α t, Q α t = a α cosω α t+ϕ α Q α normal coordinate, ω α normal frequency a α arbitrary amplitude, ϕ α arbitrary phase the complete solution in components xt = s x α t = α=1 x i t = s α=1 s A α Q α t α=1 A α i Q α t the solution contains s arbitrary amplitudes a α and phases ϕ α which have to be found from the initial conditions x0 and ẋ0 In the language of vector algebra the transformation from vector x to Q is a linear transformation x = Û Q, x i = α U iα Q α and U iα A α i with x = x 1. x s, Q = Q 1. Q s, Û = A A s 1. A 1 s..... A s s for that transformation both quadratic forms for T and V become diagonal: s L = L α, L α = 1 M Q α α 1 K αq α 4.3 α=1 M α = A α, ˆmA α, K α = A α, ˆkA α 65
5 the Lagrangian has the form of s non-interacting oscillations we get the one-dimensional equations of motion for s independent normal coordinates Q α Note that normal frequencies are real: M α Qα +K α Q α = 0, α = 1,...,s ω α = K α M α = Aα, ˆmA α A α, ˆkA α 0 So, each normal coordinate Q α corresponds to an oscillation with one frequency ω α : Q α +ω αq α = 0 the oscillation is called normal mode Thesolutionintheoriginalvariablesx i isalinearsuperpositionofoscillationswithdifferent frequencies 66
6 4. Orthogonality of normal modes Assume two eigenvectors A α,a β are given with ω α ω β the corresponding normal mode vectors are x α = A α Q α and x β = A β Q β there is no reason that the scalar products of two eigenvectors or normal mode vectors have to vanish: A β, A α 0, x β, x α 0 Consider the eigenvector equations [see 4.] ω α,β ˆmA α,β = ˆkA α,β multiply from the left with A β,α ωα ωβ A β, ˆmA α = A α, ˆmA β = A β, ˆkA α A α, ˆkA β 4.4 and subtract A β, ˆmA α ωβ ω α A α, ˆmA β = A β, ˆkA α A α, ˆkA β it is known that for each real matrix ˆn: A,ˆnB = ˆn T A,B in our case ˆk and ˆm are real and symmetric ˆk T = ˆk, ˆm T = ˆm using A,B = B,A the difference vanishes ωα A β, ˆmA α ωβ ˆmA α, A β = A β, ˆkA α ˆkA α, A β 0 we get ωα ωβ A β, ˆmA α = 0 Therefore: A β, ˆmA α = 0 and [see 4.4] A β, ˆkA α = 0 67
7 Conclusions: The normal mode vectors x α and x β with ω α ω β are orthogonal to each other, if their scalar product is defined using the so called metric tensors ˆm or ˆk x α and x β are orthogonal in the metric of mass or elasticity x β, ˆmx α = 0, x β, ˆkx α = 0 In case of degenerate frequencies: e.g. the normal mode vectors x 1 and x with ω 1 = ω the linear combination ax 1 +bx is also a solution with the same frequency the space of solution is the plane given by the vectors x 1 and x choose a pair of independent vectors which satisfy the orthogonality condition in the mass metric togetherwithallothernormalmodevectorswithnon-degeneratefrequenciestheybuildthe basisforthenormalcoordinatesleadingtothediagonallagrangianofdecoupledoscillators
8 4.3 Thedouble pendulumin the fieldof constant gravity an example s =, use as generalized coordinates θ 1 and θ for the notations, see the figure Cartesian coordinates and the generalized coordinates are related via x 1 = l 1 sinθ 1, x = l 1 sinθ 1 +l sinθ z 1 = l 1 cosθ 1, z = l 1 cosθ 1 +l cosθ the Lagrangian can be found as follows see Landau/Lifshitz, 5, problem 1 L = 1 m 1 +m l1 θ m l θ +m l 1 l θ1 θ cosθ 1 θ +m 1 +m gl 1 cosθ 1 +m gl cosθ Consider the special case m 1 = m = m, l 1 = l, l = l L = 4ml θ ml θ +ml θ1 θ cosθ 1 θ +4mgl cosθ 1 +mgl cosθ linearize the Lagrangian θ i 1 L = 1 ml 8 θ 1 +4 θ 1 θ + θ 1 mgl4θ 1 +θ identify the matrices of mass and elasticity ˆm = ml, ˆk = mgl the equations of motion 8 θ 1 + θ +4ω 0θ 1 = 0, θ 1 + θ +ω 0θ = 0, ω 0 = 69 g l
9 Search the solution in the form θ1 x = = A cosωt+ϕ = θ A1 A cosωt+ϕ the eigenvalue equation ω ˆm+ˆkA = 0 [ 8 ml ω +mgl ] A1 A = 0 divide the matrix equation by ml and introduce ω0 = g l [ ] ω +ω0 A A = 0 8ω +4ω 0 ω ω ω +ω 0 A1 A = 0 solve the characteristic equation 8ω det +4ω0 ω ω ω +ω0 = 0 ω 1, = 3 5 ω 0 determine the eigenvectors A 1 and A from 8ω i +4ω 0 ω i ω i ω i +ω 0 A i 1 A i = 0 the A i are defined up to a normalization for A 1 we get the two equations 8ω 1 +4ω 0A 1 1 ω 1A 1 = 0, ω 1A ω 1 +ω 0A 1 = 0 from the second equation we obtain same result from the first equation choose A 1 1 = 1 A 1 = the normal mode vectors x 1, = A 1, Q 1, = A 1 = ω 1 A 1 ω0 ω1 1 = 5 1A , analogously A = ±1 ±1 a 5 1 1, cosω 1, t+ϕ 1, = Q 5 1 1, 70
10 complete solution x = x 1 +x Q = 1 Q 5 1Q Q with the angles θ 1 t = a 1 cosω 1 t+ϕ 1 a cosω t+ϕ θ1 θ t = 5 1a 1 cosω 1 t+ϕ a cosω t+ϕ the constants a 1,a,ϕ 1,ϕ are found from the initial conditions Let us also check that the Lagrangian becomes diagonal using the normal coordinates Q α 1 L = M Q α α 1 K αq α and ωα = K α M α we find identify α=1 mgl4θ1 +θ = [ 5 5 1Q 1 + ] 5+1Q ml 8 θ 1 +4 θ 1 θ + θ = [ Q Q ] θ mgl ml M 1, = 5 5±1ml, K 1, = 5 5 1mgl 5 1 ω1, g = 5±1 l = 5 1 g 4 l = 3 5 ω0 M 1, and K 1, can be found also from M α = A α, ˆmA α, K α = A α, ˆkA α In a plane with orthogonal axes θ 1,θ the eigenvectors A 1 and A give the directions of the new axes of normal coordinates Q 1 and Q the axes Q 1 and Q are not orthogonal to each other A 1,A = 1, = but the vectors A 1 and A are orthogonal in the metric of mass or elasticity A 1, ˆk mgl A = 1, = 1, =
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationMath 1302, Week 8: Oscillations
Math 302, Week 8: Oscillations T y eq Y y = y eq + Y mg Figure : Simple harmonic motion. At equilibrium the string is of total length y eq. During the motion we let Y be the extension beyond equilibrium,
More information! 4 4! o! +! h 4 o=0! ±= ± p i And back-substituting into the linear equations gave us the ratios of the amplitudes of oscillation:.»» = A p e i! +t»»
Topic 6: Coupled Oscillators and Normal Modes Reading assignment: Hand and Finch Chapter 9 We are going to be considering the general case of a system with N degrees of freedome close to one of its stable
More informationAssignment 6. Using the result for the Lagrangian for a double pendulum in Problem 1.22, we get
Assignment 6 Goldstein 6.4 Obtain the normal modes of vibration for the double pendulum shown in Figure.4, assuming equal lengths, but not equal masses. Show that when the lower mass is small compared
More informationRotational motion of a rigid body spinning around a rotational axis ˆn;
Physics 106a, Caltech 15 November, 2018 Lecture 14: Rotations The motion of solid bodies So far, we have been studying the motion of point particles, which are essentially just translational. Bodies with
More informationMatrix Representation
Matrix Representation Matrix Rep. Same basics as introduced already. Convenient method of working with vectors. Superposition Complete set of vectors can be used to express any other vector. Complete set
More informationIntroduction to Matrix Algebra
Introduction to Matrix Algebra August 18, 2010 1 Vectors 1.1 Notations A p-dimensional vector is p numbers put together. Written as x 1 x =. x p. When p = 1, this represents a point in the line. When p
More informationCoupled Oscillators Monday, 29 October 2012 normal mode Figure 1:
Coupled Oscillators Monday, 29 October 2012 In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid. Physics 111 θ 1 l m k θ 2
More informationIn which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid.
Coupled Oscillators Wednesday, 26 October 2011 In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid. Physics 111 θ 1 l k θ 2
More informationMath Bootcamp An p-dimensional vector is p numbers put together. Written as. x 1 x =. x p
Math Bootcamp 2012 1 Review of matrix algebra 1.1 Vectors and rules of operations An p-dimensional vector is p numbers put together. Written as x 1 x =. x p. When p = 1, this represents a point in the
More informationLinear Algebra & Geometry why is linear algebra useful in computer vision?
Linear Algebra & Geometry why is linear algebra useful in computer vision? References: -Any book on linear algebra! -[HZ] chapters 2, 4 Some of the slides in this lecture are courtesy to Prof. Octavia
More informationAppendix A Equations of Motion in the Configuration and State Spaces
Appendix A Equations of Motion in the Configuration and State Spaces A.1 Discrete Linear Systems A.1.1 Configuration Space Consider a system with a single degree of freedom and assume that the equation
More information10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections )
c Dr. Igor Zelenko, Fall 2017 1 10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections 7.2-7.4) 1. When each of the functions F 1, F 2,..., F n in right-hand side
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More information2. Matrix Algebra and Random Vectors
2. Matrix Algebra and Random Vectors 2.1 Introduction Multivariate data can be conveniently display as array of numbers. In general, a rectangular array of numbers with, for instance, n rows and p columns
More informationSeminar 6: COUPLED HARMONIC OSCILLATORS
Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached
More information1.1 To observe the normal modes of oscillation of a two degree of freedom system.
I. BJECTIVES. To observe the normal modes of oscillation of a two degree of freedom system.. To determine the natural frequencies and mode shapes of the system from solution of the Eigenvalue problem..3
More informationMulti Degrees of Freedom Systems
Multi Degrees of Freedom Systems MDOF s http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 9, 07 Outline, a System
More informationReview of Linear Algebra
Review of Linear Algebra Definitions An m n (read "m by n") matrix, is a rectangular array of entries, where m is the number of rows and n the number of columns. 2 Definitions (Con t) A is square if m=
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More information22.3. Repeated Eigenvalues and Symmetric Matrices. Introduction. Prerequisites. Learning Outcomes
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationIf the symmetry axes of a uniform symmetric body coincide with the coordinate axes, the products of inertia (Ixy etc.
Prof. O. B. Wright, Autumn 007 Mechanics Lecture 9 More on rigid bodies, coupled vibrations Principal axes of the inertia tensor If the symmetry axes of a uniform symmetric body coincide with the coordinate
More informationSolutions to homework assignment #3 Math 119B UC Davis, Spring = 12x. The Euler-Lagrange equation is. 2q (x) = 12x. q(x) = x 3 + x.
1. Find the stationary points of the following functionals: (a) 1 (q (x) + 1xq(x)) dx, q() =, q(1) = Solution. The functional is 1 L(q, q, x) where L(q, q, x) = q + 1xq. We have q = q, q = 1x. The Euler-Lagrange
More informationMatrices and Linear Algebra
Contents Quantitative methods for Economics and Business University of Ferrara Academic year 2017-2018 Contents 1 Basics 2 3 4 5 Contents 1 Basics 2 3 4 5 Contents 1 Basics 2 3 4 5 Contents 1 Basics 2
More informationIn the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as
2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,
More information5 Linear Algebra and Inverse Problem
5 Linear Algebra and Inverse Problem 5.1 Introduction Direct problem ( Forward problem) is to find field quantities satisfying Governing equations, Boundary conditions, Initial conditions. The direct problem
More informationOscillations and Waves
Oscillations and Waves Somnath Bharadwaj and S. Pratik Khastgir Department of Physics and Meteorology IIT Kharagpur Module : Oscillations Lecture : Oscillations Oscillations are ubiquitous. It would be
More informationAA 242B / ME 242B: Mechanical Vibrations (Spring 2016)
AA 242B / ME 242B: Mechanical Vibrations (Spring 206) Solution of Homework #3 Control Tab Figure : Schematic for the control tab. Inadequacy of a static-test A static-test for measuring θ would ideally
More information11 a 12 a 13 a 21 a 22 a b 12 b 13 b 21 b 22 b b 11 a 12 + b 12 a 13 + b 13 a 21 + b 21 a 22 + b 22 a 23 + b 23
Chapter 2 (3 3) Matrices The methods used described in the previous chapter for solving sets of linear equations are equally applicable to 3 3 matrices. The algebra becomes more drawn out for larger matrices,
More information1 Review of simple harmonic oscillator
MATHEMATICS 7302 (Analytical Dynamics YEAR 2017 2018, TERM 2 HANDOUT #8: COUPLED OSCILLATIONS AND NORMAL MODES 1 Review of simple harmonic oscillator In MATH 1301/1302 you studied the simple harmonic oscillator:
More informationA = 3 B = A 1 1 matrix is the same as a number or scalar, 3 = [3].
Appendix : A Very Brief Linear ALgebra Review Introduction Linear Algebra, also known as matrix theory, is an important element of all branches of mathematics Very often in this course we study the shapes
More informationVectors To begin, let us describe an element of the state space as a point with numerical coordinates, that is x 1. x 2. x =
Linear Algebra Review Vectors To begin, let us describe an element of the state space as a point with numerical coordinates, that is x 1 x x = 2. x n Vectors of up to three dimensions are easy to diagram.
More informationA VERY BRIEF LINEAR ALGEBRA REVIEW for MAP 5485 Introduction to Mathematical Biophysics Fall 2010
A VERY BRIEF LINEAR ALGEBRA REVIEW for MAP 5485 Introduction to Mathematical Biophysics Fall 00 Introduction Linear Algebra, also known as matrix theory, is an important element of all branches of mathematics
More informationLecture 7. Econ August 18
Lecture 7 Econ 2001 2015 August 18 Lecture 7 Outline First, the theorem of the maximum, an amazing result about continuity in optimization problems. Then, we start linear algebra, mostly looking at familiar
More informationRepeated Eigenvalues and Symmetric Matrices
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationENGI 9420 Lecture Notes 2 - Matrix Algebra Page Matrix operations can render the solution of a linear system much more efficient.
ENGI 940 Lecture Notes - Matrix Algebra Page.0. Matrix Algebra A linear system of m equations in n unknowns, a x + a x + + a x b (where the a ij and i n n a x + a x + + a x b n n a x + a x + + a x b m
More informationLinear Algebra Review. Vectors
Linear Algebra Review 9/4/7 Linear Algebra Review By Tim K. Marks UCSD Borrows heavily from: Jana Kosecka http://cs.gmu.edu/~kosecka/cs682.html Virginia de Sa (UCSD) Cogsci 8F Linear Algebra review Vectors
More informationQuantum Computing Lecture 2. Review of Linear Algebra
Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces
More informationMathematical Methods wk 2: Linear Operators
John Magorrian, magog@thphysoxacuk These are work-in-progress notes for the second-year course on mathematical methods The most up-to-date version is available from http://www-thphysphysicsoxacuk/people/johnmagorrian/mm
More informationM2A2 Problem Sheet 3 - Hamiltonian Mechanics
MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,
More informationLinear Algebra, part 2 Eigenvalues, eigenvectors and least squares solutions
Linear Algebra, part 2 Eigenvalues, eigenvectors and least squares solutions Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Main problem of linear algebra 2: Given
More information1 Matrices and matrix algebra
1 Matrices and matrix algebra 1.1 Examples of matrices A matrix is a rectangular array of numbers and/or variables. For instance 4 2 0 3 1 A = 5 1.2 0.7 x 3 π 3 4 6 27 is a matrix with 3 rows and 5 columns
More informationWeek Quadratic forms. Principal axes theorem. Text reference: this material corresponds to parts of sections 5.5, 8.2,
Math 051 W008 Margo Kondratieva Week 10-11 Quadratic forms Principal axes theorem Text reference: this material corresponds to parts of sections 55, 8, 83 89 Section 41 Motivation and introduction Consider
More informationBasic Concepts in Matrix Algebra
Basic Concepts in Matrix Algebra An column array of p elements is called a vector of dimension p and is written as x p 1 = x 1 x 2. x p. The transpose of the column vector x p 1 is row vector x = [x 1
More informationIntroductory Physics. Week 2015/05/29
2015/05/29 Part I Summary of week 6 Summary of week 6 We studied the motion of a projectile under uniform gravity, and constrained rectilinear motion, introducing the concept of constraint force. Then
More informationKnowledge Discovery and Data Mining 1 (VO) ( )
Knowledge Discovery and Data Mining 1 (VO) (707.003) Review of Linear Algebra Denis Helic KTI, TU Graz Oct 9, 2014 Denis Helic (KTI, TU Graz) KDDM1 Oct 9, 2014 1 / 74 Big picture: KDDM Probability Theory
More informationWaves & Normal Modes. Matt Jarvis
Waves & Normal Modes Matt Jarvis February 27, 2016 Contents 1 Oscillations 2 1.0.1 Simple Harmonic Motion - revision................... 2 2 Normal Modes 5 2.1 The coupled pendulum..............................
More informationLecture 10: Eigenvectors and eigenvalues (Numerical Recipes, Chapter 11)
Lecture 1: Eigenvectors and eigenvalues (Numerical Recipes, Chapter 11) The eigenvalue problem, Ax= λ x, occurs in many, many contexts: classical mechanics, quantum mechanics, optics 22 Eigenvectors and
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationConsider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.
Chapter 4 Energy and Stability 4.1 Energy in 1D Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be T = 1 2 mẋ2 and the potential energy
More informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More informationMath Matrix Algebra
Math 44 - Matrix Algebra Review notes - (Alberto Bressan, Spring 7) sec: Orthogonal diagonalization of symmetric matrices When we seek to diagonalize a general n n matrix A, two difficulties may arise:
More informationCHAPTER 12 OSCILLATORY MOTION
CHAPTER 1 OSCILLATORY MOTION Before starting the discussion of the chapter s concepts it is worth to define some terms we will use frequently in this chapter: 1. The period of the motion, T, is the time
More informationHamilton-Jacobi theory
Hamilton-Jacobi theory November 9, 04 We conclude with the crowning theorem of Hamiltonian dynamics: a proof that for any Hamiltonian dynamical system there exists a canonical transformation to a set of
More informationLinear Algebra & Geometry why is linear algebra useful in computer vision?
Linear Algebra & Geometry why is linear algebra useful in computer vision? References: -Any book on linear algebra! -[HZ] chapters 2, 4 Some of the slides in this lecture are courtesy to Prof. Octavia
More informationLinear Algebra 1 Exam 1 Solutions 6/12/3
Linear Algebra 1 Exam 1 Solutions 6/12/3 Question 1 Consider the linear system in the variables (x, y, z, t, u), given by the following matrix, in echelon form: 1 2 1 3 1 2 0 1 1 3 1 4 0 0 0 1 2 3 Reduce
More informationLinear Algebra (Review) Volker Tresp 2017
Linear Algebra (Review) Volker Tresp 2017 1 Vectors k is a scalar (a number) c is a column vector. Thus in two dimensions, c = ( c1 c 2 ) (Advanced: More precisely, a vector is defined in a vector space.
More informationChap. 15: Simple Harmonic Motion
Chap. 15: Simple Harmonic Motion Announcements: CAPA is due next Tuesday and next Friday. Web page: http://www.colorado.edu/physics/phys1110/phys1110_sp12/ Examples of periodic motion vibrating guitar
More informationSymmetric and anti symmetric matrices
Symmetric and anti symmetric matrices In linear algebra, a symmetric matrix is a square matrix that is equal to its transpose. Formally, matrix A is symmetric if. A = A Because equal matrices have equal
More informationMATRIX ALGEBRA. or x = (x 1,..., x n ) R n. y 1 y 2. x 2. x m. y m. y = cos θ 1 = x 1 L x. sin θ 1 = x 2. cos θ 2 = y 1 L y.
as Basics Vectors MATRIX ALGEBRA An array of n real numbers x, x,, x n is called a vector and it is written x = x x n or x = x,, x n R n prime operation=transposing a column to a row Basic vector operations
More informationL = 1 2 a(q) q2 V (q).
Physics 3550 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6, 11.1 11.3 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More informationAppendix A: Matrices
Appendix A: Matrices A matrix is a rectangular array of numbers Such arrays have rows and columns The numbers of rows and columns are referred to as the dimensions of a matrix A matrix with, say, 5 rows
More informationChapter 6: Vector Analysis
Chapter 6: Vector Analysis We use derivatives and various products of vectors in all areas of physics. For example, Newton s 2nd law is F = m d2 r. In electricity dt 2 and magnetism, we need surface and
More informationClassical Mechanics. Luis Anchordoqui
1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body
More informationECEEN 5448 Fall 2011 Homework #5 Solutions
ECEEN 5448 Fall 211 Homework #5 Solutions Professor David G. Meyer December 8, 211 1. Consider the 1-dimensional time-varying linear system ẋ t (u x) (a) Find the state-transition matrix, Φ(t, τ). Here
More informationOR MSc Maths Revision Course
OR MSc Maths Revision Course Tom Byrne School of Mathematics University of Edinburgh t.m.byrne@sms.ed.ac.uk 15 September 2017 General Information Today JCMB Lecture Theatre A, 09:30-12:30 Mathematics revision
More informationLinear Algebra (Review) Volker Tresp 2018
Linear Algebra (Review) Volker Tresp 2018 1 Vectors k, M, N are scalars A one-dimensional array c is a column vector. Thus in two dimensions, ( ) c1 c = c 2 c i is the i-th component of c c T = (c 1, c
More informationMath 1060 Linear Algebra Homework Exercises 1 1. Find the complete solutions (if any!) to each of the following systems of simultaneous equations:
Homework Exercises 1 1 Find the complete solutions (if any!) to each of the following systems of simultaneous equations: (i) x 4y + 3z = 2 3x 11y + 13z = 3 2x 9y + 2z = 7 x 2y + 6z = 2 (ii) x 4y + 3z =
More information. The following is a 3 3 orthogonal matrix: 2/3 1/3 2/3 2/3 2/3 1/3 1/3 2/3 2/3
Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. An n n matrix
More informationEXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. 1. Determinants
EXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. Determinants Ex... Let A = 0 4 4 2 0 and B = 0 3 0. (a) Compute 0 0 0 0 A. (b) Compute det(2a 2 B), det(4a + B), det(2(a 3 B 2 )). 0 t Ex..2. For
More information7.1 Creation and annihilation operators
Chapter 7 Second Quantization Creation and annihilation operators. Occupation number. Anticommutation relations. Normal product. Wick s theorem. One-body operator in second quantization. Hartree- Fock
More informationPhysics III: Final Solutions (Individual and Group)
Massachusetts Institute of Technology MITES 7 Physics III First and Last Name: Physics III: Final Solutions (Individual and Group Instructor: Mobolaji Williams Teaching Assistant: Rene García (Tuesday
More information21.55 Worksheet 7 - preparation problems - question 1:
Dynamics 76. Worksheet 7 - preparation problems - question : A coupled oscillator with two masses m and positions x (t) and x (t) is described by the following equations of motion: ẍ x + 8x ẍ x +x A. Write
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2016 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationMATHEMATICAL PHYSICS
MATHEMATICAL PHYSICS Third Year SEMESTER 1 015 016 Classical Mechanics MP350 Prof. S. J. Hands, Prof. D. M. Heffernan, Dr. J.-I. Skullerud and Dr. M. Fremling Time allowed: 1 1 hours Answer two questions
More informationPhys 201. Matrices and Determinants
Phys 201 Matrices and Determinants 1 1.1 Matrices 1.2 Operations of matrices 1.3 Types of matrices 1.4 Properties of matrices 1.5 Determinants 1.6 Inverse of a 3 3 matrix 2 1.1 Matrices A 2 3 7 =! " 1
More informationDIAGONALIZATION OF THE STRESS TENSOR
DIAGONALIZATION OF THE STRESS TENSOR INTRODUCTION By the use of Cauchy s theorem we are able to reduce the number of stress components in the stress tensor to only nine values. An additional simplification
More informationSection 1.6. M N = [a ij b ij ], (1.6.2)
The Calculus of Functions of Several Variables Section 16 Operations with Matrices In the previous section we saw the important connection between linear functions and matrices In this section we will
More informationSolution Set Five. 2 Problem #2: The Pendulum of Doom Equation of Motion Simple Pendulum Limit Visualization...
: Solution Set Five Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein November 4, 2015 Contents 1 Problem #1: Coupled Mass Oscillator System. 2 1.1 Normal Modes.......................................
More informationAddition of Angular Momenta
Addition of Angular Momenta What we have so far considered to be an exact solution for the many electron problem, should really be called exact non-relativistic solution. A relativistic treatment is needed
More informationj=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.
Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. Let u = [u
More informationCS 246 Review of Linear Algebra 01/17/19
1 Linear algebra In this section we will discuss vectors and matrices. We denote the (i, j)th entry of a matrix A as A ij, and the ith entry of a vector as v i. 1.1 Vectors and vector operations A vector
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationSolutions 2: Simple Harmonic Oscillator and General Oscillations
Massachusetts Institute of Technology MITES 2017 Physics III Solutions 2: Simple Harmonic Oscillator and General Oscillations Due Wednesday June 21, at 9AM under Rene García s door Preface: This problem
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
More informationq n. Q T Q = I. Projections Least Squares best fit solution to Ax = b. Gram-Schmidt process for getting an orthonormal basis from any basis.
Exam Review Material covered by the exam [ Orthogonal matrices Q = q 1... ] q n. Q T Q = I. Projections Least Squares best fit solution to Ax = b. Gram-Schmidt process for getting an orthonormal basis
More informationLinear Algebra Review (Course Notes for Math 308H - Spring 2016)
Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,
More informationENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations
ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential
More informationVectors and matrices: matrices (Version 2) This is a very brief summary of my lecture notes.
Vectors and matrices: matrices (Version 2) This is a very brief summary of my lecture notes Matrices and linear equations A matrix is an m-by-n array of numbers A = a 11 a 12 a 13 a 1n a 21 a 22 a 23 a
More information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
More informationIntroduction to Matrices
214 Analysis and Design of Feedback Control Systems Introduction to Matrices Derek Rowell October 2002 Modern system dynamics is based upon a matrix representation of the dynamic equations governing the
More informationBackground Mathematics (2/2) 1. David Barber
Background Mathematics (2/2) 1 David Barber University College London Modified by Samson Cheung (sccheung@ieee.org) 1 These slides accompany the book Bayesian Reasoning and Machine Learning. The book and
More informationHomework 2. Solutions T =
Homework. s Let {e x, e y, e z } be an orthonormal basis in E. Consider the following ordered triples: a) {e x, e x + e y, 5e z }, b) {e y, e x, 5e z }, c) {e y, e x, e z }, d) {e y, e x, 5e z }, e) {
More informationMathematical foundations - linear algebra
Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar
More informationNIU PHYS 500, Fall 2006 Classical Mechanics Solutions for HW6. Solutions
NIU PHYS 500, Fall 006 Classical Mechanics Solutions for HW6 Assignment: HW6 [40 points] Assigned: 006/11/10 Due: 006/11/17 Solutions P6.1 [4 + 3 + 3 = 10 points] Consider a particle of mass m moving in
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2006 Christopher J. Cramer. Lecture 5, January 27, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2006 Christopher J. Cramer Lecture 5, January 27, 2006 Solved Homework (Homework for grading is also due today) We are told
More informationDifferential Equations
Differential Equations A differential equation (DE) is an equation which involves an unknown function f (x) as well as some of its derivatives. To solve a differential equation means to find the unknown
More information