Dominant eigenvalue using the inverse power method

22 次查看(过去 30 天)
I wrote the code for the inverse power method and the latter is giving me the least dominant eigenvalue. Here is my code
function [v,lamda] = IPM(B,tol)
tic;
A=inv(B);
n=size(A,1);
v=rand(n,1);
v=v/norm(v);
res=1;
while (res > tol)
W= A*v;
lamda=max(abs(W));
v= W/lamda;
res=norm(A*v-lamda*v);
toc
end
I want to invert it to get the dominant eigenvalue. Please help.
  1 个评论
Geoff Hayes
Geoff Hayes 2014-5-3
It is my understanding that the Power Method returns the dominant eigenvalue, and that the Inverse Power Method can be used to find the smallest (least dominant) eigenvalue OR it can be used to approximate an eigenvalue close to some number (which implies the least dominant eigenvalue if the chosen number is zero…which is the case in the above code).

请先登录,再进行评论。

回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Linear Algebra 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by