A为非奇异矩阵,则存在正交阵Q和上三角阵R,使得A有:A=QR,且当R的对角元均为正时,分解是唯一的.
定义:mt_decompose_qr(M:Array,Q:Array,R:Array)
说明:对矩阵进行QR分解,M=Q*R
定义:mt_decompose_qr (m:array,q:array,r:array,E:array);
说明:对矩阵进行QRe分解,M:*e=Q:*R.该分解使得R矩阵的对角元按绝对值从大倒小排列
参数:
M:待分解矩阵
Q:变参返回正交矩阵Q
R:变参返回上三角矩阵R
e:变参返回列置换矩阵e
例题
H := corr(rand(100,4));
mt_decompose_qr(H,Q,R); //QR分解
return array(Q,R);