求最小二乘法的算法

[Breeze]

就是通过二次最小平方回归所得到的二次曲线，即 y=cx2+bx+a, 用一组数据求出c, b, a,那位大下给个算法，多谢了

[Breeze]

我是要二次曲线拟合阿

[improgrammer]

是要源码还是解题思路？

[improgrammer]

思路：
1、求下列和：A=sum(x*x*x*x), B=sum(x*x*x), C=sum(x*x), D=sum(x), E=sum(y), F=sum(x*y), G=sum(x*x*y)

2、解线性方程组：FF * XX = YY;其中，FF={C,D,N;B,C,D;A,B,C}, XX={c,b,a}, YY={E,F,G}

解线性方程组的算法另外再求吧，OK？

[jixingzhong]

//最小二乘法直线拟合
BOOL CalculateLineKB(CFoldPointList *m_FoldList,double &k,double &b)
{
//最小二乘法直线拟合
//m_FoldList为关键点(x,y)的链表
//拟合直线方程(Y=kX+b)


if(m_FoldList==NULL)return FALSE;
long lCount=m_FoldList->GetCount();
if(lCount<2)return FALSE;
CFoldPoint *pFold;
double mX,mY,mXX,mXY,n;
mX=mY=mXX=mXY=0;
n=lCount;
POSITION pos=m_FoldList->GetHeadPosition();
while(pos != NULL)
{
pFold=m_FoldList->GetNext(pos);
mX+=pFold->X;
mY+=pFold->Y;
mXX+=pFold->X*pFold->X;
mXY+=pFold->X*pFold->Y;
}
if(mX*mX-mXX*n==0)return FALSE;
k=(mY*mX-mXY*n)/(mX*mX-mXX*n);
b=(mY-mX*k)/n;
return TRUE;
}

[jixingzhong]

VC 程序，楼主应该能看懂的吧，不是很深奥 ...