1 Orthogonalisation in finite precision arithmetic
|
|
- Clinton McCoy
- 5 years ago
- Views:
Transcription
1 1 Orthogonlistion in finite precision rithmetic We investigte the differences nd similrities between the following four wys to compute the QR-decomposition of given rectngulr mtrix A C m n in Mtlb: (CGS) (MGS) (HOU) (GIV) Clssicl Grm-Schmidt Modified Grm-Schmidt Householder reflections Givens rottions It turns out tht in finite precision rithmetic, the results my differ considerbly. 1.1 The QR-decomposition Let n m. A QR-decomposition of A C m n with columns 1,..., n is mtrix fctoristion A = QR (1) where Q is n m n mtrix whose columns q 1,..., q n re orthonorml nd hve the dditionl property tht spn{q 1,..., q l } = spn{ 1,..., l } for ll l {1,..., n}. (2) This implies tht R is n n n upper tringulr mtrix. A QR-decomposition exists for every A C m n. If A hs rnk n, there exists exctly one QR-decomposition of A in which the digonl entries of R re positive. See Theorem 7.2 in Trefethen & Bu. The computer generted results ˆQ nd ˆR of QR-decomposition of mtrix A re generlly not exct. Even though ˆR my be exctly upper tringulr, the entries of ˆR my well differ from those of R. Moreover, the columns of ˆQ re lmost surely not exctly orthogonl. This so-clled loss of orthogonlity results from finite precision rithmetic, the fct tht computers necessrily compute with only finitely mny numbers. 1.2 Tringulr Orthogonliztion: Clssicl versus Modified Grm Schmidt The Clssicl Grm-Schmidt (CGS) lgorithm is chrcterized by the following structure. Jorgen Pedersen Grm ( ) nd Erhrd Schmidt ( ) Assume tht the first l 1 columns q 1,..., q l 1 of the unitry fctor Q together with the (l 1) (l 1) leding principle submtrix of R re lredy known, then the l-th column q l of Q nd the l-th column of R re determined by the reltions r ll q l = l r 1l q 1 r 2l q 2... r l 1,l q l, where r jl = q j l (3) 1
2 for ll j {1,..., l 1}. Thus, j is orthogonlly projected on subspce of dimension l 1, being the spn of q 1,..., q l 1, nd the result is subtrcted from j. Equivlently one my sy tht j is projected onto the orthogonl complement of the spn of q 1,..., q l 1, which hs dimension m l + 1. Remrk: Formul (3) for this projection ssumes tht q 1..., q l 1 form n orthonorml bsis of its spn. Precisely this ssumption will not be fulfilled in finite precision rithmetic. This introduces, prt from rounding errors, lso modeling error. The Modified Grm-Schmidt (MGS) lgorithm is subtle vrint tht is in exct rithmetic equivlent with CGS. However, in finite precision rithmetic it generlly yields more ccurte results. Insted of (3), MGS evlutes r ll q l = (I q l 1 q l 1 )... (I q 2q 2)(I q 1 q 1) l. (4) The mthemticl equivlence with (3) follows from the fct tht if q 1 q 2 =, then (I q 2 q 2)(I q 1 q 1) = I q 1 q 1 q 2 q 2 + q 2 q 2q 1 q 1 = I q 1 q 1 q 2 q 2. This lso explins tht when implemented on computer, in which cse the computed columns of Q re not exctly orthonorml, the methods produce different results. Intuition: In expression (4), l is projected on spce of dimension m 1, the result is lso projected on spce of dimension m 1, nd so on: in totl there re l 1 consecutive projections on subspce of dimension m 1. Some intuition for the success of MGS is the following. Assume tht (I q 1 q1 ) l introduces n error due to finite precision rithmetic. This error is (smll) vector with components tht re rndomly distributed rndom over ll directions q 1,..., q m tht re going to be computed during the process. But m 1 of those m components will be projected to zero in the consecutive steps. Admittedly, pproximtely projected to zero, but this is better thn not doing nything with them t ll, s in (3). Importnt: Expressions of the form (I vv )w should never be implemented s (i) compute the m 2 entries of the m m mtrix I vv nd (ii) multiply with w. Much less expensive is to compute w v(v w), which does not involve ny mtrix Two equivlent orders of computtion in CGS nd MGS It is common prctice in both CGS nd MGS tht the column q l of Q is computed before ny computtions re done in which the columns l+1,..., n of A re involved, becuse these re not needed for the computtion of q l. Just becuse of tht, it is possible to choose nother order of computtion. Indeed, s soon s q 1 is computed, the complete first row of R cn lredy be computed. Moreover, ll columns 2,..., n of A cn be rid of their q 1 -components. For both CGS nd MGS this results in the sme vectors â 2,..., â n. Next, q 2 cn be computed by normlising â 2, nd â 3,..., â n cn be rid of their q 2 components. For CGS this results in vectors â C 3,..., âc n nd for MGS in â M 3,..., âm n, nd so on. Recll tht CGS nd MGS my produce different results only from the third vector onwrds. Remrk: Trefethen & Bu confuse the reder by giving the CGS lgorithm (Algorithm 7.1) in the usul order of computtion, nd MGS (Algorithm 8.1) in the lterntive order. This suggests tht this spect explins the better results of MGS in comprison with CGS, which is flse. The rel gin lies in the difference between (3) nd (4)! 2
3 1.2.2 Exercise () Implement Algorithms 7.1 nd 8.1 in Mtlb s [Q,R]=CGS(A) nd [Q,R]=MGS(A). (b) Give implementtions of CGS nd MGS in which the other order of computtion is used. (c) Verify by exmple tht both versions of CGS produce the sme result. (d) Verify by exmple tht both versions of MGS produce the sme result. Such verifictions re best done with smll exmples in full precision, mening, using Mtlb s formt long. You my wish to store output using the diry commnd of Mtlb Exercise In Mtlb, the commnd A=hilb(n) genertes the Hilbert mtrix of dimensions n n. Use both CGS nd MGS to compute QR-decomposition of the n n Hilbert mtrix for n {4, 5,..., 1}. For ech QR-decomposition, compute using Mtlb the norm I Q Q, which equls zero if nd only if Q hs orthonorml columns. Use the semilogy commnd nd plot both the 7-vectors with norms for CGS nd MGS in logrithmic scle. Use hold to get both grphs in the sme plot. With semilogy(v, -*r ) you cn drw the grph for v in red nd with sterisks, wheres with semilogy(v, -o ) it is in blue (defult) with circles. With legend you cn mke legend nd reposition it using the mouse. See lso grid title xlbel ylbel to pimp the picture even more. Export nd store it s jpg-file. Keep the vribles in memory, the ones from MGS you will use gin in Exercise 1.4.2! 1.3 Tringulr Orthogonliztion versus Orthogonl Tringulriztion As rgued on pge 61 in Trefethen & Bu you cn interpret the orthogonliztion of A using Grm-Schmidt s tringulr orthogonliztion. Indeed, s soon s q 1 is determined, the first row of R cn be computed nd ll columns fo A cn be ridded of their q 1 -components. This corresponds to right-multipliction of A by n upper tringulr mtrix, nd results in n itertive process, A = A, A 1 = A R 1, A 2 = A 1 R 2,..., A n = A n 1 R n, (5) where R j is the mtrix tht removes the q j -component from columns j +1,..., n of the mtrix it cts upon (from the right). Thus, fter n steps, the unitry fctor remins: A n = AR 1 R 2 R n = Q, nd hence, R 1 R 2... R n = R 1. (6) Alterntively, one my try to multiply the mtrix A itertively from the left with unitry mtrix until the upper tringulr fctor remins: nd thus,  = A,  1 = Q 1   2 = Q 2  1...,  n = Q n  n 1, (7) Q n Q 2 Q 1 A = R, nd hence, Q n Q 2 Q 1 = Q. (8) This ltter strtegy will generlly produce better results in finite precision rithmetic. To understnd this, consider the instructive (though not very relistic) sitution in which A 3
4 itself contins errors of some kind, but tht for some mysterious reson in (6) nd (8) ll computtions re performed exctly. Writing the inexct strt mtrix s A + E, the excte mtrix A plus n error mtrix E, then with (6) nd (8) we find tht (A + E)R 1 R 2 R n = AR 1 + ER 1 nd Q n Q 2 Q 1 (A + E) = Q A + Q E. (9) Mesured in unitrily invrint norm (such s 2 nd F ) we know however tht Q E E (1) wheres ER 1 cn in principle be much lrger thn E. Indeed, ER 1 is the solution X of the upper tringulr system XR = E, whose solution cn be much lrger thn the dt due to smll pivots! Remrk: The ssumption tht ll opertions on A + E re performed exctly is unrelistic, nd the negtive effects of tringulr orthogonliztion will in prctice itertively ccumulte! Exercise Let A=hilb(8) nd choose E=.1*rnd(8), mtrix with smll entries compred to those of A. Compute using the Mtlb commnd [Q,R]=qr(A) QR-decomposition of A. Mtlb uses for this the method of the next section. () How orthonorml is this Q in comprison with the one computed in Exercise 1.2.3? Use now the Mtlb commndo [QE,RE]=qr(A+E) to compute QR-decomposition of A + E. Motivted by the result in () we will pretend tht this QR-decomposition is the exct QRdecomposition of A + E, nd thus, tht its fctors re the ones from (9). (b) Compute norm(qe *E)/norm(E) en norm(e*inv(re))/norm(e). Use diry to store your results. Do they confirm the intuition given bove? Add your comments to your diry-file (simply by opening it in text editor). 1.4 Orthogonl Tringulriztion: Householder reflections It is not difficult to find the mtrix of reflection S in hyperplne V, tht mps given x C m on x e 1, where e 1 is the first cnonicl bsis vector of C m. Indeed, the difference v = x x e 1 should be orthogonl to the plne V of reflection, nd thus P (x) = (I vv v v )x = x v (v x) (v v) en S(x) = (I 2 vv v v )x = x 2v (v x) (v v) (11) re respectively equl to the orthogonl projection of x onto V nd the reflection of x in V. See lso Exercise 6.1 in Trefethen & Bu, which introduces reflections. Given ny m n mtrix A with 1 = Ae 1 we cn use (11) to compute the (unitry) mtrix H 1 for which 1... b b 2n H 1 A =.... (12) b m2... b mn 4
5 Similrly, we cn lso determine the mtrix of the reflection tht mps 1 on 1 e 1, for which the hyperplne of reflection is orthogonl to 1 + x 1. This second reflection comes in hndy if 1 en 1 e 1 re lmost equl: s is well known, subtrcting two lmost equl numbers leds to loss of significnt digits, so-clled cncelltion errors. Alston Householder ( ) For this reson, Alston Householder proposed to use the sign of the first entry of x to determine whether to reflect x on x e 1 or on x e 1 : choose v = x + sgn(e 1 x) x e 1 nd reflect in v. Numericlly this is the stble choice. For the purpose of computing QR-decomposition of mtrix A using itertive left-multipliction with unitry mtrices, this mkes no difference. We cn now continue to find the mtix H 2 tht reflects the first column of the (m 1) (n 1) bottom right submtrix of H 1 A onto multiple of the first cnonicl bsis vector, nd so on, to compute QR-decomposition of A. The unitry mtrix Q cn then be explicitly computed s product of ll reflection mtrices. One my lso choose to simply store ll norml vectors to the reflection plnes. It depends on the ctul ppliction of the QR-decomposition if you relly need the fctor Q explicitly. See Algorithms 1.2 nd 1.3 in Trefethen & Bu Exercise Implement QR-decomposition using Householder reflections s [Q,R]=House(A). Compute the unitry mtrix Q by hving ech reflection ct on both A nd I until A hs been tringulrized. Indeed, t tht moment, I hs been trnsformed into Q Exercise Use your code House to compute QR-decomposition of the Hilbert mtrices of size n n for n {4, 5,..., 1}. Agin compute I Q Q. Just s in Exercise 1.2.3, mke grph, but now compre MGS with Householder QR, nd store it s jpg-file. Agin, do not throw wy your computtionl results just yet. 1.5 Orthogonl Tringulriztion: Givens rottions Aprt from reflections, lso rottions cn be used s unitry trnsformtions with which mtrix A cn be tringulrized. Given nonzero vector in R 2 it is esy to write down the orthogonl mtrix of the rottion bout the origin tht mps this vector on multiple of e 1, R(, b) [ b ] [ ] = 2 + b 2 where R(, b) = [ b 2 b b ]. (13) 5
6 Note: From the bove we cn lredy conclude tht [ ] [ A = = R(, b) 2 + b 2 b ] 1 = UΣV (14) which is singulr vlue decomposition of the 2 1 mtrix ( b)! A vector x = (x 1, x 2, x 3, x 4 ) R 4 cn be isometriclly mpped onto multiple of e 1 by product of three rottions. The first is in the (x 3, x 4 )-plne, the second in the (x 2, x 3 )-plne, nd the third in the (x 1, x 2 )-plne. Schemticlly, x 1 x 2 x 3 x 4 R 34(x 3,x 4) x 1 x 2 y 3 R 23(x 2,y 3) x 1 y 2 R 12(x 1,y 2) y 1. (15) Next, we show how to employ such rottions to bring 4 3 mtrix in upper tringulr form. This we do column-wise s follows. Apply the three rottions mpping the first column of X onto multiple of e 1 ; X = R 34(, ) R 23(, ) R 12(, ) = Y. The product Y = R 12 R 23 R 34 X hs three zeros below (1, 1). Note in prticulr tht the second nd third column of Y re generlly not the sme s those of X! Apply two rottions to crete zeros on the two positions below (2, 2) vn Y ; Y = R 34(, ) R 23(, ) = Z. Observe tht this cn be done without destroying previously introduced zero entries. Apply the rottion tht introduces nother zero entry, t position (3, 4); Z = R 34(, ) = U. The finl result U, fter pplying six rottions, is upper tringulr. Wllce Givens ( ) 6
7 A plne rottion in the context of computing QR-decomposition is nmed fter Wllce Givens. In de context of similrity trnsforms, they re clled fter Crl Jcobi Exercise Investigve the commnd givens nd use it to write code [Q,R]=PlneRots(A) tht computes QR-decomposition of A using Givens rottions. In this implementtion, compute Q by pplying the rottions to both A nd I until A is tringulrized Exercise Use your code PlneRots to compute QR-decomposition of the Hilbert mtrices of size n n for n {4, 5,..., 1}. Compute I Q Q. Just s in Exercises nd 1.4.2, mke grph compring Householder QR nd Givens QR. 1.6 Hnding in Hnd in the following the m-files of both vrints of CGS; the m-files vn both vrints of MGS; the m-file for House; the m-file for PlneRots; the jpg s of your figures for Exercises nd nd the diry prts for Exercises nd 1.3.1, in one file, plese. Thus, ten files in totl. Plese send them s seprte files in one e-mil (i.e. without zipping) to Jn: J.H.Westerdiep@uv.nl Dedline for this ssignment is Sundy October 1, 23:59 hrs. Wrning: This ssignment not only counts towrds your finl grde, it lso decides bout pssing for the course. Therefore, ech student hs to hnd in his or her work by him/herself. Copying (or cosmeticlly dpting with the purpose of hiding plgirism) codes or prts of codes of fellow students is strictly prohibited nd my led to hering by the FNWI Plgirism Bord. However, coopertion with fellow students on globl level is llowed nd encourged. In tht cse it is dviced (nd good prctice) to explicitly stte in your work with whom you hve cooperted, s to explin possible superficil similrity between your work. 7
Numerical Linear Algebra Assignment 008
Numericl Liner Algebr Assignment 008 Nguyen Qun B Hong Students t Fculty of Mth nd Computer Science, Ho Chi Minh University of Science, Vietnm emil. nguyenqunbhong@gmil.com blog. http://hongnguyenqunb.wordpress.com
More informationMath 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of
More informationHere we study square linear systems and properties of their coefficient matrices as they relate to the solution set of the linear system.
Section 24 Nonsingulr Liner Systems Here we study squre liner systems nd properties of their coefficient mtrices s they relte to the solution set of the liner system Let A be n n Then we know from previous
More informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk bout solving systems of liner equtions. These re problems tht give couple of equtions with couple of unknowns, like: 6 2 3 7 4
More informationLecture Note 9: Orthogonal Reduction
MATH : Computtionl Methods of Liner Algebr 1 The Row Echelon Form Lecture Note 9: Orthogonl Reduction Our trget is to solve the norml eution: Xinyi Zeng Deprtment of Mthemticl Sciences, UTEP A t Ax = A
More information7.2 The Definite Integral
7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
More informationHW3, Math 307. CSUF. Spring 2007.
HW, Mth 7. CSUF. Spring 7. Nsser M. Abbsi Spring 7 Compiled on November 5, 8 t 8:8m public Contents Section.6, problem Section.6, problem Section.6, problem 5 Section.6, problem 7 6 5 Section.6, problem
More informationMATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1
MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further
More informationSUMMER KNOWHOW STUDY AND LEARNING CENTRE
SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18
More information1 The Lagrange interpolation formula
Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x
More informationImproper Integrals, and Differential Equations
Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted
More informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk out solving systems of liner equtions. These re prolems tht give couple of equtions with couple of unknowns, like: 6= x + x 7=
More informationHandout: Natural deduction for first order logic
MATH 457 Introduction to Mthemticl Logic Spring 2016 Dr Json Rute Hndout: Nturl deduction for first order logic We will extend our nturl deduction rules for sententil logic to first order logic These notes
More information1 Linear Least Squares
Lest Squres Pge 1 1 Liner Lest Squres I will try to be consistent in nottion, with n being the number of dt points, nd m < n being the number of prmeters in model function. We re interested in solving
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More informationChapter 0. What is the Lebesgue integral about?
Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous
More informationChapter 14. Matrix Representations of Linear Transformations
Chpter 4 Mtrix Representtions of Liner Trnsformtions When considering the Het Stte Evolution, we found tht we could describe this process using multipliction by mtrix. This ws nice becuse computers cn
More information20 MATHEMATICS POLYNOMIALS
0 MATHEMATICS POLYNOMIALS.1 Introduction In Clss IX, you hve studied polynomils in one vrible nd their degrees. Recll tht if p(x) is polynomil in x, the highest power of x in p(x) is clled the degree of
More informationMATRICES AND VECTORS SPACE
MATRICES AND VECTORS SPACE MATRICES AND MATRIX OPERATIONS SYSTEM OF LINEAR EQUATIONS DETERMINANTS VECTORS IN -SPACE AND -SPACE GENERAL VECTOR SPACES INNER PRODUCT SPACES EIGENVALUES, EIGENVECTORS LINEAR
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More information1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),
1. Guss-Jcobi qudrture nd Legendre polynomils Simpson s rule for evluting n integrl f(t)dt gives the correct nswer with error of bout O(n 4 ) (with constnt tht depends on f, in prticulr, it depends on
More informationIntroduction To Matrices MCV 4UI Assignment #1
Introduction To Mtrices MCV UI Assignment # INTRODUCTION: A mtrix plurl: mtrices) is rectngulr rry of numbers rrnged in rows nd columns Exmples: ) b) c) [ ] d) Ech number ppering in the rry is sid to be
More informationBest Approximation in the 2-norm
Jim Lmbers MAT 77 Fll Semester 1-11 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the -norm Suppose tht we wish to obtin function f n (x) tht is liner combintion
More informationThings to Memorize: A Partial List. January 27, 2017
Things to Memorize: A Prtil List Jnury 27, 2017 Chpter 2 Vectors - Bsic Fcts A vector hs mgnitude (lso clled size/length/norm) nd direction. It does not hve fixed position, so the sme vector cn e moved
More informationRecitation 3: More Applications of the Derivative
Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech
More informationOrthogonal Polynomials and Least-Squares Approximations to Functions
Chpter Orthogonl Polynomils nd Lest-Squres Approximtions to Functions **4/5/3 ET. Discrete Lest-Squres Approximtions Given set of dt points (x,y ), (x,y ),..., (x m,y m ), norml nd useful prctice in mny
More informationChapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY
Chpter 3 MTRIX In this chpter: Definition nd terms Specil Mtrices Mtrix Opertion: Trnspose, Equlity, Sum, Difference, Sclr Multipliction, Mtrix Multipliction, Determinnt, Inverse ppliction of Mtrix in
More informationWeek 10: Line Integrals
Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.
More informationModule 6: LINEAR TRANSFORMATIONS
Module 6: LINEAR TRANSFORMATIONS. Trnsformtions nd mtrices Trnsformtions re generliztions of functions. A vector x in some set S n is mpped into m nother vector y T( x). A trnsformtion is liner if, for
More informationImproper Integrals. Type I Improper Integrals How do we evaluate an integral such as
Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More information( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that
Arc Length of Curves in Three Dimensionl Spce If the vector function r(t) f(t) i + g(t) j + h(t) k trces out the curve C s t vries, we cn mesure distnces long C using formul nerly identicl to one tht we
More informationBefore we can begin Ch. 3 on Radicals, we need to be familiar with perfect squares, cubes, etc. Try and do as many as you can without a calculator!!!
Nme: Algebr II Honors Pre-Chpter Homework Before we cn begin Ch on Rdicls, we need to be fmilir with perfect squres, cubes, etc Try nd do s mny s you cn without clcultor!!! n The nth root of n n Be ble
More informationMath 426: Probability Final Exam Practice
Mth 46: Probbility Finl Exm Prctice. Computtionl problems 4. Let T k (n) denote the number of prtitions of the set {,..., n} into k nonempty subsets, where k n. Argue tht T k (n) kt k (n ) + T k (n ) by
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationChapters 4 & 5 Integrals & Applications
Contents Chpters 4 & 5 Integrls & Applictions Motivtion to Chpters 4 & 5 2 Chpter 4 3 Ares nd Distnces 3. VIDEO - Ares Under Functions............................................ 3.2 VIDEO - Applictions
More informationGoals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite
Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite
More informationMatrices, Moments and Quadrature, cont d
Jim Lmbers MAT 285 Summer Session 2015-16 Lecture 2 Notes Mtrices, Moments nd Qudrture, cont d We hve described how Jcobi mtrices cn be used to compute nodes nd weights for Gussin qudrture rules for generl
More informationBernoulli Numbers Jeff Morton
Bernoulli Numbers Jeff Morton. We re interested in the opertor e t k d k t k, which is to sy k tk. Applying this to some function f E to get e t f d k k tk d k f f + d k k tk dk f, we note tht since f
More informationReview of basic calculus
Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below
More informationScientific notation is a way of expressing really big numbers or really small numbers.
Scientific Nottion (Stndrd form) Scientific nottion is wy of expressing relly big numbers or relly smll numbers. It is most often used in scientific clcultions where the nlysis must be very precise. Scientific
More informationThe Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.
Clculus Li Vs The Fundmentl Theorem of Clculus. The Totl Chnge Theorem nd the Are Under Curve. Recll the following fct from Clculus course. If continuous function f(x) represents the rte of chnge of F
More informationChapter 3 Polynomials
Dr M DRAIEF As described in the introduction of Chpter 1, pplictions of solving liner equtions rise in number of different settings In prticulr, we will in this chpter focus on the problem of modelling
More informationIntroduction to Determinants. Remarks. Remarks. The determinant applies in the case of square matrices
Introduction to Determinnts Remrks The determinnt pplies in the cse of squre mtrices squre mtrix is nonsingulr if nd only if its determinnt not zero, hence the term determinnt Nonsingulr mtrices re sometimes
More informationMATH 174A: PROBLEM SET 5. Suggested Solution
MATH 174A: PROBLEM SET 5 Suggested Solution Problem 1. Suppose tht I [, b] is n intervl. Let f 1 b f() d for f C(I; R) (i.e. f is continuous rel-vlued function on I), nd let L 1 (I) denote the completion
More informationREPRESENTATION THEORY OF PSL 2 (q)
REPRESENTATION THEORY OF PSL (q) YAQIAO LI Following re notes from book [1]. The im is to show the qusirndomness of PSL (q), i.e., the group hs no low dimensionl representtion. 1. Representtion Theory
More informationa < a+ x < a+2 x < < a+n x = b, n A i n f(x i ) x. i=1 i=1
Mth 33 Volume Stewrt 5.2 Geometry of integrls. In this section, we will lern how to compute volumes using integrls defined by slice nlysis. First, we recll from Clculus I how to compute res. Given the
More informationAbstract inner product spaces
WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the
More informationLinear Algebra 1A - solutions of ex.4
Liner Algebr A - solutions of ex.4 For ech of the following, nd the inverse mtrix (mtritz hofkhit if it exists - ( 6 6 A, B (, C 3, D, 4 4 ( E i, F (inverse over C for F. i Also, pick n invertible mtrix
More informationDuality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.
Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we
More informationStudent Activity 3: Single Factor ANOVA
MATH 40 Student Activity 3: Single Fctor ANOVA Some Bsic Concepts In designed experiment, two or more tretments, or combintions of tretments, is pplied to experimentl units The number of tretments, whether
More informationChapter 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More informationFormal languages, automata, and theory of computation
Mälrdlen University TEN1 DVA337 2015 School of Innovtion, Design nd Engineering Forml lnguges, utomt, nd theory of computtion Thursdy, Novemer 5, 14:10-18:30 Techer: Dniel Hedin, phone 021-107052 The exm
More informationConvert the NFA into DFA
Convert the NF into F For ech NF we cn find F ccepting the sme lnguge. The numer of sttes of the F could e exponentil in the numer of sttes of the NF, ut in prctice this worst cse occurs rrely. lgorithm:
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More informationAPPROXIMATE INTEGRATION
APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be
More information5.7 Improper Integrals
458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the
More informationContents. Outline. Structured Rank Matrices Lecture 2: The theorem Proofs Examples related to structured ranks References. Structure Transport
Contents Structured Rnk Mtrices Lecture 2: Mrc Vn Brel nd Rf Vndebril Dept. of Computer Science, K.U.Leuven, Belgium Chemnitz, Germny, 26-30 September 2011 1 Exmples relted to structured rnks 2 2 / 26
More informationState space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies
Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response
More informationInfinite Geometric Series
Infinite Geometric Series Finite Geometric Series ( finite SUM) Let 0 < r < 1, nd let n be positive integer. Consider the finite sum It turns out there is simple lgebric expression tht is equivlent to
More informationLecture 1. Functional series. Pointwise and uniform convergence.
1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is
More informationN 0 completions on partial matrices
N 0 completions on prtil mtrices C. Jordán C. Mendes Arújo Jun R. Torregros Instituto de Mtemátic Multidisciplinr / Centro de Mtemátic Universidd Politécnic de Vlenci / Universidde do Minho Cmino de Ver
More information13.3 CLASSICAL STRAIGHTEDGE AND COMPASS CONSTRUCTIONS
33 CLASSICAL STRAIGHTEDGE AND COMPASS CONSTRUCTIONS As simple ppliction of the results we hve obtined on lgebric extensions, nd in prticulr on the multiplictivity of extension degrees, we cn nswer (in
More informationACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019
ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS MATH00030 SEMESTER 208/209 DR. ANTHONY BROWN 7.. Introduction to Integrtion. 7. Integrl Clculus As ws the cse with the chpter on differentil
More informationParse trees, ambiguity, and Chomsky normal form
Prse trees, miguity, nd Chomsky norml form In this lecture we will discuss few importnt notions connected with contextfree grmmrs, including prse trees, miguity, nd specil form for context-free grmmrs
More informationReview of Gaussian Quadrature method
Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge
More informationChapter 3. Vector Spaces
3.4 Liner Trnsformtions 1 Chpter 3. Vector Spces 3.4 Liner Trnsformtions Note. We hve lredy studied liner trnsformtions from R n into R m. Now we look t liner trnsformtions from one generl vector spce
More informationMAA 4212 Improper Integrals
Notes by Dvid Groisser, Copyright c 1995; revised 2002, 2009, 2014 MAA 4212 Improper Integrls The Riemnn integrl, while perfectly well-defined, is too restrictive for mny purposes; there re functions which
More informationp-adic Egyptian Fractions
p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction
More informationMTH 5102 Linear Algebra Practice Exam 1 - Solutions Feb. 9, 2016
Nme (Lst nme, First nme): MTH 502 Liner Algebr Prctice Exm - Solutions Feb 9, 206 Exm Instructions: You hve hour & 0 minutes to complete the exm There re totl of 6 problems You must show your work Prtil
More informationCHAPTER 4 MULTIPLE INTEGRALS
CHAPTE 4 MULTIPLE INTEGAL The objects of this chpter re five-fold. They re: (1 Discuss when sclr-vlued functions f cn be integrted over closed rectngulr boxes in n ; simply put, f is integrble over iff
More informationf(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral
Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one
More informationMath 270A: Numerical Linear Algebra
Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner
More informationBases for Vector Spaces
Bses for Vector Spces 2-26-25 A set is independent if, roughly speking, there is no redundncy in the set: You cn t uild ny vector in the set s liner comintion of the others A set spns if you cn uild everything
More informationNumerical Integration. 1 Introduction. 2 Midpoint Rule, Trapezoid Rule, Simpson Rule. AMSC/CMSC 460/466 T. von Petersdorff 1
AMSC/CMSC 46/466 T. von Petersdorff 1 umericl Integrtion 1 Introduction We wnt to pproximte the integrl I := f xdx where we re given, b nd the function f s subroutine. We evlute f t points x 1,...,x n
More informationLesson 1.6 Exercises, pages 68 73
Lesson.6 Exercises, pges 68 7 A. Determine whether ech infinite geometric series hs finite sum. How do you know? ) + +.5 + 6.75 +... r is:.5, so the sum is not finite. b) 0.5 0.05 0.005 0.0005... r is:
More informationWe partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.
Mth 255 - Vector lculus II Notes 4.2 Pth nd Line Integrls We begin with discussion of pth integrls (the book clls them sclr line integrls). We will do this for function of two vribles, but these ides cn
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More informationElementary Linear Algebra
Elementry Liner Algebr Anton & Rorres, 1 th Edition Lecture Set 5 Chpter 4: Prt II Generl Vector Spces 163 คณ ตศาสตร ว ศวกรรม 3 สาขาว ชาว ศวกรรมคอมพ วเตอร ป การศ กษา 1/2555 163 คณตศาสตรวศวกรรม 3 สาขาวชาวศวกรรมคอมพวเตอร
More informationSelf-similarity and symmetries of Pascal s triangles and simplices mod p
Sn Jose Stte University SJSU ScholrWorks Fculty Publictions Mthemtics nd Sttistics Februry 2004 Self-similrity nd symmetries of Pscl s tringles nd simplices mod p Richrd P. Kubelk Sn Jose Stte University,
More informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution
More informationLecture 2e Orthogonal Complement (pages )
Lecture 2e Orthogonl Complement (pges -) We hve now seen tht n orthonorml sis is nice wy to descrie suspce, ut knowing tht we wnt n orthonorml sis doesn t mke one fll into our lp. In theory, the process
More information308K. 1 Section 3.2. Zelaya Eufemia. 1. Example 1: Multiplication of Matrices: X Y Z R S R S X Y Z. By associativity we have to choices:
8K Zely Eufemi Section 2 Exmple : Multipliction of Mtrices: X Y Z T c e d f 2 R S X Y Z 2 c e d f 2 R S 2 By ssocitivity we hve to choices: OR: X Y Z R S cr ds er fs X cy ez X dy fz 2 R S 2 Suggestion
More informationDiscrete Least-squares Approximations
Discrete Lest-squres Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve
More informationset is not closed under matrix [ multiplication, ] and does not form a group.
Prolem 2.3: Which of the following collections of 2 2 mtrices with rel entries form groups under [ mtrix ] multipliction? i) Those of the form for which c d 2 Answer: The set of such mtrices is not closed
More informationECO 317 Economics of Uncertainty Fall Term 2007 Notes for lectures 4. Stochastic Dominance
Generl structure ECO 37 Economics of Uncertinty Fll Term 007 Notes for lectures 4. Stochstic Dominnce Here we suppose tht the consequences re welth mounts denoted by W, which cn tke on ny vlue between
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationLine and Surface Integrals: An Intuitive Understanding
Line nd Surfce Integrls: An Intuitive Understnding Joseph Breen Introduction Multivrible clculus is ll bout bstrcting the ides of differentition nd integrtion from the fmilir single vrible cse to tht of
More informationAPPLICATIONS OF THE DEFINITE INTEGRAL
APPLICATIONS OF THE DEFINITE INTEGRAL. Volume: Slicing, disks nd wshers.. Volumes by Slicing. Suppose solid object hs boundries extending from x =, to x = b, nd tht its cross-section in plne pssing through
More information13: Diffusion in 2 Energy Groups
3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups
More informationfractions Let s Learn to
5 simple lgebric frctions corne lens pupil retin Norml vision light focused on the retin concve lens Shortsightedness (myopi) light focused in front of the retin Corrected myopi light focused on the retin
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More informationRiemann Integrals and the Fundamental Theorem of Calculus
Riemnn Integrls nd the Fundmentl Theorem of Clculus Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University September 16, 2013 Outline Grphing Riemnn Sums
More informationMatrix Solution to Linear Equations and Markov Chains
Trding Systems nd Methods, Fifth Edition By Perry J. Kufmn Copyright 2005, 2013 by Perry J. Kufmn APPENDIX 2 Mtrix Solution to Liner Equtions nd Mrkov Chins DIRECT SOLUTION AND CONVERGENCE METHOD Before
More informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More informationEE263 homework 8 solutions
EE263 Prof S Boyd EE263 homework 8 solutions 37 FIR filter with smll feedbck Consider cscde of 00 one-smple delys: u z z y () Express this s liner dynmicl system x(t + ) = Ax(t) + Bu(t), y(t) = Cx(t) +
More informationNumerical Methods I. Olof Widlund Transcribed by Ian Tobasco
Numericl Methods I Olof Widlund Trnscribed by In Tobsco Abstrct. This is prt one of two semester course on numericl methods. The course ws offered in Fll 011 t the Cournt Institute for Mthemticl Sciences,
More information