DOCSLIB.ORG

Arxiv:2012.03523V3 [Math.NT] 16 May 2021 Favne Cetfi Optn Eerh(SR Spr Ft of Part (CM4)

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Arxiv:2012.03523V3 [Math.NT] 16 May 2021 Favne Cetfi Optn Eerh(SR Spr Ft of Part (CM4) WRONSKIAN´ ALGEBRA AND BROADHURST–ROBERTS QUADRATIC RELATIONS YAJUN ZHOU To the memory of Dr. W (1986–2020) ABSTRACT. Throughalgebraicmanipulationson Wronskian´ matrices whose entries are reducible to Bessel moments, we present a new analytic proof of the quadratic relations conjectured by Broadhurst and Roberts, along with some generalizations. In the Wronskian´ framework, we reinterpret the de Rham intersection pairing through polynomial coefficients in Vanhove’s differential operators, and compute the Betti inter- section pairing via linear sum rules for on-shell and off-shell Feynman diagrams at threshold momenta. From the ideal generated by Broadhurst–Roberts quadratic relations, we derive new non-linear sum rules for on-shell Feynman diagrams, including an infinite family of determinant identities that are compatible with Deligne’s conjectures for critical values of motivic L-functions. CONTENTS 1. Introduction 1 2. W algebra 5 2.1. Wronskian´ matrices and Vanhove operators 5 2.2. Block diagonalization of on-shell Wronskian´ matrices 8 2.3. Sumrulesforon-shellandoff-shellBesselmoments 11 3. W⋆ algebra 20 3.1. Cofactors of Wronskians´ 21 3.2. Threshold behavior of Wronskian´ cofactors 24 3.3. Quadratic relations for off-shell Bessel moments 30 4. Broadhurst–Roberts quadratic relations 37 4.1. Quadratic relations for on-shell Bessel moments 37 4.2. Arithmetic applications to Feynman diagrams 43 4.3. Alternative representations of de Rham matrices 49 References 53 arXiv:2012.03523v3 [math.NT] 16 May 2021 Date: May 18, 2021. Keywords: Bessel moments, Feynman integrals, Wronskian´ matrices, Bernoulli numbers Subject Classification (AMS 2020): 11B68, 33C10, 34M35 (Primary) 81T18, 81T40 (Secondary) * This research was supported in part by the Applied Mathematics Program within the Department of Energy (DOE) Office of Advanced Scientific Computing Research (ASCR) as part of the Collaboratory on Mathematics for Mesoscopic Modeling of Materials (CM4). i BROADHURST–ROBERTSQUADRATICRELATIONS 1 1. INTRODUCTION In perturbative quantum field theory, one frequently encounters Feynman diagrams in the form of on-shell Bessel moments [33, 29, 3, 32, 34, 9, 12, 10, 13, 14, 25] ∞ a b n IKM(a,b;n) := [I0(t)] [K0(t)] t dt (1.1) Z0 = 1 π t cosθ = t coshu for certain non-negative integers a,b,n Z 0, where I0(t) π 0 e dθ and K0(t) 0∞ e− du are modified Bessel functions of zeroth order.∈ These≥ Bessel moments satisfy a wealth of algebraic relations, R R which had been discovered by Broadhurst and collaborators through numerical experimentations [9, 12, 10, 13, 14], before becoming formally proven truths [50, 48, 49, 53, 54, 25]. The Broadhurst–Mellit determinant formulae [9, Conjectures 4 and 7] were among the earliest verified [49] non-linear algebraic relations for Bessel moments: the k k matrices Mk := (IKM(a,2k +1 a;2b × − − 1))1 a,b k and Nk := (IKM(a,2k + 2 a;2b 1))1 a,b k have determinants ≤ ≤ − − ≤ ≤ k k j j (k+1)2/2 k+1 k+1 j (2 j) − π 2π (2 j 1) − detMk = , detNk = − , (1.2) ∏ 2 j+1 Γ + ∏ j j=1 (2 j + 1) ((k 1)/2) j=1 (2 j) p Γ = x 1 t where (x) : 0∞ t − e− dt for x > 0. = T These determinant formulae are refined by the Broadhurst–Roberts quadratic relations Bm PmDmPm R + b 1 a m 3 [10, 13, 14] for the period matrix Pm := (( 1) − π − 2 IKM(a,m + 2 a;2b 1)) m+1 , in which 1 a,b 2 m+1 m+1 − m−+1 m+−1 ≤ ≤ the Betti matrix Bm Q 2 × 2 and the de Rham matrix Dm Q 2 × 2 are filled with rational numbers. Fresán, Sabbah∈ and Yu have recently proved an equiv∈alent form [25, Theorem 1.4] of the Broadhurst–Roberts quadratic relations, with an algebro-geometric tour de force that draws on their profound insights into certain exponential motives [22], whose associated Galois representations are intimately related to moments of Kloosterman sums [23, 24] (analogs of on-shell Bessel moments in positive characterisitic). In this paper, we present an analytic proof of the Broadhurst–Roberts quadratic relations, by revis- iting the Wronskian´ method in our verification of the Broadhurst–Mellit determinant formulae [49]. Our Wronskian-based´ proof also produces a new family of quadratic relations, which embody the Broadhurst–Roberts relations for on-shell Bessel moments as special cases. n n n 0 In §2, writing D f (u) := d f (u)/du for n Z>0, and D f (u) := f (u), we will consider a Wronskian´ i 1 ∈ matrix W[ f1(u),..., fm(u)] = (D − f j(u))1 i, j m whose entries are reducible to on-shell Bessel moments ≤ ≤ (more precisely, entries of Pm and Pm 2) as u 1−, and whose determinant detW := W[ f1(u),..., fm(u)] = i 1 − → det(D − f j(u))1 i, j m is explicitly computable through our proof of the Broadhurst–Mellit determinant ≤ ≤ ⋆ ⋆ formulae [49]. In §3, our key observation is the adjugate relation W W = WW = WIm, where the ⋆ ⋆ T m m matrix adjugate W := adjW of W[ f1(u),..., fm(u)] factorizes into W = SW U, with S Q × and m m T ∈ U/W Q(u) × . In §4, we explore the u 1− limit of the matrix equation SW UW = WIm, and establish∈ quadratic relations among on-shell Bessel→ moments. For a generic parameter u (0,1), the individual elements of the Wronskian´ matrix W[ f1(u),..., fm(u)] in §2 will be representable by∈ off-shell Bessel moments [49, Definition 2.1] ∞ a b n IKM(a + 1,b;n u) := I0(√ut)[I0(t)] [K0(t)] t dt, (1.3) | Z0 ∞ a b n eIKM(a,b + 1;n u) := K0(√ut)[I0(t)] [K0(t)] t dt, (1.4) | 0 Z e 2 YAJUN ZHOU and differentiated Bessel moments [49, Definition 2.3] ∞ a b n+1 IKM´ (a + 1,b;n u) := + I1(√ut)[I0(t)] [K0(t)] t dt, (1.5) | Z0 ∞ a b n+1 IKM´ (a,b + 1;n u) := K1(√ut)[I0(t)] [K0(t)] t dt, (1.6) | − Z0 where I1(x) = +d I0(x)/d x and K1(x) = d K0(x)/d x. Accordingly, our efforts in §3 will culminate in a quadratic relation involving these generalized− Bessel moments and some explicitly computable rational functions (in lieu of the rational numbers in the Betti matrix B and the de Rham matrix D of the Broadhurst–Roberts formulation), as stated below. 5 3( 1)m Theorem 1.1 (Wronskian´ algebra). For m Z>0,u 0, − 2− , normalize off-shell Feynman diagrams as ∈ ∈ IKM(1,m+1;1 u)+(m+1)IKM(1,m+1;1 u) = (|m+2)π(m+1)/2 | , j 1, (u) := eIKM( j,m+2 j;1 u) e m + (1.7) m, j (m+1)/2+−1 j | , j Z 2, 2 1 , F π − ∈ ∩ IKM( j m ,m+2+ m j;1 u) e − 2 2 − | m + (m⌊+1)⌋/2+ m/2⌊+1 ⌋j , j Z 2 2,m , π ⌊ ⌋ − ∈ ∩ e and define a polynomial + m 1 2 m(u) := u 2 ∏ (u n ), (1.8) L n Z [1,m+1] − n ∈m+∩1 (mod 2) ≡ 5 3( 1)m where x is the greatest integer less than or equal to x. For u 0, − − , the Wro´nskian matrix ⌊ ⌋ ∈ 2 Wm(u) := W[ m,1(u),..., m,m(u)] has determinant F F m + 1 ( 1) m/4 [(m + 1)!]m 1 := = ⌊ ⌋ Wm(u) detWm(u) −m(m 1)/2 m+1 m/2 (1.9) m + 2 2 ∏ nn m(u) − n=1 |L | and satisfies a quadratic relation 1 T = [Vm(u)]− Wm(u)SmWm(u) (1.10) m(u) |L | for matrices V (u) = ( 1)m+1VT (u) Q(u)m m and S = ( 1)m+1ST Qm m. m − m ∈ × m − m ∈ × Explicitly, the Vanhove matrix V (u) Q(u)m m is upper-left triangular, with entries1 m ∈ × m n a b+1 = a+n+m+1 n a D − − ℓm,n(u) (Vm(u))a,b : ∑ ( 1) − , (1.11) n=a+b 1 − b 1 m(u) − − L 2 where the polynomials ℓm, j(u) Z[u] are defined through the Vanhove–Verrill relation (cf. [43, §3.6.1] and [44, §3.6.1]): ∈ m m +1 k 2 + j = m + 1 k m 1 α1 + αn αn+1 ∑ ℓm, j(u)D : uϑ ∑ u − ∑ (ϑ k) − ∏ αn(αn m 2)(ϑ k n) − , − − − − j=0 k=1 n Z [1,k],αk+1=1 n=1 ∈αn∩Z [1,m+1] b α∈n+1∩ αn 2 b b ≤ − (1.12) 1 n n! Throughout this article, we write := for the standard binomial coefficient with n Z 0,m Z [0,n]. By m m!(n m)! ∈ ≥ ∈ ∩ n = − B convention, we also define 0 when n Z 0,k Zr [0,n]. Empty sums like ∑n=A( ) with A > B are treated as zero. k ∈ ≥ ∈ ··· 2These polynomials can be constructed by a finite enumeration over integer parameters α ,...,α in (1.12). See Table I 1 n for a partial list of the Vanhove–Verrill polynomials. BROADHURST–ROBERTSQUADRATICRELATIONS 3 = m j TABLE I. A partial list of Vanhove’s operators Lm ∑ j=0 ℓm, j(u)D m Lm e 1 (u 4)uD1 + (u 2)D0 − − 2 (eu 9)(u 1)uD2 + (3u2 20u + 9)D1 + (u 3)D0 − − − − 3 (u 16)(u 4)u2 D3 + 6u(u2 15u + 32)D2 + (7u2 68u + 64)D1 + (u 4)D0 − − − − − 4 (u 25)(u 9)(u 1)u2 D4 + 2u(5u3 140u2 + 777u 450)D3 + (25u3 518u2 + 1839u 450)D2 + (3u 5)(5u 57)D1 + (u 5)D0 − − − − − − − − − − 5 (u 36)(u 16)(u 4)u3 D5 + 5u2(3u3 140u2 + 1568u 3456)D4 + u(65u3 2408u2 + 19836u 27648)D3 − − − − − − − +6(15u3 406u2 + 2078u 1152)D2 + (31u2 516u + 1020)D1 + (u 6)D0 − − − − 6 (u 49)(u 25)(u 9)(u 1)u3 D6 + 3u2(7u4 504u3 + 9870u2 51664u + 33075)D5 − − − − − − +2u(70u4 4179u3 + 64749u2 247993u + 99225)D4 + 2(175u4 8232u3 + 92394u2 217012u + 33075)D3 − − − − +(301u3 10218u2 + 69561u 62020)D2 + 3(21u2 428u + 1071)D1 + (u 7)D0 − − − − 7 (u 64)(u 36)(u 16)(u 4)u4 D7 + 28u3(u4 105u3 + 3276u2 32800u + 73728)D6 − − − − − − +2u2(133u4 11928u3 + 307950u2 2432512u + 3981312)D5 + 30u(35u4 2590u3 + 52420u2 297856u + 294912)D4 − − − − +3(567u4 32820u3 + 474944u2 1616896u + 589824)D3 + 6(161u3 6645u2 + 57096u 70912)D2 − − − − +(127u2 3076u + 9312)D1 + (u 8)D0 − − 8 (u 81)(u 49)(u 25)(u 9)(u 1)u4 D8 + 12u3(3u5 440u4 + 20482u3 345620u2 + 1762035u 1190700)D7 − − − − − − − − +6u2(77u5 9856u4 + 391358u3 5455164u2 + 21884733u 10716300)D6 − − − +6u(441u5 48048u4 + 1569150u3 17075168u2 + 48951657u 14288400)D5 − − − +3(2317u5 207295u4 + 5264259u3 40651797u2 + 68940108u 7144200)D4 − − − +6(1295u4 89980u3 + 1614657u2 7325514u + 4441878)D3 + (3025u3 148126u2 + 1547883u 2474334)D2 − − − − +(255u2 7172u + 25461)D1 + (u 9)D0 − − A B = 1 m m = Sm Sm 3 for (ϑ f )(u) D [uf (u)]; the matrix Sm Q × can be partitioned into Sm C D for ∈ Sm Sm m +2 a 2 − + m + b + + + + ⌊ ⌋ s m 2 a 2 1 s a+b+m+1 [1 (m 1)δa,1][1 (m 1)δb,1] ∑ ( 1) m+1 −+s ⌊ ⌋ − 1 + ( 1) = − 2 b 1 (SA ) := − s 1 − (1.13) m a,b 1 m a b 1+( 1)m m 2 − + − 1 ( 1) 2 2 − 4 −⌊ 2 ⌋− (a 1)!(m + 2 a)! − j k − − where a,b Z 1, m + 1 , ∈ ∩ 2 B m+1 C (S )a,b = (1) (S )b ,a m ′ − m ′ m +2 a 2 − + m + s m+1 +s + + ⌊ ⌋ s m 2 a 2 1 + b′ 2 a+b +m [1 (m 1)δa,1] ∑ ( 1) m+1 −+s ⌊ ⌋+ − ( 1) + 1 + ( 1) ′ = − 2 b′ 1 − b′ 1 := − s 1 (1.14) 1 m a 1+( 1)m b 1+(1)m h m+1 i 2 (a 1)m+ − + ′ + − − ( 1) − 2 − 4 2 4 − 2 (a 1)!(m + 2 a)! − j k j k − − m m 1 where a Z 1, + 1 ,b Z 1, − , and ∈ ∩ 2 ′ ∈ ∩ 2 + m+1 + + + m 1 2 1 + a′ + b′ m 1 m 1 D 1 ( 1) ( 4) − [1 ( 1) ][1 ( 1) ] 2 2 (S )a ,b := − − − − (1.15) m ′ ′ + 2 a′ + b′ + + 2 m 1 ! ( 1) 2 2 a′ 1 b′ 1 2 − m 1 where a′,b′ Z 1, − .
Load more

Recommended publications
  • Linear Algebra: Example Sheet 2 of 4
    Michaelmas Term 2014 SJW Linear Algebra: Example Sheet 2 of 4 1. (Another proof of the row rank column rank equality.) Let A be an m × n matrix of (column) rank r. Show that r is the least integer for which A factorises as A = BC with B 2 Matm;r(F) and C 2 Matr;n(F). Using the fact that (BC)T = CT BT , deduce that the (column) rank of AT equals r. 2. Write down the three types of elementary matrices and find their inverses. Show that an n × n matrix A is invertible if and only if it can be written as a product of elementary matrices. Use this method to find the inverse of 0 1 −1 0 1 @ 0 0 1 A : 0 3 −1 3. Let A and B be n × n matrices over a field F . Show that the 2n × 2n matrix IB IB C = can be transformed into D = −A 0 0 AB by elementary row operations (which you should specify). By considering the determinants of C and D, obtain another proof that det AB = det A det B. 4. (i) Let V be a non-trivial real vector space of finite dimension. Show that there are no endomorphisms α; β of V with αβ − βα = idV . (ii) Let V be the space of infinitely differentiable functions R ! R. Find endomorphisms α; β of V which do satisfy αβ − βα = idV . 5. Find the eigenvalues and give bases for the eigenspaces of the following complex matrices: 0 1 1 0 1 0 1 1 −1 1 0 1 1 −1 1 @ 0 3 −2 A ; @ 0 3 −2 A ; @ −1 3 −1 A : 0 1 0 0 1 0 −1 1 1 The second and third matrices commute; find a basis with respect to which they are both diagonal.
    [Show full text]
  • Linear Algebra Review James Chuang
    Linear algebra review James Chuang December 15, 2016 Contents 2.1 vector-vector products ............................................... 1 2.2 matrix-vector products ............................................... 2 2.3 matrix-matrix products ............................................... 4 3.2 the transpose .................................................... 5 3.3 symmetric matrices ................................................. 5 3.4 the trace ....................................................... 6 3.5 norms ........................................................ 6 3.6 linear independence and rank ............................................ 7 3.7 the inverse ...................................................... 7 3.8 orthogonal matrices ................................................. 8 3.9 range and nullspace of a matrix ........................................... 8 3.10 the determinant ................................................... 9 3.11 quadratic forms and positive semidefinite matrices ................................ 10 3.12 eigenvalues and eigenvectors ........................................... 11 3.13 eigenvalues and eigenvectors of symmetric matrices ............................... 12 4.1 the gradient ..................................................... 13 4.2 the Hessian ..................................................... 14 4.3 gradients and hessians of linear and quadratic functions ............................. 15 4.5 gradients of the determinant ............................................ 16 4.6 eigenvalues
    [Show full text]
  • A New Algorithm to Obtain the Adjugate Matrix Using CUBLAS on GPU González, H.E., Carmona, L., J.J
    ISSN: 2277-3754 ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 7, Issue 4, October 2017 A New Algorithm to Obtain the Adjugate Matrix using CUBLAS on GPU González, H.E., Carmona, L., J.J. Information Technology Department, ININ-ITTLA Information Technology Department, ININ proved to be the best option for most practical applications. Abstract— In this paper a parallel code for obtain the Adjugate The new transformation proposed here can be obtained from Matrix with real coefficients are used. We postulate a new linear it. Next, we briefly review this topic. transformation in matrix product form and we apply this linear transformation to an augmented matrix (A|I) by means of both a minimum and a complete pivoting strategies, then we obtain the II. LU-MATRICIAL DECOMPOSITION WITH GT. Adjugate matrix. Furthermore, if we apply this new linear NOTATION AND DEFINITIONS transformation and the above pivot strategy to a augmented The problem of solving a linear system of equations matrix (A|b), we obtain a Cramer’s solution of the linear system of Ax b is central to the field of matrix computation. There are equations. That new algorithm present an O n3 computational several ways to perform the elimination process necessary for complexity when n . We use subroutines of CUBLAS 2nd A, b R its matrix triangulation. We will focus on the Doolittle-Gauss rd and 3 levels in double precision and we obtain correct numeric elimination method: the algorithm of choice when A is square, results. dense, and un-structured. nxn Index Terms—Adjoint matrix, Adjugate matrix, Cramer’s Let us assume that AR is nonsingular and that we rule, CUBLAS, GPU.
    [Show full text]
  • Field of U-Invariants of Adjoint Representation of the Group GL (N, K)
    Field of U-invariants of adjoint representation of the group GL(n,K) K.A.Vyatkina A.N.Panov ∗ It is known that the field of invariants of an arbitrary unitriangular group is rational (see [1]). In this paper for adjoint representation of the group GL(n,K) we present the algebraically system of generators of the field of U-invariants. Note that the algorithm for calculating generators of the field of invariants of an arbitrary rational representation of an unipotent group was presented in [2, Chapter 1]. The invariants was constructed using induction on the length of Jordan-H¨older series. In our case the length is equal to n2; that is why it is difficult to apply the algorithm. −1 Let us consider the adjoint representation AdgA = gAg of the group GL(n,K), where K is a field of zero characteristic, on the algebra of matri- ces M = Mat(n,K). The adjoint representation determines the representation −1 ρg on K[M] (resp. K(M)) by formula ρgf(A)= f(g Ag). Let U be the subgroup of upper triangular matrices in GL(n,K) with units on the diagonal. A polynomial (rational function) f on M is called an U- invariant, if ρuf = f for every u ∈ U. The set of U-invariant rational functions K(M)U is a subfield of K(M). Let {xi,j} be a system of standard coordinate functions on M. Construct a X X∗ ∗ X X∗ X∗ X matrix = (xij). Let = (xij) be its adjugate matrix, · = · = det X · E. Denote by Jk the left lower corner minor of order k of the matrix X.
    [Show full text]
  • Formalizing the Matrix Inversion Based on the Adjugate Matrix in HOL4 Liming Li, Zhiping Shi, Yong Guan, Jie Zhang, Hongxing Wei
    Formalizing the Matrix Inversion Based on the Adjugate Matrix in HOL4 Liming Li, Zhiping Shi, Yong Guan, Jie Zhang, Hongxing Wei To cite this version: Liming Li, Zhiping Shi, Yong Guan, Jie Zhang, Hongxing Wei. Formalizing the Matrix Inversion Based on the Adjugate Matrix in HOL4. 8th International Conference on Intelligent Information Processing (IIP), Oct 2014, Hangzhou, China. pp.178-186, 10.1007/978-3-662-44980-6_20. hal-01383331 HAL Id: hal-01383331 https://hal.inria.fr/hal-01383331 Submitted on 18 Oct 2016 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Distributed under a Creative Commons Attribution| 4.0 International License Formalizing the Matrix Inversion Based on the Adjugate Matrix in HOL4 Liming LI1, Zhiping SHI1?, Yong GUAN1, Jie ZHANG2, and Hongxing WEI3 1 Beijing Key Laboratory of Electronic System Reliability Technology, Capital Normal University, Beijing 100048, China [email protected], [email protected] 2 College of Information Science & Technology, Beijing University of Chemical Technology, Beijing 100029, China 3 School of Mechanical Engineering and Automation, Beihang University, Beijing 100083, China Abstract. This paper presents the formalization of the matrix inversion based on the adjugate matrix in the HOL4 system.
    [Show full text]
  • Linearized Polynomials Over Finite Fields Revisited
    Linearized polynomials over finite fields revisited∗ Baofeng Wu,† Zhuojun Liu‡ Abstract We give new characterizations of the algebra Ln(Fqn ) formed by all linearized polynomials over the finite field Fqn after briefly sur- veying some known ones. One isomorphism we construct is between ∨ L F n F F n n( q ) and the composition algebra qn ⊗Fq q . The other iso- morphism we construct is between Ln(Fqn ) and the so-called Dickson matrix algebra Dn(Fqn ). We also further study the relations between a linearized polynomial and its associated Dickson matrix, generalizing a well-known criterion of Dickson on linearized permutation polyno- mials. Adjugate polynomial of a linearized polynomial is then intro- duced, and connections between them are discussed. Both of the new characterizations can bring us more simple approaches to establish a special form of representations of linearized polynomials proposed re- cently by several authors. Structure of the subalgebra Ln(Fqm ) which are formed by all linearized polynomials over a subfield Fqm of Fqn where m|n are also described. Keywords Linearized polynomial; Composition algebra; Dickson arXiv:1211.5475v2 [math.RA] 1 Jan 2013 matrix algebra; Representation. ∗Partially supported by National Basic Research Program of China (2011CB302400). †Key Laboratory of Mathematics Mechanization, AMSS, Chinese Academy of Sciences, Beijing 100190, China. Email: [email protected] ‡Key Laboratory of Mathematics Mechanization, AMSS, Chinese Academy of Sciences, Beijing 100190, China. Email: [email protected] 1 1 Introduction n Let Fq and Fqn be the finite fields with q and q elements respectively, where q is a prime or a prime power.
    [Show full text]
  • Lecture 5: Matrix Operations: Inverse
    Lecture 5: Matrix Operations: Inverse • Inverse of a matrix • Computation of inverse using co-factor matrix • Properties of the inverse of a matrix • Inverse of special matrices • Unit Matrix • Diagonal Matrix • Orthogonal Matrix • Lower/Upper Triangular Matrices 1 Matrix Inverse • Inverse of a matrix can only be defined for square matrices. • Inverse of a square matrix exists only if the determinant of that matrix is non-zero. • Inverse matrix of 퐴 is noted as 퐴−1. • 퐴퐴−1 = 퐴−1퐴 = 퐼 • Example: 2 −1 0 1 • 퐴 = , 퐴−1 = , 1 0 −1 2 2 −1 0 1 0 1 2 −1 1 0 • 퐴퐴−1 = = 퐴−1퐴 = = 1 0 −1 2 −1 2 1 0 0 1 2 Inverse of a 3 x 3 matrix (using cofactor matrix) • Calculating the inverse of a 3 × 3 matrix is: • Compute the matrix of minors for A. • Compute the cofactor matrix by alternating + and – signs. • Compute the adjugate matrix by taking a transpose of cofactor matrix. • Divide all elements in the adjugate matrix by determinant of matrix 퐴. 1 퐴−1 = 푎푑푗(퐴) det(퐴) 3 Inverse of a 3 x 3 matrix (using cofactor matrix) 3 0 2 퐴 = 2 0 −2 0 1 1 0 −2 2 −2 2 0 1 1 0 1 0 1 2 2 2 0 2 3 2 3 0 Matrix of Minors = = −2 3 3 1 1 0 1 0 1 0 2 3 2 3 0 0 −10 0 0 −2 2 −2 2 0 2 2 2 1 −1 1 2 −2 2 Cofactor of A (퐂) = −2 3 3 .∗ −1 1 −1 = 2 3 −3 0 −10 0 1 −1 1 0 10 0 2 2 0 adj A = CT = −2 3 10 2 −3 0 2 2 0 0.2 0.2 0 1 1 A-1 = ∗ adj A = −2 3 10 = −0.2 0.3 1 |퐴| 10 2 −3 0 0.2 −0.3 0 4 Properties of Inverse of a Matrix • (A-1)-1 = A • (AB)-1 = B-1A-1 • (kA)-1 = k-1A-1 where k is a non-zero scalar.
    [Show full text]
  • Contents 2 Matrices and Systems of Linear Equations
    Business Math II (part 2) : Matrices and Systems of Linear Equations (by Evan Dummit, 2019, v. 1.00) Contents 2 Matrices and Systems of Linear Equations 1 2.1 Systems of Linear Equations . 1 2.1.1 Elimination, Matrix Formulation . 2 2.1.2 Row-Echelon Form and Reduced Row-Echelon Form . 4 2.1.3 Gaussian Elimination . 6 2.2 Matrix Operations: Addition and Multiplication . 8 2.3 Determinants and Inverses . 11 2.3.1 The Inverse of a Matrix . 11 2.3.2 The Determinant of a Matrix . 13 2.3.3 Cofactor Expansions and the Adjugate . 16 2.3.4 Matrices and Systems of Linear Equations, Revisited . 18 2 Matrices and Systems of Linear Equations In this chapter, we will discuss the problem of solving systems of linear equations, reformulate the problem using matrices, and then give the general procedure for solving such systems. We will then study basic matrix operations (addition and multiplication) and discuss the inverse of a matrix and how it relates to another quantity known as the determinant. We will then revisit systems of linear equations after reformulating them in the language of matrices. 2.1 Systems of Linear Equations • Our motivating problem is to study the solutions to a system of linear equations, such as the system x + 3x = 5 1 2 : 3x1 + x2 = −1 ◦ Recall that a linear equation is an equation of the form a1x1 + a2x2 + ··· + anxn = b, for some constants a1; : : : an; b, and variables x1; : : : ; xn. ◦ When we seek to nd the solutions to a system of linear equations, this means nding all possible values of the variables such that all equations are satised simultaneously.
    [Show full text]
  • Lecture No. 6
    Preliminaries: linear algebra State-space models System interconnections Control Theory (035188) lecture no. 6 Leonid Mirkin Faculty of Mechanical Engineering Technion | IIT Preliminaries: linear algebra State-space models System interconnections Outline Preliminaries: linear algebra State-space models System interconnections in terms of state-space realizations Preliminaries: linear algebra State-space models System interconnections Outline Preliminaries: linear algebra State-space models System interconnections in terms of state-space realizations Preliminaries: linear algebra State-space models System interconnections Linear algebra, notation and terminology A matrix A Rn×m is an n m table of real numbers, 2 × 2 3 a a a m 11 12 ··· 1 6 a21 a22 a1m 7 A = 6 . ··· . 7 = [aij ]; 4 . 5 · an an anm 1 2 ··· accompanied by their manipulation rules (addition, multiplication, etc). The number of linearly independent rows (or columns) in A is called its rank and denoted rank A. A matrix A Rn×m is said to 2 have full row (column) rank if rank A = n (rank A = m) − be square if n = m, tall if n > m, and fat if n < m − be diagonal (A = diag ai ) if it is square and aij = 0 whenever i = j − f g 6 be lower (upper) triangular if aij = 0 whenever i > j (i < j) − 0 0 be symmetric if A = A , where the transpose A = [aji ] − · The identity matrix I = diag 1 Rn×n (if the dimension is clear, just I ). n · It is a convention that ·A0 = I fforg every2 nonzero square A. Preliminaries: linear algebra State-space models System interconnections Block matrices n×m P P Let A R and let i=1 ni = n and i=1 mi = n for some ni ; mi N.
    [Show full text]
  • 1111: Linear Algebra I
    1111: Linear Algebra I Dr. Vladimir Dotsenko (Vlad) Michaelmas Term 2015 Dr. Vladimir Dotsenko (Vlad) 1111: Linear Algebra I Michaelmas Term 2015 1 / 10 Row expansion of the determinant Our next goal is to prove the following formula: i1 i2 in det(A) = ai1C + ai2C + ··· + ainC ; in other words, if we multiply each element of the i-th row of A by its cofactor and add results, the number we obtain is equal to det(A). Let us first handle the case i = 1. In this case, once we write the corresponding minors with signs, the formula reads 11 12 n+1 1n det(A) = a11A - a12A + ··· + (-1) a1nA : Let us examine the formula for det(A). It is a sum of terms corresponding to pick n elements representing each row and each column of A. In particular, it involves picking an element from the first row. This can be one of the elements a11; a12;:::; a1n. What remains, after we picked the element a1k , is to pick other elements; now we have to avoid the row 1 and the column k. This means that the elements we need to pick are precisely those involved in the determinant A1k , and we just need to check that the signs match. Dr. Vladimir Dotsenko (Vlad) 1111: Linear Algebra I Michaelmas Term 2015 2 / 10 Row expansion of the determinant 1 In fact, it is quite easy to keep track of all signs. The column in the k two-row notation does not add inversions in the first row, and adds k - 1 inversions in the second row, since k appears before 1; 2;:::; k - 1.
    [Show full text]
  • Arxiv:1912.06967V1 [Math.RA] 15 Dec 2019 Values Faui Eigenvector Unit a of Ler Teacher
    MORE ON THE EIGENVECTORS-FROM-EIGENVALUES IDENTITY MALGORZATA STAWISKA Abstract. Using the notion of a higher adjugate of a matrix, we generalize the eigenvector-eigenvalue formula surveyed in [4] to arbitrary square matrices over C and their possibly multiple eigenvalues. 2010 Mathematics Subject Classification. 15A15, 15A24, 15A75 Keywords: Matrices, characteristic polynomial, eigenvalues, eigenvec- tors, adjugates. Dedication: To the memory of Witold Kleiner (1929-2018), my linear algebra teacher. 1. Introduction Very recently the following relation between eigenvectors and eigen- values of a Hermitian matrix (called “the eigenvector-eigenvalue iden- tity”) gained attention: Theorem 1.1. ([4]) If A is an n × n Hermitian matrix with eigen- values λ1(A), ..., λn(A) and i, j = 1, ..., n, then the jth component vi,j of a unit eigenvector vi associated to the eigenvalue λi(A) is related to the eigenvalues λ1(Mj), ..., λn−1(Mj) of the minor Mj of A formed by removing the jth row and column by the formula n n−1 2 |vi,j| Y (λi(A) − λk(A)) = Y(λi(A) − λk(Mj)). arXiv:1912.06967v1 [math.RA] 15 Dec 2019 k=1;k=6 i k=1 The survey paper ([4]) presents several proofs and some generaliza- tions of this identity, along with its many variants encountered in the modern literature, while also commenting on some sociology-of-science aspects of its jump in popularity. It should be pointed out that for sym- metric matrices over R this formula was established already in 1834 by Carl Gustav Jacob Jacobi as formula 30 in his paper [9] (in the process of proving that every real quadratic form has an orthogonal diagonal- ization).
    [Show full text]
  • Systems of 2 Autonomous, Linear, Odes
    In each of the following pairs of problems, we first find a particular solution using the method of variation of parameters using the given fundamental set of solutions fx1; x2g, and then we find a solution of the specified initial value problem. For the method of variation of parameters we use Theorem 4.7.1 in the special case that P (t) = A is a constant matrix. If X(t) is a fundamental matrix for the homogeneous equation x0 = Ax, then Z −1 xp(t) = X(t) X (t)g(t) dt is a particular solution to the inhomogeneous equation x0 = Ax + g(t). −1 1 Note that X (t) = W (t) adj(X(t)), where adj(X(t)) is the adjugate matrix of X(t), and W (t) = det X(t) is the Wronskian. To solve the initial value problem once we have xp(t) we solve for c1 and c2 so that x(0) = c1x1(0) + c2x2(0) + xp(0) equals the prescribed initial value. Problems 2 and 6 −3 1 1 1 1 2 x0 = x+ ; x = e−2t ; x = e−4t ; x(0) = : 1 −3 4t 1 1 2 −1 −1 The Wronskian is W (t) = det X(t) = −2e−6t, so we have Z 1 Z −e−4t −e−4t 1 u(t) = X−1(t)g(t) dt = − e6t dt 2 −e−2t e−2t 4t 1 Z e2t(1 + 4t) 1 e2t(−1 + 4t) = dt = : 2 e4t(1 − 4t) 4 e4t( 1 − 2t) Thus, our particular solution will be 1 −1 + 4t + 1 − 2t 1 t x (t) = X(t)u(t) = = : p 4 −1 + 4t − 1 + 2t 2 −1 + 3t To find the solution of the initial value problem, we need to find c1 and c2 so that 1 1 1 0 2 c + c + = : 1 1 2 −1 2 −1 −1 We must have c1 = 3=4 and c2 = 5=4, so our solution is 3 5 1 3e−2t + 5e−4t + 2t x(t) = x (t) + x (t) + x (t) = : 4 1 4 2 p 4 3e−2t − 5e−4t − 2 + 6t 1 Problems 3 and 7 −1 −2 e−t 1 −1 −1 x0 = x + ; x = e−t ;
    [Show full text]