Chapter 7 Singular Value Decomposition
|
|
- Rosaline Bridges
- 6 years ago
- Views:
Transcription
1 EE448/58 Vesion. John Stensy Chapte 7 Singla Vale Decomposition This chapte coves the singla vale decomposition (SVD) of an m n matix. This algoithm has applications in many aeas inclding nmeical analysis and signal/image pocessing. Seveal SVD applications ae discssed in Chapte 8. It is not ncommon to find in the nmeical linea algea liteate statements that poclaim the SVD algoithm as one of the most impotant diagnostic tools availale fo the analysis of linea algeaic systems. Theoem 7- (SVD) et A C m n e an m n complex-valed matix with ank. Then thee exists nitay matices U C m m and V C n n sch that U * AV R m n, (7-) whee... > ae called the singla vales of the matix A (ememe that > is the lagest and > is the smallest and these singla vales appea on the diagonal of a diagonal matix). oof: Given A C m n of ank, the matix A * A is n n Hemitian and (at least) non-negative definite (we have aged this efoe). As we know, the eigenvales of A * A ae eal-valed; of them ae positive, and n- of them ae zeo. Denote these eigenvales as... > (with sqae-oots denoted as... > ). et v, v,..., v, v +,..., v n denote the n othonomal n eigenvectos of A * A. Then, we can wite CH7.DC age 7-
2 EE448/58 Vesion. John Stensy A A vi R S T i vi, i,, + i n (7-) whee v j, vi vi v j, i j i j. (7-3) ote that v i, i, is an othonomal asis of R(A*), and v i, + i n, is an othonomal asis of K(A). ow, define the vectos i Av i/ i, i, (7-4) whee i i. ote that the vectos i, i, ae nit length m vectos since v A Av v v i Av i i ii i i i / i, (7-5) i i i and they ae othonomal. ow, se Gam-Schmidt (o any othe pocede) to fill ot the i and fom a set of m othonomal m vectos,,...,, +,..., m with i j, i j i j, (7-6) i, j m. ote that i, i, is an othonomal asis of R(A), and i, + i m, is an othonomal asis of K(A*). ow, define the two matices CH7.DC age 7-
3 EE448/58 Vesion. John Stensy U V m v v v n. (7-7) Unitay matices U and V ae m m and n n, espectively (U and V ae nitay: U * U - and V * V - ). ow, se nitay U and V to compte the m n matix Σ U AV m A v v v n Av Av Av n Av Av Av n mav mav mav n. (7-8) ow, let s ij Av i j, i m, j n, e the elements of m n matix Σ. We age that only the fist diagonal elements of Σ ae non-zeo. Fist, conside the elements s ij, i, j, in the ppe-left lock; with the aid of (7-), these elements ae compted to e vi A Av j vi j v j i, i j s ij i A v j S i i i j R T, i, j. (7-9) The second lock to e consideed ae those elements s ij, i, + j n, in the ppe-ight {n-} lock; with the aid of (7-), these elements ae compted to e vi A Av j v s i ij i A v j v j i + j n i i c h,,. (7-) The thid lock to e consideed ae those elements s ij, + i m, j, in the lowe-left {m} lock; with the aid of (7-4), these elements ae compted to e CH7.DC age 7-3
4 EE448/58 Vesion. John Stensy s ij i A v j i ( j j ), + i m, j. (7-) The foth, and last, lock to e consideed ae those elements s ij, + i m, + j n, in the lowe-ight {m-} {n-} lock. ow, fo i, the vectos Av i ae non-zeo and othogonal; this implies they ae independent. Since ank[a], we mst have Av j, + j n, dependent on Av k, k. Hence, de to the othonomal nate of the j, the elements in the lowe-ight {m-} {n-} lock ae s ij i A v j i a linea comination of Av k, k i a linea comination of k, k g g, + i m, + j n. (7-) The te nate of m n matix Σ is clea. Fom (7-9) thogh (7-), we see that Σ U AV, (7-3) a matix whose only non-zeo enties ae the fist diagonal elements. ata's SVD Fnction Given an m n matix A, ata can calclate the SVD of A. The ata syntax is [U S V] svd(a). (7-4) CH7.DC age 7-4
5 EE448/58 Vesion. John Stensy Example % Ente matix A > A [6 3;4 ; ] A % Fist we calclate the svd sing the asic theoy > B A'*A B % Calclate eignevales and eigenvectos of A T A. > [V D] eig(b) V D 7 % V is an othogonal matix sed in svd pocede % Calclate U sed in svd. Use othonomal asis of K[A T ] to %"flsh ot" U (ecall Fige 5. of the class notes) > U [A*V(:,)/sqt(D(,)) nll(a')] U % Use matices U and V to tansfom A to svd fom >U'*A*V ans % Boy, that was a lot of wok! et's do it the easy way - let % the ata svd fnction do the wok! >[U S V] svd(a) U S V % ote that S S. ata's svd fnction podced the same % decomposition (howeve, y sing diffeent U and V) CH7.DC age 7-5
6 EE448/58 Vesion. John Stensy Bases fo ange[a], ke[a*], ange[a*] and ke[a] Fom SVD In the poof of Theoem 7-, we sed the fact that,,..., ae independent and span ange[a]. ow, since k, k m, ae othonomal, the vectos +, +,..., m ae othogonal to ange[a]; that is, they span ke[a*] (ecall Fige 5.). Hence, the m m matix U can e patitioned as U [ ] m span ange[a] span ke[a ] (7-5) In the poof of Theoem 7-, we sed the fact that Av j, + j n. Fom this, it is easy to see that v +, v +,..., v n span ke[a]. This implies that v, v,..., v spans ange[a*]. Hence, the n n matix V e patitioned as V [ v ] v 44 v 443 v + v v n (7-6) span ange[a ] span ke[a] SVD Related ojection peatos The SVD povides a method of finding the othogonal pojection onto the ange and kenel of any m n matix A. Fom (7-5) and (7-6), we conclde that RA thogonal ojection on Range[A] (7-7) KA thogonal ojection on Ke[A] v + v + v v + v + v n n (7-8) CH7.DC age 7-6
7 EE448/58 Vesion. John Stensy RA* thogonal ojection on Range[A ] v v v v v v (7-9) KA* thogonal ojection on Ke[A ] m m (7-) Singla Vales of A - ; -om of A - Sppose A UΣV * is invetale. Then all singla vales k >, k n. Also, we know that A - VΣ U *, (7-) whee Σ / / / n. (7-) Clealy, / k, k n, ae the singla vales of A -. Since n is the smallest singla vale of A, / n is the lagest singla vale of A -. As a eslt, we can wite A / n. (7-3) SVD Expansion of atix A n n, we have et A e an m n matix of ank. Since A UΣV *, whee U m m, Σ is m n and V CH7.DC age 7-7
8 EE448/58 Vesion. John Stensy A m so that v v v n (7-4) A i i v i. (7-5) i ote that m n matix A can e epesented y i, v i and i, fo i. If matix A is lage, t has small ank, then matix A can e moe efficiently stoed y saving the i, v i and i, i. This techniqe has een sed in image compession and othe aeas that deal with lage dimensional, low-ank, matices. meical Rank of a atix y SVD Analysis The nmeical detemination of the ank of m n matix A is not tivial. Rank detemination can e complicated y impecise epesentation of A (i.e., fzzy data ) and onding eos that occ in all pactical calclations. Given small ε >, the ε-ank of matix A is defined as ank( A, ε) min ank( B). (7-6) AB ε That is, look at all m n matices B that ae within ε of A. Compte the ank of each εneigho of A. The smallest of these anks is ank(a,ε). Conside an application whee m n matix A {a ij } is detemined expeimentally with each a ij coect to within ±.. Then it might make sense to detemine ank(a,.). If m n matix A contains floating point nmes, then it is easonale to egad A as nmeically ank CH7.DC age 7-8
9 EE448/58 Vesion. John Stensy deficient if ank(a,ε) < min(m,n), whee ε max( m, n) A eps. The ε-ank of a matix can e detemined, in most cases, y SVD analysis. In the nmeical analysis commnity, SVD is geneally egaded as the most eliale method fo the detemination of nmeical ank. The next theoem shows how SVD is sed to detemine nmeical ank. Theoem 7- et the SVD of m n matix A e A UΣV*, whee U is m m, V is n n, and U, V ae nitay. et k e an intege, and define k Ak i i v i. (7-7) i Then we have min ank( B) k A B A A k k +. (7-8) atix A k is a ank k matix that is k+ in distance fom A, and no othe ank k matix is close to A then is A k. Eqivalently, matix A k comes as close to A as does any ank k matix. oof: Since U * A k V diag(,,..., k,,..., ), it follows that ank(a k ) k and U (A - A k ) V diag(, 4,... 4, 3, k+,..., ρ), (7-9) k zeos whee ρ min(m, n). Also, it follows that A A k k +, since the -nom is invaiant to nitay tansfomations. ow, let B any m n matix of ank k. Then (n - k) is the nllity of B, and we can find an othonomal set x, x,..., x n-k of n vectos fo which CH7.DC age 7-9
10 EE448/58 Vesion. John Stensy ke[b] span{x, x,..., x n-k}. (7-3) In an n-dimensional space, thee ae at most n linealy independent vectos so that span{x, x,..., x n-k} span{v, v,..., v k+}. (7-3) et Z span{x, x,..., x n-k} span{v, v,..., v k+} and Z. ote that F k A Z i i vi Z i vi Z i i vi Z i H G I + i K J ( ) + ( ), (7-3) i i k+ assming that k+ (othewise omit the second sm on the ight-hand side of (7-3)). ow, since Z {v k+,..., v }, the second sm (on the ight-hand side of (7-3)) is zeo and k+ A Z i ( vi Z) i. (7-33) i Finally, since BZ (ecall that Z span{x, x,..., x n-k} K(B)), we have k+ A B ( A B) Z AZ i vi Z ( ) k+. (7-34) i So a ank k matix can e no close to A than k+, and we have finished o poof. et m n matix A have fll ank; that is, ank[a] min(n,m). Theoem 7. says that the smallest singla vale of matix A is the -nom distance of A to the set of all m n, ank deficient, matices. We can calclate the ε-ank of A fom knowledge of the singla vales. If ε ank(a,ε), CH7.DC age 7-
11 EE448/58 Vesion. John Stensy then we have ε +, (7-35) ε ε ρ whee ρ min(m, n). So, to detemine ε ank, place ε in the odeing of singla vales, and ead ε as the index of the fist singla vale lage than ε. SVD and the inimm om, east Sqaes olem Conside again the geneal AX polem ( may, o may not e, in ange[a]). Hee, we want to find the minimm nom X that minimizes AX. SVD analysis can solve this polem. Theoem 7-3 Sppose U * AV Σ is the SVD of m n matix A, whee ank[a]. et (7-7) descie colmn patitions of matices U and V. Then fo any the minimm nom X that minimizes AX is given y F I X i v i, (7-36) i HG i KJ whee i, i, ae the non-zeo singla vales of A. In addition, the minimm eo satisfies m AX i i + d i. (7-37) oof: Fo any X we have CH7.DC age 7-
12 EE448/58 Vesion. John Stensy CH7.DC age 7- AX AX A X m m U U U VV U U Σα, m m + +, m m i i i i i,, d i d i i m +, (7-38) whee α [α α... α n ] T V * X. Clealy, the "optimm" X that minimizes AX eqies that α i i i /, i. If we set αi fo + i n, then the eslting X has the minimm -nom. With this α, the optimm X ecomes
13 EE448/58 Vesion. John Stensy CH7.DC age 7-3 X v v v n i i i i F HG I KJ Vα / /, (7-39) and the smallest eo satisfies (7-37). sedo Invese Fom SVD ote that (7-36) can e witten as / / / m X v v v v i i i i n ( / ) Fom this epesentation, we can conclde that A + V Σ + U *, (7-4) / / / Σ +. ote that Σ + is otained fom Σ y simply inveting the non-zeo singla vales. Eqations (7-4) (7-4)
14 EE448/58 Vesion. John Stensy (7-4) and (7-4) ae sed y ata to calclate the psedo invese A +. Comptational Sensitivity of the sedo Invese The psedo invese A + can e difficlt to compte accately. It can e vey sensitive to ond-off eos (and othe small changes in A); A + does not always depend continosly on the enties of A. This fact can e seen easily fom inspection of (7-4) and (7-4). Fo some A with vey small non-zeo singla vales, small changes in the elements of A can change the ank of A and the nme of non-zeo singla vales. This can have a damatic inflence on the compted A +. Example Conside the matix A and the petation/eo δa given y A, δa ε. It is easy to compte A +, (A + δa) + / ε. ote that a small petation matix δa (i.e., small ε) added to A can case an extemely lage change in the psedo invese (A + δa) +. In pactice, we mst compae the small singla vales of A with some small toleance that eflects the eos in the data and/o comptation. We shold set to zeo those singla vales of A that fall elow some pescied toleance. This zeoes ot the coesponding /(singla vale) tems fond in (7-4). This is the appoach taken y ata's pinv fnction. The pinv fnction accepts a theshold vale as an agment (the theshold is called tol) when it ses (7-4) and (7-4) to CH7.DC age 7-4
15 EE448/58 Vesion. John Stensy calclate A +. Fnction pinv(a,tol) sets to zeo those singla vales that ae smalle than tol. Example Conside the following ata scipt % Ente atix A and petation da A [ ; ; ] A da [ ; E-6; ] da E-6 % Calclate pinv(a+da, E-5), se tol E-5 pinv(a+da,e-5) ans %Calclate pinv(a+da, E-6), se tol E-6 pinv(a+da,e-6) ans %Calclate pinv(a+da, E-7), se tol E-7 pinv(a+da,e-7) ans % Wow! What a change! Image Compession and SVD In some applications an image can e epesented as an m n matix A of eal nmes. ow, the matix A is "highly singla", and it contains a lot of edndant infomation (think of all the "white space" occpied y "nothing" in a typical image). Using SVD, it can e epesented in the seies min( m, n) A i i v i. (7-43) i ften, only a few of the i ae significant, and the expansion can e tncated to k tems with only some loss of claity. CH7.DC age 7-5
16 EE448/58 Vesion. John Stensy Fo example, conside the it map image illstated y Fig. 7-a. This image eqies ,536 eal nmes to epesent. The SVD of this image eveals singla vales that fall off apidly. Fige 7- is a log-scale plot of the singla vales of this image. As it tns ot 86 /.5 -. By setting k, 87 k 56, we etain only those singla vales that ae within 99.5% of (we tncate the smallest.5%). Bt we edce the stoage eqiements to ,8 eal Singla Vales k Fige 7-: og plot of singla vales. nmes; this podces a edction of appoximately 33% in stoage eqiements. Howeve, not mch was lost in the compessed image, as can e seen fom Fig. 7-. By setting k, 5 k 56, we edce the stoage eqiements to ,65 eal nmes; this podces a edction of appoximately 6% in stoage eqiements. Howeve, little was lost in the compessed image, as can e seen fom Fig. 7-c. By setting k, 6 k 56, we edce the stoage eqiements to ,85 eal nmes; this podces a edction of appoximately 8% in stoage eqiements. Howeve, now the degadation in the compessed image is ovios, as can e seen fom Fig. 7-d. By setting k, 6 k 56, we edce the stoage eqiements to ,695 eal nmes; this podces a edction of appoximately 88% in stoage eqiements. Howeve, now the degadation in the compessed image is vey ovios, as can e seen fom Fig. 7-e. Finally, y setting k, 6 k 56, we edce the stoage eqiements to ,565 eal nmes; this podces a edction of appoximately 96% in stoage eqiements. Howeve, now the image is not ecognizale, as can e seen fom Fig. 7-f. CH7.DC age 7-6
17 EE448/58 Vesion. John Stensy Fige 7-a: iginal image. Fige 7-: Expansion tncated to 86 tems. Fige 7-c: Expansion tncated to 5 tems. Fige 7-d: Expansion tncated to 5 tems. Fige 7-e: Expansion tncated to 5 tems. Fige 7-f: Expansion tncated to 5 tems. CH7.DC age 7-7
Chapter 5 Linear Equations: Basic Theory and Practice
Chapte 5 inea Equations: Basic Theoy and actice In this chapte and the next, we ae inteested in the linea algebaic equation AX = b, (5-1) whee A is an m n matix, X is an n 1 vecto to be solved fo, and
More informationPROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.
POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and
More informationA Crash Course in (2 2) Matrices
A Cash Couse in ( ) Matices Seveal weeks woth of matix algeba in an hou (Relax, we will only stuy the simplest case, that of matices) Review topics: What is a matix (pl matices)? A matix is a ectangula
More informationMAC Module 12 Eigenvalues and Eigenvectors
MAC 23 Module 2 Eigenvalues and Eigenvectos Leaning Objectives Upon completing this module, you should be able to:. Solve the eigenvalue poblem by finding the eigenvalues and the coesponding eigenvectos
More informationA STUDY OF HAMMING CODES AS ERROR CORRECTING CODES
AGU Intenational Jounal of Science and Technology A STUDY OF HAMMING CODES AS ERROR CORRECTING CODES Ritu Ahuja Depatment of Mathematics Khalsa College fo Women, Civil Lines, Ludhiana-141001, Punjab, (India)
More information15.081J/6.251J Introduction to Mathematical Programming. Lecture 6: The Simplex Method II
15081J/6251J Intoduction to Mathematical Pogamming ectue 6: The Simplex Method II 1 Outline Revised Simplex method Slide 1 The full tableau implementation Anticycling 2 Revised Simplex Initial data: A,
More informationVector d is a linear vector function of vector d when the following relationships hold:
Appendix 4 Dyadic Analysis DEFINITION ecto d is a linea vecto function of vecto d when the following elationships hold: d x = a xxd x + a xy d y + a xz d z d y = a yxd x + a yy d y + a yz d z d z = a zxd
More information3.1 Random variables
3 Chapte III Random Vaiables 3 Random vaiables A sample space S may be difficult to descibe if the elements of S ae not numbes discuss how we can use a ule by which an element s of S may be associated
More informationworking pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50
woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,
More informationFigure 1. We will begin by deriving a very general expression before returning to Equations 1 and 2 to determine the specifics.
Deivation of the Laplacian in Spheical Coodinates fom Fist Pinciples. Fist, let me state that the inspiation to do this came fom David Giffiths Intodction to Electodynamics textbook Chapte 1, Section 4.
More informationMATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form
MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE ANDRAS VASY We conside second ode constant coefficient scala linea PDEs on R n. These have the fom Lu = f L = a ij xi xj + b i xi + c i whee a ij b i and
More information18.06 Problem Set 4 Solution
8.6 Poblem Set 4 Solution Total: points Section 3.5. Poblem 2: (Recommended) Find the lagest possible numbe of independent vectos among ) ) ) v = v 4 = v 5 = v 6 = v 2 = v 3 =. Solution (4 points): Since
More informationAPPENDIX. For the 2 lectures of Claude Cohen-Tannoudji on Atom-Atom Interactions in Ultracold Quantum Gases
APPENDIX Fo the lectues of Claude Cohen-Tannoudji on Atom-Atom Inteactions in Ultacold Quantum Gases Pupose of this Appendix Demonstate the othonomalization elation(ϕ ϕ = δ k k δ δ )k - The wave function
More informationIntegral Control via Bias Estimation
1 Integal Contol via Bias stimation Consie the sstem ẋ = A + B +, R n, R p, R m = C +, R q whee is an nknown constant vecto. It is possible to view as a step istbance: (t) = 0 1(t). (If in fact (t) vaies
More informationComplex Eigenvalues. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB
Pepaed by Vince Zaccone Fo ampus Leaning ssistance Sevices at USB omplex Numbes When solving fo the oots of a quadatic equation, eal solutions can not be found when the disciminant is negative. In these
More informationOn the Quasi-inverse of a Non-square Matrix: An Infinite Solution
Applied Mathematical Sciences, Vol 11, 2017, no 27, 1337-1351 HIKARI Ltd, wwwm-hikaicom https://doiog/1012988/ams20177273 On the Quasi-invese of a Non-squae Matix: An Infinite Solution Ruben D Codeo J
More informationStanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012
Stanfod Univesity CS59Q: Quantum Computing Handout 8 Luca Tevisan Octobe 8, 0 Lectue 8 In which we use the quantum Fouie tansfom to solve the peiod-finding poblem. The Peiod Finding Poblem Let f : {0,...,
More informationAs is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.
Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.
More informationScattering in Three Dimensions
Scatteing in Thee Dimensions Scatteing expeiments ae an impotant souce of infomation about quantum systems, anging in enegy fom vey low enegy chemical eactions to the highest possible enegies at the LHC.
More informationWeighted least-squares estimators of parametric functions of the regression coefficients under a general linear model
Ann Inst Stat Math (2010) 62:929 941 DOI 10.1007/s10463-008-0199-8 Weighted least-squaes estimatos of paametic functions of the egession coefficients unde a geneal linea model Yongge Tian Received: 9 Januay
More informationf(y) signal norms system gain bounded input bounded output (BIBO) stability For what G(s) and f( ) is the closed-loop system stable?
Lecte 5 Inpt otpt stabilit Cose Otline o How to make a cicle ot of the point + i, and diffeent was to sta awa fom it... Lecte -3 Lecte 4-6 Modelling and basic phenomena (lineaization, phase plane, limit
More informationGeometry of the homogeneous and isotropic spaces
Geomety of the homogeneous and isotopic spaces H. Sonoda Septembe 2000; last evised Octobe 2009 Abstact We summaize the aspects of the geomety of the homogeneous and isotopic spaces which ae most elevant
More informationChapter 3: Theory of Modular Arithmetic 38
Chapte 3: Theoy of Modula Aithmetic 38 Section D Chinese Remainde Theoem By the end of this section you will be able to pove the Chinese Remainde Theoem apply this theoem to solve simultaneous linea conguences
More informationCOORDINATE TRANSFORMATIONS - THE JACOBIAN DETERMINANT
COORDINATE TRANSFORMATIONS - THE JACOBIAN DETERMINANT Link to: phsicspages home page. To leave a comment o epot an eo, please use the auilia blog. Refeence: d Inveno, Ra, Intoducing Einstein s Relativit
More informationAuchmuty High School Mathematics Department Advanced Higher Notes Teacher Version
The Binomial Theoem Factoials Auchmuty High School Mathematics Depatment The calculations,, 6 etc. often appea in mathematics. They ae called factoials and have been given the notation n!. e.g. 6! 6!!!!!
More informationQuantum Fourier Transform
Chapte 5 Quantum Fouie Tansfom Many poblems in physics and mathematics ae solved by tansfoming a poblem into some othe poblem with a known solution. Some notable examples ae Laplace tansfom, Legende tansfom,
More informationChapter 2: Introduction to Implicit Equations
Habeman MTH 11 Section V: Paametic and Implicit Equations Chapte : Intoduction to Implicit Equations When we descibe cuves on the coodinate plane with algebaic equations, we can define the elationship
More informationBASIC ALGEBRA OF VECTORS
Fomulae Fo u Vecto Algeba By Mi Mohammed Abbas II PCMB 'A' Impotant Tems, Definitions & Fomulae 01 Vecto - Basic Intoduction: A quantity having magnitude as well as the diection is called vecto It is denoted
More information6 Matrix Concentration Bounds
6 Matix Concentation Bounds Concentation bounds ae inequalities that bound pobabilities of deviations by a andom vaiable fom some value, often its mean. Infomally, they show the pobability that a andom
More informationDymore User s Manual Two- and three dimensional dynamic inflow models
Dymoe Use s Manual Two- and thee dimensional dynamic inflow models Contents 1 Two-dimensional finite-state genealized dynamic wake theoy 1 Thee-dimensional finite-state genealized dynamic wake theoy 1
More informationJORDAN CANONICAL FORM AND ITS APPLICATIONS
JORDAN CANONICAL FORM AND ITS APPLICATIONS Shivani Gupta 1, Kaajot Kau 2 1,2 Matheatics Depatent, Khalsa College Fo Woen, Ludhiana (India) ABSTRACT This pape gives a basic notion to the Jodan canonical
More informationChem 453/544 Fall /08/03. Exam #1 Solutions
Chem 453/544 Fall 3 /8/3 Exam # Solutions. ( points) Use the genealized compessibility diagam povided on the last page to estimate ove what ange of pessues A at oom tempeatue confoms to the ideal gas law
More informationPHYS 705: Classical Mechanics. Small Oscillations
PHYS 705: Classical Mechanics Small Oscillations Fomulation of the Poblem Assumptions: V q - A consevative system with depending on position only - The tansfomation equation defining does not dep on time
More informationLecture 16 Root Systems and Root Lattices
1.745 Intoduction to Lie Algebas Novembe 1, 010 Lectue 16 Root Systems and Root Lattices Pof. Victo Kac Scibe: Michael Cossley Recall that a oot system is a pai (V, ), whee V is a finite dimensional Euclidean
More informationNew problems in universal algebraic geometry illustrated by boolean equations
New poblems in univesal algebaic geomety illustated by boolean equations axiv:1611.00152v2 [math.ra] 25 Nov 2016 Atem N. Shevlyakov Novembe 28, 2016 Abstact We discuss new poblems in univesal algebaic
More information. Using our polar coordinate conversions, we could write a
504 Chapte 8 Section 8.4.5 Dot Poduct Now that we can add, sutact, and scale vectos, you might e wondeing whethe we can multiply vectos. It tuns out thee ae two diffeent ways to multiply vectos, one which
More informationProblem Set #10 Math 471 Real Analysis Assignment: Chapter 8 #2, 3, 6, 8
Poblem Set #0 Math 47 Real Analysis Assignment: Chate 8 #2, 3, 6, 8 Clayton J. Lungstum Decembe, 202 xecise 8.2 Pove the convese of Hölde s inequality fo = and =. Show also that fo eal-valued f / L ),
More informationInformation Retrieval Advanced IR models. Luca Bondi
Advanced IR models Luca Bondi Advanced IR models 2 (LSI) Pobabilistic Latent Semantic Analysis (plsa) Vecto Space Model 3 Stating point: Vecto Space Model Documents and queies epesented as vectos in the
More informationMDS Array Codes with Optimal Rebuilding
0 IEEE Intenational Symposium on Infomation Theoy Poceedings MDS Aay Codes with Optimal Reuilding Itzhak Tamo, Zhiying Wang and Jehoshua Buck Electical Engineeing Depatment, Califonia Institute of Technology,
More informationRoot-Exponential Accuracy for Coarse
Root-Exponential Accuacy fo Coase 1 Quantization of Finite Fame Expansions Felix Kahme, Rayan Saab, and Rachel Wad Abstact In this note, we show that by quantizing the N-dimensional fame coefficients of
More informationLecture 04: HFK Propagation Physical Optics II (Optical Sciences 330) (Updated: Friday, April 29, 2005, 8:05 PM) W.J. Dallas
C:\Dallas\0_Couses\0_OpSci_330\0 Lectue Notes\04 HfkPopagation.doc: Page of 9 Lectue 04: HFK Popagation Physical Optics II (Optical Sciences 330) (Updated: Fiday, Apil 9, 005, 8:05 PM) W.J. Dallas The
More informationOn the reconstruction of the coronal magnetic field from coronal Hanle / Zeeman observations
On the econstction of the coonal magnetic field fom coonal Hanle / Zeeman obsevations M. Kama & B. Inheste Max-Planck Institte fo Sola System Reseach GERMANY Linda 005 Coonal Magnetic Field Magnetic field
More informationConjugate Gradient Methods. Michael Bader. Summer term 2012
Gadient Methods Outlines Pat I: Quadatic Foms and Steepest Descent Pat II: Gadients Pat III: Summe tem 2012 Pat I: Quadatic Foms and Steepest Descent Outlines Pat I: Quadatic Foms and Steepest Descent
More informationMultiple Criteria Secretary Problem: A New Approach
J. Stat. Appl. Po. 3, o., 9-38 (04 9 Jounal of Statistics Applications & Pobability An Intenational Jounal http://dx.doi.og/0.785/jsap/0303 Multiple Citeia Secetay Poblem: A ew Appoach Alaka Padhye, and
More informationChapter 12: Kinematics of a Particle 12.8 CURVILINEAR MOTION: CYLINDRICAL COMPONENTS. u of the polar coordinate system are also shown in
ME 01 DYNAMICS Chapte 1: Kinematics of a Paticle Chapte 1 Kinematics of a Paticle A. Bazone 1.8 CURVILINEAR MOTION: CYLINDRICAL COMPONENTS Pola Coodinates Pola coodinates ae paticlaly sitable fo solving
More informationCALCULUS II Vectors. Paul Dawkins
CALCULUS II Vectos Paul Dawkins Table of Contents Peface... ii Vectos... 3 Intoduction... 3 Vectos The Basics... 4 Vecto Aithmetic... 8 Dot Poduct... 13 Coss Poduct... 21 2007 Paul Dawkins i http://tutoial.math.lama.edu/tems.aspx
More informationMODULE 5a and 5b (Stewart, Sections 12.2, 12.3) INTRO: In MATH 1114 vectors were written either as rows (a1, a2,..., an) or as columns a 1 a. ...
MODULE 5a and 5b (Stewat, Sections 2.2, 2.3) INTRO: In MATH 4 vectos wee witten eithe as ows (a, a2,..., an) o as columns a a 2... a n and the set of all such vectos of fixed length n was called the vecto
More informationWhen two numbers are written as the product of their prime factors, they are in factored form.
10 1 Study Guide Pages 420 425 Factos Because 3 4 12, we say that 3 and 4 ae factos of 12. In othe wods, factos ae the numbes you multiply to get a poduct. Since 2 6 12, 2 and 6 ae also factos of 12. The
More informationInternet Appendix for A Bayesian Approach to Real Options: The Case of Distinguishing Between Temporary and Permanent Shocks
Intenet Appendix fo A Bayesian Appoach to Real Options: The Case of Distinguishing Between Tempoay and Pemanent Shocks Steven R. Genadie Gaduate School of Business, Stanfod Univesity Andey Malenko Gaduate
More informationMATHS FOR ENGINEERS ALGEBRA TUTORIAL 8 MATHEMATICAL PROGRESSIONS AND SERIES
MATHS FOR ENGINEERS ALGEBRA TUTORIAL 8 MATHEMATICAL PROGRESSIONS AND SERIES O completio of this ttoial yo shold be able to do the followig. Eplai aithmetical ad geometic pogessios. Eplai factoial otatio
More informationChromatic number and spectral radius
Linea Algeba and its Applications 426 2007) 810 814 www.elsevie.com/locate/laa Chomatic numbe and spectal adius Vladimi Nikifoov Depatment of Mathematical Sciences, Univesity of Memphis, Memphis, TN 38152,
More informationSkps Media
Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps
More informationTopic 4a Introduction to Root Finding & Bracketing Methods
/8/18 Couse Instucto D. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: cumpf@utep.edu Topic 4a Intoduction to Root Finding & Backeting Methods EE 4386/531 Computational Methods in EE Outline
More informationThe Strain Compatibility Equations in Polar Coordinates RAWB, Last Update 27/12/07
The Stain Compatibility Equations in Pola Coodinates RAWB Last Update 7//7 In D thee is just one compatibility equation. In D polas it is (Equ.) whee denotes the enineein shea (twice the tensoial shea)
More informationarxiv: v1 [math.co] 1 Apr 2011
Weight enumeation of codes fom finite spaces Relinde Juius Octobe 23, 2018 axiv:1104.0172v1 [math.co] 1 Ap 2011 Abstact We study the genealized and extended weight enumeato of the - ay Simplex code and
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More informationScientific Computing II
Scientific Computing II Conjugate Gadient Methods Michael Bade Summe 2014 Conjugate Gadient Methods, Summe 2014 1 Families of Iteative Solves elaxation methods: Jacobi-, Gauss-Seidel-Relaxation,... Ove-Relaxation-Methods
More informationCylindrical and Spherical Coordinate Systems
Clindical and Spheical Coodinate Sstems APPENDIX A In Section 1.2, we leaned that the Catesian coodinate sstem is deined b a set o thee mtall othogonal saces, all o which ae planes. The clindical and spheical
More informationAvailable online through ISSN
Intenational eseach Jounal of Pue Algeba -() 01 98-0 Available online though wwwjpainfo ISSN 8 907 SOE ESULTS ON THE GOUP INVESE OF BLOCK ATIX OVE IGHT OE DOAINS Hanyu Zhang* Goup of athematical Jidong
More informationBrief summary of functional analysis APPM 5440 Fall 2014 Applied Analysis
Bief summay of functional analysis APPM 5440 Fall 014 Applied Analysis Stephen Becke, stephen.becke@coloado.edu Standad theoems. When necessay, I used Royden s and Keyzsig s books as a efeence. Vesion
More informationSurveillance Points in High Dimensional Spaces
Société de Calcul Mathématique SA Tools fo decision help since 995 Suveillance Points in High Dimensional Spaces by Benad Beauzamy Januay 06 Abstact Let us conside any compute softwae, elying upon a lage
More informationA Bijective Approach to the Permutational Power of a Priority Queue
A Bijective Appoach to the Pemutational Powe of a Pioity Queue Ia M. Gessel Kuang-Yeh Wang Depatment of Mathematics Bandeis Univesity Waltham, MA 02254-9110 Abstact A pioity queue tansfoms an input pemutation
More informationLecture 10. Vertical coordinates General vertical coordinate
Lectue 10 Vetical coodinates We have exclusively used height as the vetical coodinate but thee ae altenative vetical coodinates in use in ocean models, most notably the teainfollowing coodinate models
More informationResearch Design - - Topic 17 Multiple Regression & Multiple Correlation: Two Predictors 2009 R.C. Gardner, Ph.D.
Reseach Design - - Topic 7 Multiple Regession & Multiple Coelation: Two Pedictos 009 R.C. Gadne, Ph.D. Geneal Rationale and Basic Aithmetic fo two pedictos Patial and semipatial coelation Regession coefficients
More informationENGI 4430 Non-Cartesian Coordinates Page xi Fy j Fzk from Cartesian coordinates z to another orthonormal coordinate system u, v, ˆ i ˆ ˆi
ENGI 44 Non-Catesian Coodinates Page 7-7. Conesions between Coodinate Systems In geneal, the conesion of a ecto F F xi Fy j Fzk fom Catesian coodinates x, y, z to anothe othonomal coodinate system u,,
More informationarxiv:gr-qc/ v1 17 Feb 2001
Scatteing poblem of scala wave in womhole geomety Sng-Won Kim Depatment of Science Edcation, Ewha Women s Univesity, Seol -75, Koea Jne 5, 7 In this pape, we stdy the scatteing poblem of the scala wave
More informationTHE JEU DE TAQUIN ON THE SHIFTED RIM HOOK TABLEAUX. Jaejin Lee
Koean J. Math. 23 (2015), No. 3, pp. 427 438 http://dx.doi.og/10.11568/kjm.2015.23.3.427 THE JEU DE TAQUIN ON THE SHIFTED RIM HOOK TABLEAUX Jaejin Lee Abstact. The Schensted algoithm fist descibed by Robinson
More informationA pathway to matrix-variate gamma and normal densities
Linea Algeba and its Applications 396 005 317 38 www.elsevie.com/locate/laa A pathway to matix-vaiate gamma and nomal densities A.M. Mathai Depatment of Mathematics and Statistics, McGill Univesity, 805
More information0606 ADDITIONAL MATHEMATICS
UNIVERSITY OF CAMBRIDGE INTERNATIONAL EXAMINATIONS Intenational Geneal Cetificate of Seconday Education MARK SCHEME fo the Octobe/Novembe 011 question pape fo the guidance of teaches 0606 ADDITIONAL MATHEMATICS
More informationPHYS 1444 Lecture #5
Shot eview Chapte 24 PHYS 1444 Lectue #5 Tuesday June 19, 212 D. Andew Bandt Capacitos and Capacitance 1 Coulom s Law The Fomula QQ Q Q F 1 2 1 2 Fomula 2 2 F k A vecto quantity. Newtons Diection of electic
More informationMath 2263 Solutions for Spring 2003 Final Exam
Math 6 Solutions fo Sping Final Exam ) A staightfowad appoach to finding the tangent plane to a suface at a point ( x, y, z ) would be to expess the cuve as an explicit function z = f ( x, y ), calculate
More informationTemporal-Difference Learning
.997 Decision-Making in Lage-Scale Systems Mach 17 MIT, Sping 004 Handout #17 Lectue Note 13 1 Tempoal-Diffeence Leaning We now conside the poblem of computing an appopiate paamete, so that, given an appoximation
More information10/04/18. P [P(x)] 1 negl(n).
Mastemath, Sping 208 Into to Lattice lgs & Cypto Lectue 0 0/04/8 Lectues: D. Dadush, L. Ducas Scibe: K. de Boe Intoduction In this lectue, we will teat two main pats. Duing the fist pat we continue the
More informationA Relativistic Electron in a Coulomb Potential
A Relativistic Electon in a Coulomb Potential Alfed Whitehead Physics 518, Fall 009 The Poblem Solve the Diac Equation fo an electon in a Coulomb potential. Identify the conseved quantum numbes. Specify
More informationEEO 401 Digital Signal Processing Prof. Mark Fowler
EEO 41 Digital Signal Pocessing Pof. Mak Fowle Note Set #31 Linea Phase FIR Design Optimum Equiipple (Paks-McClellan) Reading: Sect. 1.2.4 1.2.6 of Poakis & Manolakis 1/2 Motivation The window method and
More informationGravity and isostasy
Gavity and isostasy Reading: owle p60 74 Theoy of gavity Use two of Newton s laws: ) Univesal law of gavitation: Gmm = m m Univesal gavitational constant G=6.67 x 0 - Nm /kg ) Second law of motion: = ma
More informationBerkeley Math Circle AIME Preparation March 5, 2013
Algeba Toolkit Rules of Thumb. Make sue that you can pove all fomulas you use. This is even bette than memoizing the fomulas. Although it is best to memoize, as well. Stive fo elegant, economical methods.
More informationComputing Approximate GCD of Multivariate Polynomials by Structure Total Least Norm 1)
MM Reseach Pepints, 388 41 MMRC, AMSS, Academia Sinica No 23, Decembe 24 Computing Appoximate GCD of Multivaiate Polynomials by Stuctue Total Least Nom 1) Lihong Zhi and Zhengfeng Yang Key Laboatoy of
More information(read nabla or del) is defined by, k. (9.7.1*)
9.7 Gadient of a scala field. Diectional deivative Some of the vecto fields in applications can be obtained fom scala fields. This is vey advantageous because scala fields can be handled moe easily. The
More informationMany Electron Atoms. Electrons can be put into approximate orbitals and the properties of the many electron systems can be catalogued
Many Electon Atoms The many body poblem cannot be solved analytically. We content ouselves with developing appoximate methods that can yield quite accuate esults (but usually equie a compute). The electons
More informationQuestion Bank. Section A. is skew-hermitian matrix. is diagonalizable. (, ) , Evaluate (, ) 12 about = 1 and = Find, if
Subject: Mathematics-I Question Bank Section A T T. Find the value of fo which the matix A = T T has ank one. T T i. Is the matix A = i is skew-hemitian matix. i. alculate the invese of the matix = 5 7
More informationRigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018
Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining
More informationOn the integration of the equations of hydrodynamics
Uebe die Integation de hydodynamischen Gleichungen J f eine u angew Math 56 (859) -0 On the integation of the equations of hydodynamics (By A Clebsch at Calsuhe) Tanslated by D H Delphenich In a pevious
More informationQuasi-Randomness and the Distribution of Copies of a Fixed Graph
Quasi-Randomness and the Distibution of Copies of a Fixed Gaph Asaf Shapia Abstact We show that if a gaph G has the popety that all subsets of vetices of size n/4 contain the coect numbe of tiangles one
More informationF-IF Logistic Growth Model, Abstract Version
F-IF Logistic Gowth Model, Abstact Vesion Alignments to Content Standads: F-IFB4 Task An impotant example of a model often used in biology o ecology to model population gowth is called the logistic gowth
More informationThe Substring Search Problem
The Substing Seach Poblem One algoithm which is used in a vaiety of applications is the family of substing seach algoithms. These algoithms allow a use to detemine if, given two chaacte stings, one is
More information0606 ADDITIONAL MATHEMATICS 0606/01 Paper 1, maximum raw mark 80
UNIVERSITY OF CAMBRIDGE INTERNATIONAL EXAMINATIONS Intenational Geneal Cetificate of Seconday Education MARK SCHEME fo the Octobe/Novembe 009 question pape fo the guidance of teaches 0606 ADDITIONAL MATHEMATICS
More informationC/CS/Phys C191 Shor s order (period) finding algorithm and factoring 11/12/14 Fall 2014 Lecture 22
C/CS/Phys C9 Sho s ode (peiod) finding algoithm and factoing /2/4 Fall 204 Lectue 22 With a fast algoithm fo the uantum Fouie Tansfom in hand, it is clea that many useful applications should be possible.
More informationAnalytic Evaluation of two-electron Atomic Integrals involving Extended Hylleraas-CI functions with STO basis
Analytic Evaluation of two-electon Atomic Integals involving Extended Hylleaas-CI functions with STO basis B PADHY (Retd.) Faculty Membe Depatment of Physics, Khalikote (Autonomous) College, Behampu-760001,
More informationMath 124B February 02, 2012
Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial
More informationPart V: Closed-form solutions to Loop Closure Equations
Pat V: Closed-fom solutions to Loop Closue Equations This section will eview the closed-fom solutions techniques fo loop closue equations. The following thee cases will be consideed. ) Two unknown angles
More informationQuestion 1: The dipole
Septembe, 08 Conell Univesity, Depatment of Physics PHYS 337, Advance E&M, HW #, due: 9/5/08, :5 AM Question : The dipole Conside a system as discussed in class and shown in Fig.. in Heald & Maion.. Wite
More informationA Low-Complexity and High-Throughput RTL Design of a BCH (15,7) Decoder
ITB J. ICT, Vol., No.,, - A Low-Compleity and High-Thoghpt RTL Design of a BCH (,) Decode Henda Setiawan Electical Engineeing Depatment, Islamic Univesity of Indonesia Jl. Kaliang Km.. Yogyakata, Indonesia,
More informationA Converse to Low-Rank Matrix Completion
A Convese to Low-Rank Matix Completion Daniel L. Pimentel-Alacón, Robet D. Nowak Univesity of Wisconsin-Madison Abstact In many pactical applications, one is given a subset Ω of the enties in a d N data
More informationCSCE 478/878 Lecture 4: Experimental Design and Analysis. Stephen Scott. 3 Building a tree on the training set Introduction. Outline.
In Homewok, you ae (supposedly) Choosing a data set 2 Extacting a test set of size > 3 3 Building a tee on the taining set 4 Testing on the test set 5 Repoting the accuacy (Adapted fom Ethem Alpaydin and
More information1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More informationMath 301: The Erdős-Stone-Simonovitz Theorem and Extremal Numbers for Bipartite Graphs
Math 30: The Edős-Stone-Simonovitz Theoem and Extemal Numbes fo Bipatite Gaphs May Radcliffe The Edős-Stone-Simonovitz Theoem Recall, in class we poved Tuán s Gaph Theoem, namely Theoem Tuán s Theoem Let
More informationComputational Methods of Solid Mechanics. Project report
Computational Methods of Solid Mechanics Poject epot Due on Dec. 6, 25 Pof. Allan F. Bowe Weilin Deng Simulation of adhesive contact with molecula potential Poject desciption In the poject, we will investigate
More informationPushdown Automata (PDAs)
CHAPTER 2 Context-Fee Languages Contents Context-Fee Gammas definitions, examples, designing, ambiguity, Chomsky nomal fom Pushdown Automata definitions, examples, euivalence with context-fee gammas Non-Context-Fee
More informationgr0 GRAPHS Hanan Samet
g0 GRPHS Hanan Samet ompute Science epatment and ente fo utomation Reseach and Institute fo dvanced ompute Studies Univesity of Mayland ollege Pak, Mayland 074 e-mail: hjs@umiacs.umd.edu opyight 1997 Hanan
More information