0
$\begingroup$

I'm working on a problem and was lead to trying to find an approximation for the square root of a matrix. I came across a way of doing this using holomorphic functional calculus. However, my first attempt at a computation resulted in failure. Did I miss something? --

Wikipedia gives the expansion of the square root matrix as follows: $$\frac{A^{1/2}}{\|A\|^{1/2}} = I - \sum_{n = 1}^{\infty}\left\lvert {1/2 \choose n}\right\rvert \left(I-\frac{A}{\|A\|}\right)^n = I - \frac 1 2\left(I-\frac{A}{\|A\|}\right) - \frac 1 8 \left(I-\frac{A}{\|A\|}\right)^2 \,...$$ In this post, I'm taking the matrix norm to be the max of the absolute values of the column-sums of matrix $A$.

I thought I might be able to get an estimate of square root matrix by taking the first two terms and using the third term as an error bound. But a simple numerical example appears to not work. I'm not sure if I'm using the expansion incorrectly, or there's some hypothesis I'm missing. In particular for $A = \begin{pmatrix} 4 & 0 \\ 0 & 16 \end{pmatrix}, \|A\| = 16$.

I get an error:$$\big\|A^{1/2}/\|A\|^{1/2} - I + (I-A/\|A\|)/2 \big\| = 1/8.$$ Error bound using the third term: $$\big\|(I-A/\|A\|)^2/8 \big\| = 9/128.$$

[Remark: This is cross posted.]

$\endgroup$
9
  • 2
    $\begingroup$ the Wikipedia formula you cite says that it needs the eigenvalues of $A$ to lie in the interval (0,2), which does not apply to your choice of $A$. $\endgroup$ Sep 26 '20 at 14:22
  • 2
    $\begingroup$ doesn't it state the requirement that the spectrum of $A$ is contained within a unit disc in the complex plane? $\endgroup$ Sep 26 '20 at 15:08
  • 1
    $\begingroup$ As another approach, see mathoverflow.net/questions/301272/… $\endgroup$
    – Mahdi
    Sep 26 '20 at 19:03
  • 1
    $\begingroup$ Why do you expect the first omitted term to bound the error ? This series does not have alternating signs, and as far as I understand this result is not true even for scalar inputs. $\endgroup$ Sep 28 '20 at 14:14
  • 1
    $\begingroup$ No, the $n+1$st term does not estimate the error in Taylor series, in general. Case in point: $1/(1-x)$. $\endgroup$ Sep 29 '20 at 10:02

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.