Take the 2-minute tour ×
Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. It's 100% free, no registration required.

Let $A$ be a symmetric positive semidefinite matrix and $I$ the identity matrix.

Given the linear equation

$$ y = A(A + \sigma^2I)^{-1} x $$

Write $A$ in terms of its eigenvectors $|u_i\rangle$, $$ A=\sum_{i=1}^n \lambda_i|u_i\rangle\langle u_i| $$ and assume $$ x = \sum_{i=1}^n \gamma_i |u_i\rangle $$

How can one prove that $$ y = \sum_{i=1}^n \frac{\gamma_i\lambda_i}{\lambda_i + \sigma^2}|u_i\rangle $$

I have been trying to use the matrix inversion lemma, but I can't get the result. Is there something fundamental that I am missing?

Thanks

share|improve this question
    
Given $A$ in the above form, and assuming the eigenvectors are normalized, $(A+\sigma^2I)^{-1}$ is $\sum\limits_{i=1}^n\frac{1}{\lambda_i+\sigma^2}\left|u_i\right>\left<u_i\right|‌​$. –  Peter Morgan Jun 8 '12 at 15:42

1 Answer 1

up vote 2 down vote accepted

First note that \begin{align} (A+\sigma^2 I) &= \sum_i \lambda_i |u_i\rangle \langle u_i| + \sigma^2\sum_i|u_i\rangle \langle u_i|\\ &=\sum_i (\lambda_i +\sigma^2)|u_i\rangle \langle u_i| \end{align} It's diagonal with respect to $|u_i\rangle$, so the inverse is simply,

\begin{align} (A+\sigma^2 I)^{-1} &=\sum_i \frac{1}{(\lambda_i +\sigma^2)}|u_i\rangle \langle u_i| \end{align} Using the vector decomposition $x = \sum_i \gamma_i|u_i\rangle$, we have \begin{align} (A+\sigma^2 I)^{-1}x = \sum_i \frac{\gamma_i}{(\lambda_i +\sigma^2)}|u_i\rangle \end{align} and so \begin{align} A(A+\sigma^2 I)^{-1}x &=\sum_i \frac{\gamma_i}{(\lambda_i +\sigma^2)}A|u_i\rangle\\ &= \sum_i \frac{\gamma_i\lambda_i}{(\lambda_i +\sigma^2)}|u_i\rangle \end{align}

share|improve this answer
    
Gosh! First line is a saver... I did not note that! Thanks Olaf –  JuanPi Jun 8 '12 at 15:50

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.