Tridiagonal reduction redux
|
|
- Corey Wiggins
- 5 years ago
- Views:
Transcription
1 Week 15: Moday, Nov 27 Wedesday, Nov 29 Tridiagoal reductio redux Cosider the problem of computig eigevalues of a symmetric matrix A that is large ad sparse Usig the grail code i LAPACK, we ca compute all the eigevalues of a -by- tridiagoal matrix i O() time; but the usual Householder-based algorithm to reduce A to tridiagoal form costs O( 3 ) Though it s possible to maitai sparsity i some cases, it is ot altogether straightforward So let s see if we ca try somethig differet What we wat is a orthogoal matrix Q such that AQ = QT, where α 1 β 1 β 1 α 2 β 2 T = β 2 α 3 β 1 β 1 Writig colum j of this matrix equatio gives us α Aq j = β j 1 q j 1 + α j q j + β j q j+1, where α j = q T j Aq j ad β j 1 = q T j 1Aq j If we rearrage this slightly, we have r j+1 = Aq j q j α j q j 1 β j 1 = Aq j (q j q T j )Aq j (q j 1 q T j 1)Aq j q j+1 = r j+1 /β j That is, q j+1 is exactly what we compute if we orthogoalize Aq j agaist the vectors q j ad q j 1 Note that Aq j is automatically orthogoal to q 1,, q j 2, sice q T k Aq j = T kj is zero for k < j 1 This process of icremetally producig a orthoormal basis by multiplyig the last vector by A ad orthogoalizig the result agaist what came before gives us the basic Laczos procedure If we start with q 1 = e 1, the Laczos procedure will give us the same T matrix 1 as the Householder 1 Well, almost the same T matrix Note that the off-diagoal elemets β j might be positive or egative i the Householder tridiagoalizatio, but they are always positive i the Laczos algorithm However, the two T matrices satisfy a similarity relatio with a diagoal matrix whose diagoal etries are ±1
2 algorithm, at least i exact arithmetic But while a step i the Householder algorithm costs O( 2 ), the cost of a Laczos step is a matrix-vector multiply (which may cost less tha O( 2 ) if A is sparse) followed by O() work i dot products ad gaxpy operatios Also, ote that the Laczos step requires oly a few vectors of itermediate storage I exact arithmetic, we could iterate util we foud β k = 0, which would correspod to havig a ivariat subspace spaed by {q 1,, q k } Except for very special choices of startig vectors, though, this would happe oly for k = Partial tridiagoalizatio ad Krylov subspaces What happes whe we ru oly a few steps of the Laczos iteratio? After m steps, we have computed m colums of Q ad m rows of T ; ad there is some hope that we might be able to use the leadig m-by-m subblock of T (ie the block Rayleigh quotiet with the first m colums of Q) i order to extract some useful iformatio about the spectrum of A The reaso why we might have such a hope is that we otice that the Krylov subspace spaed by q 1,, q m is also spaed by the first m iterates of a power method: K m (A, q 1 ) = spa{q 1, q 2,, q m } = spa{q 1, Aq 1,, A m 1 q 1 } As we have previously discussed, the power basis is ot a particularly useful basis for umerical computatios, but it is useful for thikig about Krylov subspaces I particular, because it cotais the mth vector produced by power iteratio startig from q 1, we expect that the largest Ritz value associated with K m (A, q 1 ) (ie the largest magitude eigevalue of the block Rayleigh quotiet (Q :,1:m ) T AQ :,1:m ) will be at least as large as the Rayleigh quotiet estimate produced by m steps of power iteratio Of course, the Krylov subspace K m (A, q 1 ) has a lot more iformatio tha just the vector from step m of a power iteratio For example, otice that if σ is ay shift, the K m (A σi, q 1 ) = K m (A, q 1 ); that is, Krylov subspaces are shift ivariat This is relevat because addig a shift to A ca chage which of the exterior eigevalues i the spectrum is domiat So if λ 1 λ 2 λ, the we expect to be able to extract good estimates to λ 1 ad λ from a few steps of Laczos iteratio uder the assumptio that power iteratio with certai shifted matrices coverges sufficietly fast
3 Laczos i iexact arithmetic The observat reader will have oticed somethig usettlig about the basic Laczos iteratio: at the heart of the Laczos iteratio is a Gram- Schmidt orthogoalizatio step, which we said i our discussio of QR factorizatios could become umerically ustable I particular, whe r k+1 = (A α k )q k + q k 1 β k 1 is small (ie whe β k is small relative to A ), cacellatio may cause the computed vector to cosisty mostly of roudoff, so that it is o loger orthogoal to the precedig vectors But a small value of β k correspods exactly to covergece to a good approximatio to a ivariat subspace! Thus, i floatig poit arithmetic, the basic Laczos iteratio teds to ru reasoably well util covergece, at which poit it effectively restarts with a radom startig vector which is made up primarily of roudoff error So we get covergece of the Ritz values of the computed T to the first few exterior eigevalues; ad the more ghost Ritz values coverge to the first few exterior eigevalues; ad we cotiue o i this fashio for as log as we are willig to let thigs ru This ot state of affairs is ot etirely satisfactory, but we ca fix it by re-orthogoalizatio: that is, we improve the orthogoality of the computed basis by explicitly orthogoalizig Aq k agaist vectors to which it should already be orthoormal i exact arithmetic Oe approach is to use complete orthogoalizatio: compute Aq j ad the orthogoalize agaist all previous vectors (eg by a sequece of Householder trasformatios) This is expesive i terms of storage (ad data movemet) ad arithmetic, but it does work A less expesive solutio, selective orthogoalizatio, ivolves orthogoalizig agaist coverged eigevector estimates (Ritz vectors) I the iterest of time, we will ot discuss reorthogoalizatio further If you are iterested, though or if you are iterested i ay umber of other aspects of Laczos iteratio that we skip over I recommed lookig at Parlett s book The Symmetric Eigevalue Problem Stewart s book Matrix Algorithms, Vol 2: Eigesystems also has a ice treatmet Why did umerical istability ot arise whe we discussed CG? It s there behid the scees, ad it explais why ruig CG for steps geerally does ot provide the exact aswer to a liear system But obody sesible uses CG as a direct method ayhow; it is used as a iterative solver, ad it turs out that the loss of orthogoality does ot matter that much to the CG iteratio cotiuig to make progress
4 Partial tridiagoalizatio ad residual bouds Suppose AQ = QT, where T ii = α i, T i,i+1 = β i If we take m steps of Laczos iteratio, we geerate Q :,1:m = [ q 1 q 2 q m ] as well as the first m coefficiets (α i ) m i=1 ad (β i ) m i=1 Let us deote Q m ad T m as the leadig m colums of Q ad the leadig m m submatrix of T respectively; the writig the first m colums of AQ ad QT gives us (1) AQ m = Q m T m + β m q m+1 e T m Here T m is a block Rayleigh quotiet T m = Q T maq m, ad the eigevalues of T m (the Ritz values) are used to approximate the eigevalues of A Now cosider the eigedecompositio T m = Y ΘY T where Y T Y = I ad Θ = diag(θ 1,, θ m ) The postmultiplyig (1) by Y gives AQ m Y = Q m Y Θ + β m q m+1 e T my The colums of Z = Q m Y = [ z 1 z m ] are approximate eigevalues correspodig to the approximate eigevalues θ 1, θ 2,, θ m For each colum, we have Az k z k θ k = β m q m+1 e T my k, which meas Az k z k θ k 2 = β m e T my k This is useful because, as we discussed before, i the symmetric case a small residual error implies a small distace to the closest eigevalue This is also useful because the residual error ca be computed with o further matrix operatios we eed oly to look at quatities that we would already compute i the process of obtaiig the tridiagoal coefficiets ad the correspodig Ritz values I particular, ote that we ca compute the residual for the Ritz pair (approximate eigepair) (z k, θ k ) without explicitly computig z k! The polyomial coectio Suppose λ 1 λ 2 λ are the eigevalues of A, with correspodig eigevectors w 1 through w I matrix form, the, A = W ΛW T
5 Now, suppose we ru m steps of Laczos iteratio (i exact arithmetic), ad let θ 1 be the largest eigevector of T m (the largest Ritz value) This is the same as sayig that θ 1 maximizes the Rayleigh quotiet over the Krylov subspace K m (A, q 1 ): θ 1 = max y R m y T T m y y T y (Q m y) T A(Q m y) = max y R m (Q m y) T (Q m y) = max z T Az z K m(a,q 1 ) z T z Now, ay vector z K m (A, q 1 ) ca be writte as z = c 0 q 1 + c 1 Aq c m 1 A m 1 q 1 = p(a)q 1, where p is a polyomial of degree at most m 1 Thus, θ 1 = d T p(λ)λp(λ)d max deg(p)<m d T p(λ) 2 d q1 T p(a)ap(a)q 1 max deg(p)<m q1 T p(a) 2 q 1 This is still somewhat awkward, so let us rewrite thigs i terms of the eigevector basis Defie d = Z T q 1 ; the i=1 θ 1 = d2 i p(λ i ) 2 λ i i=1 d2 i p(λ i) 2 = max deg(p)<m Let s give a ame to the fuctio i this maximizatio: φ(p) = i=1 d2 i p(λ i ) 2 λ i i=1 d2 i p(λ i) 2 = λ 1 d2 i p(λ i ) 2 (λ 1 λ i ) i=1 d2 i p(λ λ i) 2 1 (λ 1 λ ) d2 i p(λ i ) 2 d 2 1p(λ 1 ) 2 We kow that θ 1 = max deg(p)<m φ(p) λ 1 ; we would also like a lower boud o θ 1 If we ca show that this lower boud approaches λ 1 at some rate with icreasig m, the we kow θ 1 will coverge to λ 1 at least as fast Note that if we foud a polyomial that was zero at λ 2,, λ ad ozero at λ 1, the we would recover λ 1 exactly But such a polyomial usually has too high a degree What we would like i order to get a boud, the, is a degree m polyomial which is ot too big o [λ 2, λ ], but is relatively large at λ 2 Havig see it s appearace i our aalysis of Chebyshev iteratio ad CG, you will ot be surprised that a good cadidate is our old fried, the rescaled Chebyshev polyomial ˆT k (x) = T k ( x λ λ 2 λ )
6 Note that ˆT k (λ j ) 1 for j > 1 ad ˆT k (λ 1 ) grows rapidly with k Note that for ay p, we have φ(p) = λ 1 d2 i p(λ i ) 2 (λ 1 λ i ) i=1 d2 i p(λ i) 2 λ 1 (λ 1 λ ) d2 i p(λ i ) 2 d 2 1p(λ 1 ) 2 Now, ote that i=1 d2 i = d 2 2 = q = 1, that ˆT m 1 (λ i ) 2 1 for i = 2,, ad that ˆT m 1 (λ 1 ) = T c m 1 (1 + 2ρ 1 ) Thus, φ(ĉ m 1 ) λ 1 (λ 1 λ ) d2 ˆT i m 1 (λ i ) 2 1 d 2 1 d 2 ˆT λ 1 (λ 1 λ ) 1 m 1 (λ 1 ) 2 d 2 1T m 1 (1 + 2ρ 1 ) 2 The etry d 1 = w T 1 q 1 is the cosie of the agle φ 1 betwee the first eigevector w 1 ad the startig vector q 1 ; so we ca write (1 d 2 1)/d 2 1 = ta(φ 1 ) 2, which gives us the fial boud λ 1 θ 1 = max φ(p) φ(ĉ m 1)λ 1 (λ 1 λ ) ta(φ 1 ) 2 deg p<m T m 1 (1 + 2ρ 1 ) 2 where ρ 1 = (λ 1 λ 2 )/(λ 2 λ ) We otice a few thigs from this boud The rate of covergece (determied by T m 1 (1 + 2ρ 1 ) 2 ) is strictly better tha the rate of covergece for power iteratio, which makes sese sice we are usig strictly more iformatio i the Laczos iteratio tha we use i the power iteratio (we maximize over a m-dimesioal subspace rather tha a 1-dimesioal space) Also, the trick ivolved is a good oe: istead of reasoig about matrices, reaso about polyomials But, as with our aalysis of CG, this boud may be very pessimistic, sice it does ot take ito accout the distributio of eigevalues i [λ 2, λ ] For example, if may eigevalues of A cluster ear zero (as happes with certai discretizatios of compact operators, for istace), the we might get much better covergece tha idicated by a basic boud that works for ay distributio of eigevalues betwee [λ 2, λ ]
Basic 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 informationLecture 8: October 20, Applications of SVD: least squares approximation
Mathematical Toolkit Autum 2016 Lecturer: Madhur Tulsiai Lecture 8: October 20, 2016 1 Applicatios of SVD: least squares approximatio We discuss aother applicatio of sigular value decompositio (SVD) of
More informationInfinite Sequences and Series
Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet
More informationMath 155 (Lecture 3)
Math 55 (Lecture 3) September 8, I this lecture, we ll cosider the aswer to oe of the most basic coutig problems i combiatorics Questio How may ways are there to choose a -elemet subset of the set {,,,
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 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 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 information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationGeometry of LS. LECTURE 3 GEOMETRY OF LS, PROPERTIES OF σ 2, PARTITIONED REGRESSION, GOODNESS OF FIT
OCTOBER 7, 2016 LECTURE 3 GEOMETRY OF LS, PROPERTIES OF σ 2, PARTITIONED REGRESSION, GOODNESS OF FIT Geometry of LS We ca thik of y ad the colums of X as members of the -dimesioal Euclidea space R Oe ca
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 informationAlgebra of Least Squares
October 19, 2018 Algebra of Least Squares Geometry of Least Squares Recall that out data is like a table [Y X] where Y collects observatios o the depedet variable Y ad X collects observatios o the k-dimesioal
More informationTheorem: Let A n n. In this case that A does reduce to I, we search for A 1 as the solution matrix X to the matrix equation A X = I i.e.
Theorem: Let A be a square matrix The A has a iverse matrix if ad oly if its reduced row echelo form is the idetity I this case the algorithm illustrated o the previous page will always yield the iverse
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(A sequence also can be thought of as the list of function values attained for a function f :ℵ X, where f (n) = x n for n 1.) x 1 x N +k x N +4 x 3
MATH 337 Sequeces Dr. Neal, WKU Let X be a metric space with distace fuctio d. We shall defie the geeral cocept of sequece ad limit i a metric space, the apply the results i particular to some special
More informationSeptember 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1
September 0 s (Edecel) Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright
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 informationApply 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 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 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 informationIntroduction to Machine Learning DIS10
CS 189 Fall 017 Itroductio to Machie Learig DIS10 1 Fu with Lagrage Multipliers (a) Miimize the fuctio such that f (x,y) = x + y x + y = 3. Solutio: The Lagragia is: L(x,y,λ) = x + y + λ(x + y 3) Takig
More informationProperties and Tests of Zeros of Polynomial Functions
Properties ad Tests of Zeros of Polyomial Fuctios The Remaider ad Factor Theorems: Sythetic divisio ca be used to fid the values of polyomials i a sometimes easier way tha substitutio. This is show by
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 informationMath 61CM - Solutions to homework 3
Math 6CM - Solutios to homework 3 Cédric De Groote October 2 th, 208 Problem : Let F be a field, m 0 a fixed oegative iteger ad let V = {a 0 + a x + + a m x m a 0,, a m F} be the vector space cosistig
More informationPolynomial Functions and Their Graphs
Polyomial Fuctios ad Their Graphs I this sectio we begi the study of fuctios defied by polyomial expressios. Polyomial ad ratioal fuctios are the most commo fuctios used to model data, ad are used extesively
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 informationLINEAR ALGEBRA. Paul Dawkins
LINEAR ALGEBRA Paul Dawkis Table of Cotets Preface... ii Outlie... iii Systems of Equatios ad Matrices... Itroductio... Systems of Equatios... Solvig Systems of Equatios... 5 Matrices... 7 Matrix Arithmetic
More informationSequences, Series, and All That
Chapter Te Sequeces, Series, ad All That. Itroductio Suppose we wat to compute a approximatio of the umber e by usig the Taylor polyomial p for f ( x) = e x at a =. This polyomial is easily see to be 3
More informationMathematics review for CSCI 303 Spring Department of Computer Science College of William & Mary Robert Michael Lewis
Mathematics review for CSCI 303 Sprig 019 Departmet of Computer Sciece College of William & Mary Robert Michael Lewis Copyright 018 019 Robert Michael Lewis Versio geerated: 13 : 00 Jauary 17, 019 Cotets
More informationCHAPTER 3. GOE and GUE
CHAPTER 3 GOE ad GUE We quicly recall that a GUE matrix ca be defied i the followig three equivalet ways. We leave it to the reader to mae the three aalogous statemets for GOE. I the previous chapters,
More information5.1 Review of Singular Value Decomposition (SVD)
MGMT 69000: Topics i High-dimesioal Data Aalysis Falll 06 Lecture 5: Spectral Clusterig: Overview (cotd) ad Aalysis Lecturer: Jiamig Xu Scribe: Adarsh Barik, Taotao He, September 3, 06 Outlie Review of
More informationSequences. Notation. Convergence of a Sequence
Sequeces A sequece is essetially just a list. Defiitio (Sequece of Real Numbers). A sequece of real umbers is a fuctio Z (, ) R for some real umber. Do t let the descriptio of the domai cofuse you; it
More informationA sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as
More informationCS276A Practice Problem Set 1 Solutions
CS76A Practice Problem Set Solutios Problem. (i) (ii) 8 (iii) 6 Compute the gamma-codes for the followig itegers: (i) (ii) 8 (iii) 6 Problem. For this problem, we will be dealig with a collectio of millio
More informationProblems from 9th edition of Probability and Statistical Inference by Hogg, Tanis and Zimmerman:
Math 224 Fall 2017 Homework 4 Drew Armstrog Problems from 9th editio of Probability ad Statistical Iferece by Hogg, Tais ad Zimmerma: Sectio 2.3, Exercises 16(a,d),18. Sectio 2.4, Exercises 13, 14. Sectio
More informationOptimization Methods MIT 2.098/6.255/ Final exam
Optimizatio Methods MIT 2.098/6.255/15.093 Fial exam Date Give: December 19th, 2006 P1. [30 pts] Classify the followig statemets as true or false. All aswers must be well-justified, either through a short
More informationRecurrence Relations
Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial(-)); } Let t be the umber of multiplicatios eeded to calculate factorial(). The
More informationThe Binomial Theorem
The Biomial Theorem Robert Marti Itroductio The Biomial Theorem is used to expad biomials, that is, brackets cosistig of two distict terms The formula for the Biomial Theorem is as follows: (a + b ( k
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 informationSequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet
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 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 informationThe Growth of Functions. Theoretical Supplement
The Growth of Fuctios Theoretical Supplemet The Triagle Iequality The triagle iequality is a algebraic tool that is ofte useful i maipulatig absolute values of fuctios. The triagle iequality says that
More informationLesson 10: Limits and Continuity
www.scimsacademy.com Lesso 10: Limits ad Cotiuity SCIMS Academy 1 Limit of a fuctio The cocept of limit of a fuctio is cetral to all other cocepts i calculus (like cotiuity, derivative, defiite itegrals
More informationReview Problems 1. ICME and MS&E Refresher Course September 19, 2011 B = C = AB = A = A 2 = A 3... C 2 = C 3 = =
Review Problems ICME ad MS&E Refresher Course September 9, 0 Warm-up problems. For the followig matrices A = 0 B = C = AB = 0 fid all powers A,A 3,(which is A times A),... ad B,B 3,... ad C,C 3,... Solutio:
More informationSection A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics
X0/70 NATIONAL QUALIFICATIONS 005 MONDAY, MAY.00 PM 4.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.
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 E-2b Lecture #8 Notes This week is all about determiats. We ll discuss how to defie them, how to calculate them, lear the allimportat property kow as multiliearity, ad show that a square matrix A
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 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 informationThe Discrete Fourier Transform
The Discrete Fourier Trasform Complex Fourier Series Represetatio Recall that a Fourier series has the form a 0 + a k cos(kt) + k=1 b k si(kt) This represetatio seems a bit awkward, sice it ivolves two
More informationMIDTERM 3 CALCULUS 2. Monday, December 3, :15 PM to 6:45 PM. Name PRACTICE EXAM SOLUTIONS
MIDTERM 3 CALCULUS MATH 300 FALL 08 Moday, December 3, 08 5:5 PM to 6:45 PM Name PRACTICE EXAM S Please aswer all of the questios, ad show your work. You must explai your aswers to get credit. You will
More informationAnalysis of Algorithms. Introduction. Contents
Itroductio The focus of this module is mathematical aspects of algorithms. Our mai focus is aalysis of algorithms, which meas evaluatig efficiecy of algorithms by aalytical ad mathematical methods. We
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 information1 Approximating Integrals using Taylor Polynomials
Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................
More informationMATH 304: MIDTERM EXAM SOLUTIONS
MATH 304: MIDTERM EXAM SOLUTIONS [The problems are each worth five poits, except for problem 8, which is worth 8 poits. Thus there are 43 possible poits.] 1. Use the Euclidea algorithm to fid the greatest
More informationAverage-Case Analysis of QuickSort
Average-Case Aalysis of QuickSort Comp 363 Fall Semester 003 October 3, 003 The purpose of this documet is to itroduce the idea of usig recurrece relatios to do average-case aalysis. The average-case ruig
More informationMachine Learning for Data Science (CS 4786)
Machie Learig for Data Sciece CS 4786) Lecture 9: Pricipal Compoet Aalysis The text i black outlies mai ideas to retai from the lecture. The text i blue give a deeper uderstadig of how we derive or get
More informationNotes for Lecture 11
U.C. Berkeley CS78: Computatioal Complexity Hadout N Professor Luca Trevisa 3/4/008 Notes for Lecture Eigevalues, Expasio, ad Radom Walks As usual by ow, let G = (V, E) be a udirected d-regular graph with
More informationMA131 - Analysis 1. Workbook 3 Sequences II
MA3 - Aalysis Workbook 3 Sequeces II Autum 2004 Cotets 2.8 Coverget Sequeces........................ 2.9 Algebra of Limits......................... 2 2.0 Further Useful Results........................
More informationLinear regression. Daniel Hsu (COMS 4771) (y i x T i β)2 2πσ. 2 2σ 2. 1 n. (x T i β y i ) 2. 1 ˆβ arg min. β R n d
Liear regressio Daiel Hsu (COMS 477) Maximum likelihood estimatio Oe of the simplest liear regressio models is the followig: (X, Y ),..., (X, Y ), (X, Y ) are iid radom pairs takig values i R d R, ad Y
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 informationMath 2784 (or 2794W) University of Connecticut
ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really
More informationv = -!g(x 0 ) Ûg Ûx 1 Ûx 2 Ú If we work out the details in the partial derivatives, we get a pleasing result. n Ûx k, i x i - 2 b k
The Method of Steepest Descet This is the quadratic fuctio from to that is costructed to have a miimum at the x that solves the system A x = b: g(x) = - 2 I the method of steepest descet, we
More informationFeedback in Iterative Algorithms
Feedback i Iterative Algorithms Charles Byre (Charles Byre@uml.edu), Departmet of Mathematical Scieces, Uiversity of Massachusetts Lowell, Lowell, MA 01854 October 17, 2005 Abstract Whe the oegative system
More informationMachine Learning Theory Tübingen University, WS 2016/2017 Lecture 12
Machie Learig Theory Tübige Uiversity, WS 06/07 Lecture Tolstikhi Ilya Abstract I this lecture we derive risk bouds for kerel methods. We will start by showig that Soft Margi kerel SVM correspods to miimizig
More informationChapter 10: Power Series
Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because
More informationMachine Learning for Data Science (CS 4786)
Machie Learig for Data Sciece CS 4786) Lecture & 3: Pricipal Compoet Aalysis The text i black outlies high level ideas. The text i blue provides simple mathematical details to derive or get to the algorithm
More information6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.
6 Itegers Modulo I Example 2.3(e), we have defied the cogruece of two itegers a,b with respect to a modulus. Let us recall that a b (mod ) meas a b. We have proved that cogruece is a equivalece relatio
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 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 informationa 3, a 4, ... are the terms of the sequence. The number a n is the nth term of the sequence, and the entire sequence is denoted by a n
60_090.qxd //0 : PM Page 59 59 CHAPTER 9 Ifiite Series Sectio 9. EXPLORATION Fidig Patters Describe a patter for each of the followig sequeces. The use your descriptio to write a formula for the th term
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 informationf(x) dx as we do. 2x dx x also diverges. Solution: We compute 2x dx lim
Math 3, Sectio 2. (25 poits) Why we defie f(x) dx as we do. (a) Show that the improper itegral diverges. Hece the improper itegral x 2 + x 2 + b also diverges. Solutio: We compute x 2 + = lim b x 2 + =
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 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 informationPosted-Price, Sealed-Bid Auctions
Posted-Price, Sealed-Bid Auctios Professors Greewald ad Oyakawa 207-02-08 We itroduce the posted-price, sealed-bid auctio. This auctio format itroduces the idea of approximatios. We describe how well this
More informationw (1) ˆx w (1) x (1) /ρ and w (2) ˆx w (2) x (2) /ρ.
2 5. Weighted umber of late jobs 5.1. Release dates ad due dates: maximimizig the weight of o-time jobs Oce we add release dates, miimizig the umber of late jobs becomes a sigificatly harder problem. For
More informationEstimation of Backward Perturbation Bounds For Linear Least Squares Problem
dvaced Sciece ad Techology Letters Vol.53 (ITS 4), pp.47-476 http://dx.doi.org/.457/astl.4.53.96 Estimatio of Bacward Perturbatio Bouds For Liear Least Squares Problem Xixiu Li School of Natural Scieces,
More information1 Generating functions for balls in boxes
Math 566 Fall 05 Some otes o geeratig fuctios Give a sequece a 0, a, a,..., a,..., a geeratig fuctio some way of represetig the sequece as a fuctio. There are may ways to do this, with the most commo ways
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 informationMath Solutions to homework 6
Math 175 - Solutios to homework 6 Cédric De Groote November 16, 2017 Problem 1 (8.11 i the book): Let K be a compact Hermitia operator o a Hilbert space H ad let the kerel of K be {0}. Show that there
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 5-4 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 facts,
More informationDirichlet s Theorem on Arithmetic Progressions
Dirichlet s Theorem o Arithmetic Progressios Athoy Várilly Harvard Uiversity, Cambridge, MA 0238 Itroductio Dirichlet s theorem o arithmetic progressios is a gem of umber theory. A great part of its beauty
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 information7.1 Convergence of sequences of random variables
Chapter 7 Limit Theorems Throughout this sectio we will assume a probability space (, F, P), i which is defied a ifiite sequece of radom variables (X ) ad a radom variable X. The fact that for every ifiite
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 informationSection 6.4: Series. Section 6.4 Series 413
ectio 64 eries 4 ectio 64: eries A couple decides to start a college fud for their daughter They pla to ivest $50 i the fud each moth The fud pays 6% aual iterest, compouded mothly How much moey will they
More informationThe picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled
1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how
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 informationFFTs in Graphics and Vision. The Fast Fourier Transform
FFTs i Graphics ad Visio The Fast Fourier Trasform 1 Outlie The FFT Algorithm Applicatios i 1D Multi-Dimesioal FFTs More Applicatios Real FFTs 2 Computatioal Complexity To compute the movig dot-product
More informationQuantum Information & Quantum Computation
CS9A, Sprig 5: Quatum Iformatio & Quatum Computatio Wim va Dam Egieerig, Room 59 vadam@cs http://www.cs.ucsb.edu/~vadam/teachig/cs9/ Admiistrivia Do the exercises. Aswers will be posted at the ed of the
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More information18.440, March 9, Stirling s formula
Stirlig s formula 8.44, March 9, 9 The factorial fuctio! is importat i evaluatig biomial, hypergeometric, ad other probabilities. If is ot too large,! ca be computed directly, by calculators or computers.
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 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 information5.3 Preconditioning. - M is easy to deal with in parallel (reduced approximate direct solver)
5.3 Precoditioig Direct solvers: Sequetial, loosig sparsity terative solvers: easy parallel ad sparse, but possibly slowly coverget Combiatio of both methods: clude precoditioer i the form - x - b, such
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More information18.S096: Homework Problem Set 1 (revised)
8.S096: Homework Problem Set (revised) Topics i Mathematics of Data Sciece (Fall 05) Afoso S. Badeira Due o October 6, 05 Exteded to: October 8, 05 This homework problem set is due o October 6, at the
More informationDiscrete Mathematics for CS Spring 2007 Luca Trevisan Lecture 22
CS 70 Discrete Mathematics for CS Sprig 2007 Luca Trevisa Lecture 22 Aother Importat Distributio The Geometric Distributio Questio: A biased coi with Heads probability p is tossed repeatedly util the first
More informationFrequency Domain Filtering
Frequecy Domai Filterig Raga Rodrigo October 19, 2010 Outlie Cotets 1 Itroductio 1 2 Fourier Represetatio of Fiite-Duratio Sequeces: The Discrete Fourier Trasform 1 3 The 2-D Discrete Fourier Trasform
More informationThe Phi Power Series
The Phi Power Series I did this work i about 0 years while poderig the relatioship betwee the golde mea ad the Madelbrot set. I have fially decided to make it available from my blog at http://semresearch.wordpress.com/.
More information