@dbb627 看到过你在知道里的回答,那个设了权限不能问,在这儿请教一下~
那个问题
function l = rqrtz(A,M)
%QR算法求矩阵全部特征值
%已知矩阵:A
%迭代步数:M
%求得的矩阵特征值:l
A = hess(A);
for i=1:M
N = size(A);
n = N(1,1);
u = A(n,n);
[q,r]=qr(A-u*eye(n,n));
A = r*q+u*eye(n,n);
l = diag(A);
end
>> rqrtz(A,50)
这个循环50次是一般就都够了,还是要根据我自己的数据在验证一下啊
谢谢了哈
那个问题
function l = rqrtz(A,M)
%QR算法求矩阵全部特征值
%已知矩阵:A
%迭代步数:M
%求得的矩阵特征值:l
A = hess(A);
for i=1:M
N = size(A);
n = N(1,1);
u = A(n,n);
[q,r]=qr(A-u*eye(n,n));
A = r*q+u*eye(n,n);
l = diag(A);
end
>> rqrtz(A,50)
这个循环50次是一般就都够了,还是要根据我自己的数据在验证一下啊
谢谢了哈