10个重要的算C语言实现源代码.docVIP

1.拉格朗日插值多项式 ，用于离散数据的拟合 C/C++ code #include stdio.h #include conio.h #include alloc.h float lagrange(float *x,float *y,float xx,int n) /*拉格朗日插值算法*/ { int i,j; float *a,yy=0.0; /*a作为临时变量，记录拉格朗日插值多项式*/ a=(float *)malloc(n*sizeof(float)); for(i=0;i=n-1;i++) { a[i]=y[i]; for(j=0;j=n-1;j++) if(j!=i) a[i]*=(xx-x[j])/(x[i]-x[j]); yy+=a[i]; } free(a); return yy; } main() { int i,n; float x[20],y[20],xx,yy; printf(Input n:); scanf(%d,n); if(n=20) {printf(Error!The value of n must in (0,20).); getch();return 1;} if(n=0) {printf(Error! The value of n must in (0,20).); getch(); return 1;} for(i=0;i=n-1;i++) { printf(x[%d]:,i); scanf(%f,x[i]); } printf(\n); for(i=0;i=n-1;i++) { printf(y[%d]:,i);scanf(%f,y[i]);} printf(\n); printf(Input xx:); scanf(%f,xx); yy=lagrange(x,y,xx,n); printf(x=%f,y=%f\n,xx,yy); getch(); } 2.牛顿插值多项式，用于离散数据的拟合 C/C++ code #include stdio.h #include conio.h #include alloc.h void difference(float *x,float *y,int n) { float *f; int k,i; f=(float *)malloc(n*sizeof(float)); for(k=1;k=n;k++) { f[0]=y[k]; for(i=0;ik;i++) f[i+1]=(f[i]-y[i])/(x[k]-x[i]); y[k]=f[k]; } return; } main() { int i,n; float x[20],y[20],xx,yy; printf(Input n:); scanf(%d,n); if(n=20) {printf(Error! The value of n must in (0,20).); getch(); return 1;} if(n=0) {printf(Error! The value of n must in (0,20).);getch(); return 1;} for(i=0;i=n-1;i++) { printf(x[%d]:,i); scanf(%f,x[i]); } printf(\n); for(i=0;i=n-1;i++) { printf(y[%d]:,i);scanf(%f,y[i]);} printf(\n); difference(x,(float *)y,n); printf(Input xx:); scanf(%f,xx); yy=y[20]; for(i=n-1;i=0;i--) yy=yy*(xx-x[i])+y[i]; printf(NewtonInter(%f)=%f,xx,yy); getch(); } 3.高斯列主元消去法,求解其次线性方程组 C/C++ code #includestdio.h #include math.h #define N 20 int main() { int n,i,j,k; int mi,tmp,mx; float a[N][N],b[N],x[N]; printf(\nInput n:); scanf(%d,n); if(nN) { printf(The input n should in(0,N)!\n); getch(); return 1; } if(n=0