最小二乘法实验报告

计算方法与实习实验报告2.已知实验数据如下：Xi1.02.53.54.0Yi3.81.5026.033.0试用形如y=a+b𝑥2的抛物线进行最小二乘拟合。程序如下：#includeiostream.hdoublekk(){doublek=0;for(inti=0;i4;i++){k+=1*1;}returnk;}doublekj(doublex[4]){doublek=0;for(inti=0;i4;i++){k+=x[i]*x[i];}returnk;}doublejj(doublex[4]){doublek=0;for(inti=0;i4;i++){k+=x[i]*x[i]*x[i]*x[i];}returnk;}doubleyk(doubley[4]){doublek=0;for(inti=0;i4;i++){k+=y[i];}returnk;}doubleyj(doublex[4],doubley[4]){doublek=0;for(inti=0;i4;i++){k+=y[i]*x[i]*x[i];}returnk;}voidmain(void){doubleX[4]={1.0,2.5,3.5,4.0};doubleY[4]={3.8,1.50,26.0,33.0};doubleA[2][3];doublea=0,b=0;A[0][0]=kk();A[0][1]=kj(X);A[1][0]=kj(X);A[1][1]=jj(X);A[0][2]=yk(Y);A[1][2]=yj(X,Y);doublek=A[1][0]/A[0][0];A[1][0]=0;A[1][1]=A[1][1]-k*A[0][1];A[1][2]=A[1][2]-k*A[0][2];b=A[1][2]/A[1][1];a=(A[0][2]-A[0][1]*b)/A[0][0];cout二次曲线拟合为:y=a+bx*xendl;}拟合结果如结果与答案几乎没有误差，可以认为算法正确，实验基本成功。心得：经过本次试验，认识了最先二乘法在曲线拟合时的优点，这激发了我认真学习最小二乘法拟合曲线的兴趣，一定要好好掌握并在生活中加以应用。