Matlab 提供了基于 Lapack 子程序 ZGEQRF/ZGEQPF 和 ZUNGQR 的 QR 分解程序。 这样的例程使对角线条目成为真实的。 对于一般对角线条件,此代码“qr_lapack”支持在 QR 分解中允许任何给定的对角复数相位,例如所有非负值,与原始矩阵 A 中的相同,或任何其他。 例如, DIAG_VEC = 诊断(A); [Q, R] = qr_lapack (A, ECON, DIAG_VEC); 此代码直接访问 Lapack 例程,而不是使用 matlab 命令 qr。 基于 Lapack 的 QR 分解后,Q 和 R 被重新计算为具有 R 对角线项的特定相位。 具体测试还有其他示例代码: * qr_diag_pos.m:QR 使对角线为正实数。 * qr_diag_phase.m:QR 产生与 A 中相同的对角相位。 * qr_diag_arg.m:另一种使
2022-04-13 15:59:34
22KB
matlab
1