我有一个类型的线性系统
我想最小化。我知道有不同的方法可以直接求解具有不同数值稳定性的系统(正规方程 + Cholesky、QR 分解、SVD 分解)。
我想问一下如何估计数值误差对每个参数的估计的影响。
可以仅根据条件编号完成吗?
是否应该使用蒙特卡罗模拟来完成?
应该如何设计这样的模拟?
我有一个类型的线性系统
我想最小化。我知道有不同的方法可以直接求解具有不同数值稳定性的系统(正规方程 + Cholesky、QR 分解、SVD 分解)。
我想问一下如何估计数值误差对每个参数的估计的影响。
可以仅根据条件编号完成吗?
是否应该使用蒙特卡罗模拟来完成?
应该如何设计这样的模拟?
我认为大多数求解器都会给你一个残差。如果你翻转你的方程来阅读:
求解器通常会为您提供残差,例如:.
现在,如果您想比较不同的求解器,您可以将其用作误差的度量。根据经验:您不会找到适合每个系统 A 的求解器。通常还有速度/稳定性权衡等。