数值分析实验.docVIP

/*拉格朗日插值法*/ #include stdio.h main() { static float Lx[10],Ly[10]; int n,i,j; float x,y,p; y=0; printf(enter n=); scanf(%d,n); // 输入数组xi,yi的个数 printf(enter xi\n); for(i=1;i=n;i++) scanf(%f,Lx[i]); printf(enter yi\n); for(i=1;i=n;i++) scanf(%f,Ly[i]); printf(enter x=); scanf(%f,x); for(i=1;i=n;i++) { p=1 ; for(j=1;j=n;j++) { if(i!=j) p = p*(x-Lx[j])/(Lx[i]-Lx[j]) ; // 计算插值基函数p } y=p*Ly[i]+y; } printf(y=%f\n,y); } /*最小二乘法*/ #includeiostream using namespace std; void main() { int m,n,i,j,k,l; float s,x[100],y[100],a[100][100],b[100],c[100][100],ct[100][100]; cout请输入m的值：; cinm; cout请输入n的值：; cinn; cout请输入x的值：endl; for(i=1;i=m;i++) cinx[i]; cout请输入y的值：endl; for(i=1;i=m;i++) ciny[i]; cout数组c的值：endl; for(i=1;i=m;i++) { c[i][1]=1; for(j=2;j=n+1;j++) c[i][j]=x[i]*c[i][j-1]; } for(i=1;i=m;i++) { for(j=1;j=n+1;j++) coutc[i][j] ; coutendl; } cout数组ct的值：endl; for(i=1;i=m;i++) for(j=1;j=n+1;j++) ct[j][i]=c[i][j]; for(j=1;j=n+1;j++) { for(i=1;i=m;i++) coutct[j][i] ; coutendl; } for(k=1;k=n+1;k++) for(j=1;j=n+1;j++) { s=0; for(i=1;i=m;i++) s=s+ct[k][i]*c[i][j]; a[k][j]=s; } cout数组a的值：endl; for(i=1;i=n+1;i++) { for(j=1;j=n+1;j++) couta[i][j] ; coutendl; } cout数组b的值：endl; for(k=1;k=n+1;k++) { s=0; for(i=1;i=m;i++) s=s+ct[k][i]*y[i]; b[k]=s; coutb[k]endl; } cout所求数组的值为:; k=1; while(k!=n+1) { if(a[k][k]!=0) { for(i=k+1;i=n+1;i++) { a[i][k]=a[i][k]/a[k][k]; b[i]=b[i]-a[i][k]*b[k]; for(j=k+1;j=n+1;j++) a[i][j]=a[i][j]-a[i][k]*a[k][j]; } } k=k+1; } for(k=n+1;k=1;k--) { s=0; for(l=k+1;l=n+1;l++) s=s+a[k][l]*b[l]; b[k]=(b[k]-s)/a[k][k]