Multiple Degree of Freedom Systems. The Millennium bridge required many degrees of freedom to model and design with.
|
|
- Gervais Arnold
- 5 years ago
- Views:
Transcription
1 Multiple Degree of Freedom Systems The Millennium bridge required many degrees of freedom to model and design with.
2 The first step in analyzing multiple degrees of freedom (DOF) is to look at DOF DOF: Minimum number of coordinates to specify the position of a system Many systems have more than DOF Examples of DOF systems car with sprung and unsprung mass (both heave) elastic pendulum (radial and angular) motions of a ship (roll and pitch) D.J. Inman /58 Mechanical Engineering at Virginia Tech
3 4. Two-Degree-of-Freedom Model (Undamped) A degree of freedom system used to base much of the analysis and conceptual development of MDOF systems on.
4 Free-Body Diagram of each mass Figure 4. k (x -x ) k x m m x x
5 Summing forces yields the equations of motion: m x ( t) k x ( t) k x ( t) x ( t) m x ( t) k x ( t) x ( t) (4.) Rearranging terms: m x ( t) ( k k ) x ( t) k x ( t) 0 m x ( t) k x ( t) k x ( t) 0 (4.)
6 Note that it is always the case that A Degree-of-Freedom system has Two equations of motion! Two natural frequencies (as we shall see)!
7 The dynamics of a DOF system consists of homogeneous and coupled equations Free vibrations, so homogeneous eqs. Equations are coupled: Both have x and x. If only one mass moves, the other follows Example: pitch and heave of a car model In this case the coupling is due to k. Mathematically and Physically If k = 0, no coupling occurs and can be solved as two independent SDOF systems
8 Initial Conditions Two coupled, second -order, ordinary differential equations with constant coefficients Needs 4 constants of integration to solve Thus 4 initial conditions on positions and velocities x (0) x, x (0) x, x (0) x, x (0) x
9 Solution by Matrix Methods The two equations can be written in the form of a single matrix equation (see pages 7-75 if matrices and vectors are a struggle for you) : x( t) x( t) x( t) x( t ), ( t ), ( t) x ( t) x x ( t) x x ( t) (4.4), (4.5) m 0 k k k M, K 0 m k k (4.6), (4.9) Mx Kx 0 m x m x ( t) ( k ( t) k k x ) x ( t) k ( t) k x x ( t) 0 ( t) 0
10 Initial Conditions IC s can also be written in vector form x x (0), and x(0) x x 0 0 x 0 0
11 The approach to a Solution: For DOF we assumed the scalar solution ae lt Similarly, now we assume the vector form: Let x( t) ue jt (4.5) j, u 0,, u unknown - M K ue 0 jt - M K u 0 (4.6) (4.7)
12 This changes the differential equation of motion into algebraic vector equation: - M K u0 (4.7) This is two algebraic equation in 3 uknowns ( vector of two elements and scalar): u u u =, and
13 The condition for solution of this matrix equation requires that the the matrix inverse does not exist: If the inv - M K exists u 0: which is the static equilibrium position. For motion to occur u 0 - M K does not exist or det - M K 0 (4.9) The determinant results in equation in one unknown (called the characteristic equation)
14 Back to our specific system: the characteristic equation is defined as det- M K 0 det m k k k k 0 m k (4.0) m m 4 (m k m k m k ) k k 0 (4.) Eq. (4.) is quadratic in so four solutions result: and and
15 Once is known, use equation (4.7) again to calculate the corresponding vectors u and u This yields vector equation for each squared frequency: and ( M K) u 0 (4.) ( M K) u 0 (4.3) Each of these matrix equations represents equations in the unknowns components of the vector, but the coefficient matrix is singular so each matrix equation results in only independent equation. The following examples clarify this.
16 Examples 4..5 & 4..6:calculating u and m =9 kg,m =kg, k =4 N/m and k =3 N/m The characteristic equation becomes =( -)( -4)=0 = and =4 or,3 rad/s,,4 rad/s Each value of yields an expression or u:
17 Computing the vectors u For =, denote u then we have u (- M K) u 0 u 7 9() 3 u () u 0 9u 3u 0 and 3u u 0 equations, unknowns but DEPENDENT! (the nd equation is -3 times the first)
18 Only the direction of vectors u can be determined, not the magnitude as it remains arbitrary u u u u 3 3 results from both equations: only the direction, not the magnitude can be determined! This is because: det( M K) 0. The magnitude of the vector is arbitrary. To see this suppose that u satisfies ( M K) u 0, so does au, a arbitrary. So ( M K) au 0 ( M K) u 0
19 Likewise for the second value of : For = 4, let u then we have u (- M K) u 0 u 7 9(4) 3 u (4) u 0 9u 3u 0 or u u 3 Note that the other equation is the same
20 What to do about the magnitude! Several possibilities, here we just fix one element: Choose: Choose: u u u u 3 3
21 Thus the solution to the algebraic matrix equation is:,3, has mode shape u 3,4, has mode shape u Here we have introduce the name mode shape to describe the vectors u and u. The origin of this name comes later 3
22 Return now to the time response: We have computed four solutions: x( t) u e, u e, u e, u e jt jt jt jt Since linear, we can combine as: x() t au e bu e cu e du e jt jt jt jt j t jt jt jt x() t ae be u ce de u A sin( t ) u A sin( t ) u (4.4) (4.6) where A, A,, and are constants of integration determined by initial conditions. Note that to go from the exponential form to to sine requires Euler s formula for trig functions and uses up the +/- sign on omega
23 Physical interpretation of all that math! Each of the TWO masses is oscillating at TWO natural frequencies and The relative magnitude of each sine term, and hence of the magnitude of oscillation of m and m is determined by the value of A u and A u The vectors u and u are called mode shapes because the describe the relative magnitude of oscillation between the two masses
24 What is a mode shape? First note that A, A, and are determined by the initial conditions Choose them so that A = = =0 Then: x(t) x (t) x (t) A Thus each mass oscillates at (one) frequency with magnitudes proportional to u the st mode shape u u sin t A u sin t
25 A graphic look at mode shapes: If IC s correspond to mode or, then the response is purely in mode or mode. Mode : x k m x =A/3 k x x =A m u 3 Mode : x k m x =-A/3 k x x =A m u 3
26 Example 4..7 given the initial conditions compute the time response cos cos cos 3 cos 3 ) ( ) ( sin sin sin 3 sin 3 ) ( ) ( 0 0 (0) mm, 0 = (0) consider t A t A t A t A t x t x t A t A t A t A t x t x x x
27 cos cos cos 3 cos sin sin sin 3 sin 3 0 mm A A A A A A A A At t=0 we have
28 4 equations in 4 unknowns: 3 A sin A sin 0 A sin A sin 0 A cos A cos 0 A cos A cos Yields: A.5 mm, A.5 mm, rad
29 The final solution is: x (t) 0.5cos t 0.5cost x (t).5cos t.5cost (4.34) These initial conditions gives a response that is a combination of modes. Both harmonic, but their summation is not. Figure 4.3
30 Solution as a sum of modes x(t) a u cos t a u cos t Determines how the first frequency contributes to the response Determines how the second frequency contributes to the response
31 Things to note Two degrees of freedom implies two natural frequencies Each mass oscillates at with these two frequencies present in the response and beats could result Frequencies are not those of two component systems k m.63, k m.73 The above is not the most efficient way to calculate frequencies as the following describes
32 Some matrix and vector reminders A T xx a b d b c c ad cb c a A x x m 0 M x Mx m x m x 0 T m T M 0 x Mx 0 for every value of x except 0 Then M is said to be positive definite
33 4. Eigenvalues and Natural Frequencies Can connect the vibration problem with the algebraic eigenvalue problem developed in math This will give us some powerful computational skills And some powerful theory All the codes have eigen-solvers so these painful calculations can be automated
34 Some matrix results to help us use available computational tools: A matrix M is defined to be symmetric if M M T A symmetric matrix M is positive definite if x T Mx 0 for all nonzero vectors x A symmetric positive definite matrix M can be factored M LL T Here L is upper triangular, called a Cholesky matrix
35 If the matrix L is diagonal, it defines the matrix square root The matrix square root is the matrix M / such that M / M / M If M is diagonal, then the matrix square root is just the root of the diagonal elements: m L M / 0 (4.35) 0 m
36 A change of coordinates is introduced to capitalize on existing mathematics For a symmetric, positive definite matrix M: m 0 m 0 M, M, M / 0 0 m m 0 m / / Let x( t) M q( t) and multiply by M : / / / / M MM t M KM t q( ) q( ) 0 (4.38) I identity K symmetric or q( t) Kq( t) 0 where K M KM / / k K is called the mass normalized stiffness and is similar to the scalar m used extensively in single degree of freedom analysis. The key here is that K is a SYMMETRIC matrix allowing the use of many nice properties and computational tools m 0
37 How the vibration problem relates to the real symmetric eigenvalue problem Assume q( t) ve in q( t) Kq( t) 0 jt ve Kve 0, v 0 or jt jt Kv v Kv lv v 0 vibration problem real symmetric eigenvalue problem (4.40) (4.4) Note that the martrix as does n K contains the same type of information in the single degree of freedom case.
38 Important Properties of the n x n Real Symmetric Eigenvalue Problem There are n eigenvalues and they are all real valued There are n eigenvectors and they are all real valued The set of eigenvectors are orthogonal The set of eigenvectors are linearly independent The matrix is similar to a diagonal matrix
39 Square Matrix Review Let a ik be the ik th element of A then A is symmetric if a ik = a ki denoted A T =A A is positive definite if x T Ax > 0 for all nonzero x (also implies each l i > 0) The stiffness matrix is usually symmetric and positive semi definite (could have a zero eigenvalue) The mass matrix is positive definite and symmetric (and so far, its diagonal)
40 Normal and orthogonal vectors x x M, y x n y M y n n, inner product is x T y x i y i i x orthogonal to y if x T y 0 x is normal if x T x if a the set of vectores is is both orthogonal and normal it is called an orthonormal set The norm of x is x x T x (4.43)
41 Normalizing any vector can be done by dividing it by its norm: x x T x has norm of (4.44) To see this compute x x T x T x x T x x x T x x T x x T x
42 Examples 4.. through 4..4 / / K M KM so K which is symmetric. 3 3-l l det( K li) det l 6l 8 0 which has roots: l and l 4
43 ( K l I) v 0 3 v 0 3 v 0 v v v v 0 v ( ) The first normalized eigenvector
44 Likewise the second normalized eigenvector is computed and shown to be orthogonal to the first, so that the set is orthonormal v, v T v ( ) 0 v T v ( ) v T v ( ()()) v i are orthonormal
45 Modes u and Eigenvectors v are different but related: u v and u v x M / q u M / v (4.37) Note M / u v
46 This orthonormal set of vectors is used to form an Orthogonal Matrix P T P P v v called a matrix of eigenvectors (normalized) T T v v v 0 v I T T 0 v v v v P is called an orthogonal matrix P KP P Kv Kv P lv l v T T T T T lv v lv v l 0 diag(, ) T T l 0 l (4.47) v v lv v P is also called a modal matrix
47 Example 4..3 compute P and show that it is an orthogonal matrix From the previous example: P v v P T P 0 0 I
48 Example 4..4 Compute the square of the frequencies by matrix manipulation T P KP rad/s and rad/s In general: T P KP diag l diag( ) (4.48) i i
49 Example 4..5 The equations of motion: m x ( k k ) x k x 0 m x k x ( k k ) x 0 3 (4.49) Figure 4.4 In matrix form these become: m 0 0 m x k k k k k k 3 x 0 (4.50)
50 Next substitute numerical values and compute P and m kg, m 4 kg, k k 3 0 N/m and k = N/m 0 M, K 0 4 / / K M KM det K li det l 5l 35 0 l l l.890 and l rad/s and.098 rad/s
51 Next compute the eigenvectors For l equation (4.4 ) becomes: Normalizing v v 9.089v yields v v v v (9.089) v v v v v 0. 09, and v , likewise v
52 Next check the value of P to see if it behaves as its suppose to: P v v T P KP T PP Yes!
53 A note on eigenvectors In the previous section, we could have chosed v to be v instead of v because one can always multiple an eigenvector by a constant and if the constant is - the result is still a normalized vector. Does this make any difference? No! Try it in the previous example
54 All of the previous examples can and should be solved by hand to learn the methods However, they can also be solved on calculators with matrix functions and with the codes listed in the last section In fact, for more then two DOF one must use a code to solve for the natural frequencies and mode shapes. Next we examine 3 other formulations for solving for modal data
55 Matlab commands To compute the inverse of the square matrix A: inv(a) or use A\eye(n) where n is the size of the matrix [P,D]=eig(A) computes the eigenvalues and normalized eigenvectors (watch the order). Stores them in the eigenvector matrix P and the diagonal matrix D (D=)
56 More commands To compute the matrix square root use sqrtm(a) To compute the Cholesky factor: L= chol(m) To compute the norm: norm(x) To compute the determinant det(a) To enter a matrix: K=[7-3;-3 3]; M=[9 0;0 ]; To multiply: K*inv(chol(M))
57 An alternate approach to normalizing mode shapes From equation (4.7) M K u 0, u 0 Now scale the mode shapes by computing such that i u i T M i u i i u i T u i w u i i i is called mass normalized and it satisfies: Mw Kw 0 w Kw, i, T i i i i i i (4.53)
58 There are 3 approaches to computing mode shapes and frequencies (i) Mu Ku (ii) u M Ku (iii) v M KM v (i) Is the Generalized Symmetric Eigenvalue Problem easy for hand computations, inefficient for computers (ii) Is the Asymmetric Eigenvalue Problem very expensive computationally (iii) Is the Symmetric Eigenvalue Problem the cheapest computationally
3. Mathematical Properties of MDOF Systems
3. Mathematical Properties of MDOF Systems 3.1 The Generalized Eigenvalue Problem Recall that the natural frequencies ω and modes a are found from [ - ω 2 M + K ] a = 0 or K a = ω 2 M a Where M and K are
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 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 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 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 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 informationDamping Matrix. Donkey2Ft
1 Damping Matrix DonkeyFt Damping in a single-degree-of-freedom (SDOF) system is well studied. Whether the system is under-damped, over-damped, or critically damped is well known. For an under-damped system,
More informationReduction to the associated homogeneous system via a particular solution
June PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 5) Linear Algebra This study guide describes briefly the course materials to be covered in MA 5. In order to be qualified for the credit, one
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 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 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 informationBasic Linear Algebra in MATLAB
Basic Linear Algebra in MATLAB 9.29 Optional Lecture 2 In the last optional lecture we learned the the basic type in MATLAB is a matrix of double precision floating point numbers. You learned a number
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 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.6: 16, 20, 24, 27, 28
.6: 6, 2, 24, 27, 28 6) If A is positive definite, then A is positive definite. The proof of the above statement can easily be shown for the following 2 2 matrix, a b A = b c If that matrix is positive
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 informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
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 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 informationMATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.
MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation
More informationMTH 2032 SemesterII
MTH 202 SemesterII 2010-11 Linear Algebra Worked Examples Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education December 28, 2011 ii Contents Table of Contents
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 informationand let s calculate the image of some vectors under the transformation T.
Chapter 5 Eigenvalues and Eigenvectors 5. Eigenvalues and Eigenvectors Let T : R n R n be a linear transformation. Then T can be represented by a matrix (the standard matrix), and we can write T ( v) =
More information1. Select the unique answer (choice) for each problem. Write only the answer.
MATH 5 Practice Problem Set Spring 7. Select the unique answer (choice) for each problem. Write only the answer. () Determine all the values of a for which the system has infinitely many solutions: x +
More informationYORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics MATH M Test #1. July 11, 2013 Solutions
YORK UNIVERSITY Faculty of Science Department of Mathematics and Statistics MATH 222 3. M Test # July, 23 Solutions. For each statement indicate whether it is always TRUE or sometimes FALSE. Note: For
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 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 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 informationDynamics of Structures
Dynamics of Structures Elements of structural dynamics Roberto Tomasi 11.05.2017 Roberto Tomasi Dynamics of Structures 11.05.2017 1 / 22 Overview 1 SDOF system SDOF system Equation of motion Response spectrum
More informationMATH 304 Linear Algebra Lecture 34: Review for Test 2.
MATH 304 Linear Algebra Lecture 34: Review for Test 2. Topics for Test 2 Linear transformations (Leon 4.1 4.3) Matrix transformations Matrix of a linear mapping Similar matrices Orthogonality (Leon 5.1
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 informationSHOCK RESPONSE OF MULTI-DEGREE-OF-FREEDOM SYSTEMS Revision F By Tom Irvine May 24, 2010
SHOCK RESPONSE OF MULTI-DEGREE-OF-FREEDOM SYSTEMS Revision F By Tom Irvine Email: tomirvine@aol.com May 4, 010 Introduction The primary purpose of this tutorial is to present the Modal Transient method
More informationIntroduction to Matrices
POLS 704 Introduction to Matrices Introduction to Matrices. The Cast of Characters A matrix is a rectangular array (i.e., a table) of numbers. For example, 2 3 X 4 5 6 (4 3) 7 8 9 0 0 0 Thismatrix,with4rowsand3columns,isoforder
More information11.2 Basic First-order System Methods
112 Basic First-order System Methods 797 112 Basic First-order System Methods Solving 2 2 Systems It is shown here that any constant linear system u a b = A u, A = c d can be solved by one of the following
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 informationMATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors.
MATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors. Orthogonal sets Let V be a vector space with an inner product. Definition. Nonzero vectors v 1,v
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 informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationRemark 1 By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.
Sec 5 Eigenvectors and Eigenvalues In this chapter, vector means column vector Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called
More informationFinal Review Written by Victoria Kala SH 6432u Office Hours R 12:30 1:30pm Last Updated 11/30/2015
Final Review Written by Victoria Kala vtkala@mathucsbedu SH 6432u Office Hours R 12:30 1:30pm Last Updated 11/30/2015 Summary This review contains notes on sections 44 47, 51 53, 61, 62, 65 For your final,
More informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
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 informationChap 3. Linear Algebra
Chap 3. Linear Algebra Outlines 1. Introduction 2. Basis, Representation, and Orthonormalization 3. Linear Algebraic Equations 4. Similarity Transformation 5. Diagonal Form and Jordan Form 6. Functions
More informationLast week we presented the following expression for the angles between two vectors u and v in R n ( ) θ = cos 1 u v
Orthogonality (I) Last week we presented the following expression for the angles between two vectors u and v in R n ( ) θ = cos 1 u v u v which brings us to the fact that θ = π/2 u v = 0. Definition (Orthogonality).
More informationCheat Sheet for MATH461
Cheat Sheet for MATH46 Here is the stuff you really need to remember for the exams Linear systems Ax = b Problem: We consider a linear system of m equations for n unknowns x,,x n : For a given matrix A
More informationEigenvalues and Eigenvectors
Contents Eigenvalues and Eigenvectors. Basic Concepts. Applications of Eigenvalues and Eigenvectors 8.3 Repeated Eigenvalues and Symmetric Matrices 3.4 Numerical Determination of Eigenvalues and Eigenvectors
More informationMore chapter 3...linear dependence and independence... vectors
More chapter 3...linear dependence and independence... vectors It is important to determine if a set of vectors is linearly dependent or independent Consider a set of vectors A, B, and C. If we can find
More informationTopic 15 Notes Jeremy Orloff
Topic 5 Notes Jeremy Orloff 5 Transpose, Inverse, Determinant 5. Goals. Know the definition and be able to compute the inverse of any square matrix using row operations. 2. Know the properties of inverses.
More informationlinearly indepedent eigenvectors as the multiplicity of the root, but in general there may be no more than one. For further discussion, assume matrice
3. Eigenvalues and Eigenvectors, Spectral Representation 3.. Eigenvalues and Eigenvectors A vector ' is eigenvector of a matrix K, if K' is parallel to ' and ' 6, i.e., K' k' k is the eigenvalue. If is
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 2018 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
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 informationRemark By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.
Sec 6 Eigenvalues and Eigenvectors Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called an eigenvalue of A if there is a nontrivial
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 informationLecture Note 1: Background
ECE5463: Introduction to Robotics Lecture Note 1: Background Prof. Wei Zhang Department of Electrical and Computer Engineering Ohio State University Columbus, Ohio, USA Spring 2018 Lecture 1 (ECE5463 Sp18)
More informationPolar Form of a Complex Number (I)
Polar Form of a Complex Number (I) The polar form of a complex number z = x + iy is given by z = r (cos θ + i sin θ) where θ is the angle from the real axes to the vector z and r is the modulus of z, i.e.
More informationMatrix Iteration. Giacomo Boffi.
http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 12, 2016 Outline Second -Ritz Method Dynamic analysis of MDOF
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 information2. Every linear system with the same number of equations as unknowns has a unique solution.
1. For matrices A, B, C, A + B = A + C if and only if A = B. 2. Every linear system with the same number of equations as unknowns has a unique solution. 3. Every linear system with the same number of equations
More informationFinal Exam, Linear Algebra, Fall, 2003, W. Stephen Wilson
Final Exam, Linear Algebra, Fall, 2003, W. Stephen Wilson Name: TA Name and section: NO CALCULATORS, SHOW ALL WORK, NO OTHER PAPERS ON DESK. There is very little actual work to be done on this exam if
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 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 informationEigenvalues, Eigenvectors and the Jordan Form
EE/ME 701: Advanced Linear Systems Eigenvalues, Eigenvectors and the Jordan Form Contents 1 Introduction 3 1.1 Review of basic facts about eigenvectors and eigenvalues..... 3 1.1.1 Looking at eigenvalues
More information1 Last time: least-squares problems
MATH Linear algebra (Fall 07) Lecture Last time: least-squares problems Definition. If A is an m n matrix and b R m, then a least-squares solution to the linear system Ax = b is a vector x R n such that
More informationQuiz ) Locate your 1 st order neighbors. 1) Simplify. Name Hometown. Name Hometown. Name Hometown.
Quiz 1) Simplify 9999 999 9999 998 9999 998 2) Locate your 1 st order neighbors Name Hometown Me Name Hometown Name Hometown Name Hometown Solving Linear Algebraic Equa3ons Basic Concepts Here only real
More informationHomework 3 Solutions Math 309, Fall 2015
Homework 3 Solutions Math 39, Fall 25 782 One easily checks that the only eigenvalue of the coefficient matrix is λ To find the associated eigenvector, we have 4 2 v v 8 4 (up to scalar multiplication)
More informationEIGENVALUES AND EIGENVECTORS 3
EIGENVALUES AND EIGENVECTORS 3 1. Motivation 1.1. Diagonal matrices. Perhaps the simplest type of linear transformations are those whose matrix is diagonal (in some basis). Consider for example the matrices
More informationMath 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination
Math 0, Winter 07 Final Exam Review Chapter. Matrices and Gaussian Elimination { x + x =,. Different forms of a system of linear equations. Example: The x + 4x = 4. [ ] [ ] [ ] vector form (or the column
More information(Mathematical Operations with Arrays) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Mathematical Operations with Arrays) Contents Getting Started Matrices Creating Arrays Linear equations Mathematical Operations with Arrays Using Script
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 informationChapter 4 & 5: Vector Spaces & Linear Transformations
Chapter 4 & 5: Vector Spaces & Linear Transformations Philip Gressman University of Pennsylvania Philip Gressman Math 240 002 2014C: Chapters 4 & 5 1 / 40 Objective The purpose of Chapter 4 is to think
More informationJUST THE MATHS UNIT NUMBER 9.8. MATRICES 8 (Characteristic properties) & (Similarity transformations) A.J.Hobson
JUST THE MATHS UNIT NUMBER 9.8 MATRICES 8 (Characteristic properties) & (Similarity transformations) by A.J.Hobson 9.8. Properties of eigenvalues and eigenvectors 9.8. Similar matrices 9.8.3 Exercises
More informationMTH Linear Algebra. Study Guide. Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education
MTH 3 Linear Algebra Study Guide Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education June 3, ii Contents Table of Contents iii Matrix Algebra. Real Life
More informationProperties of Matrices and Operations on Matrices
Properties of Matrices and Operations on Matrices A common data structure for statistical analysis is a rectangular array or matris. Rows represent individual observational units, or just observations,
More informationNOTES ON LINEAR ALGEBRA CLASS HANDOUT
NOTES ON LINEAR ALGEBRA CLASS HANDOUT ANTHONY S. MAIDA CONTENTS 1. Introduction 2 2. Basis Vectors 2 3. Linear Transformations 2 3.1. Example: Rotation Transformation 3 4. Matrix Multiplication and Function
More informationMath 4A Notes. Written by Victoria Kala Last updated June 11, 2017
Math 4A Notes Written by Victoria Kala vtkala@math.ucsb.edu Last updated June 11, 2017 Systems of Linear Equations A linear equation is an equation that can be written in the form a 1 x 1 + a 2 x 2 +...
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 informationReview of Linear Algebra
Review of Linear Algebra Dr Gerhard Roth COMP 40A Winter 05 Version Linear algebra Is an important area of mathematics It is the basis of computer vision Is very widely taught, and there are many resources
More informationTheory of Vibrations in Stewart Platforms
Theory of Vibrations in Stewart Platforms J.M. Selig and X. Ding School of Computing, Info. Sys. & Maths. South Bank University London SE1 0AA, U.K. (seligjm@sbu.ac.uk) Abstract This article develops a
More informationMidterm 1 NAME: QUESTION 1 / 10 QUESTION 2 / 10 QUESTION 3 / 10 QUESTION 4 / 10 QUESTION 5 / 10 QUESTION 6 / 10 QUESTION 7 / 10 QUESTION 8 / 10
Midterm 1 NAME: RULES: You will be given the entire period (1PM-3:10PM) to complete the test. You can use one 3x5 notecard for formulas. There are no calculators nor those fancy cellular phones nor groupwork
More information1. What is the determinant of the following matrix? a 1 a 2 4a 3 2a 2 b 1 b 2 4b 3 2b c 1. = 4, then det
What is the determinant of the following matrix? 3 4 3 4 3 4 4 3 A 0 B 8 C 55 D 0 E 60 If det a a a 3 b b b 3 c c c 3 = 4, then det a a 4a 3 a b b 4b 3 b c c c 3 c = A 8 B 6 C 4 D E 3 Let A be an n n matrix
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 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 information( ) Chapter 3: Free Vibration of the Breakwater. 3.1 Introduction
Chapter : Free Vibration of the Breakwater. Introduction Before conducting forcing analyses on the breakwater, a free vibration study is necessary. This chapter will describe the analysis conducted to
More informationVAR Model. (k-variate) VAR(p) model (in the Reduced Form): Y t-2. Y t-1 = A + B 1. Y t + B 2. Y t-p. + ε t. + + B p. where:
VAR Model (k-variate VAR(p model (in the Reduced Form: where: Y t = A + B 1 Y t-1 + B 2 Y t-2 + + B p Y t-p + ε t Y t = (y 1t, y 2t,, y kt : a (k x 1 vector of time series variables A: a (k x 1 vector
More informationLinear Algebra in Actuarial Science: Slides to the lecture
Linear Algebra in Actuarial Science: Slides to the lecture Fall Semester 2010/2011 Linear Algebra is a Tool-Box Linear Equation Systems Discretization of differential equations: solving linear equations
More informationStudy Guide for Linear Algebra Exam 2
Study Guide for Linear Algebra Exam 2 Term Vector Space Definition A Vector Space is a nonempty set V of objects, on which are defined two operations, called addition and multiplication by scalars (real
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 informationB553 Lecture 5: Matrix Algebra Review
B553 Lecture 5: Matrix Algebra Review Kris Hauser January 19, 2012 We have seen in prior lectures how vectors represent points in R n and gradients of functions. Matrices represent linear transformations
More information14 Singular Value Decomposition
14 Singular Value Decomposition For any high-dimensional data analysis, one s first thought should often be: can I use an SVD? The singular value decomposition is an invaluable analysis tool for dealing
More informationMath (P)Review Part I:
Lecture 1: Math (P)Review Part I: Linear Algebra Computer Graphics CMU 15-462/15-662, Fall 2017 Homework 0.0 (Due Monday!) Exercises will be a bit harder / more rigorous than what you will do for the rest
More informationMATH.2720 Introduction to Programming with MATLAB Vector and Matrix Algebra
MATH.2720 Introduction to Programming with MATLAB Vector and Matrix Algebra A. Vectors A vector is a quantity that has both magnitude and direction, like velocity. The location of a vector is irrelevant;
More informationMAT1302F Mathematical Methods II Lecture 19
MAT302F Mathematical Methods II Lecture 9 Aaron Christie 2 April 205 Eigenvectors, Eigenvalues, and Diagonalization Now that the basic theory of eigenvalues and eigenvectors is in place most importantly
More informationSTRUCTURAL DYNAMICS BASICS:
BASICS: STRUCTURAL DYNAMICS Real-life structures are subjected to loads which vary with time Except self weight of the structure, all other loads vary with time In many cases, this variation of the load
More informationLecture 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 informationav 1 x 2 + 4y 2 + xy + 4z 2 = 16.
74 85 Eigenanalysis The subject of eigenanalysis seeks to find a coordinate system, in which the solution to an applied problem has a simple expression Therefore, eigenanalysis might be called the method
More informationMath 3313: Differential Equations Second-order ordinary differential equations
Math 3313: Differential Equations Second-order ordinary differential equations Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Mass-spring & Newton s 2nd law Properties
More informationChapter Two Elements of Linear Algebra
Chapter Two Elements of Linear Algebra Previously, in chapter one, we have considered single first order differential equations involving a single unknown function. In the next chapter we will begin to
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More informationThis property turns out to be a general property of eigenvectors of a symmetric A that correspond to distinct eigenvalues as we shall see later.
34 To obtain an eigenvector x 2 0 2 for l 2 = 0, define: B 2 A - l 2 I 2 = È 1, 1, 1 Î 1-0 È 1, 0, 0 Î 1 = È 1, 1, 1 Î 1. To transform B 2 into an upper triangular matrix, subtract the first row of B 2
More information