4
$\begingroup$

Is there some simple upper bound on $||(B^{-1}+A^{-1})^{-1}||$, where $A,B$ are $n \times n$ symmetric matrices?

$\endgroup$
4
  • 3
    $\begingroup$ What norm are you denoting by $||\cdot ||$? $\endgroup$
    – Noah Stein
    Commented May 21, 2012 at 18:22
  • $\begingroup$ Of course, for posdef matrices suchs bounds follow trivially using the operator harmonic, geometric, arithmetic mean inequalities. So I guess, the broader class of symmetric or Hermitian matrices in your question is deliberate.... $\endgroup$
    – Suvrit
    Commented May 21, 2012 at 18:41
  • 1
    $\begingroup$ My conjecture is that for any symmetric norm, your lhs $\le \| |a^{-1}b^{-1}|^{-1/2} \|$, where $|x|$ denotes the operator absolute value. Corresponding to this, one can then derive (i guess) an arithmetic-mean-style upper bound. $\endgroup$
    – Suvrit
    Commented May 21, 2012 at 18:51
  • $\begingroup$ (I threw in a factor of $2$ into your lhs to replace it by the more symmetric harmonic mean, while writing out the above inequality) $\endgroup$
    – Suvrit
    Commented May 21, 2012 at 18:51

3 Answers 3

7
$\begingroup$

You can use the surprising identity $(A^{-1}+B^{-1})^{-1}=A(A+B)^{-1}B$, and take the norms of the three factors separately.

$\endgroup$
2
  • $\begingroup$ Depending on the norm you are using, this might not give a very tight bound. Suppose $A=B$ and has one eigenvalue $\epsilon$ and the rest $1$. Then the spectral norm of $A$ is $1$, the spectral norm of $(A+B)^{-1})=1/(2\epsilon$, and the spectral norm of $B$ is $1$, while the actual spectral norm of $(A^{-1}+B^{-1})^{-1}$ is $1/2$. But if $\epsilon <<1$, then $1/(2\epsilon)<<1/2$. $\endgroup$
    – Will Sawin
    Commented May 21, 2012 at 18:45
  • $\begingroup$ I'm afraid this throws me back to the Woodbury formula, which I was coming from, but it's still a nice answer! $\endgroup$ Commented May 21, 2012 at 19:41
2
$\begingroup$

To expand on my first comment, if $A, B > 0$ are symmetric positive definite matrices. Then, it is known that

$$\left(\frac{A^{-1}+B^{-1}}{2}\right)^{-1} \le A\sharp B \le \frac{A+B}{2},$$ where the inequalities are in the Löwner partial order, and $A\sharp B := A^{1/2}(A^{1/2}BA^{1/2})^{-1/2}A^{1/2}$ denotes the matrix geometric mean.

These operator inequalities are of course, stronger than corresponding norm inequalities (based on unitarily invariant norms).

For the case where you don't have positive matrices, I think the conjecture mentioned in my second argument can be expanded into a proof --- maybe if I get time, I'll try to expand that.

$\endgroup$
0
2
$\begingroup$

Probably not unless $A$ and $B$ are positive-definite, since if $B$ is very close to $-A$ then $B^{-1}+A^{-1}$ is very small and so its inverse is very large. In fact, depending on the norm, they probably need to be close only on one shared or almost-shared eigenvector.

For spectral norm of positive-definite matrices, we have a nice answer. The highest eigenvalue of $(A^{-1}+B^{-1})^{-1}$ is the lowest eigenvalue of $A^{-1}+B^{-1}$, which one can find by minimizing $x^T(A^{-1}+B^{-1})x$ with respect to $x^Tx=1$. But the minimum for $A^{-1}$ is its lowest eigenvalue, $1/||A||$, and the minimum for $B^{-1}$ is similarly $1/||B||$. Thus:

$x^T( A^{-1}+B^{-1})^{-1} x= x^T A^{-1} x+ x^T B^{-1} x\geq 1/||A||+1/||B||$

So the spectral norm of the harmonic sum is bounded by the harmonic sum of the spectral norms!

$\endgroup$
1
  • $\begingroup$ "So the spectral norm of the harmonic sum is bounded by the harmonic sum of the spectral norms!" Cosmic harmony... :) $\endgroup$ Commented May 21, 2012 at 19:39

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .