计算数值方法实验报告
double s; int i,j,k,n,index; cout<<\请输入已知点的个数n=\cin>>n; cout< 计算数值方法实验报告 b[i]=sumY[i]; } k=1; //用高斯消元法解方程组 do{ for(j=k+1;j<=index+1;j++) l[j][k]=a[j][k]/a[k][k]; for(i=k+1;i<=index+1;i++){ for(j=k+1;j<=index+1;j++) a[i][j]=a[i][j]-l[i][k]*a[k][j]; b[i]=b[i]-l[i][k]*b[k]; } if(k==index+1) break; k++; }while(1); x[index+1]=b[index+1]/a[index+1][index+1]; for(i=index;i>=1;i--) { s=0; for(j=i+1;j<=index+1;j++) s=s+a[i][j]*x[j]; x[i]=(b[i]-s)/a[i][i]; } cout<<\拟合系数为:\输出拟合系数 for(i=1;i<=index+1;i++) cout< 22 计算数值方法实验报告 五、结果讨论和分析 分析： 数据拟合的具体作法是：对给定的数据（xi ,yi）（i=0,1,…,m），在取定的函数类中，求p(x)属于此函数类，使误差ri=p(xi)- yi (i=0,1,…,m)的平方和最小，即 ∑ri=∑（∑p(xi)-yi）=min 22从几何意义上讲，就是寻求与给定点（xi ,yi）（i=0,1,…,m）的距离平方和为最小的曲线y=p(x)。 实验地点 23 综合楼六层606室 指导教师 王峥