Eigenvalues of Matrix Addition

[vkillion]

Hello,

I have a linear algebra problem that I need help with.

Basically, I need to get the eigenvalues and eigenvectors of several (sometimes tens of thousands) very large matrices (6^n x 6^n, where n>= 3, to be specific). Currently, we are just using MATLAB's eig() function to get them. I am trying to find optimizations for the simulations to cut down on computing time. There are three matrices that we use.

H_constant - generated before the loop. Real and symmetric about the diagonal. Does not change after initial calculation.

H_location - generated during each iteration. Diagonal.

H_final - H_constant + H_location. Therefore, it is also real and symmetric about the diagonal.

It is H_final that we need the eigenvalues and eigenvectors of. My theory is that we calculate the eigenvalues and eigenvectors of H_constant (which won't change after the initial calculation) once. We use this result with the eigenvalues of H_location (the diagonal), to get the eigenvalues and eigenvectors of H_final1. This would reduce our computation from tens of thousands of eig() calls to 1 eig() call and tens of thousands of very simple operations. I don't remember enough of my linear algebra to prove such a theory.

I hope I was able to explain the problem well enough. I hope someone is able to help me with this problem.

Thank you,

Vincent

 

[marcusl]

The eigenvalues of a sum of matrices C=A+B equal the sum of their eigenvalues, that is, c_n = a_n+b_n, only in the most special of cases. A and B diagonal is one such case. In general your proposed approach is invalid.

 

[vkillion]

Thank you for your response.

I knew it wouldn't be as easy as adding them together. I wonder though if there isn't a closed-form solution to this, maybe there are some approximations we can make.

 

[HallsofIvy]

IF A and B have the same eigenvector, v, then $(A+ B)v= Av+ Bv= \lambda_A v+ \lambda_B v= (\lambda_A+ \lambda_B) v$
but that is a very special situation.

 

[blue_raver22]

Hello Vincent,

Thank you for reaching out with your linear algebra problem. I understand that you need to find the eigenvalues and eigenvectors of very large matrices, and you are currently using MATLAB's eig() function to do so. You are looking for optimizations to reduce computing time, and you have a theory about using the eigenvalues of H_constant and H_location to get the eigenvalues of H_final.

Your theory is actually correct and can be proven using some basic linear algebra principles. First, let's define what eigenvalues and eigenvectors are. Eigenvalues are the values that, when multiplied by the corresponding eigenvectors, give the same matrix back. In other words, they represent the scaling factor for a particular direction in the matrix. Eigenvectors are the vectors that, when multiplied by the corresponding eigenvalue, give the same vector back. They represent the direction in which the matrix is stretched or compressed.

Now, let's look at your matrices. H_constant is real and symmetric, which means it has real eigenvalues and orthogonal eigenvectors. This makes it easier to compute the eigenvalues and eigenvectors. H_location is diagonal, which means it has the eigenvalues on the diagonal and the eigenvectors are just the standard basis vectors. H_final, being the sum of H_constant and H_location, will have the same eigenvectors as H_constant and the eigenvalues will be the sum of the eigenvalues of H_constant and H_location.

So, your theory of using the eigenvalues of H_constant and H_location to get the eigenvalues of H_final is correct. This will significantly reduce your computation time from tens of thousands of eig() calls to just one eig() call and some simple operations. I would recommend implementing this approach and testing it to see the improvement in computing time.

I hope this helps and good luck with your simulations!

Best,