曲线拟合C语言程序

1、#in clude<stdio.h>#in clude<math.h>void n ihe();void gs();void mai n()int i,j,m, n;float o50;float x50 , y50,a5050;printf(”输入数据节点数n =”,n);scan f("%d",&n);for( i=1;i<=n ;i+)prin tf(" i =%dn",i);printf("各节点的数据xi n");for(i=1;i<=n ;i+)prin tf("x%d

2、 = ",i);scan f("%f", &xi);printf("各节点的数据yi n”);for(i=1;i<=n ;i+)prin tf("y%d = ",i); scan f("%f", &yi);prin tf("n");printf("拟合的多项式次数m =", m);sca nf("%d",&m);void nihe(float x50, float y50, int m ,int n)int i,j,k=0,c

3、=1,w=1;float f,a5050, o50;dof=0;for(i=1;i<=n ;i+)k);f=f+pow( xi , k)*pow( xiacc=f ;ac+1c-1=f;ac-1c+1=f;C+；k+;while(k<=m);k=1;c=1;dof=0;for(i=1;i<=n ;i+) f=f+pow( xi , k);ac+1c=f;acc+1=f;c+;k+;k+;while(k<=m+1);k=0;c=1;dof=0;for(i=1;i<=n ;i+)f=f+yi*pow( xi , k);k+;acm+2=f;c+;while(k<

4、=m);for(i=1;i<=m+1;i+)for(j=1;j<=m+2;j+)prin tf("%ft",aij);prin tf("n");printf(”n");void gs(float a5050, float o50, int m)int i,j,c=1,z=1,t=1,w=1, n;float l,max,p,q;n=m+1;domax=azt;for(i=z+1,j=t;i<=n ;i+)if(fabs (aij)>fabs(max)max=aij;w=i;elsew=z;if(c<=n-1)for

5、(i=z,j=t;j<=n+1;j+)q=aij;aij=awj;awj=q;for(i=1;i<=n ;i+)for(j=1;j<=n+1;j+)prin tf("%ft",aij);prin tf("n ”);printf(”n");for(i=z, w=t;i<=n; i+)l=ai+1w/aww;for(j=1;j<=n+1;j+)ai+1j=ai+1j-l*awj;z+;t+;C+;while(c<=n);if(a n n=0)printf(” 矩阵奇异 n”); return;elsei=n;oi=(ai n+1/aii);for(i=n-1;i>=1;i-)p=0;for(j=i;j< n;j+)p=p+ai