Gernot Hoffmann. Faddejev Algorithm for Eigenvalues and Eigenvectors
|
|
- Chester Singleton
- 6 years ago
- Views:
Transcription
1 Gerot Hoffma Faddejev Algorithm for Eigevalues ad Eigevectors Cotets. Itroductio. Faddejev algorithm / Geeral. Example / Diagoal matrix with distict eigevalues 4. Example / Matrix with complex eigevalues 5. Example / Matrix with equal eigevalues 6 4. Rotatio matrices / Geeral 7 4. Rotatio matrices / Eigevalues ad eigevectors 8 5. Newto-Bairstow algorithm 9 6. Refereces 0
2 . Itroductio May tass i physics, mechaics ad computer graphics result i the so-called eigevalue problem for matrices. Give is a quadratic matrix with x elemets, which are here real umbers. The matrix is for the geeral case ot symmetric. () Ax x ( ) ( A I) x 0 ( ) det( A I) 0 ( 4) ( A I) 0 Equatio () is a homogeeous system of liear equatios for uows. No-trivial solutios ca be expected if the determiat () is zero. This delivers primarily a polyomial fuctio f(). The eigevalues are foud as the roots of f()0, which is called the characteristic equatio. Equatio () results for each eigevalue i equatio (4) for the eigevector x. A eigevector ca be multiplied by ay factor. A uit eigevector has the Euclidia orm oe. Complex eigevectors are here ormalized for oe real value oe. The author is programmig a library for liear algebra ad related topics i PostScript, for simplicity called the Matrix Library [6], [7]. PostScript is geerally based o sigle float precisio - a rather limitig fact. Programmig umerically algorithms is hampered by lacig structures as i Fortra, Pascal or C. PostScript is used because of excellet graphics ad typographical quality for vector graphics. It is ot expected that highly complex tass of liear algebra have to be hadled. O this occasio, the author remembered the otherwise hardly ow algorithm by Leverrier, improved by Faddejev []. This is here simply called Faddejev algorithm. Just a few matrix operatios deliver the determiat, the iverse matrix ad the coefficiets of the characteristic equatio. Oce the roots are foud - here by Newto-Bairstow [] - the eigevectors are calculated by usig itermediate matrices from the previous process. Solvig the eigevalue problem via the characteristic equatio is owadays ot cosidered as professioal. Oe of the better alteratives is the R-Algorithm []. Programmig this i PostScript for o-symmetric matrices with possible complex solutios is for the author practically impossible. O the other had, for small matrices, without ay symmetry, for real or complex distict eigevalues ad eigevectors, the Faddejev method is straightforward. Practical problems arise for simple cases, e.g. the calculatio of the eigevalues for a x idetity matrix. The three eigevalues are all. A vector base of three liearly idepedet eigevectors ca be costructed, but it caot be foud by umerical calculatios. Faddejev s boo cotais a major bug, cocerig the eigevector extractio. The problem could be solved - slow but safe, if the tas is ot affected by geeral problems lie multiple eigevalues, cocerig the liearly idepedet vector base. This documet shows a couple of examples.
3 . Faddejev algorithm / Geeral Some simple matrix operatios deliver the ecessary iformatios for the calculatio of the determiat, the iverse matrix ad the coefficiets of the characteristic equatio. Oce the roots are foud (eigevalues), it should be possible to extract the respective eigevectors from a matrix, which is easily calculated, based o previous itermediate results. Ufortuately, the situatio is ot as simple as explaied by Faddejev. Faddejev algorithm for a matrix A: I is a idetity matrix. tr meas trace : tr( A) i a ii A A p tr( A ) / A AB p tr( A )/ B A p I A B A p I AB p tr( A )/ B A p I Always valid: B 0 Determiat ad iverse matrix: det A ( ) p If p 0 : A B / p Characteristic equatio: f( ) (- ) [ p p... p p ] 0 The eigevalues are the roots of the characteristic equatio. They ca be multiple. For istace 5 ad triple :,,,, 4 5 The matrix cotais eigevectors for i colums: I+ B B Faddejev assumes for differet eigevalues, that each colum of ca be used as the respective eigevector. This is ot always true. Horer recursio for the eigevector matrix ( i to -):,0 I i, i, + Bi,
4 . Faddejev algorithm / Diagoal matrix with distict eigevalues A diagoal matrix with diagoal values,, delivers the eigevalues,,. The roots of the characteristic equatio are foud i a upredictable order. The matrix cotais ideed for each eigevalue the respective eigevector: it is the oly oe o-zero colum. I a program oe caot chec for o-zero. Alteratively, the colum with the largest orm is selected. Faddejev says, that ay colum ca be used. This is correct, if the matrix does ot cosist of diagoal elemets or diagoally arraged blocs. Faddejev algorithm, arraged for maual calculatios ( ) : p tr( A )/ p tr( A )/ p tr( A )/ B A p I B A p I B A p A A A AB A AB I p 6 p p Determiat ad iverse matrix: det A ( ) p A B / p ( 6 / ) / / Characteristic equatio: p f( ) (- ) [ p... p p ] [ 6 + 6] 0 Eigevalues, assumed to be foud by this arbitrary order: Eigevector matrix ad ormalized eigevectors: I+ B B
5 . Faddejev algorithm / Matrix with complex eigevalues For complex eigevalues the eigevectors are complex as well. The orm for searchig the largest eigevector is based o the Euclidia orm of the absolute values of the complex umbers (squareroot of the sum of the squares of all real ad imagiary umbers). For istace i oe ca use the first or the secod colum. Normalizatio is based o value for the real part of oe compoet. Faddejev algorithm, arraged for maual calculatios ( ) : p tr( A )/ p tr( A )/ p tr( A )/ B A p I B A p I B A p A A A AB A AB I p p 4 p Determiat ad iverse matrix: det A ( ) p / / A B / p ( / ) 0 / / Characteristic equatio: f( ) (-) [ p p... p p ] [ + 4 ] 0 Eigevalues, assumed to be foud by this arbitrary order: + j j Eigevector matrix ad ormalized eigevectors: I+ B B j 0 0 j j j j j
6 . Faddejev algorithm / Matrix with equal eigevalues The example below shows i fact two cases. For all values the eigevalues are the same, ad two of them are equal. For,, the eigevectors caot be foud safely. For 0 two liearly idepedet eigevectors exist, but they caot be calculated. For 0 both eigevectors are equal. They ca be calculated, but they do ot form a liearly idepedet vector base. Faddejev algorithm, arraged for maual calculatios ( ) : p tr( A )/ p tr( A )/ p tr( A )/ B A p I B A p I B A p A A A AB A AB I p 4 p 5 p Determiat ad iverse matrix: det A ( ) p A B / p ( / ) / Characteristic equatio: f( ) (- ) [ p p... p p] [ ] 0 Eigevalues, two are equal: Eigevector matrix ad ormalized eigevectors: I+ B B 0 + ( ) ( A I) ( A I) ( A I)
7 4. Rotatio matrices / Geeral Geeral rotatio matrices are built by multiplyig sigle axis matrices. The resultig compoud matrices are uder all circumstaces orthoormal: the determiat is oe, colums are orthoormal, ad rows as well. The iverse matrix is easily foud by trasposig. A arbitrary D rotatio by three Euler agles ca be replaced by a sigle axis rotatio. The rotatio axis is calculated either by the eigevector for the eigevalue or the formula below. The axis is udefied for 0 ad 80. Sigle axis rotatio matrices: Coordiate trasformatios about x-axis, y-axis ad z-axis: x Xx x Yx x Z x X Y Z 0 0 cos( α) si( α) 0 i( α) cos( α) cos( β) 0 i( β) 0 0 si( β) 0 cos( β) cos( γ) si( γ) i( γ) cos( γ) y 4 y z 4 z x 4 x Coordiate trasformatio, Aircraft Euler agles: x4 A4x Compoud rotatio matrix, usig ψyaw, θpitch ad φroll : A X () φ Y () θ Z ( ψ) 4 4 A 4 cos( θ)cos( ψ) cos( θ)si( ψ) i( θ) cos( φ)si( ψ) + si( φ)si( θ)cos( ψ) cos( φ)cos( ψ) + si( φ)si( θ)si( ψ) si( φ)cos( θ) si( φ)si( ψ) + cos( φ)si( θ)cos( ψ) si( φ)cos( ψ) + cos( φ)si( θ)si( ψ) cos( φ)cos( θ) The Euler theorem says: Each rotatio i D by a matrix R ca be replaced by a equivalet sigle axis rotatio with a agle η about a axis. The trasformatios were itroduced above as coordiate trasformatios. T 4 For a rotatio we have to use the iverse (trasposed) matrix R A. Agle ad axis: η arccos[ 05. ( r + r + r )] 05. r r r r si( η) r r A rotatio matrix has always oe eigevalue For si( η). is the ormalized eigevector for. 0 the axis is udefied. This happes for η 0 ( o rotatio) ad η 80. 7
8 4. Rotatio matrices / Eigevalues ad eigevectors A sigle rotatio about the z-axis shows already some fudametal features of rotatio matrices. For better uderstadability the stadard eigevalue equatios are show as well. Faddejev algorithm for a rotatio matrix: Rotatio about the z-axis, usig c cos( ψ) ad s si( ψ) p + c p ( + c) p ( + c) c s 0 0 s ( + c) 0 c c c ( + c) s 0 s c 0 ( + c) c 0 0 Determiat ad iverse matrix: det A ( ) p c s A B / p c 0 A 0 0 Characteristic equatio: T f( ) (- ) [ p p... p p] [ ( + c) + ( + c) ] 0 Eigevalues, two are cojugate complex: c+ js c js Eigevector matrix ad ormalized eigevectors: I ( + c) + c s( ) + B B (( ) ( + c) + c c ( c) 0 0 c js( c) s( c) js j( c) ( c) js ( c) + js js( c) 0 ( c) + js j( c) 0 0 j s 0 s + js( c) s( c) + js s+ j( c) ( c) + js ( c) js s + js( c) 0 ( c) js s+ j( c) 0 0 j s 0 ( A I) c c s c 0 s c c js ( A I) s c 0 s js j ( c) js 0 s 0 c js (A I) s c 0 s js j ( c) + js 0 s 0 8
9 5. Newto-Bairstow algorithm The Newto-Bairstow algorithm [] fids all roots of a polyomial with real coefficiets. Polyomial with degree with real coefficiets: fx ( ) a x +a x +... a x + a 0 Fid all roots x...x for fx ( ) 0 Fid a quadratic divisor: px ( ) x rx Solve quadratic equatio, usig the Hoffma algorithm [8]. Results are two, oe or o roots. This algorithm is protected agaist divisio overflow. Two roots are either real or cojugate complex. Fid the reduced polyomial with degree m : gx ( ) bm x bx+ b0 Repeat steps util m. The iteratio is executed by a double Horer algorithm [7]. 9
10 5. Refereces [] D.K.Faddejev ud W.N.Faddejeva Numerische Methode der lieare Algebra R.Oldebourg Verlag Müche, Wie / 979 [] H.R.Schwarz Numerische Mathemati B.G.Teuber Stuttgart / 99 [] R.Zurmühl ud S.Fal Matrize ud ihre Aweduge, Teil Spriger- Verlag Berli... / 984 [4] J.Hosche ud D.Lasser Grudlage der geometrische Dateverarbeitug B.G.Teuber Stuttgart / 99 [5] G.Hoffma Euler Agles ad Projectios [6] G.Hoffma PostScript Tutor [7] G.Hoffma Matrix Library Reame as *.eps [8] G.Hoffma Solutios for a uadratic Equatio This doc Jauary 6, 0: replaced by docs-hoffma.de/ Gerot Hoffma September 9 / 007 Jauary 6 / 0 Website Load Browser / Clic here 0
Apply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.
Eigevalue-Eigevector Istructor: Nam Su Wag eigemcd Ay vector i real Euclidea space of dimesio ca be uiquely epressed as a liear combiatio of liearly idepedet vectors (ie, basis) g j, j,,, α g α g α g α
More informationChimica Inorganica 3
himica Iorgaica Irreducible Represetatios ad haracter Tables Rather tha usig geometrical operatios, it is ofte much more coveiet to employ a ew set of group elemets which are matrices ad to make the rule
More informationCHAPTER I: Vector Spaces
CHAPTER I: Vector Spaces Sectio 1: Itroductio ad Examples This first chapter is largely a review of topics you probably saw i your liear algebra course. So why cover it? (1) Not everyoe remembers everythig
More informationInverse Matrix. A meaning that matrix B is an inverse of matrix A.
Iverse Matrix Two square matrices A ad B of dimesios are called iverses to oe aother if the followig holds, AB BA I (11) The otio is dual but we ofte write 1 B A meaig that matrix B is a iverse of matrix
More informationSession 5. (1) Principal component analysis and Karhunen-Loève transformation
200 Autum semester Patter Iformatio Processig Topic 2 Image compressio by orthogoal trasformatio Sessio 5 () Pricipal compoet aalysis ad Karhue-Loève trasformatio Topic 2 of this course explais the image
More informationEigenvalues and Eigenvectors
5 Eigevalues ad Eigevectors 5.3 DIAGONALIZATION DIAGONALIZATION Example 1: Let. Fid a formula for A k, give that P 1 1 = 1 2 ad, where Solutio: The stadard formula for the iverse of a 2 2 matrix yields
More informationMATH10212 Linear Algebra B Proof Problems
MATH22 Liear Algebra Proof Problems 5 Jue 26 Each problem requests a proof of a simple statemet Problems placed lower i the list may use the results of previous oes Matrices ermiats If a b R the matrix
More informationUnit 4: Polynomial and Rational Functions
48 Uit 4: Polyomial ad Ratioal Fuctios Polyomial Fuctios A polyomial fuctio y px ( ) is a fuctio of the form p( x) ax + a x + a x +... + ax + ax+ a 1 1 1 0 where a, a 1,..., a, a1, a0are real costats ad
More informationMatrix Algebra 2.2 THE INVERSE OF A MATRIX Pearson Education, Inc.
2 Matrix Algebra 2.2 THE INVERSE OF A MATRIX MATRIX OPERATIONS A matrix A is said to be ivertible if there is a matrix C such that CA = I ad AC = I where, the idetity matrix. I = I I this case, C is a
More informationCHAPTER 5. Theory and Solution Using Matrix Techniques
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 information1 Last time: similar and diagonalizable matrices
Last time: similar ad diagoalizable matrices Let be a positive iteger Suppose A is a matrix, v R, ad λ R Recall that v a eigevector for A with eigevalue λ if v ad Av λv, or equivaletly if v is a ozero
More informationAPPENDIX F Complex Numbers
APPENDIX F Complex Numbers Operatios with Complex Numbers Complex Solutios of Quadratic Equatios Polar Form of a Complex Number Powers ad Roots of Complex Numbers Operatios with Complex Numbers Some equatios
More informationComplex Number Theory without Imaginary Number (i)
Ope Access Library Joural Complex Number Theory without Imagiary Number (i Deepak Bhalchadra Gode Directorate of Cesus Operatios, Mumbai, Idia Email: deepakm_4@rediffmail.com Received 6 July 04; revised
More informationStochastic Matrices in a Finite Field
Stochastic Matrices i a Fiite Field Abstract: I this project we will explore the properties of stochastic matrices i both the real ad the fiite fields. We first explore what properties 2 2 stochastic matrices
More informationDeterminants of order 2 and 3 were defined in Chapter 2 by the formulae (5.1)
5. Determiats 5.. Itroductio 5.2. Motivatio for the Choice of Axioms for a Determiat Fuctios 5.3. A Set of Axioms for a Determiat Fuctio 5.4. The Determiat of a Diagoal Matrix 5.5. The Determiat of a Upper
More informationNotes The Incremental Motion Model:
The Icremetal Motio Model: The Jacobia Matrix I the forward kiematics model, we saw that it was possible to relate joit agles θ, to the cofiguratio of the robot ed effector T I this sectio, we will see
More informationThe Method of Least Squares. To understand least squares fitting of data.
The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve
More informationWhere do eigenvalues/eigenvectors/eigenfunctions come from, and why are they important anyway?
Where do eigevalues/eigevectors/eigeuctios come rom, ad why are they importat ayway? I. Bacgroud (rom Ordiary Dieretial Equatios} Cosider the simplest example o a harmoic oscillator (thi o a vibratig strig)
More informationMathematics 3 Outcome 1. Vectors (9/10 pers) Lesson, Outline, Approach etc. This is page number 13. produced for TeeJay Publishers by Tom Strang
Vectors (9/0 pers) Mathematics 3 Outcome / Revise positio vector, PQ = q p, commuicative, associative, zero vector, multiplicatio by a scalar k, compoets, magitude, uit vector, (i, j, ad k) as well as
More informationx c the remainder is Pc ().
Algebra, Polyomial ad Ratioal Fuctios Page 1 K.Paulk Notes Chapter 3, Sectio 3.1 to 3.4 Summary Sectio Theorem Notes 3.1 Zeros of a Fuctio Set the fuctio to zero ad solve for x. The fuctio is zero at these
More informationAfter the completion of this section the student should recall
Chapter III Liear Algebra September 6, 7 6 CHAPTER III LINEAR ALGEBRA Objectives: After the completio of this sectio the studet should recall - the cocept of vector spaces - the operatios with vectors
More informationAppendix F: Complex Numbers
Appedix F Complex Numbers F1 Appedix F: Complex Numbers Use the imagiary uit i to write complex umbers, ad to add, subtract, ad multiply complex umbers. Fid complex solutios of quadratic equatios. Write
More information, then cv V. Differential Equations Elements of Lineaer Algebra Name: Consider the differential equation. and y2 cos( kx)
Cosider the differetial equatio y '' k y 0 has particular solutios y1 si( kx) ad y cos( kx) I geeral, ay liear combiatio of y1 ad y, cy 1 1 cy where c1, c is also a solutio to the equatio above The reaso
More information5.1. The Rayleigh s quotient. Definition 49. Let A = A be a self-adjoint matrix. quotient is the function. R(x) = x,ax, for x = 0.
40 RODICA D. COSTIN 5. The Rayleigh s priciple ad the i priciple for the eigevalues of a self-adjoit matrix Eigevalues of self-adjoit matrices are easy to calculate. This sectio shows how this is doe usig
More informationPolynomials with Rational Roots that Differ by a Non-zero Constant. Generalities
Polyomials with Ratioal Roots that Differ by a No-zero Costat Philip Gibbs The problem of fidig two polyomials P(x) ad Q(x) of a give degree i a sigle variable x that have all ratioal roots ad differ by
More informationSignal Processing in Mechatronics
Sigal Processig i Mechatroics Zhu K.P. AIS, UM. Lecture, Brief itroductio to Sigals ad Systems, Review of Liear Algebra ad Sigal Processig Related Mathematics . Brief Itroductio to Sigals What is sigal
More informationCALCULATION OF FIBONACCI VECTORS
CALCULATION OF FIBONACCI VECTORS Stuart D. Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithaca.edu ad Dai Novak Departmet of Mathematics, Ithaca College
More informationM 340L CS Homew ork Set 6 Solutions
1. Suppose P is ivertible ad M 34L CS Homew ork Set 6 Solutios A PBP 1. Solve for B i terms of P ad A. Sice A PBP 1, w e have 1 1 1 B P PBP P P AP ( ).. Suppose ( B C) D, w here B ad C are m matrices ad
More informationM 340L CS Homew ork Set 6 Solutions
. Suppose P is ivertible ad M 4L CS Homew ork Set 6 Solutios A PBP. Solve for B i terms of P ad A. Sice A PBP, w e have B P PBP P P AP ( ).. Suppose ( B C) D, w here B ad C are m matrices ad D is ivertible.
More information(VII.A) Review of Orthogonality
VII.A Review of Orthogoality At the begiig of our study of liear trasformatios i we briefly discussed projectios, rotatios ad projectios. I III.A, projectios were treated i the abstract ad without regard
More informationMatrices and vectors
Oe Matrices ad vectors This book takes for grated that readers have some previous kowledge of the calculus of real fuctios of oe real variable It would be helpful to also have some kowledge of liear algebra
More informationSymmetric Matrices and Quadratic Forms
7 Symmetric Matrices ad Quadratic Forms 7.1 DIAGONALIZAION OF SYMMERIC MARICES SYMMERIC MARIX A symmetric matrix is a matrix A such that. A = A Such a matrix is ecessarily square. Its mai diagoal etries
More informationChapter 2. Periodic points of toral. automorphisms. 2.1 General introduction
Chapter 2 Periodic poits of toral automorphisms 2.1 Geeral itroductio The automorphisms of the two-dimesioal torus are rich mathematical objects possessig iterestig geometric, algebraic, topological ad
More informationLemma Let f(x) K[x] be a separable polynomial of degree n. Then the Galois group is a subgroup of S n, the permutations of the roots.
15 Cubics, Quartics ad Polygos It is iterestig to chase through the argumets of 14 ad see how this affects solvig polyomial equatios i specific examples We make a global assumptio that the characteristic
More informationFor a 3 3 diagonal matrix we find. Thus e 1 is a eigenvector corresponding to eigenvalue λ = a 11. Thus matrix A has eigenvalues 2 and 3.
Closed Leotief Model Chapter 6 Eigevalues I a closed Leotief iput-output-model cosumptio ad productio coicide, i.e. V x = x = x Is this possible for the give techology matrix V? This is a special case
More informationMatrix Algebra from a Statistician s Perspective BIOS 524/ Scalar multiple: ka
Matrix Algebra from a Statisticia s Perspective BIOS 524/546. Matrices... Basic Termiology a a A = ( aij ) deotes a m matrix of values. Whe =, this is a am a m colum vector. Whe m= this is a row vector..2.
More information1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y
Questio (a) A square matrix A= A is called positive defiite if the quadratic form waw > 0 for every o-zero vector w [Note: Here (.) deotes the traspose of a matrix or a vector]. Let 0 A = 0 = show that:
More informationTEACHER CERTIFICATION STUDY GUIDE
COMPETENCY 1. ALGEBRA SKILL 1.1 1.1a. ALGEBRAIC STRUCTURES Kow why the real ad complex umbers are each a field, ad that particular rigs are ot fields (e.g., itegers, polyomial rigs, matrix rigs) Algebra
More information4755 Mark Scheme June Question Answer Marks Guidance M1* Attempt to find M or 108M -1 M 108 M1 A1 [6] M1 A1
4755 Mark Scheme Jue 05 * Attempt to fid M or 08M - M 08 8 4 * Divide by their determiat,, at some stage Correct determiat, (A0 for det M= 08 stated, all other OR 08 8 4 5 8 7 5 x, y,oe 8 7 4xy 8xy dep*
More informationChapter Vectors
Chapter 4. Vectors fter readig this chapter you should be able to:. defie a vector. add ad subtract vectors. fid liear combiatios of vectors ad their relatioship to a set of equatios 4. explai what it
More informationRandom Matrices with Blocks of Intermediate Scale Strongly Correlated Band Matrices
Radom Matrices with Blocks of Itermediate Scale Strogly Correlated Bad Matrices Jiayi Tog Advisor: Dr. Todd Kemp May 30, 07 Departmet of Mathematics Uiversity of Califoria, Sa Diego Cotets Itroductio Notatio
More informationREVISION SHEET FP1 (MEI) ALGEBRA. Identities In mathematics, an identity is a statement which is true for all values of the variables it contains.
the Further Mathematics etwork wwwfmetworkorguk V 07 The mai ideas are: Idetities REVISION SHEET FP (MEI) ALGEBRA Before the exam you should kow: If a expressio is a idetity the it is true for all values
More information11. FINITE FIELDS. Example 1: The following tables define addition and multiplication for a field of order 4.
11. FINITE FIELDS 11.1. A Field With 4 Elemets Probably the oly fiite fields which you ll kow about at this stage are the fields of itegers modulo a prime p, deoted by Z p. But there are others. Now although
More informationIntroduction to Optimization Techniques. How to Solve Equations
Itroductio to Optimizatio Techiques How to Solve Equatios Iterative Methods of Optimizatio Iterative methods of optimizatio Solutio of the oliear equatios resultig form a optimizatio problem is usually
More informationZeros of Polynomials
Math 160 www.timetodare.com 4.5 4.6 Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered with fidig the solutios of polyomial equatios of ay degree
More information( ) 2 + k The vertex is ( h, k) ( )( x q) The x-intercepts are x = p and x = q.
A Referece Sheet Number Sets Quadratic Fuctios Forms Form Equatio Stadard Form Vertex Form Itercept Form y ax + bx + c The x-coordiate of the vertex is x b a y a x h The axis of symmetry is x b a + k The
More informationBasic Iterative Methods. Basic Iterative Methods
Abel s heorem: he roots of a polyomial with degree greater tha or equal to 5 ad arbitrary coefficiets caot be foud with a fiite umber of operatios usig additio, subtractio, multiplicatio, divisio, ad extractio
More informationCALCULUS BASIC SUMMER REVIEW
CALCULUS BASIC SUMMER REVIEW NAME rise y y y Slope of a o vertical lie: m ru Poit Slope Equatio: y y m( ) The slope is m ad a poit o your lie is, ). ( y Slope-Itercept Equatio: y m b slope= m y-itercept=
More informationTMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods
TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A
More information8. Applications To Linear Differential Equations
8. Applicatios To Liear Differetial Equatios 8.. Itroductio 8.. Review Of Results Cocerig Liear Differetial Equatios Of First Ad Secod Orders 8.3. Eercises 8.4. Liear Differetial Equatios Of Order N 8.5.
More informationMATH 205 HOMEWORK #2 OFFICIAL SOLUTION. (f + g)(x) = f(x) + g(x) = f( x) g( x) = (f + g)( x)
MATH 205 HOMEWORK #2 OFFICIAL SOLUTION Problem 2: Do problems 7-9 o page 40 of Hoffma & Kuze. (7) We will prove this by cotradictio. Suppose that W 1 is ot cotaied i W 2 ad W 2 is ot cotaied i W 1. The
More informationRecursive Algorithms. Recurrences. Recursive Algorithms Analysis
Recursive Algorithms Recurreces Computer Sciece & Egieerig 35: Discrete Mathematics Christopher M Bourke cbourke@cseuledu A recursive algorithm is oe i which objects are defied i terms of other objects
More informationREVISION SHEET FP1 (MEI) ALGEBRA. Identities In mathematics, an identity is a statement which is true for all values of the variables it contains.
The mai ideas are: Idetities REVISION SHEET FP (MEI) ALGEBRA Before the exam you should kow: If a expressio is a idetity the it is true for all values of the variable it cotais The relatioships betwee
More information2 Geometric interpretation of complex numbers
2 Geometric iterpretatio of complex umbers 2.1 Defiitio I will start fially with a precise defiitio, assumig that such mathematical object as vector space R 2 is well familiar to the studets. Recall that
More informationComplex Numbers Solutions
Complex Numbers Solutios Joseph Zoller February 7, 06 Solutios. (009 AIME I Problem ) There is a complex umber with imagiary part 64 ad a positive iteger such that Fid. [Solutio: 697] 4i + + 4i. 4i 4i
More information6.003 Homework #3 Solutions
6.00 Homework # Solutios Problems. Complex umbers a. Evaluate the real ad imagiary parts of j j. π/ Real part = Imagiary part = 0 e Euler s formula says that j = e jπ/, so jπ/ j π/ j j = e = e. Thus the
More informationFirst, note that the LS residuals are orthogonal to the regressors. X Xb X y = 0 ( normal equations ; (k 1) ) So,
0 2. OLS Part II The OLS residuals are orthogoal to the regressors. If the model icludes a itercept, the orthogoality of the residuals ad regressors gives rise to three results, which have limited practical
More informationNotes on iteration and Newton s method. Iteration
Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f
More informationOn Nonsingularity of Saddle Point Matrices. with Vectors of Ones
Iteratioal Joural of Algebra, Vol. 2, 2008, o. 4, 197-204 O Nosigularity of Saddle Poit Matrices with Vectors of Oes Tadeusz Ostrowski Istitute of Maagemet The State Vocatioal Uiversity -400 Gorzów, Polad
More informationPAPER : IIT-JAM 2010
MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure
More informationEllipsoid Method for Linear Programming made simple
Ellipsoid Method for Liear Programmig made simple Sajeev Saxea Dept. of Computer Sciece ad Egieerig, Idia Istitute of Techology, Kapur, INDIA-08 06 December 3, 07 Abstract I this paper, ellipsoid method
More informationLECTURE 8: ORTHOGONALITY (CHAPTER 5 IN THE BOOK)
LECTURE 8: ORTHOGONALITY (CHAPTER 5 IN THE BOOK) Everythig marked by is ot required by the course syllabus I this lecture, all vector spaces is over the real umber R. All vectors i R is viewed as a colum
More informationOrthogonal transformations
Orthogoal trasformatios October 12, 2014 1 Defiig property The squared legth of a vector is give by takig the dot product of a vector with itself, v 2 v v g ij v i v j A orthogoal trasformatio is a liear
More informationLinearly Independent Sets, Bases. Review. Remarks. A set of vectors,,, in a vector space is said to be linearly independent if the vector equation
Liearly Idepedet Sets Bases p p c c p Review { v v vp} A set of vectors i a vector space is said to be liearly idepedet if the vector equatio cv + c v + + c has oly the trivial solutio = = { v v vp} The
More informationTHE EIGENVALUE PROBLEM
November 5, 2013 APPENDIX D HE EIGENVALUE PROLEM Eigevalues ad Eigevectors are Properties of the Equatios that Simulate the ehavior of a Real Structure D.1 INRODUCION he classical mathematical eigevalue
More informationWHEN IS THE (CO)SINE OF A RATIONAL ANGLE EQUAL TO A RATIONAL NUMBER? 1. My Motivation Some Sort of an Introduction
WHEN IS THE (CO)SINE OF A RATIONAL ANGLE EQUAL TO A RATIONAL NUMBER? JÖRG JAHNEL 1. My Motivatio Some Sort of a Itroductio Last term I taught Topological Groups at the Göttige Georg August Uiversity. This
More informationWhy learn matrix algebra? Vectors & Matrices with statistical applications. Brief history of linear algebra
R Vectors & Matrices with statistical applicatios x RXX RXY y RYX RYY Why lear matrix algebra? Simple way to express liear combiatios of variables ad geeral solutios of equatios. Liear statistical models
More informationMatrix Algebra 2.3 CHARACTERIZATIONS OF INVERTIBLE MATRICES Pearson Education, Inc.
2 Matrix Algebra 2.3 CHARACTERIZATIONS OF INVERTIBLE MATRICES 2012 Pearso Educatio, Ic. Theorem 8: Let A be a square matrix. The the followig statemets are equivalet. That is, for a give A, the statemets
More informationAddition: Property Name Property Description Examples. a+b = b+a. a+(b+c) = (a+b)+c
Notes for March 31 Fields: A field is a set of umbers with two (biary) operatios (usually called additio [+] ad multiplicatio [ ]) such that the followig properties hold: Additio: Name Descriptio Commutativity
More informationNANYANG TECHNOLOGICAL UNIVERSITY SYLLABUS FOR ENTRANCE EXAMINATION FOR INTERNATIONAL STUDENTS AO-LEVEL MATHEMATICS
NANYANG TECHNOLOGICAL UNIVERSITY SYLLABUS FOR ENTRANCE EXAMINATION FOR INTERNATIONAL STUDENTS AO-LEVEL MATHEMATICS STRUCTURE OF EXAMINATION PAPER. There will be oe 2-hour paper cosistig of 4 questios.
More informationNICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =
AN INTRODUCTION TO SCHRÖDER AND UNKNOWN NUMBERS NICK DUFRESNE Abstract. I this article we will itroduce two types of lattice paths, Schröder paths ad Ukow paths. We will examie differet properties of each,
More informationPROPERTIES OF AN EULER SQUARE
PROPERTIES OF N EULER SQURE bout 0 the mathematicia Leoard Euler discussed the properties a x array of letters or itegers ow kow as a Euler or Graeco-Lati Square Such squares have the property that every
More informationKinetics of Complex Reactions
Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet
More informationMEI Casio Tasks for Further Pure
Task Complex Numbers: Roots of Quadratic Equatios. Add a ew Equatio scree: paf 2. Chage the Complex output to a+bi: LpNNNNwd 3. Select Polyomial ad set the Degree to 2: wq 4. Set a=, b=5 ad c=6: l5l6l
More informationPatterns in Complex Numbers An analytical paper on the roots of a complex numbers and its geometry
IB MATHS HL POTFOLIO TYPE Patters i Complex Numbers A aalytical paper o the roots of a complex umbers ad its geometry i Syed Tousif Ahmed Cadidate Sessio Number: 0066-009 School Code: 0066 Sessio: May
More informationApplication of Jordan Canonical Form
CHAPTER 6 Applicatio of Jorda Caoical Form Notatios R is the set of real umbers C is the set of complex umbers Q is the set of ratioal umbers Z is the set of itegers N is the set of o-egative itegers Z
More information6. Kalman filter implementation for linear algebraic equations. Karhunen-Loeve decomposition
6. Kalma filter implemetatio for liear algebraic equatios. Karhue-Loeve decompositio 6.1. Solvable liear algebraic systems. Probabilistic iterpretatio. Let A be a quadratic matrix (ot obligatory osigular.
More informationState Space Representation
Optimal Cotrol, Guidace ad Estimatio Lecture 2 Overview of SS Approach ad Matrix heory Prof. Radhakat Padhi Dept. of Aerospace Egieerig Idia Istitute of Sciece - Bagalore State Space Represetatio Prof.
More information3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials
Math 60 www.timetodare.com 3. Properties of Divisio 3.3 Zeros of Polyomials 3.4 Complex ad Ratioal Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered
More informationMath 475, Problem Set #12: Answers
Math 475, Problem Set #12: Aswers A. Chapter 8, problem 12, parts (b) ad (d). (b) S # (, 2) = 2 2, sice, from amog the 2 ways of puttig elemets ito 2 distiguishable boxes, exactly 2 of them result i oe
More informationPhysics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.
Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;
More informationLINEAR RECURSION RELATIONS - LESSON FOUR SECOND-ORDER LINEAR RECURSION RELATIONS
LINEAR RECURSION RELATIONS - LESSON FOUR SECOND-ORDER LINEAR RECURSION RELATIONS BROTHER ALFRED BROUSSEAU St. Mary's College, Califoria Give a secod-order liear recursio relatio (.1) T. 1 = a T + b T 1,
More informationa for a 1 1 matrix. a b a b 2 2 matrix: We define det ad bc 3 3 matrix: We define a a a a a a a a a a a a a a a a a a
Math S-b Lecture # Notes This wee is all about determiats We ll discuss how to defie them, how to calculate them, lear the allimportat property ow as multiliearity, ad show that a square matrix A is ivertible
More informationPresentation of complex number in Cartesian and polar coordinate system
a + bi, aεr, bεr i = z = a + bi a = Re(z), b = Im(z) give z = a + bi & w = c + di, a + bi = c + di a = c & b = d The complex cojugate of z = a + bi is z = a bi The sum of complex cojugates is real: z +
More informationCourse : Algebraic Combinatorics
Course 18.312: Algebraic Combiatorics Lecture Notes # 18-19 Addedum by Gregg Musier March 18th - 20th, 2009 The followig material ca be foud i a umber of sources, icludig Sectios 7.3 7.5, 7.7, 7.10 7.11,
More informationA TYPE OF PRIMITIVE ALGEBRA*
A TYPE OF PRIMITIVE ALGEBRA* BT J. H. M. WEDDERBURN I a recet paper,t L. E. Dickso has discussed the liear associative algebra, A, defied by the relatios xy = yo(x), y = g, where 8 ( x ) is a polyomial
More informationWe will conclude the chapter with the study a few methods and techniques which are useful
Chapter : Coordiate geometry: I this chapter we will lear about the mai priciples of graphig i a dimesioal (D) Cartesia system of coordiates. We will focus o drawig lies ad the characteristics of the graphs
More informationECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations
ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows
More informationG r a d e 1 1 P r e - C a l c u l u s M a t h e m a t i c s ( 3 0 S )
G r a d e 1 1 P r e - C a l c u l u s M a t h e m a t i c s ( 3 0 S ) Grade 11 Pre-Calculus Mathematics (30S) is desiged for studets who ited to study calculus ad related mathematics as part of post-secodary
More informationLinear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy
Liear Differetial Equatios of Higher Order Basic Theory: Iitial-Value Problems d y d y dy Solve: a( ) + a ( )... a ( ) a0( ) y g( ) + + + = d d d ( ) Subject to: y( 0) = y0, y ( 0) = y,..., y ( 0) = y
More information(I.C) Matrix algebra
(IC) Matrix algebra Before formalizig Gauss-Jorda i terms of a fixed procedure for row-reducig A, we briefly review some properties of matrix multiplicatio Let m{ [A ij ], { [B jk ] p, p{ [C kl ] q be
More informationAH Checklist (Unit 3) AH Checklist (Unit 3) Matrices
AH Checklist (Uit 3) AH Checklist (Uit 3) Matrices Skill Achieved? Kow that a matrix is a rectagular array of umbers (aka etries or elemets) i paretheses, each etry beig i a particular row ad colum Kow
More information(3) If you replace row i of A by its sum with a multiple of another row, then the determinant is unchanged! Expand across the i th row:
Math 50-004 Tue Feb 4 Cotiue with sectio 36 Determiats The effective way to compute determiats for larger-sized matrices without lots of zeroes is to ot use the defiitio, but rather to use the followig
More informationR is a scalar defined as follows:
Math 8. Notes o Dot Product, Cross Product, Plaes, Area, ad Volumes This lecture focuses primarily o the dot product ad its may applicatios, especially i the measuremet of agles ad scalar projectio ad
More informationLinear Programming and the Simplex Method
Liear Programmig ad the Simplex ethod Abstract This article is a itroductio to Liear Programmig ad usig Simplex method for solvig LP problems i primal form. What is Liear Programmig? Liear Programmig is
More informationt distribution [34] : used to test a mean against an hypothesized value (H 0 : µ = µ 0 ) or the difference
EXST30 Backgroud material Page From the textbook The Statistical Sleuth Mea [0]: I your text the word mea deotes a populatio mea (µ) while the work average deotes a sample average ( ). Variace [0]: The
More informationAlgorithms and Data Structures 2014 Exercises and Solutions Week 13
Algorithms ad Data Structures 204 Exercises ad Solutios Week 3 Toom-Cook (cotiued) Durig the last lecture, two polyomials A(x) a 0 + a x ad B(x) b 0 + b x both of degree were multiplied, first by evaluatig
More informationAssignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.
PHYSICS 34 QUANTUM PHYSICS II (25) Assigmet 2 Solutios 1. With respect to a pair of orthoormal vectors ϕ 1 ad ϕ 2 that spa the Hilbert space H of a certai system, the operator  is defied by its actio
More informationSection Let A =. Then A has characteristic equation λ. 2 4λ + 3 = 0 or (λ 3)(λ 1) = 0. Hence the eigenvalues of A are λ 1 = 3 and λ 2 = 1.
Sectio 63 4 3 Let A The A has characteristic equatio λ 2 4λ + 3 or (λ 3)(λ ) Hece the eigevalues of A are λ 3 ad λ 2 λ 3 The correspodig eigevectors satisfy (A λ I 2 )X, or 3 3 or equivaletly x 3y Hece
More informationCALCULATING FIBONACCI VECTORS
THE GENERALIZED BINET FORMULA FOR CALCULATING FIBONACCI VECTORS Stuart D Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithacaedu ad Dai Novak Departmet
More information