大学生数学建模：数学建模插值与拟合.pptVIP

插值与拟合 主讲人：刘 婷 Email:lt1210@ 最小二乘法原理: 例. 在研究某单分子化学反应速度时, 得到下列数据: 因此 a , b 应满足法方程组: 通过计算确定某些经验公式类型的方法: 特别, 当数据点分布近似一条直线时, 问题为确定 a, b 令 满足: 使 得 解此线性方程组 即得 a, b 称为法方程组 (注意其特点) 例. 为了测定刀具的磨损速度, 每隔 1 小时测一次刀 具的厚度, 得实验数据如下: 找出一个能使上述数据大体适合的经验公式. 解: 通过在坐标纸上描点可看出它们 大致在一条直线上, 列表计算: 故可设经验公式为 27.0 26.8 26.5 26.3 26.1 25.7 25.3 24.8 0 1 2 3 4 5 6 7 0 1 2 3 4 5 6 7 得法方程组 解得 故所求经验公式为 0 0 27.0 0 7 49 24.8 137.6 28 140 208.5 717.0 为衡量上述经验公式的优劣, 计算各点偏差如下: 称为均方误差, 对本题均方误差 它在一定程度上反映了经验函数的好坏. 偏差平方和为 27.0 26.8 26.5 26.3 26.1 25.7 25.3 24.8 0 1 2 3 4 5 6 7 27.125 26.518 25.911 25.303 26.821 26.214 25.607 25.000 －0.125 －0.018 0.189 －0.003 －0.021 0.086 0.093 －0.200 57.6 41.9 31.0 22.7 16.6 12.2 8.9 6.5 3 6 9 12 15 18 21 24 1 2 3 4 5 6 7 8 其中? 表示从实验开始算起的时间, y 表示时刻 ? 反应 物的量. 试根据上述数据定出经验公式 解: 由化学反应速度的理论知, 经验公式应取 其中k , m 为待定常数. 对其取对数得 (线性函数) (书中取的是常用对数) 经计算得 解得: 所求经验公式为 其均方误差为 观测数据: 用最小二乘法确定a, b 线性最小二乘拟合 f(x)=a1r1(x)+ …+amrm(x)中函数的选取 1. 通过机理分析建立数学模型来确定 f(x)； + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + f=a1+a2x f=a1+a2x+a3x2 f=a1+a2x+a3x2 f=a1+a2/x f=aebx f=ae-bx 2. 将数据 (xi,yi) i=1, …n 作图，通过直观判断确定 f(x)： 1、线性最小二乘拟合 2、非线性最小二乘拟合 用MATLAB作线性最小二乘拟合 1. 作多项式f(x)=a1xm+ …+amx+am+1拟合,可利用已有程序: a=polyfit(x,y,m) 2. 对超定方程组 可得最小二乘意义下的解。 ，用 3.多项式在x处的值y可用以下命令计算：y=polyval（a，x） 输出拟合多项式系数 a=[a1, …am , am+1] (数组)） 输入同长度 的数组X，Y 拟合多项 式次数 即要求 出二次多项式: 中 的 使得: 例 对下面一组数据作二次多项式拟合 1）输入以下命令： x=0:0.1:1; y=[-0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2]; A=polyfit(x,y,2) z=pol