MAC Module 12 Eigenvalues and Eigenvectors
|
|
- Emory Webster
- 5 years ago
- Views:
Transcription
1 MAC 23 Module 2 Eigenvalues and Eigenvectos
2 Leaning Objectives Upon completing this module, you should be able to:. Solve the eigenvalue poblem by finding the eigenvalues and the coesponding eigenvectos of an n x n matix. Find the algebaic multiplicity and the geometic multiplicity of an eigenvalue. 2. Find a basis fo each eigenspace of an eigenvalue. 3. Detemine whethe a matix A is diagonalizable. 4. Find a matix P, P -, and D that diagonalize A if A is diagonalizable. 5. Find an othogonal matix P with P - = P T and D that diagonalize A if A is symmetic and diagonalizable. 6. Detemine the powe and the eigenvalues of a matix, A k. Rev.F9 2
3 Eigenvalues and Eigenvectos The majo topics in this module: Eigenvalues, Eigenvectos, Eigenspace, Diagonalization and Othogonal Diagonalization Rev.9 3
4 What ae Eigenvalues and Eigenvectos? If A is an n x n matix and λ is a scala fo which Ax = λx has a nontivial solution x Rⁿ, then λ is an eigenvalue of A and x is a coesponding eigenvecto of A. Ax = λx is called the eigenvalue poblem fo A. Note that we can ewite the equation Ax = λx = λi n x as follows: λi n x - Ax = o (λi n - A)x =. x = is the tivial solution. But ou solutions must be nonzeo vectos called eigenvectos that coespond to each of the distinct eigenvalues. Rev.F9 4
5 What ae Eigenvalues and Eigenvectos? Since we seek a nontivial solution to (λi n - A)x = (λi - A)x =, λi - A must be singula to have solutions x. This means that the det(λi - A) =. The det(λi - A) = p(λ) = is the chaacteistic equation, whee det(λi - A) = p(λ) is the chaacteistic polynomial. The deg(p(λ)) = n and the n oots of p(λ), λ, λ 2,,λ n, ae the eigenvalues of A. The polynomial p(λ) always has n oots, so the zeos always exist; but some may be complex and some may be epeated. In ou examples, all of the oots will be eal. Fo each λ i we solve fo x i = p i the coesponding eigenvecto, and Ap i = λ i p i fo each distinct eigenvalue. Rev.F9 5
6 Rev.F9 How to Solve the Eigenvalue Poblem, Ax = λx? Example : Find the eigenvalues and the coesponding eigenvectos of A. 3 2 A 5 Step : Find the chaacteistic equation of A and solve fo its eigenvalues. p() I A ( 3) (2)(5) 2 3 ( 5) ( 2) Thus, the eigenvalues ae 5, 2 2. Each eigenvalue has algebaic multiplicity. Let D which is a diagonal matix. 6
7 How to Solve the Eigenvalue Poblem, Ax = λx?(cont.) Step 2: Use Gaussian elimination with back-substitution to solve the (λi - A) x = fo λ and λ 2. The augmented matix fo the system with λ = -5: 5I A ~ 2 2 The second column is not a leading column, so x 2 = t is a fee vaiable, and x = x 2 = t. Thus, the solution coesponding to λ = -5 is Rev.F9 x x x t t : t 5 2 2,t. 7.
8 How to Solve the Eigenvalue Poblem, Ax = λx? (Cont.) Since t is a fee vaiable, thee ae infinitely many eigenvectos. Fo convenience, we choose t =, and x p as the eigenvecto fo λ = -5. If we want p to be a unit vecto, we will choose t so that / 2 Howeve, t = is fine in this p poblem. / 2. The augmented matix fo the system with λ 2 = 2: 2I A : Rev.F :
9 How to Solve the Eigenvalue Poblem, Ax = λx? (Cont.) Again, the second column is not a leading column, so x 2 = t is a fee vaiable, and x = - 2x 2 / 5 = - 2t / 5. Thus, the solution coesponding to λ = 2 is x x x 2 2 t 5 t t 2 5,t. Fo convenience, we choose t = 5 and x p as the eigenvecto fo λ = 2. Alight, we have finished solving the eigenvalue poblem fo 3 2 A. 5 Rev.F9 9
10 Eigenspaces Fo each eigenvalue λ, thee is an eigenspace E λ with a basis fomed fom the linealy independent eigenvectos fo λ. The dim(e λ ) is the geometic multiplicity of λ, which is the numbe of linealy independent eigenvectos associated with λ. We will see that the geometic multiplicity equals the algebaic multiplicity fo each eigenvalue. B { p } is a basis fo E, the eigenspace of, and B 2 { p 2 } is a basis fo E 2, the eigenspace of 2. So, E span(b ) span({ p }) span E 2 span(b 2 ) span({ p 2 }) span 2 5 and dim(e ) dim(e 2 ). T T, Rev.F9
11 How to Detemine if a Matix A is Diagonalizable? If P = [ p p 2 ], then AP = PD even if A is not diagonalizable. Since AP p 2 p () (2) 5() 2(5) p p2 5 2 A p p2 A p Ap PD. Rev.F9
12 How to Detemine if a Matix A is Diagonalizable? So, ou two distinct eigenvalues both have algebaic multiplicity and geometic multiplicity. This ensues that p and p 2 ae not scala multiples of each othe; thus, p and p 2 ae linealy independent eigenvectos of A. Since A is 2 x 2 and thee ae two linealy independent eigenvectos fom the solution of the eigenvalue poblem, A is diagonalizable and P - AP = D. We can now constuct P, P - and D. Let P Then, P Rev.F9 p p2 5 / 7 2 / 7 / 7 / and D
13 How to Detemine if a Matix A is Diagonalizable? (Cont.) Note that, if we multiply both sides on the left by P, then AP A Rev.F9 5 P AP p p / 7 2 / 7 / 7 / 7 A p Ap / 7 2 / 7 / 7 / p p2 2 5 p p / 7 / 7 / 7 4 / () 2(2) 5() 2(5) PD becomes D.
14 How to Detemine if a Matix A is Diagonalizable? (Cont.) Example 2: Find the eigenvalues and eigenvectos fo A. 3 4 A 3 Step : Find the eigenvalues fo A. Recall: The deteminant of a tiangula matix is the poduct of the elements at the diagonal. Thus, the chaacteistic equation of A is λ = has algebaic multiplicity and λ 2 = 3 has algebaic multiplicity 2.. Rev.F9 p() det(i A) I A ( 3) 2 ( ). 4
15 How to Detemine if a Matix A is Diagonalizable? (Cont.) Step 2: Use Gaussian elimination with back-substitution to solve (λi - A) x =. Fo λ =, the augmented matix fo the system is I A : :. Column 3 is not a leading column and x 3 = t is a fee vaiable. The geometic multiplicity of λ = is one, since thee is only one fee vaiable. x 2 = and x = 2x 2 = Rev.F9 5
16 How to Detemine if a Matix A is Diagonalizable? (Cont.) The eigenvecto coesponding to λ = is x x x 2 x 3 t t ou choice fo the eigenvecto.. If we choose t, then p B { p } is a basis fo the eigenspace, E, with dim(e ). The dimension of the eigenspace is because the eigenvalue has only one linealy independent eigenvecto. Thus, the geometic multiplicity is and the algebaic multiplicity is fo λ =. is Rev.F9 6
17 How to Detemine if a Matix A is Diagonalizable? (Cont.) The augmented matix fo the system with λ 2 = 3 is 3I A 4 2 : : : Column is not a leading column and x = t is a fee vaiable. Since thee is only one fee vaiable, the geometic multiplicity of λ 2 is one. Rev.F9 7
18 How to Detemine if a Matix A is Diagonalizable? (Cont.) x 2 = x 3 = and the eigenvecto coesponding to λ 2 = 3 is x t x x 2 t,we choose t, and p 2 is x 3 ou choice fo the eigenvecto. B 2 { p 2 } is a basis fo the eigenspace, E 2, with dim(e 2 ). The dimension of the eigenspace is because the eigenvalue has only one linealy independent eigenvecto. Thus, the geometic multiplicity is while the algebaic multiplicity is 2 fo λ 2 = 3. This means thee will not be enough linealy independent eigenvectos fo A to be diagonalizable. Thus, A is not diagonalizable wheneve the geometic multiplicity is less than the algebaic multiplicity fo any eigenvalue. Rev.F9 8
19 How to Detemine if a Matix A is Diagonalizable? (Cont.) This time, AP A p 2 p2 2 p2 p p2 p2 p p2 p2 A p Ap 2 Ap PD. P does not exist since the columns of P ae not linealy independent. It is not possible to solve fo D P AP, so A is not diagonalizable. Rev.F9 9
20 A is Diagonalizable iff A is Simila to a Diagonal Matix Fo A, an n x n matix, with chaacteistic polynomial oots, 2,..., n, then AP A p p2... p n A p Ap 2... Ap n p 2 p2... n pn PD p p2... p n O n fo eigenvalues λ i of A with coesponding eigenvectos p i. P is invetible iff the eigenvectos that fom its columns ae linealy independent iff dim(e i ) geometic multiplicity algbaic multiplicity fo each distinct i. Rev.F9 2
21 A is Diagonalizable iff A is Simila to a Diagonal Matix (Cont.) This gives us n linealy independent eigenvectos fo P, so P - exists. Theefoe, A is diagonalizable since P AP P PD D O n. The squae matices S and T ae simila iff thee exists a nonsingula P such that S = P - TP o PSP - = T. Since A is simila to a diagonal matix, A is diagonalizable. Rev.F9 2
22 Anothe Example Example 4: Solve the eigenvalue poblem Ax = λx and find the eigenspace, algebaic multiplicity, and geometic multiplicity fo each eigenvalue. A Step : Wite down the chaacteistic equation of A and solve fo its eigenvalues. p() I A ()( ) Rev.F9 22
23 p() Anothe Example (Cont.) ( )( 4)( 5) 8( ) ( 2 5 4)( 5) 8( ) ( ) ( )( 2 2) ( )( 2)( ) ( 2)( ) 2. So the eigenvalues ae 2, 2. Since the facto (λ - 2) is fist powe, λ = 2 is not a epeated oot. λ = 2 has an algebaic multiplicity of. On the othe hand, the facto (λ +) is squaed, λ 2 = - is a epeated oot, and it has an algebaic multiplicity of 2. Rev.F9 23
24 Anothe Example (Cont.) Step 2: Use Gaussian elimination with back-substitution to solve (λi - A) x = fo λ and λ 2. Fo λ = 2, the augmented matix fo the system is 2I A ~ / ~ : / 2 3 / 2 / 2. In this case, x 3 =, x 2 =, and x = -/2() + = + =. Rev.F9 24
25 25 Rev.F9 Anothe Example (Cont.) Thus, the eigenvecto coesponding to λ = 2 is x x x 2 x 3,. If we choose p, then B is a basis fo the eigenspace of 2. E span({ p }) and dim(e ), so the geometic multiplicity is. A x 2 x o (2I A) x
26 Anothe Example (Cont.) Fo λ 2 = -, the augmented matix fo the system is ()I A ~ ~ x 3 = t, x 2 = s, and x = -s + 2t. Thus, the solution has two linealy independent eigenvectos fo λ 2 = - with x x x 2 x 3 s 2t s t s t 2,s,t. Rev.F9 26
27 27 Rev.F9 Anothe Example (Cont.) If we choose p 2, and p 3 2, then B 2, 2 is a basis fo E 2 span({ p 2, p 3 }) and dim(e 2 ) 2, so the geometic multiplicity is 2.
28 Anothe Example (Cont.) Since the geometic multiplicity is equal to the algebaic multiplicity fo each distinct eigenvalue, we found thee linealy independent eigenvectos. The matix A is diagonalizable since P = [p p 2 p 3 ] is nonsingula. Thus, we have AP PD as follows : Rev.F9 28
29 We can find P - as follows: Anothe Example (Cont.) P I : : : 2. So, P 2. Rev.F9 29
30 Anothe Example (Cont.) Note that A and D ae simila matices. AP PD gives us A APP PDP. Thus, PDP A Rev.F9 3
31 Anothe Example (Cont.) Also, D = P - AP = D So, A and D ae simila with D = P - AP and A = PD P -. Rev.F9 3
32 The Othogonal Diagonalization of a Symmetic Matix If P is an othogonal matix, its invese is its tanspose, P - = P T. Since P T P T p T p 2 p p2 p3 p T 3 p p p p 2 p p 3 p 2 p p2 p 2 p2 p 3 p 3 p p3 p 2 p3 p 3 p i p j T T T p p p p2 p p3 p T T T 2 p p2 p2 p2 p3 T T p 3 p p3 p2 p T 3p 3 I because p i p j fo i j and p i p j fo i j,2,3. So, P P T. Rev.F9 32
33 The Othogonal Diagonalization of a Symmetic Matix (Cont.) A is a symmetic matix if A = A T. Let A be diagonalizable so that A = PDP -. But A = A T and A T (PDP ) T (PDP T ) T (P T ) T D T P T PDP T A. This shows that fo a symmetic matix A to be diagonalizable, P must be othogonal. If P - P T, then A A T. The eigenvectos of A ae mutually othogonal but not othonomal. This means that the eigenvectos must be scaled to unit vectos so that P is othogonal and composed of othonomal columns. Rev.F9 33
34 The Othogonal Diagonalization of a Symmetic Matix (Cont.) Example 5: Detemine if the symmetic matix A is diagonalizable; if it is, then find the othogonal matix P that othogonally diagonalizes the symmetic matix A. Let A () 3 ( 2) 5 5 2, then det(i A) ( 2)( 5)2 ( 2) ( 4)( 2 4 2) ( 4)( 2)( 6) Thus, 4, 2 2, 3 6. Since we have thee distinct eigenvalues, we will see that we ae guaanteed to have thee linealy independent eigenvectos. Rev.F9 34
35 The Othogonal Diagonalization of a Symmetic Matix (Cont.) Since λ = 4, λ 2 = -2, and λ 3 = 6, ae distinct eigenvalues, each of the eigenvalues has algebaic multiplicity. An eigenvalue must have geometic multiplicity of at least one. Othewise, we will have the tivial solution. Thus, we have thee linealy independent eigenvectos. We will use Gaussian elimination with back-substitution as follows: Rev.F9 35
36 The Othogonal Diagonalization of a Symmetic Matix (Cont.) Rev.F9 Fo 4, I A 6 : : x 2 s, x 3, x s. x x x 2 x 3 s o p / 2 / 2. 36
37 The Othogonal Diagonalization of a Symmetic Matix (Cont.) Rev.F9 Fo 2 2, 2 I A 7 7 : 7 x 3 s, x 2, x. x x x 2 x 3 s o p 2. 37
38 The Othogonal Diagonalization of a Symmetic Matix (Cont.) Fo 3 6, 3 I A 8 : x 2 s, x 3, x s. x x x 2 x 3 s o p 3 / 2 / 2. Rev.F9 38
39 The Othogonal Diagonalization of a Symmetic Matix (Cont.) As we can see the eigenvectos of A ae distinct, so {p, p 2, p 3 } is linealy independent, P - exists fo P =[p p 2 p 3 ] and AP PD PDP. Thus A is diagonalizable. Since A = A T (A is a symmetic matix) and P is othogonal with appoximate scaling of p, p 2, p 3, P - = P T. PP PP T / 2 / 2 / 2 / 2 / 2 / 2 / 2 / 2 I. Rev.F9 39
40 The Othogonal Diagonalization of a Symmetic Matix (Cont.) As we can see the eigenvectos of A ae distinct, so {p, p 2, p 3 } is linealy independent, P - exists fo P =[p p 2 p 3 ] and AP PD PDP. Thus A is diagonalizable. Since A = A T (A is a symmetic matix) and P is othogonal with appoximate scaling of p, p 2, p 3, P - = P T. PP PP T / 2 / 2 / 2 / 2 / 2 / 2 / 2 / 2 I. Rev.F9 4
41 The Othogonal Diagonalization of a Symmetic Matix (Cont.) Note that A and D ae simila matices. PD P - = PDP T / 2 / 2 / 2 / / 2 / 2 / 2 / 2 4 / 2 6 / 2 4 / 2 6 / 2 2 / 2 / 2 / 2 / A. Rev.F9 4
42 The Othogonal Diagonalization of a Symmetic Matix (Cont.) Also, D = P - AP = P T AP / 2 / 2 / 2 / / 2 / 2 / 2 / 2 4 / 2 4 / / 2 6 / 2 / 2 / 2 / 2 / So, A and D ae simila with D = P T AP and A = PDP T. Rev.F9 42
43 The Othogonal Diagonalization of a Symmetic Matix (Cont.) A T = (PD P - ) T = (PD P T ) T = (P T ) T D T P T = P D T P T T / 2 / 2 4 / 2 / 2 / 2 / 2 2 / 2 / 2 6 / 2 / 2 / 2 / 2 Rev.F / 2 / 2 / 2 / 2 4 / 2 6 / 2 / 2 / / 2 6 / / 2 / 2 2 = A. This shows that if A is a symmetic matix, P must be othogonal with P - = P T. T 43
44 How to Detemine the Powe and the Eigenvalues of a Matix, A k? Fom Example, the diagonal matix fo matix A is : D and A 3 PDP PDP PDP PD 3 P 2 5 (5) 3 64 / / / 7 29 / 7 (2) 3 P AP D A PDP, 5 / 7 2 / 7 / 7 / / 7 2 / 7 / 7 / 7. Fo A 3,the eigenvalues, ae 3 25 aand In geneal, the powe of a matix, A k PD k P. and the eigenvalues ae i k, whee i is on the main diagonal of D. Rev.F9 44
45 We have leaned to: What have we leaned?. Solve the eigenvalue poblem by finding the eigenvalues and the coesponding eigenvectos of an n x n matix. Find the algebaic multiplicity and the geometic multiplicity of an eigenvalue. 2. Find a basis fo each eigenspace of an eigenvalue. 3. Detemine whethe a matix A is diagonalizable. 4. Find a matix P, P -, and D that diagonalize A if A is diagonalizable. 5. Find an othogonal matix P with P - = P T and D that diagonalize A if A is symmetic and diagonalizable. 6. Detemine the powe and the eigenvalues of a matix, A k. Rev.F9 45
46 Cedit Some of these slides have been adapted/modified in pat/whole fom the following textbook: Anton, Howad: Elementay Linea Algeba with Applications, 9th Edition Rev.F9 46
MAC Module 12 Eigenvalues and Eigenvectors
MAC 23 Module 2 Eigenvalues and Eigenvectors Learning Objectives Upon completing this module, you should be able to:. Solve the eigenvalue problem by finding the eigenvalues and the corresponding eigenvectors
More informationMAC Module 12 Eigenvalues and Eigenvectors. Learning Objectives. Upon completing this module, you should be able to:
MAC Module Eigenvalues and Eigenvectors Learning Objectives Upon completing this module, you should be able to: Solve the eigenvalue problem by finding the eigenvalues and the corresponding eigenvectors
More informationA Crash Course in (2 2) Matrices
A Cash Couse in ( ) Matices Seveal weeks woth of matix algeba in an hou (Relax, we will only stuy the simplest case, that of matices) Review topics: What is a matix (pl matices)? A matix is a ectangula
More informationPHYS 705: Classical Mechanics. Small Oscillations
PHYS 705: Classical Mechanics Small Oscillations Fomulation of the Poblem Assumptions: V q - A consevative system with depending on position only - The tansfomation equation defining does not dep on time
More informationComplex Eigenvalues. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB
Pepaed by Vince Zaccone Fo ampus Leaning ssistance Sevices at USB omplex Numbes When solving fo the oots of a quadatic equation, eal solutions can not be found when the disciminant is negative. In these
More informationMATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form
MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE ANDRAS VASY We conside second ode constant coefficient scala linea PDEs on R n. These have the fom Lu = f L = a ij xi xj + b i xi + c i whee a ij b i and
More information18.06 Problem Set 4 Solution
8.6 Poblem Set 4 Solution Total: points Section 3.5. Poblem 2: (Recommended) Find the lagest possible numbe of independent vectos among ) ) ) v = v 4 = v 5 = v 6 = v 2 = v 3 =. Solution (4 points): Since
More informationChapter 5 Linear Equations: Basic Theory and Practice
Chapte 5 inea Equations: Basic Theoy and actice In this chapte and the next, we ae inteested in the linea algebaic equation AX = b, (5-1) whee A is an m n matix, X is an n 1 vecto to be solved fo, and
More informationPart V: Closed-form solutions to Loop Closure Equations
Pat V: Closed-fom solutions to Loop Closue Equations This section will eview the closed-fom solutions techniques fo loop closue equations. The following thee cases will be consideed. ) Two unknown angles
More informationMODULE 5a and 5b (Stewart, Sections 12.2, 12.3) INTRO: In MATH 1114 vectors were written either as rows (a1, a2,..., an) or as columns a 1 a. ...
MODULE 5a and 5b (Stewat, Sections 2.2, 2.3) INTRO: In MATH 4 vectos wee witten eithe as ows (a, a2,..., an) o as columns a a 2... a n and the set of all such vectos of fixed length n was called the vecto
More informationMU+CU+KU=F MU+CU+KU=0
MEEN 67 Handout # MODAL ANALYSIS OF MDOF Systems with VISCOUS DAMPING ^ Symmetic he motion of a n-dof linea system is descibed by the set of 2 nd ode diffeential equations MU+CU+KU=F t whee U (t) and F
More informationAppendix A. Appendices. A.1 ɛ ijk and cross products. Vector Operations: δ ij and ɛ ijk
Appendix A Appendices A1 ɛ and coss poducts A11 Vecto Opeations: δ ij and ɛ These ae some notes on the use of the antisymmetic symbol ɛ fo expessing coss poducts This is an extemely poweful tool fo manipulating
More information3.6 Applied Optimization
.6 Applied Optimization Section.6 Notes Page In this section we will be looking at wod poblems whee it asks us to maimize o minimize something. Fo all the poblems in this section you will be taking the
More informationOn the Quasi-inverse of a Non-square Matrix: An Infinite Solution
Applied Mathematical Sciences, Vol 11, 2017, no 27, 1337-1351 HIKARI Ltd, wwwm-hikaicom https://doiog/1012988/ams20177273 On the Quasi-invese of a Non-squae Matix: An Infinite Solution Ruben D Codeo J
More informationof the contestants play as Falco, and 1 6
JHMT 05 Algeba Test Solutions 4 Febuay 05. In a Supe Smash Bothes tounament, of the contestants play as Fox, 3 of the contestants play as Falco, and 6 of the contestants play as Peach. Given that thee
More informationConjugate Gradient Methods. Michael Bader. Summer term 2012
Gadient Methods Outlines Pat I: Quadatic Foms and Steepest Descent Pat II: Gadients Pat III: Summe tem 2012 Pat I: Quadatic Foms and Steepest Descent Outlines Pat I: Quadatic Foms and Steepest Descent
More informationPolitical Science 552
Political Science 55 Facto and Pincial Comonents Path : Wight s Rules 4 v 4 4 4u R u R v 4. Path may ass though any vaiable only once on a single tavese. Path may go backwads, but not afte going fowad.
More informationLecture 15, 16: Diagonalization
Lecture 15, 16: Diagonalization Motivation: Eigenvalues and Eigenvectors are easy to compute for diagonal matrices. Hence, we would like (if possible) to convert matrix A into a diagonal matrix. Suppose
More information8 Separation of Variables in Other Coordinate Systems
8 Sepaation of Vaiables in Othe Coodinate Systems Fo the method of sepaation of vaiables to succeed you need to be able to expess the poblem at hand in a coodinate system in which the physical boundaies
More informationEM Boundary Value Problems
EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do
More informationPROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.
POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and
More information-Δ u = λ u. u(x,y) = u 1. (x) u 2. (y) u(r,θ) = R(r) Θ(θ) Δu = 2 u + 2 u. r = x 2 + y 2. tan(θ) = y/x. r cos(θ) = cos(θ) r.
The Laplace opeato in pola coodinates We now conside the Laplace opeato with Diichlet bounday conditions on a cicula egion Ω {(x,y) x + y A }. Ou goal is to compute eigenvalues and eigenfunctions of the
More informationRigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018
Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining
More informationRecall : Eigenvalues and Eigenvectors
Recall : Eigenvalues and Eigenvectors Let A be an n n matrix. If a nonzero vector x in R n satisfies Ax λx for a scalar λ, then : The scalar λ is called an eigenvalue of A. The vector x is called an eigenvector
More informationTHE CONE THEOREM JOEL A. TROPP. Abstract. We prove a fixed point theorem for functions which are positive with respect to a cone in a Banach space.
THE ONE THEOEM JOEL A. TOPP Abstact. We pove a fixed point theoem fo functions which ae positive with espect to a cone in a Banach space. 1. Definitions Definition 1. Let X be a eal Banach space. A subset
More informationDymore User s Manual Two- and three dimensional dynamic inflow models
Dymoe Use s Manual Two- and thee dimensional dynamic inflow models Contents 1 Two-dimensional finite-state genealized dynamic wake theoy 1 Thee-dimensional finite-state genealized dynamic wake theoy 1
More informationA Power Method for Computing Square Roots of Complex Matrices
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 13, 39345 1997 ARTICLE NO. AY975517 A Powe Method fo Computing Squae Roots of Complex Matices Mohammed A. Hasan Depatment of Electical Engineeing, Coloado
More informationWeighted least-squares estimators of parametric functions of the regression coefficients under a general linear model
Ann Inst Stat Math (2010) 62:929 941 DOI 10.1007/s10463-008-0199-8 Weighted least-squaes estimatos of paametic functions of the egession coefficients unde a geneal linea model Yongge Tian Received: 9 Januay
More informationA pathway to matrix-variate gamma and normal densities
Linea Algeba and its Applications 396 005 317 38 www.elsevie.com/locate/laa A pathway to matix-vaiate gamma and nomal densities A.M. Mathai Depatment of Mathematics and Statistics, McGill Univesity, 805
More information(n 1)n(n + 1)(n + 2) + 1 = (n 1)(n + 2)n(n + 1) + 1 = ( (n 2 + n 1) 1 )( (n 2 + n 1) + 1 ) + 1 = (n 2 + n 1) 2.
Paabola Volume 5, Issue (017) Solutions 151 1540 Q151 Take any fou consecutive whole numbes, multiply them togethe and add 1. Make a conjectue and pove it! The esulting numbe can, fo instance, be expessed
More informationI. CONSTRUCTION OF THE GREEN S FUNCTION
I. CONSTRUCTION OF THE GREEN S FUNCTION The Helmohltz equation in 4 dimensions is 4 + k G 4 x, x = δ 4 x x. In this equation, G is the Geen s function and 4 efes to the dimensionality. In the vey end,
More informationScientific Computing II
Scientific Computing II Conjugate Gadient Methods Michael Bade Summe 2014 Conjugate Gadient Methods, Summe 2014 1 Families of Iteative Solves elaxation methods: Jacobi-, Gauss-Seidel-Relaxation,... Ove-Relaxation-Methods
More information9.1 The multiplicative group of a finite field. Theorem 9.1. The multiplicative group F of a finite field is cyclic.
Chapte 9 Pimitive Roots 9.1 The multiplicative goup of a finite fld Theoem 9.1. The multiplicative goup F of a finite fld is cyclic. Remak: In paticula, if p is a pime then (Z/p) is cyclic. In fact, this
More information5.61 Physical Chemistry Lecture #23 page 1 MANY ELECTRON ATOMS
5.6 Physical Chemisty Lectue #3 page MAY ELECTRO ATOMS At this point, we see that quantum mechanics allows us to undestand the helium atom, at least qualitatively. What about atoms with moe than two electons,
More informationDIAGONALIZATION. In order to see the implications of this definition, let us consider the following example Example 1. Consider the matrix
DIAGONALIZATION Definition We say that a matrix A of size n n is diagonalizable if there is a basis of R n consisting of eigenvectors of A ie if there are n linearly independent vectors v v n such that
More informationQuantum Mechanics I - Session 5
Quantum Mechanics I - Session 5 Apil 7, 015 1 Commuting opeatos - an example Remine: You saw in class that Â, ˆB ae commuting opeatos iff they have a complete set of commuting obsevables. In aition you
More informationAn Application of Fuzzy Linear System of Equations in Economic Sciences
Austalian Jounal of Basic and Applied Sciences, 5(7): 7-14, 2011 ISSN 1991-8178 An Application of Fuzzy Linea System of Equations in Economic Sciences 1 S.H. Nassei, 2 M. Abdi and 3 B. Khabii 1 Depatment
More informationMany Electron Atoms. Electrons can be put into approximate orbitals and the properties of the many electron systems can be catalogued
Many Electon Atoms The many body poblem cannot be solved analytically. We content ouselves with developing appoximate methods that can yield quite accuate esults (but usually equie a compute). The electons
More informationPreliminary Exam: Quantum Physics 1/14/2011, 9:00-3:00
Peliminay Exam: Quantum Physics /4/ 9:-: Answe a total of SIX questions of which at least TWO ae fom section A and at least THREE ae fom section B Fo you answes you can use eithe the blue books o individual
More informationNew problems in universal algebraic geometry illustrated by boolean equations
New poblems in univesal algebaic geomety illustated by boolean equations axiv:1611.00152v2 [math.ra] 25 Nov 2016 Atem N. Shevlyakov Novembe 28, 2016 Abstact We discuss new poblems in univesal algebaic
More informationΔt The textbook chooses to say that the average velocity is
1-D Motion Basic I Definitions: One dimensional motion (staight line) is a special case of motion whee all but one vecto component is zeo We will aange ou coodinate axis so that the x-axis lies along the
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More information3D-Central Force Problems II
5.73 ectue # - 1 3D-Cental Foce Poblems II ast time: [x,p] vecto commutation ules: genealize fom 1-D to 3-D conugate position and momentum components in Catesian coodinates Coespondence Pinciple Recipe
More informationx 1 b 1 Consider the midpoint x 0 = 1 2
1 chapte 2 : oot-finding def : Given a function f(), a oot is a numbe satisfying f() = 0. e : f() = 2 3 = ± 3 question : How can we find the oots of a geneal function f()? 2.1 bisection method idea : Find
More information6 Matrix Concentration Bounds
6 Matix Concentation Bounds Concentation bounds ae inequalities that bound pobabilities of deviations by a andom vaiable fom some value, often its mean. Infomally, they show the pobability that a andom
More informationJORDAN CANONICAL FORM AND ITS APPLICATIONS
JORDAN CANONICAL FORM AND ITS APPLICATIONS Shivani Gupta 1, Kaajot Kau 2 1,2 Matheatics Depatent, Khalsa College Fo Woen, Ludhiana (India) ABSTRACT This pape gives a basic notion to the Jodan canonical
More informationPractice Integration Math 120 Calculus I Fall 2015
Pactice Integation Math 0 Calculus I Fall 05 Hee s a list of pactice eecises. Thee s a hint fo each one as well as an answe with intemediate steps... ( + d. Hint. Answe. ( 8 t + t + This fist set of indefinite
More informationQuestion Bank. Section A. is skew-hermitian matrix. is diagonalizable. (, ) , Evaluate (, ) 12 about = 1 and = Find, if
Subject: Mathematics-I Question Bank Section A T T. Find the value of fo which the matix A = T T has ank one. T T i. Is the matix A = i is skew-hemitian matix. i. alculate the invese of the matix = 5 7
More informationPractice Integration Math 120 Calculus I D Joyce, Fall 2013
Pactice Integation Math 0 Calculus I D Joyce, Fall 0 This fist set of indefinite integals, that is, antideivatives, only depends on a few pinciples of integation, the fist being that integation is invese
More information= e2. = 2e2. = 3e2. V = Ze2. where Z is the atomic numnber. Thus, we take as the Hamiltonian for a hydrogenic. H = p2 r. (19.4)
Chapte 9 Hydogen Atom I What is H int? That depends on the physical system and the accuacy with which it is descibed. A natual stating point is the fom H int = p + V, (9.) µ which descibes a two-paticle
More informationDiagonalization of Matrix
of Matrix King Saud University August 29, 2018 of Matrix Table of contents 1 2 of Matrix Definition If A M n (R) and λ R. We say that λ is an eigenvalue of the matrix A if there is X R n \ {0} such that
More informationNon-Linear Dynamics Homework Solutions Week 2
Non-Linea Dynamics Homewok Solutions Week Chis Small Mach, 7 Please email me at smach9@evegeen.edu with any questions o concens eguading these solutions. Fo the ececises fom section., we sketch all qualitatively
More informationWhen two numbers are written as the product of their prime factors, they are in factored form.
10 1 Study Guide Pages 420 425 Factos Because 3 4 12, we say that 3 and 4 ae factos of 12. In othe wods, factos ae the numbes you multiply to get a poduct. Since 2 6 12, 2 and 6 ae also factos of 12. The
More informationResearch Design - - Topic 17 Multiple Regression & Multiple Correlation: Two Predictors 2009 R.C. Gardner, Ph.D.
Reseach Design - - Topic 7 Multiple Regession & Multiple Coelation: Two Pedictos 009 R.C. Gadne, Ph.D. Geneal Rationale and Basic Aithmetic fo two pedictos Patial and semipatial coelation Regession coefficients
More informationPHYS 301 HOMEWORK #10 (Optional HW)
PHYS 301 HOMEWORK #10 (Optional HW) 1. Conside the Legende diffeential equation : 1 - x 2 y'' - 2xy' + m m + 1 y = 0 Make the substitution x = cos q and show the Legende equation tansfoms into d 2 y 2
More informationInformation Retrieval Advanced IR models. Luca Bondi
Advanced IR models Luca Bondi Advanced IR models 2 (LSI) Pobabilistic Latent Semantic Analysis (plsa) Vecto Space Model 3 Stating point: Vecto Space Model Documents and queies epesented as vectos in the
More informationKEPLER S LAWS OF PLANETARY MOTION
EPER S AWS OF PANETARY MOTION 1. Intoduction We ae now in a position to apply what we have leaned about the coss poduct and vecto valued functions to deive eple s aws of planetay motion. These laws wee
More informationSchur product of matrices and numerical radius (range) preserving maps. Dedicated to Professor Roger Horn on the occasion of his sixty fifth birthday.
Schu poduct of matices and numeical adius (ange) peseving maps Chi-Kwong Li 1 and Edwad Poon 2 Dedicated to Pofesso Roge Hon on the occasion of his sixty fifth bithday. Abstact Let F (A) be the numeical
More informationf h = u, h g = v, we have u + v = f g. So, we wish
Answes to Homewok 4, Math 4111 (1) Pove that the following examples fom class ae indeed metic spaces. You only need to veify the tiangle inequality. (a) Let C be the set of continuous functions fom [0,
More informationChapter 3: Theory of Modular Arithmetic 38
Chapte 3: Theoy of Modula Aithmetic 38 Section D Chinese Remainde Theoem By the end of this section you will be able to pove the Chinese Remainde Theoem apply this theoem to solve simultaneous linea conguences
More information10/04/18. P [P(x)] 1 negl(n).
Mastemath, Sping 208 Into to Lattice lgs & Cypto Lectue 0 0/04/8 Lectues: D. Dadush, L. Ducas Scibe: K. de Boe Intoduction In this lectue, we will teat two main pats. Duing the fist pat we continue the
More informationEigenvalues and Eigenvectors 7.2 Diagonalization
Eigenvalues and Eigenvectors 7.2 Diagonalization November 8 Goals Suppose A is square matrix of order n. Provide necessary and sufficient condition when there is an invertible matrix P such that P 1 AP
More informationChapter 2: Introduction to Implicit Equations
Habeman MTH 11 Section V: Paametic and Implicit Equations Chapte : Intoduction to Implicit Equations When we descibe cuves on the coodinate plane with algebaic equations, we can define the elationship
More informationMethod for Approximating Irrational Numbers
Method fo Appoximating Iational Numbes Eic Reichwein Depatment of Physics Univesity of Califonia, Santa Cuz June 6, 0 Abstact I will put foth an algoithm fo poducing inceasingly accuate ational appoximations
More informationA Relativistic Electron in a Coulomb Potential
A Relativistic Electon in a Coulomb Potential Alfed Whitehead Physics 518, Fall 009 The Poblem Solve the Diac Equation fo an electon in a Coulomb potential. Identify the conseved quantum numbes. Specify
More informationTopic 4a Introduction to Root Finding & Bracketing Methods
/8/18 Couse Instucto D. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: cumpf@utep.edu Topic 4a Intoduction to Root Finding & Backeting Methods EE 4386/531 Computational Methods in EE Outline
More informationSpring 2019 Exam 2 3/27/19 Time Limit: / Problem Points Score. Total: 280
Math 307 Spring 2019 Exam 2 3/27/19 Time Limit: / Name (Print): Problem Points Score 1 15 2 20 3 35 4 30 5 10 6 20 7 20 8 20 9 20 10 20 11 10 12 10 13 10 14 10 15 10 16 10 17 10 Total: 280 Math 307 Exam
More informationReading Assignment. Problem Description for Homework #9. Read Chapters 29 and 30.
Reading Assignment Read Chaptes 29 and 30. Poblem Desciption fo Homewok #9 In this homewok, you will solve the inhomogeneous Laplace s equation to calculate the electic scala potential that exists between
More informationLinear Algebra Math 221
Linea Algeba Math Open Book Eam Open Notes Sept Calculatos Pemitted Sho all ok (ecept #). ( pts) Gien the sstem of equations a) ( pts) Epess this sstem as an augmented mati. b) ( pts) Bing this mati to
More informationAs is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.
Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.
More informationLecture 16 Root Systems and Root Lattices
1.745 Intoduction to Lie Algebas Novembe 1, 010 Lectue 16 Root Systems and Root Lattices Pof. Victo Kac Scibe: Michael Cossley Recall that a oot system is a pai (V, ), whee V is a finite dimensional Euclidean
More informationLecture 3 Eigenvalues and Eigenvectors
Lecture 3 Eigenvalues and Eigenvectors Eivind Eriksen BI Norwegian School of Management Department of Economics September 10, 2010 Eivind Eriksen (BI Dept of Economics) Lecture 3 Eigenvalues and Eigenvectors
More informationSolutions to Problems : Chapter 19 Problems appeared on the end of chapter 19 of the Textbook
Solutions to Poblems Chapte 9 Poblems appeae on the en of chapte 9 of the Textbook 8. Pictue the Poblem Two point chages exet an electostatic foce on each othe. Stategy Solve Coulomb s law (equation 9-5)
More informationExceptional regular singular points of second-order ODEs. 1. Solving second-order ODEs
(May 14, 2011 Exceptional egula singula points of second-ode ODEs Paul Gaett gaett@math.umn.edu http://www.math.umn.edu/ gaett/ 1. Solving second-ode ODEs 2. Examples 3. Convegence Fobenius method fo solving
More informationEigenvalues and Eigenvectors 7.1 Eigenvalues and Eigenvecto
7.1 November 6 7.1 Eigenvalues and Eigenvecto Goals Suppose A is square matrix of order n. Eigenvalues of A will be defined. Eigenvectors of A, corresponding to each eigenvalue, will be defined. Eigenspaces
More informationMA557/MA578/CS557. Lecture 14. Prof. Tim Warburton. Spring
MA557/MA578/CS557 Lectue 4 Sping 3 Pof. Tim Wabuton timwa@math.unm.edu Matlab Notes To ceate a symbolic vaiable (say theta) use the command: theta = sym( theta ); Example manipulation: Matlab cont To evaluate
More informationConstruction and Analysis of Boolean Functions of 2t + 1 Variables with Maximum Algebraic Immunity
Constuction and Analysis of Boolean Functions of 2t + 1 Vaiables with Maximum Algebaic Immunity Na Li and Wen-Feng Qi Depatment of Applied Mathematics, Zhengzhou Infomation Engineeing Univesity, Zhengzhou,
More information. Using our polar coordinate conversions, we could write a
504 Chapte 8 Section 8.4.5 Dot Poduct Now that we can add, sutact, and scale vectos, you might e wondeing whethe we can multiply vectos. It tuns out thee ae two diffeent ways to multiply vectos, one which
More informationHua Xu 3 and Hiroaki Mukaidani 33. The University of Tsukuba, Otsuka. Hiroshima City University, 3-4-1, Ozuka-Higashi
he inea Quadatic Dynamic Game fo Discete-ime Descipto Systems Hua Xu 3 and Hioai Muaidani 33 3 Gaduate School of Systems Management he Univesity of suuba, 3-9- Otsua Bunyo-u, oyo -0, Japan xuhua@gssm.otsua.tsuuba.ac.jp
More information3.1 Random variables
3 Chapte III Random Vaiables 3 Random vaiables A sample space S may be difficult to descibe if the elements of S ae not numbes discuss how we can use a ule by which an element s of S may be associated
More informationarxiv: v1 [math.co] 1 Apr 2011
Weight enumeation of codes fom finite spaces Relinde Juius Octobe 23, 2018 axiv:1104.0172v1 [math.co] 1 Ap 2011 Abstact We study the genealized and extended weight enumeato of the - ay Simplex code and
More information556: MATHEMATICAL STATISTICS I
556: MATHEMATICAL STATISTICS I CHAPTER 5: STOCHASTIC CONVERGENCE The following efinitions ae state in tems of scala anom vaiables, but exten natually to vecto anom vaiables efine on the same obability
More information1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More informationA POWER IDENTITY FOR SECOND-ORDER RECURRENT SEQUENCES
A OWER IDENTITY FOR SECOND-ORDER RECURRENT SEQUENCES V. E. Hoggatt,., San ose State College, San ose, Calif. and D. A. Lind, Univesity of Viginia, Chalottesville, V a. 1. INTRODUCTION The following hold
More informationCMSC 425: Lecture 5 More on Geometry and Geometric Programming
CMSC 425: Lectue 5 Moe on Geomety and Geometic Pogamming Moe Geometic Pogamming: In this lectue we continue the discussion of basic geometic ogamming fom the eious lectue. We will discuss coodinate systems
More informationChromatic number and spectral radius
Linea Algeba and its Applications 426 2007) 810 814 www.elsevie.com/locate/laa Chomatic numbe and spectal adius Vladimi Nikifoov Depatment of Mathematical Sciences, Univesity of Memphis, Memphis, TN 38152,
More informationWhat to Expect on the Placement Exam
What to Epect on the Placement Eam Placement into: MTH o MTH 44 05 05 The ACCUPLACER placement eam is an adaptive test ceated by the College Boad Educational Testing Sevice. This document was ceated to
More informationQuantum Fourier Transform
Chapte 5 Quantum Fouie Tansfom Many poblems in physics and mathematics ae solved by tansfoming a poblem into some othe poblem with a known solution. Some notable examples ae Laplace tansfom, Legende tansfom,
More informationMath 124B February 02, 2012
Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial
More informationMultiple Experts with Binary Features
Multiple Expets with Binay Featues Ye Jin & Lingen Zhang Decembe 9, 2010 1 Intoduction Ou intuition fo the poect comes fom the pape Supevised Leaning fom Multiple Expets: Whom to tust when eveyone lies
More informationA Converse to Low-Rank Matrix Completion
A Convese to Low-Rank Matix Completion Daniel L. Pimentel-Alacón, Robet D. Nowak Univesity of Wisconsin-Madison Abstact In many pactical applications, one is given a subset Ω of the enties in a d N data
More informationVector d is a linear vector function of vector d when the following relationships hold:
Appendix 4 Dyadic Analysis DEFINITION ecto d is a linea vecto function of vecto d when the following elationships hold: d x = a xxd x + a xy d y + a xz d z d y = a yxd x + a yy d y + a yz d z d z = a zxd
More informationLacunary I-Convergent Sequences
KYUNGPOOK Math. J. 52(2012), 473-482 http://dx.doi.og/10.5666/kmj.2012.52.4.473 Lacunay I-Convegent Sequences Binod Chanda Tipathy Mathematical Sciences Division, Institute of Advanced Study in Science
More informationCALCULUS II Vectors. Paul Dawkins
CALCULUS II Vectos Paul Dawkins Table of Contents Peface... ii Vectos... 3 Intoduction... 3 Vectos The Basics... 4 Vecto Aithmetic... 8 Dot Poduct... 13 Coss Poduct... 21 2007 Paul Dawkins i http://tutoial.math.lama.edu/tems.aspx
More informationPermutations and Combinations
Pemutations and Combinations Mach 11, 2005 1 Two Counting Pinciples Addition Pinciple Let S 1, S 2,, S m be subsets of a finite set S If S S 1 S 2 S m, then S S 1 + S 2 + + S m Multiplication Pinciple
More informationConjugate Gradient (CG) Method
Optimization II Conugate Gaient CG Metho Anothe metho oes not equie explicit secon eivatives, an oes not even stoe appoximation to Hessian matix CG geneates sequence of conugate seach iections, implicitly
More informationx x2 2 B A ) v(0, t) = 0 and v(l, t) = 0. L 2. This is a familiar heat equation initial/boundary-value problem and has solution
Hints to homewok 7 8.2.d. The poblem is u t ku xx + k ux fx u t A u t B. It has a souce tem and inhomogeneous bounday conditions but none of them depend on t. So as in example 3 of the notes we should
More information2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum
2. Electostatics D. Rakhesh Singh Kshetimayum 1 2.1 Intoduction In this chapte, we will study how to find the electostatic fields fo vaious cases? fo symmetic known chage distibution fo un-symmetic known
More informationworking pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50
woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,
More informationQIP Course 10: Quantum Factorization Algorithm (Part 3)
QIP Couse 10: Quantum Factoization Algoithm (Pat 3 Ryutaoh Matsumoto Nagoya Univesity, Japan Send you comments to yutaoh.matsumoto@nagoya-u.jp Septembe 2018 @ Tokyo Tech. Matsumoto (Nagoya U. QIP Couse
More information