计算方法(8)--第五章-插值法(2).pptVIP

* * * * * 第(3)步 f [xi , xi+1] f [xi-1 , xi] 计算步骤 例5.5 给出函数表 试求f (x) 在区间[0,3]上的三次样条插值函数 用节点处二阶导数表示的三次样条插值函数 6f [xi-1 , xi , xi+1 ] 计算步骤 ? 分段插值常用命令：?? interp1: 一维分段插值 interp2: 二维分片插值 interp3: 三维分块插值 spline: 三次样条插值 一维分段插值 y0=interp1(x,y,x0,插值方法选项) plot(x,y,x0,y0):绘插值曲线 ‘插值方法选项’有以下四种 （1）nearest: 最接近点插值 （2）linear: 分段线性插值 （3）cubic: 分段三次插值 （4）spline:分段三次样条插值 本章重点 1、拉格朗日插值 2、牛顿插值 3、三次样条插值的定义 习题五 P133-----11, 12（2） 试判断下面函数是否为三次样条函数： ① ② ② ② * Hermite插值是一个最简单的密切插值多项式 * * * * * * * * * * * * * * * Hermite插值的提出 假设函数y=f(x)是 在[a,b]上有一定光滑性的函数,在[a,b] 上有n+1个互异点xo…xn, f(x)在这些点上取值yo…...yn.求一个确定的函数p(x)在上面n+1个点上满足p(xi)=yi i=0,1,…,n.这是最简单的插值问题,如果除了知道f(x)在插值节点上的取值外,还知道f(x)在插值节点xi上的 1≤mi≤n阶导数，如何来构造插值函数呢? Hermite插值就是既满足插值节点xi的函数值条件又满足微商条件的插值函数。 5.4 埃尔米特(Hermite)插值 一、一般情形的 Hermite插值问题 一、一般情形的 Hermite插值问题 实际应用中常采用分段低次插值。 （1）分段线性插值 （2）分段二次插值与分段三次插值 （3）分段Hermite插值 (4) 分段三次样条插值 一、分段低次插值 5.5三次样条插值 1.分段线性插值 在每个子区间[xi ,xi+1](i=0,1,2,???,n-1)上是一次插值多项式; 分段线性插值曲线图： 分段线性插值函数的余项 注意: h随分段增多而减少，因此用分段插值提高精 度是很好的途径. 定理：设f(x)在[a,b]上有二阶连续导数f″(x) ， 且| f″(x)| ≤m2, 记： h = max |xi+1-xi|,就有估计： |R(x)| =|f(x)- ?(x) |≤m2h2/ 8 , x∈[a, b]。 2.分段二次插值与分段三次插值 例: 在[-4,4]上给出等距节点函数表，若用分段二次插值计算ex的近似值，要使截断误差不超过10-6,问使用函数表的步长h 应为多少? 解：设 xi-1≤x≤xi+1, 则有 xi-1=xi-h, xi+1=xi+h, x=xi+th (-1≤t≤1) 过三点 xi-1，xi，xi+1的二次插值误差为： 三次样条插值函数的定义 一般的分段插值存在着一个缺点,就是会导致插值函数在子区间的端点(衔接处)不光滑,即导数不连续,但对于一些实际问题却不许在插值节点处有一阶和二阶导数的间断,例如高速飞机的机翼外形,内燃机进排气门的凸轮曲线以及船体放样等等。为了满足这些要求,人们引入了样条插值的概念。 所谓“样条”(spline)是工程绘图中的一种工具，它是有弹性的细长木条。绘图时,往往用样条在样点用压铁顶住，样条在自然弹性弯曲下形成的光滑曲线称之为样条曲线,此曲线不仅具有一阶连续导数，而且还有连续的曲率（即具有二阶连续导数）。可以证明，此曲线是分段的三次曲线。 样条函数就是对这样的曲线进行数学模拟得到的。它除了要求给出各个结点处的函数值外,只需提供两个边界点处导数信息,便可满足对光滑性的要求。 二、三次样条插值 定义5.4 三次样条插值函数的构造 (一)用节点处的一阶导数表示的三次样条插值函数 第(1)步 第(2)步 * Hermite插值是一个最简单的密切插值多项式 * * * * * * * * * * * * * * * * * * * *