数值计算方法课件 第四章 多项式插值与函数逼近VIP

性质2（最佳逼近性质）设是任一被插函数，是带有斜率边界条件的三次插值样条函数，是与有相同分割的任一三次样条函数，则有证明：证明方法与性质1类似（自己证明）性质3（误差估计）设函数，是区间的一个分割，是关于的带有Ⅰ型(斜率边界)或Ⅱ型(二阶导数边界)边界条件的插值函数，则有误差估计其中是分割比，并且系数与是最优估计。性质说明：三次样条插值函数本身连同它的一、二、三阶导数分别收敛到及其相应导数，具有强收敛性。广义多项式§6 函数逼近/* Approximation of Function */ 一、函数逼近问题的提法假设是定义在某区间上的函数，现寻求另一个构造简单、计算量小的函数来近似地代替：为区间上的一个线性无关函数系为一组实常数。就是我们前面讨论的多项式逼近若线性无关函数系取常用的函数系：幂函数系：三角函数系：指数函数系：多项式插值的收敛性问题：注：结论并不是对所有函数都成立函数逼近构造思想：要求构造函数在整个区间上与已知函数的误差尽可能小误差度量标准：其中为权函数（2）（1）对于给定的函数系，寻求一组系数使得函数满足（1）（2）一致逼近逼近二、最佳平方逼近/*Best Approximation in Quadratic Norm*/ 假设，是[a,b] 上的一个线性无关函数系,且, 为[a,b] 上的一个权函数如果存在一组系数使得广义多项式满足称函数为在[a,b] 上关于权函数的最佳平方逼近或最小二乘逼近；特别，若，则称是在[a,b] 上的最佳平方逼近. 由定义可以看出，最佳平方逼近问题实际上是个多元极值问题记由极值的必要条件即：记将代入前式：令称矩阵是关于函数系的Gram( 格拉姆)矩阵易证Gram 矩阵为实对称正定矩阵：上述方程组存在唯一解设由上述方程组的解确定的广义多项式为：对于任意广义多项式下面证明即记设给定节点，则其最佳平方逼近唯一存在，且可以由前述Gram 组成的方程组求解构造。注：前述Gram 组成的方程组通常称为法方程组最佳平方逼近可以通过求解法方程组而得到Gram 矩阵是实对称正定矩阵一、三次样条函数的力学背景在工程技术和数学应用中经常遇到这样一类数据处理问题：在平面上给定了一组有序的离散点列，要求用一条光滑曲线把这些点按次序连接起来。. . . . . . . . 压铁弹性木条. 数据点形象地称之为样条曲线在力学上，通常均匀细木条可以看作弹性细梁，压铁看作是作用在梁上的集中载荷，“样条曲线”就模拟为弹性细梁在外加集中载荷作用下的弯曲变形曲线。设细梁刚度系数为，弯矩为，样条曲线的曲率为由力学知识：当时（即“小挠度”的情况）上述微分方程简化为：是线性函数因此，“样条曲线”可近似认为是三次多项式二、三次样条函数定义及求法设在区间上给定一个分割，定义在上的函数如果满足下列条件：(1) 在每个小区间内是三次多项式(2) 在整个区间上，为二阶连续可导函数，即在每个节点处则称为三次样条函数假设现在已知函数在节点处的函数值：如果三次样条函数满足则称为插值于的三次样条函数，简称三次样条插值函数。如何求的三次样条插值函数: 4n 个未知数3n-1 个条件线性插值函数1、M连续方程与的表达式记因为在每一个子区间上都是线性函数两边积分两边再积分一次？代入插值条件：在整个区间上，的表达式为：未知数n+1 个的计算方法：由由由其中写成方程组的形式：上述方程组称为的M连续方程n-1 个方程n+1 个未知数三弯矩方程2、m连续方程与的表达式记在区间上采用Hermite 插值的计算方法：对两边求导（微分2次）两者相等得到方程组：其中同前写成方程组的形式：n-1 个方程n+1 个未知数三转角方程上述方程组称为的m连续方程M、m连续方程的求解：需要补充附加条件3、边界条件/*boundary conditions */ 已知端点的斜率：已知端点的二阶导数：设是以为周期的周期函数，对附加周期性条件：即要求三次样条插值函数在端点处函数值、一阶导数值和二阶导数值相同。M连续方程在各类边界条件下的求解方法对于第一类边界条件由得从而得到方程组（三对角）：可用追赶法求解对于第二类边界条件类似地可以得到方程组（三对角）：上述两种情况得到的方程组，可以写成统一形式：其中时为第二类边界条件时为第一类边界条件对于第三类边界条件得到方程其中第三类边界条件对应的方程组：对三对角算法经过修改后可以求解上述方程组不是三对角方程组注：三次样条与分段Hermite 插值的根本区别在于S(x) 自身光滑，不需要知道f 的导数值（除了在2个端点处的函数值）；而Hermite 插值依赖于f 在许多插值节点的导数值。f(x) H(x) S(x) 3 2 6 2 0 0 1 3 1 0 例1：已知函数在的数据表：解：求在区间上的三次样条插值函数。三弯矩方程组用追赶法求