数值分析 三四章数分析 三四章.doc

数值分析第二次上机作业 第三章 函数逼近与快速傅里叶变换 三次曲线拟合的程序及计算结果: #include stdio.h #include conio.h #include math.h #include process.h #define N 5//N个点 #define T 3 //T次拟合 #define W 1//权函数 #define PRECISION 0.00001 float pow_n(float a,int n) { int i; if(n==0) return(1); float res=a; for(i=1;in;i++) {res*=a;} return(res); } void mutiple(float a[][N],float b[][T+1],float c[][T+1]) { float res=0; int i,j,k; for(i=0;iT+1;i++) for(j=0;jT+1;j++) {res=0; for(k=0;kN;k++) {res+=a[i][k]*b[k][j]; c[i][j]=res; } } } void matrix_trans(float a[][T+1],float b[][N]) { int i,j; for(i=0;iN;i++) {for(j=0;jT+1;j++) {b[j][i]=a[i][j]; } } } void init(float x_y[][2],int n) { int i; printf(请输入%d个已知点：\n,N); for(i=0;in;i++) { printf((x%d y%d):,i,i); scanf(%f %f,x_y[i][0],x_y[i][1]); } } void get_A(float matrix_A[][T+1],float x_y[][2],int n) { int i,j; for(i=0;iN;i++) { for(j=0;jT+1;j++) { matrix_A[i][j]=W*pow_n(x_y[i][0],j); } } } void print_array(float array[][T+1],int n) { int i,j; for(i=0;in;i++) { for(j=0;jT+1;j++) { printf(%-g,array[i][j]); } printf(\n); } } void convert(float argu[][T+2],int n) { int i,j,k,p,t; float rate,temp; for(i=1;in;i++) { for(j=i;jn;j++) { if(argu[i-1][i-1]==0) { for(p=i;pn;p++) { if(argu[p][i-1]!=0) break; } if(p==n) { printf(方程组无解!\n); exit(0); } for(t=0;tn+1;t++) { temp=argu[i-1][t]; argu[i-1][t]=argu[p][t]; argu[p][t]=temp; } } rate=argu[j][i-1]/argu[i-1][i-1]; for(k=i-1;kn+1;k++) { argu[j][k]-=argu[i-1][k]*rate; if(fabs(argu[j][k])=PRECISION) argu[j][k]=0; } } } } void compute(float argu[][T+2],int n,float root[]) { int i,j; float temp; for(i=n-1;i=0;i--) { temp=argu[i][n]; for(j=n-1;ji;j--) { temp-=argu[i][j]*root[j]; } root[i]=temp/argu[i][i]; } } void get_y(float trans_A[][N],float x_y[][2],float y[],int n) { int i,j; float temp; for(i=0;in;i++) { temp=0; for(j=0;jN;j++) { temp+=trans_A[i][j]*x_y[j][1]; } y[i]=temp; } } void cons_formula(float coef_A[][T+1],float y[],float coef_form[][T+2]) { int i,j; for(i=0;iT+1;i++) { for(j=0;jT+2;j++) { if(j==T+1) coef_form[i][j]=y[i]; el