MHB Approximation of eigenvalue with inverse iteration method

AI Thread Summary
The discussion focuses on using the inverse iteration method to approximate the eigenvalue of a given matrix, starting with an initial approximation of 1.2 and a vector of (1, 1, 1). Participants clarify that the power iteration method is not suitable for this case, as it only finds the largest eigenvalue. Instead, they suggest using the inverse iteration method, which requires an initial eigenvalue approximation, and the Rayleigh quotient iteration method, which refines both the eigenvalue and eigenvector. It is recommended to apply the Rayleigh quotient at each step for better accuracy. Overall, the consensus leans towards using the Rayleigh quotient iteration for improved results.
mathmari
Gold Member
MHB
Messages
4,984
Reaction score
7
Hey! :giggle:

We have the matrix $\begin{pmatrix}2 & 1/2 & 1 \\ 1/2 & 3/2 & 1/2 \\ 1 & 1/2 & 2\end{pmatrix}$.
We take as initial approximation of $\lambda_2$ the $1.2$. We want to calculate this value approximately using the inverse iteration (2 steps) using as starting vector $x^{(0)}=\begin{pmatrix}1 \\ 1\\ 1\end{pmatrix}$.

At the inverse iteration method do we have to use at each step the Rayleigh-Quotient or only at the beginning and then just the power iteration ?
I think to get a better approximation that we have to use the Rayleigh-Quotient at each step. Is that correct?

:unsure:
 
Mathematics news on Phys.org
Hey mathmari!

It seems to me that we cannot use the power iteration method, since it can only find the eigenvalue with the greatest magnitude.
Instead we can use the inverse iteration method, which finds an eigenvector given an approximation of the eigenvalue, which does not improve the eigenvalue.
Or we can use the Rayleigh quotient iteration method, which is an extension of the inverse iteration method. It improves the eigenvalue as well as the eigenvector.

Since the initial eigenvalue approximation is still a bit off, it seems to me that it's best to use the Rayleigh quotient iteration for every iteration.
Or else we should try some variants and see what works best. :unsure:
 
Klaas van Aarsen said:
It seems to me that we cannot use the power iteration method, since it can only find the eigenvalue with the greatest magnitude.
Instead we can use the inverse iteration method, which finds an eigenvector given an approximation of the eigenvalue, which does not improve the eigenvalue.
Or we can use the Rayleigh quotient iteration method, which is an extension of the inverse iteration method. It improves the eigenvalue as well as the eigenvector.

Since the initial eigenvalue approximation is still a bit off, it seems to me that it's best to use the Rayleigh quotient iteration for every iteration.
Or else we should try some variants and see what works best. :unsure:

Ok! Thank you for your answer! 🤩
 
Suppose ,instead of the usual x,y coordinate system with an I basis vector along the x -axis and a corresponding j basis vector along the y-axis we instead have a different pair of basis vectors ,call them e and f along their respective axes. I have seen that this is an important subject in maths My question is what physical applications does such a model apply to? I am asking here because I have devoted quite a lot of time in the past to understanding convectors and the dual...
Insights auto threads is broken atm, so I'm manually creating these for new Insight articles. In Dirac’s Principles of Quantum Mechanics published in 1930 he introduced a “convenient notation” he referred to as a “delta function” which he treated as a continuum analog to the discrete Kronecker delta. The Kronecker delta is simply the indexed components of the identity operator in matrix algebra Source: https://www.physicsforums.com/insights/what-exactly-is-diracs-delta-function/ by...

Similar threads

Replies
4
Views
1K
Replies
21
Views
4K
Replies
2
Views
2K
Replies
9
Views
3K
Replies
5
Views
2K
Replies
1
Views
10K
Replies
10
Views
2K
Back
Top