差值与拟合3讲解.pptVIP

插值与拟合 数学建模工作室 杨高飞 2019/5/8 插值与拟合 ? 插值和拟合的定义 ? 常用插值方法 ? 拟合的原理 ? matlab 中的拟合 ? matlab 拟合工具箱使用 1 2 3 4 5 6 7 5 10 15 20 25 30 35 1 2 3 4 5 6 7 5 10 15 20 25 30 35 插值法：求过已知有限个数据点的近似函数。 拟合法：已知有限个数据点，求近似函数，不要求过已知 数据点，只要求在某种意义下它在这些点上的总偏差最小 。 在实际应用中两种方法的选择 1 、如果给定的数据是少量的且被认为是严格精 确的，那么宜选择插值方法。 2 、如果给定的数据时大量的测试和统计的结果 ，并不是必须严格遵守的，而是起定性地控制 作用的，那么宜选用数据拟合的方法。 常用的插值方法 ? 拉格朗日多项式插值 ? 牛顿插值 ? 分段线性插值 ? Hermite 插值 ? 三次样条插值 插值的维数 ? interp1 ： 一维插值； ? interp2 ： 二维插值； ? interp3 ： 三维插值； ? interpn ： n 维插值； matlab 中插值的函数 ? 拉格朗日多项式插值 Matlab 中没有现成的 Lagrange 插值函数，必须编写 一个 M 文件实现 Lagrange 插值 ? 牛顿插值（实际中比较少运用） ? Hermite 插值 Matlab 中没有现成的 Hermite 插值函数，必须编写 一个 M 文件实现插值。 ? 三次样条插值 拉格朗日插值法 ? 怎样构造插值函数 Ln(x)? ? 从 n=1,n=2 推广到一般情况 . 当 n=1 ，线性插值 构造方程： 化简： 令： 则： 其中： 0 1 0 1 0 0 x x y y x x y y ? ? ? ? ? 1 ) ( , 0 ) ( ; 0 ) ( , 1 ) ( ) ( ) ( ) ( ) ( , ) ( ) ( 1 1 0 1 1 0 0 0 1 1 0 0 0 1 0 1 1 0 1 0 0 1 0 1 1 0 1 0 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? x l x l x l x l x l y x l y x L x x x x x l x x x x x l x x x x y x x x x y y x L 同理：当给出 n+1 个点时： 函数为： 其中： ) , 0 ( , ) )...( )( )...( )( ( ) )...( )( ).....( )( ( ) ( ) ( ) ( 1 1 1 0 1 1 1 0 0 n k x x x x x x x x x x x x x x x x x x x x x l x l y x L n k k k k k k k n k k k n k k k n ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? matlab 中插值的函数 ? 分段线性插值 用 Matlab 实现分段线性插值不需要编制函数程序， Matlab 中有现成的一维插值函数 interp1 。 y=interp1(x0,y0,x,method) method 指定插值的方法，默认为线性插值。其值可为： nearest 最近项插值 linear 线性插值 spline 立方样条插值 cubic 立方插值。 所有的插值方法要求 x0 是单调的 matlab 各种插值方法的比较 ? nearest ：执行速度最快，输出结果为直角转折，插 值效果较差； ? linear ：默认值，在样本点上斜率变化很大； ? spline ：最花时间，但输出结果也最平滑，十分逼近 理论值 ? cubic ：最占内存，输出结果与 spline 差不多。 matlab 中插值的函数 ? 三次样条插值 Matlab 中三次样条插值也有现成的函数： y=interp1(x0,y0,x,spline) ； y=spline(x0,y0,x) ；