MATLAB用QR方法算矩阵的特征值代码

以下是MATLAB用QR方法算矩阵的特征值的代码：

function [eig_val, iter] = qr_eig(A, tol)
% A: 待求特征值的矩阵
% tol: 精度要求
% eig_val: 矩阵A的特征值
% iter: 迭代次数

n = size(A, 1);
iter = 0;
eig_val = zeros(n, 1);

while max(abs(diag(A) - eig_val)) > tol
eig_val = diag(A);
[Q, R] = qr(A - eig_val * eye(n));
A = R * Q + eig_val * eye(n);
iter = iter + 1;
end

end

其中，A为待求特征值的矩阵，tol为精度要求，eig_val为矩阵A的特征值，iter为迭代次数。