拉格朗日差值多项式(C语言实例).pdfVIP

(1)拉格朗日差值多项式; #includestdio.h #includemath.h #defineMAX 100 voidmain() { int i,j,k,m,n,N,mi; floattmp,mx; floatX[MAX][MAX],Y[MAX],x[MAX],y[MAX],a[MAX]; printf(\n 输入拟合多项式的次数:\n); scanf(%d,m); printf(\n 输入给定点的个数n 及坐标(x,y):\n); scanf(%d,N); printf(\n); for(i 0;iN;i++) scanf(%f,%f,x[i],y[i]); for(i 0;i m;i++) { for(j i;j m;j++) { tmp 0; for(k 0;kN;k++) tmp tmp+pow(x[k],(i+j)); X[i][j] tmp; X[j][i] X[i][j]; } } for(i 0;i m;i++) { tmp 0; for(k 0;kN;k++) tmp tmp+y[k]*pow(x[k],i); Y[i] tmp; } for(j 0;jm;j++) { for(i j+1,mi j,mx fabs(X[j][j]);i m;i++) if(fabs(X[i][j])mx) { mi i; mx fabs(X[i][j]); } if(jmi) { tmp Y[j]; Y[j] Y[mi]; Y[mi] tmp; for(k j;k m;k++) { tmp X[j][k]; X[j][k] X[mi][k]; X[mi][k] tmp; } } for(i j+1;i m;i++) { tmp -X[i][j]/X[j][j]; Y[i]+ Y[j]*tmp; for(k j;k m;k++) X[i][k]+ X[j][k]*tmp; } } a[m] Y[m]/X[m][m]; for(i m-1;i 0;i--) { a[i] Y[i]; for(j i+1;j m;j++) a[i]- X[i][j]*a[j]; a[i]/ X[i][i]; } printf(\n 所求的二次多项式为:\n); printf(P(x) %f,a[0]); for(i 1;i m;i++) printf(+(%f)*x^%d,a[i],i); } 运行结果： 输入拟合多项式的次数: 5 输入给定点的个数n 及坐标(x,y): 3 1,2 5,3 4,2 所求的二次多项式为: P(x) 1.980417+(0.282759)*x^1+(-0.299937)*x^2+(0.022071)*x^3+(0.016624)*x^4+(-0.0 01934)*x^5Pressanykey tocontinue 总结：拉格朗日计算公式中，只需要知道各个点即可