计算方法-5数据拟合方法.pptVIP

1 * * * * * * * * * * * * * * * * * * * * * * * * 计 算 方 法 主 讲： 孙 在 Email: sunzai@cjlu.edu.cn 2009 第5章 数据拟合方法 §5–1 最小二乘原理 §5–2 超定方程组的最小二乘解 §5–3 其它函数的拟合 §5–4 Bezier曲线及应用 学习要点 最小二乘法的概念； 曲线的最小二乘拟合方法； 5.1 最小二乘原理 实验测量得到多组数据，而数据之间存在一定的函数关系，为了求得逼近函数与已给函数，同时使得其偏差按某种方法度量能达到最小，可以采用最小二乘法。 例1 测得铜导线在温度tj时的电阻rj，求出电阻r与温度t的近似表达式。 j 1 2 3 4 5 6 7 tj/℃ 19.1 25.0 30.1 36.0 40.0 45.0 50.0 rj/Ω 76.30 77.80 79.25 80.80 82.35 83.90 85.10 电阻r和温度t的数据关系 可以看到电阻与温度呈现类似直线的线性关系。 设此直线方程为： r=a+bt a,b待定。 误差为： Rj=a+btj-rj j=1,2…7 选择a，b使得Rj的平方和最小，即 用这种方法来求a,b的原理就是一种最小二乘原理 线性相关 定义5.1　设x1,x2,…xn为互不相同的点，函数f0(x), f1(x) ..f0(x), fn (x) , 是n+1个已知函数 如果存在不全为零的常数c0,c1,…cn使得 则称函数 f0(x), f1(x) ..f0(x), fn (x)关于点x1,x2,…xn是线性相关的，否则称线性无关。 定义5.2　给定数据(xj,yj),j=1,2,…n。假设拟合函数的形式为 这里 为已知线性无关函数。求系数 使得 取最小值，称 为拟合函数，或经验公式。 是 的m+1元二次多项式，求其极值方法为： 设 内积标识 ,称为点集函数。记 称为u和v的内积。有三个性质： 利用内积的定义，前式可以表达为： 上方程组称为正规方程组，是关于 的m+1阶线性方程组。可以用各种消去法、矩阵分解、和迭代法求解，然后得到拟合函数： 这个方程组也是对称矩阵，也可以写为： 对例1中： r=a+bt 则： 正规方程组为： 可以求得a=70.57，b=0.2915，代入得到电阻与温度的拟合曲线方程： r=70.57+0.2915t r=70.57+0.2915t 两个问题？ （1）正规方程组是否有解； （2）由正规方程组得到的系数 是否能使 取得最小值 例2 设某一实验中测得的两个吧变量x和y的一组数据如下，求一代数多项式曲线，使其能最好地拟合实验数据。 i 1 2 3 4 5 6 7 8 9 xi 1 3 4 5 6 7 8 9 10 yi 10 5 4 2 1 1 2 3 4 可以看到x与y呈现类似抛物线关系。 设拟合曲线方程为： y=a+bx+cx2 正规方程组为： 分别计算各项： 代入正规方程组为： 可以求得a=13.4609，b=-3.60585，c=0.267616 代入得到拟合曲线方程： y=13.4609- 3.60585x+0.267616x2 y=13.4609- 3.60585x+0.267616x2 5.2 超定方程组的最小二乘解 对线性方程组 Ax=b 若mn，此方程组称为超定方程组 通常对任意一组xj（j=1..n）： 不会全为零。可以寻求一组xi*（i=1..n）,使得 取最小值。 即求 用矩阵表达即为： ATAx=ATb 当A为列满秩时，有唯一解，此解称为超定方程组的最小二乘解 例3 用最小二乘法求下列超定方程组的近似解： 得方程组： 解得：x1=0.792 72， x2=-1.464 1 4.1, 4.3 作业 * * * * * *