GSL 系列 6 — 线性代数 5 — 完全正交分解

写在前面

若无特别说明，本篇代码均来自头文件 gsl_linalg.h

完全正交分解

完全正交分解可以看作是QR分解的推广，对于矩阵 $A$ ($M\times N$)，有如下分解：
$$AP=Q\left(\begin{array}{cc}{R}_{11}& 0\\ 0& 0\end{array}\right)Z^T$$
其中：
$P$：转置矩阵，$N\times N$
$Q$：正交矩阵，$M\times M$
$R_{11}$：上三角矩阵，$r\times r,r=rank(A)$
$Z$：正交矩阵， N ×