曲线拟合程序

时间:2015-01-22 16:29:30
【文件属性】:

文件名称:曲线拟合程序

文件大小:3KB

文件格式:CPP

更新时间:2015-01-22 16:29:30

曲线拟合

应用最小二乘法原理求出拟合多项式。 部分代码 int xianxingfangchengzu(double **a,int n,double *b,double *p,double dt)//用高斯列主元法来求解法方程组 { int i,j,k,l; double c,t; for(k=1;k<=n;k++) { c=0.0; for(i=k;i<=n;i++) if(fabs(a[i-1][k-1])>fabs(c)) { c=a[i-1][k-1]; l=i; }if(fabs(c)<=dt) return(0); if(l!=k) { for(j=k;j<=n;j++) { t=a[k-1][j-1]; a[k-1][j-1]=a[l-1][j-1]; a[l-1][j-1]=t; } t=b[k-1]; b[k-1]=b[l-1]; b[l-1]=t; } c=1/c; for(j=k+1;j<=n;j++) { a[k-1][j-1]=a[k-1][j-1]*c; for(i=k+1;i<=n;i++) a[i-1][j-1]-=a[i-1][k-1]*a[k-1][j-1]; } b[k-1]*=c; for(i=k+1;i<=n;i++) b[i-1]-=b[k-1]*a[i-1][k-1]; } for(i=n;i>=1;i--) for(j=i+1;j<=n;j++) b[i-1]-=b[j-1]*a[i-1][j-1]; cout.precision(12); for(i=0;i


网友评论

  • 简单了点,用处不大。
  • 可以用,但是计算误差的时候似乎有点问题
  • 算法只是实例,太简单
  • 可以使用,使用挺方便