【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

>>阅读原文<<