第13讲 插值解析.ppt

* 数学建模与数学实验 插 值 一、什么是插值 所谓插值就是在已知离散数据的基础上补插连续函数，使得这条连续曲线（或曲面）通过全部已知的离散数据点，利用插值方法可通过函数在有限个点处的取值状况，估算出函数在其他点处的近似值。 * 拉格朗日插值 分段线性插值 三次样条插值 一 维 插 值 一、插值的定义 二、插值的方法 三、用Matlab解插值问题 * 二维插值 一、二维插值定义 二、网格节点插值法 三、用Matlab解插值问题 最邻近插值 分片线性插值 双线性插值 网格节点数据的插值 散乱数据的插值 一维插值问题的数学描述为： 已知某一函数 y = g(x)（ g(x)的解析表达式可能十分复杂, 也可以是未知的）在区间[a, b]上n +1个互异点 xj 处的函数值 y j ， j = 0,1,…,n ，还可能知道g(x)在[a, b]上有若干阶导数，如何求出g(x)在[a, b]上任一点x*的近似值。 二、一维插值问题的数学描述 １．数学描述 ? ? ? ? ? ? 一维插值方法的基本思想是：根据g(x)在区间[a,b]上 n +1个互异点 x j （称为节点）的函数值 y j ，j = 0,1,…,n，求 一个足够光滑、简单便于计算的函数f (x)（称为插值函数） 作为g(x)的近似表达式，使得 f (x j ) = y j ， j = 0，1，…，n。 （1） 然后计算f (x)在区间[a,b](称为插值区间)上点x（称为插值 点）的值作为原函数g(x)（称为被插函数）在此点的近似 值。求插值函数f (x)的方法称为插值方法，（1）式称为插值条件。 代数多项式比较简单，常用多项式作为插值函数。 2．基本思想 ? ? ? ? ? ? 常用的一维插值方法有：拉格朗日（Lagrange）多项式插值、分段线性插值、三次样条插值、牛顿（Newton）插值、Hermite插值、最近邻插值和B样条插值等。 3．常用一维插值方法 一、Lagrange（拉格朗日）插值 1．拉格朗日线性插值 插值多项式： 直线的两点式表达式 分别称为节点x0和x1的一次 插值基函数 插值函数为基函数的线性组合，组合系数就是对应节点上的函数值 在 上，过点 2. 拉格朗日二次插值 二次插值多项式： 对于 ，可假定： 基函数满足： 基函数如何确定？ 从而 过点 构造一组插值基函数（n次多项式） 缺点：当 n 比较大时，插值多项式Ln (x)的收敛性与稳定性变 差，逼近效果不理想。 li(x)是 n次多项式，且满足 由插值基函数li(x)构造 n次拉格朗日插值多项式如下： 3. 拉格朗日n 次插值 拉格朗日n次插值的误差估计 拉格朗日插值产生的截断误差为Rn(x)，则 * 拉格朗日多项式插值的这种振荡现象叫 Runge现象 【例1】Runge（龙格）现象：用函数 在区间[-1, 1] 上产生11个等距节点，然后调用自编lagrange函数作拉格朗日插值。 * 分段线性插值 ? ? ? ? ? ? xj xj-1 xj+1 x0 xn x o y 是一个细的、可弯曲的木制或塑料条，在飞机或轮船等的设计制造过程中为描绘出光滑的外形曲线(放样)所用的工具 三次样条本质上是一段一段的三次多项式拼合而成的曲线，在拼接处,不仅函数是连续的,且一阶和二阶导数也是连续的，1946年,Schoenberg将样条引入数学,即所谓的样条函数 1. 什么是样条？ 三次样条插值的数学描述 三次样条插值 2. 三次样条插值函数 3. 三次样条插值原理 在n个小区间构造S(x)，共有n个三次多项式，需确定4n个参数 在所有节点上 n+1个方程 在除端点外的节点上 3(n-1)个方程 这样就得到4n - 2个方程，为保证待定参数的唯一性，还差两个方程。为此，常用的方法是对边界节点除函数值外附加要求，这就是所谓的边界条件。根据实际问题的不同，三次样条插值常用到下列三类边界条件。 周期边界条件： 当y = g(x)是以 b-a= x0 - xn 为周期的周期函数时，要求S(x) 也是周期函数，故端点处要满足 此条件称为周期条件。 m边界条件： 即给定端点处的一阶导数值。 M边界条件： 即给定端点处的二阶导数值。 * 比分段线性插值更光滑。 ? ? ? ? ? ? ? ? ? x y xi-1 xi a b * 用MATLA