曲线拟合最小二乘法.docVIP

数学与计算科学学院 实 验 报 告 实验项目名称 曲线拟合最小二乘法 所属课程名称 数值计算 实 验 类 型 验证实验 实 验 日 期 班 级 学 号 姓 名 成 绩 一、实验概述： 【实验目的】 掌握最小二乘法， 熟练掌握求函数的近似最小二乘拟合多项式 【实验原理】 在科学实验数据处理中，往往要根据一组给定的实验数据，求出自变量x与因变量y的函数关系，这是为待定参数，由于观测数据总有误差，且待定参数ai的数量比给定数据点的数量少(即n＜m)，因此它不同于插值问题.这类问题不要求通过点，而只要求在给定点上的误差的平方和最小.当时，即(3.1)这里是线性无关的函数族，假定在上给出一组数据，以及对应的一组权，这里为权系数，要求使最小，其中(3.2) 这就是最小二乘逼近，得到的拟合曲线为y=s(x)，这种方法称为曲线拟合的最小二乘法.　　(.2)中实际上是关于的多元函数，求I的最小值就是求多元函数I的极值，由极值必要条件，可得　　　　　　(.3)根据内积定义(见第三章)引入相应带权内积记号(3.4)则(3.3)可改写为　这是关于参数的线性方程组，用矩阵表示为 (3.5)(3.5)称为法方程.当线性无关，且在点集上至多只有n个不同零点，则称在X上满足Haar条件，此时(3.5)的解存在唯一(证明见［3］).记(3.5)的解为　　　　　　　　 从而得到最小二乘拟合曲线　　　　　　　　　　(3.6)可以证明对，有　　　　　故(3.6)得到的即为所求的最小二乘解.它的平方误差为　　　　　　　　　　　　　　　　　　(3.7)均方误差为　　　　　　　　　　　　在最小二乘逼近中，若取，则，表示为　　　　　　　　　　　　　(3.8)此时关于系数的法方程(3.5)是病态方程，通常当n≥3时都不直接取作为基. y=2.014000e+000+2.250000e+000x p171习题3： y=-1.071071e+001+1.404405e+001x p280数值试验六-1-拟合多项式：y=4.957055e-001+2.463190e+000x 【实验结论（结果）小结结论#includestdio.h #includemath.h #define N 5 void main() { long double x[N]={-1,-0.5,0,0.5,1}; long double y[N]={-0.22,0.88,2.00,3.13,4.28}; long double u11=0,u12=0,u21=0,u22=0,c1=0,c2=0; long double A,B; int i; for(i=0;iN;i++) { u12+=x[i]; u22+=x[i]*x[i]; c1+=y[i]; c2+=x[i]*y[i]; } u21=u12; u11=N; A=(c1*u22-c2*u12)/(u11*u22-u21*u12); B=(c1*u21-c2*u11)/(u12*u21-u22*u11); printf(y=%le+%lex\n,A,B); } p171习题3： #includestdio.h #includemath.h #define N 8 void main() { long double x[N]={1,2,3,4,5,6,7,8}; long double y[N]={15.3,20.5,27.4,36.6,49.1,65.6,87.8,117.6}; long double u11=0,u12=0,u21=0,u22=0,c1=0,c2=0; long double A,B; int i; for(i=0;iN;i++) { u12+=x[i]; u22+=x[i]*x[i]; c1+=y[i]; c2+=x[i]*y[i]; } u21=u12; u11=N; A=(c1*u22-c2*u12)/(u11*u22-u21*u12); B=(c1*u21-c2*u11)/(u12*u21-u22*u11); printf(y=%le+%lex\n,A,B); } p280数值试验六-1-拟合多项式实验： #includestdi