实验3__曲线拟合的最小二乘法2..docVIP

实验三 曲线拟合的最小二乘法 实验目的： 1.最小二乘法的基本原理； 2．用多项式作最小二乘曲线拟合原理的基础上，通过编程实现一组实验数据 的最小二乘拟合曲线。 2、实验要求: 认真分析题目的条件和要求，复习相关的理论知识，选择适当的解决方案和算法； 编写上机实验程序，作好上机前的准备工作； 上机调试程序，算各种方案，记录计算的结果（包括必要的中间结果）； 分析和解释计算结果； 按照要求书写实验报告； 3、实验内容： 1) 给定数据如下： x ： 0.15，0.4，0.6 ，1.01 ，1.5 ，2.2 ，2.4， 2.7， 2.9， 3.5 ，3.8 ，4.4， 4.6 ，5.1 ，6.6， 7.6； y ： 4.4964，5.1284，5.6931 ，6.2884 ，7.0989 ，7.5507 ，7.5106， 8.0756， 7.8708，8.2403 ，8.5303 ，8.7394， 8.9981 ，9.1450 ，9.5070，9.9115； 试作出幂函数拟合数据。 2) 已知一组数据： x ： 0，0.1，0.2 ，0.3 ，0.4 ，0.5 ，0.6，0.7，0.8，0.9 ，1 y : -0.447，1.978，3.28 ，6.16 ，7.08 ，7.34 ，7.66，9.56，9.48，9.30 ，11.2； 试用最小二乘法求多项式函数，使与此组数据相拟合。 4、题目： 曲线拟合的最小二乘法 5、原理：由已知的离散数据点选择与实验点误差最小的曲线 称为曲线拟合的最小二乘法。 若记 上式可改写为这个方程成为法方程，可写成距阵形式 其中 。 它的平方误差为： ? 6、设计思想： 从几何意义上讲，就是寻求与给定点 (i=0,1,…,m)的距离平方和为最小的曲线。函数 称为拟合 函数或最小二乘解，求拟合函数 的方法称为曲线拟合的最小二乘法。#includestdio.h #includemath.h void main() { double a0[16][2],a1[2][16],A[2][2],Y[2];//a1是a0的转置，A表示 a0*a1,Y表示a1*y[j] int i,j,k; double x[16]={0.15,0.4,0.6,1.01,1.5,2.2,2.4,2.7,2.9,3.5,3.8,4.4, 4.6,5.1,6.6, 7.6}; double y[16]={4.4964,5.1284,5.6931,6.2884,7.0989,7.5507,7.5106, 8.0756,7.8708,8.2403,8.5303,8.7394,8.9981,9.1450,9.5070,9.9115}; double m0,m1,n; //m1是幂函数X的系数，n是X的指数 for(i=0;i=15;i++) { a0[i][0]=1; a0[i][1]=log(x[i]); y[i]=log(y[i]); printf(输入X的值： ); printf(%f ,a0[i][1]); printf(得到的对应的函数值：); printf(%f \n ,y[i]);} printf(\n); for(i=0;i=15;i++) for(j=0;j=1;j++) {a1[j][i]=a0[i][j];}//以上正确 A[0][0]=0;A[0][1]=0;A[1][0]=0;A[1][1]=0;Y[0]=0;Y[1]=0; for(i=0;i=1;i++) { for(j=0;j=1;j++) {for(k=0;k=15;k++) {A[i][j]+=a1[i][k]*a0[k][j];} } } for(i=0;i=1;i++) { for(j=0;j=15;j++) {Y[i]+=a1[i][j]*y[j];} } m0=(Y[0]*A[1][1]-Y[1]*A[0][1])/(A[0][0]*A[1][1]-A[1][0]