Power Method Eigenvectors Code
显示 更早的评论
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
2021-1-28
What do you mean by "look smoother"?
Benjamin Boettcher
2021-1-28
编辑：Benjamin Boettcher
2021-1-28
回答（1 个）
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;
类别
在 帮助中心 和 File Exchange 中查找有关 Mathematics 的更多信息
产品
2021-1-28
2021-1-28
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
