【matlab】 QR分解 求矩阵的特征值

"QR_H.m" function [Q,R] = QR_tao(A) %输入矩阵A %输出正交矩阵Q和上三角矩阵R [n,n]=size(A); E = eye(n); X = zeros(n,1); R = zeros(n); P1 = E; for k=1:n-1 s = -sign(A(k,k))*norm(A(k:n,k)); R(k,k) = -s; if

>>阅读原文<<