数值分析 曲线拟合

《数值分析 曲线拟合》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第三章实验作业由物理实验得到下列一组数据：xi11.522.533.544.555.56yi33.4079.50122.65159.05189.15214.15238.65252.50267.55280.50296.65xi6.577.58yi301.65310.40318.15325.15用抛物线做曲线拟合。1.N-S图;C语言程序：#include"stdio.h"#include"conio.h"#include"stdlib.h"#include"math.h"#defineN15//N个节点#defineM2//M次拟合#def

2、ineK2*Mvoidzhuyuan(intk,intn,floata[M+1][M+2]){intt,i,j; floatx,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;}}voidxiaoyuan(intn,floata[M+1][M+2]){intk,i,j; for(i=0;i

3、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]; }}voidhuidai(intn,floata[M+1][M+2],floatx[M+1]){inti,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

4、]/a[i][i]; }}voidmain(){floatx_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];inti,j,n;printf("请输入%d个已知点：",N);for(i=0;i

5、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;i

6、拟合曲线方程为:y(x)=%g",m[0]);for(i=1;i<=n;i++){printf("+%g",m[i]);for(j=0;j