我准备在几天内进行数字考试并且有一个旧的考试要求在一个MATLAB函数中实现对称矩阵的最大和最小特征值的搜索,当然不使用 eig (A)。现在既不是功率方法也不是逆功率方法是单独实现的问题,但在我无法得到的问题中有一个规范:
您应该尽可能少地访问A的元素。
是否存在如何获取以下代码片段的技巧,使用幂迭代方法,尽可能高效地计算最小值?
电源方法的附加代码摘录:
while error > tol
i = i + 1;
y_k = z_k/norm(z_k);
z_new = A*y_k;
mu_k = transpose(y_k)*z_new;
%mu_small = transpose(y_small)*z_newsmal;
error = abs(mu_k - mu_old);
mu_old = mu_k;
z_k = z_new;
end