Hessenberg QR迭代计算特征值matlab代码

以下是使用Hessenberg QR迭代计算特征值的MATLAB代码：

function [eig_val, eig_vec] = hessenberg_qr_iteration(A, max_iter)
% Hessenberg QR迭代计算特征值和特征向量
n = size(A, 1);
eig_vec = eye(n);
for k = 1:max_iter
    % Hessenberg变换
    for i = n:-1:2
        x = A(i-1:i, i);
        [c, s] = givens(x(1), x(2));
        G = blkdiag(eye(i-2), c, s);
        A([i-1, i], :) = G * A([i-1, i], :);
        A(:, [i-1, i]) = A(:, [i-1, i]) * G';
        eig_vec(:, [i-1, i]) = eig_vec(:, [i-1, i]) * G';
    end
    % QR分解
    for i = 1:n-1
        [c, s] = givens(A(i, i), A(i+1, i));
        G = [c, s; -s, c];
        A(i:i+1, i:n) = G * A(i:i+1, i:n);
        A(:, i:i+1) = A(:, i:i+1) * G';
        eig_vec(:, i:i+1) = eig_vec(:, i:i+1) * G';
    end
    % 判断是否收敛
    if max(abs(diag(A, -1))) < eps
        break;
    end
end
eig_val = diag(A);
end

function [c, s] = givens(a, b)
% 计算Givens矩阵的元素
if b == 0
    c = 1;
    s = 0;
else
    if abs(b) > abs(a)
        t = -a / b;
        s = 1 / sqrt(1 + t^2);
        c = s * t;
    else
        t = -b / a;
        c = 1 / sqrt(1 + t^2);
        s = c * t;
    end
end
end

函数`hessenberg_qr_iteration`实现了Hessenberg QR迭代算法，可以计算矩阵的特征值和特征向量。函数`givens`计算Givens矩阵的元素。在迭代过程中，先将矩阵A通过Hessenberg变换转化为上Hessenberg矩阵，再通过QR分解求解特征值，直到满足收敛条件。在求解特征向量时，需要将每次变换得到的Givens矩阵乘到特征向量矩阵上，以保证特征向量的正确性。