Power method to solve eigenvalues
#matlab code
A = [];
x = [];%%input A and intial x
k = 1;
[xp p] = max(abs(x));
x = x/x§;
while k <= 14
y = A*x;
u = y§;
[yp p] = max(abs(y));
if yp == 0
display(‘A has ev 0’);
break
end
err = norm(x-y/yp,inf);
x = y/yp;
if err<0.00001
display(u);
display(x);
break
end
k=k+1;
end
Power Method for dominate eigenvalue
最新推荐文章于 2022-01-28 12:13:16 发布