MATLAB QR 算法求矩阵特征值代码示例

% 使用 QR 算法求解矩阵特征值

% 输入矩阵 A A = [1 2 3; 4 5 6; 7 8 9];

% 初始化 Q 和 R Q = eye(size(A)); R = A;

% 迭代计算 for i = 1:100 [Q, R] = qr(R * Q); end

% 输出特征值 eigenvalues = diag(R * Q)

% 验证结果 eig(A)