{"title":"On certain matrix algebras related to quasi-Toeplitz matrices","authors":"Dario A. Bini, Beatrice Meini","doi":"10.1007/s11075-024-01855-3","DOIUrl":null,"url":null,"abstract":"<p>Let <span>\\(A_\\alpha \\)</span> be the semi-infinite tridiagonal matrix having subdiagonal and superdiagonal unit entries, <span>\\((A_\\alpha )_{11}=\\alpha \\)</span>, where <span>\\(\\alpha \\in \\mathbb C\\)</span>, and zero elsewhere. A basis <span>\\(\\{P_0,P_1,P_2,\\ldots \\}\\)</span> of the linear space <span>\\(\\mathcal {P}_\\alpha \\)</span> spanned by the powers of <span>\\(A_\\alpha \\)</span> is determined, where <span>\\(P_0=I\\)</span>, <span>\\(P_n=T_n+H_n\\)</span>, <span>\\(T_n\\)</span> is the symmetric Toeplitz matrix having ones in the <i>n</i>th super- and sub-diagonal, zeros elsewhere, and <span>\\(H_n\\)</span> is the Hankel matrix with first row <span>\\([\\theta \\alpha ^{n-2}, \\theta \\alpha ^{n-3}, \\ldots , \\theta , \\alpha , 0, \\ldots ]\\)</span>, where <span>\\(\\theta =\\alpha ^2-1\\)</span>. The set <span>\\(\\mathcal {P}_\\alpha \\)</span> is an algebra, and for <span>\\(\\alpha \\in \\{-1,0,1\\}\\)</span>, <span>\\(H_n\\)</span> has only one nonzero anti-diagonal. This fact is exploited to provide a better representation of symmetric quasi-Toeplitz matrices <span>\\(\\mathcal{Q}\\mathcal{T}_S\\)</span>, where, instead of representing a generic matrix <span>\\(A\\in \\mathcal{Q}\\mathcal{T}_S\\)</span> as <span>\\(A=T+K\\)</span>, where <i>T</i> is Toeplitz and <i>K</i> is compact, it is represented as <span>\\(A=P+H\\)</span>, where <span>\\(P\\in \\mathcal {P}_\\alpha \\)</span> and <i>H</i> is compact. It is shown experimentally that the matrix arithmetic obtained this way is much more effective than that implemented in the toolbox <span>CQT-Toolbox</span> of <i>Numer. Algo.</i> 81(2):741–769, 2019.</p>","PeriodicalId":54709,"journal":{"name":"Numerical Algorithms","volume":"32 1","pages":""},"PeriodicalIF":1.7000,"publicationDate":"2024-06-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Numerical Algorithms","FirstCategoryId":"100","ListUrlMain":"https://doi.org/10.1007/s11075-024-01855-3","RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0
Abstract
Let \(A_\alpha \) be the semi-infinite tridiagonal matrix having subdiagonal and superdiagonal unit entries, \((A_\alpha )_{11}=\alpha \), where \(\alpha \in \mathbb C\), and zero elsewhere. A basis \(\{P_0,P_1,P_2,\ldots \}\) of the linear space \(\mathcal {P}_\alpha \) spanned by the powers of \(A_\alpha \) is determined, where \(P_0=I\), \(P_n=T_n+H_n\), \(T_n\) is the symmetric Toeplitz matrix having ones in the nth super- and sub-diagonal, zeros elsewhere, and \(H_n\) is the Hankel matrix with first row \([\theta \alpha ^{n-2}, \theta \alpha ^{n-3}, \ldots , \theta , \alpha , 0, \ldots ]\), where \(\theta =\alpha ^2-1\). The set \(\mathcal {P}_\alpha \) is an algebra, and for \(\alpha \in \{-1,0,1\}\), \(H_n\) has only one nonzero anti-diagonal. This fact is exploited to provide a better representation of symmetric quasi-Toeplitz matrices \(\mathcal{Q}\mathcal{T}_S\), where, instead of representing a generic matrix \(A\in \mathcal{Q}\mathcal{T}_S\) as \(A=T+K\), where T is Toeplitz and K is compact, it is represented as \(A=P+H\), where \(P\in \mathcal {P}_\alpha \) and H is compact. It is shown experimentally that the matrix arithmetic obtained this way is much more effective than that implemented in the toolbox CQT-Toolbox of Numer. Algo. 81(2):741–769, 2019.
期刊介绍:
The journal Numerical Algorithms is devoted to numerical algorithms. It publishes original and review papers on all the aspects of numerical algorithms: new algorithms, theoretical results, implementation, numerical stability, complexity, parallel computing, subroutines, and applications. Papers on computer algebra related to obtaining numerical results will also be considered. It is intended to publish only high quality papers containing material not published elsewhere.