matlab求QR分解,QL分解,LQ分解和RQ分解
Matlab的命令里只有QR分解,如何求QL分解和LQ,RQ分解?
LQ分解:
[Q 1,R 1]= qr (H');
H=(Q 1*R 1)'=R 1' *Q 1'=L*Q;
所以下三角阵L=R 1',正交阵Q=Q 1'
QL分解:
[Q 1,R 1]= qr ( inv(H' ) );
H'=(Q 1*R 1) -1= R 1 -1* Q 1 -1;
H= (R 1 -1* Q 1 -1) ' =( Q 1 -1)' * (R 1 -1)'= Q 1* (R 1 -1)';
所以正交阵Q= Q 1,下三角阵L= (R 1 -1)'.
RQ分解:
[Q 1, R 1]=qr(inv(H));
H -1= Q 1* R 1;
H=( R 1) -1*( Q 1) -1= R 1 -1* Q 1';
所以上三角阵R= R 1 -1,正交阵Q= Q 1'.
转载:
快去成为你想要的样子!