有解析解，你干嘛非整个数值解？
【 在 vociferous 的大作中提到: 】
: a*x0+b*x1+c*x2+d = y      线性方程拟合a,b,c,d 有高效的C源代码方案吗？
--
FROM 117.136.81.*

【 在 qyli 的大作中提到: 】
: 有解析解，你干嘛非整个数值解？

有误差存在的情况下，哪有解析解？
--
FROM 117.152.203.*

最小二乘相当于超定方程求解吧，A*B=Y，A的每行都是一组[x0，x1，x2，1]，B是[a,b,c,d]的转置，Y是A中各行对应的y组成的列向量。  B=inv(A^t *A) * (A^t * Y），A^t是A的转置。全主元高斯消元的求逆c代码好像很容易搜到。
【 在 snoopyzhao 的大作中提到: 】
:
: 有误差存在的情况下，哪有解析解？
--
FROM 223.104.39.*

是的，高斯消元求逆是可以的。
【 在 toutouqi 的大作中提到: 】
: 最小二乘相当于超定方程求解吧，A*B=Y，A的每行都是一组[x0，x1，x2，1]，B是[a,b,c,d]的转置，Y是A中各行对应的y组成的列向量。  B=inv(A^t *A) * (A^t * Y），A^t是A的转置。全主元高斯消元的求逆c代码好像很容易搜到。
--
FROM 50.83.12.*

线性方程就没有不是解析解的。只是变量很多或者取值病态会有数值误差放大不够稳定的问题。
【 在 snoopyzhao 的大作中提到: 】
: 【 在 toutouqi 的大作中提到: 】
: : 最小二乘相当于超定方程求解吧，A*B=Y，A的每行都是一组[x0，x1，x2，1]，B是[a,b,c,d]的转置，Y是A中各行对应的y组成的列向量。  B=inv(A^t *A) * (A^t * Y），A^t是A的转置。全主元高斯消元的求逆c代码好像很容易搜到。
:
: ...................
--
FROM 61.148.245.*