! 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»»
|
|
- Clare McDonald
- 5 years ago
- Views:
Transcription
1 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 equilibrium points. We are going to show that for small displacements from equilibrium, the system acts like N independent SHOs, usually with N different characteristic frequencies. One or more of these SHOs cn be excited, depending on the initial conditions. If only one oscillation frequency is excited, the N degrees of freedom move synchronously at a common mode frequency. The ratios of the displacements for the different degrees of freedom, known as mode displacement ratios are an intrinsic characteristic of the normal mode that is oscillating. Any motion can be described as a linear combination of modes. We will now develop this formally.. The Double Pendulum Revisited Recall our solution for the double pendulum. The Lagrangian for this system is: L = ml h _ + _ + cos ( ) _ _ i + mgl ( cos + cos ) From this, by considering small displacements (throwing out high order terms and making small angle approximations) we got two coupled differential equations: + =! o + = g! o We looked explicitly for solutions where both masses oscillate with the same frequency: ; (t) = ; e ±i!t Note that we did not nd a general solution to the problem where the oscillations of ; are not necessarily at the same frequency. By substituting our solution into the EOM, we found two linear equations (! o! )! =0! +! o! =0 which can have non-trivial solutions only if the determinant of the coefcients vanishes. So, we got two roots corresponding to two different frequencies:»» (! o! )! =0! (! o! ) The determinant vanishing gives us:! o!! 4 =0
2 ! 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»» = A p e i! t We are now going to develop this more formally for arbitrary systems with N generalized co-ordinates.. A Review of Matrix Algebra The formalism we are about to develop requires some matrix algebra, so lets do a quick review. A collection of basic facts: Matrix multiplication is not commutative: if we consider two matrices A and B:! o C = AB 6= BA We will use Einstein summation notation. In this convention, repeated indices are summed over. So, we would write the product of two matrices as: c ij = a ik b kj where k is summed over, and c ij refers to the ij th element ofthematrix C. The transpose of a matrix is dened by interchanging rows and columns. We write it as then The inverse of a matrix is dened such that if so we can easily see that ~A ~x 0 = A~x A ~x 0 = ~x A A = AA = μ where μ is the identity matrix 0 0 : In the particular case where A = A, ~ then the matrix is said to be orthogonal. The adjoint is formed by taking the complex conjugate of A ~ A y = A unitary matrix is dened as one where ~A Λ A y A = and further, if A y = A, the matrix is self-adjoint or hermitian.
3 .3 Small Oscillations Now return to the problem at hand. Consider a conservative system no explicit time dependence to L, and V = V (q). We will be looking at systems with small displacements from equilibrium, so that we can linearie the EOM (as we discussed in the last topic). As we sawfrom expanding the Lagrangian, Taylor's theorem guarantees that most systems are linear if displacements are small enough. Our condition for static equilibrium states that the generalized forces Q i = io This must be true for all dimensions (no saddle points). In general we can write the Lagrangian as: L = X m ij (~q) _q i _q j V (~q) Why does the mass term depend on q? This will be the case if the co-ordinates are not Cartesian. Take the case of the double pendulum. Here we had the term cos ( ) _ _, so we have a coefcient that depends on the q's. Our stability =0 If we change co-ordinates in a singular manner: ~q! 0 =0 implying the static solution is co-ordinate independent. Lets look at small perturbations about the equilibrium position: ~q! ~q o + ~ (t) Expand the potential about the equilibrium position, ~q V (~q) =V(~q o )+ i j + h:o:t: jo The kinetic energy T = m ij _q i _q j = m ij _ i _ j 3
4 Expand the m ij mij (~q) k o k + ::::: Now T is already second order in the _ 's! since we want lowest order, keep only rst term in the m ij. Let m ij be elements of a square, real, symmetric matrix T, and be the same for a matrix V. How do we know V is symmetric? It doesn't matter which order we take the partial derivatives in, so exchange of indeces leaves the value unchanged. So, we can express the Lagrangian as o The equations of motion are: or in matrix notation: L = T ij _ i _ j V ij i j = d~ dt T d ~ V dt T ij j + V ij j =0 T + V =0 alinear ODE We have N equations, each may contain all the co-ordinates. Note: T must be positive denite, because the kinetic energy must be positive: d~ dt T d dt > 0 In matrix language, T being positive denite (having positive eigenvalues) guarantees the kinetic energy will be positive no matter what is (in other words, any matrix that obeys the above equation always is called positive denite). As for our double pendulum example, we will seek solutions for which all the co-ordinates have the same oscillation frequency the normal modes. = a k e i! kt where we should have asmany normal modes as we have degrees of freedom. So, the EOM are V = T =! T V ij a j! T ij a j =0 and we have N linear homogeneous equations for the a's. These have non-trivial solutions only if the determinant of the matrix formed by the coefcents vanishes: V! T V! T ::: V! T V! T ::: ::: ::: ::: =0 4
5 We can write this as a modied eigenvalue equation. Let! = jv Tj =0 which is a polynocial equation for. From the EOM, we get a characteristic (secular) equation for k th eigenvector a k for the k th normal mode (ie. for each root of for each eigenvalue substitute that root into the equations to get the ratio of the a i 's, and these form the components of the eigenvector. We therefore have N eigenvalues/eigenvectors. Our modied eigenvalue equation is written V a k = k T a k (called a modifed eigenvalue eqn. since usually V a = a). We can show that the k must be real. Take the transposed complex conjugate equation: a y l V = Λ l ay l T and dot it into a k a y l V a k Λ l ay l T a k =0=( k Λ l)a y l T a k and for l = k ( k Λ k) a y k T a k =0 let a k = ff k + i k a y T a k k = ~ff k T ff k + ~ k T k + i ~ff k T k ~ k T ff k Now the last term vanishes because T is a symmetric matrix. So this proves a y T a k k is real (= twice the kinetic energy so additionally is >0). For ( k Λ k) a y T a k k = 0, we must have then k Λ k =0=) k real. Now k = ay k V a k a y k T a k The denominator is T > 0, so if V is positive denite, the numerator is > 0 and k are both real and positive, so we get stable oscillations and stable equilibria. Note however that if there exists an a k such that a y k V a k < 0=) <0;! <0 and! is imaginary. The solution is then a k _ e j!jt, and we have exponential growth. So this proves that V is a minimum for stable equilibrium. ie. a y k V a k = V > 0 for >0 In other words - we took our equilibrium to be at V = 0, so for stable equilibrium any other valid state must have larger potential energy. 5
6 Suppose that all the 's are distinct (ie. there are no degenerate roots) now k Λ l 6= 0ifk6=l ( k Λ l ) ~a k T a l =0 =)~a k T a l = 0 for k 6= l We therefore have an orthogonality condition for non-degenerate roots. We can determine the ratio of elements of a k, but the normalization is undetermined. We usually normalize by requiring for each a k! we have N such equations. A = ~a k T a k = 0 and we can write the normalization condition In general, let a a ::: a n j j j j j j j j j j j j ~A T A = μ.4 Example: Two coupled oscillators C A Consider two coupled oscillators with the Lagrangian given by: L = _x + _y! o x + y + ffxy where ffxy is a coupling term. This type of Lagrangian could apply to two pendula of unit mass coupled by a spring, or electrical circuits coupled by mutual inductance. Note that if ff =0,wehave two independent oscillators with angular frequency! o. 6
7 The equations of motion are: @x =0 Welook for the normal mode solutions where x; y~e i!t :! o! x ffy =0! o! y ffx =0 we get nontrival solutions to these linear equations only if so! o! ff ff! o! =0! o! = ±ff! ; =! o ff We nd the ratios of the oscillation amplitudes for the eigenvectors by plugging the eigenfrequencies back into the linear equations:! =! o ff! A x = A y the eigenvectors are therefore:! =! o + ff! A x = A y! = p! o ff x y = p for the symmetric mode (where both pendula oscillate in the same direction), and p x! =! o + ff = p y This is a common char- If we call the eigenvector for the antisymmetric mode, where they oscillate out of phase. acteristic that the antisymmetric mode has higher frequency. solutions ; then x = + p y = p For weak coupling, variations in x and y are the superposition of oscillations with nearly identical frequencies. We get a beat frequency phenomenon. To see this, look at the limit of weak coupling where ff<<! o ;! ß! o ff 7! o
8 so in this case! ß! o + ff x _ e i!ot ff! o! o + e i!ot + = p e i!ot e i! bt + e i! bt ff! o where we let! b = ff!o, so x / e i!ot cos! b t and y _ e i!ot ff! o e i!ot = p e i!ot e i! bt e i! bt + ff! o so y / p e i!ot sin! b t So we can see that the beats are out of phase, and energy is exchanged between the two oscillators..5 Normal Co-ordinates We now illustrate the concept of normal co-ordinates, which are just the a k 's. We can expand our generalized coordinates (the 's) in terms of the eigenvectors (as we did above for x; y (t) =C k a k e i! kt where c k is a complex amplitude factor, and a k is a vector with real components. So the physical coordinates are a superposition of the normal modes. The complex amplitudes are then determined by the initial conditions, but the frequencies and phases of each mode, and the relative amplitudes of the motion of the 's for that mode are characteristic of the mode. 8
9 A transformation into the normal coordinate system diagonalizes T and V, so that we can write for the k th mode k =Re c k a k e i!tλ (note here k isn't summed over), where k is the normal co-ordinate for that mode, and L = m k _ k! k k (summed over k) in our previous example (ignoring the phase factor for now); and x / ( + ) y / ( ) x / Re [(cos! o t + i sin! o t) cos(! b t)] / cos! o t cos! b t y / Re [(cos! o t + i sin! o t)( i sin(! b t))] / sin! o t sin! b t so x = x max at t = 0, and y =0att=0. But, we can shift the phase (time) by multiplying by a complex constant..6 Degeneracy So far we have assumed that all roots of the secular equation are distinct. Nothing guarantees of course that the roots will be distinct, and we can have multiple l = k. We can then get two normal eigenvectors which will not, in general, be orthogonal. In this case we have to construct an orthogonal set. We can do this by taking a linear combination of the degenerate eigenvectors (which are linearly independent but not orthogonal). Lets look at an example for the caseof two degenerate roots corresponding to the normalized eigenvectors a ;a (remember we normalize them so that ~a T a = ). We construct a third vector a 3 = c a + c a We want to make ~a T a 3 =0. also, ~a 3 T a =c +c (~a T a )=0 c c = ~a T a ~a 3 T a 3 ==c +c +c c (~a T a )= (Note: ~a T a = ~a T a since T is symmetric). We have two independent equations that we can solve for c and c. We then use a 3 and a as our orthogonal eigenvectors. They are automatically orthogonal to the other (non-degenerate) vectors. This is called the Gram-Schmidt process. We can follow the same procedure for higher multiplicity roots also. Note that even though the 's are degenerate, the modes describe different motions of the system (which happen to have the same frequency). 9
10 .7 Molecular Vibrations Consider a diatomic molecule a physical system that can be approximated by two masses coupled by a spring: The Euler-Lagrange equations give us and we get the eigenfrequencies from the roots are L = m _x + m _x k (x x ) m ẍ +k(x x )=0 m ẍ +k(x x )=0 m! k k k m! k =0 m! k m! k k =0! =0; k μ where μ = m m m +m. What does the root! =0correspond to? It is associated with a nonoscilliatory motion of the system, ie. translation of the center of mass.! = k=μ corresponds to vibration about the center of mass. ur normal coordinates are for! =0;x =x. x =x clearly corresponds to a translation of the system. We can remove this by considering only motion about the CM, thereby reducing the number of degrees of freedom. This will help with larger numbers of atoms..7. The linear, symmetric, triatomic molecule As seen in the previous example, we would like to rst suppress the center of mass motion. To do this, we make the total momentum = 0. The CM at rest implies m x + m x + m 3 x 3 =0 x = m m (x +x 3 ) 0
11 We can use this to eliminate one or our coordinates, reducing the dimension of the problem by one L = m _x + _x 3 + m _x Λ k (x x ) +(x x 3 ) eliminate x, and do a lot of algebra, and we nd L = m M _Q a + m Q m _ s km Q a k Q s m where Q a = x + x p 3 Q s = x x p 3 and M = m + m is the total mass. These motions correspond to our normal modes which we could have guessed physically: q Our "masses/spring constant" are m =k and m M m = km m k 0 m 0 are so the normal frequencies! 0 =! (symmetric) = k m! (antismmetric) = km m m There is also a transverse mode of oscillation. Again, we reduce the degrees of freedom by removing the CM motion, but in this case also rotation: m (y + y 3 )=m y suppress translation Now we also want to suppress rotation in the plane of the molecule, so use conservation of angular momentum r a = r a0 + u a
12 where u a is the displacement from equilibrium, and the angular momentum L = X a = d dt m a r a v a ß X m a r a0 _u a X m a r a0 u a and since we can choose the origin arbitrarily, X m a r a0 u a =0 For our molecule, this implies: =0 y = y 3 (supress rotation) Now we write the Lagrangian. Let f be the deviation of from ß.» y y f = l + y 3 y l The potential energy V = k0 l f and k = k 0 l L = m We use our two constraints _y + _y 3 + m _y k0 (y + y 3 +y ) y = y 3 ; y m m = y to remove y and y 3 : L = m _y + m m _y 4 k0 y + m m y
13 From the E-L equations: so m + m ÿ m = 4k + m 0 y m! M =k 0 m m In general, if we have N atoms, we have 3N degrees of freedom. 3N-6 correspond to vibration (we subtract 6 3 for translation and 3 for rotation). For co-linear atoms, we have 3N-5 dof, since rotation about the lineof the atoms is not signicant. For a linear molecule, n- are in line, so n-4 bring the atom out of line. Here both have the same frequency ( mutually perpendicular planes). 3
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 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 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 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 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 informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
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 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 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 informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationSmall oscillations and normal modes
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
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 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 informationNormal modes. where. and. On the other hand, all such systems, if started in just the right way, will move in a simple way.
Chapter 9. Dynamics in 1D 9.4. Coupled motions in 1D 491 only the forces from the outside; the interaction forces cancel because they come in equal and opposite (action and reaction) pairs. So we get:
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 informationAdvanced Vibrations. Elements of Analytical Dynamics. By: H. Ahmadian Lecture One
Advanced Vibrations Lecture One Elements of Analytical Dynamics By: H. Ahmadian ahmadian@iust.ac.ir Elements of Analytical Dynamics Newton's laws were formulated for a single particle Can be extended to
More informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
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 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 informationQuantum Mechanics- I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras
Quantum Mechanics- I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 6 Postulates of Quantum Mechanics II (Refer Slide Time: 00:07) In my last lecture,
More informationNORMAL MODES, WAVE MOTION AND THE WAVE EQUATION. Professor G.G.Ross. Oxford University Hilary Term 2009
NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION Professor G.G.Ross Oxford University Hilary Term 009 This course of twelve lectures covers material for the paper CP4: Differential Equations, Waves and
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 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 informationLecture 7: Vectors and Matrices II Introduction to Matrices (See Sections, 3.3, 3.6, 3.7 and 3.9 in Boas)
Lecture 7: Vectors and Matrices II Introduction to Matrices (See Sections 3.3 3.6 3.7 and 3.9 in Boas) Here we will continue our discussion of vectors and their transformations. In Lecture 6 we gained
More informationLinear Algebra (part 1) : Matrices and Systems of Linear Equations (by Evan Dummit, 2016, v. 2.02)
Linear Algebra (part ) : Matrices and Systems of Linear Equations (by Evan Dummit, 206, v 202) Contents 2 Matrices and Systems of Linear Equations 2 Systems of Linear Equations 2 Elimination, Matrix Formulation
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 informationPROBLEMS In each of Problems 1 through 12:
6.5 Impulse Functions 33 which is the formal solution of the given problem. It is also possible to write y in the form 0, t < 5, y = 5 e (t 5/ sin 5 (t 5, t 5. ( The graph of Eq. ( is shown in Figure 6.5.3.
More informationAssignment 11 (C + C ) = (C + C ) = (C + C) i(c C ) ] = i(c C) (AB) = (AB) = B A = BA 0 = [A, B] = [A, B] = (AB BA) = (AB) AB
Arfken 3.4.6 Matrix C is not Hermition. But which is Hermitian. Likewise, Assignment 11 (C + C ) = (C + C ) = (C + C) [ i(c C ) ] = i(c C ) = i(c C) = i ( C C ) Arfken 3.4.9 The matrices A and B are both
More informationOPERATORS AND MEASUREMENT
Chapter OPERATORS AND MEASUREMENT In Chapter we used the results of experiments to deduce a mathematical description of the spin-/ system. The Stern-Gerlach experiments demonstrated that spin component
More informationOur goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always
October 5 Relevant reading: Section 2.1, 2.2, 2.3 and 2.4 Our goal is to solve a general constant coecient linear second order ODE a d2 y dt + bdy + cy = g (t) 2 dt where a, b, c are constants and a 0.
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 informationQuantum Field Theory Notes. Ryan D. Reece
Quantum Field Theory Notes Ryan D. Reece November 27, 2007 Chapter 1 Preliminaries 1.1 Overview of Special Relativity 1.1.1 Lorentz Boosts Searches in the later part 19th century for the coordinate transformation
More informationChapter 2. Linear Algebra. rather simple and learning them will eventually allow us to explain the strange results of
Chapter 2 Linear Algebra In this chapter, we study the formal structure that provides the background for quantum mechanics. The basic ideas of the mathematical machinery, linear algebra, are rather simple
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 information[Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty.]
Math 43 Review Notes [Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty Dot Product If v (v, v, v 3 and w (w, w, w 3, then the
More information1 Matrices and vector spaces
Matrices and vector spaces. Which of the following statements about linear vector spaces are true? Where a statement is false, give a counter-example to demonstrate this. (a) Non-singular N N matrices
More informationLab #2 - Two Degrees-of-Freedom Oscillator
Lab #2 - Two Degrees-of-Freedom Oscillator Last Updated: March 0, 2007 INTRODUCTION The system illustrated in Figure has two degrees-of-freedom. This means that two is the minimum number of coordinates
More informationMathematical Formulation of the Superposition Principle
Mathematical Formulation of the Superposition Principle Superposition add states together, get new states. Math quantity associated with states must also have this property. Vectors have this property.
More informationMath 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.
Math 350 Fall 2011 Notes about inner product spaces In this notes we state and prove some important properties of inner product spaces. First, recall the dot product on R n : if x, y R n, say x = (x 1,...,
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 informationChapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are real numbers. 1
More informationPY 351 Modern Physics Short assignment 4, Nov. 9, 2018, to be returned in class on Nov. 15.
PY 351 Modern Physics Short assignment 4, Nov. 9, 2018, to be returned in class on Nov. 15. You may write your answers on this sheet or on a separate paper. Remember to write your name on top. Please note:
More information1 Matrices and Systems of Linear Equations
Linear Algebra (part ) : Matrices and Systems of Linear Equations (by Evan Dummit, 207, v 260) Contents Matrices and Systems of Linear Equations Systems of Linear Equations Elimination, Matrix Formulation
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Revised February, 1982 NOTES ON MATRIX METHODS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Revised February, 198 NOTES ON MATRIX METHODS 1. Matrix Algebra Margenau and Murphy, The Mathematics of Physics and Chemistry, Chapter 10, give almost
More informationLinear Second-Order Differential Equations LINEAR SECOND-ORDER DIFFERENTIAL EQUATIONS
11.11 LINEAR SECOND-ORDER DIFFERENTIAL EQUATIONS A Spring with Friction: Damped Oscillations The differential equation, which we used to describe the motion of a spring, disregards friction. But there
More information18.06 Problem Set 8 - Solutions Due Wednesday, 14 November 2007 at 4 pm in
806 Problem Set 8 - Solutions Due Wednesday, 4 November 2007 at 4 pm in 2-06 08 03 Problem : 205+5+5+5 Consider the matrix A 02 07 a Check that A is a positive Markov matrix, and find its steady state
More informationElementary maths for GMT
Elementary maths for GMT Linear Algebra Part 2: Matrices, Elimination and Determinant m n matrices The system of m linear equations in n variables x 1, x 2,, x n a 11 x 1 + a 12 x 2 + + a 1n x n = b 1
More informationConceptual Questions for Review
Conceptual Questions for Review Chapter 1 1.1 Which vectors are linear combinations of v = (3, 1) and w = (4, 3)? 1.2 Compare the dot product of v = (3, 1) and w = (4, 3) to the product of their lengths.
More informationMATH 583A REVIEW SESSION #1
MATH 583A REVIEW SESSION #1 BOJAN DURICKOVIC 1. Vector Spaces Very quick review of the basic linear algebra concepts (see any linear algebra textbook): (finite dimensional) vector space (or linear space),
More information22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes
Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration
More informationProblem set 2. Math 212b February 8, 2001 due Feb. 27
Problem set 2 Math 212b February 8, 2001 due Feb. 27 Contents 1 The L 2 Euler operator 1 2 Symplectic vector spaces. 2 2.1 Special kinds of subspaces....................... 3 2.2 Normal forms..............................
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 informationComputational Methods CMSC/AMSC/MAPL 460. Eigenvalues and Eigenvectors. Ramani Duraiswami, Dept. of Computer Science
Computational Methods CMSC/AMSC/MAPL 460 Eigenvalues and Eigenvectors Ramani Duraiswami, Dept. of Computer Science Eigen Values of a Matrix Recap: A N N matrix A has an eigenvector x (non-zero) with corresponding
More information1 Fundamental physical postulates. C/CS/Phys C191 Quantum Mechanics in a Nutshell I 10/04/07 Fall 2007 Lecture 12
C/CS/Phys C191 Quantum Mechanics in a Nutshell I 10/04/07 Fall 2007 Lecture 12 In this and the next lecture we summarize the essential physical and mathematical aspects of quantum mechanics relevant to
More informationAN INTRODUCTION TO LAGRANGE EQUATIONS. Professor J. Kim Vandiver October 28, 2016
AN INTRODUCTION TO LAGRANGE EQUATIONS Professor J. Kim Vandiver October 28, 2016 kimv@mit.edu 1.0 INTRODUCTION This paper is intended as a minimal introduction to the application of Lagrange equations
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More informationQuantum Mechanics- I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras
Quantum Mechanics- I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 4 Postulates of Quantum Mechanics I In today s lecture I will essentially be talking
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2018 8 Lecture 8 8.1 Matrices July 22, 2018 We shall study
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of n-dof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More informationPHYS 705: Classical Mechanics. Rigid Body Motion Introduction + Math Review
1 PHYS 705: Classical Mechanics Rigid Body Motion Introduction + Math Review 2 How to describe a rigid body? Rigid Body - a system of point particles fixed in space i r ij j subject to a holonomic constraint:
More informationconventions and notation
Ph95a lecture notes, //0 The Bloch Equations A quick review of spin- conventions and notation The quantum state of a spin- particle is represented by a vector in a two-dimensional complex Hilbert space
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
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 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 informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationMath Linear Algebra II. 1. Inner Products and Norms
Math 342 - Linear Algebra II Notes 1. Inner Products and Norms One knows from a basic introduction to vectors in R n Math 254 at OSU) that the length of a vector x = x 1 x 2... x n ) T R n, denoted x,
More informationJim Lambers MAT 610 Summer Session Lecture 1 Notes
Jim Lambers MAT 60 Summer Session 2009-0 Lecture Notes Introduction This course is about numerical linear algebra, which is the study of the approximate solution of fundamental problems from linear algebra
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 informationSecond quantization: where quantization and particles come from?
110 Phys460.nb 7 Second quantization: where quantization and particles come from? 7.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system? 7.1.1.Lagrangian Lagrangian
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationSection 9 Variational Method. Page 492
Section 9 Variational Method Page 492 Page 493 Lecture 27: The Variational Method Date Given: 2008/12/03 Date Revised: 2008/12/03 Derivation Section 9.1 Variational Method: Derivation Page 494 Motivation
More informationVector Space Basics. 1 Abstract Vector Spaces. 1. (commutativity of vector addition) u + v = v + u. 2. (associativity of vector addition)
Vector Space Basics (Remark: these notes are highly formal and may be a useful reference to some students however I am also posting Ray Heitmann's notes to Canvas for students interested in a direct computational
More informationIntroduction to Group Theory
Chapter 10 Introduction to Group Theory Since symmetries described by groups play such an important role in modern physics, we will take a little time to introduce the basic structure (as seen by a physicist)
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 informationPhysics 70007, Fall 2009 Answers to Final Exam
Physics 70007, Fall 009 Answers to Final Exam December 17, 009 1. Quantum mechanical pictures a Demonstrate that if the commutation relation [A, B] ic is valid in any of the three Schrodinger, Heisenberg,
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
More information1. Matrix multiplication and Pauli Matrices: Pauli matrices are the 2 2 matrices. 1 0 i 0. 0 i
Problems in basic linear algebra Science Academies Lecture Workshop at PSGRK College Coimbatore, June 22-24, 2016 Govind S. Krishnaswami, Chennai Mathematical Institute http://www.cmi.ac.in/~govind/teaching,
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. Linear Equations and Matrices. Copyright 2005, W.R. Winfrey
Copyright 2005, W.R. Winfrey Topics Preliminaries Systems of Linear Equations Matrices Algebraic Properties of Matrix Operations Special Types of Matrices and Partitioned Matrices Matrix Transformations
More informationIntroduction to Determinants
Introduction to Determinants For any square matrix of order 2, we have found a necessary and sufficient condition for invertibility. Indeed, consider the matrix The matrix A is invertible if and only if.
More informationCHAPTER 3. Matrix Eigenvalue Problems
A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 3 A COLLECTION OF HANDOUTS ON SYSTEMS OF ORDINARY DIFFERENTIAL
More informationLinear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space
Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................
More informationDiagonalization by a unitary similarity transformation
Physics 116A Winter 2011 Diagonalization by a unitary similarity transformation In these notes, we will always assume that the vector space V is a complex n-dimensional space 1 Introduction A semi-simple
More informationAnalytical Mechanics for Relativity and Quantum Mechanics
Analytical Mechanics for Relativity and Quantum Mechanics Oliver Davis Johns San Francisco State University OXPORD UNIVERSITY PRESS CONTENTS Dedication Preface Acknowledgments v vii ix PART I INTRODUCTION:
More informationMath 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Math 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008 Exam 2 will be held on Tuesday, April 8, 7-8pm in 117 MacMillan What will be covered The exam will cover material from the lectures
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationSAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear Algebra
1.1. Introduction SAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear algebra is a specific branch of mathematics dealing with the study of vectors, vector spaces with functions that
More informationFREE VIBRATION RESPONSE OF UNDAMPED SYSTEMS
Lecture Notes: STRUCTURAL DYNAMICS / FALL 2011 / Page: 1 FREE VIBRATION RESPONSE OF UNDAMPED SYSTEMS : : 0, 0 As demonstrated previously, the above Equation of Motion (free-vibration equation) has a solution
More informationLecture 11. Linear systems: Cholesky method. Eigensystems: Terminology. Jacobi transformations QR transformation
Lecture Cholesky method QR decomposition Terminology Linear systems: Eigensystems: Jacobi transformations QR transformation Cholesky method: For a symmetric positive definite matrix, one can do an LU decomposition
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 informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
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 informationQM and Angular Momentum
Chapter 5 QM and Angular Momentum 5. Angular Momentum Operators In your Introductory Quantum Mechanics (QM) course you learned about the basic properties of low spin systems. Here we want to review that
More informationChapter 7. Homogeneous equations with constant coefficients
Chapter 7. Homogeneous equations with constant coefficients It has already been remarked that we can write down a formula for the general solution of any linear second differential equation y + a(t)y +
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 informationSAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear Algebra
SAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to 1.1. Introduction Linear algebra is a specific branch of mathematics dealing with the study of vectors, vector spaces with functions that
More information5.1 Classical Harmonic Oscillator
Chapter 5 Harmonic Oscillator 5.1 Classical Harmonic Oscillator m l o l Hooke s Law give the force exerting on the mass as: f = k(l l o ) where l o is the equilibrium length of the spring and k is the
More information