这个问题不是很深。假设我有一个小的矩形矩阵, 行数和列数之间-, 分别。
给定一个右手边, 我想解决最小二乘问题, 即我想找到具有范数最小残差和在这些最小化器中具有最小规范。
由于系统很小,而且有一个很好的广义条件数,我考虑求解正规方程相反,比如说,通过 Cholesky 分解。
在这个阶段,效率不如获得准确的解决方案重要。根据您的经验,您是否会争辩说这是一个不错的选择,而不是计算 SVD?
这个问题不是很深。假设我有一个小的矩形矩阵, 行数和列数之间-, 分别。
给定一个右手边, 我想解决最小二乘问题, 即我想找到具有范数最小残差和在这些最小化器中具有最小规范。
由于系统很小,而且有一个很好的广义条件数,我考虑求解正规方程相反,比如说,通过 Cholesky 分解。
在这个阶段,效率不如获得准确的解决方案重要。根据您的经验,您是否会争辩说这是一个不错的选择,而不是计算 SVD?
如果需要,您可以使用“半正规方程”;这里的关键是QR分解中的因子与 Cholesky 三角形相同(直至输入符号的变化), 自从
.
如果您不能或不愿意存储正交矩阵,无论是显式存储还是将其分解为 Householder 向量,这将特别方便。
简而言之,你可以解方程以获得您的最小二乘解决方案。
尽管如此,没有什么比奇异值分解所提供的诊断能力更好的了。如果你负担得起,你应该使用它,除非你提前知道你的最小二乘问题总是条件良好的。
当然,去反转. 该矩阵是对称且正定的,因此您也可以使用 CG 来做到这一点。或者,您可以使用 QR 等算法通过计算原始矩形矩阵的分解来解决原始问题。