DOCSLIB.ORG

On Some Pentadiagonal Matrices: Their Determinants and Inverses

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
On Some Pentadiagonal Matrices: Their Determinants and Inverses Annales Univ. Sci. Budapest., Sect. Comp. 51 (2020) 3950 ON SOME PENTADIAGONAL MATRICES: THEIR DETERMINANTS AND INVERSES Carlos M. da Fonseca (Safat, Kuwait) László Losonczi (Debrecen, Hungary) Dedicated to the 70th birthday of Professor Antal Járai Communicated by László Szili (Received February 7, 2020; accepted June 2, 2020) Abstract. In this paper we consider pentadiagonal (n+1)×(n+1) matri- ces with two super-diagonals and two sub-diagonals at distances k and 2k from the main diagonal. We give an explicit formula for their determinants provided that (n + 1)=3 ≤ k ≤ n=2: We consider the Toeplitz and imper- fect Toeplitz versions of such matrices, and give explicit formulas for their determinants. We also show that the inverse matrix can be obtained as the product of an upper triangle matrix with two super-diagonals, a diagonal and a lower triangle matrix with two sub-diagonals. 1. Introduction Let n; k; ` be given positive integers with 1 ≤ k < ` ≤ n and denote by Mn the set of n × n complex matrices. Consider the pentadiagonal matrix A0 = (aij) 2 Mn+1 where for i; j = 0; 1; : : : ; n, 8 Lj if j − i = −`; > > lj if j − i = −k; > < di if j − i = 0; aij = if > ri j − i = k; > if > Ri j − i = `; :> 0 otherwise: Key words and phrases: Toeplitz matrix, pentadiagonal, tridiagonal matrices. 2010 Mathematics Subject Classication: Primary 15A15, 15B9, Secondary 65F40, 15A42. 40 C. M. da Fonseca and L. Losonczi Notice that the numbering of entries starts with zero. We also denote the matrix A0 by An+1;k;`(L; l; d; r; R) and denote by D0 = Dn+1;k;`(L; l; d; r; R) its determinant where the diagonal vectors L; l; d; r; R are dened by L = (L0;:::;Ln−`); l = (l0; : : : ; ln−k); d = (d0; : : : ; dn); R = (R0;:::;Rn−`); r = (r0; : : : ; rn−k): We will refer to the diagonal vectors L; l; d; r; R as `th sub-diagonal, kth sub- diagonal, (main) diagonal, kth super-diagonal, `th super-diagonal, respectively. If L; R are zero vectors then our matrix becomes a tridiagonal one denoted by An+1;k(l; d; r) and its determinant by Dn+1;k(l; d; r). If 1 ≤ k = ` ≤ n, then the two super-diagonals and the two sub-diagonals slip together and our matrix becomes a tridiagonal matrix An+1;k;k(L; l; d; r; R) = An+1;k(L + l; d; r + R) and its determinant D0 = Dn+1;k(L + l; d; r + R): We shall call An+1;k;`(L; l; d; r; R) a (general) k; `-pentadiagonal matrix while An+1;k(l; d; r) will be termed as k-tridiagonal. In [6] we developed a method to reduce the determinant of k; `-pentadiagonal matrices to tridiagonal determinants provided that k + ` ≥ n + 1. If k ≥ ≥ (n + 1)=3, then by suitable matrix multiplications we can reduce a k; 2k- pentadiagonal matrix to a tridiagonal one. A similar reduction of tridagonals leads to diagonal matrices. In this way we get an explicit formula for the determinant of a general k; 2k-pentadiagonal matrix. In case of Toeplitz pentadiagonal matrices the diagonal vectors are con- stant vectors, i.e., L = (L; : : : ; L), R = (R; : : : ; R) 2 Rn+1−`, l = (l; : : : ; l), r = (r; : : : ; r) 2 Rn+1−k, d = (d; : : : ; d) 2 Rn+1 and for the matrix and its determinant the notations An+1;k;`(L; l; d; r; R) and Dn+1;k;`(L; l; d; r; R) will be used. For imperfect Toeplitz matrices (the term used by Marr and Vineyard [9]) the main diagonal is n+1 (d − α; : : : ; d − α; d; : : : ; d; d − β; : : : ; d − β) 2 R | {z } | {z } | {z } k n+1−2k k while the other diagonals are the same as above. To indicate the imperfectness we add a superscript (α,β) to the notation of Toeplitz pentadiagonal matrices. Thus an imperfect -pentadiagonal matrix (α,β) has the k; 2k An+1;k;2k(L; l; d; r; R) form On some pentadiagonal matrices 41 0 d − α r R 1 . B .. C B d − α r C B .. .. C B . R C B C B d − α C B C B .. C B l d . C B C B l d C B . C B .. .. .. C B C B d r C B C B .. C B . d − β C B C B L d − β C B . C @ .. .. .. A L l d − β As far as we can tell, the rst particular k; `-pentadiagonal Toeplitz matrix was studied in 1928 by Jen® Egerváry and Otto Szász [3]. One of the matrices they ¯ considered was An+1;k;`(δ; ,¯ −λ, , δ) with k+` = n+1 and jj = jδj = 1. Their result has been largely ignored in the pure and applied matrix theory as we can see e.g. in [2, 10]. For some other papers on the topic, the reader is refereed to [1,6, 11, 12]. An interesting graph theoretical approach can be found in [7,5]. An exhaustive list of recent references is given in the survey [4]. Our paper was inspired by Marr and Vineyard [9] who have shown that the product of two 1-tridiagonal Toeplitz matrices is an imperfect Toeplitz matrix (α,β) which is related to the corresponding Toeplitz matrix by a two-step An+1;1;2 recursion. In this way they succeeded to express the determinant Dn+1;1;2 of 1; 2-pentadiagonal Toeplitz matrix in terms of Chebyshev polynomials of second kind. A similar approach was used in [13] to nd a formula for the inverse of a 1; 2-pentadiagonal Toeplitz matrix. 2. Reduction of general k; 2k-pentadiagonal determinants to tridiagonal ones Our starting point is the matrix A0 = An+1;k;2k(L; l; d; r; R) where 1 ≤ k < < 2k ≤ n and k + 2k ≤ n + 1, i.e., (n + 1)=3 ≤ k ≤ n=2 < (n + 1)=2. Let us dene the matrices B = (bij);C = (cij),F = (fij);G = (gij) of Mn+1 by 42 C. M. da Fonseca and L. Losonczi 8 1 if i=j; > <> lj b = − if i=j+k; j = 0; : : : ; k−1; ij d > j > : 0 otherwise; 8 1 if i=j; > > r < − i if i=0; : : : ; k−1; j =i+k; cij = d > i > : 0 otherwise; 8 1 if i=j; > <> Lj f = − if i=j+2k; j =0; : : : ; n−2k; ij d > j > : 0 otherwise; 8 1 if i=j; > <> Ri g = − if i=0; : : : ; n−2k; j =i+2k; ij d > i > : 0 otherwise: Performing four matrix multiplications (for some details see [6], Theo- rem 3, (i)) we get the matrix L ∗ where A1 = F BA0CG = E1 A1 E1 2 Mk is a diagonal matrix with diagonal elements d0; : : : ; dk−1 and ∗ (1) (1) (1) (1) (1) (1) (1) (1) (1) (1) A1 =An+1−k;k;k(L ; l ; d ; r ;R )=An+1−k;k(L +l ; d ; r +R ) is tridiagonal, where the iterated diagonal vectors are L r L r L(1) = − 0 0 ;:::; − n−2k n−2k ; d0 dn−2k (1) l = (lk; : : : ; ln−k) ; (1) l0r0 lk−1rk−1 L0R0 Ln−2kRn−2k d = dk − ; : : : ; d2k−1 − ; d2k − ; : : : ; dn − ; d0 dk−1 d0 dn−2k (1) r = (rk; : : : ; rn−k) ; R l R l R(1) = − 0 0 ;:::; − n−2k n−2k : d0 dn−2k On some pentadiagonal matrices 43 3. Further reduction to diagonal matrix The matrix A1 can be reduced to a diagonal matrix with two matrix mul- tiplications. Dene the matrices U = (uij);V = (vij) 2 Mn+1 by 8 1 if i = j; > > > < lj − Lj−krj−k=dj−k uij = − if i = j + k; j = k; : : : ; n − k; > dj − lj−krj−k=dj−k > > :> 0 otherwise; 8 1 if i = j; > > > < ri − Ri−kli−k=di−k vij = − if i = k; : : : ; n − k; j = i + k; > dj − li−kri−k=di−k > > :> 0 otherwise. Multiplying A1 from the left by U, then from the right by V we get a diagonal matrix A2 = UA1V with diagonal (3.1) d0; : : : ; dk−1; l0r0 lk−1rk−1 dk − ; : : : ; d2k−1 − ; d0 dk−1 L0r0 R0l0 lk − rk − L0R0 d0 d0 d2k − − ; l0r0 d0 dk − d0 Ln−2krn−2k Rn−2kln−2k ln−k − rn−k − Ln−2kRn−2k dn−2k dn−2k : : : ; dn − − ; ln−2krn−2k dn−2k dn−k − dn−2k where in the above lines k; k; n + 1 − 2k elements are listed. Because the determinants of the matrices B; C; F; G; U; V all are 10s, the determinant of our general pentadiagonal matrix An+1;k;2k is the product of the diagonal elements of A2: To calculate this product we multiply the elements of the rst two groups pairwise, to obtain d0dk − l0r0; : : : ; dk−1d2k−1 − lk−1rk−1: 44 C. M. da Fonseca and L. Losonczi Simplifying the general term of the third group we get, for j = 0; : : : ; n − 2k, Lj rj Rjlj lk+j − d rk+j − LjRj j dj d2k+j − − = dj ljrj dk+j − dj d d d −d l r −L R d −l r d +l R l +r L r = j j+k j+2k j j+k j+k j j j+k j j j+2k j j j+k j j j+k : djdj+k−ljrj Thus we have proved Theorem 3.1. Assuming (n+1)=3≤ k ≤n=2 the determinant Dn+1;k;2k(L; l; d; r; R) of the general k; 2k-pentadiagonal matrix is (3.2) n−2k Q (djdj+kdj+2k −djlj+krj+k −LjRjdj+k −ljrjdj+2k +ljRjlj+k +rjLjrj+k)· j=0 k−1 Q · (djdj+k−ljrj) j=n+1−2k Remark 3.1.
Load more

Recommended publications
  • Solving Systems of Linear Equations by Gaussian Elimination
    Chapter 3 Solving Systems of Linear Equations By Gaussian Elimination 3.1 Mathematical Preliminaries In this chapter we consider the problem of computing the solution of a system of n linear equations in n unknowns. The scalar form of that system is as follows: a11x1 +a12x2 +... +... +a1nxn = b1 a x +a x +... +... +a x = b (S) 8 21 1 22 2 2n n 2 > ... ... ... ... ... <> an1x1 +an2x2 +... +... +annxn = bn > Written in matrix:> form, (S) is equivalent to: (3.1) Ax = b, where the coefficient square matrix A Rn,n, and the column vectors x, b n,1 n 2 2 R ⇠= R . Specifically, a11 a12 ... ... a1n a21 a22 ... ... a2n A = 0 1 ... ... ... ... ... B a a ... ... a C B n1 n2 nn C @ A 93 94 N. Nassif and D. Fayyad x1 b1 x2 b2 x = 0 1 and b = 0 1 . ... ... B x C B b C B n C B n C @ A @ A We assume that the basic linear algebra property for systems of linear equa- tions like (3.1) are satisfied. Specifically: Proposition 3.1. The following statements are equivalent: 1. System (3.1) has a unique solution. 2. det(A) =0. 6 3. A is invertible. In this chapter, our objective is to present the basic ideas of a linear system solver. It consists of two main procedures allowing to solve efficiently (3.1). 1. The first, referred to as Gauss elimination (or reduction) reduces (3.1) into an equivalent system of linear equations, which matrix is upper triangular. Specifically one shows in section 4 that Ax = b Ux = c, () where c Rn and U Rn,n is given by: 2 2 u11 u12 ..
    [Show full text]
  • Arxiv:2009.05100V2
    THE COMPLETE POSITIVITY OF SYMMETRIC TRIDIAGONAL AND PENTADIAGONAL MATRICES LEI CAO 1,2, DARIAN MCLAREN 3, AND SARAH PLOSKER 3 Abstract. We provide a decomposition that is sufficient in showing when a symmetric tridiagonal matrix A is completely positive. Our decomposition can be applied to a wide range of matrices. We give alternate proofs for a number of related results found in the literature in a simple, straightforward manner. We show that the cp-rank of any irreducible tridiagonal doubly stochastic matrix is equal to its rank. We then consider symmetric pentadiagonal matrices, proving some analogous results, and providing two different decom- positions sufficient for complete positivity. We illustrate our constructions with a number of examples. 1. Preliminaries All matrices herein will be real-valued. Let A be an n n symmetric tridiagonal matrix: × a1 b1 b1 a2 b2 . .. .. .. . A = .. .. .. . bn 3 an 2 bn 2 − − − bn 2 an 1 bn 1 − − − bn 1 an − We are often interested in the case where A is also doubly stochastic, in which case we have ai = 1 bi 1 bi for i = 1, 2,...,n, with the convention that b0 = bn = 0. It is easy to see that− if a− tridiagonal− matrix is doubly stochastic, it must be symmetric, so the additional hypothesis of symmetry can be dropped in that case. We are interested in positivity conditions for symmetric tridiagonal and pentadiagonal matrices. A stronger condition than positive semidefiniteness, known as complete positivity, arXiv:2009.05100v2 [math.CO] 10 Mar 2021 has applications in a variety of areas of study, including block designs, maximin efficiency- robust tests, modelling DNA evolution, and more [5, Chapter 2], as well as recent use in mathematical optimization and quantum information theory (see [14] and the references therein).
    [Show full text]
  • Inverse Eigenvalue Problems Involving Multiple Spectra
    Inverse eigenvalue problems involving multiple spectra G.M.L. Gladwell Department of Civil Engineering University of Waterloo Waterloo, Ontario, Canada N2L 3G1 [email protected] URL: http://www.civil.uwaterloo.ca/ggladwell Abstract If A Mn, its spectrum is denoted by σ(A).IfA is oscillatory (O) then σ(A∈) is positive and discrete, the submatrix A[r +1,...,n] is O and itsspectrumisdenotedbyσr(A). Itisknownthatthereisaunique symmetric tridiagonal O matrix with given, positive, strictly interlacing spectra σ0, σ1. It is shown that there is not necessarily a pentadiagonal O matrix with given, positive strictly interlacing spectra σ0, σ1, σ2, but that there is a family of such matrices with positive strictly interlacing spectra σ0, σ1. The concept of inner total positivity (ITP) is introduced, and it is shown that an ITP matrix may be reduced to ITP band form, or filled in to become TP. These reductions and filling-in procedures are used to construct ITP matrices with given multiple common spectra. 1Introduction My interest in inverse eigenvalue problems (IEP) stems from the fact that they appear in inverse vibration problems, see [7]. In these problems the matrices that appear are usually symmetric; in this paper we shall often assume that the matrices are symmetric: A Sn. ∈ If A Sn, its eigenvalues are real; we denote its spectrum by σ(A)= ∈ λ1, λ2,...,λn ,whereλ1 λ2 λn. The direct problem of finding σ{(A) from A is} well understood.≤ At≤ fi···rst≤ sight it appears that inverse eigenvalue T problems are trivial: every A Sn with spectrum σ(A) has the form Q Q ∈ ∧ where Q is orthogonal and = diag(λ1, λ2,...,λn).
    [Show full text]
  • The Rule of Hessenberg Matrix for Computing the Determinant of Centrosymmetric Matrices
    CAUCHY –Jurnal Matematika Murni dan Aplikasi Volume 6(3) (2020), Pages 140-148 p-ISSN: 2086-0382; e-ISSN: 2477-3344 The Rule of Hessenberg Matrix for Computing the Determinant of Centrosymmetric Matrices Nur Khasanah1, Agustin Absari Wahyu Kuntarini 2 1,2Department of Mathematics, Faculty of Science and Technology UIN Walisongo Semarang Email: [email protected] ABSTRACT The application of centrosymmetric matrix on engineering takes their part, particularly about determinant rule. This basic rule needs a computational process for determining the appropriate algorithm. Therefore, the algorithm of the determinant kind of Hessenberg matrix is used for computing the determinant of the centrosymmetric matrix more efficiently. This paper shows the algorithm of lower Hessenberg and sparse Hessenberg matrix to construct the efficient algorithm of the determinant of a centrosymmetric matrix. Using the special structure of a centrosymmetric matrix, the algorithm of this determinant is useful for their characteristics. Key Words : Hessenberg; Determinant; Centrosymmetric INTRODUCTION One of the widely used studies in the use of centrosymmetric matrices is how to get determinant from centrosymmetric matrices. Besides this special matrix has some applications [1], it also has some properties used for determinant purpose [2]. Special characteristic centrosymmetric at this entry is evaluated at [3] resulting in the algorithm of centrosymmetric matrix at determinant. Due to sparse structure of this entry, the evaluation of the determinant matrix has simpler operations than full matrix entries. One special sparse matrix having rules on numerical analysis and arise at centrosymmetric the determinant matrix is the Hessenberg matrix. The role of Hessenberg matrix decomposition is the important role of computing the eigenvalue matrix.
    [Show full text]
  • MAT TRIAD 2019 Book of Abstracts
    MAT TRIAD 2019 International Conference on Matrix Analysis and its Applications Book of Abstracts September 8 13, 2019 Liblice, Czech Republic MAT TRIAD 2019 is organized and supported by MAT TRIAD 2019 Edited by Jan Bok, Computer Science Institute of Charles University, Prague David Hartman, Institute of Computer Science, Czech Academy of Sciences, Prague Milan Hladík, Department of Applied Mathematics, Charles University, Prague Miroslav Rozloºník, Institute of Mathematics, Czech Academy of Sciences, Prague Published as IUUK-ITI Series 2019-676ø by Institute for Theoretical Computer Science, Faculty of Mathematics and Physics, Charles University Malostranské nám. 25, 118 00 Prague 1, Czech Republic Published by MATFYZPRESS, Publishing House of the Faculty of Mathematics and Physics, Charles University in Prague Sokolovská 83, 186 75 Prague 8, Czech Republic Cover art c J. Na£eradský, J. Ne²et°il c Jan Bok, David Hartman, Milan Hladík, Miroslav Rozloºník (eds.) c MATFYZPRESS, Publishing House of the Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic, 2019 i Preface This volume contains the Book of abstracts of the 8th International Conference on Matrix Anal- ysis and its Applications, MAT TRIAD 2019. The MATTRIAD conferences represent a platform for researchers in a variety of aspects of matrix analysis and its interdisciplinary applications to meet and share interests and ideas. The conference topics include matrix and operator theory and computation, spectral problems, applications of linear algebra in statistics, statistical models, matrices and graphs as well as combinatorial matrix theory and others. The goal of this event is to encourage further growth of matrix analysis research including its possible extension to other elds and domains.
    [Show full text]
  • Some Bivariate Stochastic Models Arising from Group Representation Theory
    Available online at www.sciencedirect.com ScienceDirect Stochastic Processes and their Applications ( ) – www.elsevier.com/locate/spa Some bivariate stochastic models arising from group representation theory Manuel D. de la Iglesiaa,∗, Pablo Románb a Instituto de Matemáticas, Universidad Nacional Autónoma de México, Circuito Exterior, C.U., 04510, Ciudad de México, Mexico b CIEM, FaMAF, Universidad Nacional de Córdoba, Medina Allende s/n Ciudad Universitaria, Córdoba, Argentina Received 14 September 2016; accepted 31 October 2017 Available online xxxx Abstract The aim of this paper is to study some continuous-time bivariate Markov processes arising from group representation theory. The first component (level) can be either discrete (quasi-birth-and-death processes) or continuous (switching diffusion processes), while the second component (phase) will always be discrete and finite. The infinitesimal operators of these processes will be now matrix-valued (eithera block tridiagonal matrix or a matrix-valued second-order differential operator). The matrix-valued spherical functions associated to the compact symmetric pair (SU(2) × SU(2); diag SU(2)) will be eigenfunctions of these infinitesimal operators, so we can perform spectral analysis and study directly some probabilistic aspects of these processes. Among the models we study there will be rational extensions of the one-server queue and Wright–Fisher models involving only mutation effects. ⃝c 2017 Elsevier B.V. All rights reserved. MSC: 60J10; 60J60; 33C45; 42C05 Keywords: Quasi-birth-and-death processes; Switching diffusions; Matrix-valued orthogonal polynomials; Wright–Fisher models 1. Introduction It is very well known that many important results of one-dimensional stochastic processes can be obtained by using spectral methods.
    [Show full text]
  • CS321 Numerical Analysis
    CS321 Numerical Analysis Lecture 5 System of Linear Equations Professor Jun Zhang Department of Computer Science University of Kentucky Lexington, KY 40506-0046 System of Linear Equations a11 x1 a12 x2 a1n xn b1 a21 x1 a22 x2 a2n xn b2 an1 x1 an2 x2 ann xn bn where aij are coefficients, xi are unknowns, and bi are right-hand sides. Written in a compact form is n aij x j bi , i 1,,n j1 The system can also be written in a matrix form Ax b where the matrix is a11 a12 a1n a a a A 21 22 2n an1 an2 ann and x [x , x ,, x ]T ,b [b ,b ,b ]T 1 2 n 1 2 n 2 An Upper Triangular System An upper triangular system a11 x1 a12 x2 a13 x3 a1n xn b1 a22 x2 a23 x3 a2n xn b2 a33 x3 a3n xn b3 an1,n1 xn1 an1,n xn bn1 ann xn bn is much easier to find the solution: bn xn ann from the last equation and substitute its value in other equations and repeat the process n 1 xi bi aij x j aii ji1 for i = n – 1, n – 2,…, 1 3 Karl Friedrich Gauss (April 30, 1777 – February 23, 1855) German Mathematician and Scientist 4 Gaussian Elimination Linear systems are solved by Gaussian elimination, which involves repeated procedure of multiplying a row by a number and adding it to another row to eliminate a certain variable For a particular step, this amounts to aik aij aij akj (k j n) akk aik bi bi bk akk th After this step, the variable xk, is eliminated in the (k + 1) and in the later equations The Gaussian elimination modifies a matrix into an upper triangular form such that aij = 0 for all i > j.
    [Show full text]
  • Quantum Phase Transitions Mediated by Clustered Non-Hermitian
    Quantum phase transitions mediated by clustered non-Hermitian degeneracies Miloslav Znojil The Czech Academy of Sciences, Nuclear Physics Institute, Hlavn´ı130, 250 68 Reˇz,ˇ Czech Republic and Department of Physics, Faculty of Science, University of Hradec Kr´alov´e, Rokitansk´eho 62, 50003 Hradec Kr´alov´e, Czech Republic e-mail: [email protected] Abstract The phenomenon of degeneracy of an N plet of bound states is studied in the framework of the − quasi-Hermitian (a.k.a. symmetric) formulation of quantum theory of closed systems. For PT − a general non-Hermitian Hamiltonian H = H(λ) such a degeneracy may occur at a real Kato’s exceptional point λ(EPN) of order N and of the geometric multiplicity alias clusterization index K. The corresponding unitary process of collapse (loss of observability) can be then interpreted as a generic quantum phase transition. The dedicated literature deals, predominantly, with the non-numerical benchmark models of the simplest processes where K = 1. In our present paper it is shown that in the “anomalous” dynamical scenarios with 1 <K N/2 an analogous approach ≤ is applicable. A multiparametric anharmonic-oscillator-type exemplification of such systems is constructed as a set of real-matrix N by N Hamiltonians which are exactly solvable, maximally non-Hermitian and labeled by specific ad hoc partitionings (N) of N. R arXiv:2102.12272v1 [quant-ph] 24 Feb 2021 1 1 Introduction The experimentally highly relevant phenomenon of a quantum phase transition during which at least one of the observables loses its observability status is theoretically elusive. In the conventional non-relativistic Schr¨odinger picture, for example, the observability of the energy represented by a self-adjoint Hamiltonian H = H(λ) appears too robust for the purpose.
    [Show full text]
  • Arxiv:2108.07942V1 [Math.CO]
    ALGORITHMIC TECHNIQUES FOR FINDING RESISTANCE DISTANCES ON STRUCTURED GRAPHS E. J. EVANS AND A. E. FRANCIS Abstract. In this paper we give a survey of methods used to calculate values of resistance distance (also known as effective resistance) in graphs. Resistance distance has played a prominent role not only in circuit theory and chemistry, but also in combinatorial matrix theory and spectral graph theory. Moreover resistance distance has applications ranging from quantifying biological struc- tures, distributed control systems, network analysis, and power grid systems. In this paper we discuss both exact techniques and approximate techniques and for each method discussed we provide an illustrative example of the technique. We also present some open questions and conjectures. 1. Introduction The resistance distance (occasionally referred to as the effective resistance) of a graph is a measure that quantifies its structural properties. Resistance distance has its origin in electrical circuit theory and its first known application to graph structure occurred in the analysis of chemical structure [22]. Resistance distance in graphs has played a prominent role not only in circuit theory and chemistry [18, 22, 30], but also in combinatorial matrix theory [4, 36] and spectral graph theory [2, 10, 13, 29]. A few specific examples of the use of resistance distance are: Spielman and Srivastava [29] have used resistance distance between nodes • of graphs to develop an algorithm to rapidly sparsify a given graph while maintaining spectral properties. Ghosh, Boyd, and Saberi [21] considered the problem of minimizing the • total resistance distance by allocating edge weights on a given graph. This problem has applications to Markov chains and continuous-time averaging arXiv:2108.07942v2 [math.CO] 13 Sep 2021 networks.
    [Show full text]
  • DETERMINANTS of MULTIDIAGONAL MATRICES∗ 1. Introduction. the First Purpose of This Paper Is to Present and Discuss Algo- Rithm
    Electronic Journal of Linear Algebra ISSN 1081-3810 A publication of the International Linear Algebra Society Volume 25, pp. 102-118, November 2012 ELA DETERMINANTS OF MULTIDIAGONAL MATRICES∗ KATARZYNA FILIPIAK† , AUGUSTYN MARKIEWICZ† , AND ANETA SAWIKOWSKA‡ Abstract. The formulas presented in [Molinari, L.G. Determinants of block tridiagonal matri- ces. Linear Algebra Appl., 2008; 429, 2221–2226] for evaluating the determinant of block tridiagonal matrices with (or without) corners are used to derive the determinant of any multidiagonal matri- ces with (or without) corners with some specified non-zero minors. Algorithms for calculation the determinant based on this method are given and properties of the determinants are studied. Some applications are presented. Key words. Multidiagonal matrix with corners, Multidiagonal matrix without corners, Deter- minant. AMS subject classifications. 15A15. 1. Introduction. The first purpose of this paper is to present and discuss algo- rithms for calculation of the determinant of multidiagonal matrices with (or without) corners. The second purpose is to compare them with some other algorithms in the case of pentadiagonal matrix with (or without) corners. We consider the determinant of a multidiagonal matrix of order n with and (p) (p) without corners, i.e., matrices A = (ai,j ) and B = (bi,j ) with ai,j = 0 if p 1 p 1 p 1 − < i j < n − and b = 0 if i j > − , where p (odd) denotes 2 | − | − 2 i,j | − | 2 the number of diagonals. In particular, tridiagonal and pentadiagonal matrices with corners are of the form: ∗Received by the editors on October 25, 2011. Accepted for publication on November 3, 2012 Handling Editor Natalia Bebiano.
    [Show full text]
  • A Numerical Algorithm for Computing the Inverse of a Toeplitz Pentadiagonal Matrix
    Journal of Applied Mathematics and Computational Mechanics 2018, 17(3), 83-95 www.amcm.pcz.pl p-ISSN 2299-9965 DOI: 10.17512/jamcm.2018.3.08 e-ISSN 2353-0588 A NUMERICAL ALGORITHM FOR COMPUTING THE INVERSE OF A TOEPLITZ PENTADIAGONAL MATRIX 1 2 1 Boutaina Talibi , Ahmed Driss Aiat Hadj , Driss Sarsri 1 National School of Applied Sciences, Abdelmalek Essaadi University Tangier, Morocco 2 Regional Center of the Trades of Education and Training (CRMEF)-Tangier, Avenue My Abdelaziz Souani, BP: 3117, Tangier, Morocco [email protected], [email protected], [email protected] Received: 2 August 2018; Accepted: 27 November 2018 Abstract. In the current paper, we present a computationally efficient algorithm for obtain- ing the inverse of a pentadiogonal toeplitz matrix. Few conditions are required, and the al- gorithm is suited for implementation using computer algebra systems. MSC 2010: 15B05 65F40 33C45 Keywords: pentadiagonal matrix, Toeplitz matrix, K-Hessenberg matrix, inverse 1. Introduction Pentadiagonal Toeplitz matrices often occur when solving partial differential equations numerically using the finite difference method, the finite element method and the spectral method, and could be applied to the mathematical representation of high dimensional, nonlinear electromagnetic interference signals [1-7]. Pentadiagonal matrices are a certain class of special matrices, and other common types of special matrices are Jordan, Frobenius, generalized Vandermonde, Hermite, centrosymmetric, and arrowhead matrices [8-12]. Typically, after the original partial differential equations are processed with these numerical methods, we need to solve the pentadiagonal Toeplitz systems of linear equations for obtaining the numerical solutions of the original partial differ- ential equations.
    [Show full text]
  • Fiedler-Comrade and Fiedler-Chebyshev Pencils 1
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by University of Essex Research Repository FIEDLER-COMRADE AND FIEDLER-CHEBYSHEV PENCILS VANNI NOFERINI† ¶ AND JAVIER PEREZ´ ‡¶ Abstract. Fiedler pencils are a family of strong linearizations for polynomials expressed in the monomial basis, that include the classical Frobenius companion pencils as special cases. We generalize the definition of a Fiedler pencil from monomials to a larger class of orthogonal polynomial bases. In particular, we derive Fiedler-comrade pencils for two bases that are extremely important in practical applications: the Chebyshev polynomials of the first and second kind. The new approach allows one to construct linearizations having limited bandwidth: a Chebyshev analogue of the pentadiagonal Fiedler pencils in the monomial basis. Moreover, our theory allows for linearizations of square matrix polynomials expressed in the Chebyshev basis (and in other bases), regardless of whether the matrix polynomial is regular or singular, and for recovery formulae for eigenvectors, and minimal indices and bases. Key words. Fiedler pencil, Chebyshev polynomial, linearization, matrix polynomial, singular matrix polynomial, eigenvector recovery, minimal basis, minimal indices AMS subject classifications. 15A22, 15A18, 15A23, 65H04, 65F15 1. Motivation. In computational mathematics, many applications require to compute the roots of a polynomial expressed in a nonstandard basis. Particularly relevant in practice [28] are the Chebyshev polynomials of the first kind, that we denote by Tk(x), see (3.1) for their formal definition. For example, suppose that we want to approximate numerically the roots of the polynomial T5(x) − 4T4(x)+4T2(x) − T1(x) (1.1) The roots of (1.1) are easy to compute analytically and they are ±1/2, ±1, and 2.
    [Show full text]