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
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 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!