Computing eigenvalues of semi-infinite quasi-Toeplitz matrices

03/12/2022
by   D. A. Bini, et al.
0

A quasi-Toeplitz (QT) matrix is a semi-infinite matrix of the form A=T(a)+E where T(a) is the Toeplitz matrix with entries (T(a))_i,j=a_j-i, for a_j-i∈ℂ, i,j≥ 1, while E is a matrix representing a compact operator in ℓ^2. The matrix A is finitely representable if a_k=0 for k<-m and for k>n, given m,n>0, and if E has a finite number of nonzero entries. The problem of numerically computing eigenpairs of a finitely representable QT matrix is investigated, i.e., pairs (λ, v) such that A v=λ v, with λ∈ℂ, v=(v_j)_j∈ℤ^+, v 0, and ∑_j=1^∞ |v_j|^2<∞. It is shown that the problem is reduced to a finite nonlinear eigenvalue problem of the kind WU(λ)β=0, where W is a constant matrix and U depends on λ and can be given in terms of either a Vandermonde matrix or a companion matrix. Algorithms relying on Newton's method applied to the equation WU(λ)=0 are analyzed. Numerical experiments show the effectiveness of this approach. The algorithms have been included in the CQT-Toolbox [Numer. Algorithms 81 (2019), no. 2, 741–769].

READ FULL TEXT
research
03/27/2023

Algorithms for square root of semi-infinite quasi-Toeplitz M-matrices

A quasi-Toeplitz M-matrix A is an infinite M-matrix that can be written ...
research
02/08/2021

Geometric means of quasi-Toeplitz matrices

We study means of geometric type of quasi-Toeplitz matrices, that are se...
research
07/05/2019

Rational Krylov and ADI iteration for infinite size quasi-Toeplitz matrix equations

We consider a class of linear matrix equations involving semi-infinite m...
research
07/23/2019

Solving quadratic matrix equations arising in random walks in the quarter plane

Quadratic matrix equations of the kind A_1X^2+A_0X+A_-1=X are encountere...
research
12/14/2022

On Eigenvalue Gaps of Integer Matrices

Given an n× n matrix with integer entries in the range [-h,h], how close...
research
07/29/2009

An Explicit Construction of Gauss-Jordan Elimination Matrix

A constructive approach to get the reduced row echelon form of a given m...
research
03/07/2021

Multilevel approximation of Gaussian random fields: Covariance compression, estimation and spatial prediction

Centered Gaussian random fields (GRFs) indexed by compacta such as smoot...

Please sign up or login with your details

Forgot password? Click here to reset