常微分方程的函数逼近与曲线拟合实验.docVIP

实验二 函数逼近与曲线拟合 【】 在某冶炼过程中，根据统计数据的含碳量与时间关系，试求含碳量与时间的拟合曲线。 (分) 0 5 10 15 20 25 30 35 40 45 50 55 0 1.27 2.16 2.86 3.44 3.87 4.15 4.37 4.51 4.58 4.02 4.64 【】1、用最小二乘法进行曲线拟合； 2、近似解析表达式为； 3、打印出拟合函数，并打印出与的误差，； 4、另外选取一个近似表达式，尝试拟合效果的比较； 5、* 绘制出曲线拟合图。 #include stdio.h #include conio.h #include stdlib.h #include math.h #define N 12//N个节点 #define M 2//M次拟合 #define K 2*M void zhuyuan (int k,int n,float a[M+1][M+2]) {int t,i,j; float x,y; x=fabs(a[k][k]);t=k; for (i=k+1;i=n;i++) if (fabs(a[i][k])x) {x=fabs(a[i][k]);t=i;} for (j=k;j=n+1;j++) {y=a[k][j];a[k][j]=a[t][j];a[t][j]=y;} } void xiaoyuan(int n,float a[M+1][M+2]) {int k,i,j; for(i=0;in;i++) {zhuyuan(i,n,a); for (j=i+1;j=n;j++) for (k=i+1;k=n+1;k++) a[j][k]=a[j][k]-a[j][i]*a[i][k]/a[i][i]; } } void huidai(int n,float a[M+1][M+2],float x[M+1]) {int i,j; x[n]=a[n][n+1]/a[n][n]; for (i=n-1;i=0;i--) { x[i]=a[i][n+1]; for (j=i+1;j=n;j++) x[i]=x[i]-a[i][j]*x[j]; x[i]=x[i]/a[i][i]; } } void main() {float x_y[N][2],A[N][K+1],B[N][M+1],AA[K+1],BB[M+1],a[M+1][M+2],m[M+1]; int i,j,n; 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]); } for(i=0;iN;i++) { A[i][0]=1; for(j=1;j=K;j++) A[i][j]=A[i][j-1]*x_y[i][0]; for(j=0;j=M;j++) B[i][j]=A[i][j]*x_y[i][1]; } for(j=0;j=K;j++) for(AA[j]=0,i=0;iN;i++) AA[j]+=A[i][j]; for(j=0;j=M;j++) for(BB[j]=0,i=0;iN;i++) BB[j]+=B[i][j]; for(i=0;iM+1;i++) {a[i][M+1]=BB[i]; for(j=0;j=M;j++) a[i][j]=AA[i+j];} n=M; printf(正规系数矩阵为:\n); for(i=0;i=n;i++) {for(j=0;j=n+1;j++) printf(%f ,a[i][j]); printf(\n);} xiaoyuan(n,a); huidai(n,a,m); printf(拟合曲线方程为:\ny(x)=%g,m[0]); for(i=1;i=n;i++) { printf( + %g,m[i]); for(j=0;ji;j++) { printf(*X); } } } p3=polyfit(x,y,3); y3=polyval(p3,x); e3=norm(y-y3); t=0:5:60; pt3=polyval(p3,t); plot(t,pt3); plot(t,pt3); title(3次拟合函数)