期刊
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS
卷 31, 期 3, 页码 1279-1302出版社
SIAM PUBLICATIONS
DOI: 10.1137/090757058
关键词
matrix square root; Brouwer's fixed point theorem; Krawczyk's method; Kronecker structures; interval analysis; circular arithmetic
We present methods to compute verified square roots of a square matrix A. Given an approximation X to the square root, obtained by a classical floating point algorithm, we use interval arithmetic to find an interval matrix which is guaranteed to contain the error of X. Our approach is based on the Krawczyk method, which we modify in two different ways in such a manner that the computational complexity for an n x n matrix is reduced to n(3). The methods are based on the spectral decomposition or, in the case that the eigenvector matrix is ill conditioned, on a similarity transformation to block diagonal form. Numerical experiments prove that our methods are computationally efficient and that they yield narrow enclosures provided X is a good approximation. This is particularly true for symmetric matrices, since their eigenvector matrix is perfectly conditioned.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据