TMA 4180 Optimeringsteori Least Squares Optimization. Harald E. Krogstad, rev. 2010
|
|
- Rhoda Reynolds
- 5 years ago
- Views:
Transcription
1 TMA 48 Optimeringsteori Lest Squres Optimiztion Hrld E. Krogstd, rev. This note ws originlly prepred for erlier versions of the course. Nocedl nd Wright hs nice introduction to Lest Squre (LS) optimiztion problems in Chpter, nd the note is now therefore only smll supplement. It re ects tht LS problems re by fr the most common cse for unconstrined optimiztion. See lso N&W, p for typicl exmples of LS problems. Lest Squre problems hve often their origin in tting models to observtions. In its simplest form, we know this from the problem of tting regression line, y = x + b, through set of dt points fx i ; y i g ; i = ; ; N. When N >, it is in generl impossible to put the line through ll points, but we try to determine n optiml line, for exmple by determining the pir f ; b g which minimizes the objective function f(; b) = NX (x i + b y i ) : () i= In this simple cse, it is esy to derive the solution nlyticlly, but in generl the solution hs to be found numericlly. Since these problems re so common, there hs been lot of work involved in dpting the generl lgorithms for unconstrined optimiztion to this specil cse. It is more e ective to use specilly dpted lgorithms insted of the more generl ones. Below we rst consider the liner LS problem, which some of you would know from liner nlysis or multiple liner regression in sttistics. The nottion my di er somewht from wht you know, but the ides re the sme. The singulr vlue decomposition is n importnt tool for nlyzing the liner problem, nd the liner problem is lso importnt for the solutions of the vrious itertion steps in non-liner LS problems. In the non-liner lgorithms we once more meet bsic tools like line serch nd trust region methods, which utilize the specil form of the grdient nd Hessin the LS problems hve. Liner Lest Squre Problems. The origin Although LS problems stem from vriety of di erent situtions, we shll consider sitution where n fctors in uence the result of n experiment. For ech experiment we know the vlues of the fctors, sy f ; :::; n g, nd we re ble to mesure the result. We re interested in nding function f ( ; :::; n ) which cn be used to compute the result without hving to repet the experiment ech time. In prctice, this my be quite di cult: A lot of fctors my be involved
2 We re not quite sure whether fctor counts or not There my be fctors unknown to us We hve no ide how the function f looks It is expensive nd/or time consuming to crry out experiments We hve only limited time nd budget to gure out solution All these points re prt of rel life, where you often nd yourself in group of collegues where you re the lest incompetent member. The rst guess if one does not know nything will often be to ssume tht there is liner dependency between the fctors nd the result. In other words, we re looking for set of coe cients or weight fctors fx ; ; x n g such tht x + x + ::: + n x n : () Note tht things re turned round here: When we pply the formul fter the x-s re determined, we put in f ; ; n g nd get n estimte of, denoted ^. The expression should in generl lso include constnt term, nd this is obtined by letting =, so tht the formul reds ^ = x + x + + n x n. In order to determine the x-s, we need to crry out set of clibrtion experiments. If we do not know nything bout the x-s, we would like to hve m > n experiments: If the number of experiments is less thn n (often the sitution in prctice), there will in generl be mny di erent sets of weight fctors tting Eqn., nd if we hve exctly n experiments, we could still stisfy Eqn. exctly, but would hve no ide how well the eqution is stis ed in other cses. Let us store the vlues of the fctors in the mn mtrix A nd the results in the m dimensionl vector b, such tht A = f i ; i ; ; in g; i = ; ; m; () b = f ; ; m g : (4) Since it would be too good to be true if we relly found vector x such tht Ax = b, we re stis ed if we re ble to nd x such tht the sum of the squred devitions is minimized: min x mx ( i x + i x + ::: + in x n i ) : (5) i= In mtrix nottion this is equivlent to nding the minimum of the function f(x) = kax bk = (Ax b) (Ax b): (6)
3 . Solution of the full rnk problem By crrying out the multiplictions in Eqn. 6, we nd tht f(x) = (Ax b) (Ax b) = x A Ax x A b + b b ; (7) which we, by now, should recognize s the qudrtic test problem. If A A hppens to be positive de nite (which requires tht m n), we lredy know tht the solution, x, stis es the eqution A Ax = A b: (8) (Convince yourself tht ll mtrix multiplictions mke sense!). The system of equtions in 8 is clled the Norml Equtions. It turns out tht A A is positive de nite if nd only if A hs full column rnk. In this cse, where m > n, this mens tht rnk(a) = n, or tht ll columns of A re linerly independent. Tht is equivlent to Ax =, x =, which in turn is equivlent to x A Ax = kaxk =, x =. Even if A hs full rnk mthemticlly, it is not necessrily resonble numericlly to solve the norml equtions. The condition number of A is the rtio between its lrgest nd smllest singulr vlue (see below), nd the condition number of A A will be the squre of tht. In order to void numericl problems, the best is to use QR-fctoriztion of A. The -norm is invrint under orthogonl trnsformtions, nd this reduces the problem to solving liner system of equtions with n unknowns which hs the upper n n sub-mtrix of R s its coe cient mtrix. A short summry of this is given in N&W p. 5, nd more extensive tretment is found in the book of Golub nd vnlon [], Chpter 5. Cses with few fctors nd lrge number of experiments re not likely to experience ny problems with the Norml Equtions.. The rnk de cient problem If the rnk of A, r = rnk (A), is less thn n, the problem will still hve solutions, but they re not unique. The null spce of A, N (A), is the set of ll vectors z such tht Az =, nd when the column rnk, rnk (A), is less thn n, N (A) contins non-zero vectors (recll the de nition of rnk!). Obviously, if x is solution to 6, then x + z, z N (A) will be solution s well since A Az = nd f (x + z) = f (x ). The existence of solutions is not obvious, but follows by considering Eqn. 8 s n eqution not in R n but in the smller orthogonl complement of the null-spce of A, N (A)?. Check out ll de nitions you do not know/recll nd try to crry out the following steps:. Show tht A b nd A Ax re members of N (A)? by tking the sclr product with vector z N (A).
4 . Show tht the mtrix A A is non-singulr when restricted to N (A)? (This mounts to show tht kayk = () y = when y N (A)? ).. Show tht there is unique x N (A)? such tht A Ax = A b (This my be little tough if you re not fmilir with liner lgebr nd liner trnsformtions. It is not tht importnt, the solution exists). The rnk de cient cse will occur if we try to explin our result by including too mny fctors. In prctice, it is smrt to strt with few fctors nd then include one dditionl fctor t time. By supervising kax bk, we get feeling of when enough is enough. Even if A is rnk de cient, the function f (x) in Eqn. 7 is still convex since A A, nd the set of globl minim, tht there re indeed optiml solutions, so = fy; f(y) = ming, is therefore closed convex set (We proved bove is de nitely non-empty). Are ll these possible solutions eqully good? Mthemticlly yes, but in prctice we re using x for predicting or estimting, nd the predictor ^ for given set of fctors,, is simply ^ = x : (9) Unnecessry lrge components in x will mgnify ny mesuring errors in the fctors. Therefore, it is resonble to use s smll x s possible. If we, on the convex set, de ne the strictly convex function g(y) = kyk, the function will hve unique minimum, nd this would be the specil solution we re looking for: x = rgmin kxk : () A t Ax=A t b Cn you prove tht x hs to be equl to the solution in N (A)?? (Recll tht ny solution my be written uniquely s x = x + z, where z N (A) nd x N (A)? ). There re severl numericl methods for nding x described in N&W, s well s in []. The most robust lthough not the chepest computtionlly is bsed on the Singulr Vlue Decomposition of A... The Singulr Vlue Decomposition The Singulr Vlue Decomposition (SVD) is door-opener to lot of prcticl lgorithms, nd very much used in modern pplictions of liner lgebr. The reduced form SVD of generl m n mtrix A is A = UV = rx i u i vi; () i= where U = 4 j j j u u ::: u r j j j 5 ; V = 4 j j j v v ::: v r j j j 5 () 4
5 hve orthogonl columns, U U = I (r) ; V V = I (r) ; () nd r = rnk (A). The r r digonl mtrix contins the so-clled singulr vlues of A. These singulr vlues re the squre root of the non-negtive eigenvlues of A A, or AA (Note tht the mtrices U nd V my be extended to full m m nd n n orthogonl mtrices nd to n m n mtrix, s shown in N&W, p We do not ssume or need this below). There re few things to check here: The mtrix A represents liner trnsformtion from R n! R m nd, similrly, A liner trnsformtion from R m! R n. v ; v ; :::; v r R n nd u ; u ; :::; u r R m. U nd V consist of the eigenvectors of AA nd A A tht correspond to non-zero eigenvlues, (AA ) U = U nd (A A) V = V. The columns of V spn n orthogonl bsis for N (A)? R n, nd the vectors in U n orthogonl bsis for the rnge of A, R (A) R m. The mtrix A + = V U (4) is clled pseudo-inverse or generlized inverse of A, nd is the most fmous member mong lrge number of generlized inverses. This prticulr generlized inverse is clled the Moore- Penrose inverse of A. Exercise: Prove tht A + A cts s the identity mtrix on N (A)?, nd AA + s the identity mtrix on R (A) (Hint: Use tht ll vectors x N (A)? my be written x = V s, s R r, nd similrly z = Ut; t R r for z R (A) ). You could consult textbook (E.g. [] or your own liner lgebr book) in order to brush this up. It is not so esy to grsp, or to remember!.. The solution in terms of the SVD Given the SVD of A, the unique solution with the smllest norm of Eqn. 8 is x = A + b = V U b = rx i= u i b i v i : (5) (Prove this yourself by showing tht A A (A + b) = A b!). In Eqn. 5, we ssume tht the singulr vlues re ordered ccording to their size, tht is, r : (6) 5
6 The dt set Figure : Simulted 6 dt points with some underlying function dded to lot of noise. This is the opposite ordering from wht we hve been using for eigenvlues. Contrry to wht you lern in the liner lgebr courses, determining the rnk of big mtrix is fr from trivil. If the singulr vlues decrese grdully to, it is quite di cult to sy where to stop (perhps when i < 6?) nd de ne the rest to be. The nice thing bout Eqn. 5 is tht we cn dd one term t time nd stop when the solution strts to get unstble. Smll singulr vlues mgnify errors in b, U or V, nd this cse is clled n ill-conditioned problem. It seems to be n inherent property of rel life problems of some size tht they ll tend to be ill-conditioned. This will be discussed lter for Inverse Problems... A worked exmple The following simple nd somewht unrelistic exmple shows how the SVD works. We hve lrge nd very noisy dt set of function of one vrible (y). The vlues of y rnge between nd, nd the ide is to extrct the function by tting polynomil. Typicl dt re shown in Fig., nd if this hs been rel dt set, we would not even think of tting something else thn constnt vlue, or t most stright line. The dt re produced by dding independent normlly distributed noise with stndrd devition equl to.6 to the function vlues. Let us nevertheless tke very bold ttitude nd try to t polynomil of order 6 (!) by just running the Mtlb functions poly t nd polyvl: p = polyfit(y,b,6); Y = polyvl(p,y); plot(y,dt,.,y,y); In this cse, the fctors re the vrious powers of vrible y, = ; y; y ; ; y 6 weights x re the coe cients in the polynomil, nd the ^ = x + x y + x y + + x 6 y 6 : (7) 6
7 Number of terms = Figure : Dt nd 6th order polynomil tted by poly t in Mtlb. As expected, the result is relly not tht gret, s Fig. shows. Also Mtlb protests nd wrns us tht the mtrices involved re close to singulr. The progrm suggests removing some dt points, or do some centering nd scling of the dt. Nothing of this would ctully help, the obvious cure is clerly to reduce the order of the polynomil. But let us not give up tht fst. We stick to 6th order nd form the fctor mtrix A of dimension 6 6 (6 fctors nd 6 experiments). A = 6 4. nd the dt vector b = (b ; b ; ; b 6 ). y y y 6 y y y 6.. y 6 y6 y6 7 5 ; (8) We re then fced with the problem of nding the best lods x = (x ; x ; ; x 6 ) to the dt, tht is, solve the LS problem, Let us consider two possibilities: x = rg min xr 6 kax bk : (9) Use of the SVD Use of the norml equtions Both re strightforwrd to progrm, the SVD even s simple s [U Sigm V] = svd(a)! The rnk (A) computed by Mtlb is equl to 6, nd this seems to be in rough greement with plot of the singulr vlues shown on Fig.. It is tempting to use the solution on the SVD-form, x = V U b = rx i= u i b i v i : () 7
8 5 5 σ k k Figure : The singulr vlues of the A-mtrix. The mchine ccurcy seems to be reched fter bout singulr vlues, but the singulr vlues never become exctly. where we include one term t time. This is displyed in Fig. 4, where the result is plotted every eight term. The grphs in the middle show some similrities, nd this is even true for the grphs t the bottom, lthough using more terms in Eqn. thn the rnk is not very meningful. The SVD solution with r = rnk (A) terms, the Norml eqution solution, nd the (up to now) unknown true function re shown in Fig. 5. The nl plot, Fig. 6 shows the SVD solution nd the t of polynomil of order r. The solutions re bout eqully good, but we hve not tried to optimize this further. By running the progrm with di erent dt, the optiml SVD solution seems, on the verge, to be slightly better thn the Norml Eqution solution. The reder should be wrned tht the, fter ll, civilized behviour of these results is in prt due to Mtlb s very crefully coded numerics! Non-Liner Lest Squre Problems In the generl non-liner lest squre cse we hve vector-vlued function, h(x) R m where x R n, nd we wnt to nd the minimum of f(x) = kh(x)k = mx h i (x) = h(x) h(x): () i= The simple Exmple. in N&W on p is of this form. All the grdients of the vrious components in h mke up wht is clled the Jcobin mtrix of h, de ned by J(x) (x) j i=;:::;m; j=;:::;n 6 4 rh (x) rh (x). rh m (x) 7 5 : () Recll tht we ssume tht rh i re n-dimensionl row vectors. By tking derivtives term by 8
9 Number of terms =8 Number of terms =6 Dt nd polynomil Dt nd polynomil Number of terms =4 Number of terms = Dt nd polynomil Dt nd polynomil Number of terms =4 Number of terms =48 Dt nd polynomil Dt nd polynomil Figure 4: SVD solutions truncted t the term indicted. All polynomils hve order 6 (The dt set here di ers from the bove). 9
10 Dt nd ll solutions Figure 5: Dt nd ll polynomils: Red: SVD solution with r terms; Blue: Norml Equtions; Green: Underlying exct function..5.5 Solutions Figure 6: Red: SVD solution of order r = rnk (A). Blue: A simple polynomil t of order r. Green: Input function.
11 term, we esily obtin (do it yourself!) tht the grdient nd the Hessin of f re rf(x) = J (x)h(x); () mx r f (x) = J (x)j(x) + h i (x)r h i (x); (4) where r h i re the Hessin mtrices for h i, i = ; ; m. Note tht if h i (x) re liner functions, sy h(x) = Ax b, then f(x) = kax bk, nd i= rf(x) = A (Ax b) = A h(x); r f(x) = A A: (5) We observe tht A corresponds to the Jcobin, J, wheres the Hessin in the generl non-liner cse gets n dditionl term B = mx h i (x)r h i (x): (6) i= If the minimiztion leds to smll vlues of h i (x), or the problem is nerly liner such tht r h i (x)-s re relly smll, the rst term in r f(x) will dominte nd This is utilized in Guss-Newtons method. method could be written or, in terms of the serch direction, p k, r f(x) J(x) J(x): (7) We remember tht the itertion step in Newton s r f(x k )(x k+ x k ) = rf(x k ) ; (8) r f(x k )p k = rf(x k ) : (9) In the Guss-Newtons method, the serch direction p k for line serch, x k+ = x k + k p k, is obtined by pproximting the Hessin by the rst prt of the sum in Eqn. 4, J(x k ) J(x k )p k = J(x k ) h(x k ): () We observe from the previous section tht this is the sme s solving the liner LS problem min p kj(x k )p + h(x k )k : () Prcticl experience hs shown tht if J(x) hs full rnk nd h(x) relly gets smll close to the solution, J(x k ) J(x k ) becomes good pproximtion to the exct Hessin nd Guss-Newton converges bout s fst s full Newton. This is illustrted for Guss-Newton on the Rosenbrock Bnn Function on Fig. 7, copied from Mtlb Optimiztion Toolbox documenttion. The guide sys tht only 48 function evlutions were needed, even when the grdient ws computed numericlly. If, on the contrry, J(x) is rnk de cient, or the residuls (the components of h) re not smll, the performnce my be very poor. This is discussed in N&W, p
12 Figure 7: This gure, copied from the documenttion of Mtlb Optimiztion Toolbox, shows tht Guss-Newton is very e cient when the lst prt of the Hessin mtrix my be ignored. Wheres Guss-Newton is line serch method, it is lso resonble use the sme pproximte Hessin in trust region setting. The qudrtic pproximtion to the objective is then m (p) = f k + J k h k p + p J k J k p; () nd with sphericl domin D with rdius, we hve the sme type of model problem s before, p = rgmin sd m (p) : () If J(x k ) hs full rnk, the solution of will be identicl to the Guss-Newton solution, p GN = Jk J k J k h k ; (4) s long s kp GN k. Otherwise, the solution hs to be on the boundry of D nd obtined from J k J k + k I p = J k h k : (5) with vlue on k tht ensures kpk =. Note tht in this cse, the solution of the constrined problem is strightforwrd (in principle): J k J k, nd J k J k + I > for ll >. Moreover, kp k!!. This trust region lgorithm is clled Levenberg-Mrqurdt s method, nd is the stndrd method for non-liner LS problems nd one of the lgorithm implemented in Optimiztion Toolbox in MATLAB. If it is possible to compute the mtrix B in resonble wy, it will in generl be preferble to use it, s discussed in []. Such methods re implemented in the NAG librry of numericl lgorithms.
13 References [] Gill, P.E. nd W. Murry: Algorithms for the solution of the nonliner lest-squre problems, SIAM J. Num. Anl. Vol. 5(978) pp [] Golub, G.H. nd C.F. VnLon: Mtrix Computtions, Johns Hopkins University Press, 989. [] Luenberger, D.G.: Liner nd Nonliner Progrmming, nd ed., Addison Westley, 984.
1 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 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 informationBest Approximation. Chapter The General Case
Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given
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 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 informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
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 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 informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationQuadratic Forms. Quadratic Forms
Qudrtic Forms Recll the Simon & Blume excerpt from n erlier lecture which sid tht the min tsk of clculus is to pproximte nonliner functions with liner functions. It s ctully more ccurte to sy tht we pproximte
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 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 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 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 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 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 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 informationCMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature
CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy
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 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 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 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 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 informationRecitation 3: Applications of the Derivative. 1 Higher-Order Derivatives and their Applications
Mth 1c TA: Pdric Brtlett Recittion 3: Applictions of the Derivtive Week 3 Cltech 013 1 Higher-Order Derivtives nd their Applictions Another thing we could wnt to do with the derivtive, motivted by wht
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 informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
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 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 informationChapter 3 Solving Nonlinear Equations
Chpter 3 Solving Nonliner Equtions 3.1 Introduction The nonliner function of unknown vrible x is in the form of where n could be non-integer. Root is the numericl vlue of x tht stisfies f ( x) 0. Grphiclly,
More information1B40 Practical Skills
B40 Prcticl Skills Comining uncertinties from severl quntities error propgtion We usully encounter situtions where the result of n experiment is given in terms of two (or more) quntities. We then need
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 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 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 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 informationMath Lecture 23
Mth 8 - Lecture 3 Dyln Zwick Fll 3 In our lst lecture we delt with solutions to the system: x = Ax where A is n n n mtrix with n distinct eigenvlues. As promised, tody we will del with the question of
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 informationChapter 5. Numerical Integration
Chpter 5. Numericl Integrtion These re just summries of the lecture notes, nd few detils re included. Most of wht we include here is to be found in more detil in Anton. 5. Remrk. There re two topics with
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 informationW. We shall do so one by one, starting with I 1, and we shall do it greedily, trying
Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)
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 informationLinearity, linear operators, and self adjoint eigenvalue problems
Linerity, liner opertors, nd self djoint eigenvlue problems 1 Elements of liner lgebr The study of liner prtil differentil equtions utilizes, unsurprisingly, mny concepts from liner lgebr nd liner ordinry
More informationNumerical 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 informationNumerical Integration
Numericl Integrtion Wouter J. Den Hn London School of Economics c 2011 by Wouter J. Den Hn June 3, 2011 Qudrture techniques I = f (x)dx n n w i f (x i ) = w i f i i=1 i=1 Nodes: x i Weights: w i Qudrture
More informationElements of Matrix Algebra
Elements of Mtrix Algebr Klus Neusser Kurt Schmidheiny September 30, 2015 Contents 1 Definitions 2 2 Mtrix opertions 3 3 Rnk of Mtrix 5 4 Specil Functions of Qudrtic Mtrices 6 4.1 Trce of Mtrix.........................
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 informationOrthogonal Polynomials
Mth 4401 Gussin Qudrture Pge 1 Orthogonl Polynomils Orthogonl polynomils rise from series solutions to differentil equtions, lthough they cn be rrived t in vriety of different mnners. Orthogonl polynomils
More informationMultivariate problems and matrix algebra
University of Ferrr Stefno Bonnini Multivrite problems nd mtrix lgebr Multivrite problems Multivrite sttisticl nlysis dels with dt contining observtions on two or more chrcteristics (vribles) ech mesured
More information4.4 Areas, Integrals and Antiderivatives
. res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order
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 informationLecture 3. Limits of Functions and Continuity
Lecture 3 Limits of Functions nd Continuity Audrey Terrs April 26, 21 1 Limits of Functions Notes I m skipping the lst section of Chpter 6 of Lng; the section bout open nd closed sets We cn probbly live
More informationThe Wave Equation I. MA 436 Kurt Bryan
1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string
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 informationARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac
REVIEW OF ALGEBRA Here we review the bsic rules nd procedures of lgebr tht you need to know in order to be successful in clculus. ARITHMETIC OPERATIONS The rel numbers hve the following properties: b b
More informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationIntermediate Math Circles Wednesday, November 14, 2018 Finite Automata II. Nickolas Rollick a b b. a b 4
Intermedite Mth Circles Wednesdy, Novemer 14, 2018 Finite Automt II Nickols Rollick nrollick@uwterloo.c Regulr Lnguges Lst time, we were introduced to the ide of DFA (deterministic finite utomton), one
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 informationLecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.
Lecture 3 3 Solving liner equtions In this lecture we will discuss lgorithms for solving systems of liner equtions Multiplictive identity Let us restrict ourselves to considering squre mtrices since one
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 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 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 informationLecture 13 - Linking E, ϕ, and ρ
Lecture 13 - Linking E, ϕ, nd ρ A Puzzle... Inner-Surfce Chrge Density A positive point chrge q is locted off-center inside neutrl conducting sphericl shell. We know from Guss s lw tht the totl chrge on
More informationLECTURE 19. Numerical Integration. Z b. is generally thought of as representing the area under the graph of fèxè between the points x = a and
LECTURE 9 Numericl Integrtion Recll from Clculus I tht denite integrl is generlly thought of s representing the re under the grph of fèxè between the points x = nd x = b, even though this is ctully only
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 informationJim Lambers MAT 169 Fall Semester Lecture 4 Notes
Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of 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 informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
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 informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationNUMERICAL INTEGRATION
NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationContinuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht
More information1 Probability Density Functions
Lis Yn CS 9 Continuous Distributions Lecture Notes #9 July 6, 28 Bsed on chpter by Chris Piech So fr, ll rndom vribles we hve seen hve been discrete. In ll the cses we hve seen in CS 9, this ment tht our
More informationNOTES ON HILBERT SPACE
NOTES ON HILBERT SPACE 1 DEFINITION: by Prof C-I Tn Deprtment of Physics Brown University A Hilbert spce is n inner product spce which, s metric spce, is complete We will not present n exhustive mthemticl
More informationMain topics for the First Midterm
Min topics for the First Midterm The Midterm will cover Section 1.8, Chpters 2-3, Sections 4.1-4.8, nd Sections 5.1-5.3 (essentilly ll of the mteril covered in clss). Be sure to know the results of the
More informationUNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3
UNIFORM CONVERGENCE Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3 Suppose f n : Ω R or f n : Ω C is sequence of rel or complex functions, nd f n f s n in some sense. Furthermore,
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties
More informationNumerical Integration
Chpter 1 Numericl Integrtion Numericl differentition methods compute pproximtions to the derivtive of function from known vlues of the function. Numericl integrtion uses the sme informtion to compute numericl
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 informationAQA Further Pure 1. Complex Numbers. Section 1: Introduction to Complex Numbers. The number system
Complex Numbers Section 1: Introduction to Complex Numbers Notes nd Exmples These notes contin subsections on The number system Adding nd subtrcting complex numbers Multiplying complex numbers Complex
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 informationand that at t = 0 the object is at position 5. Find the position of the object at t = 2.
7.2 The Fundmentl Theorem of Clculus 49 re mny, mny problems tht pper much different on the surfce but tht turn out to be the sme s these problems, in the sense tht when we try to pproimte solutions we
More informationChapter 5 : Continuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 216 Néhémy Lim Chpter 5 : Continuous Rndom Vribles Nottions. N {, 1, 2,...}, set of nturl numbers (i.e. ll nonnegtive integers); N {1, 2,...}, set of ll
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 informationZ b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...
Chpter 7 Numericl Methods 7. Introduction In mny cses the integrl f(x)dx cn be found by finding function F (x) such tht F 0 (x) =f(x), nd using f(x)dx = F (b) F () which is known s the nlyticl (exct) solution.
More informationOverview of Calculus I
Overview of Clculus I Prof. Jim Swift Northern Arizon University There re three key concepts in clculus: The limit, the derivtive, nd the integrl. You need to understnd the definitions of these three things,
More informationMORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)
MORE FUNCTION GRAPHING; OPTIMIZATION FRI, OCT 25, 203 (Lst edited October 28, 203 t :09pm.) Exercise. Let n be n rbitrry positive integer. Give n exmple of function with exctly n verticl symptotes. Give
More information12 TRANSFORMING BIVARIATE DENSITY FUNCTIONS
1 TRANSFORMING BIVARIATE DENSITY FUNCTIONS Hving seen how to trnsform the probbility density functions ssocited with single rndom vrible, the next logicl step is to see how to trnsform bivrite probbility
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct
More informationMath 61CM - Solutions to homework 9
Mth 61CM - Solutions to homework 9 Cédric De Groote November 30 th, 2018 Problem 1: Recll tht the left limit of function f t point c is defined s follows: lim f(x) = l x c if for ny > 0 there exists δ
More informationg i fφdx dx = x i i=1 is a Hilbert space. We shall, henceforth, abuse notation and write g i f(x) = f
1. Appliction of functionl nlysis to PEs 1.1. Introduction. In this section we give little introduction to prtil differentil equtions. In prticulr we consider the problem u(x) = f(x) x, u(x) = x (1) where
More informationEuler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )
Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s
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 dierentil eqution (ODE) du f(t) dt with initil condition u() : Just
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 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 informationECON 331 Lecture Notes: Ch 4 and Ch 5
Mtrix Algebr ECON 33 Lecture Notes: Ch 4 nd Ch 5. Gives us shorthnd wy of writing lrge system of equtions.. Allows us to test for the existnce of solutions to simultneous systems. 3. Allows us to solve
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 informationThe First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).
The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples
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 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 information