Matlab的命令里只有QR分解,如何求QL分解和LQ,RQ分解?
LQ分解:
[Q1,R1]= qr (H');
H=(Q1*R1)'=R1' *Q1'=L*Q;
所以下三角阵L=R1',正交阵Q=Q1'
QL分解:
[Q1,R1]= qr ( inv(H' ) );
H'=(Q1*R1)-1= R1-1* Q1-1;
H= (R1-1* Q1-1) ' =( Q1-1)'
* (R1-1)'= Q1* (R1-1)';
所以正交阵Q= Q1,下三角阵L= (R1-1)'.
RQ分解:
[Q1, R1]=qr(inv(H));
H-1= Q1* R1;
H=( R1)-1*( Q1)-1= R1-1*
Q1';
所以上三角阵R= R1-1,正交阵Q= Q1'.
代码:matlab求QR分解,QL分解,LQ分解和RQ分解。
最新推荐文章于 2023-12-04 21:24:27 发布