H := corr(rand(100,4));
mt_decompose_qr(H,Q,R);    //QR分解
return array(Q,R);