This report shows on real data that the direct methods such as LDL decomposition and Gaussian elimination for solving linear systems with ill-conditioned matrices provide inaccurate results due to divisions by very small numbers, which in turn results in peaking phenomena and large estimation errors. Richardson iteration provides accurate results without peaking phenomena since division by small numbers is absent in the Richardson approach. In addition, two preconditioners are considered and compared in the Richardson iteration: 1) the simplest and robust preconditioner based on the maximum row sum matrix norm and 2) the optimal one based on calculation of the eigenvalues. It is shown that the simplest preconditioner is more robust for ill-conditioned case and therefore it is recommended for many applications.
翻译:本报告在真实数据中显示,直接方法,如LDL分解法和Gaussian消除法等用于用条件差的矩阵解决线性系统的直接方法,由于极小的分化而得出不准确的结果,这反过来又导致现象的高峰和大的估计误差。Richardson的迭代法提供了准确的结果,而没有达到高峰现象,因为在Richardson方法中没有按小数进行分解。此外,在Richardson重复法中考虑和比较了两个先决条件:1)基于最大行和总总矩阵规范的最简单和可靠的先决条件;2)基于计算精精度值的最佳先决条件。这表明,最简单的先决条件对于条件差的情况更为有力,因此建议许多应用。