matlab的qr算法代码

以下是MATLAB中使用QR算法求解特征值和特征向量的代码：

% 定义一个矩阵A
A = [1, 2, 3; 4, 5, 6; 7, 8, 9];

% 使用QR算法求解特征值和特征向量
[Q, R] = qr(A);
for i = 1:100
    [Q, R] = qr(A);
    A = R * Q;
end

% 输出特征值和特征向量
eig(A)

其中，qr函数用于计算QR分解，eig函数用于计算矩阵的特征值和特征向量。在QR算法中，我们通过多次QR分解来逼近矩阵的特征值和特征向量。在上述代码中，我们通过100次QR分解来逼近矩阵的特征值和特征向量，最终输出特征值。