知识库 > 数学方法 > 数值计算 > 矩阵的分解

QR分解    

  • 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);