最小二乘曲线拟合程序.doc

《最小二乘曲线拟合程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、如：根据数据表，用最小二乘法求三次多项式拟合这组数据。数据表x012345678y12.521355.543#include"stdio.h"#include"conio.h"#include"stdlib.h"#include"math.h"#defineN9#defineM2#defineK2*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[

2、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、ti,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];}}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("Pleaseinqut%d",N);for(i=0;i

4、",i,i);scanf("%f%f",&x_y[i][0],&x_y[i][1]);}for(i=0;i

5、=0;i

6、("*X");}}}输入数据表中的数据，运行结果如下：