用最小二乘法做直线拟合真的是非常简单,我随便写了一下并经过测试,需要使用的朋友随便改几个地方(例如变量名等)就可以用在自己的应用程序中了。在这里给出拟合直线Y=a0+a1X 以及Y=aX的例子。
首先是关于最小二乘法的数学基础,就不多说了,无非就是求偏导之类的。
在VC中的具体代码如下(先是拟合直线Y=a0+a1X的):
函数声明:
void linearReg_Y_a0_a1X(double &a_0, double &a_1, const CArraydouble,double & X_Con_Lin_Reg, const CArraydouble,double & Y_Volts_Lin_Reg);
函数定义:
void CYourDlg::linearReg_Y_a0_a1X(double &a0, double &a1, const CArraydouble,double & X_Con_Lin_Reg, const CArraydouble,double & Y_Volts_Lin_Reg)
{
double sumXiYi(0.0), sumXi(0.0), sumYi(0.0)...[ 查看全文 ]