曲线拟合C语言程序.doc

《曲线拟合C语言程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、#include#includevoidnihe();voidgs();voidmain(){inti,j,m,n;floato[50];floatx[50],y[50],a[50][50];printf("输入数据节点数n=",n);scanf("%d",&n);for(i=1;i<=n;i++){printf("i=%d",i);}printf("各节点的数据x[i]");for(i=1;i<=n;i++){printf("x[%d]=",i);scanf

2、("%f",&x[i]);}printf("各节点的数据y[i]");for(i=1;i<=n;i++){printf("y[%d]=",i);scanf("%f",&y[i]);}printf("");printf("拟合的多项式次数m=",m);scanf("%d",&m);}voidnihe(floatx[50],floaty[50],intm,intn){inti,j,k=0,c=1,w=1;floatf,a[50][50],o[50];;do{f=0;for(i=1;i<=n;i++

3、){f=f+pow(x[i],k)*pow(x[i],k);}a[c][c]=f;a[c+1][c-1]=f;a[c-1][c+1]=f;c++;k++;}while(k<=m);k=1;c=1;do{f=0;for(i=1;i<=n;i++){f=f+pow(x[i],k);}a[c+1][c]=f;a[c][c+1]=f;c++;k++;k++;}while(k<=m+1);k=0;c=1;do{f=0;for(i=1;i<=n;i++){f=f+y[i]*pow(x[i],k);}k++;a[c

4、][m+2]=f;c++;}while(k<=m);for(i=1;i<=m+1;i++){for(j=1;j<=m+2;j++){printf("%ft",a[i][j]);}printf("");}printf("------------------------------------------------------------------------");}voidgs(floata[50][50],floato[50],intm){inti,j,c=1,z=1,t=1,w=1,n

5、;floatl,max,p,q;n=m+1;do{max=a[z][t];for(i=z+1,j=t;i<=n;i++){if(fabs(a[i][j])>fabs(max)){max=a[i][j];w=i;}else{w=z;}}if(c<=n-1){for(i=z,j=t;j<=n+1;j++){q=a[i][j];a[i][j]=a[w][j];a[w][j]=q;}}for(i=1;i<=n;i++){for(j=1;j<=n+1;j++){printf("%ft",a[i][j]);}p

6、rintf("");}printf("------------------------------------------------------------------------");for(i=z,w=t;i<=n;i++){l=a[i+1][w]/a[w][w];for(j=1;j<=n+1;j++){a[i+1][j]=a[i+1][j]-l*a[w][j];}}z++;t++;c++;}while(c<=n);if(a[n][n]==0){printf("矩阵奇异");retu

7、rn;}else{i=n;o[i]=(a[i][n+1]/a[i][i]);for(i=n-1;i>=1;i--){p=0;for(j=i;j