Bidiagonal Factorization of Banded Recursion Matrices for Mixed-Type Multiple Orthogonal Polynomials
Abstract
Given a banded matrix TN with p subdiagonals and q superdiagonals arising from the Gauss--Borel factorization MN = LN-1UN-1 of a moment matrix, this paper constructs explicitly its bidiagonal factorization \[ TN = L1 ·s Lp\, Uq ·s U1. \] Bidiagonal factorizations of this type are central to the study of oscillatory banded matrices and to the spectral Favard theorem for multiple orthogonal polynomials The factorization is obtained via Christoffel transformations of the moment matrix. Provided that the perturbed moment matrices MN,(b,0) and MN,(0,a) admit a Gauss--Borel factorization, each bidiagonal factor is a quotient of the corresponding Gauss--Borel factors: \[ Ub = UN,(b,0)-1UN,(b-1,0), La = LN,(0,a-1)LN,(0,a)-1. \] Explicit Christoffel-type formulas for the entries of the bidiagonal factors are then derived in terms of certain tau-determinants evaluated at the origin: \[ Ub,n = -τBb-1,n\,τBb,n+1 τBb-1,n+1\,τBb,n, La,n+1 = -τAa-1,n+2\,τAa,n τAa-1,n+1\,τAa,n+1. \] As an illustration, the theory is applied to the recurrence matrices of multiple Hahn orthogonal polynomials. For two weights the tetradiagonal case is handled via contiguous hypergeometric relations; for three weights, i.e. the pentadiagonal case, the direct hypergeometric representations are required. In both cases fully explicit bidiagonal factorizations are obtained.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.