Power Method Eigenvectors Code

4 次查看（过去 30 天）
Benjamin Boettcher
Benjamin Boettcher 2021-1-28
评论： Benjamin Boettcher 2021-1-28
I am attempting to make a code for a manual computation of the Power Method.
I think I have it working properly but wanted to make sure. Any thoughts?
function [lambda,x,k] = power_method(A,x,tol,maxit);
% Initialize
n = length(A);
k = 0;
y = zeros(1,n);
tolerance = 1;
while k < maxit && (tolerance > tol)
yk = A * x; % y(k) = A* x(k)
xnext = yk / norm(yk); % X(k+1) = yk / 2norm(yk)
lambda = dot (xnext,(A * xnext)); % lambda = X(k+1) * (A * X(k+1))
k = k + 1;
tolerance = norm(yk - lambda * x);
err = tolerance;
x = xnext;
end
  2 个评论
James Tursa
James Tursa 2021-1-28
What do you mean by "look smoother"?
Benjamin Boettcher
Benjamin Boettcher 2021-1-28
编辑：Benjamin Boettcher 2021-1-28
It seems to work, I guess I just wanted it verified. and see if it could be made any better.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

James Tursa
James Tursa 2021-1-28
The only obvious thing I would note is that the following line doesn't accomplish anything because you don't use err in your code, so it can be removed:
err = tolerance;

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by