Numerical Analysis: the power method with shifts

AI Thread Summary
The discussion centers on determining the optimal shift value, beta, for the power method applied to a symmetric matrix A to achieve the fastest convergence to its largest eigenvalue, lambda_1. It is noted that large positive or negative values of beta hinder the method's effectiveness. The convergence rate improves when the ratio of the largest to the second-largest eigenvalues is significant, but this does not always correspond to the shifted eigenvalues. A proposed solution suggests that beta should be set to -(1/2)(lambda_2 + lambda_n) to optimize convergence, although the reasoning behind this conclusion requires further clarification. The thread concludes with a request for assistance in deriving this result.
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!
 
Thread 'PLL - How to find all the gains of a PI corrector and fix Ki ? MATLAB'
Hello everyone, it's not really for my homework (because I'm not at school anymore) but I'm new and I don't know how to start a new forum. I've never done automatic before and I need help. I have a 2nd ordre PLL on MATLAB, with an NCO and a PI corrector; the PLL works via the Mid/End algorithm (that's part is ok). I need to find the link between the theoretical formulas and their values in the code to find the value of the gain Ki, and possibly that of Kp if it's not correct. Unfortunately...
Thread 'Have I solved this structural engineering equation correctly?'
Hi all, I have a structural engineering book from 1979. I am trying to follow it as best as I can. I have come to a formula that calculates the rotations in radians at the rigid joint that requires an iterative procedure. This equation comes in the form of: $$ x_i = \frac {Q_ih_i + Q_{i+1}h_{i+1}}{4K} + \frac {C}{K}x_{i-1} + \frac {C}{K}x_{i+1} $$ Where: ## Q ## is the horizontal storey shear ## h ## is the storey height ## K = (6G_i + C_i + C_{i+1}) ## ## G = \frac {I_g}{h} ## ## C...
Back
Top