Difference Equations
|
|
- Arron Sullivan
- 5 years ago
- Views:
Transcription
1 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 5a =7 (1) true for each nonnegatve nteger. An equaton of ths type s called a dfference equaton, and our man am of ths artcle s to explan how to solve t (.e. fnd a formula for computng the ndvdual terms of the sequence). We wll deal only wth the smplest case of such equatons, namely those lnear n all the a s; furthermore, we allow the a s to have only constant coeffcents (such as our above example - note that the a -terms are usually collected on the left hand sde of the equaton). If the rght hand sde of the equaton s zero, t s called homogeneous, otherwse (as n our example) the equaton s non-homogeneous. We wll be able to solve a non-homogeneous equaton only when the rght hand sde s a polynomal n, further multpled by a constant rased to the power of, e.g. ( 4). Note that (1) s of ths type - the constant s smply equal to 1. Homogeneous case - Introducton Havng an equaton of the a + +5a +1 3a =0 () type, we attempt to solve t by assumng that the soluton has the followng form (hopng that our guess wll prove correct): a tral = λ where λ s a constant whose exact value remans to be establshed. Substtutng ths tral soluton nto () yelds λ + +5λ +1 3λ =0 1
2 or, equvalently (dvdng the prevous equaton by λ ), the followng so called characterstc equaton (ts LHS s the correspondng characterstc polynomal): λ +5λ 3=0 The last equaton has clearly only two roots, namely λ 1 = 1 and λ = 3. Ths mples that each a =( 1 ) and a =( 3) consttutes a dfferent soluton to () - we wll call them (stretchng the usual termnology a bt) basc solutons. One can show that the fully general soluton of () can be bult by takng a lnear combnaton of the two basc solutons, thus: a general = A ( 1 ) + B ( 3) where A and B are two arbtrary numbers. To fx the values of A and B and construct a specfc soluton to (), we need to be gven numercal values of exactly two members of the a -sequence. In the case of an equaton wth two basc solutons, there are two common possbltes:.1 Intal condtons s the name by whch we refer to the stuaton when the frst two members of the sequence (usually a 0 and a 1 ) are gven, e.g. a 0 = 1 and a 1 =. Ths clearly leads to two ordnary lnear equatons for A and B, namely (n the case of the prevous example): A + B = 1 A 3B = whch can be easly solved to yeld: A = 7 and B = 5 7. The specfc soluton whch meets not only () but also the two ntal condtons s then: a specfc = 7 µ ( 3) basedonwhchwecaneaslyfnd the value of any a, e.g. a 10 = ( 3) = 51. Note that n ths case (of gven ntal condtons), we can also fnd a 10 (or any other specfc a ) by the usual recursve procedure,.e. a = 5 a a 0 = 13 a 3 = 5 a + 3 a 1 = 77 4.
3 (no formula necessary). As a second, more nterestng example, let us try to fnd a formula for the determnant of the followng tr-dagonal matrx: α β β α β β α β... 0 a det. 0 0 β α β β α where s the number of rows (and columns). By expandng t along the frst row, one gets: a = α a 1 β a wth a 1 = α and a = α β. Clearly, a + α a +1 + β a =0 s an equvalent way of presentng the same set of equatons. The roots of the correspondng characterstc poynomal are λ 1 = α + p α 4β and λ = α p α 4β The general soluton to (??) sthus a = A λ 1 + B λ The ntal condtons mply A λ 1 + B λ = α A λ 1 + B λ = α β whch yelds A = λ 1 and B = λ λ 1 λ λ 1 λ The fnal, specfc solutonsthus a = λ+1 1 λ +1 λ 1 λ whch can be wrtten n a more explct form of a = 1 [/] X +1 j+1 α j (α 4β ) j (3) j=0 where [/] mples the nteger part of /. 3
4 . Boundary condtons We wll now return to solvng (), but ths tme the frst and the last value of the sequence are gven, e.g. a 0 = 97 3 and a = 3. Smlarly to dealng wth ntal condtons, ths leads to two ordnary lnear equatons for A and B, namely A + B = 97 A 10 + B ( 3)10 = whch yeld A =3and B = 1 3. The correspondng specfc soluton then reads a specfc =3 µ ( 3) = µ 5 1 ( 3) 1 Note that n ths case the recursve procedure would not work, and the formula soluton s thus the only way to solve the problem. 3 Homogeneous case - Complcatons A natural extenson of our prevous examples s to relate more than three consecutve values of the a sequence, e.g. a +3 3a + a +1 +3a =0 (4) It s not to dffcult to fgure out that the correspondng characterstc equaton s λ 3 3λ λ +3=0 ths tme havng three roots, namely λ 1 =3,λ =1and λ 3 = 1. The general soluton to (4) s thus a = A 3 + B + C ( 1) To fnd a specfc soluton, three dstnct values of the a sequence must be gven; f these are a 0,a 1 and a, we can stll call t an ntal-value problem (and have the opton of solvng t recursvely), but there s no clear-cut verson of boundary condtons. In general, t s mportant to realze that there are several equvalent ways of presentng (4), such as a 3a 1 a +3a 3 =0 (the only dfference beng that now one would say: true for 3). 4
5 In ths context, one should watch out for skpped ndces, e.g. a 3a 1 +a 3 =0 (5) whose characterstc polynomal s λ 3 3λ + (no λ term!) wth roots of 1, 1+ 3 and 1 3. The general soluton to (5) s thus a = A + B (1 + 3) + C (1 3) (6) 3.1 Double and multple roots An obvous complcaton arses when two (or more) roots of the characterstc polynomal are equal to each other, e.g. a +3 4a + 3a a =0 (7) whose characterstc polynomal has the followng roots:, 3 and 3. Clearly, ( ) and 3 are stll two basc solutons of ths dfference equaton, but where s the thrd? One can easly verfy that, n a case lke ths, multplyng 3 by creates yet another possble soluton to (7) - check t out! The fully general soluton s then constructed as a lnear combnaton of these: a = A ( ) + B 3 + C 3 (8) Convertng t to a specfc soluton for a gven set of ntal values, e.g. a 0 =, a 1 =0and a = 1 s done n the usual way: A + B = A +3B +3C = 0 4A +9B +18C = 1 The correspondng soluton s: A = ,B= 5 and C = 15, whch leads to a spec. = 17 5 ( ) Smlarly, n a case of a trple root (say, equal to 3), we would construct the correspondng three basc solutons by frst takng ( 3), then multplyng ths by, and fnally multplyng t by. The general pattern should now be obvous. Thus, for example, f the characterstc polynomal has the followng set of roots:, 3, 3, 4, 4 and 4, the general soluton s: a = A + B 3 + C 3 + D ( 4) + E ( 4) + F ( 4) Note that some textbooks, nstead of multplyng a multple root by,, 3, etc. use the followng alternate selecton:, ( 1), ( 1)( ),... The latter approach has some conceptual and perhaps even computatonal advantages, but here we prefer usng the former. The two schemes are fully equvalent. 5
6 3. Complex roots Snce ths secton deals wth complex numbers, we need the symbol to denote the purely magnary unt. We wll thus have to swtch from the old a to usng a n. Later on, we wll return to our regular notaton. As we all know, roots of a polynomal can easly turn out to be complex (comng n complex-conjugate pars). Formally, our prevous soluton remans correct, for example: a n+ 4a n+1 +13a n =0 (9) has a characterstc polynomal wth roots of +3 and +3. We can stll wrte a n = A ( + 3) n + A ( 3) n where A and A are complex conjugates of each other (to yeld a real answer), but t s often more convenent to avod complex numbers and use, nstead of (+3) n and ( 3) n, ther real and purely magnary parts (qute legtmately, as we can always replace any two basc solutons by ther lnear combnatons). Ths s acheved more easly by convertng +3 to ts polar representaton, thus: cos(arctan 3 )+sn(arctan 3 ) The real and purely magnary parts of ( + 3) n are then, respectvely: 13 n/ cos(n arctan 3 ) and 13 n/ sn(n arctan 3 ) The general soluton to (9) can then be wrtten (avodng complex numbers) thus: a n =13 n/ A sn(n arctan 3 )+B cos(n arctan 3 ) Luckly, many practcal stuatons usually manage to steer clear of the complex case. Double and multple complex roots are dealt wth n the standard manner (multplyng the usual basc soluton by powers of n). So s the constructon of a specfc soluton. 4 Non-homogeneous case We now consder dfference equatons wth a non-zero rght hand sde (.e. havng at least one term not multpled by an a ; such terms can appear on ether sde of the equaton, but we wll always transfer them to the RHS). In general, the RHS can be any expresson nvolvng and a handful of constant parameters, but we wll explctly treat only the case of the RHS havng the form of P () θ (10) 6
7 where P () s a polynomal n, and θ s a specfc number. One can show that a general soluton to a dfference equaton wth a RHS (10) can be bult as follows: Take the general soluton to the correspondng homogeneous equaton, and add to t a so called partcular soluton to the complete (non-homogeneous) equaton. Note that the partcular soluton wll have no arbtrary constants (the old A, B, C, etc.) - the frst part of the soluton takes care of these. To fnd a partcular soluton s, for our partcular RHS (10), rather easy (at least n prncple); all we need to know s that t wll have the followng form: Q() θ m (11) where Q() s a polynomal of the same degree as P (), but wth undetermned (meanng: yet-to-be-determned) coeffcents, and m s the multplcty of θ as a root of the (homogeneous) characterstc polynomal. Typcally, θ s not a root of the characterstc polynomal, whch means that m =0and the last factor of (11) dsappears. But then, f θ happens to be a smple root, m =1;f t s a double root, m =, etc. 4.1 Specal case of θ =1 In ths case, the RHS s smply a polynomal n. One has to realze that, mplctly, t s stll multpled by 1. Thus, to establsh the value of m, one has to check whether 1 s a root of the characterstc polynomal (and ts multplcty). As an example, let us solve Eq. 1. Its characterstc polynomal has the followng roots: 5 and 1. The RHS of (1) s a lnear polynomal, whch means that θ = 1, whose multplcty as a root of the characterstc polynomal s m =1. The form of the partcular soluton s then a part. =(q 0 + q 1 ) (1) Note that we need a complete lnear polynomal (.e. a polynomal wth both lnear and constant coeffcents), even though P () had only a lnear term. Also,weneedtorealzethatq 0 and q 1 (unlke A and B) arenot arbtrary numbers; to fnd the actual partcular soluton, we must frst substtute (1) nto (1) - ths can be done effcently by Maple - and then solve for the correct values of q 0 and q 1 by matchng the coeffcents of all powers of. Substtutng (1) nto the LHS of (1) yelds: 7q 0 +11q 1 +14q 1 To match ths to the RHS of (1), namely to 7, we need (matchng lnear terms): 14q 1 =7 whch mmedately yelds q 1 = 1, and (matchng absolute terms): 7q 0 +11q 1 =0 7
8 whch mples that q 0 = The partcular soluton s thus a part. = whch leads to the followng fully general soluton of (1): a gen. = A ( 5 ) + B Note that only at ths pont we would be n a poston to deal wth ether ntal or boundary condtons (a common mstake s to omt the partcular soluton when fndng A and B). For example, gven that a 0 =3and a 1 =6, we get whch yelds A = More examples A + B = 3 5 A + B = and B = 49, and the followng specfc soluton: a spec. = ( 5 ) In ths secton we concentrate on buldng partcular solutons only (extendng them to general soluton s qute trval). Another frequent specal case of (10) arses when the polynomal P () s just a constant (of zero degree), e.g. a +5 3a +4 4a a + +84a a =5 3 (13) To fnd the correct form of a partcular soluton, we only need to know the multplcty of θ =3as a root of the characterstc polynomal λ 5 3λ 4 4λ 3 +13λ +84λ +36 (14) Ths can be fgured out wthout solvng the correspondng equaton! Frst, we substtute 3 for λ n (14). Snce we get zero, 3 s clearly a root. To fnd ts multplcty, we repeatedly dfferentate (14) wth respect to λ, then make the same substtuton, tll we reach a non-zero value. The order of the frst non-zero dervatve yelds the multplcty of the root. In the case of our example, the frst dervatve stll evaluates to zero, but the second one s equal to 350. Ths mples that the correspondng multplcty m s equal to. The partcular soluton wll thus have the followng form: a = q 3 (15) 8
9 where the correct value of q s yet to be found. Substtutng (15) nto the LHS of (13), best done wth the help of Maple, yelds 3150q 3 To make ths equal to the RHS of (9), q = 5 The partcular soluton s then a part. = = For our fnal example, we just replace the RHS of (13) by 3+ ( ) (16).e. a quadratc polynomal, multpled by θ, where θ = 1. Frst we need to know the multplcty of 1 as a root of (14) - the characterstc polynomal has not changed! The correspondng substtuton yelds 0 for the polynomal tself, and for ts frst dervatve. The correspondng m s thus equal to 1. The partcular soluton wll have a form of q 0 + q 1 + q ( ) When substtuted nto the LFS of (13), ths yelds (the amount of algebra s formdable - Maple now becomes ndspensable): µ 441q0 1105q q 16 88q q q 16 1 ( ) We can make the last expresson equal to (16) by makng q = 16 q 1 = = and q 0 = q 1 73q 441 = The partcular soluton thus reads: a part. = (wehavecancelledoutthefactorof 8). 4.3 Superposton prncple ( ) 3 133, Fnally, when the RHS of a non-homegeneous equaton s a sum of two or more terms of type (10), t s easy to show that the correspondng partcular soluton wll be the correspondng sum (superposton) of the(two or more)ndvdual partcular solutons, constructed separately for each term (gnorng the rest). Thus, replacng the RHS of our prevous example (yet one more tme) by 3+ 5 (17) 9
10 we fst have to realze that ths does not have the form of (10), but can be wrttenasassumof 3+ (18) and 5 (19) whch, ndvdually, we do know how to deal wth. The frst of these s just a quadratc polynomal (mplctly, θ =1), and we can also readly verfy that 1 s not a root of (14), mplyng that m =0. The frst partcular soluton wll thus have the form of q 0 + q 1 + q Ths, substtuted nto the LHS of (13), yelds: 108q + (108q 1 +7q ) + (108q 0 +36q 1 78q ) mplyng that q = 1 108,q 1 = = 1 16 and q 0 = 3 36q 1+78q The correspondng partcular soluton s a part(1) = = Smlarly, to fnd a partcular soluton to match (19), we frst have to fnd the multplcty of θ =as a root of (14). Substtutng for λ agan results n a non-zero value, agan mplyng that m =0. The form of the second partcular soluton s thus smply q Ths, substtuted nto the LHS of (13), results n yeldng q = Thus, a part() 80q = 4 The partcular soluton whch solves the whole equaton s thus a superposton (sum) of the two, namely: a part. = To buld the fully general soluton, we frst need to fnd all roots of (14). We already know that 3 s a double root and 1 a smple root; t s not dffcult to verfy that s also a double root - note that replacng n (17) by ( ) would have made the prevous example a lot more dffcult! 10
11 Ths mples that a general = A 3 +B 3 +C ( 1 ) +D ( ) +E ( ) Only at ths pont one could start constructng a specfc soluton, based on fve gven values of a. Ths would lead to a lnear set of ordnary equatons for A, B, C, D and E (farly routne). Hopefully, no explct example s necessary. 11
1 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 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 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 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 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 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 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 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 informationThe KMO Method for Solving Non-homogenous, m th Order Differential Equations
The KMO Method for Solvng Non-homogenous, m th Order Dfferental Equatons Davd Krohn Danel Marño-Johnson John Paul Ouyang March 14, 2013 Abstract Ths paper shows a smple tabular procedure for fndng the
More information= z 20 z n. (k 20) + 4 z k = 4
Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5
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 informationSection 3.6 Complex Zeros
04 Chapter Secton 6 Comple Zeros When fndng the zeros of polynomals, at some pont you're faced wth the problem Whle there are clearly no real numbers that are solutons to ths equaton, leavng thngs there
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More information12. The Hamilton-Jacobi Equation Michael Fowler
1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
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 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 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 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 informationExample: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,
The greatest common dvsor of two ntegers a and b (not both zero) s the largest nteger whch s a common factor of both a and b. We denote ths number by gcd(a, b), or smply (a, b) when there s no confuson
More information9 Characteristic classes
THEODORE VORONOV DIFFERENTIAL GEOMETRY. Sprng 2009 [under constructon] 9 Characterstc classes 9.1 The frst Chern class of a lne bundle Consder a complex vector bundle E B of rank p. We shall construct
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 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 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 informationFundamental loop-current method using virtual voltage sources technique for special cases
Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,
More informationSolution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method
Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another well-known teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222
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 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 information8.6 The Complex Number System
8.6 The Complex Number System Earler n the chapter, we mentoned that we cannot have a negatve under a square root, snce the square of any postve or negatve number s always postve. In ths secton we want
More informationLinear Regression Analysis: Terminology and Notation
ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented
More informationa b a In case b 0, a being divisible by b is the same as to say that
Secton 6.2 Dvsblty among the ntegers An nteger a ε s dvsble by b ε f there s an nteger c ε such that a = bc. Note that s dvsble by any nteger b, snce = b. On the other hand, a s dvsble by only f a = :
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 informationAffine transformations and convexity
Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton
More informationTHE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens
THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of
More informationWeek 5: Neural Networks
Week 5: Neural Networks Instructor: Sergey Levne Neural Networks Summary In the prevous lecture, we saw how we can construct neural networks by extendng logstc regresson. Neural networks consst of multple
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 informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
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 informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationTHE SUMMATION NOTATION Ʃ
Sngle Subscrpt otaton THE SUMMATIO OTATIO Ʃ Most of the calculatons we perform n statstcs are repettve operatons on lsts of numbers. For example, we compute the sum of a set of numbers, or the sum of the
More informationCHAPTER 5: Lie Differentiation and Angular Momentum
CHAPTER 5: Le Dfferentaton and Angular Momentum Jose G. Vargas 1 Le dfferentaton Kähler s theory of angular momentum s a specalzaton of hs approach to Le dfferentaton. We could deal wth the former drectly,
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 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 informationSELECTED PROOFS. DeMorgan s formulas: The first one is clear from Venn diagram, or the following truth table:
SELECTED PROOFS DeMorgan s formulas: The frst one s clear from Venn dagram, or the followng truth table: A B A B A B Ā B Ā B T T T F F F F T F T F F T F F T T F T F F F F F T T T T The second one can be
More informationCOMPLEX NUMBERS AND QUADRATIC EQUATIONS
COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not
More informationCausal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim
Causal Damonds M. Aghl, L. Bombell, B. Plgrm Introducton The correcton to volume of a causal nterval due to curvature of spacetme has been done by Myrhem [] and recently by Gbbons & Solodukhn [] and later
More informationMAE140 - Linear Circuits - Winter 16 Final, March 16, 2016
ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have
More informationarxiv: v1 [math.ho] 18 May 2008
Recurrence Formulas for Fbonacc Sums Adlson J. V. Brandão, João L. Martns 2 arxv:0805.2707v [math.ho] 8 May 2008 Abstract. In ths artcle we present a new recurrence formula for a fnte sum nvolvng the Fbonacc
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 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 informationDeterminants Containing Powers of Generalized Fibonacci Numbers
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 19 (2016), Artcle 1671 Determnants Contanng Powers of Generalzed Fbonacc Numbers Aram Tangboonduangjt and Thotsaporn Thanatpanonda Mahdol Unversty Internatonal
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 informationDISCRIMINANTS AND RAMIFIED PRIMES. 1. Introduction A prime number p is said to be ramified in a number field K if the prime ideal factorization
DISCRIMINANTS AND RAMIFIED PRIMES KEITH CONRAD 1. Introducton A prme number p s sad to be ramfed n a number feld K f the prme deal factorzaton (1.1) (p) = po K = p e 1 1 peg g has some e greater than 1.
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 informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More information1 Generating functions, continued
Generatng functons, contnued. Generatng functons and parttons We can make use of generatng functons to answer some questons a bt more restrctve than we ve done so far: Queston : Fnd a generatng functon
More informationand problem sheet 2
-8 and 5-5 problem sheet Solutons to the followng seven exercses and optonal bonus problem are to be submtted through gradescope by :0PM on Wednesday th September 08. There are also some practce problems,
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 informationSolutions to Problem Set 6
Solutons to Problem Set 6 Problem 6. (Resdue theory) a) Problem 4.7.7 Boas. n ths problem we wll solve ths ntegral: x sn x x + 4x + 5 dx: To solve ths usng the resdue theorem, we study ths complex ntegral:
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 informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationIterative General Dynamic Model for Serial-Link Manipulators
EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 1. Introducton Iteratve General Dynamc Model for Seral-Lnk Manpulators In ths set of notes, we are gong to develop a method for computng a general
More informationBallot Paths Avoiding Depth Zero Patterns
Ballot Paths Avodng Depth Zero Patterns Henrch Nederhausen and Shaun Sullvan Florda Atlantc Unversty, Boca Raton, Florda nederha@fauedu, ssull21@fauedu 1 Introducton In a paper by Sapounaks, Tasoulas,
More informationModule 2. Random Processes. Version 2 ECE IIT, Kharagpur
Module Random Processes Lesson 6 Functons of Random Varables After readng ths lesson, ou wll learn about cdf of functon of a random varable. Formula for determnng the pdf of a random varable. Let, X be
More informationMath 261 Exercise sheet 2
Math 261 Exercse sheet 2 http://staff.aub.edu.lb/~nm116/teachng/2017/math261/ndex.html Verson: September 25, 2017 Answers are due for Monday 25 September, 11AM. The use of calculators s allowed. Exercse
More informationDigital Signal Processing
Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over
More informationA NOTE ON CES FUNCTIONS Drago Bergholt, BI Norwegian Business School 2011
A NOTE ON CES FUNCTIONS Drago Bergholt, BI Norwegan Busness School 2011 Functons featurng constant elastcty of substtuton CES are wdely used n appled economcs and fnance. In ths note, I do two thngs. Frst,
More informationFACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP
C O L L O Q U I U M M A T H E M A T I C U M VOL. 80 1999 NO. 1 FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP BY FLORIAN K A I N R A T H (GRAZ) Abstract. Let H be a Krull monod wth nfnte class
More informationCOS 521: Advanced Algorithms Game Theory and Linear Programming
COS 521: Advanced Algorthms Game Theory and Lnear Programmng Moses Charkar February 27, 2013 In these notes, we ntroduce some basc concepts n game theory and lnear programmng (LP). We show a connecton
More informationApproximate Smallest Enclosing Balls
Chapter 5 Approxmate Smallest Enclosng Balls 5. Boundng Volumes A boundng volume for a set S R d s a superset of S wth a smple shape, for example a box, a ball, or an ellpsod. Fgure 5.: Boundng boxes Q(P
More information1 Generating functions, continued
Generatng functons, contnued. Exponental generatng functons and set-parttons At ths pont, we ve come up wth good generatng-functon dscussons based on 3 of the 4 rows of our twelvefold way. Wll our nteger-partton
More information5 The Rational Canonical Form
5 The Ratonal Canoncal Form Here p s a monc rreducble factor of the mnmum polynomal m T and s not necessarly of degree one Let F p denote the feld constructed earler n the course, consstng of all matrces
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 informationLecture Space-Bounded Derandomization
Notes on Complexty Theory Last updated: October, 2008 Jonathan Katz Lecture Space-Bounded Derandomzaton 1 Space-Bounded Derandomzaton We now dscuss derandomzaton of space-bounded algorthms. Here non-trval
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 informationPRIMES 2015 reading project: Problem set #3
PRIMES 2015 readng project: Problem set #3 page 1 PRIMES 2015 readng project: Problem set #3 posted 31 May 2015, to be submtted around 15 June 2015 Darj Grnberg The purpose of ths problem set s to replace
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 informationHMMT February 2016 February 20, 2016
HMMT February 016 February 0, 016 Combnatorcs 1. For postve ntegers n, let S n be the set of ntegers x such that n dstnct lnes, no three concurrent, can dvde a plane nto x regons (for example, S = {3,
More informationComplex Numbers. x = B B 2 4AC 2A. or x = x = 2 ± 4 4 (1) (5) 2 (1)
Complex Numbers If you have not yet encountered complex numbers, you wll soon do so n the process of solvng quadratc equatons. The general quadratc equaton Ax + Bx + C 0 has solutons x B + B 4AC A For
More informationPolynomials. 1 What is a polynomial? John Stalker
Polynomals John Stalker What s a polynomal? If you thnk you already know what a polynomal s then skp ths secton. Just be aware that I consstently wrte thngs lke p = c z j =0 nstead of p(z) = c z. =0 You
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
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 informationGraph Reconstruction by Permutations
Graph Reconstructon by Permutatons Perre Ille and Wllam Kocay* Insttut de Mathémathques de Lumny CNRS UMR 6206 163 avenue de Lumny, Case 907 13288 Marselle Cedex 9, France e-mal: lle@ml.unv-mrs.fr Computer
More information18.781: Solution to Practice Questions for Final Exam
18.781: Soluton to Practce Questons for Fnal Exam 1. Fnd three solutons n postve ntegers of x 6y = 1 by frst calculatng the contnued fracton expanson of 6. Soluton: We have 1 6=[, ] 6 6+ =[, ] 1 =[,, ]=[,,
More informationMA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials
MA 323 Geometrc Modellng Course Notes: Day 13 Bezer Curves & Bernsten Polynomals Davd L. Fnn Over the past few days, we have looked at de Casteljau s algorthm for generatng a polynomal curve, and we have
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 informationSubset Topological Spaces and Kakutani s Theorem
MOD Natural Neutrosophc Subset Topologcal Spaces and Kakutan s Theorem W. B. Vasantha Kandasamy lanthenral K Florentn Smarandache 1 Copyrght 1 by EuropaNova ASBL and the Authors Ths book can be ordered
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 informationRadar Trackers. Study Guide. All chapters, problems, examples and page numbers refer to Applied Optimal Estimation, A. Gelb, Ed.
Radar rackers Study Gude All chapters, problems, examples and page numbers refer to Appled Optmal Estmaton, A. Gelb, Ed. Chapter Example.0- Problem Statement wo sensors Each has a sngle nose measurement
More informationψ = i c i u i c i a i b i u i = i b 0 0 b 0 0
Quantum Mechancs, Advanced Course FMFN/FYSN7 Solutons Sheet Soluton. Lets denote the two operators by  and ˆB, the set of egenstates by { u }, and the egenvalues as  u = a u and ˆB u = b u. Snce the
More informationRemarks on the Properties of a Quasi-Fibonacci-like Polynomial Sequence
Remarks on the Propertes of a Quas-Fbonacc-lke Polynomal Sequence Brce Merwne LIU Brooklyn Ilan Wenschelbaum Wesleyan Unversty Abstract Consder the Quas-Fbonacc-lke Polynomal Sequence gven by F 0 = 1,
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 informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
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 informationLecture 2: Gram-Schmidt Vectors and the LLL Algorithm
NYU, Fall 2016 Lattces Mn Course Lecture 2: Gram-Schmdt Vectors and the LLL Algorthm Lecturer: Noah Stephens-Davdowtz 2.1 The Shortest Vector Problem In our last lecture, we consdered short solutons to
More informationRandom Walks on Digraphs
Random Walks on Dgraphs J. J. P. Veerman October 23, 27 Introducton Let V = {, n} be a vertex set and S a non-negatve row-stochastc matrx (.e. rows sum to ). V and S defne a dgraph G = G(V, S) and a drected
More informationProblem Solving in Math (Math 43900) Fall 2013
Problem Solvng n Math (Math 43900) Fall 2013 Week four (September 17) solutons Instructor: Davd Galvn 1. Let a and b be two nteger for whch a b s dvsble by 3. Prove that a 3 b 3 s dvsble by 9. Soluton:
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationSL n (F ) Equals its Own Derived Group
Internatonal Journal of Algebra, Vol. 2, 2008, no. 12, 585-594 SL n (F ) Equals ts Own Derved Group Jorge Macel BMCC-The Cty Unversty of New York, CUNY 199 Chambers street, New York, NY 10007, USA macel@cms.nyu.edu
More information