用matlab写一个带原地移位的反幂法的代码

以下是一个用MATLAB编写的带原地移位的反幂法的代码：

function [eig_val, eig_vec] = inv_power_shift(A, mu, tol, max_iter)
% INV_POWER_SHIFT 使用带原地移位的反幂法求矩阵的最小特征值和对应的特征向量
% A: 待求特征值和特征向量的矩阵
% mu: 移位参数
% tol: 精度要求
% max_iter: 最大迭代次数
% eig_val: 最小特征值
% eig_vec: 对应的特征向量

% 初始化特征向量
n = size(A, 1);
eig_vec = ones(n, 1);

% 反幂法迭代计算
for iter = 1 : max_iter
% 计算移位后的矩阵
B = A - mu * eye(n);

% LU分解
[L, U, P] = lu(B);

% 前代和回代求解线性方程组
y = L \ (P * eig_vec);
eig_vec = U \ y;

% 正则化
eig_vec = eig_vec / norm(eig_vec);

% 计算特征值
eig_val = eig_vec' * A * eig_vec;

% 判断是否达到精度要求
if abs(eig_val - mu) < tol
break;
end

% 更新移位参数
mu = eig_val;
end