A rapidly convergent descent method for minimization
|
|
- Luke Cook
- 6 years ago
- Views:
Transcription
1 A rapdly convergent descent method for mnmzaton By R. Fletcher and M. J. D. Powell A powerful teratve descent method for fndng a local mnmum of a functon of several varables s descrbed. A number of theorems are proved to show that t always converges and that t converges rapdly. Numercal tests on a varety of functons confrm these theorems. The method has been used to solve a system of one hundred non-lnear smultaneous equatons. 1. Introducton We are concerned n ths paper wth the general problem of fndng an unrestrcted local mnmum of a functon J[x u x 2, x n ) of several varables x x, x 2,., x n. We suppose that the functon of nterest can be calculated at all ponts. It s convenent to group functons nto two man classes accordng to whether the gradent vector g, = Wbx, s defned analytcally at each pont or must be estmated from the dfferences of values of/. The method descrbed n ths paper s applcable to the case of a defned gradent. For the other case a useful method and general dscusson are gven by Rosenbrock (196). Methods usng the gradent nclude the classcal method of steepest descents (Courant, 1943; Curry, 1944; and Householder, 1953), Levenberg's modfcaton of damped steepest descents (1944), a somewhat smlar varaton due to Booth (1957), the conjugate gradent method of Hestenes and Stefel (1952), smlar methods of Martn and Tee (1961), the "Partan" method of Shah, Buehler and Kempthorne (1961), and a method due to Powell (1962). In ths paper we descrbe a powerful method wth rapd convergence whch s based upon a procedure descrbed by Davdon (1959). Davdon's work has been lttle publczed, but n our opnon consttutes a consderable advance over current alternatves. We have made both a smplfcaton by whch certan orthogonalty condtons whch are mportant to the rate of attanng the soluton are preserved, and also an mprovement n the crteron of convergence. Because, near the mnmum, the second-order terms n the Taylor seres expanson domnate, the only methods whch wll converge quckly for a general functon are those whch wll guarantee tofnd the mnmum of a general quadratc speedly. Only the latter four methods of the last paragraph do ths, and the procedures of Hestenes and Stefel and of Martn and Tee are not applcable to a general functon. Of course the generalzed Newton-Raphson method (Householder, 1953) has fast convergence eventually, but t requres second dervatves of the functon to be evaluated, and frequently fals to converge from a poor approxmaton to the mnmum. The method descrbed has quadratc convergence and s superor to "Partan" and to Powell's method, both n that t makes use of nformaton determned by prevous teratons and also n that each teraton s quck and smple to carry out. 163 Furthermore, t yelds the curvature of the functon at the mnmum, so excellent tests for convergence and estmates of varance can be made. The method s gven an elegant theoretcal bass, and proofs of stablty and of the rate of convergence are ncluded. The results of numercal tests wth a varety of functons are also gven. These confrm that the method s probably the most powerful general procedure for fndng a local mnmum whch s known at the present tme. 2. taton It s convenent to descrbe the method n terms of the Drac bra-ket notaton (Drac, 1958) appled to real vectors. In ths notaton the column vector (x u x 2., x n ) s wrtten as x>. The row vector wth these same elements s denoted by <x. The scalar product of <x and j> s wrtten (x\y} and we may note that <x\y~) s= 2x,-y, = j,x; s= <j/ x>. / The constructon.x><j>, however, denotes a lnear operator wth matrx elements D } = x,y>j so that.x><j> #.F><*. A general lnear operator or matrx wll be denoted by a captal letter n bold type. It then follows that say H\xy s a column vector, (x\h s a row vector and (x\h\yy s a scalar. We reserve / to denote the functon of nterest, x> to denote ts arguments and g> to denote ts gradent. Hence the standard quadratc form n n dmensons / = /o + a,x, = 1 becomes n ths notaton and also \g> = = I j= 1 G\x>. 3. The method If we consder the quadratc form (1) then, gven the matrx G u = 'b 2 fl'bx ( 'bxj, we can calculate the dsplacement between the pont \x} and the mnmum \x o y as x o >- x>=-g-' g>. (3) In ths method the matrx G~' s not evaluated drectly; (1) (2)
2 Descent method for nstead a matrx H s used whch may ntally be chosen to be any postve defnte symmetrc matrx. Ths matrx s modfed after the th teraton usng the nformaton ganed by movng down the drecton \sfy = - J5TV> (4) n accordance wth (3). The modfcaton s such that y, the step to the mnmum down the lne s effectvely an egenvector of the matrx H +l G. Ths ensures that as the procedure converges H tends to G~ l evaluated at the mnmum. It s convenent to take the unt matrx ntally for H so that the frst drecton s down the lne of steepest descent. Let the current pont be x'> wth gradent g'> and matrx H'. The teraton can then be stated as follows. Set s<> = - H^g'y. Obtan a.' such that f(\x'} + a' s'>) s a mnmum wth respect to A along \x'} + A s'> and a' >. We wll prove that a' can always be chosen to be postve. Set a ( > = a-y>. (5) Set x' + '> = x' Evaluate /( x /+I» and orthogonal to a'y, that s notng that > s Set Set where and y\h'\yy Set / = + 1 and repeat. There are two obvous and very useful ways of termnatng the procedure, and they arse because ^'> tends to the correcton to \x*y. One s to stop when the predcted absolute dstance from the mnmum (s^1)* s less than a prescrbed amount, and the other s to fnsh when every component of s'> s less than a prescrbed accuracy. Two addtonal safeguards have been found necessary n automatc computer programs. The frst s to work through at least n (the number of varables) teratons, and the second s to apply the tests to CT'> as well as to \s'y. The method of obtanng the mnmum along a lne s not central to the theory. The suggested procedure gven n the Appendx, whch uses cubc nterpolaton, s based on that gven by Davdon, and has been found satsfactory. We shall now show that the process s stable, and demonstrate that f/( x» s the quadratc form (1) then (6) (7) 164 mnmzaton the procedure termnates n n teratons. We shall also explan the theoretcal justfcaton for the manner n whch the matrx H s modfed. 4. Stablty It s usual for descent methods to be stable because one ensures that the functon to be mnmzed s decreased by each step. It wll be shown n ths Secton that the drecton of search \s'}, defned by equaton (4), s downhll, so «' can always be chosen to be postve. Because g'> s the drecton of steepest ascent, the drecton \s'} wll be downhll f and only f s postve. We wsh the drecton of search to be downhll for all possble \g'y so we must prove that H' s postve defnte. Because H has been chosen to be postve defnte an nductve argument wll be used. In the proof t s assumed that H l s postve defnte and consequently that a' s postve. It s proved that, for any x>, (x\h +x \xy >. We may defne \py = (H''Y\xy and?>= (#')* /> as the square root of a postve defnte matrx exsts. From (7) <y\h'\yy on account of Schwartz's nequalty. But <a' /> = <a' g'+ 1 > from (6) from (4) and (5) Hence <x ^l+l x> > for all non-trval x>. Therefore H + s postve defnte and the procedure s stable. 5. Quadratc convergence In ths Secton t s assumed that / s the quadratc form (1) and that /has a well defned mnmum. It s proved that n ths case the method fnds the mnmum n n teratons. The method of proof s to show that o- >, CT'>,..., o*> are lnearly ndependent egenvectors of H k + > G wth egenvalue unty. Therefore t wll follow that H"G s the unt matrx. By defnton from (2) = G\a>y. (8)
3 Also from (8) The equatons Descent method for mnmzaton Ths result can be proved from the orthogonalty condtons (1) because these mply that S'GS = A, where S s the matrx of vectors <x'>, and A s a dagonal matrx - H'\y'> by usng (7) wth elements <a' (? <7'>. = //' y> + = k>. (9) <o \G\o>> = <<j<k (1) and H k G\a>> = a'> < / < k (11) wll now be consdered. It s clear from (9) that they are true f k 1. It wll be proved that f they are true for k they are true for k + 1. From (2) = \a> Hence by defnton G = Therefore G- =SA 1 S / and as A s a dagonal matrx ths reduces to Therefore from the defnton of A' and equaton (8), equaton (15) s proved. The form of the term B' can be deduced because equaton (9) must be vald. For a quadratc we must have Therefore from (1) and (6) Hence from (11) = < < k. > = so from (4) and (5) -a'<a' G a*> =. Therefore Also from (8), (11) and (13) = < < k. (12) (13) Therefore usng the above result and equatons (7), (11) and (13) y = H k G\a l y <<k. (14) Equatons (9), (13) and (14) prove the nducton. Equaton (1) proves that the vectors a >, ICT 1 ),..., <r"~'> are lnearly ndependent and therefore H" = G~ l. That the mnmum s found by n teratons s proved by equaton (12). g"> must be orthogonal to CT >, CT'>,..., O"-'> whch s only possble f g"> s dentcally zero. 6. Improvng the matrx H The matrx H' s modfed by addng to t two terms A' and B'. A' s the factor whch makes H tend to G~ l n the sense that for a quadratc G-< = n-i A 1. (15) /=o 165 Therefore as A'G\&y = a'> the equaton B'G\a l y = - H'G a'> = - //' /> must be satsfed. Ths mples that the smplest form for B> s. = _ H'\yy <z\ and as B' s to be symmetrc ths gves. = _ <y\h'\y'y' Although Davdon's method nvolves these relatons, some of the other deas used by hm can cause H not to tend to G~' even n the quadratc case. The effect n the non-quadratc case would depend upon the functon n queston but mght well lead to slower convergence. 7. Numercal results comparson wth other procedures As a comparson wth other methods we use the functon gven by Rosenbrock ftxy, x 2 ) = 1(x 2 - x]) 2 + (1 - x,) 2 startng at ( 1-2, 1-). Ths functon s dffcult to mnmze on account of ts havng a steep sded valley followng the curve x\ = x 2. Eghteen teratons were requred to reach the mnmum, each one requrng the mnmum to be calculated n only one drecton. Table 1 shows how ths procedure compares wth the classcal steepest descent method and Powell's method, one of the procedures wth quadratc convergence. The table takes nto account that the latter method requres mnma to be found n three drectons for each teraton. It wll be seen that ths method s consderably more effcent than that of Powell, both of these beng far more effcent than steepest descents.
4 Table 1 A comparson n two dmensons Descent method for mnmzaton Table 2 A quadratc functon EQUIVALENT n STEEPEST DESCENTS /(*!. X2) ' POWELL'S METHOD fx u x 2 ) x x lo" 9 OUR METHOD f{x\,x 2 ) x 1-8 ITERATION X f H \ > A l : I o-' 5 1 A smlar comparson was made wth the functon gven by Powell:,, x 2, x 3, x 4 ) = (x, + 1x 2 ) 2 + 5(x 3 x 4 ) 2 (x 2-1(x, - x 4 )< startng at (3, 1,,. 1). In sx teratons the method reduced /from 215 to 2-5 x 1~ 8. Powell's method took the equvalent of seven teratons to reach 9 x 1~ 3, whereas steepest descents only reached 6-36 n seven teratons. The method also brought out the sngularty of G at the mnmum of/, the elements of H becomng ncreasngly large. To compare ths varaton of Davdon's method wth hs orgnal method the smple quadratc /(*1, *2> = A~ 2 *1*2 + 2x1 was used. The complete progress of the method descrbed s gven n Table 2, showng that t does termnate n two teratons and that H does converge to G~ { whch for ths functon s It wll be notced also, as proved, that C"" 1 = "LA!. In Davdon's method, although a value of/of smlar order of magntude had been reached n two teratons, Hhad only reached /-95-47\ \-47-48A Ths was due to one of the alternatves allowed by Davdon. Also hs procedure for termnatng the process was unsatsfactory, and the computaton had to be stopped manually. A non-quadratc test n three dmensons was also 166 made, by usng a functon wth a steep sded helcal valley. Ths functon /(x,, x 2, x 3 ) = 1{[x 3 - lo(x,, x 2 )] 2 + [r(x {, x 2 ) - I] 2 } + x\ where 2TT6(X U X 2 ) = arctan (x 2 /x,), x t > = 77 + arctan (x 2 /x), x t < and / (*!, x 2 ) = (x 2 + x 2 )* has a helcal valley n the x 3 drecton wth ptch 1 and radus 1. It s only consdered for TT/2 <2TT<3 77/2 that s 2-5<x 3 <7-5. It has a mnmum at the pont (1,, ). Both methods were started from ( 1,, ) and H set to the unt matrx. The method gven n ths paper converged n eghteen teratons, whereas Davdon's method requred only ten. However, on account of the more complcated nature of Davdon's teratons, the mnmum often beng sought along more than one drecton n a sngle teraton, the tme taken by the two procedures was almost dentcal. The progress of ths method on the functon s gven n Table Numercal results functons of a large number of varables Tests were also made to fnd out whether the method s sutable for fndng the mnmum of a functon of a large number of varables. In these tests the Stretch computer was used to solve non-lnear smultaneous equatons n up to a hundred varables. The equatons were n S Ajj sn ocj + Bjj cos a,- = E, = 1, 2,..., n
5 n Table 3 A functon wth a steep-sded helcal valley * o- 5 X o- 5 so that the functon to be mnmzed was Descent method for mnmzaton / 2-5 x IO x X x IO- 4 3 x IO- 6 7 x IO- 8 /= {E, - {A lt sn ccj + B,j cos a,)} 2. = j The matrx elements of A and B were generated as random ntegers between 1 and +1, and the values of the varables <x h = 1, 2,..., n were generated randomly between n and -n. For these values the rght-hand sdes of the equatons, E h were worked out. The method of ths paper was appled to fnd optmum values of «,- startng from (a, -fo-ls,) where the 8,'s were also generated as random numbers between IT and 7T. In each run the crteron for convergence was that every a should be found to accuracy 1. The method was entrely successful. Table 4 shows that the number of tmes/and ts dervatves had to be calculated was approxmately lnear n the number of varables. The total tme taken for all the runs was ffteen mnutes, ten mnutes of whch was spent on the fnal case. That a dfferent mnmum was found on fve occasons was not surprsng because f A = B t may be shown that there are up to 2" real solutons to the equatons such that a, < n. Ths abundance of mnma emphaszes the power of the method because n every case t converged to a reasonable soluton. The progress of these tests s nterestng. For the frst n teratons the changes n the functon were smlar to those experenced wth the method of steepest descents, that s a substantal change occurred ntally due to descendng nto a nearby valley, after whch convergence was slow. However, after n teratons had been completed 167 Table 4 Applcaton to a functon of many varables n NO. OF TIMES /EVALUATED WHETHER EXPECTED MINIMUM FOUND a good approxmaton to the fnal matrx H had been accumulated, after whch the functon was decreased substantally at each teraton. For example n the hundred-varable tral the functon to be mnmzed was decreased from to n the frst ten teratons, and to after one hundred teratons. After 12 teratons t was down to 1342, and after 14 to 147. The functon was reduced to -44 by 16 teratons, and the mnmum was found on the 162 nd. The second ffty-varable tral was even more strkng. Ten teratons reduced the functon from 2538 to 4264, ffty teratons reduced t to 3526, and a further ten teratons reduced t to 27. The concluson to be drawn from ths behavour s that for many applcatons of the method a substantal number of the teratons requred wll be spent on settng up the nverse of the matrx of second dervatves. Therefore, f a good postve defnte approxmaton to H can be calculated ntally, as s the case when the method s beng appled to solvng smultaneous equatons, then ths approxmaton should be chosen for H. 9. Concluson The numercal examples show clearly that the type of method gven by Davdon s consderably superor to other methods prevously avalable. The smplfcatons we have made enable programs to be wrtten more easly, and they do not seem to mpar the speed of convergence. It s obvously practcable to apply ths method to fnd a local mnmum of a general functon of a large number of varables whose frst dervatves can be evaluated quckly, even f only poor ntal approxmatons to a soluton are known.
6 Descent method for mnmzaton 1. Acknowledgement where w = (z 2 <,g x \s l y (g y \s'y)* One of us, R.F., s ndebted to Dr. C. M. Reeves for hs constant help and encouragement, and also to the D.S.I.R. for the provson of a research studentshp. Appendx The mnmum on a lne A smple algorthm s gven for estmatng the parameter a.'. A pont y> s chosen on \x 1 } + A j'> wth A >. Let f x, \g x y,f y and \g y y denote the values of the functon and gradent at the ponts \x'} and />. Then an estmate of a' can be formed by nterpolatng cubcally, usng the functon values f x and f y and the components of the gradents along 5'>. Ths s gven by a' "A References <gyw> + W Z = 1 - <gy\s'> - 2W and z = ^ (f x - f y ) + <g x \s>> + <g, * >. A sutable choce of the pont \y'y s gven by h / ' \' = MINIMUM OF (l, /o s the predcted lower bound of/( x», for example zero n least-squares calculatons. Ths value of -q ensures that the choce of /> s reasonable. It s necessary to check that /( *'> + a' s'» s less than both f x and f y. If t s not, the nterpolaton must be repeated over a smaller range. Davdon suggests one should ensure that the mnmum s located between \x! ) and y> by testng the sgn of <^ ^'> and comparng f x and f y before nterpolatng. The reader s referred to Davdon's report for more extensve detals of ths stage. BOOTH, A. D. (1957). Numercal Methods, London: Butterworths. COURANT, R. (1943). "Varatonal methods for the soluton of problems of equlbrum and vbratons," Bull. Amer. Math. Soc, Vol. 49, p. 1. CURRY, H. D. (1944). "The method of steepest descent for non-lnear mnmzaton problems;" Qu. App. Maths., Vol. 2, p DAVIDON, W. C. (1959). "Varable metrc method for mnmzaton," A.E.C. Research and Development Report, ANL-599 (Rev.). DIRAC, P. A. M. (1958). The Prncples of Quantum Mechancs, Oxford: O.U.P. HESTENES, M. R., and STIEFEL, E. (1952). "Methods of conjugate gradents for solvng lnear systems," /. Res. N.B.S., Vol. 49, p. 49. HOUSEHOLDER, A. S. (1953). Prncples of Numercal Analyss, New York: McGraw-Hll. LEVENBERG, K. (1944).."A method for the soluton of certan non-lnear problems n least squares," Qu. App. Maths., Vol. 2 p MARTIN, D. W., and TEE, G. J. (1961). "Iteratve methods for lnear equatons wth symmetrc postve defnte matrx," The Computer Journal, Vol. 4, p POWELL, M. J. D. (1962). "An teratve method for fndng statonary values of a functon of several varables," The Computer Journal, Vol. 5, p ROSENBROCK, H. H. (I96). "An automatc method for fndng the greatest or least value of a functon," The Computer Journal, Vol. 3, p SHAH, B. V., BUEHLER, R. J., and KEMPTHORNE, O. (1961). "The method of parallel tangents (Partan) for fndng an optmum," Offce of Naval Research Report, NR (. 2). Book Revew (.contnued from p. 143) here. It seems entrely wrong that modern source languages should be domnated by the sequental regmes of early machne codes, and any move towards conventonal mathematcal forms s to be welcomed. Agan, the approach here seems rather tentatve and some major benefts are lost. I thnk t s preferable to make sequental codng subordnate to defntons rather than the other way round: here les the key to the very mportant problem of ntegratng the translator wth a realstc operatng system. Amongst the other problems tackled are the handlng of complex varables,- recurrence relatons, and drect transfer of control to parts of'the program not smlarly accessble n ALGOL. It may be regarded as a trbute to ALGOL that an attempt has been made to graft such a system onto the same tree. At the same tme t- must accept a measure of responsblty for the fact that the above deas were not more fully developed and n use three years ago. J. K. ILIFFE. 168
1 GSW Iterative Techniques for y = Ax
1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus
More informationSingular Value Decomposition: Theory and Applications
Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationModule 9. Lecture 6. Duality in Assignment Problems
Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept
More informationFoundations of Arithmetic
Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an
More informationStanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011
Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationIV. Performance Optimization
IV. Performance Optmzaton A. Steepest descent algorthm defnton how to set up bounds on learnng rate mnmzaton n a lne (varyng learnng rate) momentum learnng examples B. Newton s method defnton Gauss-Newton
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationVector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence.
Vector Norms Chapter 7 Iteratve Technques n Matrx Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematcs Unversty of Calforna, Berkeley Math 128B Numercal Analyss Defnton A vector norm
More informationTopic 5: Non-Linear Regression
Topc 5: Non-Lnear Regresson The models we ve worked wth so far have been lnear n the parameters. They ve been of the form: y = Xβ + ε Many models based on economc theory are actually non-lnear n the parameters.
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationYong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )
Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often
More informationA New Refinement of Jacobi Method for Solution of Linear System Equations AX=b
Int J Contemp Math Scences, Vol 3, 28, no 17, 819-827 A New Refnement of Jacob Method for Soluton of Lnear System Equatons AX=b F Naem Dafchah Department of Mathematcs, Faculty of Scences Unversty of Gulan,
More informationFeature Selection: Part 1
CSE 546: Machne Learnng Lecture 5 Feature Selecton: Part 1 Instructor: Sham Kakade 1 Regresson n the hgh dmensonal settng How do we learn when the number of features d s greater than the sample sze n?
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationInexact Newton Methods for Inverse Eigenvalue Problems
Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationρ some λ THE INVERSE POWER METHOD (or INVERSE ITERATION) , for , or (more usually) to
THE INVERSE POWER METHOD (or INVERSE ITERATION) -- applcaton of the Power method to A some fxed constant ρ (whch s called a shft), x λ ρ If the egenpars of A are { ( λ, x ) } ( ), or (more usually) to,
More informationFor now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.
Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationChapter - 2. Distribution System Power Flow Analysis
Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationLecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem.
prnceton u. sp 02 cos 598B: algorthms and complexty Lecture 20: Lft and Project, SDP Dualty Lecturer: Sanjeev Arora Scrbe:Yury Makarychev Today we wll study the Lft and Project method. Then we wll prove
More informationOn a direct solver for linear least squares problems
ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationReview of Taylor Series. Read Section 1.2
Revew of Taylor Seres Read Secton 1.2 1 Power Seres A power seres about c s an nfnte seres of the form k = 0 k a ( x c) = a + a ( x c) + a ( x c) + a ( x c) k 2 3 0 1 2 3 + In many cases, c = 0, and the
More information8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS
SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationSupplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso
Supplement: Proofs and Techncal Detals for The Soluton Path of the Generalzed Lasso Ryan J. Tbshran Jonathan Taylor In ths document we gve supplementary detals to the paper The Soluton Path of the Generalzed
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationSIO 224. m(r) =(ρ(r),k s (r),µ(r))
SIO 224 1. A bref look at resoluton analyss Here s some background for the Masters and Gubbns resoluton paper. Global Earth models are usually found teratvely by assumng a startng model and fndng small
More informationU.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017
U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that
More informationDepartment of Chemical and Biological Engineering LECTURE NOTE II. Chapter 3. Function of Several Variables
LECURE NOE II Chapter 3 Functon of Several Varables Unconstraned multvarable mnmzaton problem: mn f ( x), x R x N where x s a vector of desgn varables of dmenson N, and f s a scalar obectve functon - Gradent
More information2.3 Nilpotent endomorphisms
s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationLeast Squares Fitting of Data
Least Squares Fttng of Data Davd Eberly Geoetrc Tools, LLC http://www.geoetrctools.co/ Copyrght c 1998-2015. All Rghts Reserved. Created: July 15, 1999 Last Modfed: January 5, 2015 Contents 1 Lnear Fttng
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationC/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1
C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned
More informationOPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming
OPTIMIATION Introducton ngle Varable Unconstraned Optmsaton Multvarable Unconstraned Optmsaton Lnear Programmng Chapter Optmsaton /. Introducton In an engneerng analss, sometmes etremtes, ether mnmum or
More informationPrimer on High-Order Moment Estimators
Prmer on Hgh-Order Moment Estmators Ton M. Whted July 2007 The Errors-n-Varables Model We wll start wth the classcal EIV for one msmeasured regressor. The general case s n Erckson and Whted Econometrc
More informationFeb 14: Spatial analysis of data fields
Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationOn the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros
Appled Mathematcal Scences, Vol. 5, 2011, no. 75, 3693-3706 On the Interval Zoro Symmetrc Sngle-step Procedure for Smultaneous Fndng of Polynomal Zeros S. F. M. Rusl, M. Mons, M. A. Hassan and W. J. Leong
More informationNorms, Condition Numbers, Eigenvalues and Eigenvectors
Norms, Condton Numbers, Egenvalues and Egenvectors 1 Norms A norm s a measure of the sze of a matrx or a vector For vectors the common norms are: N a 2 = ( x 2 1/2 the Eucldean Norm (1a b 1 = =1 N x (1b
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationCSCE 790S Background Results
CSCE 790S Background Results Stephen A. Fenner September 8, 011 Abstract These results are background to the course CSCE 790S/CSCE 790B, Quantum Computaton and Informaton (Sprng 007 and Fall 011). Each
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationBeyond Zudilin s Conjectured q-analog of Schmidt s problem
Beyond Zudln s Conectured q-analog of Schmdt s problem Thotsaporn Ae Thanatpanonda thotsaporn@gmalcom Mathematcs Subect Classfcaton: 11B65 33B99 Abstract Usng the methodology of (rgorous expermental mathematcs
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationj) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1
Random varables Measure of central tendences and varablty (means and varances) Jont densty functons and ndependence Measures of assocaton (covarance and correlaton) Interestng result Condtonal dstrbutons
More informationFormulas for the Determinant
page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use
More informationSalmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2
Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to
More informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
More informationTime-Varying Systems and Computations Lecture 6
Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy
More informationQuantum Mechanics I - Session 4
Quantum Mechancs I - Sesson 4 Aprl 3, 05 Contents Operators Change of Bass 4 3 Egenvectors and Egenvalues 5 3. Denton....................................... 5 3. Rotaton n D....................................
More information= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.
Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationQuantum Mechanics for Scientists and Engineers. David Miller
Quantum Mechancs for Scentsts and Engneers Davd Mller Types of lnear operators Types of lnear operators Blnear expanson of operators Blnear expanson of lnear operators We know that we can expand functons
More informationOn the correction of the h-index for career length
1 On the correcton of the h-ndex for career length by L. Egghe Unverstet Hasselt (UHasselt), Campus Depenbeek, Agoralaan, B-3590 Depenbeek, Belgum 1 and Unverstet Antwerpen (UA), IBW, Stadscampus, Venusstraat
More informationMATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS
MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS These are nformal notes whch cover some of the materal whch s not n the course book. The man purpose s to gve a number of nontrval examples
More information1 Convex Optimization
Convex Optmzaton We wll consder convex optmzaton problems. Namely, mnmzaton problems where the objectve s convex (we assume no constrants for now). Such problems often arse n machne learnng. For example,
More informationwhere the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt
Physcs 543 Quantum Mechancs II Fall 998 Hartree-Fock and the Self-consstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes.
More informationElectronic Quantum Monte Carlo Calculations of Energies and Atomic Forces for Diatomic and Polyatomic Molecules
RESERVE HIS SPACE Electronc Quantum Monte Carlo Calculatons of Energes and Atomc Forces for Datomc and Polyatomc Molecules Myung Won Lee 1, Massmo Mella 2, and Andrew M. Rappe 1,* 1 he Maknen heoretcal
More informationProblem Set 9 Solutions
Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationChapter 3 Differentiation and Integration
MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v10157-01-00-y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION
More informationLimited Dependent Variables
Lmted Dependent Varables. What f the left-hand sde varable s not a contnuous thng spread from mnus nfnty to plus nfnty? That s, gven a model = f (, β, ε, where a. s bounded below at zero, such as wages
More information