下载吧在线观看:用过最小二乘法曲线拟合的高手请进，看看为何不准？？

以下这个最小二乘法的源代码是网上搜到的最多的例子，
bool CalculateCurveParameter(float *X,float *Y,long M,long N,float *A)
{
//X,Y -- X,Y两轴的坐标
//M -- 结果变量组数
//N -- 采样数目
//A -- 结果参数

register long i,j,k;
float Z,D1,D2,C,P,G,Q;
float B[50],T[50],S[50];

if(M>N)M=N;
for(i=0;i<M;i++)
A[i]=0;
Z=0;
B[0]=1;
D1=N;
P=0;
C=0;
for(i=0;i<N;i++)
{
P=P+(*X)-Z;
C=C+(*Y);
}
C=C/D1;
P=P/D1;
A[0]=C*B[0];
if(M>1)
{
T[1]=1;
T[0]=-P;
D2=0;
C=0;
G=0;
for(i=0;i<N;i++)
{
Q=(*X)-Z-P;
D2=D2+Q*Q;
C=(*Y)*Q+C;
G=((*X)-Z)*Q*Q+G;
}
C=C/D2;
P=G/D2;
Q=D2/D1;
D1=D2;
A[1]=C*T[1];
A[0]=C*T[0]+A[0];
}
for(j=2;j<M;j++)
{
S[j]=T[j-1];
S[j-1]=-P*T[j-1]+T[j-2];
if(j>=3)
{
for(k=j-2;k>=1;k--)
S[k]=-P*T[k]+T[k-1]-Q*B[k];
}
S[0]=-P*T[0]-Q*B[0];
D2=0;
C=0;
G=0;
for(i=0;i<N;i++)
{
Q=S[j];
for(k=j-1;k>=0;k--)
Q=Q*((*X)-Z)+S[k];
D2=D2+Q*Q;
C=(*Y)*Q+C;
G=((*X)-Z)*Q*Q+G;
}
C=C/D2;
P=G/D2;
Q=D2/D1;
D1=D2;
A[j]=C*S[j];
T[j]=S[j];
for(k=j-1;k>=0;k--)
{
A[k]=C*S[k]+A[k];
B[k]=T[k];
T[k]=S[k];
}
}
return true;
}
在这里我们可以知道M为阶数，A为系数。我想验证一下这个函数，输入了两个数组X[4]={1,2,3,4},Y【4】＝{1,2,3,4}，阶数为M＝1，N＝4，按照道理它应该给我们生成a【0】＝0，a【1】＝1，即y＝x，但是它居然生成了a【0】＝1，a【1】＝0，即y＝1的函数，我很茫然。求救各位高手告诉我为什么，谢谢。
其他输入不变，当我输入x＝1，2，3，4；y＝2，4，6，8，得到了a【0】＝2，a【1】＝0，基本正确，但是输入y＝2x＋1的值x＝1，2，3，4；y＝3，5，7，9时，出现了a【0】＝3，a【1】＝0的情况。
当我把幂M改成2时，输入y＝（x）2的几个值：x＝1，2，3，4；y＝1，4，9，16。得到了a【0】＝ －1.0012 a[1]= -1.0012 a[2]= 0
看来好像还是有些出入哦。