matlab实现QR分解求矩阵特征值

以下是使用Matlab实现QR分解求矩阵特征值的代码示例：

function [Q] = QR(A, it_max)
    A1 = A;
    for N = 1:it_max
        [Q1, R1] = QR_tao(A1);
        A2 = R1 * Q1;
        A1 = A2;
    end
    Q = A1;
end

% 测试
A = [13, -3, -1, 0; -3, 6, 0, -2; -1, 0, 10, -1; 0, -2, -1, 7];
[Q] = QR(A, 50);
eig(A)

在上述代码中，`QR`函数实现了QR分解求矩阵特征值的过程。`A`是输入的矩阵，`it_max`是迭代次数。在每次迭代中，通过调用`QR_tao`函数进行QR分解，然后更新矩阵`A1`。最后，返回更新后的矩阵`A1`作为特征值矩阵`Q`。

在测试部分，我们定义了一个矩阵`A`，然后调用`QR`函数进行QR分解求特征值，并使用`eig`函数验证结果。