※ 引述《j0958322080 (Tidus)》之铭言：
: ※ [本文转录自 C_and_CPP 看板 #1QFvpcqU ]
: 作者: j0958322080 (Tidus) 看板: C_and_CPP
: 标题: [问题] 解最小平方法的问题 Ax~b
: 时间: Sun Dec 24 20:26:10 2017
: ※ 引述《j0958322080 (Tidus)》之铭言：
: 开发平台(Platform): (Ex: Win10, Linux, ...)
: win10
: 编译器(Ex: GCC, clang, VC++...)+目标环境(跟开发平台不同的话需列出)
: GCC
: 额外使用到的函数库(Library Used): (Ex: OpenGL, ...)
: NO
: 问题(Question)：
: 我想要去FIT一条四次方的曲线，其中 x 的值为50000左右，
: 依照理论我会用到x^4，这样整个矩阵A*A^T的最大值与最小值会差到40次方，
: 我自己写了一个程式用 LU 分解去计算反矩阵，求得的反矩阵跟 EXCEL 的结果完全一样，
: 可是我发现那两个矩阵(A*A^T)和(A*A^T)^-1在 EXCEL 里面乘起来不是单位矩阵，
: 而且有些非对角线元素甚至达到10^8，这样的结果不知道是否会与我想要的解差很多??
: 因为目前只有想到用反矩阵解，不知道有没有什么比较好的算法可以解的比较精确??
: 喂入的资料(Input)：
: [ 1 x x^2 x^3 x^4] = A
这是 Pseudo-inverse matrix 的 Condition number 问题
在数值分析的书应该会提到.
去找 Golomb Matrix Computation 应该有讲到