用最小二乘法做直线拟合真的是非常简单,我随便写了一下并经过测试,需要使用的朋友随便改几个地方(例如变量名等)就可以用在自己的应用程序中了。在这里给出拟合直线Y=a0+a1X 以及Y=aX的例子。
我写了一个曲线拟合的类,这个其实就是其中的一部分。目前这个类可以拟合如下几种曲线:Y=aX,Y=aX+b,Y=aX^2+bX,Y=X^a+bX,Y=alnX+b,Y=aX^b,Y=aX^2+bX+c。
首先是关于最小二乘法的数学基础,就不多说了,无非就是求偏导之类的。
在VC中的具体代码如下(先是拟合直线Y=a0+a1X的):
函数声明:
void
linearReg_Y_a0_a1X(
double
&
a_0,
double
&
a_1,
const
CArray
<
double
,
double
>
&
X_Con_Lin_Reg,
const
CArray
<
double
,
double
>
&
Y_Volts_Lin_Reg);
函数定义:
void
CYourDlg::linearReg_Y_a0_a1X(
double
&
a0,
double
&
a1,
const
CArray
<
double
,
double
>
&
X_Con_Lin_Reg,
const
CArray
<
double
,