给我一个MINRES迭代法的Matlab代码

MINRES (Modified Newton-Raphson Method for Non-Hermitian Linear Systems) 是一种用于求解非对称线性系统Ax=b的迭代方法。下面是一个简单的MATLAB代码示例：

function [x, iter] = minres(A, b, tol, maxiter)
% A: 系统矩阵
% b: 方程右侧向量
% tol: 迭代终止误差阈值
% maxiter: 最大迭代次数

% 初始化
x = zeros(size(b));
r = b - A*x; % 初始残差
rprev = r;
H = eye(size(A)); % 使用单位矩阵作为初始Hessians的近似
iter = 0;

while norm(r) > tol && iter < maxiter
    % 计算下一个搜索方向 p
    [v, d] = eig(H, r'*r);
    index = find(abs(d) == min(d(abs(d) > 0)));
    if ~isempty(index)
        p = v(:, index(1)) * sign(diag(H)*d(index(1)));
    else % 如果所有特征值都是负的，则使用原点方向
        p = r / norm(r);
    end
    
    % 更新迭代向量
    alpha = r'*r / p'*A*p;
    x = x + alpha*p;
    
    % 更新残差和Hessian
    r = r - alpha*A*p;
    H = H - alpha*(p'*A)*pinv(p'*A)*p;
    
    % 检查收敛
    rprev = r;
    iter = iter + 1;
end
end

这个函数会返回迭代得到的解`x`以及迭代次数`iter`。请注意，在实际应用中，MINRES可能需要处理复数情况，并且MATLAB的`pinv`函数在这里用于计算伪逆。