Numerical Analysis: the power method with shifts

Click For Summary
SUMMARY

The discussion focuses on determining the optimal shift value, beta, for the power method applied to a symmetric matrix A with distinct eigenvalues. The optimal beta is identified as beta = -(1/2)(lambda_2 + lambda_n), which enhances convergence to the largest eigenvalue, lambda_1. The analysis emphasizes that the convergence rate is influenced by the ratio of the largest to the second largest eigenvalues, particularly when they are centered around zero. The discussion concludes with a successful resolution of the problem after exploring various cases of eigenvalue distributions.

PREREQUISITES
  • Understanding of symmetric matrices and eigenvalues
  • Familiarity with the power method for eigenvalue computation
  • Knowledge of convergence rates in numerical analysis
  • Basic concepts of matrix shifts and their effects on eigenvalues
NEXT STEPS
  • Study the power method for eigenvalue problems in detail
  • Explore the impact of matrix shifts on convergence rates
  • Learn about the Rayleigh quotient and its application in eigenvalue problems
  • Investigate numerical stability in iterative methods for eigenvalue computation
USEFUL FOR

Students and professionals in numerical analysis, particularly those focusing on eigenvalue computations and iterative methods. This discussion is beneficial for anyone looking to optimize the power method for symmetric matrices.

sarahr
Messages
13
Reaction score
0

Homework Statement



Consider a symmetric matrix, A, n x n with distinct eigenvalues lambda_1 > lambda_2 > ... > lambda_n (note: i didnt miss anything here typing this, there are no absolute values here). What value of the shift beta will give fastest convergence to lamba_1 and its eigenvalue when the power method is applied to A + betaI ?

Homework Equations





The Attempt at a Solution



First, I know that when beta is very large positive or
negative: the power method works badly or not at all for computing
lambda_1.

Second, the rate of convergence is best when the ratio
between the largest and second largest eigenvalues (in magnitude) is
large. But.. these are not necessarily the shifted versons of lambda_1
and lambda_2 (lambda_1+ beta, lambda_2 + beta). For some shifts,
the shifted version lambda_n at the other end might be the largest or
second largest in magnitude.

I've been trying to look at all the different cases (like the largest eig & second largest eig being both positive, both negative, both centered around zero, etc). I was thinking that the the ratio was largest when the largest eig and the second largest eig are equally centered around zero.. but not completely sure how i should state beta..

any further ideas?
 
Physics news on Phys.org
nevermind! i figured it out! :)
 
Answer?

Does anyone know how to solve this question?

I think the answer is beta = -(1/2)(lambda_2 + lambda_n), however I have no idea how to reach this. Any help appreciated!
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
4
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
1
Views
4K
  • · Replies 6 ·
Replies
6
Views
2K