对称正定矩阵分解，计算结果错误，实在找不到原因了

chiangtp · 发表于 2015-3-29 00:52:23

1. Single precision (~7 precision digits) 的計算環境下, 反算驗證, 兩個答案都滿足

a.f90 (1.12 KB, 下载次数: 6)

2. 矩陣A須以 double precision (15-16 precision digits) 表達, 方得見兩個相異的答案

aa.f90 (2.61 KB, 下载次数: 6)

3. 系統(矩陣A)有點 ill-conditioned, Condition number (based on L2 norm) C~7000

ill-conditioned.pdf (120.86 KB, 下载次数: 8)

chiangtp · 发表于 2015-4-1 14:19:37

請心直口快的pasuka參考/評價: 數值方法與程式(林聰悟-林佳慧)

[线性代数] 对称正定矩阵分解，计算结果错误，实在找不到原因了